Abstract
We report a time course of SARS-CoV-2 RNA concentrations in primary sewage sludge during the Spring COVID-19 outbreak in a northeastern U.S. metropolitan area. SARS-CoV-2 RNA was detected in all environmental samples and, when adjusted for the time lag, the virus RNA concentrations were highly correlated with the COVID-19 epidemiological curve (R2=0.99) and local hospital admissions (R2=0.99). SARS-CoV-2 RNA concentrations were a seven-day leading indicator ahead of compiled COVID-19 testing data and led local hospital admissions data by three days. Decisions to implement or relax public health measures and restrictions require timely information on outbreak dynamics in a community.
Introduction
The most common metric followed to track the progression of the COVID-19 epidemic within communities is derived from testing symptomatic cases and evaluating the number of positive tests over time.1 However, tracking positive tests is a lagging indicator for the epidemic progression.2, 3 Testing is largely prompted by symptoms, which may take up to five days to present4, and individuals can shed virus prior to exhibiting symptoms. There is a pressing need for additional methods for early sentinel surveillance and real-time estimations of community disease burden so that public health authorities may modulate and plan epidemic responses accordingly.
SARS-CoV-2 RNA is present in the stool of COVID-19 patients5-7 and has recently been documented in raw wastewater.8-10 Thus, monitoring raw wastewater (sewage) within a community’s collection system can potentially provide information on the prevalence and dynamics of infection for entire populations.11 When municipal raw wastewater discharges into treatment facilities, solids are settled and collected into a matrix called (primary) sewage sludge, which has been shown to contain a broad diversity of human viruses including commonly circulating coronavirus strains.12 Primary sludge provides a well-mixed and concentrated sample that may be advantageous for monitoring SARS-CoV-2. As viral shedding can occur before cases are detected, we hypothesize that the time course of SARS-CoV-2 RNA concentrations in primary sewage sludge is a leading indicator of outbreak dynamics within a community served by the treatment plant.
Results
During the COVID-19 outbreak from March 19, 2020 to May 1, 2020 in the New Haven, Connecticut (CT), USA metropolitan area (Figure 1A), daily primary sludge samples were collected from the wastewater treatment facility which serves approximately 200,000 residents. SARS-CoV-2 viral RNA concentrations were quantitatively compared with local hospital admission data and community COVID-19 compiled testing data. SARS-CoV-2 viral RNA was detectable in all samples tested and ranged from 1.7 × 103 virus RNA copies mL-1 to 4.6 × 105 virus RNA copies mL-1. The lower concentration in this range corresponds to a qRT-PCR cycle threshold (CT) value of 38.75 and can be considered a detection threshold for this method and sludge matrix. Overall, 96.5% of all CT values were less than 38 and values between 38 and 40 were reported as positive if detection occurred with virus nucleocapsid N1 and N2 primer sets and both replicates. Replicated samples demonstrated similar SARS-CoV-2 RNA concentration values and comparisons between SARS-CoV-2 primer values with the human ribonuclease P (RP) gene primer values indicate that dynamic temporal concentration changes in SARS-CoV-2 RNA were from actual change in virus concentration (Figure 1B). Concentration comparisons between replicates produced an average slope of 0.98 ± 0.17 st. dev., while replicates comparing N and RP primer values result in a flat line, with average slope of 0.05 ± 0.05 st. dev. (Figure 1B). Overall, the SARS-CoV-2 RNA concentration time series results consistently reveal a curve that increases over times, peaks, and then decreases to May 1 (Figure 1C).
(A) Epidemiology curve for COVID-19 new cases (based on testing) from the cities in the wastewater catchment area (New Haven, East Haven, and Hamden, CT) served by the East Shore Water Pollution Abatement Facility (ESWPAF)1; (B) Example comparison of SARS-CoV-2 virus RNA replicates (N1 primer set) and comparison between N1 primer concentrations versus human RP concentrations; (C) SARS-CoV-2 virus RNA concentration time course for all replicates and primers considered.
New Haven COVID-19 epidemic suggest that these data may provide useful epidemiological insights (Figure 1C). SARS-CoV-2 RNA sludge concentrations were quantitatively compared with data that are commonly used to track the community progression of COVID-19 including hospital admissions (Figure 2A,B) and COVID-19 compiled testing data for the four municipalities (New Haven, East Haven, Hamden, and Woodbridge, CT) served by the ESWPAF (Figure 2C). All three measures traced the initial wave of the SARS-CoV-2 outbreak in the New Haven metropolitan area. Applying Locally Weighted Scatterplot Smoothing (LOWESS) to the data and rescaling enables comparison (Figure 2A,B). The virus RNA time course peaked 3 days earlier than hospital admissions (April 9 versus April 12) and a cross correlation analysis revealed a correlation coefficient (R=0.996) between smoothed RNA and hospital data when the latter was shifted 3 days forward.
(A) Average sludge SARS-CoV-2 RNA concentration time course data (o) and average hospital admissions data (o) with LOWESS smoothing; (B) rescaled smoothed SARS-CoV-2 virus RNA concentrations (---) and hospital admissions (---); (C) smoothed sludge SARS-CoV-2 virus RNA concentration (---) with smoothed COVID-19 epidemiology curve (---); (D) regression between smoothed virus RNA and new COVID-19 cases (ascending), slope=1,240 virus RNA copies/new case, R2=0.99; (E) regression between smoothed virus RNA and new COVID-19 cases (descending), slope=1,305 virus RNA copies/new case, R2=0.97.
Comparing the LOWESS-smoothed daily virus RNA and COVID-19 testing data reveals that the peak in COVID-19 cases (April 16) occurred 7 days following the peak in SARS-CoV-2 RNA concentration (April 9) (Figure 2C). A cross correlation analysis resulted in a maximum correlation coefficient (R=0.994) when the COVID-19 testing data is adjusted 7 days forward. Regressing the adjusted COVID-19 testing data and SARS-CoV-2 virus RNA concentrations yield a value of 1,305 virus RNA copies/mL per new COVID-19 case report as the outbreak ascended (Figure 2D), and 1,240 virus RNA copies/mL per new COVID-19 case report as the outbreak descended (Figure 2E).
Discussion
We produced a SARS-CoV-2 RNA concentration time course in primary sewage sludge during a COVID-19 outbreak in the New Haven, CT metropolitan area. Approximately 200,000 people are served by the treatment facility and COVID-19 total documented cases by testing rose from less than 29 to 2,609 during the March 19 to May 1, 2020 surveillance period. Our results demonstrate: (1) the utility of SARS-CoV-2 primary sludge monitoring to accurately track outbreaks in a community and (2) primary sludge SARS-CoV-2 RNA concentrations can be a leading indicator over other commonly used epidemiology approaches including summarized COVID-19 test results and hospital admissions.
This study uniquely utilized primary sewage sludge instead of raw wastewater for virus RNA measurements. Sewage sludge is comprised of solids that are removed during primary sedimentation steps and typically gravity thickened. As a result, primary sludge has a greater solids content (2 to 5%) than raw wastewater (0.01 to 0.05%). Due to the elevated solids content and the high case load observed during the outbreak (~1,200 per 100,000 population), the concentrations of SARS-CoV-2 RNA reported here ranged from two to three orders of magnitude greater than raw wastewater SARS-CoV-2 values previously reported8, 10, 13. Monitoring primary sludge is broadly applicable. Wastewater treatment plants with primary and secondary treatment are standard in many regions of the world, and treatment facilities are rapidly expanding in urban areas of lower and middle-income countries.14 Within the US, approximately 16,000 treatment plants serve more than 250,000,000 people.15
SARS-CoV-2 RNA concentrations in sewage sludge were a leading indicator of community outbreak dynamics over hospitalization and compiled COVID-19 testing data. SARS-CoV-2 RNA concentrations led hospital admissions by 3 days and COVID-19 cases by 7 days. Hospital admissions to Yale New Haven Hospital from the four towns served by the wastewater treatment facility both rose and fell more slowly than the observed RNA concentrations. The rates of SARS-CoV-2 RNA increases and decreases were similar to the COVID-19 testing data, and reflect the very concentrated outbreak that took place within New Haven and its neighboring towns. The strong correlation between the virus RNA and COVID-19 case data provides an approach to independently evaluate the local SARS-CoV-2 testing rate and to estimate the number of new cases. While the results here demonstrate the potential utility of such an approach, we note that these relationships could be treatment plant specific, as primary sludge process trains are not uniform, and that compiled COVID-19 testing is prompted by symptoms, thus underestimating the true number of cases.16
In conclusion, the SARS-CoV-2 RNA concentrations in primary sewage sludge were tracked throughout a COVID-19 epidemic and compared with traditional outbreak epidemiological indicators. SARS-CoV-2 RNA concentration in primary sludge closely followed the epidemiology curves established by compiled COVID-19 testing data and hospital admissions, but was a leading indicator by seven and three days, respectively. Our study could have substantial policy implications. Jurisdictions can use primary sludge SARS-CoV-2 concentrations to preempt community outbreak dynamics or provide an additional basis for easing restrictions, especially when there are limitations in clinical testing. Raw wastewater and sludge-based surveillance is particularly useful for low and middle-income countries where clinical testing capacity is limited.
Methods
Sample collection
Primary sewage sludge was collected from the East Shore Water Pollution Abatement Facility (ESWPAF) in New Haven, CT, USA. Samples were taken daily from March 19 to May 1, 2020 between 8 and 10 am EDT and stored at −80°C prior to analysis. The first sampling dates were prior to widespread individual testing in the region, and prior to the start of stay at home restrictions implemented throughout the State of Connecticut, USA (Figure 1A). From the sampling start and end dates, cities served by the ESWPAF experienced greater than 85 times increase in confirmed COVID-19 cases (by testing) from less than 29 cases to 2,609.17 The plant serves an estimated population of 200,000 people with average treated flows of 1.75 m3/s. Sludge collected from ESWPAF is primary sludge, sampled at the outlet of a gravity thickener, ranging in solids content from 2.6% to 5%.
Viral RNA quantitative testing
To quantify SARS-CoV-2 virus RNA concentrations in primary sludge, 2.5 mL of well mixed sludge were added directly to a commercial kit optimized for isolation of total RNA from soil (RNeasey PowerSoil Total RNA kit, Qiagen). Isolated RNA pellets were dissolved in 50 μL of ribonuclease free water and total RNA measured by spectrophotometry (NanoDrop, ThermoFisher Scientific). SARS-CoV-2 RNA was quantified by one-step quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) using the U.S. Center for Disease Control N1 and N2 primers sets.18 For control, analysis was also conducted for the human Ribonuclease P (RP) gene. Samples were analyzed using the Bio-Rad iTaq Universal Probes One-Step kit in 20 μL reactions run at 50°C for 10 min, 95°C for 1 min followed by 40 cycles of 95°C for 10 s and 60°C for 30 s per the manufacturer’s recommendations. SARS-CoV-2 RNA concentrations were determined using a standard curve as previously described18 and presented as virus RNA copies. The SARS-CoV-2 concentration results were normalized to the total RNA extracted to reduce day to day variations in treatment plant flow, sludge solids content, and RNA extraction efficiency. All samples were diluted 5x for use as template to ensure the removal of inhibition. Performing qRT-PCR on undiluted and 5x diluted samples typically resulted in the expected 5x decrease in concentration. Sewage sludge from 2018 was used as a control and no SARS-CoV-2 detection was observed from either N1 or N2 primers. These control sludges were consistently positive for the human RP gene.
Epidemiological Analysis
Daily COVID-19 admissions to the Yale New Haven Hospital were compiled from hospital records and confirmed by laboratory testing. Numbers of laboratory-confirmed COVID-19 cases in the towns served by the ESWPAF (New Haven, East Haven, Hamden, and Woodbridge CT) were compiled from daily reports published by the CT Department of Public Health.17 SARS-CoV-2 RNA, COVID-19 cases, and hospital admissions time series were smoothed using Locally Weighted Scatterplot Smoothing (LOWESS) to better visualize the time trend in the data.19
Data Availability
Once published, data will be available on the COVID tracker website.
Acknowledgements
We wish to thank the East Shore Water Pollution Abatement Facility, New Haven, CT, USA for primary sludge sampling assistance. AZ is supported by the Yale Institute of Global Health and NDG is supported by a gift from the Huffman Family Donor Advised Fund.
References
