Transcriptome signatures of wastewater effluent exposure in larval zebrafish vary with seasonal mixture composition in an effluent-dominated stream

Wastewater treatment plant (WWTP) effluent-dominated streams provide critical habitat for aquatic and terrestrial or-ganisms but also continually expose them to complex mixtures of pharmaceuticals that can potentially impair growth, behavior, and reproduction. Currently, few biomarkers are available that relate to pharmaceutical-specific mecha-nisms of action. In the experiment reported in this paper, zebrafish (Danio rerio) embryos at two developmental stages were exposed to water samples from three sampling sites (0.1 km upstream of the outfall, at the effluent outfall, and 0.1 km below the outfall) during base-flow conditions from two months (January and May) of a temperate-region effluent-dominated stream containing a complex mixture of pharmaceuticals and other contaminants of emerging con-cern. RNA-sequencing identified potential biological impacts and biomarkers of WWTP effluent exposure that extend past traditional markers of endocrine disruption. Transcriptomics revealed changes to a wide range of biological func-tions and pathways including cardiac, neurological, visual, metabolic, and signaling pathways. These transcriptomic changes varied by developmental stage and displayed sensitivity to variable chemical composition and concentration of effluent, thus indicating a need for stage-specific biomarkers. Some transcripts are known to be associated with genes related to pharmaceuticals that were present in the collected samples. Although traditional biomarkers of endo-crine disruption were not enriched in either month, a high estrogenicity signal was detected upstream in May and im-plicates the presence of unidentified chemical inputs not captured by the targeted chemical analysis. This work reveals associations between bioeffects of exposure, stage of development, and the composition of chemical mixtures in effluent-dominated surface water. The work underscores the importance of measuring effects beyond the endocrine system when assessing the impact of bioactive chemicals in WWTP effluent and identifies a need for non-targeted chemical analysis when bioeffects are not explained by the targeted analysis.