Abstract
Since the start of the COVID-19 pandemic in 2020, wastewater surveillance programs were established, or upscaled, in many countries around the world and have proven to be a cost-effective way of monitoring infectious disease pathogens. Many of these programs use RT-qPCR, and quantify the viral concentrations in samples based on standard curves, by including preparations of a reference material with known nucleic acid or virus concentrations in the RT-qPCR analyses. In high-throughput monitoring programs it is possible to combine data from multiple previous runs, circumventing the need for duplication and resulting in decreased costs and prolonged periods during which the reference material is obtained from the same batch. However, over time, systematic shifts in standard curves are likely to occur. This would affect the reliability and usefulness of wastewater surveillance as a whole. We aim to find an optimal combination of standard curve data to compensate for run-to-run measurement variance while ensuring enough flexibility to capture systematic longitudinal shifts. Based on more than 4000 observations obtained with the CDC N1 and N2 assays, taken as a part of the National Sewage Surveillance program at the Dutch National Institute for Public Health and the Environment, we show that seasonal and long-term shifts in RT-qPCR efficiency and sensitivity occur. We find that in our setting, using five days of standard-curve data to quantify, results in the least error prone curve or best approximation. This results in differences up to 100% in quantified viral loads when averaged out over a nationwide program of >300 treatment plants. Results show that combining standard curves from a limited set of runs can be a valid approach to quantification without obscuring the trends in the viral load of interest.