- 1Department of Civil and Environmental Engineering, Michigan State University, East Lansing, MI, United States
- 2Great Lakes Institute for Environmental Research, University of Windsor, Windsor, ON, Canada
- 3Great Lakes Center for Fresh Waters and Human Health, Bowling Green State University, Bowling Green, OH, United States
- 4Great Lakes Water Authority, Detroit, MI, United States
Introduction: Wastewater surveillance has proven to be a valuable approach to monitoring the spread of SARS-CoV-2, the virus that causes Coronavirus disease 2019 (COVID-19). Recognizing the benefits of wastewater surveillance as a tool to support public health in tracking SARS-CoV-2 and other respiratory pathogens, numerous wastewater virus sampling and concentration methods have been tested for appropriate applications as well as their significance for actionability by public health practices.
Methods: Here, we present a 34-week long wastewater surveillance study that covers nearly 4 million residents of the Detroit (MI, United States) metropolitan area. Three primary concentration methods were compared with respect to recovery of SARS-CoV-2 from wastewater: Virus Adsorption-Elution (VIRADEL), polyethylene glycol precipitation (PEG), and polysulfone (PES) filtration. Wastewater viral concentrations were normalized using various parameters (flow rate, population, total suspended solids) to account for variations in flow. Three analytical approaches were implemented to compare wastewater viral concentrations across the three primary concentration methods to COVID-19 clinical data for both normalized and non-normalized data: Pearson and Spearman correlations, Dynamic Time Warping (DTW), and Time Lagged Cross Correlation (TLCC) and peak synchrony.
Results: It was found that VIRADEL, which captures free and suspended virus from supernatant wastewater, was a leading indicator of COVID-19 cases within the region, whereas PEG and PES filtration, which target particle-associated virus, each lagged behind the early alert potential of VIRADEL. PEG and PES methods may potentially capture previously shed and accumulated SARS-CoV-2 resuspended from sediments in the interceptors.
Discussion: These results indicate that the VIRADEL method can be used to enhance the early-warning potential of wastewater surveillance applications although drawbacks include the need to process large volumes of wastewater to concentrate sufficiently free and suspended virus for detection. While lagging the VIRADEL method for early-alert potential, both PEG and PES filtration can be used for routine COVID-19 wastewater monitoring since they allow a large number of samples to be processed concurrently while being more cost-effective and with rapid turn-around yielding results same day as collection.
1. Introduction
Wastewater surveillance has been widely adopted by researchers and health agencies as an effective tool for tracking Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) in wastewater amid the Coronavirus Disease 2019 (COVID-19) pandemic (1–13). SARS-CoV-2 was first identified in Wuhan, Hubei, China, and was designated a Public Health Emergency of International Concern on January 30th, 2020, by the World Health Organization (WHO). COVID-19 was later declared a pandemic on March 11th, 2020 (who.int). Numerous studies have demonstrated that SARS-CoV-2 can be shed from the gastrointestinal tract of infected individuals and its viral RNA can persist and be detected in wastewater (14–18). To increase the sensitivity of the assay used to detect viral RNA in wastewater, samples are routinely concentrated prior to quantification (19–21).
Methods used in published studies to recover and concentrate SARS-CoV-2 viral RNA from wastewater encompass a wide range of techniques including Virus Adsorption-Elution (VIRADEL), polyethylene glycol precipitation (PEG), ultrafiltration, ultracentrifugation, concentrating pipette, filtration and so forth. Some of the methods, such as VIRADEL, exclude large solids and focus on free and suspended viral particles in supernatant wastewater. Other methods, such as PEG precipitation and filtration, target particulate matter and the associated viruses that are sorbed onto solids. Notably, this fraction may preferentially settle within the sewer when flow is reduced and likewise is susceptible to resuspension when flows are elevated (3, 22).
The recovery efficiencies of concentration methods are variable, differing between method, virus type and conditioning of the wastewater sample. Notably, VIRADEL was found to be effective for concentrating viruses from water samples with recovery efficiencies of more than 90% for poliovirus (23, 24), 54.4% for murine norovirus (MNV) (25), 51% for echovirus (26), 35% for enteric virus (27), and 4.7% for adenovirus (28). Likewise, PEG was found to be effective for concentrating viruses in water samples, with recovery efficiencies of 89.5% for echovirus (29), 86% for hepatitis A virus (30), 68% for poliovirus (30), and 56.7% (31) and 26.4% (32) for SARS-CoV-2. Filtration was reported to recover virus from wastewater samples with recovery efficiencies ranging from 26.7 to 65.7% for murine hepatitis virus (33), and 90% for human betacoronavirus OC43 (34).
Applying different concentration methods can achieve different goals. For instance, use of VIRADEL to concentrate SARS-CoV-2 can provide early warnings of impending COVID-19 cases (1, 3, 13). PEG precipitation is an economical and widely adopted method that allows a large number of samples to be processed concurrently and it is suitable for routine COVID-19 wastewater monitoring (22, 35). Likewise, filtration presents a cost-effective and simple approach commonly applied to recover cells and viral particles from environmental samples for nucleic acid extraction (36), which has also been applied to recovery of SARS-CoV-2 from wastewater (12, 35, 37, 38).
Here we present a comparison of three primary concentration methods (VIRADEL, PEG and filtration) to detect SARS-CoV-2 viral RNA in wastewater, in relation to COVID-19 cases amid the transition from Delta to Omicron Variants of Concerns (VOCs) circulating in the Detroit, MI metropolitan area. Similarities and correlations were examined among the three concentration methods with both normalized and non-normalized data. The lead/lag time of each method in relation to the total COVID-19 cases was also assessed. The results presented in this study will assist researchers and public health practitioners to select appropriate primary concentration methods for the recovery of SARS-CoV-2 from wastewater for different wastewater surveillance practices.
2. Materials and methods
Untreated wastewater samples were collected weekly from the Water Resource Recovery Facility (WRRF) of the Great Lakes Water Authority (GLWA) located in Detroit, MI, United States, between October 1, 2021, and May 31, 2022. The WRRF serves the needs of Detroit and 76 area communities with a service area of more than 2,450 square kilometers serving nearly 4 million people. WRRF collects and treats stormwater, as well as residential, industrial, and commercial waste, depending on service areas, with its semi-combined sewershed system. WRRF receives wastewater via three main interceptors including the Detroit River Interceptor (DRI), the North Interceptor-East Arm (NIEA), and the Oakwood-Northwest-Wayne County Interceptor (ONWI) (Figure 1), serving the City of Detroit as well as the three largest Michigan counties by population: Wayne, Oakland, and Macomb. Composite samples collected over 24-h were used to compare the polyethylene glycol (PEG) precipitation and filtration methods, however, the larger volumes required by the virus adsorption-elution (VIRADEL) method necessitated a targeted approach with samples collected between 15:30 to 18:00 each afternoon. The samples were collected from the three interceptors at the point of discharge into the WRRF and maintained chilled on ice during transport to the lab for primary concentration and sample analysis.
2.1. Virus adsorption-elution method
The United States Environmental Protection Agency virus adsorption-elution (VIRADEL) method employing electropositive or electronegative filters was reported to recover and concentrate viruses from wastewater samples previously (1–4, 13, 20, 35). Electronegative filters require preconditioning such as adjusting the pH, prior to downstream concentration processes. Electropositive filters do not require any preconditioning (20, 35). In this study, depending on the quantity of suspended solids in the wastewater, 10 to 50 L of untreated wastewater (grab sample) was passed through NanoCeram electropositive cartridge filters (Argonide, Sanford, FL, United States) at a rate less than 11 L/min using a previously described method (1–4). Flow meter readings were tracked at the beginning and end of each sampling event to measure the total volume of wastewater passing through the filters. Following sampling, the NanoCeram filters were transported on ice to the lab for sample analysis within 24 h. The elution process releases viral particles captured by the filters (20). Viruses were eluted using 1.5% beef extract containing 0.05 M glycine, based on a previously described method (1–4). Subsequently, the eluates containing viruses were flocculated by adjusting the pH, following multiple centrifugations and resuspension of particles in sodium phosphate. Afterwards, supernatants containing viruses were separated by adjusting the pH and centrifugation. Finally, the supernatants containing viruses were passed through 0.45 μm and 0.22 μm Millipore filters (MilliporeSigma, Burlington, MA, United States), which were followed by aliquoting and storage of the final aliquots at -80°C for downstream molecular analysis (1–4, 20). Bacteriophage Phi6 was applied as a proxy virus to estimate the recovery efficiency during virus concentration (3, 29, 39). Figure 2 demonstrates the workflow of the VIRADEL method.
2.2. Polyethylene glycol precipitation method
From a 24-h composite sample of untreated wastewater collected in a 1 L Nalgene bottle, 100 mL samples were mixed with 0.2 M sodium chloride and 8% polyethylene glycol (w/v). Samples were mixed gently on a magnetic stirrer at 4°C for 2 h, followed by centrifugation at 4700 × g for 45 min at 4°C. The supernatant was removed, and the pellet was resuspended in the remaining liquid (approximately 2-3 mL). The final concentrate volumes were between 1 to 6 mL. All sample concentrates were then subjected to downstream analysis including RNA extraction and RT-ddPCR (Figure 3).
2.3. Filtration method
Composite samples of raw wastewater collected as for the PEG method were concentrated by filtering 50-120 mL through 0.22 μm Sterivex PES cartridge filters (MilliporeSigma, Burlington, MA, United States) using a 50 mL syringe fitted into a caulking gun. Immediately following filtration, the filters were sealed and flash-frozen through immersion in liquid nitrogen as described previously (38). Subsequently, filters were subjected to downstream processes including RNA extraction and RT-qPCR (Figure 4).
2.4. RNA extraction, RT-ddPCR, RT-qPCR
Following VIRADEL and PEG methods, viral RNA was extracted using the QIAamp Viral RNA kit (Qiagen, Germantown, MD, United States), following the manufacturer’s protocol modified by use of 140 μL elution buffer to extract the viral RNA (1–4). RT-ddPCR was performed on a QX200 AutoDG Droplet Digital PCR system (Bio-Rad, Hercules, CA, United States), using the One-step RT-ddPCR Advanced Kit for Probes (Bio-Rad, Hercules, CA, United States) as described previously (2, 3). United States Centers for Disease Control and Prevention (US CDC) primers and probes that target the N1 and N2 genes of SARS-CoV-2 were used (2, 3, 13). N1 N2 gene Duplex Assay Reaction Mixture was reported previously (2, 3, 13). Following the preparation of the Duplex Mixture and oil droplets generation, samples were run on a C1000 Touch Thermal Cycler (Bio-Rad, Hercules, CA, United States) using the thermocycling conditions which were reported previously (2, 3, 13). Subsequently, the measurement of fluorescence was performed on a QX200 Droplet Reader (Bio-Rad, Hercules, CA, United States). For each RT-ddPCR run, positive controls (PTCs), negative controls (NTCs), and process negative controls were included, which were described previously (3). All samples were run in triplicate. The Limit of Detection (LOD) and Limit of Blank (LOB) for RT-ddPCR were described and determined previously (2, 3, 13).
Following the filtration method, filters were thawed, and RNA was extracted from the filters using the AllPrep PowerViral DNA/RNA kit (Qiagen, Germantown, MD, United States) modified by addition of 5% 2-mercaptoethanol (v/v). RNA was eluted in 50 μL of RNAse free water. Samples were not treated with DNase upon extraction. Assays for SARS-CoV-2 targeted regions of the nucleocapsid (N) gene using US CDC primers and probes for the N1 and N2 regions (40). Reagents were supplied by Integrated DNA Technologies (Coralville, IA, United States). Reactions contained 5 μL of RNA template mixed with 10 μL of 2 × RT-qPCR master mix (Takyon TM Dry One-Step RT Probe MasterMix No Rox, Eurogentec, Liège, Belgium) and primers and probes in a final reaction volume of 20 μL. Reaction inhibition was assessed using VetMAX XENO Internal Positive Control RNA (Applied Biosystems Corp., Waltham, MA, United States). Due to repeated incidence of inhibition with wastewater samples processed by filtration, template was diluted 1:5 in all reactions. Technical triplicates were run for detection of gene targets. Thermal cycling was performed using a MA6000 qPCR thermocycler (Sansure Biotech, Changsha, China). RT was performed at 48°C for 10 min, followed by polymerase activation at 95°C for 3 min, and 50 cycles of denaturation, annealing/extension at 95°C for 10 s, then 60°C for 45 s, respectively. The EDX SARS-CoV-2 synthetic RNA standard (Exact Diagnostics, Fort Worth, TX, United States) was used to create a 7-point standard curve to quantify N1 and N2 gene targets. No template controls yielded no amplification, and we report a limit of detection of 5 gene copies of N1 and N2 per reaction containing 5 μL of template RNA for RT-qPCR.
2.5. COVID-19 clinical data
Publicly available clinical data were accessed on August 22, 2022, for the period between October 1, 2021, and May 31, 2022, for the city of Detroit, as well as Wayne, Macomb, and Oakland counties (Figure 5A).1 Clinical data with a 7-day moving average (3, 41, 42) was used for further statistical analysis (Figure 5B). COVID-19 clinical data were only available per city/county for the Detroit metropolitan area. Each interceptor received wastewater from portions of each city/county. Therefore, only the total SARS-CoV-2 concentrations can be correlated to the total COVID-19 cases of each city/county (3, 13).
Figure 5. (A) COVID-19 cases in the City of Detroit, as well as Wayne, Macomb, and Oakland counties; (B) 7-day moving average of the COVID-19 cases.
2.6. Data analysis and visualization
Data were tracked and organized using Microsoft Excel version 16.66.1. R version 4.1.3 was applied to perform data analysis including Pearson and Spearman correlations, Dynamic Time Warping (DTW), Time Lagged Cross Correlation (TLCC) and peak synchrony, depending primarily on ggplot2 package for visualization, and packages including dtw, synchrony, dplyr, and ggpubr. Missing data from samples were filled using linear interpolation for further analysis (3, 13, 43). For VIRADEL samples, 128 genes concentrations were measured for both N1 and N2 genes between 10/1/21 and 5/31/22. For PEG samples, 88 gene concentrations were measured for both N1 and N2 genes between 10/1/21 and 5/31/22. For filtration samples, 66 gene concentrations were measured for both N1 and N2 genes between 10/1/21 and 5/31/22. To perform correlation analysis between weekly gene concentrations and daily clinical cases, linear interpolation was conducted to generate daily data based on weekly measurements. The number of interpolated daily gene concentrations were 179, 199, and 210 for VIRADEL, PEG, and filtration, respectively.
To account for the changing flow in wastewater, dilution events, and variability in the solids portion of the wastewater, four approaches (flow rate, flow rate/population, TSS, flow rate×TSS) of normalizing the N1 and N2 gene concentrations (gc/L) were implemented using Eq. (1), Eq. (2), Eq. (3), and Eq. (4) (3, 44, 45). TSS, or “Total Suspended Solids,” is an estimate of the entire solids in wastewater in contrast to the liquid fraction or dissolved matter (45). In addition, other parameters, including sanitary percentage and Biological Oxygen Demand (BOD), proved ineffective for normalizing N1 and N2 gene concentrations for the Detroit area and other areas, thus, they were not considered in the current study (3, 9). SARS-CoV-2 gene concentrations measured in the wastewater following VIRADEL, PEG, and filtration methods are reported as gene copies per L (gc/L). The units after normalization using flow rate, flow rate/population, TSS, and flow rate×TSS, are gene copies per day (gc/day), gene copies per day per person (gc/day/person), gene copies per mg TSS (gc/mg TSS), and gene copies per L per pounds/day {gc/[L(pounds/day)]}, respectively.
is the normalized concentration of SARS-CoV-2 in gc/day. is the normalized concentration of SARS-CoV-2 in gc/day/person. is the normalized concentration of SARS-CoV-2 in gc/mg TSS. is the normalized concentration of SARS-CoV-2 in gc/[L(pounds/day)]. is the volume of wastewater flowing into WWTP interceptors during sampling events (MGD). is the conversion factor between L and MGD. is the conversion factor between mg and pounds. is the total population in the Detroit metropolitan area served by WRRF’s interceptors including ONWI, NIEA, and DRI. represents the total suspended solids (mg/L).
2.6.1. Correlations among N1 and N2 gene concentrations by VIRADEL, PEG, and filtration
Multiple studies investigated the applications of both Pearson and Spearman correlations on analyzing the relationship between wastewater viral concentrations of SARS-CoV-2 and COVID-19 clinical cases as well as the relationship among wastewater viral concentrations by different genes or methods (9, 46, 47). In this study, Pearson and Spearman correlations were performed among N1 and N2 gene concentrations {gc/L, gc/day, gc/day/person, gc/mg TSS, gc/(L[pounds/day)]} by VIRADEL, PEG, and filtration methods. The Pearson correlation measures how two time series among VIRADEL, PEG, and filtration gene concentrations covary during the study period, and indicate their linear relationships. The Spearman correlation coefficient is a simple and straightforward approach to analyze the degree of associations between two time series (48).
2.6.2. Dynamic time warping
One commonly used algorithm for quantifying the similarities/dissimilarities between time series data is the Euclidean distance (ED), but numerous studies demonstrated that ED is insensitive to time shifting and patterns between time series since it compares the data points of time series in a settled sequence and cannot consider time shifting or patterns (49, 50). Dynamic time warping (DTW) is a well-established algorithm that circumvents the limitations of ED and compares two time series by computing dynamic distances between them considering regional distortions, time shifting, and the optimal warping that best aligns the time series between each other (50, 51). Therefore, similar patterns that occur at different times between time series can be considered as matching, thus, the similarity of time series can be evaluated considering their time shifting and shapes by DTW algorithm (50). The DTW algorithm was proposed previously (51).
The outcome of DTW analysis indicates two time series with the most similar patterns by calculating the minimum overall dissimilarity or the DTW minimum distance where data points on one time series best align data points on another time series (51). Multiple studies investigated the similarities between time series using DTW algorithm (50, 52, 53). However, to our knowledge, this is the first study to apply DTW algorithm to compare the similarities between wastewater gene concentrations data by three concentration methods (VIRADEL, PEG, and filtration), as well as comparing the similarities between wastewater gene concentrations data and COVID-19 clinical data. In this study, package dtw and related packages in R (version 4.1.3) were implemented to calculate DTW for the normalized {gc/day, gc/day/person, gc/mg TSS, and gc/[L(pounds/day)]} and non-normalized (gc/L) data to analyze the similarities/dissimilarities between VIRADEL, PEG, and filtration methods.
One limitation is that the minimum DTW distance can be affected by the scaling factor of time series data. For instance, the minimum DTW distance between PEG (gc/day/person) and COVID-19 cases can be smaller than the distance between VIRADEL (gc/day/person) and cases, indicating that PEG presents higher similarity to cases than VIRADEL. However, this was affected by the population factor which is a constant number but is not dynamic time series data. Using flow/population normalization including a constant factor intentionally changed the similarities among time series data. Therefore, the minimum DTW distance with flow/population normalized data was not considered for further discussions.
2.6.3. Time lagged cross correlation and peak synchrony
To estimate the leading or lagging relationships between wastewater viral concentrations by three concentration methods (VIRADEL, PEG, and filtration) and total COVID-19 cases, TLCC and peak synchrony were performed where the total COVID-19 cases were shifted over time and correlated with wastewater viral concentrations for each concentration method. TLCC refers to correlations between two time series shifted relatively in time. It can identify the direction and relationship between two time series, for instance, a leader-follower relationship, where the leader time series develop a pattern which is repeated by the follower time series (54). TLCC is widely applied in analyzing time series especially delay, lead/lag time, and lagged cross correlation and so forth (44, 54–56). TLCC is an effective approach to estimate the dynamic relationships between two time series and demonstrate how they shift over time (44).
In this study, TLCC is measured by gradually shifting total COVID-19 cases between -20 days (lagging) and + 20 days (leading), and constantly calculating the Pearson’s correlation coefficients between two time series for each shifting. Peak synchrony occurs when the peak correlation is observed. For instance, if the peak correlation is observed at the center where the lag time or offset is 0 day, this condition indicates that the time series are most synchronized at day 0 demonstrating no shifting or lag time. However, the peak correlation can be at a different offset if one time series is leading or lagging another one. R package “synchrony,” “devtools,” and related packages were implemented to calculate the TLCC and peak synchrony between gene concentrations (both normalized and non-normalized data, by VIRADEL, PEG, and filtration methods) and 7-day moving average total COVID-19 cases.
3. Results
3.1. SARS-CoV-2 viral RNA concentrations in wastewater derived by three concentration methods spanning the transition between delta and omicron VOCs
RT-ddPCR (VIRADEL and PEG samples) and RT-qPCR (filtration samples) targeting the N1 and N2 genes was used to quantify SARS-CoV-2 RNA concentrations in wastewater samples collected at GLWA’s WRRF over 34 weeks. The study period captured the third major resurgence of COVID-19 cases in the region corresponding to the transition from SARS-CoV-2 Delta (B.1.617.2) variant to Omicron (B.1.1.529) variant (3, 44).
Filtered samples yielded N1 and N2 gene concentrations higher than those of VIRADEL but lower than those of PEG, for both normalized and non-normalized data (Table 1). Filtered samples yielded mean N1 and N2 gene concentrations of 3.22E+04 and 1.50E+04 gc/L, respectively. VIRADEL samples yielded mean N1 and N2 gene concentrations of 1.61E+03 and 1.63E+03 gc/L, respectively. PEG samples yielded mean N1 and N2 gene concentrations of 1.61E+05 and 1.50E+05 gc/L, respectively. The overall observed trends of the VIRADEL total N1 and N2 gene concentrations increased steeply from early December 2021 and reached a peak in late December 2021 (Figure 6A), which heralded the major wave of COVID-19 cases in late December 2021 and early January 2022. Likewise, VIRADEL N1 and N2 gene concentrations increased in early April 2022, which preceded a resurgence of COVID-19 cases later in mid-April 2022.
Table 1. Total N1 and N2 gene concentrations measured in wastewater samples by VIRADEL, PEG, and filtration methods.
Figure 6. N1 and N2 gene concentrations (gc/L) by three concentration methods: (A) VIRADEL, (B) PEG, (C) Filtration, plotted against total COVID-19 cases.
Previous reports have demonstrated that the VIRADEL method can serve as a leading indicator of COVID-19 cases (1, 3, 13). By contrast, PEG measured N1 and N2 gene concentrations were more variable and increased significantly in January 2022, lagging the major wave of COVID-19 infections (Figure 6B). PEG N1 and N2 gene concentrations increased simultaneously with the surge of COVID-19 cases in mid-April 2022, into May 2022. N1 and N2 gene concentrations yielded by the filtration approach increased in early November 2021 and decreased in early December 2021. Thereafter, gene concentrations rapidly increased starting in mid-December 2021, peaking in mid-January 2022, which later significantly decreased to a low level in February 2022 (Figure 6C). Notably, the peak in SARS-CoV-2 measured in wastewater by this approach was staggered, lagging the major wave of COVID-19 cases.
3.2. Correlations and similarity analysis among three concentration methods
3.2.1. Correlations of N1 and N2 gene concentrations among three concentration methods
Multiple studies have applied Pearson and Spearman correlations to analyze the relationships between wastewater SARS-CoV-2 gene concentrations and COVID-19 cases (3, 9, 46), as well as the relationships among gene concentrations of SARS-CoV-2 in wastewater (47, 57). In this study, we tested the Pearson and Spearman correlations among N1 and N2 gene concentrations by VIRADEL, PEG, and filtration with normalized and non-normalized data (Table 2). A value of p that is less than 0.05 is considered statistically significant. For the non-normalized data (gc/L), the highest correlation was observed between PEG and filtration with N2 gene concentration (Pearson’s r = 0.67, Spearman’s r = 0.6). The lowest correlation was found between VIRADEL and PEG for N2 gene concentration (Pearson’s r = 0.12, Spearman’s r = 0.34). For non-normalized data (gc/L), the correlations between PEG and filtration were stronger than the correlations between VIRADEL and filtration, which in turn was stronger than the correlations between VIRADEL and PEG. For normalized data, the highest correlation was found between PEG and filtration for N1 (Pearson’s r = 0.73, Spearman’s r = 0.66) and N2 (Pearson’s r = 0.76, Spearman’s r = 0.64) gene concentrations in gc/[L(pounds/day)]. Significant correlations (Pearson coefficient > 0.63, Spearman coefficient > 0.6) were observed between PEG and filtration in gc/L, gc/mg TSS and gc/(L(pounds/day)) (Table 2). VIRADEL has stronger correlation to filtration than to PEG for both normalized and non-normalized data.
Table 2. Correlation coefficients among gene concentrations by VIRADEL, PEG, and filtration methods.
Normalizations using flow rate and flow rate/population reduced the correlations of gene concentrations among VIRADEL, PEG, and filtration compared to the correlations using the non-normalized data (gc/L) (Table 2). For instance, both Pearson and Spearman correlation coefficients between PEG and filtration were reduced from 0.67 (N2, Pearson, gc/L) and 0.6 (N2, Spearman, gc/L) to 0.45 (N2, Pearson, gc/day) and 0.5 (N2, Spearman, gc/day), respectively (Table 2). Conversely, normalizations using TSS and flow rate×TSS enhanced the correlations of gene concentrations among the three methods. For instance, higher correlation coefficients (Pearson’s r ranged from 0.73 (N1 gene) to 0.76 (N2 gene), Spearman’s r ranged from 0.64 (N2 gene) to 0.66 (N1 gene), all p < 0.05) were observed between PEG and filtration gene concentrations after normalization using flow rate×TSS compared to the correlation coefficients for non-normalized data (gc/L) (Pearson’s r ranged from 0.63 (N1 gene) to 0.67 (N2 gene), Spearman’s r = 0.6 (both N1 and N2 gene), all p < 0.05).
3.2.2. Dynamic time warping of N1 and N2 gene concentrations among three concentration methods
Detecting patterns and comparing similarities of gene concentration time series data are critical for comparing the concentration methods. Dynamic time warping (DTW) identifies the most similar patterns and the optimal warping match between two time series by calculating the minimum DTW distance (51, 53, 58). Shorter DTW distances indicate higher degree of similarity in patterns/shapes between two time series (59, 60). Table 3 presents the DTW minimum distances among the N1 and N2 gene concentrations by VIRADEL, PEG, and filtration methods. Smallest DTW distances were observed between VIRADEL and filtration for both non-normalized and normalized data, which indicated that VIRADEL has a higher degree of similarity with filtration than with PEG. Largest DTW distances were observed between VIRADEL and PEG for both non-normalized and normalized data, indicating that VIRADEL and PEG have the least similarity. This finding was consistent with the sampling and concentration mechanisms since VIRADEL targets free and suspended viral particles in the dissolved phase of wastewater, whereas PEG targets particle-associated viruses, some of which may represent previously shed and accumulated viruses in the sewer stream (3, 22).
Table 3. Dynamic time warping (DTW) minimum distances among gene concentrations by VIRADEL, PEG, and filtration methods.
Normalization using flow rate decreased the similarity among methods. For instance, the DTW distance between VIRADEL and filtration increased significantly after normalizing using flow rate (gc/day), indicating that the similarity between VIRADEL and filtration was reduced after normalization (Table 3). Conversely, normalization using TSS and flow rate×TSS strengthened the similarity among methods. For instance, the DTW distances decreased in gc/mg TSS and gc/(L(pounds/day)) comparing to the DTW distance in gc/L among the methods, indicating the improvement of similarity among methods after normalization (Table 3).
3.3. Similarity and TLCC analysis between three concentration methods and COVID-19 cases
3.3.1. Dynamic time warping between three concentration methods and COVID-19 cases
Wastewater surveillance data for COVID-19 primarily contain temporal data of viral gene concentrations and clinical cases. DTW analysis were performed between gene concentrations derived from the three concentration methods (VIRADEL, PEG, and filtration) and the 7-day moving average of total COVID-19 cases for both normalized and non-normalized data. For non-normalized data (gc/L), the smallest DTW distance was found between VIRADEL and total COVID-19 cases (Table 4). This indicates that VIRADEL (gc/L) has the highest similarity to total COVID-19 cases among the three concentration methods tested. The largest DTW distance was found between PEG (gc/L) and total COVID-19 cases, indicating the PEG method for concentration yields the least similarity to clinical cases. Normalizing gene concentration data using flow (gc/day) demonstrated similar findings. Conversely, normalization using TSS and flow×TSS can significantly increase the similarity between PEG and total COVID-19 cases but concurrently decrease the similarity between VIRADEL and total COVID-19 cases. Specifically, for normalized data (gc/mg TSS, gc/L(pounds/day)), the smallest DTW distance was identified between PEG and total COVID-19 cases, indicating the PEG has the highest similarity to total COVID-19 cases. The largest DTW distance was identified between VIRADEL and COVID-19 cases, indicating that VIRADEL has the lowest similarity to total COVID-19 cases.
Table 4. Dynamic time warping (DTW) minimum distances between gene concentrations by VIRADEL, PEG, as well as filtration methods and total COVID-19 cases.
3.3.2. Time lagged cross correlation and peak synchrony between three concentration methods and COVID-19 cases
The relative timing of the wastewater gene concentrations {gc/L, gc/day, gc/day/person, gc/mg TSS, and gc/[L(pounds/day)]} of VIRADEL, PEG and filtration were compared to the total COVID-19 cases using TLCC and peak synchrony. To evaluate if wastewater viral concentrations of the three methods lead or lag COVID-19 cases, the total COVID-19 case data were shifted by a period of −20 (lagging) to +20 days (leading) and the Pearson’s correlation coefficients were calculated between cases and wastewater viral gene concentration for each shift. The leading or lagging metric varied for each method, which was determined by comparing the strongest Pearson’s correlation coefficient.
For the VIRADEL method, both N1 and N2 gene concentrations (gc/L) were strongly correlated with COVID-19 cases, covering shifting windows between −20 and + 20 days (Figure 7A). The highest correlation coefficient was observed when offset is +12 days (Figure 7A), indicating that SARS-CoV-2 gene concentrations (gc/L) in wastewater by the VIRADEL method lead COVID-19 cases by 12 days, which concurred with previous findings of a 35-day lead time of gene concentrations preceding total COVID-19 cases prior to the Omicron surge (3). For both non-normalized and normalized data, VIRADEL always led COVID-19 cases with a variety of lead times (Table 5).
Figure 7. Pearson correlation coefficients for TLCC and peak synchrony between wastewater viral concentrations and COVID-19 cases with offsets between −20 (lagging) and + 20 (leading) days for the three methods, including (A) VIRADEL, (B) PEG, and (C) Filtration.
Table 5. Lead/lag time between wastewater viral concentrations by VIRADEL, PEG, as well as filtration methods and total COVID-19 cases.
For the PEG method (gc/L), the strongest correlation coefficients were observed with an offset of -12 days, indicating that SARS-CoV-2 gene concentrations by the PEG method lagged reported COVID-19 cases by 12 days during the study period (Figure 7B).
For the filtration method (gc/L), the highest correlation coefficient was observed with an offset of -7 days for the N1 gene and -11 days for the N2 gene, indicating that SARS-CoV-2 gene concentrations in wastewater lagged reported COVID-19 cases for 7 days (N1) and 11 days (N2), respectively (Figure 7C). Likewise, similar observations were found for normalized data where the filtration method yielded data that lagged clinical cases (Table 5). Table 5 summarized the lead/lag time between VIRADEL, PEG, and filtration methods and total COVID-19 cases. The length of the leading or lagging time differed with dissimilar normalizations. However, the leading or lagging pattern of each method did not change, where VIRADEL measurements were always leading COVID-19 cases, whereas PEG and filtration measurements routinely lagged COVID-19 cases.
4. Discussion
There is an ongoing effort to optimize methods to recover and concentrate SARS-CoV-2 from wastewater in support of actionable public health outcomes (33, 61). In this study, three concentration methods were evaluated for concentrating SARS-CoV-2 from wastewater, spanning the transition between Delta and Omicron variants circulating in the Detroit, MI metropolitan area. The three methods share common characteristics, especially downstream where they follow similar procedures of nucleic acid extraction and quantification such as RT-ddPCR or RT-qPCR. Likewise, their recovery efficiencies are reported as comparable (2, 3, 22, 34, 62).
4.1. VIRADEL: opportunities and obstacles
Several studies have previously adopted VIRADEL as the concentration method for SARS-CoV-2 in wastewater (1–4, 13). An attribute of the VIRADEL method is the ability to process large volumes (10 – 50 L) of wastewater, thus facilitating capture of free and suspended viral particles that are arguably most representative of viruses shed by recently infected individuals (3, 35). This establishes VIRADEL as a concentration method capable to provide early warning that leads case reporting (1, 3), which was also verified by TLCC analysis in this study (Table 5). Limiting widescale adoption of VIRADEL is labor-intensive preparation of sampling units which require extensive washing and disinfection prior to use. VIRADEL (63) also requires access to large volumes of wastewater which may not be available to all researchers. Further, the required large volumes may necessitate use of grab samples which typically yield higher variability than composite samples which is the sampling method of choice for many wastewater surveillance efforts (64). VIRADEL requires trained personnel for comparatively laborious work with limited samples (n = 15) processed over a relatively long time (4-6 h). VIRADEL also requires multiple large centrifuges as well as expensive and at times, supply chain-limited consumables. Therefore, VIRADEL may not be an ideal choice for routine wastewater surveillance projects in common microbiology laboratories. However, it was clear from the comparative analysis conducted that VIRADEL has clear potential to be implemented as a tool to provide early warning to inform public health actions (1, 3, 13).
4.2. PEG: opportunities and obstacles
Apart from requiring access to a centrifuge, the consumables required are widely available and relatively inexpensive, lending itself as one of the most broadly applied concentration methods for routine wastewater surveillance (3, 22, 31, 33, 62, 63). On the other hand, PEG is restricted to processing smaller volumes of wastewater (usually 0.05 to 2 L) and only a portion of the sample pellet is used to recover and extract RNA, which can be affected by the variation of samples and representation of all viruses in wastewater (3, 22, 33, 35).
Unlike VIRADEL, PEG targets particle-associated viruses consistent with reports that identify solids as the phase yielding highest SARS-CoV-2 concentrations in wastewater (63). While a fraction of these particles will represent recently deposited SARS-CoV-2, the majority may represent previously shed and accumulated viruses in the sewer stream and later resuspended during flow fluctuations, thus providing a mechanism for the method to yield data lagging clinical COVID-19 cases. Though the exact mechanism of PEG is not well understood, several studies proposed that it captures viruses that are sorbed to larger precipitates and solids, consistent with a high quantity of TSS in wastewater (3, 22). In this study, through the DTW analysis, PEG yielded data were normalized using TSS and flow × TSS, which increased the degree of similarity between PEG and total COVID-19 cases (Table 4). This demonstrated that PEG yielded data were largely affected by the presence of TSS. VIRADEL, instead, captured free and suspended viruses in the supernatant wastewater. Thus, normalizing the VIRADEL data using TSS and flow × TSS decreased the similarity between VIRADEL and cases (Table 4).
Through the TLCC analysis, this study also demonstrated that PEG gene concentrations lagged COVID-19 cases (Table 5), which embraced the aforementioned sampling mechanism of PEG (22). PEG method did not provide an early warning (leading window) for COVID-19 cases which concurred with our previous findings, whereas VIRADEL provided early warnings ahead of clinical cases while PEG lagged clinical cases for the Detroit area (3).
However, several studies using PEG provided early warnings of impending COVID-19 cases (65). Notably, in these studies, PEG was applied to different types of samples such as primary sewage sludge, which is a different sample matrix from untreated wastewater samples, thus needing more investigation on the impact of sample types on early warnings (65). Kumar et al. (66) identified early warnings using PEG in the early stage of the pandemic in August 2020 in India (66). PEG and other concentration methods [such as ultrafiltration (17, 67) and adsorption-precipitation (68)] identified early warnings in the early stage of the pandemic when testing capacities were largely limited, and societal responses to the pandemic and clinical data reporting were significantly delayed (3, 69). In addition, earlier prevalent COVID-19 variants including Alpha, Beta and Gamma were reported with longer incubation times than Delta and Omicron variants, leading to prolonged early warning potentials of wastewater surveillance in the early stages of the pandemic (70).
Though PEG was reported to provide early warnings, it may have a shorter early warning window than VIRADEL due to the fundamental disparity of their targets, that being newly contributed free and suspended viral particles versus particle-attached virus, some of which may be considered previously shed and accumulated and subsequently resuspended from sediment (3). In the current study, PEG was shown to lag clinical cases while VIRADEL was leading clinical cases for both normalized and non-normalized data (Table 5). Overall, the early warning potential of PEG needs further investigations in terms of sample types, sampling mechanisms and locations, stage of the epidemic, among other factors.
4.3. Filtration: opportunities and obstacles
Filtration is commonly applied to recover and concentrate viral RNA in water samples. It achieves generally good recovery efficiencies, is relatively inexpensive using commonly available lab equipment and simple protocols and provides consistent performance and inclusive measurement since it captures viruses from both solids and liquid fractions by nature of forcing free viral particles across trapped solids (33, 37). However, filtration has several drawbacks. First, the number of available filtration units restricts the number of samples that can be processed concurrently (33). Meanwhile, clogging of filters can occur due to high variations of turbidity in wastewater. While this can be offset in part by use of a caulking gun to exert more pressure on the sample being filtered, in reality, volumes are limited to ~0.1 L. Additionally, filtration measurements lagged the COVID-19 clinical cases in the current study, thus, its ability to provide early warnings for impending cases is called into question. The recovery efficiencies also differ with different filters (33).
4.4. Future directions
The mechanism and implications of primarily collecting viruses attached to solids that may have settled and resuspended before sampling, such as by the PEG, needs further investigations. Notably, multiple studies have reported that the integrity of SARS-CoV-2 RNA was higher when sorbed to suspended solids, organic matter, and large bio-solids which provide protection from predation and inactivation. This can be explained by the hydrophobicity of SARS-CoV-2 viral particles, leading to their adherence to wastewater solids and longer persistence compared to free viruses in the supernatant wastewater (71–73).
The implications of seasonal variations in SARS-CoV-2 persistence in wastewater needs further investigations. Seasonal variations of wastewater temperature and pH are reported to affect the persistence of viral RNA in wastewater (74). However, SARS-CoV-2 RNA was shown to be highly stable at 4°C aqueous environment or in a wide pH range at room temperature (75, 76). Multiple studies reported the detectability and persistence of SARS-CoV-2 RNA in untreated wastewater solids samples. For instance, researchers found that SARS-CoV-2 RNA was consistently detected for 29 days and 64 days at 4°C and -20°C, respectively in wastewater solids pelleted by centrifugation (77). Another study indicated that only minimal reduction of SARS-CoV-2 RNA was observed for wastewater solids samples after 100 days (78). Additionally, researchers established models to indicate that viral RNA can be detected in wastewater even with long sewer travel time (100 h), especially with lower average wastewater temperature in northern cities such as Detroit (74). A recent study also indicated that biofilms could mediate the fate of SARS-CoV-2 in wastewater, especially leading the viral RNA to prolonged presence (79).
The effect of varying sampling volumes needs further investigation. Some studies indicated that a larger sampling volume can increase the sensitivity of the sampling method, suggesting that it will detect lower levels of viral RNA in wastewater samples (80). Similarly, researchers suggested that processing of larger sample volumes may help to lower the method detection limits (74). But at the same time, keeping the required samples sizes low can reduce shipping costs between sampling location and the analytical laboratory as well as reduce space for storage (74). Other researchers indicated that detection sensitivity can be improved by increasing the sample volume from 100 mL to 500 mL wastewater for testing SARS-CoV-2 (6).
However, other researchers presented that large-volume sampling did not significantly enhance the sensitivity of methods (81). For instance, Zheng et al. (81) found that wastewater concentration methods (they used ultracentrifugation) using less volume of wastewater was preferable than larger volume of wastewater in terms of sensitivity for testing SARS-CoV-2. The study revealed that when using the same concentration methods, no significant difference was observed in the viral RNA concentrations between experiments conducted with a larger volume of wastewater and those conducted with a smaller volume (81). Some studies indicated that a larger sampling volume may also dilute the wastewater sample, which can lead to a lower viral RNA concentration (82).
Overall, the sampling volume for wastewater surveillance of SARS-CoV-2 using different concentration methods will depend on several factors, including the sensitivity of the method, the concentration of viral RNA in the wastewater, and the size of the population being monitored. It is critical to consider and address these factors when analyzing wastewater surveillance data and more in-depth research on how the sampling volume affect statistical results are needed.
The time of sampling may potentially affect results in sewershed sampling. The effect of sampling time in large interceptors, like the ones sampled in this study, is less significant, since the interceptor wastewater is mixed at the pumping stations. A few studies have reported gene concentration varying on an hourly basis (83, 84) although the temporal variability of SARS-CoV-2 concentrations in wastewater remains ambiguous (83, 85). It has been suggested that composite samples may circumvent the within-day variation of viral concentrations (83). Whereas both the PEG and filtration methods used composite samples, the large volume required for VIRADEL necessitated separate sampling which was conducted over a period of several hours to help reduce temporal variability. Further, considering the vast sewersheds and population of nearly 4 million people that GLWA’s three interceptors serve, the concentrations of SARS-CoV-2 in wastewater may be highly diluted and within-day variations can be negligible. Future studies are called to examine within-day variation of SARS-CoV-2.
Admittedly, there are caveats to the current study that should be considered and discussed. The study period was limited to the transition between Delta and Omicron VOCs that occurred between fall 2021 and winter 2022. With each successive resurgence of COVID-19, differences are reported related to disease trajectory including incubation time, shedding dynamics and disease severity (18, 86). For instance, the incubation time was shorter during the Omicron surge compared to the previous variants, inevitably reducing the early warning potentials of wastewater surveillance in the later stage of the pandemic (3, 86). Further, the changing viral shedding dynamics, viral decay kinetics, and shedding duration of the Omicron variant are not well understood and many uncertainties remain (18, 87). As such, the lead and lag times reported here cannot be extrapolated to past or future SARS-CoV-2 variants. In addition, sampling frequency was limited to weekly samples and thus less informative for establishing time series or less likely to depict accurately the actual fluctuations of wastewater viral concentrations (cdc.gov). Feng et al. (88) proposed a minimum of two samples collected weekly to establish the time series data of wastewater viral concentrations for continuous trend analysis. Some researchers have even suggested daily or very frequent sampling, if the laboratory is capable of handling increased numbers of samples, considering rapid resurgence of COVID-19 cases (89). Indeed, the filtration method has been used to analyze samples 5 days weekly since the emergence of the Omicron VOC as part of Ontario’s Wastewater Surveillance Initiative in the Windsor-Essex region located across the international border with Detroit (Q. Geng, R. Corchis-Scott, R.M. McKay, unpublished). While SARS-CoV-2 signal intensity derived from this approach does not provide a clear early warning of clinical cases, preliminary analysis supports its use as a leading indicator of COVID-19-related hospitalizations in the region (Q. Geng, R. Corchis-Scott, R.M. McKay, unpublished). This is important considering that clinical testing capacity across North America was overwhelmed by infections attributed to Omicron and is thus no longer a reliable indicator of disease prevalence (90).
5. Conclusion
This study is among the first to implement, evaluate, and compare commonly applied wastewater virus concentration methodologies to recover and concentrate SARS-CoV-2 from wastewater amid the transition between Delta and Omicron VOCs. Analytical approaches, including Pearson and Spearman correlations, Dynamic Time Warping (DTW), and Time Lagged Cross Correlation (TLCC) and peak synchrony, were performed to analyze the relations among three methods as well as the relations between methods and COVID-19 cases. To our knowledge, this is the only study to implement Dynamic Time Warping to compare wastewater surveillance time series data and successfully identify the similarities/dissimilarities among the methods and between methods and clinical data. The analytical approach used can be applied to different sample processing and concentration methods under various pandemic scenarios to evaluate method efficacy for different public health practices.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
LZ: methodology, investigation, data acquisition, data curation, formal analysis, visualization, writing–original draft, and writing–review and editing. QG and RC-S: methodology, investigation, data acquisition, and writing–review and editing. RM: methodology, investigation, data acquisition, resources, supervision, funding acquisition, and writing–review and editing. JN: funding acquisition and writing–review and editing. IX: conceptualization, funding acquisition, methodology, investigation, project administration, resources, supervision, and writing–review and editing.
Funding
This study was funded by the Michigan Department of Health and Human Services (MDHHS) and the Great Lakes Water Authority (GLWA). RMM acknowledges support from Ontario Genomics (COVID-19 Regional Genomics Initiative) with additional support provided by the Canada Foundation for Innovation–Exceptional Opportunities Fund (COVID-19 program) and National Institutes of Health [1P01ES028939-01], and National Science Foundation [OCE-1840715] awards to the Bowling Green State University Great Lakes Center for Fresh Waters and Human Health.
Acknowledgments
The authors acknowledge the support from operators, the laboratory team and managers of the Wastewater Resource Recovery Facility at the Great Lakes Water Authority. They also thank CDM Smith for supporting this work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
1. Miyani, B, Zhao, L, Spooner, M, Buch, S, Gentry, Z, Mehrotra, A, et al. Early warnings of COVID-19 second wave in Detroit. J Environ Eng (New York). (2021) 147:6021004. doi: 10.1061/(ASCE)EE.1943-7870.0001907
2. Li, Y, Miyani, B, Zhao, L, Spooner, M, Gentry, Z, Zou, Y, et al. Surveillance of SARS-CoV-2 in nine neighborhood sewersheds in Detroit Tri-County area, United States: assessing per capita SARS-CoV-2 estimations and COVID-19 incidence. Sci Total Environ. (2022) 851:158350. doi: 10.1016/j.scitotenv.2022.158350
3. Zhao, L, Zou, Y, Li, Y, Miyani, B, Spooner, M, Gentry, Z, et al. Five-week warning of COVID-19 peaks prior to the omicron surge in Detroit, Michigan using wastewater surveillance. Sci Total Environ. (2022) 844:157040:157040. doi: 10.1016/j.scitotenv.2022.157040
4. Miyani, B, Fonoll, X, Norton, J, Mehrotra, A, and Xagoraraki, I. SARS-CoV-2 in Detroit Wastewater. J Environ Eng (New York). (2020) 146:06020004. doi: 10.1061/(ASCE)EE.1943-7870.0001830
5. Sherchan, SP, Shahin, S, Ward, LM, Tandukar, S, Aw, TG, Schmitz, B, et al. First detection of SARS-CoV-2 RNA in wastewater in North America: a study in Louisiana, USA. Sci Total Environ. (2020) 743:140621. doi: 10.1016/j.scitotenv.2020.140621
6. Ahmed, W, Bertsch, PM, Angel, N, Bibby, K, Bivins, A, Dierens, L, et al. Detection of SARS-CoV-2 RNA in commercial passenger aircraft and cruise ship wastewater: a surveillance tool for assessing the presence of COVID-19 infected travellers. J Travel Med. (2020) 27:116. doi: 10.1093/jtm/taaa116
7. Ahmed, W, Bivins, A, Simpson, SL, Bertsch, PM, Ehret, J, Hosegood, I, et al. Wastewater surveillance demonstrates high predictive value for COVID-19 infection on board repatriation flights to Australia. Environ Int. (2022) 158:106938. doi: 10.1016/j.envint.2021.106938
8. Kaya, D, Niemeier, D, Ahmed, W, and Kjellerup, BV. Evaluation of multiple analytical methods for SARS-CoV-2 surveillance in wastewater samples. Sci Total Environ. (2022) 808:152033. doi: 10.1016/j.scitotenv.2021.152033
9. Ai, Y, Davis, A, Jones, D, Lemeshow, S, Tu, H, He, F, et al. Wastewater SARS-CoV-2 monitoring as a community-level COVID-19 trend tracker and variants in Ohio, United States. Sci Total Environ. (2021) 801:149757. doi: 10.1016/j.scitotenv.2021.149757
10. Xiao, A, Wu, F, Bushman, M, Zhang, J, Imakaev, M, Chai, PR, et al. Metrics to relate COVID-19 wastewater data to clinical testing dynamics. Water Res. (2022) 212:118070. doi: 10.1016/j.watres.2022.118070
11. Bivins, A, and Bibby, K. Wastewater surveillance during mass COVID-19 vaccination on a college campus. Environ Sci Technol Lett. (2021) 8:792–8. doi: 10.1021/acs.estlett.1c00519
12. Chik, AHS, Glier, MB, Servos, M, Mangat, CS, Pang, XL, Qiu, Y, et al. Comparison of approaches to quantify SARS-CoV-2 in wastewater using RT-qPCR: results and implications from a collaborative inter-laboratory study in Canada. J Environ Sci. (2021) 107:218–29. doi: 10.1016/j.jes.2021.01.029
13. Zhao, L, Zou, Y, David, RE, Withington, S, McFarlane, S, Faust, RA, et al. Simple methods for early warnings of COVID-19 surges: lessons learned from 21 months of wastewater and clinical data collection in Detroit, Michigan, United States. Sci Total Environ. (2022):161152. doi: 10.1016/j.scitotenv.2022.161152
14. Jones, DL, Baluja, MQ, Graham, DW, Corbishley, A, McDonald, JE, Malham, SK, et al. 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. Sci Total Environ. (2020) 749:141364. doi: 10.1016/j.scitotenv.2020.141364
15. He, X, Lau, EHY, Wu, P, Deng, X, Wang, J, Hao, X, et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat Med. (2020) 26:672–5. doi: 10.1038/s41591-020-0869-5
16. Haramoto, E, Malla, B, Thakali, O, and Kitajima, M. First environmental surveillance for the presence of SARS-CoV-2 RNA in wastewater and river water in Japan. Sci Total Environ. (2020) 737:140405. doi: 10.1016/j.scitotenv.2020.140405
17. Medema, G, Heijnen, L, Elsinga, G, Italiaander, R, and Brouwer, A. 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. Environ Sci Technol Lett. (2020) 7:511–6. doi: 10.1021/acs.estlett.0c00357
18. Boucau, J, Marino, C, Regan, J, Uddin, R, Choudhary, MC, Flynn, JP, et al. Duration of shedding of Culturable virus in SARS-CoV-2 omicron (BA.1) infection. N Engl J Med. (2022) 387:275–7. doi: 10.1056/NEJMc2202092
19. Farkas, K, Hillary, LS, Malham, SK, McDonald, JE, and Jones, DL. Wastewater and public health: the potential of wastewater surveillance for monitoring COVID-19. Curr Opin Environ Sci Health. (2020) 17:14–20. doi: 10.1016/j.coesh.2020.06.001
20. Xagoraraki, I, Yin, Z, and Svambayev, Z. Fate of viruses in water systems. J Environ Eng. (2014) 140:04014020. doi: 10.1061/(ASCE)EE.1943-7870.0000827
21. Xagoraraki, I. Can we predict viral outbreaks using wastewater surveillance? J Environ Eng. (2020) 146:01820003. doi: 10.1061/(ASCE)EE.1943-7870.0001831
22. Flood, MT, D’Souza, N, Rose, JB, and Aw, TG. Methods evaluation for rapid concentration and quantification of SARS-CoV-2 in raw wastewater using droplet digital and quantitative RT-PCR. Food Environ Virol. (2021) 13:303–15. doi: 10.1007/s12560-021-09488-8
23. Wallis, C, Henderson, M, and Melnick, JL. Enterovirus Concentration on Cellulose Membranes. Appl Microbiol. (1972) 23:476–80. doi: 10.1128/am.23.3.476-480.1972
24. Jakubowski, W, Hill, WF, and Clarke, NA. Comparative study of four microporous filters for concentrating viruses from drinking water. Appl Microbiol. (1975) 30:58–65. doi: 10.1128/am.30.1.58-65.1975
25. Lee, H, Kim, M, Paik, SY, Lee, CH, Jheong, WH, Kim, J, et al. Evaluation of electropositive filtration for recovering norovirus in water. J Water Health. (2011) 9:27–36. doi: 10.2166/wh.2010.190
26. Hill, VR, Polaczyk, AL, Kahler, AM, Cromeans, TL, Hahn, D, and Amburgey, JE. Comparison of hollow-Fiber ultrafiltration to the USEPA VIRADEL technique and USEPA method 1623. JEQ. (2009) 38:822–5. doi: 10.2134/jeq2008.0152
27. Black, LE, Brion, GM, and Freitas, SJ. Multivariate logistic regression for predicting Total Culturable virus presence at the intake of a potable-water treatment plant: novel application of the atypical coliform/Total coliform ratio. Appl Environ Microbiol. (2007) 73:3965–74. doi: 10.1128/AEM.02780-06
28. Francy, DS, Stelzer, EA, Brady, AMG, Huitger, C, Bushon, RN, Ip, HS, et al. Comparison of filters for concentrating microbial indicators and pathogens in Lake water samples. Appl Environ Microbiol. (2013) 79:1342–52. doi: 10.1128/AEM.03117-12
29. Ye, Y, Ellenberg, RM, Graham, KE, and Wigginton, KR. Survivability, partitioning, and recovery of enveloped viruses in untreated municipal wastewater. Environ Sci Technol. (2016) 50:5077–85. doi: 10.1021/acs.est.6b00876
30. Michael-Kordatou, I, Karaolia, P, and Fatta-Kassinos, D. Sewage analysis as a tool for the COVID-19 pandemic response and management: the urgent need for optimised protocols for SARS-CoV-2 detection and quantification. J Environ Chem Eng. (2020) 8:104306. doi: 10.1016/j.jece.2020.104306
31. Sapula, SA, Whittall, JJ, Pandopulos, AJ, Gerber, C, and Venter, H. An optimized and robust PEG precipitation method for detection of SARS-CoV-2 in wastewater. Sci Total Environ. (2021) 785:147270. doi: 10.1016/j.scitotenv.2021.147270
32. Barril, PA, Pianciola, LA, Mazzeo, M, Ousset, MJ, Jaureguiberry, MV, Alessandrello, M, et al. Evaluation of viral concentration methods for SARS-CoV-2 recovery from wastewaters. Sci Total Environ. (2021) 756:144105. doi: 10.1016/j.scitotenv.2020.144105
33. Ahmed, W, Bertsch, PM, Bivins, A, Bibby, K, Farkas, K, Gathercole, A, et al. Comparison of virus concentration methods for the RT-qPCR-based recovery of murine hepatitis virus, a surrogate for SARS-CoV-2 from untreated wastewater. Sci Total Environ. (2020) 739:139960. doi: 10.1016/j.scitotenv.2020.139960
34. Pecson, BM, Darby, E, Haas, CN, Amha, YM, Bartolo, M, Danielson, R, et al. Reproducibility and sensitivity of 36 methods to quantify the SARS-CoV-2 genetic signal in raw wastewater: findings from an interlaboratory methods evaluation in the U.S. Environ Sci (Camb). (2021) 7:504–20. doi: 10.1039/d0ew00946f
35. Lu, D, Huang, Z, Luo, J, Zhang, X, and Sha, S. Primary concentration–the critical step in implementing the wastewater based epidemiology for the COVID-19 pandemic: a mini-review. Sci Total Environ. (2020) 747:141245. doi: 10.1016/j.scitotenv.2020.141245
36. McKindles, KM, Manes, MA, DeMarco, JR, McClure, A, McKay, RM, Davis, TW, et al. Dissolved microcystin release coincident with lysis of a bloom dominated by Microcystis spp. in Western Lake Erie attributed to a novel Cyanophage. Appl Environ Microbiol. (2020) 86:20. doi: 10.1128/AEM.01397-20
37. Gonçalves, J, Koritnik, T, Mioč, V, Trkov, M, Bolješič, M, Berginc, N, et al. Detection of SARS-CoV-2 RNA in hospital wastewater from a low COVID-19 disease prevalence area. Sci Total Environ. (2021) 755:143226. doi: 10.1016/j.scitotenv.2020.143226
38. Corchis-Scott, R, Geng, Q, Seth, R, Ray, R, Beg, M, Biswas, N, et al. Averting an outbreak of SARS-CoV-2 in a university residence hall through wastewater surveillance. Microbiol Spectr. (2021) 9:21. doi: 10.1128/Spectrum.00792-21
39. Kantor, RS, Nelson, KL, Greenwald, HD, and Kennedy, LC. Challenges in measuring the recovery of SARS-CoV-2 from wastewater. Environ Sci Technol. (2021) 55:3514–9. Available from:. doi: 10.1021/acs.est.0c08210
40. Lu, X, Wang, L, Sakthivel, SK, Whitaker, B, Murray, J, Kamili, S, et al. US CDC real-time reverse transcription PCR panel for detection of severe acute respiratory syndrome coronavirus 2. Emerg Infect Dis. (2020) 26:1654–65. doi: 10.3201/eid2608.201246
41. Menkir, TF, Chin, T, Hay, JA, Surface, ED, de Salazar, PM, Buckee, CO, et al. Estimating internationally imported cases during the early COVID-19 pandemic. Nat Commun. (2021) 12:311. doi: 10.1038/s41467-020-20219-8
42. Barua, VB, Juel, MAI, Blackwood, AD, Clerkin, T, Ciesielski, M, Sorinolu, AJ, et al. Tracking the temporal variation of COVID-19 surges through wastewater-based epidemiology during the peak of the pandemic: a six-month long study in Charlotte, North Carolina. Sci Total Environ. (2022) 814:152503:152503. doi: 10.1016/j.scitotenv.2021.152503
43. Lepot, M, Aubin, JB, and Clemens, F. Interpolation in time series: an introductive overview of existing methods, their performance criteria and uncertainty assessment. Water (Basel). (2017) 9:796. doi: 10.3390/w9100796
44. Hopkins, L, Persse, D, Caton, K, Ensor, K, Schneider, R, McCall, C, et al. Citywide wastewater SARS-CoV-2 levels strongly correlated with multiple disease surveillance indicators and outcomes over three COVID-19 waves. Sci Total Environ. (2023) 855:158967. doi: 10.1016/j.scitotenv.2022.158967
45. Terry, IB, and Meschke, JS. Variability in normalization methods of COVID-19 wastewater surveillance. ProQuest dissertations and theses, Ann Arbor. (2022). Available at: http://ezproxy.msu.edu/login?url=https://www.proquest.com/dissertations-theses/variability-normalization-methods-covid-19/docview/2695037740/se-2?accountid=12598.
46. D'Aoust, PM, Mercier, E, Montpetit, D, Jia, JJ, Alexandrov, I, Neault, N, et al. Quantitative analysis of SARS-CoV-2 RNA from wastewater solids in communities with low COVID-19 incidence and prevalence. Water Res X. (2021) 188:116560. doi: 10.1016/j.watres.2020.116560
47. Vadde, KK, Al-Duroobi, H, Phan, DC, Jafarzadeh, A, Moghadam, SV, Matta, A, et al. Assessment of concentration, recovery, and normalization of SARS-CoV-2 RNA from two wastewater treatment plants in Texas and correlation with COVID-19 cases in the community. ACS ES&T Water. (2022) 2022:54. doi: 10.1021/acsestwater.2c00054
48. Ye, J, Xiao, C, Esteves, RM, and Rong, C. Time series similarity evaluation based on Spearman’s correlation coefficients and distance measures, pp. 319–331. (2015).
49. Keogh, E, and Ratanamahatana, CA. Exact indexing of dynamic time warping. Knowl Inf Syst. (2005) 7:358–86. doi: 10.1007/s10115-004-0154-9
50. Izakian, H, Pedrycz, W, and Jamal, I. Fuzzy clustering of time series data using dynamic time warping distance. Eng Appl Artif Intell. (2015) 39:235–44. doi: 10.1016/j.engappai.2014.12.015
51. Giorgino, T. Computing and visualizing dynamic time warping alignments in R: the dtw package. J Stat Softw. (2009) 31:1–24.
52. Rakthanmanon, T, Campana, B, Mueen, A, Batista, G, Westover, B, Zhu, Q, et al. Searching and mining trillions of time series subsequences under dynamic time warping. In: Proceedings of the 18th ACM SIGKDD international conference on knowledge discovery and data mining - KDD’12. New York, USA: ACM Press, p. 262. (2012).
53. Jeong, YS, Jeong, MK, and Omitaomu, OA. Weighted dynamic time warping for time series classification. Pattern Recogn. (2011) 44:2231–40. doi: 10.1016/j.patcog.2010.09.022
54. Shen, C. Analysis of detrended time-lagged cross-correlation between two nonstationary time series. Phys Lett A. (2015) 379:680–7. doi: 10.1016/j.physleta.2014.12.036
55. Mei, DC, Du, LC, and Wang, CJ. The effects of time delay on stochastic resonance in a Bistable system with correlated noises. J Stat Phys. (2009) 137:625–38. doi: 10.1007/s10955-009-9864-4
56. Li, H, Futch, SH, and Syvertsen, JP. Cross-correlation patterns of air and soil temperatures, rainfall and Diaprepes abbreviatus root weevil in citrus. Pest Manag Sci. (2007) 63:1116–23. doi: 10.1002/ps.1431
57. Lancaster, E, Byrd, K, Ai, Y, and Lee, J. Socioeconomic status correlations with confirmed COVID-19 cases and SARS-CoV-2 wastewater concentrations in small-medium sized communities. Environ Res. (2022) 215:114290. doi: 10.1016/j.envres.2022.114290
58. Berndt, DJ, and Clifford, J. Using dynamic time warping to find patterns in time series. KDD workshop. Seattle, WA, pp. 359–370. (1994).
59. Liu, L, Li, W, and Jia, H. Method of time series similarity measurement based on dynamic time warping. Comput Mater Contin. (2018) 57:97–106. doi: 10.32604/cmc.2018.03511
60. Guan, X, Huang, C, Liu, G, Meng, X, and Liu, Q. Mapping Rice cropping Systems in Vietnam Using an NDVI-based time-series similarity measurement based on DTW distance. Remote Sens. (2016) 8:19. doi: 10.3390/rs8010019
61. Kitajima, M, Ahmed, W, Bibby, K, Carducci, A, Gerba, CP, Hamilton, KA, et al. SARS-CoV-2 in wastewater: state of the knowledge and research needs. Sci Total Environ. (2020) 2020:139076. doi: 10.1016/j.scitotenv.2020.139076
62. Hata, A, Hara-Yamamura, H, Meuchi, Y, Imai, S, and Honda, R. Detection of SARS-CoV-2 in wastewater in Japan during a COVID-19 outbreak. Sci Total Environ. (2021) 758:143578. doi: 10.1016/j.scitotenv.2020.143578
63. Torii, S, Furumai, H, and Katayama, H. Applicability of polyethylene glycol precipitation followed by acid guanidinium thiocyanate-phenol-chloroform extraction for the detection of SARS-CoV-2 RNA from municipal wastewater. Sci Total Environ. (2021) 756:143067. doi: 10.1016/j.scitotenv.2020.143067
64. Bivins, A, Kaya, D, Ahmed, W, Brown, J, Butler, C, Greaves, J, et al. Passive sampling to scale wastewater surveillance of infectious disease: lessons learned from COVID-19. Sci Total Environ. (2022) 835:155347. doi: 10.1016/j.scitotenv.2022.155347
65. D’Aoust, PM, Graber, TE, Mercier, E, Montpetit, D, Alexandrov, I, Neault, N, et al. Catching a resurgence: increase in SARS-CoV-2 viral RNA identified in wastewater 48 h before COVID-19 clinical tests and 96 h before hospitalizations. Sci Total Environ. (2021) 770:145319. doi: 10.1016/j.scitotenv.2021.145319
66. Kumar, M, Joshi, M, Patel, AK, and Joshi, CG. Unravelling the early warning capability of wastewater surveillance for COVID-19: a temporal study on SARS-CoV-2 RNA detection and need for the escalation. Environ Res. (2021) 196:110946. doi: 10.1016/j.envres.2021.110946
67. Hasan, SW, Ibrahim, Y, Daou, M, Kannout, H, Jan, N, Lopes, A, et al. Detection and quantification of SARS-CoV-2 RNA in wastewater and treated effluents: surveillance of COVID-19 epidemic in the United Arab Emirates. Sci Total Environ. (2021) 764:142929. doi: 10.1016/j.scitotenv.2020.142929
68. Randazzo, W, Truchado, P, Cuevas-Ferrando, E, Simón, P, Allende, A, and Sánchez, G. SARS-CoV-2 RNA in wastewater anticipated COVID-19 occurrence in a low prevalence area. Water Res. (2020) 181:115942. doi: 10.1016/j.watres.2020.115942
69. Bibby, K, Bivins, A, Wu, Z, and North, D. Making waves: plausible lead time for wastewater based epidemiology as an early warning system for COVID-19. Water Res. (2021) 202:117438. doi: 10.1016/j.watres.2021.117438
70. Wu, Y, Kang, L, Guo, Z, Liu, J, Liu, M, and Liang, W. Incubation period of COVID-19 caused by unique SARS-CoV-2 strains. JAMA Netw Open. (2022) 5:e2228008. doi: 10.1001/jamanetworkopen.2022.28008
71. Abu Ali, H, Yaniv, K, Bar-Zeev, E, Chaudhury, S, Shagan, M, Lakkakula, S, et al. Tracking SARS-CoV-2 RNA through the wastewater treatment process. ACS ES&T Water. (2021) 1:1161–7. doi: 10.1021/acsestwater.0c00216
72. Gundy, PM, Gerba, CP, and Pepper, IL. Survival of coronaviruses in water and wastewater. Food Environ Virol. (2009) 1:10. doi: 10.1007/s12560-008-9001-6
73. Panchal, D, Prakash, O, Bobde, P, and Pal, S. SARS-CoV-2: sewage surveillance as an early warning system and challenges in developing countries. Environ Sci Pollut Res. (2021) 28:22221–40. doi: 10.1007/s11356-021-13170-8
74. Hart, OE, and Halden, RU. Computational analysis of SARS-CoV-2/COVID-19 surveillance by wastewater-based epidemiology locally and globally: feasibility, economy, opportunities and challenges. Sci Total Environ. (2020) 730:138875. doi: 10.1016/j.scitotenv.2020.138875
75. Ahmed, W, Bertsch, PM, Bibby, K, Haramoto, E, Hewitt, J, Huygens, F, et al. Decay of SARS-CoV-2 and surrogate murine hepatitis virus RNA in untreated wastewater to inform application in wastewater-based epidemiology. Environ Res. (2020) 191:110092. doi: 10.1016/j.envres.2020.110092
76. Chin, AWH, Chu, JTS, Perera, MRA, Hui, KPY, Yen, HL, Chan, MCW, et al. Stability of SARS-CoV-2 in different environmental conditions. Lancet Microb. (2020) 1:e10. doi: 10.1016/S2666-5247(20)30003-3
77. Hokajärvi, AM, Rytkönen, A, Tiwari, A, Kauppinen, A, Oikarinen, S, Lehto, KM, et al. The detection and stability of the SARS-CoV-2 RNA biomarkers in wastewater influent in Helsinki, Finland. Sci Total Environ. (2021) 770:145274. doi: 10.1016/j.scitotenv.2021.145274
78. Simpson, A, Topol, A, White, BJ, Wolfe, MK, Wigginton, KR, and Boehm, AB. Effect of storage conditions on SARS-CoV-2 RNA quantification in wastewater solids. PeerJ. (2021) 9:e11933. doi: 10.7717/peerj.11933
79. Li, J, Ahmed, W, Metcalfe, S, Smith, WJM, Choi, PM, Jackson, G, et al. Impact of sewer biofilms on fate of SARS-CoV-2 RNA and wastewater surveillance. Nat Water. (2023) 1:272–80. doi: 10.1038/s44221-023-00033-4
80. McMinn, BR, Korajkic, A, Kelleher, J, Herrmann, MP, Pemberton, AC, Ahmed, W, et al. Development of a large volume concentration method for recovery of coronavirus from wastewater. Sci Total Environ. (2021) 774:145727. doi: 10.1016/j.scitotenv.2021.145727
81. Zheng, X, Deng, Y, Xu, X, Li, S, Zhang, Y, Ding, J, et al. Comparison of virus concentration methods and RNA extraction methods for SARS-CoV-2 wastewater surveillance. Sci Total Environ. (2022) 824:153687. doi: 10.1016/j.scitotenv.2022.153687
82. Bertels, X, Demeyer, P, van den Bogaert, S, Boogaerts, T, van Nuijs, ALN, Delputte, P, et al. Factors influencing SARS-CoV-2 RNA concentrations in wastewater up to the sampling stage: a systematic review. Sci Total Environ. (2022) 820:153290. doi: 10.1016/j.scitotenv.2022.153290
83. Bivins, A, North, D, Wu, Z, Shaffer, M, Ahmed, W, and Bibby, K. Within-and between-day variability of SARS-CoV-2 RNA in municipal wastewater during periods of varying COVID-19 prevalence and positivity. Acs Es&T Water. (2021) 1:2097–108. doi: 10.1021/acsestwater.1c00178
84. Li, B, Di, DYW, Saingam, P, Jeon, MK, and Yan, T. Fine-scale temporal dynamics of SARS-CoV-2 RNA abundance in wastewater during a COVID-19 lockdown. Water Res. (2021) 197:117093. doi: 10.1016/j.watres.2021.117093
85. Ahmed, W, Bivins, A, Bertsch, PM, Bibby, K, Gyawali, P, Sherchan, SP, et al. Intraday variability of indicator and pathogenic viruses in 1-h and 24-h composite wastewater samples: implications for wastewater-based epidemiology. Environ Res. (2021) 193:110531. doi: 10.1016/j.envres.2020.110531
86. Baker, MA, Rhee, C, Tucker, R, Badwaik, A, Coughlin, C, Holtzman, MA, et al. Rapid control of hospital-based severe acute respiratory syndrome coronavirus 2 omicron clusters through daily testing and universal use of N95 respirators. Rev Infect Dis. (2022) 75:e296–9. doi: 10.1093/cid/ciac113
87. Kandel, C, Lee, Y, Taylor, M, Llanes, A, McCready, J, Crowl, G, et al. Viral dynamics of the SARS-CoV-2 omicron variant among household contacts with 2 or 3 COVID-19 vaccine doses. J Infect. (2022) 85:666–70. doi: 10.1016/j.jinf.2022.10.027
88. Feng, S, Roguet, A, McClary-Gutierrez, JS, Newton, RJ, Kloczko, N, Meiman, JG, et al. Evaluation of sampling, analysis, and normalization methods for SARS-CoV-2 concentrations in wastewater to assess COVID-19 burdens in Wisconsin communities. ACS ES&T Water. (2021) 1:1955–65. doi: 10.1021/acsestwater.1c00160
89. Zhu, Y, Oishi, W, Maruo, C, Saito, M, Chen, R, Kitajima, M, et al. Early warning of COVID-19 via wastewater-based epidemiology: potential and bottlenecks. Sci Total Environ. (2021) 767:145124. doi: 10.1016/j.scitotenv.2021.145124
Keywords: wastewater surveillance, SARS-CoV-2, COVID-19, virus adsorption-elution, polyethylene glycol precipitation, filtration, lead/lag time, dynamic time warping
Citation: Zhao L, Geng Q, Corchis-Scott R, McKay RM, Norton J and Xagoraraki I (2023) Targeting a free viral fraction enhances the early alert potential of wastewater surveillance for SARS-CoV-2: a methods comparison spanning the transition between delta and omicron variants in a large urban center. Front. Public Health. 11:1140441. doi: 10.3389/fpubh.2023.1140441
Edited by:
David Champredon, Public Health Agency of Canada (PHAC), CanadaReviewed by:
Niti B. Jadeja, University of Virginia, United StatesHyatt Green, SUNY College of Environmental Science and Forestry, United States
Madhvi Joshi, Gujarat Biotechnology Research Centre (GBRC), India
Copyright © 2023 Zhao, Geng, Corchis-Scott, McKay, Norton and Xagoraraki. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Irene Xagoraraki, eGFnb3JhcmFAbXN1LmVkdQ==