Within-Day Variability of SARS-CoV-2 RNA in Municipal Wastewater Influent During Periods of Varying COVID-19 Prevalence and Positivity
Summary (3 min read)
INTRODUCTION
- When infected with severe acute coronavirus 2 (SARS-CoV-2), the β-coronavirus which causes coronavirus disease 2019 (COVID-19), humans, both symptomatic and asymptomatic 1 , shed the virus and its RNA, in various body fluids [2] [3] [4] including: sputum, saliva, urine, and feces.
- Surveillance strategies and sampling methods for WBE remain diverse, with community-level temporal trends monitored via both primary solids 7, 8 and primary influent 9 .
- Concentrations of SARS-CoV-2 RNA in wastewater and wastewater solids correlate with COVID-19 cases 8, 17, 18 and positivity rates 19 .
- WBE studies of wastewater solids have reported improved correlation with clinical case trends 26 and no effect 8 associated with PMMoVnormalization, while a study of wastewater influent found that PMMoV-normalization decreased correlation with clinical case trends 18 . .
(which was not certified by peer review)
- The copyright holder for this preprint this version posted March 24, 2021.
- The first goal was to characterize the within-day variation in influent SARS-CoV-2 RNA concentrations, PMMoV RNA concentrations, and process control recovery efficiency.
- The copyright holder for this preprint this version posted March 24, 2021.
- During a methods comparison, significantly different recovery efficiencies were observed between a variety of process controls for a single method 42 .
- BRSV RNA also demonstrated a similar increase over seven days at 4°C .
Primary Influent Sampling Locations
- The experiments described herein were conducted at two WWTPs located in two communities, identified as community A and community B, in northern Indiana, USA.
- Records from the Environmental Protection Agency's (EPA) Enforcement and Compliance History Online (ECHO) system indicate the design flow for each WWTP is 20 million gallons per day (MGD During each experiment, 500 mL primary influent grab samples were collected at 2-hour intervals and immediately stored at 4°C.
- Samples were then transported on ice to the laboratory and again stored at 4°C until concentrated as described below within 24 hours.
24-h Sampling Experiments
- While at WWTP B, a 24-hour time-based composite sample was only prepared using the grab samples from 1 December to 2 December.
- The average hourly flow rates were recorded during each experiment and subsequently used to calculate average flow rates for each 2-hour interval and for the entire 24-hour experiment.
Direct Extractions
- In addition to the primary influent samples concentrated by the adsorption-extraction method, a paired subset of 16 samples, eight collected from WWTP A and 8 from WWTP B (December 2020 experiments), were extracted by adding 500 µL of influent directly into a Garnet PowerBead tube and extracting the nucleic acids as described above.
- The purpose of these direct extractions was to directly estimate the virus RNA concentration recovery efficiency by comparing the direct extraction enumerations and the adsorption-extraction enumerations.
RT-ddPCR
- RNA in sample extracts was detected and quantified by RT-ddPCR performed on the BioRad QX200 Droplet Digital PCR System with thermal cycling performed on the C1000 Touch Thermal Cycler (BioRad, Hercules, CA, USA).
- RNA reverse transcription and PCR amplification was performed in a single reaction using the One-Step RT-ddPCR Advanced Kit for Probes (BioRad, Hercules, CA, USA) per the manufacturer's instructions.
- Primer and probe sequences, concentrations, and thermal cycling conditions for each RT-ddPCR assay are summarized in Table S1 .
- Each RT-ddPCR experiment included notemplate controls, positive controls, and the pertinent negative extraction controls as described in further detail below.
- The RNA copy number for each RT-ddPCR reaction was estimated by manual thresholding performed in QuantaSoft Version 1.7.4 (BioRad, Hercules, CA, USA) such that the negative controls, both no-template and extraction, were negative for each assay.
RT-ddPCR Assays
- SARS-CoV-2 RNA was detected and quantified using the CDC N1 assay targeting the nucleocapsid gene 32 .
- The N1 copy number in each sample was measured in triplicate RT-ddPCR reactions using the premixed primers and probe (Table S1 ) from the 2019-nCoV RUO Kit (IDT, Coralville, IA, USA).
- At each step in the dilution series, 12 RT-ddPCR replicates were assayed.
- A cumulative Gaussian distribution was fit to the observed proportion of positive technical replicates along the dilution series, and the 95% LOD was estimated as the 95th percentile of the resulting distribution.
RNA Persistence Experiments
- In addition to the 24-h influent sampling experiments, a daily composite sample was collected from WWTP A on 23 June 2020, seeded with BRSV, and used to investigate the stability of RNA during storage, pasteurization, and freeze-thaw cycles.
- To assess persistence during storage, the composite sample was aliquoted into 50 mL centrifuge tubes.
- The tubes were incubated at either 4°C or 25°C with two tubes combined into a single 100 mL sample and processed every 24 hours from time zero to seven days.
- Persistence through pasteurization was assessed by pasteurizing two 50 mL aliquots in centrifuge tubes at 60°C for 90 minutes, with a brief vortex mix at 45 minutes, and then combining the two aliquots into a single 100 mL sample.
- After each thaw, two tubes were combined into one 100 mL sample for processing.
Process Control, Molecular Control & Concentration Recovery Efficiency
- Across all experiments, 83 primary influent samples were concentrated by adsorption-extraction and assayed for SARS-CoV-2 and PMMoV RNA: 58 from WWTP A and 25 from WWTP B (one 24-hour event did not include a composite sample).
- The observed recovery efficiency in samples from WWTP A was greater than in samples from WWTP B (p <0.0001); however, the coefficient of variation (CV) in samples from WWTP A was also greater (169%) than WWTP B (83%).
- Interestingly, the mean recovery for wastewater seeded with Hep G was greater than for PCRgrade water seeded with Hep G (37%, n = 2).
Within-Day Variation in Primary Influent
- During both 24-hour sampling intervals at WWTP A , hourly flow rates peaked from roughly 9:00 to midnight.
- Summary statistics for each parameter measured during the 24-hour sampling events are listed in Table S2 .
- PMMoV concentrations in primary influent at WWTP A were comparable between sampling events (p >0.9999), while PMMoV concentrations were greater during the May than December sampling in primary influent at WWTP B (p <0.0001).
- During the May sampling event higher recoveries were observed during periods of both high and low flow.
WWTP Influent SARS-CoV-2 RNA and Clinical Surveillance
- Due to the agreement between BRSV and SARS-CoV-2 RNA recovery, the effect of recoverycorrection on within-day trends, and the large variance associated with recovery efficiency, recovery-corrected SARS-CoV-2 RNA concentrations in primary influent were compared to county-level COVID-19 cases and positivity rates during the two weeks prior to each 24-hour sampling period.
- SARS-CoV-2 RNA concentrations did not consistently increase with increasing average daily COVID-19 cases.
- The results of the current study indicate that when sub-county level COVID-19 clinical surveillance data are not available, positivity may offer a better metric for comparison with wastewater data.
- Differences in the primary influent concentrations during each of these periods were often not meaningfully different after accounting for variation in concentration and recovery.
- These observations indicate that quantitative relationships between wastewater data and SARS-CoV-2 infection prevalence, particularly those premised on material balance, are likely to remain constrained by variability and uncertainty 21 .
Did you find this useful? Give us your feedback
Citations
116 citations
36 citations
29 citations
29 citations
29 citations
References
11,055 citations
9,365 citations
3,943 citations
3,305 citations
1,813 citations
Related Papers (5)
Frequently Asked Questions (10)
Q2. How many recoveries of BRSV were observed during a study in Virginia?
electronegative membrane filtration methods have demonstrated 252 recoveries ranging from approximately 10% to less than 1% during virus concentration from 253 sewage42.
Q3. What is the way to compare COVID-19 with clinical surveillance data?
The results of the current study indicate that when 437 sub-county level COVID-19 clinical surveillance data are not available, positivity may offer a 438 better metric for comparison with wastewater data.
Q4. What is the significance of the process control in wastewater surveillance?
The large variation in recovery efficiency observed 376 during 24-hour sampling periods draws further attention to the importance the consistent use of 377 process controls in wastewater surveillance despite their limitations
Q5. How many efficiencies were recovered using the adsorption-extraction method?
In the current study SARS-CoV-2 and PMMoV were recovered at different mean efficiencies 243 (1.9% and 19.2%, respectively) using the adsorption-extraction method.
Q6. What were the effects of lockdowns on the flow patterns of the primary influent?
flow patterns during all the sampling events 455 were likely affected by changes in human behavior patterns associated with lockdowns and 456 interrupted domestic and working routines.
Q7. What was the recovery adjustment for the composite samples?
Even with recovery adjustment, composite samples still yielded lower 365 SARS-CoV-2 RNA concentrations than the 24-hour average from grab samples.
Q8. What is the effect of the process control on the recovery efficiency of SARS-CoV-2?
The sporadic use of process or molecular 382 controls observed in the WBE literature greatly limits the ability to compare SARS-CoV-2 RNA 383 measurements within and between WWTPs56.
Q9. What is the correlation between the adsorption of SARS-CoV-2 to surfaces?
In a study of 262 SARS-CoV-2 adsorption to surfaces in solution, electrostatic adhesion correlated with both 263 solution ionic strength and surface chemistry49.
Q10. How was the recovery efficiency of the BRSV 218?
For a subset of 217 four samples where the solids fraction was removed prior to adsorption-extraction, the BRSV 218 mean recovery efficiency was 25% (95%CI: 22 – 28).