skip to main content


Title: Emerging investigator series: meta-analyses on SARS-CoV-2 viral RNA levels in wastewater and their correlations to epidemiological indicators
Background: recent applications of wastewater-based epidemiology (WBE) have demonstrated its ability to track the spread and dynamics of COVID-19 at the community level. Despite the growing body of research, quantitative synthesis of SARS-CoV-2 RNA levels in wastewater generated from studies across space and time using diverse methods has not been performed. Objective: the objective of this study is to examine the correlations between SARS-CoV-2 RNA levels in wastewater and epidemiological indicators across studies, stratified by key covariates in study methodologies. In addition, we examined the association of proportions of positive detections in wastewater samples and methodological covariates. Methods: we systematically searched the Web of Science for studies published by February 16th, 2021, performed a reproducible screening, and employed mixed-effects models to estimate the levels of SARS-CoV-2 viral RNA quantities in wastewater samples and their correlations to the case prevalence, the sampling mode (grab or composite sampling), and the wastewater fraction analyzed ( i.e. , solids, solid–supernatant mixtures, or supernatants/filtrates). Results: a hundred and one studies were found; twenty studies (671 biosamples and 1751 observations) were retained following a reproducible screening. The mean positivity across all studies was 0.68 (95%-CI, [0.52; 0.85]). The mean viral RNA abundance was 5244 marker copies per mL (95%-CI, [0; 16 432]). The Pearson correlation coefficients between the viral RNA levels and case prevalence were 0.28 (95%-CI, [0.01; 0.51]) for daily new cases or 0.29 (95%-CI, [−0.15; 0.73]) for cumulative cases. The fraction analyzed accounted for 12.4% of the variability in the percentage of positive detections, followed by the case prevalence (9.3% by daily new cases and 5.9% by cumulative cases) and sampling mode (0.6%). Among observations with positive detections, the fraction analyzed accounted for 56.0% of the variability in viral RNA levels, followed by the sampling mode (6.9%) and case prevalence (0.9% by daily new cases and 0.8% by cumulative cases). While the sampling mode and fraction analyzed both significantly correlated with the SARS-CoV-2 viral RNA levels, the magnitude of the increase in positive detection associated with the fraction analyzed was larger. The mixed-effects model treating studies as random effects and case prevalence as fixed effects accounted for over 90% of the variability in SARS-CoV-2 positive detections and viral RNA levels. Interpretations: positive pooled means and confidence intervals in the Pearson correlation coefficients between the SARS-CoV-2 viral RNA levels and case prevalence indicators provide quantitative evidence that reinforces the value of wastewater-based monitoring of COVID-19. Large heterogeneities among studies in proportions of positive detections, viral RNA levels, and Pearson correlation coefficients suggest a strong demand for methods to generate data accounting for cross-study heterogeneities and more detailed metadata reporting. Large variance was explained by the fraction analyzed, suggesting sample pre-processing and fractionation as a direction that needs to be prioritized in method standardization. Mixed-effects models accounting for study level variations provide a new perspective to synthesize data from multiple studies.  more » « less
Award ID(s):
2047470
NSF-PAR ID:
10405113
Author(s) / Creator(s):
;  ; ; ; ; ; ;
Date Published:
Journal Name:
Environmental Science: Water Research & Technology
Volume:
8
Issue:
7
ISSN:
2053-1400
Page Range / eLocation ID:
1391 to 1407
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Background

    Estimating the cumulative incidence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is essential for setting public health policies. We leveraged deidentified Massachusetts newborn screening specimens as an accessible, retrospective source of maternal antibodies for estimating statewide seroprevalence in a nontest-seeking population.

    Methods

    We analyzed 72 117 newborn specimens collected from November 2019 through December 2020, representing 337 towns and cities across Massachusetts. Seroprevalence was estimated for the Massachusetts population after correcting for imperfect test specificity and nonrepresentative sampling using Bayesian multilevel regression and poststratification.

    Results

    Statewide seroprevalence was estimated to be 0.03% (90% credible interval [CI], 0.00–0.11) in November 2019 and rose to 1.47% (90% CI: 1.00–2.13) by May 2020, following sustained SARS-CoV-2 transmission in the spring. Seroprevalence plateaued from May onward, reaching 2.15% (90% CI: 1.56–2.98) in December 2020. Seroprevalence varied substantially by community and was particularly associated with community percent non-Hispanic Black (β = .024; 90% CI: 0.004–0.044); i.e., a 10% increase in community percent non-Hispanic Black was associated with 27% higher odds of seropositivity. Seroprevalence estimates had good concordance with reported case counts and wastewater surveillance for most of 2020, prior to the resurgence of transmission in winter.

    Conclusions

    Cumulative incidence of SARS-CoV-2 protective antibody in Massachusetts was low as of December 2020, indicating that a substantial fraction of the population was still susceptible. Maternal seroprevalence data from newborn screening can inform longitudinal trends and identify cities and towns at highest risk, particularly in settings where widespread diagnostic testing is unavailable.

     
    more » « less
  2. Abstract

    Wastewater surveillance has proven to be an effective tool to monitor the transmission and emergence of infectious agents at a community scale. Workflows for wastewater surveillance generally rely on concentration steps to increase the probability of detection of low-abundance targets, but preconcentration can substantially increase the time and cost of analyses while also introducing additional loss of target during processing. To address some of these issues, we conducted a longitudinal study implementing a simplified workflow for SARS-CoV-2 detection from wastewater, using a direct column-based extraction approach. Composite influent wastewater samples were collected weekly for 1 year between June 2020 and June 2021 in Athens-Clarke County, Georgia, USA. Bypassing any concentration step, low volumes (280 µl) of influent wastewater were extracted using a commercial kit, and immediately analyzed by RT-qPCR for the SARS-CoV-2 N1 and N2 gene targets. SARS-CoV-2 viral RNA was detected in 76% (193/254) of influent samples, and the recovery of the surrogate bovine coronavirus was 42% (IQR: 28%, 59%). N1 and N2 assay positivity, viral concentration, and flow-adjusted daily viral load correlated significantly with per-capita case reports of COVID-19 at the county-level (ρ = 0.69–0.82). To compensate for the method’s high limit of detection (approximately 106–107 copies l−1 in wastewater), we extracted multiple small-volume replicates of each wastewater sample. With this approach, we detected as few as five cases of COVID-19 per 100 000 individuals. These results indicate that a direct-extraction-based workflow for SARS-CoV-2 wastewater surveillance can provide informative and actionable results.

     
    more » « less
  3. Kretzschmar, Mirjam E. (Ed.)
    Background Development of an effective antiviral drug for Coronavirus Disease 2019 (COVID-19) is a global health priority. Although several candidate drugs have been identified through in vitro and in vivo models, consistent and compelling evidence from clinical studies is limited. The lack of evidence from clinical trials may stem in part from the imperfect design of the trials. We investigated how clinical trials for antivirals need to be designed, especially focusing on the sample size in randomized controlled trials. Methods and findings A modeling study was conducted to help understand the reasons behind inconsistent clinical trial findings and to design better clinical trials. We first analyzed longitudinal viral load data for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) without antiviral treatment by use of a within-host virus dynamics model. The fitted viral load was categorized into 3 different groups by a clustering approach. Comparison of the estimated parameters showed that the 3 distinct groups were characterized by different virus decay rates ( p -value < 0.001). The mean decay rates were 1.17 d −1 (95% CI: 1.06 to 1.27 d −1 ), 0.777 d −1 (0.716 to 0.838 d −1 ), and 0.450 d −1 (0.378 to 0.522 d −1 ) for the 3 groups, respectively. Such heterogeneity in virus dynamics could be a confounding variable if it is associated with treatment allocation in compassionate use programs (i.e., observational studies). Subsequently, we mimicked randomized controlled trials of antivirals by simulation. An antiviral effect causing a 95% to 99% reduction in viral replication was added to the model. To be realistic, we assumed that randomization and treatment are initiated with some time lag after symptom onset. Using the duration of virus shedding as an outcome, the sample size to detect a statistically significant mean difference between the treatment and placebo groups (1:1 allocation) was 13,603 and 11,670 (when the antiviral effect was 95% and 99%, respectively) per group if all patients are enrolled regardless of timing of randomization. The sample size was reduced to 584 and 458 (when the antiviral effect was 95% and 99%, respectively) if only patients who are treated within 1 day of symptom onset are enrolled. We confirmed the sample size was similarly reduced when using cumulative viral load in log scale as an outcome. We used a conventional virus dynamics model, which may not fully reflect the detailed mechanisms of viral dynamics of SARS-CoV-2. The model needs to be calibrated in terms of both parameter settings and model structure, which would yield more reliable sample size calculation. Conclusions In this study, we found that estimated association in observational studies can be biased due to large heterogeneity in viral dynamics among infected individuals, and statistically significant effect in randomized controlled trials may be difficult to be detected due to small sample size. The sample size can be dramatically reduced by recruiting patients immediately after developing symptoms. We believe this is the first study investigated the study design of clinical trials for antiviral treatment using the viral dynamics model. 
    more » « less
  4. Abstract Background

    SARS-CoV-2 is an RNA virus responsible for the coronavirus disease 2019 (COVID-19) pandemic. Viruses exist in complex microbial environments, and recent studies have revealed both synergistic and antagonistic effects of specific bacterial taxa on viral prevalence and infectivity. We set out to test whether specific bacterial communities predict SARS-CoV-2 occurrence in a hospital setting.

    Methods

    We collected 972 samples from hospitalized patients with COVID-19, their health care providers, and hospital surfaces before, during, and after admission. We screened for SARS-CoV-2 using RT-qPCR, characterized microbial communities using 16S rRNA gene amplicon sequencing, and used these bacterial profiles to classify SARS-CoV-2 RNA detection with a random forest model.

    Results

    Sixteen percent of surfaces from COVID-19 patient rooms had detectable SARS-CoV-2 RNA, although infectivity was not assessed. The highest prevalence was in floor samples next to patient beds (39%) and directly outside their rooms (29%). Although bed rail samples more closely resembled the patient microbiome compared to floor samples, SARS-CoV-2 RNA was detected less often in bed rail samples (11%). SARS-CoV-2 positive samples had higher bacterial phylogenetic diversity in both human and surface samples and higher biomass in floor samples. 16S microbial community profiles enabled high classifier accuracy for SARS-CoV-2 status in not only nares, but also forehead, stool, and floor samples. Across these distinct microbial profiles, a single amplicon sequence variant from the genusRothiastrongly predicted SARS-CoV-2 presence across sample types, with greater prevalence in positive surface and human samples, even when compared to samples from patients in other intensive care units prior to the COVID-19 pandemic.

    Conclusions

    These results contextualize the vast diversity of microbial niches where SARS-CoV-2 RNA is detected and identify specific bacterial taxa that associate with the viral RNA prevalence both in the host and hospital environment.

     
    more » « less
  5. Santos, Ricardo (Ed.)
    We developed and implemented a framework for examining how molecular assay sensitivity for a viral RNA genome target affects its utility for wastewater-based epidemiology. We applied this framework to digital droplet RT-PCR measurements of SARS-CoV-2 and Pepper Mild Mottle Virus genes in wastewater. Measurements were made using 10 replicate wells which allowed for high assay sensitivity, and therefore enabled detection of SARS-CoV-2 RNA even when COVID-19 incidence rates were relatively low (~10 −5 ). We then used a computational downsampling approach to determine how using fewer replicate wells to measure the wastewater concentration reduced assay sensitivity and how the resultant reduction affected the ability to detect SARS-CoV-2 RNA at various COVID-19 incidence rates. When percent of positive droplets was between 0.024% and 0.5% (as was the case for SARS-CoV-2 genes during the Delta surge), measurements obtained with 3 or more wells were similar to those obtained using 10. When percent of positive droplets was less than 0.024% (as was the case prior to the Delta surge), then 6 or more wells were needed to obtain similar results as those obtained using 10 wells. When COVID-19 incidence rate is low (~ 10 −5 ), as it was before the Delta surge and SARS-CoV-2 gene concentrations are <10 4 cp/g, using 6 wells will yield a detectable concentration 90% of the time. Overall, results support an adaptive approach where assay sensitivity is increased by running 6 or more wells during periods of low SARS-CoV-2 gene concentrations, and 3 or more wells during periods of high SARS-CoV-2 gene concentrations. 
    more » « less