Background The spread of SARS-CoV-2 has been studied at unprecedented levels worldwide. In jurisdictions where molecular analysis was performed on large scales, the emergence and competition of numerous SARS-CoV-2lineages have been observed in near real-time. Lineage identification, traditionally performed from clinical samples, can also be determined by sampling wastewater from sewersheds serving populations of interest. Variants of concern (VOCs) and SARS-CoV-2 lineages associated with increased transmissibility and/or severity are of particular interest. Method Here, we consider clinical and wastewater data sources to assess the emergence and spread of VOCs in Canada retrospectively. Results We show that, overall, wastewater-based VOC identification provides similar insights to the surveillance based on clinical samples. Based on clinical data, we observed synchrony in VOC introduction as well as similar emergence speeds across most Canadian provinces despite the large geographical size of the country and differences in provincial public health measures. Conclusion In particular, it took approximately four months for VOC Alpha and Delta to contribute to half of the incidence. In contrast, VOC Omicron achieved the same contribution in less than one month. This study provides significant benchmarks to enhance planning for future VOCs, and to some extent for future pandemics caused by other pathogens, by quantifying the rate of SARS-CoV-2 VOCs invasion in Canada.

The search for better tools for interpreting and understanding wastewater surveillance has continued since the beginning of the coronavirus disease 2019 (COVID-19) pandemic. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued to mutate, thus complicating the interpretation of surveillance results. We assessed the Omicron variants (BA.1, BA.2, and BA.5) associated with wastewater-derived SARS-CoV-2 RNA trends by estimating the effective reproduction number (Reff) using an epidemic model that integrates explicitly the SARS-CoV-2 N2 gene concentration detected in wastewater through rt-qPCR quantitative analysis. The model inferred COVID-19 cases based on wastewater data and compared them with the ones reported by clinical surveillance. The variant of the SARS-CoV-2 associated with the wastewater-derived viral RNA was monitored through wastewater whole-genome sequencing. Three major waves between January and September 2022 were associated with the Omicron subvariants (BA.1, BA.2, and BA.5). This work showed that disease trends can be monitored using estimates of the effective reproduction number which is simple and easy to understand.

Distributions of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and fecal viral biomarkers between solid and liquid phases of wastewater are largely unknown. Herein, distributions of SARS-CoV-2, Pepper Mild Mottle Virus (PMMoV), and F-RNA bacteriophage group II (FRNAPH-II) were determined by viral RNA RT-qPCR. Comparison of viral recovery using three conventional fractionation methods included membrane filtration, a combination of mid-speed centrifugation and membrane filtration, and high-speed centrifugation. SARS-CoV-2 partitioned to the solids fraction in greater abundance compared to liquid fractions in a combination of mid-speed centrifugation and membrane filtration and high-speed centrifugation, but not in membrane filtration method in a particular assay, while fecal biomarkers (PMMoV and FRNAPH-II) exhibited the reciprocal relationship. The wastewater fractionation method had minimal effects on the solids-liquids distribution for all viral and phage markers tested; however, viral RNA load was significantly greater in solid–liquid fractions viral RNA loads compared with the than whole-wastewater PEG precipitation. A RNeasy PowerWater Kit with PCR inhibitor removal resulted in greater viral RNA loads and lesser PCR inhibition compared to a QIAamp Viral RNA Mini Kit without PCR inhibitor removal. These results support the development of improved methods and interpretation of WBE of SARS-CoV-2.•Distribution of SARS-CoV-2 to liquid and solid portions was addressed.•Addressing PCR inhibition is important in wastewater-based epidemiology.•Fraction methods have minimal effect.

Conventional wastewater surveillance (WS) relies on highly trained personnel, advanced instrumentation, and significant resources, making the development and use of simple, rapid, and sensitive alternative technologies valuable for reducing costs, time, and labor intensity. For the first time, this study investigated the use of two well-developed rapid systems, including the GeneXpert and LuminUltra, in parallel with a conventional WS reference methodology for the assessment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in three cities: Saskatoon, Prince Albert, and North Battleford, Saskatchewan, Canada. RNA extractions from wastewater samples were carried out for the conventional reference and LuminUltra methods, while GeneXpert was used for both raw and concentrated wastewater samples. Bland-Altman plots showed a combination of systematic bias and random error between these real-time reverse transcription-quantitative polymerase chain reaction (RT-qPCR)-based systems. Additionally, results indicated the reasonable performance of GeneXpert in viral detection with a sensitivity rate of >98%, as compared to the conventional reference methodology of 100% and LuminUltra with >65%. A Spearman correlation test showed meaningful relationships between the GeneXpert and conventional reference methodology viral level results across all cities, indicating GeneXpert's reliability for accurate viral detection and disease prevalence determination, specifically in limited-resource communities, with a shorter processing time and cost-effectiveness in analysis.

Wastewater surveillance (WS) helps to improve the understanding of the spread of communicable diseases in communities. WS can assist public health decision-makers in the design and implementation of timely mitigation measures. There is an increased need to use reliable, cost-effective, simple, and rapid WS systems, given traditional analytical (or 'gold-standard') programs are instrument/time-intensive, and dependent on highly skilled personnel. This study investigated the application of the portable GeneXpert platform for WS of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), influenza A virus (IAV), influenza B virus (IBV), and respiratory syncytial virus (RSV). The GeneXpert system with the Xpert Xpress-SARS-CoV-2/Flu/RSV test kit uses reverse transcription-quantitative polymerase chain reaction (RT-qPCR) to analyze wastewater samples. From September 2022 through January 2023, wastewater samples were collected from the influents of municipal wastewater treatment plants (MWTPs) of Saskatoon, Prince Albert, and North Battleford in the province of Saskatchewan, Canada. Both raw and concentrated wastewater samples were subjected to the GeneXpert analysis. Results showed that the Saskatoon wastewater viral loads were significantly correlated to Saskatchewan's influenza and COVID-19 clinical cases, with a lead time of 10 days for IAV and a lag time of 4 days for SARS-CoV-2. Additionally, the GeneXpert analysis of the three cities' wastewater samples showed that the raw WS could capture the dynamics of SARS-CoV-2 and IAV due to their correlation with concentrated WS. Interestingly, IBV loads were not detected in any wastewater samples, while the Saskatoon and Prince Albert wastewater samples collected following the 2023 holiday season (end of December and beginning of January) were positive for RSV. This study indicates that the GeneXpert has excellent potential for use in the development of an early warning system for transmissible disease in municipalities and limited-resource communities while simultaneously providing stakeholders with an efficient WS methodology.

This study investigated passive adsorption of nine antipsychotic drugs by the benthic oligochaete, Lumbriculus variegatus from spiked freshwater sediments. Based on a combination of diffusive gradients in thin-films (DGT) passive sampling and modelling, specifically the DGT-induced fluxes in sediments (DIFS) model and fractions-transfer models, an uptake and diffusion-induced model was established for oligochaetes. Lamotrigine and bupropion showed positive fluxes from sediment to overlying water. Oligochaetes accumulated both lamotrigine and carbamazepine. The simulation model for transfer showed that after 30 days antipsychotic compounds in sediment were present predominantly in stable or bound fractions, rather than labile fractions. Dissolved organic carbon and degradation were important drivers of dissipation rates of antipsychotic compounds in sediments and labile pool size. Uptake by oligochaetes and diffusion-induced models matched measured concentrations well (R2= 0.83–0.95 and R2= 0.99–1.00, respectively). Processes of uptake by oligochaetes were comparable to DIFS results, where fast desorption to the labile pool in the short response time and adsorption processes, except carbamazepine and lamotrigine, in sediments dominated kinetics. It suggests a link between bioavailability of organic pollutants for sediment-dwelling organisms and desorption kinetics involved with fractions transfer within sediments, which has important implications for the environmental risk assessment of these chemicals.

Perfluoroethylcyclohexanesulfonate (PFECHS) is an emerging perfluoroalkyl substance (PFAS) that has been considered a potential replacement for perfluorooctanesulfonic acid (PFOS). However, there is little information characterizing the toxic potency of PFECHS to zebrafish embryos and its potential for effects in aquatic environments. This study assessed toxic potency of PFECHS in vivo during both acute (96-hour postfertilization) and chronic (21-day posthatch) exposures and tested concentrations of PFECHS from 500 ng/L to 2 mg/L. PFECHS was less likely to cause mortalities than PFOS for both the acute and chronic experiments based on previously published values for PFOS exposure, but exposure resulted in a similar incidence of deformities. Exposure to PFECHS also resulted in significantly increased abundance of transcripts of peroxisome proliferator activated receptor alpha (pparα), cytochrome p450 1a1 (cyp1a1), and apolipoprotein IV (apoaIV) at concentrations nearing those of environmental relevance. Overall, these results provide further insight into the safety of an emerging PFAS alternative in the aquatic environment and raise awareness that previously considered "safer" alternatives may show similar effects as legacy PFASs.

Perfluoroethylcyclohexane sulphonate (PFECHS) is an emerging, replacement perfluoroalkyl substance (PFAS) with little information available on the toxic effects or potencies with which to characterize its potential impacts on aquatic environments. This study aimed to characterize effects of PFECHS using in vitro systems, including rainbow trout liver cells (RTL-W1 cell line) and lymphocytes separated from whole blood. It was determined that exposure to PFECHS caused minor acute toxic effects for most endpoints and that little PFECHS was concentrated into cells with a mean in vitro bioconcentration factor of 81 +/- 25 L/kg. However, PFECHS was observed to affect the mitochondrial membrane and key molecular receptors, such as the peroxisome proliferator receptor, cytochrome p450-dependent monooxygenases, and receptors involved in oxidative stress. Also, glutathione-Stransferase was significantly down-regulated at a near environmentally relevant exposure concentration of 400 ng/L. These results are the first to report bioconcentration of PFECHS, as well as its effects on the peroxisome proliferator and glutathione-S-transferase receptors, suggesting that even with little bioconcentration, PFECHS has potential to cause adverse effects.

N-(1,3-Dimethylbutyl)-N '-phenyl-p-phenylenediamine-quinone (6PPD-Q) is a recently identified contaminant that originates from the oxidation of the tire antidegradant 6PPD. 6PPD-Q is acutely toxic to select salmonids at environmentally relevant concentrations, while other fish species display tolerance to concentrations that surpass those measured in the environment. The reasons for these marked differences in sensitivity are presently unknown. The objective of this research was to explore potential toxicokinetic drivers of species sensitivity by characterizing biliary metabolites of 6PPD-Q in sensitive and tolerant fishes. For the first time, we identified an O-glucuronide metabolite of 6PPD-Q using high-resolution mass spectrometry. The semiquantified levels of this metabolite in tolerant species or life stages, including white sturgeon (Acipenser transmontanus), chinook salmon (Oncorhynchus tshawytscha), westslope cutthroat trout (Oncorhynchus clarkii lewisi), and nonfry life stages of Atlantic salmon (Salmo salar), were greater than those in sensitive species, including coho salmon (Oncorhynchus kisutch), brook trout (Salvelinus fontinalis), and rainbow trout (Oncorhynchus mykiss), suggesting that tolerant species might detoxify 6PPD-Q more effectively. Thus, we hypothesize that differences in species sensitivity are a result of differences in basal expression of biotransformation enzyme across various fish species. Moreover, the semiquantification of 6PPD-Q metabolites in bile extracted from wild-caught fish might be a useful biomarker of exposure to 6PPD-Q, thereby being valuable to environmental monitoring and risk assessment.

Wastewater surveillance (WWS) is useful to better understand the spreading of coronavirus disease 2019 (COVID-19) in communities, which can help design and implement suitable mitigation measures. The main objective of this study was to develop the Wastewater Viral Load Risk Index (WWVLRI) for three Saskatchewan cities to offer a simple metric to interpret WWS. The index was developed by considering relationships between reproduction number, clinical data, daily per capita concentrations of virus particles in wastewater, and weekly viral load change rate. Trends of daily per capita concentrations of SARS-CoV-2 in wastewater for Saskatoon, Prince Albert, and North Battleford were similar during the pandemic, suggesting that per capita viral load can be useful to quantitatively compare wastewater signals among cities and develop an effective and comprehensible WWVLRI. The effective reproduction number (Rt) and the daily per capita efficiency adjusted viral load thresholds of 85 × 106 and 200 × 106 N2 gene counts (gc)/population day (pd) were determined. These values with rates of change were used to categorize the potential for COVID-19 outbreaks and subsequent declines. The weekly average was considered 'low risk' when the per capita viral load was 85 × 106 N2 gc/pd. A 'medium risk' occurs when the per capita copies were between 85 × 106 and 200 × 106 N2 gc/pd. with a rate of change <100 %. The start of an outbreak is indicated by a 'medium-high' risk classification when the week-over-week rate of change was >100 %, and the absolute magnitude of concentrations of viral particles was >85 × 106 N2 gc/pd. Lastly, a 'high risk' occurs when the viral load exceeds 200 × 106 N2 gc/pd. This methodology provides a valuable resource for decision-makers and health authorities, specifically given the limitation of COVID-19 surveillance based on clinical data.

ABSTRACT N-(1,3-Dimethylbutyl)-N’-Phenyl-P-Phenylenediamine-Quinone (6PPD-Q) is a recently identified contaminant that originates from the oxidation of the tire anti-degradant 6PPD. 6PPD-Q is acutely toxic to select salmonids at environmentally relevant concentrations, while other fish species display tolerance to concentrations surpassing those measured in the environment. The reasons for these marked differences in sensitivity are presently unknown. The objective of this research was to explore potential toxicokinetic drivers of species sensitivity by characterizing biliary metabolites of 6PPD-Q in sensitive and tolerant fishes. For the first time, we identified an O -glucuronide metabolite of 6PPD-Q using high-resolution mass spectrometry. The semi-quantified levels of this metabolite in tolerant species or life stages, including white sturgeon ( Acipenser transmontanus ), chinook salmon ( Oncorhynchus tshawytscha ), westslope cutthroat trout ( Oncorhynchus clarkia lewisi ) and non-fry life stages of Atlantic salmon ( Salmo salar ), were greater than those in sensitive species, including coho salmon ( Oncorhynchus kisutch ), brook trout ( Salvelinus fontinalis ), and rainbow trout ( Oncorhynchus mykiss ), suggesting that tolerant species might more effectively detoxify 6PPD-Q. Thus, we hypothesize that differences in species sensitivity are a result of differences in basal expression of biotransformation enzyme across various fish species. Moreover, the semi-quantification of 6PPD-Q metabolites in bile extracted from wild-caught fish might be a useful biomarker of exposure to 6PPD-Q, thereby being invaluable to environmental monitoring and risk assessment.

In vitro biotransformation assays with primary trout hepatocytes (RT-HEP) or liver subcellular fractions (RT-S9) have been proposed as valuable tools to help scientists and regulators better understand the toxicokinetics of chemicals. While both assays have been applied successfully to a diversity of neutral organic chemicals, only the RT-S9 assay has been applied to a large number of ionizable organic chemicals. Here, a combination of an in vitro biotransformation assay with RT-HEP with an active transport assay based on the permanent rainbow trout liver cell line RTL-W1 was used to qualitatively predict the potential hepatic clearance of nine psychotropic drugs with various degrees of ionization. Predictions were compared with rates of clearance measured in isolated perfused rainbow trout livers, and the importance of active transport was verified in the presence of the active transport inhibitor cyclosporin A. For the first time, it was demonstrated that a combination of biotransformation and active transport assays is powerful for the prediction of rates of hepatic clearance of ionizable chemicals. Ultimately, it is expected that this approach will allow for use of fewer animals while at the same time improving our confidence in the use of data from in vitro assays in chemical risk assessment.

Forms of organic contaminants is an important driver of bioavailable fraction and desorption kinetics of pollutants binding to sediments. To determine fluxes and resupply of nine environmentally-relevant antipsychotic drugs, which are emerging pollutants that can have adverse effects on aquatic organisms, interface passive samplers of diffusive gradients in thin films (DGT) were deployed for 21 days, in situ at the sediment-water interface in submerged sandy riverbank sediments. At each deployment time, samples of sediment were collected and subjected to consecutive extraction of pore water, as well as rapidly-desorbing (labile), stabledesorbing, and bound residue fractions. Concentrations of antipsychotic drugs decreased with sediment depth with the greatest concentrations observed in the top 2 cm. Positive fluxes of antipsychotic drugs were observed from sediment to surface water. The dynamic fraction transfer model indicated that the labile fraction can be resupplied with a lag time (> 21 d). When results were further interpreted using the DGT-induced fluxes in soils and sediments (DIFS) model, partial resupply of antipsychotic drugs from sediment particles to porewater was demonstrated. Desorption occurred within the entirety of the observed 15 cm depth of sediment. Fastest rates of resupply were found for carbamazepine and lamotrigine. Size of the labile pool estimated by the DIFS model did not fully explain the observed resupply, while a first-order three-compartment kinetic model for the fast-desorbing fraction can be used to supplement DIFS predictions with estimations of labile pool size.

Monitoring the communal incidence of COVID-19 is important for both government and residents of an area to make informed decisions. However, continuous reliance on one means of monitoring might not be accurate because of biases introduced by government policies or behaviours of residents. Wastewater surveillance was employed to monitor concentrations of SARS-CoV-2 RNA in raw influent wastewater from wastewater treatment plants serving three Canadian Prairie cities with different population sizes. Data obtained from wastewater are not directly influenced by government regulations or behaviours of individuals. The means of three weekly samples collected using 24 h composite auto-samplers were determined. Viral loads were determined by RT-qPCR, and whole-genome sequencing was used to charaterize variants of concern (VOC). The dominant VOCs in the three cities were the same but with different proportions of sub-lineages. Sub-lineages of Delta were AY.12, AY.25, AY.27 and AY.93 in 2021, while the major sub-lineage of Omicron was BA.1 in January 2022, and BA.2 subsequently became a trace-level sub-variant then the predominant VOC. When each VOC was first detected varied among cities; However, Saskatoon, with the largest population, was always the first to present new VOCs. Viral loads varied among cities, but there was no direct correlation with population size, possibly because of differences in flow regimes. Population is one of the factors that affects trends in onset and development of local outbreaks during the pandemic. This might be due to demography or the fact that larger populations had greater potential for inter- and intra-country migration. Hence, wastewater surveillance data from larger cities can typically be used to indicate what to expect in smaller communities.

Highlights• Viral load of SARS-CoV-2 RNA was measured in three wastewater treatment plants serving communities with different population sizes• Whole-genome sequencing of RNA from wastewater revealed AY.12, AY.25, AY.27, AY.93 as common AY-lineages• Omicron became the dominant VOC within a short time in the three cities• Traces of Omicron-sub-lineage BA.2 were detected through sequencing at the sub-consensus level.

There are no standardized protocols for quantifying severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in wastewater to date, especially for population normalization. Here, a pipeline was developed, applied, and assessed to quantify SARS-CoV-2 and key variants of concern (VOCs) RNA in wastewater at Saskatoon, Canada. Normalization approaches using recovery ratio and extraction efficiency, wastewater parameters, or population indicators were assessed by comparing to daily numbers of new cases. Viral load was positively correlated with daily new cases reported in the sewershed. Wastewater surveillance (WS) had a lead time of approximately 7 days, which indicated surges in the number of new cases. WS revealed the variant α and δ driving the third and fourth wave, respectively. The adjustment with the recovery ratio and extraction efficiency improved the correlation between viral load and daily new cases. Normalization of viral concentration to concentrations of the artificial sweetener acesulfame K improved the trend of viral load during the Christmas and New Year holidays when populations were dynamic and variable. Acesulfame K performed better than pepper mild mottle virus, creatinine, and ammonia for population normalization. Hence, quality controls to characterize recovery ratios and extraction efficiencies and population normalization with acesulfame are promising for precise WS programs supporting decision-making in public health.

Dynamic processes of organic contaminants in sediments can have important toxicological implications in aquatic systems. The current study used diffusive gradients in thin-films (DGT) devices in sandy sediments spiked with nine antipsychotics and in field sandy sediments. Samplers were deployed for 1 to 30 days to determine the flux of these compounds to DGT devices and the exchange rates between the porewater and sediment solid phase. The results showed a continuous removal of antipsychotics to a binding gel and induced a mobile flux from the DGT device to the adjacent sediment solution. A dynamic model, DGT-induced fluxes in soils and sediments, was used to derive rate constants of resupply of antipsychotics from solid phase to aqueous phase (response time, Tc) and distribution coefficients for labile antipsychotics. The largest labile pool was found for lamotrigine and carbamazepine in spiked sediments. Carbamazepine, clozapine, citalopram, and lamotrigine were resupplied rapidly by sediments with Tc (25-30 min). Tc values of bupropion and amitriptyline were the longest (≈5 h), which exhibited slow desorption rates in sediments. In field sediments, high resupply was found for carbamazepine and lamotrigine, which did not show higher labile pool. The Tc values were obviously higher in the filed sediments (52-171 h). Although the adsorption process is dominant for most studied antipsychotics in both spiked sediments and field sediments, the kinetic resupply of antipsychotic compounds may not be accurately estimated by laboratory-controlled incubation experiments. More studies are needed to explore the mechanisms of desorption kinetics by using in situ DGT technique in the field.

Previously unreported natural and synthetic organo-bromine compounds (NSOBCs) have been found to contribute more than 99% of total organic bromine (TOB) in environmental matrices. We recently developed a novel untargeted method (data-independent precursor isolation and characteristic fragment, DIPIC-Frag) and identified ∼2000 NSOBCs in two sediments from Lake Michigan. In this study, this method was used to investigate the distributions of these NSOBCs in 23 surficial samples and 24 segments of a sediment core from Lake Michigan. NSOBCs were detected in all 23 surficial samples and exhibited 10- to 100-fold variations in peak abundance among locations. The pattern of distributions of NSOBCs was correlated with depth of the water column (r2 = 0.61, p < 0.001). Hierarchical cluster analysis showed that sediments in close proximity exhibited similar profiles of NSOBCs. Distributions of NSOBCs in 24 segments of a sediment core dated from 1766 to 2008 were investigated, and samples from similar depths exhibited similar profiles of NSOBCs. NSOBCs were grouped into four clusters (soft-cluster analysis) with different temporal trends of abundances. 515 and 768 of the NSOBCs were grouped into cluster 1 and cluster 3 with increasing temporal trends, especially since 1950, indicating that abundances of these compounds might have been affected by human activities.

The majority of halogenated organic compounds present in the environment remain unidentified. To address this data gap, we recently developed an untargeted method (data-independent precursor isolation and characteristic fragment; DIPIC-Frag) for identification of unknown organo-bromine compounds. In this study, the method was adapted to enable untargeted screening of natural and synthetic organo-iodine compounds (NSOICs) in sediments. A total of 4,238 NSOIC peaks were detected in sediments from Lake Michigan. Precursor ions and formulas were determined for 2,991 (71%) of the NSOIC peaks. These compounds exhibited variations in abundances (<10(3) to ∼10(7)), m/z values (206.9304-996.9474), retention times (1.0-29.7 min), and number of iodine atoms (1-4). Hierarchical cluster analysis showed that sediments in closer proximity exhibited similar profiles of NSOICs. NSOICs were screened in 10 samples of sediment from the Arctic Ocean to compare the profiles of NSOICs between freshwater and marine sediments. A total of 3,168 NSOIC peaks were detected, and profiles of NSOICs in marine sediments were clearly distinct from Lake Michigan. The coexistence of brominated and iodinated analogues indicated that some NSOICs are of natural origin. Different ratios of abundances of iodinated compounds to brominated analogues were observed and proposed as a marker to distinguish sources of NSOICs.