The impact of signal variability on epidemic growth rate estimation from wastewater surveillance data

Ewan Colman, Rowland R. Kao
doi: https://doi.org/10.1101/2023.03.07.23286904

Abstract

Background Testing samples of waste water for markers of infectious disease became a widespread method of surveillance during the COVID-19 pandemic. While these data generally correlate well with other indicators of national prevalence, samples that cover localised regions tend to be highly variable over short time scales.

Methods We introduce a procedure for estimating the realtime growth rate of pathogen prevalence using time series data from wastewater sampling. The number of copies of a target gene found in a sample is modelled as time-dependent random variable whose distribution is estimated using maximum likelihood. The output depends on a hyperparameter that controls the sensitivity to variability in the underlying data. We apply this procedure to data reporting the number of copies of the N1 gene of SARS-CoV-2 collected at water treatment works across Scotland between February 2021 and February 2023.

Results The real-time growth rate of the SARS-CoV-2 prevalence is estimated at 121 wastewater sampling sites covering a diverse range of locations and population sizes. We find that the sensitivity of the fitting procedure to natural variability determines its reliability in detecting the early stages of an epidemic wave. Applying the procedure to hospital admissions data, we find that changes in the growth rate are detected an average of 2 days earlier in wastewater than in hospital admissions data.

Conclusion We provide a robust method to generate reliable estimates of epidemic growth from highly variable data. Applying this method to samples collected at wastewater treatment works provides highly responsive situational awareness to inform public health.

Infectious disease surveillance is a fundamental pillar of public health management. Epidemiological data collected with regularity across multiple locations provides an awareness of the progress of a disease as it spreads through a population. Robust surveillance systems to support this will prevent adverse outcomes through early action, without undue risk of raising false alarms. An emerging technology that promises improve surveillance capabilities pathogen detection in wastewater [1]. This has been shown to provide early detection of epidemic trends compared to other surveillance sources [2–6]. The technology has potential to detect a range human diseases [7–13] and reach otherwise overlooked communities [14, 15].

A central question for policy makers is whether incidence of a disease is increasing or decreasing. Since an upward trend in new infections may be a reason to introduce restrictions or adjust provisions to hospitals, it is important to have reliable methods for determining epidemic trajectories in real time from the available data [16–18]. For wastewater based epidemiology using RT-qPCR this is particularly challenging as there are factors beyond changes in the underlying population prevalence that affect the observed signal an observed increase in the gene copy count between samples taken on different days does not necessarily indicate an increase in the number of infected people [19, 20].

To mitigate the effects of natural variability, smoothing can be applied over time and space [18, 21–24]. This approach can, however, have the unintended consequence of distorting some important features of the data, such as inflection points where there is a sudden change in the epidemic trajectory, or a sharp increase in one location that subsequently spreads to the wider region. We propose that an approach based on fitting an epidemiologically realistic model is more appropriate for these data as we would like the outcome to represent the underlying epidemic trajectory as accurately as possible. The question is how can we filter out the high levels of variability observed at the local scale, while retaining the ability to quickly detect changes in epidemic trajectory?

This paper introduces a method to interpret the highly variable wastewater sample values and applies it to data from the Covid-19 pandemic in Scotland. We start by discussing the various causes of variability affecting data collected from wastewater samples and describe the distribution of of values we expect to observe at a given prevalence. We then introduce a function to measure the likelihood of any given epidemic trajectory in relation to the observed wastewater data and describe a procedure to find the best fitting function from a class of exponential models. We compare results across different sampling locations, optimise the process for early detection of epidemic waves, and finally compare results to hospital admissions data.

1 Methods

1.1 Data

From February 2021 to the time of writing (February 2023), samples from Waste Water Treatment Plant (WWTP) were regularly taken by the Scottish Environment Protection Agency (SEPA) to detect fragments of SARS-Cov-2 virus RNA [25, 26]. Samples from sewerage influent were taken at 122 sampling locations across Scotland approximately two to three days each week depending on the location [27]. Each Sample was collected using a refrigerated autosampler that obtained a fixed volume of influent every hour over a 24-h period (08:00 to 08:00). Composite 24-h samples were mixed and concentrated before viral RNA was extracted using commercial kits. SARS-CoV-2 N1 gene average concentrations (gene copies/l) are obtained using RT-qPCR.

WWTP data were provided by SEPA via a publicly available portal [28]. For each sample taken the data provides the date, the WWTP, and the N1 reported value obtained from RT-qPCR. Additional data for fluid flow for a limited number of WWTPs were provided though restricted access from SEPA. Since inflow to each WWTP comes from households and other premises connected to the sewerage network, the recorded values relate to SARSCoV-2 infections in a specific geographical catchment area. Scottish water provided shapefiles giving the border of the region corresponding to each WWTP.

Restricted access to hospital admission data was provided by Public Heath Scotland. The hospital admission database includes information on individuals admitted to hospital in Scotland including primary reason for admission, other contributing conditions, date of admission, and the place of residence provided at the datazone level (each datazone contains approximately 700 residents). We take all admissions for which Covid-19 was listed as primary or other contributing reason for admission, and aggregate to the datazone level before using the proportions of overlap between the datazones and the WWTP catchment shapefiles we aggregate to the geographical area of the each WWTP catchment.

1.2 Variability

Our goal is to use data from WWTPs to observe changes in the growth rate of infections in the population. Since there are many factors that may affect the quantity of viral RNA in a sample, we cannot interpret every observed change in the reported value to be a true indicator of a change in the number of people infected. We often see the reported value fall to a fraction of what it was in the previous sample, only to jump to a high value again in the next sample. Rapid changes like this are too fast to be explained by transmission dynamics based on our current understanding of SARS-CoV-2 epidemiology.

Table 1 lists some of the causes of variability suggested in the literature. While each of these will have some effect, we generally lack sufficient data to control each of these factors; for the purpose of obtaining growth estimates it is more practical to accept that there is some inherent randomness in the wastewater sampling process. Our approach is to model the variabilty we expect to see at a given prevalence. We hypothesize that the variability between samples is predominantly the result of variability between the amount that infected individuals contribute to the wastewater sample, and this enables us to describe mathematically the distribution of outcomes we expect to see from a sample given for a given population prevalence. This distribution will provide a basis for modelling the underlying trend in prevalence.

View this table:
Table 1: Factors influencing temporal variation in the observed RNA count in wastewater samples.

Naively, it is natural to assume that the N1 gene copies contained in a 24-hour sample is contributed equally from every individual in the catchment who is infected and currently shedding the virus. In reality this may not be the case; firstly, there is some variability there in the amount each individual dispenses to waste water (related to the amount of RNA fragments passing through the gastro-intestinal tract) [34–36], and secondly, as illustrated in Figure 1, fluids and fecal matter do not dilute perfectly throughout the sewerage influent, adding additional variability to the concentration of viral RNA received by the sampler. In an extreme case, the sampler my by chance pick up a large concentration of viral RNA from one highly infected individual giving an atypically large value [33].

Figure 1:
Figure 1:

Flow of viral RNA fragments. RNA fragments shown as red circles. The amount from each household may vary across time and depend on severity of outcomes. When samples are taken, this variability is expected to affect the outcome of taking a sample. Composite sampling may reduce variability but the quantity of viral RNA detected may still be dominated by RNA originating from a single infection.

Let us suppose that on a specific day there are n individuals who are infected and contributing to some wastewater sample. We let xi represent the quantity of viral RNA contributed by individual i. We model the values xi for i = 1, 2, …, n as an i.i.d random variables from a distribution with mean µ and standard deviation s. It follows from the central limit theorem that the cumulative quantity of viral RNA detected Embedded Image is modelled as a random variable that follows a Normal distribution: Embedded Image If we assume that samples drawn from this distribution at any two times are independent, the variability over short time scales can be explained by a sufficiently high value of ns2.

This generalises to other data series including the number of reported cases or the number of hospital admissions. In both cases, each infected individual contributes either 0 or 1 to the total count depending whether they report a positive test (for case count data) or have symptoms that lead to hospitalisation (for admissions data). The total count can therefore be modelled as an Binomial distribution and, for sufficiently large numbers contributing to the signal, approximated by the Normal distribution. The methods described in this paper are therefore applicable to a wide range of sources of epidemiological surveillance.

1.3 Likelihood of underlying prevalence model

We allow n to be a function of time, n(t), representing the number of people who are shedding virus in wastewater on day t. Similarly v(t) is the total quantity of viral RNA in the sampler on day t (assuming a sample was taken). The function f (t) = n(t)µ transforms Eq. (2) to Embedded Image where D = s2/µ is a property of the distribution of individual contributions, xi, known as the index of dispersion. Note that we have assumed that this distribution is timeindependent.

Our approach to estimating the epidemic trajectory from WWTP data is to find the function f (t) and parameter D that achieves the maximum likelihood for the given data. Representing the data as two vectors t = {t1, t2, …, tk}, and y = {y1, y2, …, yk}, where ti is the time the ith sample was taken and yi is the reported quantity of viral RNA, the best fitting parameter values are found by solving Embedded Image where Embedded Image The function f (t) in Eq. (4), can be any time series over the period where data were collected. Since our goal is to remove unwanted variability and reveal the underlying epidemic trajectory, we choose to constrain f to a class of functions that are feasible given the epidemiology of the disease.

1.4 Epidemic trajectory model

We limit our choice of the epidemic trajectory f to the class of exponential functions with growth rates that change at discrete points in time. Supposing there are m times when the exponential growth rate changes, hereafter referred to as change points, we let c be a vector of length m representing the times of each change point, and r be a vector of length m + 1 where the ith entry is the growth during the interval from ci to the next change point. Thus we have Embedded Image where A1 is a free parameter (the intercept), and the other Ak are determined using the formula Embedded Image to ensure that f is continuous.

1.5 Optimization procedure

An estimate for the epidemic trajectory and natural variability is found by solving Eq. 4 over the class of functions of the form expressed in Eq. (6). The free parameters are the initial quantity of wastewater RNA (A1), the number of times the growth rate changes (m), the times of each of these changes (c), the corresponding growth rates (r), and the index of dispersion (D). If we want to fit to data over a period containing several change points, the space of possible parameter combinations becomes very large and computationally expensive to optimize by brute force. We tried several off-the-shelf methods for solving the for optimization but these struggled to consistently produce useful output. We then developed our own heuristic procedure for solving Eq. (4) which we detail here.

The outline of the procedure, illustrated in Figure 2, is as follows: A simple exponential curve is fitted to the first few data points using a hill-climbing algorithm. The remaining data points are added sequentially. Each time a new data point is added, the hill-climbing algorithm is re-applied to optimize the parameters. A candidate change point is considered between the most recent change point (or the start of the data if no change points have been added) and the time of the newly added data point. The hill-climbing algorithm is re-applied to optimize the parameters (including the time of the candidate change point and the growth rate after the change). If the inclusion of the change point yields an increase in the likelihood function, Eq. (5), that exceeds a given threshold value then the new change point is added to the parameter set of f. This is repeated until all the data have been added.

Figure 2:
Figure 2:

Model fitting. The model is fitted to the data by iteratively adding data points and optimizing the parameters using hill-climbing at each iteration. In each iteration a candidate change point (CP) is added. If there is an improvement in the likelihood resulting from the added exceeds some threshold value then the change point becomes part of the accepted model.

The choice of the threshold value has consequences for the outcome. If set close to zero it will produce a model that is highly sensitive to variability in the underlying data. Higher thresholds will find longer-term trends in the underlying data while being less sensitive to the natural fluctuations of the waste water signal. The choice of threshold value can only be determined when the goal of the analysis has been decided. In this paper the goal is to detect waves of SARS-CoV-2 infection as early as possible.

Here we describe the hill-climbing part of the algorithm. The first step selects a discrete set of possible values for each parameter. The dispersion parameter (D) can take values from 0.5% to 10% of the maximum quantity of viral RNA across the series at intervals of 0.5%, the time of each change point (c) can be any day between the midpoint of its current value and the previous change point (or start of the data), and the midpoint of its current value and the next change point (or end of the data at that iteration of the algorithm).

The initial quantity of wastewater RNA (A1) exponential growth rates (r) are not included as parameters explicitly, instead we use the gene copy value at each change point and use them to calculate growth rates after the optimization. The range for the gene copy values is 100 equally spaced values from 0 to 1.5 multiplied by the maximum value in the series. The hill-climb proceeds by perturbing parameters one at a time to adjacent values within the ranges described. It then chooses to accept the perturbation which yielded the best improvement. This continues until no improvement can be made.

2 Results

2.1 Model fitting

We use the 10 largest WWTPS (by population) to explore the effects of the threshold on the outcome. As shown in Fig. 3, higher threshold values produce fewer changes in the modelled trajectory and correspondingly higher values for the dispersion (the measure of non-epidemic variability). The threshold value of 8 was provides outcomes that are neither too sensitive to natural variability yet still responsive to real changes in the epidemic trajectory. We applied the procedure with this threshold to 121 WWTPs, 1 was omitted for containing no positive values for the quantity of viral RNA. The quality of the outcome is measured by the mean log likelihood; the maximum likelihood result divided by the number of data points. This varied across sites, however, it did not vary by the number of samples in the data (Pearson’s r,p = 0.36) and there was no significant effect of the catchment area population size (Pearson’s r,p = 0.16), implying that the procedure can be applied successfully in geographical regions as small as 1500 people (Figure S1).

Figure S1:
Figure S1:

Model fit does not appear to be affected by the number of samples or the population of the WWTP catchment. The log likelihood of model using the maximal likelihood parameter values is divided by the number of data points to give the mean log likelihood. Each point in these plots represents a WWTP.

Figure 3:
Figure 3:

Relationship between the sensitivity threshold and the number of trajectory changes. A. The top panel shows the raw data reporting the quantity of viral RNA in the collected samples. The panels below show the outcome of model fitting using maximum likelihood estimation for two values of the threshold hyperparameter. B The dependence of the number of change points on the threshold hyperparameter for the 10 largest WWTPs by catchment population. C. Relationship between the number of change points and the index of dispersion. With fewer change points, variability around the mean must necessarily be larger in order to explain the distribution of observed quantities.

2.2 Real-time estimation of the epidemic growth rate

The real-time growth rate at time t is the rate value (rm+1, in the function f given by Eq. (6)) obtained after fitting the model to time series up to and including time t. Figure 4A shows the estimated growth rates over time for one WWTP for three different threshold values. This example illustrates a general trend that we observe across all the treatment sites: that the responsiveness of the real-time estimate to changes in the data is dependent on the chosen threshold value.

Figure 4:
Figure 4:

Real-time growth rate estimates. A. Estimates for the Seafield WWTP at three threshold values. The solid line in these figures shows the estimated growth rate using the data available at each the given point in time. Higher thresholds are less responsive to high frequency variability in the data. Regions shaded red indicate periods uninterrupted positive growth. B. The number of periods uninterrupted positive growth. Low threshold values produce a high frequency of “false alarms” - changes to the trajectory that are not sustained - higher threshold values are more stable but periods of positive growth tend to begin later.

As measure of the responsiveness to changes in the epidemic trajectory we count the number of periods of uninterrupted positive growth. We see that periods of growth tend to start earlier for lower thresholds, potentially providing an earlier warning that prevalence will rise, but also causing a larger frequency of “false alarms” where the direction of growth will change for a short period before changing again. From Figure 4 we see that the number of waves predicted is highly sensitive to the threshold value when it is below 6, while threshold values of 8 or greater provide relatively consistent outcomes. We therefore suggest that values from 6 to 8 provide high responsiveness while remaining insensitive to high frequency fluctuations in the data.

2.3 Comparison to hospital admissions data

Throughout the Covid-19 epidemic the number of reported cases has provided a useful metric to help forecast hospital demand [37]. As the number of reported test results is highly dependent on test-seeking behaviour, wastewater can potentially be a less biased way to provide the same early warning. To address this we compare the growth rates obtained from wastewater data to those obtained from hospital admission data to see if one time series predicts the other.

To obtain growth rates from hospital admission data we applied our method to time series data reporting the number of people admitted to hospital for COVID-19 who live in the catchment of each of the 10 largest WWTPs (Figure 5A). We found that a threshold value of 8 provided a real-time growth rate that captured the major changes to the epidemic trajectory without being overly responsive to fluctuations in the data. It predicted a similar number of waves to the waste water output for the same threshold value (within 1 across all WWTPs tested) and therefore provided outputs suitable for comparison. Figure 5B shows an example of this for one WWTP.

Figure 5:
Figure 5:

Comparison to admissions data. A. Model fitted to data from hospital admissions of individuals living in the catchment of the Seafield WWTP using identical methods to those applied to wastewater data in previous figures. B. The real time growth rate estimated using samples from the Seafield WWTP and admissions for the corresponding catchment area. C. Pearson’s r correlation between the time-series produced in panel B and for the 10 largest WWTP catchments in Scotland over a range of lags. Circles are added at the bottom of the figure to show the value of the lag that yields the maximum correlation. The mean correlation over the 10 sites is shown as a dashed line. Positive lag implies that wastewater growth rates respond faster than those obtained from hospital admissions data.

We calculated the correlation between the time series over a range of lags. A lag of length l means that the growth rate value on day t in the wastewater time series is paired with the growth rate value on day t + l in the hospital admissions time series. The optimal lag is the value of l that give the highest correlation between the two series. A positive optimal lag implies that the signal coming from wastewater data responds quicker to changes in the epidemic trajectory than the signal coming from admissions data. We find that the optimal lag is positive for 6 of the 10 sites. Taking the mean of the Pearson’s r values over the 10 sites we find that a lag of 2 days achieved the largest correlation (Figure 5C).

2.4 Environmental factors

Rainfall entering the sewerage system is expected to dilute the concentration of viral material in the fluid leading to lower observed values. While this dilution effect will vary depending on the amount of rainfall, it can be difficult to observe directly as its magnitude is small in comparison to the variation in viral prevalence. To test whether fluid flow has an effect on the reported quantities of viral RNA, we compare fluid flow val-ues to the z-scores of the observed viral RNA quantities with respect to the modelled distribution (using a sensitivity threshold value of 8). The z-score is calculated by subtracting the model mean from the observed value and dividing by the standard deviation, in principal removing the effects of underlying epidemic trajectory.

Of the 15 sites with more than 100 recorded flow values, correlation between fluid flow and higher than expected viral RNA was found in 5 (Pearson’s r < 0, p < 0.05). In these cases the effect size is small in comparison to the overall amount of variability. For example, applying linear regression to data from Meadowhead, we find that that each additional 100 megalitres of flow per day reduces the mean sample value by 0.41 standard deviations (Figure S2).

Figure S2:
Figure S2:

The affect of fluid flow is small. We show the 10 largest WWTPs (by population) for which flow data are available. Each point in these plots represents a sample. The vertical axis shows the z-score for the sample with respect to the distribution that maximizes the model likelihood

3 Discussion

The goal of this work was to provide a methodology to interpret waste water data and help inform policy decisions. A major challenge to decision makers is deciding when a change in the reported data indicates a genuine change in the underlying epidemic trajectory, particularly when the data are subject to high levels of natural variability. Our approach accounts for the variability between samples observed in the reported RT-qPCR result and uses it as the basis for a model fitting procedure. The main outcome is an estimate of the growth rate of the epidemic that updates each time new data become available. Applying this procedure to data from sampling locations across Scotland, we found that samples collected from wastewater detect changes to the epidemic trajectory approximately two days earlier than can be observed in hospital admissions.

Despite the recent growth of waste water surveillance as a non-intrusive source of population health data, there are no clear guidelines around how these data should be used to inform policy. We note that deep analysis is not always needed; Polio, for example, is rare in the UK but confers a profound risk if it is not controlled if the virus is detected then action should be taken [38]. For highly prevalent diseases, on the other hand, deciding when to take action is less straight-forward; there is always a risk of mistaking a random fluctuation for a real spike in incidence, and this may be a reason for some to distrust this valuable source of data. Our approach of estimating the growth rate of an underlying model resolves this problem.

We highlight here some limitations of our analysis and areas for potential improvement. Firstly, the method has only been tested on SARS-CoV-2 data in Scotland during a period when prevalence was high. While the method can be applied in other settings we recommend first validating the outputs against other sources of surveillance data in those settings. Secondly, we made a specific choice to explore the exponential function as the prevalence curve. This assumption is reasonable in most cases as infectious diseases dynamics are typically driven by branching processes parameterized by the effective reproductive number. There may, however, be situations where other functions should be considered for example if we are to take into account the duration of shedding [21, 35]. Thirdly, adjusting the observed quantities of viral RNA using a normalisation factor (such as the fluid flow which we have shown to affect the model output) could potentially yield some improvement to the results [20, 27, 27, 39, 40]. Similarly, adjusting the hospital admissions data prior analysis to remove known periodic patterns may yield an improvement.

One advantage of our method is its applicability to data from individual sites distributed over a range of geographical locations. Better spatial coverage achieves earlier warnings of the start of an epidemic wave and better targeting of control measures [41]. As wastewater surveillance grows, we may see samples being taken at more WWTPs or even within the sewerage network, reaching ever smaller populations [42,43]. While our results demonstrate that surveillance for Covid-19 can be useful over catchment populations as low as 1500, we don’t know precisely where the limitations of the method lie. Future work should explore these limits, not just in terms of spatial resolution, but also the sampling frequency and underlying prevalence.

We hope that this work contributes to the wider goal of improving the way we mitigate the risks of disease through wastewater surveillance. While we are confident that the methods developed here are applicable to a range of current settings, they exist to be improved upon and adapted to the specific diseases and challenges in wastewater monitoring. The ongoing challenge is to build on the work of many environmental and biological scientists that have made these data sources possible, to push the continued development of mathematical and statistical tools, and help support decision makers and the population by providing the best possible scientific analysis.

Data Availability

Environmental data analysed in the present study are available at https://informatics.sepa.org.uk/RNAmonitoring/. Medical data were provided with restricted access subject to a data sharing agreement

https://informatics.sepa.org.uk/RNAmonitoring/

Code availability

https://github.com/EwanColman/

Estimating-epidemic-growth-rate-from-wastwater-data

Acknowledgments

We are grateful to the Scottish Environmental Protection Agency, Scottish Water, and Public Health Scotland for making data for this study available. We thank our partners in the Scottish Government Wastewater Monitoring Data Integration Subgroup for their insightful feedback. Financial support was received from The Wellcome Trust and the National COVID-19 Wastewater Epidemiology Surveillance Programme.

Footnotes

  • ↵* ecolman{at}ed.ac.uk

References

  1. [1].↵
    Joshua I. Levy, Kristian G. Andersen, Rob Knight, and Smruthi Karthikeyan. Wastewater surveillance for public health. Science, 379(6627):26–27, 2023.
    OpenUrl
  2. [2].↵
    Jordan Peccia, Alessandro Zulli, Doug E Brackney, Nathan D Grubaugh, Edward H Kaplan, Arnau Casanovas-Massana, Albert I Ko, Amyn A Malik, Dennis Wang, Mike Wang, et al. Measurement of sars-cov-2 rna in wastewater tracks community infection dynamics. Nature biotechnology, 38(10):1164–1167, 2020.
    OpenUrlPubMed
  3. [3].
    Candice L. Swift, Mirza Isanovic, Karlen E. Correa Velez, and R. Sean Norman. Sars-cov-2 concentration in wastewater consistently predicts trends in covid-19 case counts by at least two days across multiple wwtp scales. Environmental Advances, 11:100347, 2023.
    OpenUrl
  4. [4].
    Veljo Kisand, Peeter Laas, Kadi Palmik-Das, Kristel Panksep, Helen Tammert, Leena Albreht, Hille Allemann, Lauri Liepkalns, Katri Vooro, Christian Ritz, Vasili Hauryliuk, and Tanel Tenson. Prediction of covid-19 positive cases, a nation-wide sarscov-2 wastewater-based epidemiology study. Water Research, 231:119617, 2023.
    OpenUrl
  5. [5].
    Liang Zhao, Yangyang Zou, Randy E. David, Scott Withington, Stacey McFarlane, Russell A. Faust, John Norton, and Irene Xagoraraki. 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. Science of The Total Environment, 864:161152, 2023.
    OpenUrl
  6. [6].↵
    Manish Kumar, Guangming Jiang, Alok Kumar Thakur, Shreya Chatterjee, Tanushree Bhattacharya, Sanjeeb Mohapatra, Tushara Chaminda, Vinay Kumar Tyagi, Meththika Vithanage, Prosun Bhattacharya, Long D. Nghiem, Dibyendu Sarkar, Christian Sonne, and Jürgen Mahlknecht. Lead time of early warning by wastewater surveillance for covid-19: Geographical variations and impacting factors. Chemical Engineering Journal, 441:135936, 2022.
    OpenUrlCrossRef
  7. [7].↵
    Ananda Tiwari, Paula Kurittu, Ahmad I. Al-Mustapha, Viivi Heljanko, Venla Johansson, Ocean Thakali, Shyam Kumar Mishra, Kirsi-Maarit Lehto, Anssi Lipponen, Sami Oikarinen, Tarja Pitkänen, WastPan Study Group, Annamari Heikinheimo, Annika Länsivaara, Rafiqul Hyder, Erja Janhonen, Anna-Maria Hokajärvi, Anniina Sarekoski, Aleksi Kolehmainen, Soile Blomqvist, Kati Räisänen, Carita Savolainen Kopra, Teemu Möttönen, Oskari Luomala, Aapo Juutinen, Ocean Thakali, and Shyam Kumar Mishra. Wastewater surveillance of antibiotic-resistant bacterial pathogens: A systematic review. Frontiers in Microbiology, 13, 2022.
  8. [8].
    Sarah Nadeau, A.J. Devaux, Claudia Bagutti, Monica Alt, Evelyn Ilg Hampe, Melanie Kraus, Eva Würfel, Katrin N. Koch, Simon Fuchs, Sarah Tschudin-Sutter, Aurélie Holschneider, Christoph Ort, Chaoran Chen, Jana S. Huisman, Timothy R. Julian, and Tanja Stadler. Influenza transmission dynamics quantified from wastewater. medRxiv, 2023.
  9. [9].
    Inés Girón-Guzmán, Azahara Díaz-Reolid, Pilar Truchado, Albert Carcereny, David García-Pedemonte, Bruno Hernáez, Albert Bosch, Rosa María Pintó, Susana Guix, Ana Allende, Antonio Alcamí, Alba Pérez-Cataluña, and Gloria Sánchez. Spanish wastewater reveals the current spread of monkeypox virus. Water Research, 231:119621, 2023.
    OpenUrl
  10. [10].
    Marlene K. Wolfe, Alexander T. Yu, Dorothea Duong, Madhura S. Rane, Bridgette Hughes, Vikram Chan-Herur, Marisa Donnelly, Shua Chai, Bradley J. White, Duc J. Vugia, and Alexandria B. Boehm. Use of wastewater for mpox outbreak surveillance in california. New England Journal of Medicine, 0(0):ull, 0.
  11. [11].
    Casey Barber, Katherine Crank, Katerina Papp, Gabriel K. Innes, Bradley W. Schmitz, Jorge Chavez, Alessandro Rossi, and Daniel Gerrity. Community-scale wastewater surveillance of candida auris during an ongoing outbreak in southern nevada. Environmental Science & Technology, 57(4):1755–1763, 2023. PMID: 36656763.
    OpenUrlPubMed
  12. [12].
    Warish Ahmed, Aaron Bivins, Mikayla Stephens, Suzanne Metcalfe, Wendy J.M. Smith, Kwanrawee Sirikanchana, Masaaki Kitajima, and Stuart L. Simpson. Occurrence of multiple respiratory viruses in wastewater in queensland, australia: Potential for community disease surveillance. Science of The Total Environment, 864:161023, 2023.
    OpenUrl
  13. [13].↵
    Kattareeya Kumthip, Pattara Khamrin, Hiroshi Ushijima, and Niwat Maneekarn. Detection of six different human enteric viruses contaminating environmental water in chiang mai, thailand. Microbiology Spectrum, 0(0):e03512–22, 2023.
    OpenUrl
  14. [14].↵
    William M. Baldwin, Robert D. Dayton, Aaron W. Bivins, Rona S. Scott, Andrew D. Yurochko, John A. Vanchiere, Terry Davis, Connie L. Arnold, Jose.E.T. Asuncion, Mohammad A.N. Bhuiyan, Brandon Snead, William Daniel, Deborah G. Smith, Nicholas E. Goeders, Christopher G. Kevil, Jennifer Carroll, and Kevin S. Murnane. Highly socially vulnerable communities exhibit disproportionately increased viral loads as measured in community wastewater. Environmental Research, 222:115351, 2023.
    OpenUrl
  15. [15].↵
    Tatiana Prado, Gloria Rey-Benito, Marize Pereira Miagostovich, Maria Inês Zanoli Sato, Veronica Beatriz Rajal, Cesar Rossas Mota Filho, Alyne Duarte Pereira, Mikaela Renata Funada Barbosa, Camille Ferreira Mannarino, and Agnes Soares da Silva. Wastewater-based epidemiology for preventing outbreaks and epidemics in latin america – lessons from the past and a look to the future. Science of The Total Environment, 865:161210, 2023.
    OpenUrl
  16. [16].↵
    Junling Ma. Estimating epidemic exponential growth rate and basic reproduction number. Infectious Disease Modelling, 5:129–141, 2020.
    OpenUrl
  17. [17].
    Laura M Guzmán-Rincón, Edward M Hill, Louise Dyson, Michael J Tildesley, and Matt J Keeling. Bayesian estimation of real-time epidemic growth rates using gaussian processes: local dynamics of sars-cov-2 in england. medRxiv, pages 2022–01, 2022.
  18. [18].↵
    Jana S Huisman, Jérémie Scire, Lea Caduff, Xavier Fernandez-Cassi, Pravin Ganesanandamoorthy, Anina Kull, Andreas Scheidegger, Elyse Stachler, Alexandria B Boehm, Bridgette Hughes, et al. Wastewater-based estimation of the effective reproductive number of sars-cov-2. Environmental health perspectives, 130(5):057011, 2022.
    OpenUrlCrossRef
  19. [19].↵
    Matthew J. Wade, Anna Lo Jacomo, Elena Armenise, Mathew R. Brown, Joshua T. Bunce, Graeme J. Cameron, Zhou Fang, Kata Farkas, Deidre F. Gilpin, David W. Graham, Jasmine M.S. Grimsley, Alwyn Hart, Till Hoffmann, Katherine J. Jackson, David L. Jones, Chris J. Lilley, John W. McGrath, Jennifer M. McKinley, Cormac McSparron, Behnam F. Nejad, Mario Morvan, Marcos Quintela-Baluja, Adrian M.I. Roberts, Andrew C. Singer, Célia Souque, Vanessa L. Speight, Chris Sweetapple, David Walker, Glenn Watts, Andrew Weightman, and Barbara Kasprzyk-Hordern. Understanding and managing uncertainty and variability for wastewater monitoring beyond the pandemic: Lessons learned from the united kingdom national covid-19 surveillance programmes. Journal of Hazardous Materials, 424:127456, 2022.
    OpenUrlCrossRef
  20. [20].↵
    Mario Morvan, Anna Lo Jacomo, Celia Souque, Matthew J Wade, Till Hoffmann, Koen Pouwels, Chris Lilley, Andrew C Singer, Jonathan Porter, Nicholas P Evens, et al. An analysis of 45 largescale wastewater sites in england to estimate sars-cov-2 community prevalence. Nature Communications, 13(1):4313, 2022.
    OpenUrl
  21. [21].↵
    Fuqing Wu, Amy Xiao, Jianbo Zhang, Katya Moniz, Noriko Endo, Federica Armas, Richard Bonneau, Megan A. Brown, Mary Bushman, Peter R. Chai, Claire Duvallet, Timothy B. Erickson, Katelyn Foppe, Newsha Ghaeli, Xiaoqiong Gu, William P. Hanage, Katherine H. Huang, Wei Lin Lee, Mariana Matus, Kyle A. McElroy, Jonathan Nagler, Steven F. Rhode, Mauricio Santillana, Joshua A. Tucker, Stefan Wuertz, Shijie Zhao, Janelle Thompson, and Eric J. Alm. Sars-cov-2 rna concentrations in wastewater foreshadow dynamics and clinical presentation of new covid-19 cases. Science of The Total Environment, 805:150121, 2022.
    OpenUrlCrossRef
  22. [22].
    Hannah D. Greenwald, Lauren C. Kennedy, Adrian Hinkle, Oscar N. Whitney, Vinson B. Fan, Alexander Crits-Christoph, Sasha Harris-Lovett, Avi I. Flamholz, Basem Al-Shayeb, Lauren D. Liao, Matt Beyers, Daniel Brown, Alicia R. Chakrabarti, Jason Dow, Dan Frost, Mark Koekemoer, Chris Lynch, Payal Sarkar, Eileen White, Rose Kantor, and Kara L. Nelson. Tools for interpretation of wastewater sars-cov-2 temporal and spatial trends demonstrated with data collected in the san francisco bay area. Water Research X, 12:100111, 2021.
    OpenUrl
  23. [23].
    Maarten Nauta, Oliver McManus, Kristina Træholt Franck, Ellinor Lindberg Marving, Lasse Dam Rasmussen, Stine Raith Richter, and Steen Ethelberg. Early detection of local sars-cov-2 outbreaks by wastewater surveillance: a feasibility study. Epidemiology &; Infection, page 1–28, 2023.
  24. [24].↵
    Guangquan Li, Hubert Denise, Peter Diggle, Jasmine Grimsley, Chris Holmes, Daniel James, Radka Jersakova, Callum Mole, George Nicholson, Camila Rangel Smith, Sylvia Richardson, William Rowe, Barry Rowlingson, Fatemeh Torabi, Matthew J. Wade, and Marta Blangiardo. A spatio-temporal framework for modelling wastewater concentration during the covid-19 pandemic. Environment International, 172:107765, 2023.
    OpenUrl
  25. [25].↵
    Livia CT Scorza, Graeme J Cameron, Roisin Murray-Williams, David Findlay, Julie Bolland, Brindusa Cerghizan, Kirsty Campbell, David Thomson, Alexander Corbishley, David Gally, et al. Sarscov-2 rna levels in scotland’s wastewater. Scientific Data, 9(1):713, 2022.
    OpenUrl
  26. [26].↵
    Stephen F. Fitzgerald, Gianluigi Rossi, Alison S. Low, Sean P. McAteer, Brian O’Keefe, David Findlay, Graeme J. Cameron, Peter Pollard, Peter T. R. Singleton, George Ponton, Andrew C. Singer, Kata Farkas, Davey Jones, David W. Graham, Marcos Quintela-Baluja, Christine Tait-Burkard, David L. Gally, Rowland Kao, and Alexander Corbishley. Site specific relationships between covid-19 cases and sars-cov-2 viral load in wastewater treatment plant influent. Environmental Science & Technology, 55(22):15276–15286, 2021. PMID: 34738785.
    OpenUrlPubMed
  27. [27].↵
    Zhou Fang, Adrian M. I. Roberts, Claus-Dieter Mayer, Anastasia Frantsuzova, Jackie M. Potts, Graeme J. Cameron, Peter T. R. Singleton, and Iona Currie. Wastewater monitoring of COVID-19: a perspective from Scotland. Journal of Water and Health, 20(12):1688–1700, 12 2022.
    OpenUrl
  28. [28].↵
    Sepa rna monitoring dashboard. https://informatics.sepa.org.uk/RNAmonitoring/.
  29. [29].
    Van Vo, Richard L. Tillett, Katerina Papp, Ching-Lan Chang, Anthony Harrington, Michael Moshi, Edwin C. Oh, and Daniel Gerrity. Detection of the Omicron BA.1 Variant of SARS-CoV-2 in Wastewater From a Las Vegas Tourist Area. JAMA Network Open, 6(2):e230550–e230550, 02 2023.
    OpenUrl
  30. [30].
    Jean-Baptiste Burnet, Henry-Michel Cauchie, Cécile Walczak, Nathalie Goeders, and Leslie Ogorzaly. Persistence of endogenous rna biomarkers of sars-cov-2 and pmmov in raw wastewater: Impact of temperature and implications for wastewater-based epidemiology. Science of The Total Environment, 857:159401, 2023.
    OpenUrl
  31. [31].
    Warish Ahmed, Stuart L. Simpson, Paul M. Bertsch, Kyle Bibby, Aaron Bivins, Linda L. Blackall, Sílvia Bofill-Mas, Albert Bosch, João Brandão, Phil M. Choi, Mark Ciesielski, Erica Donner, Nishita D’Souza, Andreas H. Farnleitner, Daniel Gerrity, Raul Gonzalez, John F. Griffith, Pradip Gyawali, Charles N. Haas, Kerry A. Hamilton, Hapuarachchige Chanditha Hapuarachchi, Valerie J. Harwood, Rehnuma Haque, Greg Jackson, Stuart J. Khan, Wesaal Khan, Masaaki Kitajima, Asja Korajkic, Giuseppina La Rosa, Blythe A. Layton, Erin Lipp, Sandra L. McLellan, Brian McMinn, Gertjan Medema, Suzanne Metcalfe, Wim G. Meijer, Jochen F. Mueller, Heather Murphy, Coleen C. Naughton, Rachel T. Noble, Sudhi Payyappat, Susan Petterson, Tarja Pitkänen, Veronica B. Rajal, Brandon Reyneke, Fernando A. Roman, Joan B. Rose, Marta Rusinõl, Michael J. Sadowsky, Laura Sala-Comorera, Yin Xiang Setoh, Samendra P. Sherchan, Kwanrawee Sirikanchana, Wendy Smith, Joshua A. Steele, Rosalie Sabburg, Erin M. Symonds, Phong Thai, Kevin V. Thomas, Josh Tynan, Simon Toze, Janelle Thompson, Andy S. Whiteley, Judith Chui Ching Wong, Daisuke Sano, Stefan Wuertz, Irene Xagoraraki, Qian Zhang, Amity G. Zimmer-Faust, and Orin C. Shanks. Minimizing errors in rt-pcr detection and quantification of sars-cov-2 rna for wastewater surveillance. Science of The Total Environment, 805:149877, 2022.
    OpenUrl
  32. [32].
    Federica Armas, Franciscus Chandra, Wei Lin Lee, Xiaoqiong Gu, Hongjie Chen, Amy Xiao, Mats Leifels, Stefan Wuertz, Eric J Alm, and Janelle Thompson. Contextualizing wastewater-based surveillance in the covid-19 vaccination era. Environment International, 171:107718, 2023.
    OpenUrl
  33. [33].↵
    Xander Bertels, Phaedra Demeyer, Siel Van den Bogaert, Tim Boogaerts, Alexander L.N. van Nuijs, Peter Delputte, and Lies Lahousse. Factors influencing sars-cov-2 rna concentrations in wastewater up to the sampling stage: A systematic review. Science of The Total Environment, 820:153290, 2022.
    OpenUrl
  34. [34].↵
    Till Hoffmann and Justin Alsing. Faecal shedding models for sars-cov-2 rna among hospitalised patients and implications for wastewater-based epidemiology. medRxiv, 2023.
  35. [35].↵
    Fuminari Miura, Masaaki Kitajima, and Ryosuke Omori. Duration of sars-cov-2 viral shedding in faeces as a parameter for wastewater-based epidemiology: Re-analysis of patient data using a shedding dynamics model. Science of The Total Environment, 769:144549, 2021.
    OpenUrl
  36. [36].↵
    Amarylle S. van Doorn, Berrie Meijer, Chris M. A. Frampton, Murray L. Barclay, and Nanne K. H. de Boer. Systematic review with meta-analysis: Sars-cov-2 stool testing and the potential for faecal-oral transmission. Alimentary Pharmacology & Therapeutics, 52(8):1276–1288, 2020.
    OpenUrlPubMed
  37. [37].↵
    Sophie Meakin, Sam Abbott, Nikos Bosse, James Munday, Hugo Gruson, Joel Hellewell, Katharine Sherratt, and Sebastian Funk. Comparative assessment of methods for short-term forecasts of covid-19 hospital admissions in england at the local level. BMC medicine, 20(1):1–15, 2022.
    OpenUrl
  38. [38].↵
    K. M. O’Reilly, N. C. Grassly, D. J. Allen, M. Bannister-Tyrrell Cameron, A. I. Carrion Martin, M. Ramsay, R. Pebody, and M. Zambon. Surveillance optimisation to detect poliovirus in the pre-eradication era: a modelling study of england and wales. Epidemiology & Infection, 148:e157, 2020.
    OpenUrl
  39. [39].
    Qingyu Zhan, Kristina M. Babler, Mark E. Sharkey, Ayaaz Amirali, Cynthia C. Beaver, Melinda M. Boone, Samuel Comerford, Daniel Cooper, Elena M. Cortizas, Benjamin B. Currall, Jonathan Foox, George S. Grills, Erin Kobetz, Naresh Kumar, Jennifer Laine, Walter E. Lamar, Alejandro M.A. Mantero, Christopher E. Mason, Brian D. Reding, Maria Robertson, Matthew A. Roca, Krista Ryon, Stephan C. Schürer, Bhavarth S. Shukla, Natasha Schaefer Solle, Mario Stevenson, John J. Tallon Jr., Collette Thomas, Tori Thomas, Dušica Vidović, Sion L. Williams, Xue Yin, and Helena M. Solo-Gabriele. Relationships between sars-cov-2 in wastewater and covid-19 clinical cases and hospitalizations, with and without normalization against indicators of human waste. ACS ES&T Water, 2(11):1992–2003, 2022.
    OpenUrl
  40. [40].
    Rasha Maal-Bared, Yuanyuan Qiu, Qiaozhi Li, Tiejun Gao, Steve E. Hrudey, Sudha Bhavanam, Norma J. Ruecker, Erik Ellehoj, Bonita E. Lee, and Xiaoli Pang. Does normalization of sars-cov-2 concentrations by pepper mild mottle virus improve correlations and lead time between wastewater surveillance and clinical data in alberta (canada): comparing twelve sars-cov-2 normalization approaches. Science of The Total Environment, 856:158964, 2023.
    OpenUrlCrossRef
  41. [41].↵
    Ewan Colman, Petter Holme, Hiroki Sayama, and Carlos Gershenson. Efficient sentinel surveillance strategies for preventing epidemics on networks. PLOS Computational Biology, 15(11):1–19, 11 2019.
    OpenUrl
  42. [42].↵
    Prakit Saingam, Bo Li, Bao Nguyen Quoc, Tanisha Jain, Andrew Bryan, and Mari K.H. Winkler. Wastewater surveillance of sarscov-2 at intra-city level demonstrated high resolution in tracking covid-19 and calibration using chemical indicators. Science of The Total Environment, 866:161467, 2023.
    OpenUrl
  43. [43].↵
    Gyuhyon Cha, Katherine E. Graham, Kevin J. Zhu, Gouthami Rao, Blake G. Lindner, Kumru Kocaman, Seongwook Woo, Isabelle D’amico, Lilia R. Bingham, Jamie M. Fischer, Camryn I. Flores, John W. Spencer, Pranav Yathiraj, Hayong Chung, Shweta Biliya, Naima Djeddar, Liza J. Burton, Samantha J. Mascuch, Joe Brown, Anton Bryksin, Ameet Pinto, Janet K. Hatt, and Konstantinos T. Konstantinidis. Parallel deployment of passive and composite samplers for surveillance and variant profiling of sars-cov-2 in sewage. Science of The Total Environment, 866:161101, 2023.
    OpenUrl
Back to top
Posted March 09, 2023.
Download PDF
Data/Code
Email
The impact of signal variability on epidemic growth rate estimation from wastewater surveillance data
Ewan Colman, Rowland R. Kao
medRxiv 2023.03.07.23286904; doi: https://doi.org/10.1101/2023.03.07.23286904
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area