Implementing an adaptive, two-tiered SARS-CoV-2 wastewater surveillance program on a university campus using passive sampling.

Building-level wastewater-based surveillance (WBS) has been increasingly applied upstream from wastewater treatment plants to conduct targeted monitoring for SARS-CoV-2. In this study, a two-tiered, trigger-based wastewater surveillance program was developed on a university campus to monitor dormitory wastewater. The objective was to determine if passive sampling with cotton gauze as a sampling medium could be used to support institution-level public health action. Two nucleocapsid gene targets (N1 and N2) of SARS-CoV-2 as well as the endogenous fecal indicator pepper mild mottle virus (PMMoV) were quantified using RT-qPCR. >500 samples were analyzed during two contrasting surveillance periods. In the Fall of 2021 community viral burden was low and a tiered sampling network was able to isolate individual clinical cases at the building-scale. In the Winter of 2022 wastewater signals were quickly elevated by the emergence of the highly transmissible SARS-CoV-2 Omicron (B.1.1.529) variant. Prevalence of SARS-CoV-2 shifted surveillance objectives from isolating cases to monitoring trends, revealing both the benefits and limitations of a tiered surveillance design under different public health situations. Normalization of SARS-CoV-2 by PMMoV was not reflective of upstream population differences, suggesting saturation of the material occurred during the exposure period. The passive sampling method detected nearly all known clinical cases and in one instance was able to identify one pre-symptomatic individual days prior to confirmation by clinical test. Comparisons between campus samplers and municipal wastewater influent suggests that the spread of COVID-19 on the campus was similar to that of the broader community. The results demonstrate that passive sampling is an effective tool that can produce semi-quantitative data capable of tracking temporal trends to guide targeted public health decision-making at an institutional level. Practitioners of WBS can utilize these results to inform surveillance program designs that prioritize efficient resource use and rapid reporting.

Assessment of seasonality and normalization techniques for wastewater-based surveillance in Ontario, Canada.

Introduction: Wastewater-based surveillance is at the forefront of monitoring for community prevalence of COVID-19, however, continued uncertainty exists regarding the use of fecal indicators for normalization of the SARS-CoV-2 virus in wastewater. Using three communities in Ontario, sampled from 2021-2023, the seasonality of a viral fecal indicator (pepper mild mottle virus, PMMoV) and the utility of normalization of data to improve correlations with clinical cases was examined. Methods: Wastewater samples from Warden, the Humber Air Management Facility (AMF), and Kitchener were analyzed for SARS-CoV-2, PMMoV, and crAssphage. The seasonality of PMMoV and flow rates were examined and compared by Season-Trend-Loess decomposition analysis. The effects of normalization using PMMoV, crAssphage, and flow rates were analyzed by comparing the correlations to clinical cases by episode date (CBED) during 2021. Results: Seasonal analysis demonstrated that PMMoV had similar trends at Humber AMF and Kitchener with peaks in January and April 2022 and low concentrations (troughs) in the summer months. Warden had similar trends but was more sporadic between the peaks and troughs for PMMoV concentrations. Flow demonstrated similar trends but was not correlated to PMMoV concentrations at Humber AMF and was very weak at Kitchener (r = 0.12). Despite the differences among the sewersheds, unnormalized SARS-CoV-2 (raw N1-N2) concentration in wastewater (n = 99-191) was strongly correlated to the CBED in the communities (r = 0.620-0.854) during 2021. Additionally, normalization with PMMoV did not improve the correlations at Warden and significantly reduced the correlations at Humber AMF and Kitchener. Flow normalization (n = 99-191) at Humber AMF and Kitchener and crAssphage normalization (n = 29-57) correlations at all three sites were not significantly different from raw N1-N2 correlations with CBED. Discussion: Differences in seasonal trends in viral biomarkers caused by differences in sewershed characteristics (flow, input, etc.) may play a role in determining how effective normalization may be for improving correlations (or not). This study highlights the importance of assessing the influence of viral fecal indicators on normalized SARS-CoV-2 or other viruses of concern. Fecal indicators used to normalize the target of interest may help or hinder establishing trends with clinical outcomes of interest in wastewater-based surveillance and needs to be considered carefully across seasons and sites.

Intersex manifestation in the rainbow darter (Etheostoma caeruleum): Are adult male fish susceptible to developing and recovering from intersex after exposure to endocrine active compounds?

For over a decade, intersex has been observed in rainbow darter (RD) (Etheostoma caeruleum) populations living downstream wastewater treatment plants (WWTPs) in the Grand River, Ontario, Canada. To further our understanding of intersex development in adult male fish, the current study addressed three objectives: i) can intersex be induced in adult male fish, ii) is there a specific window of exposure when adult male fish are more susceptible to developing intersex, and iii) can pre-exposed adult male fish recover from intersex? To assess intersex induction in adult male fish, wild male RD were exposed in the laboratory for 22 weeks (during periods of spawning, gonadal regression, and gonadal recrudescence) to environmentally relevant concentrations of 17α-ethinylestradiol (EE2) including nominal 0, 1, and 10 ng/L. Intersex rates and severity at 10 ng/L EE2 were similar to those observed historically in adult male populations living downstream WWTPs in the Grand River and confirmed previous predictions that 1-10 ng/L EE2 would cause these adverse effects. To assess a window of sensitivity in developing intersex, male RD were exposed to nominal 0, 1 or 10 ng/L EE2 for 4 weeks during three different periods of gonadal development, including (i) spawning, (ii) early recrudescence and (iii) late recrudescence. These short-term exposures revealed that intersex incidence and severity were greater when RD were exposed while gonads were fully developed (during spawning) compared to periods of recrudescence. To assess if RD recover from intersex, wild fish were collected downstream WWTPs in the Grand River and assessed for intersex both before and after a 22-week recovery period in clean water that included gonadal regression and recrudescence. Results showed that fish did not recover from intersex, with intersex rates and severity similar to those both before and after the transition to clean water. This study further advances our knowledge on intersex manifestation in adult male fish including their sensitivity to endocrine active compounds during different periods of their annual reproductive cycle and their limited ability to recover from intersex after onset of the condition.

Omicron COVID-19 Case Estimates Based on Previous SARS-CoV-2 Wastewater Load, Regional Municipality of Peel, Ontario, Canada.

We determined correlations between SARS-CoV-2 load in untreated water and COVID-19 cases and patient hospitalizations before the Omicron variant (September 2020-November 2021) at 2 wastewater treatment plants in the Regional Municipality of Peel, Ontario, Canada. Using pre-Omicron correlations, we estimated incident COVID-19 cases during Omicron outbreaks (November 2021-June 2022). The strongest correlation between wastewater SARS-CoV-2 load and COVID-19 cases occurred 1 day after sampling (r = 0.911). The strongest correlation between wastewater load and COVID-19 patient hospitalizations occurred 4 days after sampling (r = 0.819). At the peak of the Omicron BA.2 outbreak in April 2022, reported COVID-19 cases were underestimated 19-fold because of changes in clinical testing. Wastewater data provided information for local decision-making and are a useful component of COVID-19 surveillance systems.

A comparative analysis of the partitioning behaviour of SARS-CoV-2 RNA in liquid and solid fractions of wastewater.

As fragments of SARS-CoV-2 RNA can be quantified and measured temporally in wastewater, surveillance of concentrations of SARS-CoV-2 in wastewater has become a vital resource for tracking the spread of COVID-19 in and among communities. However, the absence of standardized methods has affected the interpretation of data for public health efforts. In particular, analyzing either the liquid or solid fraction has implications for the interpretation of how viral RNA is quantified. Characterizing how SARS-CoV-2 or its RNA fragments partition in wastewater is a central part of understanding fate and behaviour in wastewater. In this study, partitioning of SARS-CoV-2 was investigated by use of centrifugation with varied durations of spin and centrifugal force, polyethylene glycol (PEG) precipitation followed by centrifugation, and ultrafiltration of wastewater. Partitioning of the endogenous pepper mild mottled virus (PMMoV), used to normalize the SARS-CoV-2 signal for fecal load in trend analysis, was also examined. Additionally, two surrogates for coronavirus, human coronavirus 229E and murine hepatitis virus, were analyzed as process controls. Even though SARS-CoV-2 has an affinity for solids, the total RNA copies of SARS-CoV-2 per wastewater sample, after centrifugation (12,000 g, 1.5 h, no brake), were partitioned evenly between the liquid and solid fractions. Centrifugation at greater speeds for longer durations resulted in a shift in partitioning for all viruses toward the solid fraction except for PMMoV, which remained mostly in the liquid fraction. The surrogates more closely reflected the partitioning of SARS-CoV-2 under high centrifugation speed and duration while PMMoV did not. Interestingly, ultrafiltration devices were inconsistent in estimating RNA copies in wastewater, which can influence the interpretation of partitioning. Developing a better understanding of the fate of SARS-CoV-2 in wastewater and creating a foundation of best practices is the key to supporting the current pandemic response and preparing for future potential infectious diseases.

An improved method for determining frequency of multiple variants of SARS-CoV-2 in wastewater using qPCR assays.

Wastewater-based surveillance has become an effective tool around the globe for indirect monitoring of COVID-19 in communities. Variants of Concern (VOCs) have been detected in wastewater by use of reverse transcription polymerase chain reaction (RT-PCR) or whole genome sequencing (WGS). Rapid, reliable RT-PCR assays continue to be needed to determine the relative frequencies of VOCs and sub-lineages in wastewater-based surveillance programs. The presence of multiple mutations in a single region of the N-gene allowed for the design of a single amplicon, multiple probe assay, that can distinguish among several VOCs in wastewater RNA extracts. This approach which multiplexes probes designed to target mutations associated with specific VOC's along with an intra-amplicon universal probe (non-mutated region) was validated in singleplex and multiplex. The prevalence of each mutation (i.e. VOC) is estimated by comparing the abundance of the targeted mutation with a non-mutated and highly conserved region within the same amplicon. This is advantageous for the accurate and rapid estimation of variant frequencies in wastewater. The N200 assay was applied to monitor frequencies of VOCs in wastewater extracts from several communities in Ontario, Canada in near real time from November 28, 2021 to January 4, 2022. This includes the period of the rapid replacement of the Delta variant with the introduction of the Omicron variant in these Ontario communities in early December 2021. The frequency estimates using this assay were highly reflective of clinical WGS estimates for the same communities. This style of qPCR assay, which simultaneously measures signal from a non-mutated comparator probe and multiple mutation-specific probes contained within a single qPCR amplicon, can be applied to future assay development for rapid and accurate estimations of variant frequencies.

Realizing the value in "non-standard" parts of the qPCR standard curve by integrating fundamentals of quantitative microbiology.

The real-time polymerase chain reaction (PCR), commonly known as quantitative PCR (qPCR), is increasingly common in environmental microbiology applications. During the COVID-19 pandemic, qPCR combined with reverse transcription (RT-qPCR) has been used to detect and quantify SARS-CoV-2 in clinical diagnoses and wastewater monitoring of local trends. Estimation of concentrations using qPCR often features a log-linear standard curve model calibrating quantification cycle (Cq) values obtained from underlying fluorescence measurements to standard concentrations. This process works well at high concentrations within a linear dynamic range but has diminishing reliability at low concentrations because it cannot explain "non-standard" data such as Cq values reflecting increasing variability at low concentrations or non-detects that do not yield Cq values at all. Here, fundamental probabilistic modeling concepts from classical quantitative microbiology were integrated into standard curve modeling approaches by reflecting well-understood mechanisms for random error in microbial data. This work showed that data diverging from the log-linear regression model at low concentrations as well as non-detects can be seamlessly integrated into enhanced standard curve analysis. The newly developed model provides improved representation of standard curve data at low concentrations while converging asymptotically upon conventional log-linear regression at high concentrations and adding no fitting parameters. Such modeling facilitates exploration of the effects of various random error mechanisms in experiments generating standard curve data, enables quantification of uncertainty in standard curve parameters, and is an important step toward quantifying uncertainty in qPCR-based concentration estimates. Improving understanding of the random error in qPCR data and standard curve modeling is especially important when low concentrations are of particular interest and inappropriate analysis can unduly affect interpretation, conclusions regarding lab performance, reported concentration estimates, and associated decision-making.

Wild fish responses to wastewater treatment plant upgrades in the Grand River, Ontario.

Municipal wastewater treatment plant (WWTP) effluent is one of several point sources of contaminants (nutrients, pharmaceuticals, estrogens, etc.) which can lead to adverse responses in aquatic life. Studies of WWTP effluent impacts on rainbow darter (Etheostoma caeruleum) collected downstream of WWTPs in the Grand River, Ontario have reported disruption at multiple levels of biological organization, including altered vitellogenin gene expression, lower levels of in vitro steroid production, and high frequency of intersex. However, major upgrades have occurred at treatment plants in the central Grand River over the last decade. Treatment upgrades to the Waterloo WWTP were initiated in 2009 but due to construction delays, the upgrades came fully on-line in 2017/2018. Responses in rainbow darter have been followed at sites associated with the outfall consistently over this entire time period. The treatment plant upgrade resulted in nitrification of effluent, and once complete there was a major reduction in effluent ammonia, selected pharmaceuticals, and estrogenicity. This study compared several key responses in rainbow darter associated with the Waterloo WWTP outfall prior to and post upgrades. Stable isotopes signatures in fish were used to track exposure to effluent and changed dramatically over time, corresponding to the effluent quality. Disruptions in in vitro steroid production and intersex in the darters that had been identified prior to the upgrades were no longer statistically different from the upstream reference sites after the upgrades. Although annual variations in water temperature and flow can potentially mask or exacerbate the effects of the WWTP effluent, major capital investments in wastewater treatment targeted at improving effluent quality have corresponded with the reduction of adverse responses in fish in the receiving environment.

Impacts on antioxidative enzymes and transcripts in darter (Etheostoma spp.) brains in the Grand River exposed to wastewater effluent.

The Grand River watershed is the largest in southern Ontario and assimilates thirty wastewater treatment plants (WWTP) with varied degrees of treatment. Many WWTPs are unable to effectively eliminate several contaminants of emerging concern (CECs) from final effluent, leading to measurable concentrations in surface waters. Exposures to CECs have reported impacts on oxidative stress measured through antioxidative enzymes (SOD, CAT, GPX). This study focuses on the effects of WWTP effluent on four Etheostoma (Darter) species endemic to the Grand River, by investigating if increased antioxidative response markers are present in darter brains downstream from the effluent outfall compared to an upstream reference site relative to the Waterloo, ON WWTP across two separate years (Oct 2020 and Oct 2021). This was assessed using transcriptional and enzyme analysis of antioxidant enzymes and an enzyme involved in serotonin synthesis, tryptophan hydroxylase (tph). In fall 2020, significant differences in transcript markers were found between sites and sexes in GSD with SOD and CAT showing increased expression downstream, in JD with both sexes showing increased SOD downstream, and an interactive effect for tph in RBD. Changes in transcripts aligned with enzyme activity where interactive effects with sex-related differences were observed in fish collected fall 2020. In contrast, transcripts measured in fall 2021 were increased upstream compared to downstream species in RBD and GSD. This study additionally displayed yearly, species and sex differences in antioxidant responses. Continued investigation on the impacts of CECs in effluent in non-target species is required to better understand WWTP effluent impacts.

The relationship between NMDA receptor function and the high ammonia tolerance of anoxia-tolerant goldfish.

Acute ammonia toxicity in vertebrates is thought to be characterized by a cascade of deleterious events resembling those associated with anoxic/ischemic injury in the central nervous system. A key event is the over-stimulation of neuronal N-methyl-D-aspartate (NMDA) receptors, which leads to excitotoxic cell death. The similarity between the responses to acute ammonia toxicity and anoxia suggests that anoxia-tolerant animals such as the goldfish (Carassius auratus Linnaeus) may also be ammonia tolerant. To test this hypothesis, the responses of goldfish were compared with those of the anoxia-sensitive rainbow trout (Oncorhynchus mykiss Walbaum) during exposure to high external ammonia (HEA). Acute toxicity tests revealed that goldfish are ammonia tolerant, with 96 h median lethal concentration (LC(50)) values of 199 µmol l(-1) and 4132 µmol l(-1) for NH(3) and total ammonia ([T(Amm)]=[NH(3)]+[NH(4)(+)]), respectively. These values were ~5-6 times greater than corresponding NH(3) and T(Amm) LC(50) values measured in rainbow trout. Further, the goldfish readily coped with chronic exposure to NH(4)Cl (3-5 mmol l(-1)) for 5 days, despite 6-fold increases in plasma [T] to ~1300 µmol l(-1) and 3-fold increases in brain [T(Amm)] to 6700 µmol l(-1). Muscle [T(Amm)] increased by almost 8-fold from ~900 µmol kg(-1) wet mass (WM) to greater than 7000 µmol kg(-1) WM by 48 h, and stabilized. Although urea excretion rates (J(Urea)) increased by 2-3-fold during HEA, the increases were insufficient to offset the inhibition of ammonia excretion that occurred, and increases in urea were not observed in the brain or muscle. There was a marked increase in brain glutamine concentration at HEA, from ~3000 µmol kg(-1) WM to 15,000 µmol kg(-1) WM after 48 h, which is consistent with the hypothesis that glutamine production is associated with ammonia detoxification. Injection of the NMDA receptor antagonists MK801 (0.5-8 mg kg(-1)) or ethanol (1-8 mg kg(-1)) increased trout survival time by 1.5-2.0-fold during exposure to 2 mmol l(-1) ammonia, suggesting that excitotoxic cell death contributes to ammonia toxicity in this species. In contrast, similar doses of MK801 or ethanol had no effect on ammonia-challenged (8-9.5 mmol l(-1) T(Amm)) goldfish survival times, suggesting that greater resistance to excitotoxic cell death contributes to the high ammonia-tolerance of the goldfish. Whole-cell recordings measured in isolated brain slices of goldfish telencephalon during in vitro exposure to 5 mmol l(-1) or 10 mmol l(-1) T(Amm) reversibly potentiated NMDA receptor currents. This observation suggested that goldfish neurons may not be completely resistant to ammonia-induced excitotoxicity. Subsequent western blot and densitometric analyses revealed that NMDA receptor NR1 subunit abundance was 40-60% lower in goldfish exposed to 3-5 mmol l(-1) T(Amm) for 5 days, which was followed by a restoration of NR1 subunit abundance after 3 days recovery in ammonia-free water. We conclude that the goldfish brain may be protected from excitotoxicity by downregulating the abundance of functional NMDA receptors during periods when it experiences increased internal ammonia.