Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Editor's Pick Research Article | Applied and Environmental Science

SARS-CoV-2 Titers in Wastewater Are Higher than Expected from Clinically Confirmed Cases

Fuqing Wu, Jianbo Zhang, Amy Xiao, Xiaoqiong Gu, Wei Lin Lee, Federica Armas, Kathryn Kauffman, William Hanage, Mariana Matus, Newsha Ghaeli, Noriko Endo, Claire Duvallet, Mathilde Poyet, Katya Moniz, Alex D. Washburne, Timothy B. Erickson, Peter R. Chai, Janelle Thompson, Eric J. Alm
Jack A. Gilbert, Editor
Fuqing Wu
aCenter for Microbiome Informatics and Therapeutics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
nDepartment of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
oDepartment of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Fuqing Wu
Jianbo Zhang
aCenter for Microbiome Informatics and Therapeutics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
nDepartment of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Amy Xiao
aCenter for Microbiome Informatics and Therapeutics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
nDepartment of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
oDepartment of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiaoqiong Gu
bSingapore-MIT Alliance for Research and Technology, National University of Singapore, Singapore
mCampus for Research Excellence and Technological Enterprise, Singapore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wei Lin Lee
bSingapore-MIT Alliance for Research and Technology, National University of Singapore, Singapore
mCampus for Research Excellence and Technological Enterprise, Singapore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Federica Armas
bSingapore-MIT Alliance for Research and Technology, National University of Singapore, Singapore
mCampus for Research Excellence and Technological Enterprise, Singapore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kathryn Kauffman
cUniversity at Buffalo, The State University of New York, Buffalo, New York, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
William Hanage
dCenter for Communicable Disease Dynamics, Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mariana Matus
eBiobot Analytics, Inc., Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Newsha Ghaeli
eBiobot Analytics, Inc., Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Noriko Endo
eBiobot Analytics, Inc., Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Claire Duvallet
eBiobot Analytics, Inc., Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mathilde Poyet
aCenter for Microbiome Informatics and Therapeutics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
nDepartment of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
oDepartment of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Katya Moniz
aCenter for Microbiome Informatics and Therapeutics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
nDepartment of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
oDepartment of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alex D. Washburne
pSelva Analytics, LLC, Bozeman, Montana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Timothy B. Erickson
fDivision of Medical Toxicology, Department of Emergency Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, USA
gHarvard Humanitarian Institute, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter R. Chai
fDivision of Medical Toxicology, Department of Emergency Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, USA
hThe Fenway Institute, Boston, Massachusetts, USA
iThe Koch Institute for Integrated Cancer Research, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Janelle Thompson
jSingapore Center for Environmental Life Sciences Engineering, Nanyang Technological University, Singapore
kAsian School of the Environment, Nanyang Technological University, Singapore
mCampus for Research Excellence and Technological Enterprise, Singapore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eric J. Alm
aCenter for Microbiome Informatics and Therapeutics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
bSingapore-MIT Alliance for Research and Technology, National University of Singapore, Singapore
eBiobot Analytics, Inc., Cambridge, Massachusetts, USA
lBroad Institute of MIT and Harvard, Cambridge, Massachusetts, USA
mCampus for Research Excellence and Technological Enterprise, Singapore
nDepartment of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
oDepartment of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Eric J. Alm
Jack A. Gilbert
University of California San Diego
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00614-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Wastewater surveillance represents a complementary approach to clinical surveillance to measure the presence and prevalence of emerging infectious diseases like the novel coronavirus SARS-CoV-2. This innovative data source can improve the precision of epidemiological modeling to understand the penetrance of SARS-CoV-2 in specific vulnerable communities. Here, we tested wastewater collected at a major urban treatment facility in Massachusetts and detected SARS-CoV-2 RNA from the N gene at significant titers (57 to 303 copies per ml of sewage) in the period from 18 to 25 March 2020 using RT-qPCR. We validated detection of SARS-CoV-2 by Sanger sequencing the PCR product from the S gene. Viral titers observed were significantly higher than expected based on clinically confirmed cases in Massachusetts as of 25 March. Our approach is scalable and may be useful in modeling the SARS-CoV-2 pandemic and future outbreaks.

IMPORTANCE Wastewater-based surveillance is a promising approach for proactive outbreak monitoring. SARS-CoV-2 is shed in stool early in the clinical course and infects a large asymptomatic population, making it an ideal target for wastewater-based monitoring. In this study, we develop a laboratory protocol to quantify viral titers in raw sewage via qPCR analysis and validate results with sequencing analysis. Our results suggest that the number of positive cases estimated from wastewater viral titers is orders of magnitude greater than the number of confirmed clinical cases and therefore may significantly impact efforts to understand the case fatality rate and progression of disease. These data may help inform decisions surrounding the advancement or scale-back of social distancing and quarantine efforts based on dynamic wastewater catchment-level estimations of prevalence.

INTRODUCTION

Improved understanding of the presence and prevalence of SARS-CoV-2 at a population level can help government, public health, and hospital officials implement appropriate policies to mitigate the exponential spread of COVID-19 and diminish the future strain on health care facilities. Despite pandemic spread of SARS-CoV-2 worldwide, broad access to testing in the United States (US) has thus far been severely limited. While it is impractical to test every US resident for SARS-CoV-2, the virus has been found in the stool of confirmed COVID-19 patients (1), making it a promising candidate for wastewater-based surveillance (WBS).

WBS can help detect the presence of pathogens across municipalities and estimate population prevalence without invasive individual testing, and inform public health officials of emerging local hot spots, permitting empirical deployment of testing centers and the efficacy of interventions. The closely related virus SARS-CoV-1 was detected in wastewater from Chinese hospitals during the 2002–2003 SARS pandemic (2), and WBS has been used for the early detection and direct mitigation of disease outbreaks in Israel, Egypt, and Sweden (3–7). We have previously used this technique to measure and map the use of pharmaceuticals across residential communities (8). Here, we describe an analytical technique to detect and quantify genetic material from SARS-CoV-2 in wastewater collected at a large treatment facility.

RESULTS

We first collected two sewage samples on 18 March 2020 at a major urban wastewater treatment facility in Massachusetts, which has two major influent streams (here named “Northern” and “Southern” influents). Samples were transported to our laboratory where we conducted sample inactivation and viral enrichment (with 80-ml samples), nucleic acid extraction, and RT-qPCR. As negative controls, we used two biobanked wastewater samples (collected on January 8 and 11) from the same treatment facility taken before the first US case was documented.

Initial testing with PCR and gel electrophoresis using primers specific for the SARS-CoV-2 S gene (9) indicated that both samples (from Southern and Northern influents) collected from the treatment facility on 18 March had a positive signal for SARS-CoV-2. We confirmed the signal by identifying a PCR product at 137 bp (see Fig. S1A in the supplemental material). A no-template control for the PCR assays was negative. Sanger sequencing of the PCR products confirmed a 97 to 98% identity match to the SARS-CoV-2 S gene (Fig. 1 and Fig. S1B).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Sanger sequencing results and alignment to SARS-CoV-2 S gene. Sequence was aligned to SARS-CoV-2 S gene (GenBank MT577641.1) using NCBI BLASTN. Highlighted gray, forward and reverse primers; highlighted yellow, aligned sequence; highlighted purple, mismatched sequence from either primer. The first 24 bases in the COVID-19 S gene are the T7 promoter added by PCR with S-F and S-R primers in Table 1. Results for the Southern influent sample were similar. Sequencing chromatograms for forward and reverse primers are provided in Fig. S1B.

FIG S1

PCR products with SARS-CoV-2 S gene primers and Sanger sequencing results. (A) cDNA was used as the template and run for 42 cycles. Lane 1, DNA ladder (1kb Plus); lane 2, negative control (no template); lane 3, 18 March sample (Southern influent); lane 4, 18 March sample (Northern influent); lane 5, sample on 8 January (combined from Southern and Northern influents); lane 6, sample on 11 January (combined from Southern and Northern influents); lane 7, 50-bp DNA ladder (first three ladders were labeled); lane 8, positive control using cDNA reverse transcribed from synthetic SARS-CoV-2 RNA; lane 9, DNA ladder (1kb Plus). (B) Sanger sequencing results of the PCR product. (Top) Full sequence of the PCR product of S gene. Gray nucleotides are the primer sequences. (Middle and bottom) Sequencing chromatograms for forward and reverse primers, respectively. Download FIG S1, EPS file, 2.9 MB.
Copyright © 2020 Wu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

We next sought to establish our viral enrichment protocol and quantify viral titer in sewage using RT-qPCR. We used the US CDC primer/probe sets (10) targeting the N1, N2, and N3 loci of the SARS-CoV-2 nucleocapsid gene to amplify cDNA reverse transcribed from viral RNA (Table 1). Standard curves for all three primer sets showed linear behavior over 6 log units (r2 ranging from 0.981 to 0.998) using DNA standards (SARS-CoV-2 nucleocapsid gene) (Fig. S2). We tested different steps in our viral enrichment process. First, we examined raw (unfiltered) sewage precipitated with polyethylene glycol 8000 (PEG), which recovers both bacterial and viral nucleic acids. Next, we looked at samples taken from 0.2-μm-filtered sewage: we considered both the material collected on the filter and the filtrate. We found the strongest and most consistent results from the PEG-precipitated viral pellet from the 0.2-μm filtrate, which was resuspended in TRIzol for RNA extraction (Table 2).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Primers and probes used in this studya

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Mean CT values for Southern and Northern samples (18 March 2020) with and without prefiltration (filtrate versus unfiltered) and solids fraction (filter)

FIG S2

Standard curves for N1, N2, and N3 primers. Serial dilutions (10-fold) of a plasmid containing SARS-CoV-2 nucleocapsid gene (N) gene were made to generate the standard curves. Value shown is averaged from three replicates. Download FIG S2, TIF file, 5.3 MB.
Copyright © 2020 Wu et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

These results suggest a simple viral enrichment and RNA extraction protocol is sufficient to achieve viral identification. Importantly, we included pasteurization (90 min at 60°C) as a first step performed before sample containers were opened, increasing the safety of the protocol. Prior work on SARS-CoV-1 indicates that a 30-min heat inactivation at 60°C is sufficient to inactivate the virus by over 6 log units (11). We increased the pasteurization time as an extra precaution for a novel virus and to account for the time required to bring sewage samples to temperature.

After demonstrating the ability of our assay to detect the presence of SARS-CoV-2 in wastewater treatment facility samples, we next quantified the viral titers of SARS-CoV-2 in the wastewater based on the standard curves for each primer set (Fig. S2). Sewage samples from two catchment areas (Southern and Northern) were tested at five sampling dates from 18 to 25 March. All the 10 samples tested were positive for SARS-CoV-2 by RT-qPCR (threshold cycle [CT] < 40 cycles) (Fig. 2). The mean SARS-CoV-2 titer estimated from all the three primer sets ranges from 57 to 303 copies per ml of sewage. In contrast, both of the biobanked samples collected from before the first known US SARS-CoV-2 case were negative for all three primer sets (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Estimated viral titer per milliliter of sewage. Estimated gene copies in sewage samples for each of the three CDC primer sets are shown for all dates. Northern (N) and Southern (S) influents are shown separately, except for biobanked January samples, which are combined.

Next, we examined how quickly viral signal degraded during storage at 4°C. All samples were received during 20 to 25 March, processed upon initial receipt, and then stored at 4°C. We reprocessed those samples starting from pasteurized sewage on 4 April. Figure 3 shows that there was variation between sample runs on different dates, but no clear trend toward lower signals at the later time point, suggesting samples can be stored at 4°C for more than a week (9 to 15 days) without significant degradation of viral RNA in pasteurized sewage samples.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Sample stability over time. Estimated copies per milliliter of sewage estimated from the same samples filtered and processed on different dates. The y axis shows samples processed between 20 and 25 March 2020, and these are plotted against samples reprocessed on 4 April 2020 (samples were pasteurized upon receipt and stored at 4°C until processing). The gray dashed line corresponds to y = x. Each point represents an average over all three primer sets for each sample.

Viral titers in sewage samples are determined by the concentrations of fecal materials in the total flow at the wastewater treatment plant. However, the sewage flow rate is not stable and is impacted by many factors including inflows from rainstorms. This factor is prominent especially during this dynamic winter-spring season in Massachusetts. To correct those variations, we use pepper mild mottle virus (PMMoV) as an internal reference for quantification of SARS-CoV-2. Previous studies have shown that PMMoV is the most abundant RNA virus in human feces and it is shed in large quantities in wastewater across the United States (12–16). PMMoV is remarkably stable in the wastewater, and its concentrations showed little seasonal variation (15, 16). Furthermore, like SARS-CoV-2, PMMoV is also a positive-sense single-stranded RNA virus, making it suitable as an internal reference to help control for sample-to-sample variability in wastewater dilution or processing. These properties allow us to calibrate the SARS-CoV-2 titers across the samples.

The qPCR results showed that PMMoV levels varied across the 12 samples, with a standard deviation of 1.34 CT units. We normalized our data in Fig. 2 based on the PMMoV concentration in each sample and its deviation from the median value. Figure 4 shows that the data after normalization are much less noisy and match the upward trend of clinical COVID-19 cases reported during 18 to 25 March (17). Thus, these data suggest that PMMoV may be useful as an internal reference to reduce sample variance.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

PMMoV normalization for sewage samples collected from 18 March to 25 March. (A) Original SARS-CoV-2 titers averaged from N1, N2, and N3 primers in the sewage sample. (B) Relative SARS-CoV-2 titers against PMMoV concentration in each of the samples. Data shown are averaged from N1, N2, and N3 primers. Influent samples from Northern (N) and Southern (S) influents are shown separately for each date.

These data suggest an order of magnitude estimate of approximately 100 viral genome copies per ml of sewage where titers in wastewater varied from 57 to 303 copies/ml, and from 21 to 506 copies/ml after PMMoV normalization. From these findings and analyses, we can draw the following deductions concerning disease prevalence. First, we note that any rigorous conclusions depend on a number of factors that are unknown, such as the timeline and load of fecal shedding, loss of viral particles in sewage lines, and the exact loss of viral RNA during experimental procedures, and thus, additional experiments will be required to calibrate these numbers. Nonetheless, we can estimate an abundance based on the lowest observed values across these samples of ∼10 copies/ml. If we assume typical stool sizes of 200 g, diluted into an average daily flow volume of 1.36 × 109 liter, and a population of 2.3 × 106 individuals each producing one stool per day, and we further assume that there is no loss of viral RNA in sewer lines and that excreted viruses are fully suspended in sewage, then we expect the viral titer in feces to be about 3,000 times higher than that in sampled raw sewage, or about 30,000 particles per ml.

Estimates of viral load in stool from positive patients are still a matter of uncertainty, where median levels of 12,882 copies/g (18) reflect a range from nondetection of SARS-CoV-2 RNA in some stool samples and levels as high as 600,000 to 30,000,000 viral genomes per ml or g of fecal material in other samples (18, 19). Estimating an average of 600,000 viral genomes per ml in stool (19) would suggest that roughly 5% of all fecal samples in the treatment facility catchment were positive for SARS-CoV-2 in the 18 to 25 March period, a number much higher than the 0.026% confirmed for the state of Massachusetts (a similar prevalence is obtained using individual counties represented by the wastewater treatment facility’s catchment or statewide estimates) on 25 March (17). Similarly, taking the extreme of stool concentration of 30 million copies/ml (18) suggests that we would have a prevalence of 0.1%, closer to, but still higher than, the number of confirmed clinical cases. Additional data on viral shedding in stool over the clinical course of the disease in patients with documented SARS-CoV-2 may be required to better interpret these findings.

DISCUSSION

Here, we established and verified a protocol to detect and quantify SARS-CoV-2 titers in wastewater samples. Starting with pasteurization, raw sewage samples were first filtered and precipitated with PEG-NaCl for viral enrichment, followed by RNA extraction, reverse transcription, and real-time PCR with SARS-CoV-2-specific primers. With this method, we detected SARS-CoV-2 in all of the 10 samples from 18 to 25 March in the influents of a wastewater treatment plant in Massachusetts. This simple protocol does not rely on expensive chemicals or materials and thus could be widely employed for viral detection in wastewater samples. Furthermore, we avoided the competitive usage of commercial kits like viral RNA extraction kits and one-step RT-PCR enzyme master mix kits that were recommended by the CDC protocols widely used in clinical testing (10).

Viral titers were quantified based on standard curves for the three primer sets (see Fig. S2 in the supplemental material). However, a complete quantification of viral titers needs further investigation of recovery rates for each step including the pasteurization, filtering, viral precipitation, RNA extraction, and reverse transcription. Recently, Ahmed et al. used the beta-coronavirus murine hepatitis virus (MHV) as a closely related surrogate for SARS-CoV-2 to quantify the recovery rates from wastewater; using various methods including PEG precipitation, they found that MHV recoveries ranged from 26.7 to 65.7% (20). To adjust the variations from the initial sample concentrations and RNA losses during these experimental procedures, we used PMMoV, the most abundant RNA virus in human stool (13, 16), as a reference. Our results showed that the PMMoV concentrations varied across the 12 samples with a standard deviation of less than 1.34 units, which also fall in the variation range of PMMoV concentrations in different seasons in a previous study (15). However, studies with more wastewater samples over a longer period are needed to understand whether the PMMoV loading is stable over time for the catchment population.

This discrepancy between confirmed cases and observed viral titers could arise from a number of factors. First, our estimates are based on numbers and assumptions that are currently subject to significant uncertainty, as discussed above, and therefore could be inaccurate. However, we note that our calculations are conservative because we (i) apply the lower limit of our observed viral titers in wastewater; (ii) assume no loss in viral titer due to degradation, sample processing, or RNA extraction; and (iii) consider high viral loads per infected stool, thus requiring fewer positive stools to generate the observed wastewater titers. The estimates of viral titer per positive stool may be too low if a small number of individuals shed very high levels of virus. While this can be determined, it will require more extensive testing of individual stool specimens from patients with documented COVID-19 disease in order to arrive at an accurate estimate. Given the importance of assessing the fraction of SARS-CoV-2 infections that present with symptoms, we note that our results are consistent with the idea that a significant fraction of cases are not detected with current testing algorithms, and that this fraction may include a large number of asymptomatic individuals (21, 22). Limitations of testing capacity and stringent criteria for testing would lead to an underestimate of clinically confirmed cases, whereas wastewater-based surveillance provides an anonymous, unbiased sample of the entire infected population, potentially explaining the magnitude of the difference we observe.

To discriminate between these possibilities, additional catchments of different sizes should be tested. By integrating data on the presence or absence of viral particles across catchments of varying size within the same geographic region, it may be possible to estimate disease prevalence independent from knowing the average viral titer in infected stool. For example, if disease prevalence is 1 in 10,000, then approximately 50% of catchment areas representing ∼7,000 individuals would be positive (although cases may cluster within households, so additional experiments or modeling would be needed to derive precise numbers). These experiments require sampling upstream areas in the sewage system, using specialized equipment to strategically capture time-integrated samples, and are under way. Once these experiments are completed in a small number of regions, they could be used to estimate disease prevalence in other regions using wastewater treatment facility data alone.

Conclusions.These data demonstrate the feasibility of measuring SARS-CoV-2 in wastewater using a method that does not require materials that are in high demand for individual clinical testing. The implications of this research are that wastewater-based surveillance can be leveraged to detect population-level prevalence of SARS-CoV-2 in cities and municipalities across the world. In a setting where in-person testing may not be available or may be overly cost-prohibitive, longitudinal analysis of wastewater can provide population-level estimates of the burden of SARS-CoV-2.

These data may help inform decisions surrounding the advancement or scale-back of social distancing and quarantine efforts based on wastewater catchment-level estimations of prevalence. Additionally, wastewater collection at the municipal or community level may allow for more granular detection of SARS-CoV-2 in cities with lower COVID-19 disease burden, thereby functioning as an early warning system to help preemptively enact public health measures prior to the widespread onset of disease.

MATERIALS AND METHODS

Sewage samples and pasteurization.Twenty-four-hour composite samples (1 liter) of raw sewage from 18 to 25 March were taken from a major urban wastewater treatment facility in Massachusetts. Those samples are from two catchment areas (Southern and Northern influents) except biobanked samples on 8 and 11 January, which were a mixture of Southern and Northern influents. In total, we collected 12 wastewater samples including 2 samples in January (8 and 11 January) and 10 samples from 18 to 25 March. Upon initial receipt, samples were placed in the biosafety cabinet with UV for 20 min and then pasteurized in a 60°C water bath for 90 min to inactivate the virus. Previous studies showed that pasteurization could effectively inactivate the virus without compromising sample quality (23, 24). Pasteurized samples were then used for viral precipitation, and the remaining samples were stored at 4°C.

Prefiltration evaluation.We tested the SARS-CoV-2 signals in the unfiltered samples, filtrate collected by filtration through a 0.22-μm polyether sulfone membrane (catalog no. SCGP00525; Millipore Sigma), and solid materials on the filter, respectively. Unfiltered sewage samples and filtrate were then precipitated with polyethylene glycol 8000 and NaCl and processed as described below. Solid materials on the filter were directly resuspended with 1.5 ml TRIzol and lysed for 10 min at room temperature. Undissolved materials were removed by centrifugation before the RNA extraction steps. The strongest and most consistent SARS-CoV-2 signal was detected from the PEG-precipitated filtrate fraction of the Northern and Southern samples, which was higher than the signal from the unfiltered samples or the solid fraction. No viral signal was detected in the solid fraction on the filter for the Southern sample, whereas only ∼10% of SARS-CoV-2 signal for N1 primer was found on the solid materials on the filter for the Northern sample (Table 2). Thus, we concluded that there was little viral RNA on the filters and used the filtrate for detection.

Viral precipitation.Pasteurized sample was first filtered through a 0.2-μm membrane (catalog no. SCGP00525; Millipore Sigma) to remove bacterial cells and debris. For the initial test, we used 80 ml filtrate of two samples on 18 March (Southern and Northern) for the viral precipitation with polyethylene glycol 8000 (8% [wt/vol]; Millipore Sigma) and NaCl (0.3 M; Millipore Sigma). Samples were shaken at room temperature for about 15 min (until the chemicals were fully dissolved) and then centrifuged at 12,000 × g for 2 h or until a pellet was visible. The viral pellet was then resuspended in 1.5 ml TRIzol reagent (ThermoFisher) for RNA extraction. After validation of this method (Fig. 1 and Table 2), we used 40 ml filtrate for the viral precipitation and processed all the 12 samples in the same way.

RNA extraction and RT-qPCR.RNA was extracted from the filtered sewage samples using the TRIzol-chloroform approach. Briefly, samples resuspended in the TRIzol reagent were thoroughly mixed with 300 μl chloroform for 1 min and incubated at room temperature for 5 min before centrifugation (16,000 × g, 15 min, 4°C). Six hundred microliters of aqueous phase was transferred to a new 1.5-ml tube and thoroughly mixed with an equal volume of isopropanol. After centrifugation at 16,000 × g for 10 min, the supernatant was discarded and the pellet was washed twice with 75% ethanol. Thirty microliters of diethyl pyrocarbonate (DEPC) water was used to recover the RNA.

cDNA was synthesized by the reverse transcription with random hexamers (IDT) based on manufacturer’s protocol (NEB catalog no. M0368). PCR includes TaqMan Fast Advanced master mix (4444557; ThermoFisher Scientific); CDC N1, N2, and N3 primer-probes (IDT); and cDNA as a template. The qPCR was carried out for 48 cycles on a Bio-Rad CFX96 real-time PCR detection system (Bio-Rad) based on the following program: polymerase activation (95°C for 2 min) and PCR (48 cycles, denaturation at 95°C for 1 s, and annealing/extension at 55°C for 30 s). Three replicates were performed for each primer, mean values were reported, and the coefficient of variation between replicates was <8%. The positive control for qPCR and construction of the standard curve for N1, N2, and N3 assays was a plasmid (initial concentration: 2 × 105 copies/μl; IDT catalog no. 10006625) containing the complete nucleocapsid gene from SARS-CoV-2.

Output from the qPCR assay for the pepper mild mottle virus (PMMoV) (primers shown in Table 1) was used to calibrate the SARS-CoV-2 titers. To adjust the SARS-CoV-2 viral titers for each sample, we first calculated the deviation of its PMMoV copies from the median of PMMoV copies in all samples, i.e., deviation factor = 10̂[k × (sample CT − median CT)], where k is the slope of the standard curve and equals −0.2991 (amplification efficiency is 99.11% for the PMMoV primer set) based on our test. We then divided the SARS-CoV-2 viral titers by this deviation. We noted that viral titers for 24 March samples (Southern and Northern) increased by 1.6- to 3.5-fold after PMMoV adjustment (Fig. 3). This is because the PMMoV concentrations for the two samples were lower than the median of PMMoV copies in the 12 sewage samples, suggesting that the two samples were more dilute than the other samples. It is interesting that there was significant precipitation in the Boston area during 23 to 24 March (https://www.timeanddate.com/weather/). The inflows from snow and heavy rainfall might explain the diluted wastewater samples on 24 March.

PCR, gel electrophoresis, and Sanger sequencing.To test SARS-CoV-2 identity in the wastewater samples and avoid potential contamination from N gene positive control (IDT catalog no. 10006625), we amplified an S gene fragment using primers shown in Table 1. PCR was performed with KAPA Hi-Fi Hot Start ready mix 2× (Fisher Scientific; KK2602), and cDNA was used as the template with the following program: 95°C for 5 min and 95°C for 30 s, 55°C for 30 s, and 72°C for 20 s for 42 cycles. PCR products were run on a 2% precast agarose E-gel (ThermoFisher; G402002). The positive control was cDNA reversed transcribed from synthetic SARS-CoV-2 RNA (Twist catalog no. 102024). PCR products were sent to Genewiz (Cambridge, MA) for purification and Sanger sequencing with forward and reverse primers (Table 1), and low-quality sequence adjacent to sequencing primers and terminal N’s were removed. Sequencing results were subjected to BLAST search with NCBI BLASTN (https://blast.ncbi.nlm.nih.gov).

ACKNOWLEDGMENTS

We thank the wastewater treatment facility who worked with us in providing the samples for analysis; Penny Chisholm (MIT) and Allison Coe (MIT) for access to equipment and other supplies, Feng Zhang (Broad Institute) for reagents; Shandrina Burns for logistical support; Katelyn Foppe for sample transportation; and Marc Lipsitch (HSPH), Karina Gin (NUS), Lee Ching Ng (NEA), and Stefan Wuertz, Mats Leifels, and Cheng Dan (NTU) for helpful discussion. We also thank Kyle Bibby (University of Notre Dame) for suggestions on PMMoV.

W.H. was supported by the National Institute of General Medical Sciences. This work was supported by the Center for Microbiome Informatics and Therapeutics and Intra-CREATE Thematic Grant (Cities) grant NRF2019-THE001-0003a to J.T. and E.J.A., and funding from the Massachusetts Consortium on Pathogen Readiness (MassCPR) and China Evergrande Group (T.B.E., P.R.C., M.M., and E.J.A.).

FOOTNOTES

    • Received July 2, 2020.
    • Accepted July 6, 2020.
  • Copyright © 2020 Wu et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Wu Y,
    2. Guo C,
    3. Tang L,
    4. Hong Z,
    5. Zhou J,
    6. Dong X,
    7. Yin H,
    8. Xiao Q,
    9. Tang Y,
    10. Qu X,
    11. Kuang L,
    12. Fang X,
    13. Mishra N,
    14. Lu J,
    15. Shan H,
    16. Jiang G,
    17. Huang X
    . 2020. Prolonged presence of SARS-CoV-2 viral RNA in faecal samples. Lancet Gastroenterol Hepatol 5:434–435. doi:10.1016/S2468-1253(20)30083-2.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Wang XW,
    2. Li J,
    3. Guo T,
    4. Zhen B,
    5. Kong Q,
    6. Yi B,
    7. Li Z,
    8. Song N,
    9. Jin M,
    10. Xiao W,
    11. Zhu X,
    12. Gu C,
    13. Yin J,
    14. Wei W,
    15. Yao W,
    16. Liu C,
    17. Li J,
    18. Ou G,
    19. Wang M,
    20. Fang T,
    21. Wang G,
    22. Qiu Y,
    23. Wu H,
    24. Chao F,
    25. Li J
    . 2005. Concentration and detection of SARS coronavirus in sewage from Xiao Tang Shan Hospital and the 309th Hospital of the Chinese People’s Liberation Army. Water Sci Technol 52:213–221. doi:10.2166/wst.2005.0266.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. Blomqvist S,
    2. Bassioni LE,
    3. Nasr EMEM,
    4. Paananen A,
    5. Kaijalainen S,
    6. Asghar H,
    7. de Gourville E,
    8. Roivainen M
    . 2012. Detection of imported wild polioviruses and of vaccine-derived polioviruses by environmental surveillance in Egypt. Appl Environ Microbiol 78:5406–5409. doi:10.1128/AEM.00491-12.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Hellmér M,
    2. Paxéus N,
    3. Magnius L,
    4. Enache L,
    5. Arnholm B,
    6. Johansson A,
    7. Bergström T,
    8. Norder H
    . 2014. Detection of pathogenic viruses in sewage provided early warnings of hepatitis A virus and norovirus outbreaks. Appl Environ Microbiol 80:6771–6781. doi:10.1128/AEM.01981-14.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Kopel E,
    2. Kaliner E,
    3. Grotto I
    . 2014. Lessons from a public health emergency—importation of wild poliovirus to Israel. N Engl J Med 371:981–983. doi:10.1056/NEJMp1406250.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Berchenko Y,
    2. Manor Y,
    3. Freedman LS,
    4. Kaliner E,
    5. Grotto I,
    6. Mendelson E,
    7. Huppert A
    . 2017. Estimation of polio infection prevalence from environmental surveillance data. Sci Transl Med 9:eaaf6786. doi:10.1126/scitranslmed.aaf6786.
    OpenUrlFREE Full Text
  7. 7.↵
    1. Duintjer Tebbens RJ,
    2. Zimmermann M,
    3. Pallansch MA,
    4. Thompson KM
    . 2017. Insights from a systematic search for information on designs, costs, and effectiveness of poliovirus environmental surveillance systems. Food Environ Virol 9:361–382. doi:10.1007/s12560-017-9314-4.
    OpenUrlCrossRef
  8. 8.↵
    1. Endo N,
    2. Ghaeli N,
    3. Duvallet C,
    4. Foppe K,
    5. Erickson TB,
    6. Matus M,
    7. Chai PR
    . 2020. Rapid assessment of opioid exposure and treatment in cities through robotic collection and chemical analysis of wastewater. J Med Toxicol 16:195–203. doi:10.1007/s13181-019-00756-5.
    OpenUrlCrossRef
  9. 9.↵
    1. Zhang F,
    2. Abudayyeh O,
    3. Gootenberg J
    . 2020. A protocol for detection of COVID-19 using CRISPR diagnostics. v.20200321.
  10. 10.↵
    CDC. 2020. Coronavirus disease 2019 (COVID-19). Centers for Disease Control and Prevention, Atlanta, GA. https://www.cdc.gov/coronavirus/2019-ncov/lab/rt-pcr-panel-primer-probes.html.
  11. 11.↵
    1. Rabenau HF,
    2. Cinatl J,
    3. Morgenstern B,
    4. Bauer G,
    5. Preiser W,
    6. Doerr HW
    . 2005. Stability and inactivation of SARS coronavirus. Med Microbiol Immunol 194:1–6. doi:10.1007/s00430-004-0219-0.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Zhang T,
    2. Breitbart M,
    3. Lee WH,
    4. Run J-Q,
    5. Wei CL,
    6. Soh SWL,
    7. Hibberd ML,
    8. Liu ET,
    9. Rohwer F,
    10. Ruan Y
    . 2006. RNA viral community in human feces: prevalence of plant pathogenic viruses. PLoS Biol 4:e3. doi:10.1371/journal.pbio.0040003.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Rosario K,
    2. Symonds EM,
    3. Sinigalliano C,
    4. Stewart J,
    5. Breitbart M
    . 2009. Pepper mild mottle virus as an indicator of fecal pollution. Appl Environ Microbiol 75:7261–7267. doi:10.1128/AEM.00410-09.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Hamza H,
    2. Rizk NM,
    3. Gad MA,
    4. Hamza IA
    . 2019. Pepper mild mottle virus in wastewater in Egypt: a potential indicator of wastewater pollution and the efficiency of the treatment process. Arch Virol 164:2707–2713. doi:10.1007/s00705-019-04383-x.
    OpenUrlCrossRef
  15. 15.↵
    1. Kitajima M,
    2. Iker BC,
    3. Pepper IL,
    4. Gerba CP
    . 2014. Relative abundance and treatment reduction of viruses during wastewater treatment processes—identification of potential viral indicators. Sci Total Environ 488–489:290–296. doi:10.1016/j.scitotenv.2014.04.087.
    OpenUrlCrossRef
  16. 16.↵
    1. Kitajima M,
    2. Sassi HP,
    3. Torrey JR
    . 2018. Pepper mild mottle virus as a water quality indicator. NPJ Clean Water 1:19. doi:10.1038/s41545-018-0019-5.
    OpenUrlCrossRef
  17. 17.↵
    Commonwealth of Massachusetts. 2020. COVID-19 updates and information. https://www.mass.gov/info-details/covid-19-updates-and-information.
  18. 18.↵
    1. Wölfel R,
    2. Corman VM,
    3. Guggemos W,
    4. Seilmaier M,
    5. Zange S,
    6. Müller MA,
    7. Niemeyer D,
    8. Jones TC,
    9. Vollmar P,
    10. Rothe C,
    11. Hoelscher M,
    12. Bleicker T,
    13. Brünink S,
    14. Schneider J,
    15. Ehmann R,
    16. Zwirglmaier K,
    17. Drosten C,
    18. Wendtner C
    . 2020. Virological assessment of hospitalized patients with COVID-2019. Nature 581:465–410. doi:10.1038/s41586-020-2196-x.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Zhang N,
    2. Gong Y,
    3. Meng F,
    4. Bi Y,
    5. Yang P,
    6. Wang F
    . 2020. Virus shedding patterns in nasopharyngeal and fecal specimens of COVID-19 patients. medRxiv doi:10.1101/2020.03.28.20043059.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Ahmed W,
    2. Bertsch PM,
    3. Bivins A,
    4. Bibby K,
    5. Farkas K,
    6. Gathercole A,
    7. Haramoto E,
    8. Gyawali P,
    9. Korajkic A,
    10. McMinn BR,
    11. Mueller JF,
    12. Simpson SL,
    13. Smith WJM,
    14. Symonds EM,
    15. Thomas KV,
    16. Verhagen R,
    17. Kitajima M
    . 2020. 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 739:139960. doi:10.1016/j.scitotenv.2020.139960.
    OpenUrlCrossRef
  21. 21.↵
    1. Li R,
    2. Pei S,
    3. Chen B,
    4. Song Y,
    5. Zhang T,
    6. Yang W,
    7. Shaman J
    . 2020. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science 368:489–493. doi:10.1126/science.abb3221.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Bai Y,
    2. Yao L,
    3. Wei T,
    4. Tian F,
    5. Jin D-Y,
    6. Chen L,
    7. Wang M
    . 2020. Presumed asymptomatic carrier transmission of COVID-19. JAMA 323:1406–1407. doi:10.1001/jama.2020.2565.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Wang T,
    2. Lien C,
    3. Liu S,
    4. Selveraj P
    . 2020. Effective heat inactivation of SARS-CoV-2. medRxiv doi:10.1101/2020.04.29.20085498.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Liu Y,
    2. Cao Z,
    3. Chen M,
    4. Zhong Y,
    5. Luo Y,
    6. Shi G,
    7. Xiang H,
    8. Luo J,
    9. Zhou H
    . 2020. Effect of heat inactivation on real-time reverse transcription PCR of the SARS-CoV-2 detection. medRxiv doi:10.1101/2020.05.19.20101469.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
SARS-CoV-2 Titers in Wastewater Are Higher than Expected from Clinically Confirmed Cases
Fuqing Wu, Jianbo Zhang, Amy Xiao, Xiaoqiong Gu, Wei Lin Lee, Federica Armas, Kathryn Kauffman, William Hanage, Mariana Matus, Newsha Ghaeli, Noriko Endo, Claire Duvallet, Mathilde Poyet, Katya Moniz, Alex D. Washburne, Timothy B. Erickson, Peter R. Chai, Janelle Thompson, Eric J. Alm
mSystems Jul 2020, 5 (4) e00614-20; DOI: 10.1128/mSystems.00614-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
SARS-CoV-2 Titers in Wastewater Are Higher than Expected from Clinically Confirmed Cases
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
SARS-CoV-2 Titers in Wastewater Are Higher than Expected from Clinically Confirmed Cases
Fuqing Wu, Jianbo Zhang, Amy Xiao, Xiaoqiong Gu, Wei Lin Lee, Federica Armas, Kathryn Kauffman, William Hanage, Mariana Matus, Newsha Ghaeli, Noriko Endo, Claire Duvallet, Mathilde Poyet, Katya Moniz, Alex D. Washburne, Timothy B. Erickson, Peter R. Chai, Janelle Thompson, Eric J. Alm
mSystems Jul 2020, 5 (4) e00614-20; DOI: 10.1128/mSystems.00614-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

COVID-19 prevalence
PMMoV
SARS-CoV-2
viral titers
wastewater

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077