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.

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.

Sex matters: Gamete-specific contribution of microRNA following parental exposure to hypoxia in zebrafish.

Oxygen availability varies among aquatic environments, and oxygen concentration has been demonstrated to drive behavioral, metabolic, and genetic adaptations in numerous aquatic species. MicroRNAs (miRNAs) are epigenetic modulators that act at the interface of the environment and the transcriptome and are known to drive plastic responses following environmental stressors. An area of miRNA that has remained underexplored is the sex specific action of miRNAs following hypoxia exposure and its effects as gene expression regulator in fishes. This study aimed to identify differences in mRNA and miRNA expression in the F1 generation of zebrafish (Danio rerio) at 1 hpf after either F0 parental male or female were exposed to 2 weeks of continuous (45 %) hypoxia. In general, F1 embryos at 1 hpf demonstrated differences in mRNA and miRNAs expression related to the stressor and to the specific sex of the F0 that was exposed to hypoxia. Bioinformatic pathway analysis of predicted miRNA:mRNA relationships indicated responses in known hypoxia signaling and mitochondrial bioenergetic pathways. This research demonstrates the importance of examining the specific male and female contributions to phenotypic variation in subsequent generations and provides evidence that there is both maternal and paternal contribution of miRNA through eggs and sperm.

Responses of microRNA and predicted mRNA and enzymatic targets in liver of two salmonids (Oncorhynchus mykiss and Salvelinus fontinalis) following air exposure.

The acute stress response is well-characterized, with rainbow trout as a teleost model for physiological and molecular responses. Air exposure, which stimulates an acute stress response, modulates liver microRNAs in rainbow trout; however, these highly conserved non-coding RNAs that bind to mRNA and repress translation, have never been measured in brook trout and it is unknown how miRNA expression responds following air exposure in this less studied salmonid. Our objective was to characterize the effects of air exposure on rainbow and brook trout liver miRNA expression, as well as the mRNA expression and enzyme activity that the miRNAs are predicted to target. Brook and rainbow trout were sampled pre- and 1-, 3-, and 24-h post- a three-minute air exposure. Plasma cortisol, glucose, and lactate were measured. Relative expression of miR-21a-5p, miR-143-3p, let-7a-5p and relative expression and enzyme activities of five predicted targets (pyruvate kinase, glucokinase, citrate synthase, cytochrome c oxidase, and catalase) were measured in liver. Rainbow and brook trout both had increases in plasma cortisol and lactate, while only rainbow trout had significant post-stress increases in plasma glucose. Furthermore, both trout species had increased miR-143-3p and miR-21a-5p relative expression 24-h post-stress. Four of the five enzymes measured had altered activity following stress. Brook trout miRNAs had inverse relative expression with relative catalase mRNA expression and cytochrome c oxidase enzyme activity, but no relationship was found in rainbow trout. Therefore, we have further characterized the transcriptional and enzymatic response to air exposure in two salmonids.

Chronic exposure to anthropogenic and climate related stressors alters transcriptional responses in the liver of zebrafish (Danio rerio) across multiple generations.

The antidepressant, venlafaxine (VFX), and climate change stressors, such as increased water temperature and decreased dissolved oxygen, are current threats to aquatic environments. This study aimed to determine how microRNAs (miRNAs) and predicted targeted transcripts were altered in livers of zebrafish exposed to these stressors, and livers of their un-exposed F1 and F2 offspring. Following a 21 day exposure to multiple stressors (1 µg/L VFX, +5 °C ambient, 50% O2), then a subsequent 21 day recovery, relative abundances of cyp3a65, hsp70, hsp90, and ppargc1a and miRNAs predicted to target them (miR-142a, miR-16c, miR-181c, and miR-129, respectively) were measured in the liver via quantitative PCR (RT-qPCR). There were significant decreases in miR-142a in the exposed F0 generation and the exposed F1 generation. While there were no changes detected in cyp3a65 relative abundance, there was a significant inverse relationship between cyp3a65 and miR-142a. Hsp70 expression significantly increased in the F1 generation, which persisted to the F2 generation and the relative abundance of hsp90 significantly increased in all generations. There was a significant reduction in miR-181c in the F1 generation, but there was no significant relationship between miR-181c and hsp90. Finally, there was a significant decrease in ppargc1a relative abundance in the F1 generation which was associated with an increase in miR-129. Combined, these results suggest that parental exposure to multiple, environmentally relevant stressors can confer transcriptional and epigenetic responses in the F1 and F2 generations, although identifying which stressor is a driving force becomes unclear.

Acute air exposure modulates the microRNA abundance in stress responsive tissues and circulating extracellular vesicles in rainbow trout (Oncorhynchus mykiss).

The hypothalamic-pituitary-interrenal axis is an important regulator of stress and metabolism in teleosts. Cortisol is secreted by the head kidney where it increases gluconeogenesis in the liver to increase circulating glucose levels. MicroRNAs (miRNAs) are small, non-coding RNA molecules that bind to the 3' untranslated region of specific mRNA to regulate their expression. MicroRNAs can also be secreted into circulation by association with extracellular vesicles (EVs) where they can influence the phenotype of other tissues. In this study, adult rainbow trout were exposed to a 3-minute acute air stress and allowed to recover for 1-, 3-, or 24-h to determine how miRNAs were altered. MicroRNAs measured in this study were chosen based on their high relative abundance in tissues that drive the stress response (miR-21a-3p, let-7a-5p, miR-143-3p) or their role in regulating DNA methylation (miR-29a-3p). In general, miRNAs increased in circulating EVs during the recovery period while decreasing in head kidney and liver at the same timepoints. Predicted targets for these miRNAs were analyzed using KEGG and DAVID functional enrichment analysis. Pathways involved in metabolism and cell signaling were predicted to be upregulated. Future studies can use these results to investigate how pathways are regulated after stress. Overall, our results indicate that miRNAs are regulated during teleost stress responses and could be supporting the cortisol-mediated changes that occur.

High Throughput Sequencing of MicroRNA in Rainbow Trout Plasma, Mucus, and Surrounding Water Following Acute Stress.

Circulating plasma microRNAs (miRNAs) are well established as biomarkers of several diseases in humans and have recently been used as indicators of environmental exposures in fish. However, the role of plasma miRNAs in regulating acute stress responses in fish is largely unknown. Tissue and plasma miRNAs have recently been associated with excreted miRNAs; however, external miRNAs have never been measured in fish. The objective of this study was to identify the altered plasma miRNAs in response to acute stress in rainbow trout (Oncorhynchus mykiss), as well as altered miRNAs in fish epidermal mucus and the surrounding ambient water. Small RNA was extracted and sequenced from plasma, mucus, and water collected from rainbow trout pre- and 1 h-post a 3-min air stressor. Following small RNA-Seq and pathway analysis, we identified differentially expressed plasma miRNAs that targeted biosynthetic, degradation, and metabolic pathways. We successfully isolated miRNA from trout mucus and the surrounding water and detected differences in miRNA expression 1-h post air stress. The expressed miRNA profiles in mucus and water were different from the altered plasma miRNA profile, which indicated that the plasma miRNA response was not associated with or immediately reflected in external samples, which was further validated through qPCR. This research expands understanding of the role of plasma miRNA in the acute stress response of fish and is the first report of successful isolation and profiling of miRNA from fish mucus or samples of ambient water. Measurements of miRNA from plasma, mucus, or water can be further studied and have potential to be applied as non-lethal indicators of acute stress in fish.

Chronic exposure to venlafaxine and increased water temperature reversibly alters microRNA in zebrafish gonads (Danio rerio).

MicroRNA (miRNA) are short, non-coding RNA that act by downregulating targeted mRNA transcripts. Only recently have they been used as endpoints in studies of aquatic toxicology. The purpose of this study was to determine the effect of an antidepressant contaminant, venlafaxine (VFX), and increased temperature on specific microRNA levels in zebrafish (Danio rerio) reproductive tissue. Adult zebrafish were exposed to one of four conditions; control, 1 µg/L VFX (VFX), 32 °C (Temp), or 1 µg/L VFX + 32 °C (VFX & Temp) for 21 days. Half of the fish were returned to control conditions for a 21-day recovery period. RT-qPCR was performed to measure relative abundances of several miRNAs known to respond to antidepressant exposure: dre-miR-22b-3p, dre-miR-301a, dre-miR-140-5p, dre-let-7d-5p, dre-miR-210-5p, and dre-miR-457b-5p. After the exposure period, dre-miR-22b-3p and dre-miR-301a showed a significant downregulation in response to all treatments. In contrast, after the recovery period, there were no significant differences in microRNA abundance. These altered microRNA are predicted to target several genes, including phosphofructokinase, and are associated with ovarian pathologies. Combined, we have shown that VFX and increased water temperature alter miRNA abundances in zebrafish reproductive tissue, an effect correlated with a functional stress response and cell cycle dysregulation.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes.

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.