Omics for aquatic ecotoxicology: control of extraneous variability to enhance the analysis of environmental effects

Evidence of interregional similarity in crayfish metabolomes at reference sites: Progress towards the metabolome as a biomonitoring tool

It has been proposed that biomonitoring may benefit from the use of metabolomics (the study of all small molecules in an organism) to detect sub-lethal organism stress through changes in the metabolite profile (i.e., the metabolome). However, to integrate the metabolome into biomonitoring programs the amount of natural variability among and within populations of indicator taxa must be established prior to generating a reference condition. This study determined variation in the metabolome among ecoregion and stream of origin in the northern crayfish (Faxonius virilis) and if that variation inhibited detection of stressor effects at sites exposed to human activities. We collected crayfish from seven minimally disturbed streams (i.e., reference streams), distributed across three level II ecoregions in central Canada and compared their metabolomes. We found ecoregion and stream origin were poor predictors of crayfish metabolomes. This result suggests crayfish metabolomes were similar, despite differing environmental conditions. Metabolomes of crayfish collected from three stream sites exposed to agricultural activity and municipal wastewater (i.e., test sites) were then compared to the crayfish metabolomes from the seven reference streams. Findings showed that crayfish metabolomes from test sites were strongly differentiated from those at all reference sites. The consistency in the northern crayfish metabolome at the studied reference streams indicates that a single reference condition may effectively detect impacts of human activities across the sampled ecoregions.

Insights on the Organ-Dependent, Molecular Sexual Dimorphism in the Zebra Mussel, Dreissena polymorpha, Revealed by Ultra-High-Performance Liquid Chromatography-Tandem Mass Spectrometry Metabolomics

The zebra mussel, Dreissena polymorpha, is extensively used as a sentinel species for biosurveys of environmental contaminants in freshwater ecosystems and for ecotoxicological studies. However, its metabolome remains poorly understood, particularly in light of the potential molecular sexual dimorphism between its different tissues. From an ecotoxicological point of view, inter-sex and inter-organ differences in the metabolome suggest variability in responsiveness, which can influence the analysis and interpretation of data, particularly in the case where males and females would be analyzed indifferently. This study aimed to assess the extent to which the molecular fingerprints of functionally diverse tissues like the digestive glands, gonads, gills, and mantle of D. polymorpha can reveal tissue-specific molecular sexual dimorphism. We employed a non-targeted metabolomic approach using liquid chromatography high-resolution mass spectrometry and revealed a significant sexual molecular dimorphism in the gonads, and to a lesser extent in the digestive glands, of D. polymorpha. Our results highlight the critical need to consider inter-sex differences in the metabolome of D. polymorpha to avoid confounding factors, particularly when investigating environmental effects on molecular regulation in the gonads, and to a lesser extent in the digestive glands.

Influence of chemical dose and exposure duration on protein synthesis in green sea turtle primary cells

Understanding the impacts of chemical exposure in marine wildlife is challenging, due to practical and ethical constraints that preclude traditional toxicology research on these animals. This study addressed some of these limitations by presenting an ethical and high throughput cell-based approach to elucidate molecular-level effects of contaminants on sea turtles. The experimental design addressed basic questions of cell-based toxicology, including chemical dose and exposure time. Primary green turtle skin cells were exposed to polychlorinated biphenyl (PCB) 153 and perfluorononanoic acid (PFNA) for 24 and 48 h, at three sub-lethal, environmentally relevant concentrations (1, 10 and 100 μg/L). Sequential window acquisition of all theoretical mass spectra (SWATH-MS) identified over 1000 differentially abundant proteins within the 1% false discovery rate (FDR) threshold. The 24 h exposure resulted in a greater number of differentially abundant proteins, compared to 48 h exposure, for both contaminants. However, there were no statistically significant dose-response relationships for the number of differentially synthesised proteins, nor differences in the proportion of increased vs decreased proteins between or within exposure times. Known in vivo markers of contaminant exposure, superoxide dismutase and glutathione S-transferase, were differentially abundant following exposure to PCB153 and PFNA. SIGNIFICANCE: Cell-based (in vitro) proteomics provides an ethical and high throughput approach to understanding the impacts of chemical contamination on sea turtles. Through investigating effects of chemical dose and exposure duration on unique protein abundance in vitro, this study provides an optimised framework for conducting cell-based studies in wildlife proteomics, and highlights that proteins detected in vitro could act as biomarkers of chemical exposure and effect in vivo.

Metabolomics for in situ monitoring of attached Crassostrea gigas and Mytilus edulis: Effects of offshore wind farms on aquatic organisms

While offshore wind power has support from countries around the world, studies show that offshore wind farms (OWFs) may affect marine organisms. Environmental metabolomics is a high-throughput method that provides a snapshot of an organism's metabolic state. To elucidate the effects of OWFs on aquatic organisms, we studied, in situ, Crassostrea gigas and Mytilus edulis attached within and outside of OWFs and their reef areas. Our results show that epinephrine, sulphaniline, and inosine 5'-monophosphate were significantly increased and L-carnitine was significantly reduced in both Crassostrea and Mytilus species from the OWFs. This may be related to immune response, oxidative stress, energy metabolism and osmotic pressure regulation of aquatic organisms. Our study shows that active selection of biological monitoring methods for risk assessment is necessary and that metabolomics of attached shellfish is useful in elucidating the metabolic pathways of aquatic organisms in OWFs.

First Identification of a Large Set of Serine Hydrolases by Activity-Based Protein Profiling in Dibutyl Phthalate-Exposed Zebrafish Larvae

Despite the involvement of several serine hydrolases (SHs) in the metabolism of xenobiotics such as dibutyl phthalate (DBP), no study has focused on mapping this enzyme class in zebrafish, a model organism frequently used in ecotoxicology. Here, we survey and identify active SHs in zebrafish larvae and search for biological markers of SH type after exposure to DBP. Zebrafish were exposed to 0, 5, and 100 µg/L DBP from 4 to 120 h post-fertilization. A significant decrease in vitellogenin expression level of about 2-fold compared to the control was found in larvae exposed to 100 µg/L DBP for 120 h. The first comprehensive profiling of active SHs in zebrafish proteome was achieved with an activity-based protein profiling (ABPP) approach. Among 49 SHs identified with high confidence, one was the carboxypeptidase ctsa overexpressed in larvae exposed to 100 µg/L DBP for 120 h. To the best of our knowledge, this is the first time that a carboxypeptidase has been identified as deregulated following exposure to DBP. The overall results indicate that targeted proteomics approaches, such as ABPP, can, therefore, be an asset for understanding the mechanism of action related to xenobiotics in ecotoxicology.

Consideration of metabolomics and transcriptomics data in the context of using avian embryos for toxicity testing

Early-life stage (ELS) avian toxicity tests have been proposed as a more ethical alternative to traditional standardized tests with adult birds. At the same time, 'omics approaches are gaining traction in the field of avian toxicology, but little has been done to characterize the metabolome and transcriptome at different life stages. The present study uses 'omics data from toxicity tests of 8 environmental chemicals in ELS and adult Japanese quail (Coturnix japonica) to address this data gap. Previous analyses of these data focused on responses to each of the individual chemicals. Here, we consider data from all studies to describe variation in the metabolome and transcriptome between life stages and across independent experiments, irrespective of chemical treatment. Of the 230 metabolites detected in liver, 163 were shared between the two life stages. However, many of the targeted bile acids that were present in the adult liver were absent from ELS samples. For the transcriptome, >90% of the 18,364 detected transcripts were common to both life stages. Based on the 213 genes solely detected in ELS liver, the neuroactive ligand-receptor interaction pathway was significantly enriched. Multivariate and hierarchical clustering analyses revealed that variability among independent experiments was higher for the adult than the ELS studies at both the metabolomic and transcriptomic levels. Our results indicate concordance of the two approaches, with less variation between independent experiments in the ELS metabolome and transcriptome than in adults, lending support for the use of ELS as an alternative toxicity testing strategy.

The onset of active gill respiration in post-embryonic zebrafish (Danio rerio) larvae triggers an increased sensitivity to neurotoxic compounds

Originally designed as a general alternative to acute fish toxicity testing (AFT), the fish embryo toxicity test (FET) has become subject to concerns with respect to neurotoxic substances. Whereas oxygen uptake in the fish embryo primarily occurs via diffusion across the skin, juvenile and adult fish rely on active ventilation of the gills. As a consequence, substances including, e.g., neurotoxicants which prevent appropriate ventilation of gills ("respiratory failure syndrome") might lead to suffocation in juvenile and adult fish, but not in skin-breathing embryos. To investigate if this respiratory failure syndrome might play a role for the higher sensitivity of juvenile and adult fish to neurotoxicants, a modified acute toxicity test using post-embryonic, early gill-breathing life-stages of zebrafish was developed with chlorpyrifos, permethrin, lindane, aldicarb, ziram and aniline as test substances. Additionally, a comparative study into bioaccumulation of lipophilic substances with log_Kow > 3.5 and swimbladder deflation as potential side effects of the respiratory failure syndrome was performed with 4 d old skin-breathing and 12 d old gill-breathing zebrafish. With respect to acute toxicity, post-embryonic 12 d larvae proved to be more sensitive than both embryos (FET) and adult zebrafish (AFT) to all test substances except for permethrin. Accumulation of chlorpyrifos, lindane and permethrin was 1.3- to 5-fold higher in 4 d old than in 12 d old zebrafish, suggesting that (intermediate) storage of substances in the yolk might reduce bioavailability and prevent metabolization, which could be a further reason for lower toxicity in 4 d than in 12 d old zebrafish. Whereas ziram and aniline showed no significant effect on the swimbladder, zebrafish exposed to chlorpyrifos, lindane and permethrin showed significantly deflated swimbladders in 12 d old larvae; in the case of aldicarb, there was a significant hyperinflation in 4 d old larvae. Swimbladder deflation in post-embryonic 12 d zebrafish larvae might be hypothesized as a reason for a lack of internal oxygen supplies during the respiratory failure syndrome, whereas in 4 d old embryos cholinergic hyperinflation of the swimbladder dominates over other effects. Regarding acute lethality, the study provides further evidence that the switch from transcutaneous to branchial respiration in post-embryonic zebrafish life-stages might be the reason for the higher sensitivity of juvenile and adult fish to neurotoxic substances.

Fish metabolome from sub-urban lakes of the Paris area (France) and potential influence of noxious metabolites produced by cyanobacteria

The recent democratization of high-throughput molecular phenotyping allows the rapid expansion of promising untargeted multi-dimensional approaches (e.g. epigenomics, transcriptomics, proteomics, and/or metabolomics). Indeed, these emerging omics tools, processed for ecologically relevant species, may present innovative perspectives for environmental assessments, that could provide early warning of eco(toxico)logical impairments. In a previous pilot study (Sotton et al., Chemosphere 2019), we explore by ¹H NMR the bio-indicative potential of metabolomics analyses on the liver of 2 sentinel fish species (Perca fluviatilis and Lepomis gibbosus) collected in 8 water bodies of the peri-urban Paris' area (France). In the present study, we further investigate on the same samples the high potential of high-throughput UHPLC-HRMS/MS analyses. We show that the LC-MS metabolome investigation allows a clear separation of individuals according to the species, but also according to their respective sampling lakes. Interestingly, similar variations of Perca and Lepomis metabolomes occur locally indicating that site-specific environmental constraints drive the metabolome variations which seem to be influenced by the production of noxious molecules by cyanobacterial blooms in certain lakes. Thus, the development of such reliable environmental metabolomics approaches appears to constitute an innovative bio-indicative tool for the assessment of ecological stress, such as toxigenic cyanobacterial blooms, and aim at being further follow up.

Characterizing toxicity pathways of fluoxetine to predict adverse outcomes in adult fathead minnows (Pimephales promelas)

Current ecotoxicity testing programs are impeded as they predominantly rely on slow and expensive animal tests measuring adverse outcomes. Therefore, new approach methodologies (NAMs) increasingly involve short-term mechanistic assays that employ molecular endpoints to predict adverse outcomes of regulatory relevance. This study aimed to elucidate the application of NAMs in adult fathead minnows using fluoxetine (FLX) as a model compound. Fish were exposed to three FLX concentrations (measured: 2.42, 10.7, and 56.7 μgL^-1) and a control. After 96 h, molecular toxicity signatures were characterized using proteomics and transcriptomics analyses in livers and brains of a sub-set of fish. The remaining fish were sampled at 21 days and assessed for liver histopathology and morphometric measurements. Fecundity was monitored throughout the study. In the livers, 56.7 μgL^-1 FLX caused enrichment of PPAR signaling in the proteome and fatty acid-related pathways in the transcriptome, potential upstream responses that led to lipid-type vacuolation of hepatocytes, observed via histopathology. Upregulated genes in the brain suggested alterations in serotonin-related signaling processes and reproductive behaviour, which may explain the observed significant decrease in fecundity. While the relationships between molecular responses and adverse outcomes remain complex, this research provided important insights into the mechanistic toxicity of FLX.

Newly Hatched Stage I American Lobster (Homarus americanus) Survival Following Exposure to Physically and Chemically Dispersed Crude Oil

Standard model species are commonly used in toxicity tests due to their biological and technical advantages but studying native species increases the specificity and relevance of results generated for the potential risk assessment to an ecosystem. Accounting for intraspecies variability and other factors, such as chemical and physical characterization of test medium, is necessary to develop a reproducible bioassay for toxicity testing with native species. In this study, larval stage I American lobster (Homarus americanus), a commercially important and native species of Atlantic Canada, was used as the test species. Toxicity tests were first conducted by exposing lobster larvae to a reference toxicant of copper sulphate (CuSO₄) and then to physically and chemically (using Corexit 9500A) dispersed oil (WAF and CEWAF, respectively). The effect on larval survival was estimated by calculating the 24-h median effect concentration (24-h EC50), and there was no difference between WAF or CEWAF exposure when the results are reported on a total petroleum hydrocarbon (TPH) basis. The 24-h EC50s ranged from 2.54 to 9.73 mg TPH/L when all trials (n = 19) are considered together. The HC5 (hazardous concentration for 5 per cent of the population) value was 2.52 mg TPH/L and similar to the EC50 value when all trials were pooled. To evaluate the reproducibility of the lobster toxicity tests, inter-trial variability was determined, and the resultant coefficients of variation (%CV) were compared to those reported for two standard test species, mysid shrimp (Americamysis bahia) and inland silverside (Menidia beryillina). This comparison showed that the %CV for the lobster toxicity tests were lower than those for the standard species tests indicating that the described larval lobster toxicity bioassay produces reliable and repeatable results.

Multi-omic approach to evaluate the response of gilt-head sea bream (Sparus aurata) exposed to the UV filter sulisobenzone

Sulisobenzone (BP-4) is one of the benzophenone type UV filters most frequently detected in aquatic ecosystems. As a suspected endocrine disrupting compound, scarce information is available yet about other molecular effects and its mechanism of action. Here, we used an integrated transcriptomic and metabolomic approach to improve the current understanding on the toxicity of BP-4 towards aquatic species. Gilt-head sea bream individuals were exposed at environmentally relevant concentrations (10 μg L^-1) for 22 days. Transcriptomic analysis revealed 371 differentially expressed genes in liver while metabolomic analysis identified 123 differentially modulated features in plasma and 118 in liver. Integration of transcriptomic and metabolomic data showed disruption of the energy metabolism (>10 pathways related to the metabolism of amino acids and carbohydrates were impacted) and lipid metabolism (5 glycerophospholipids and the expression of 3 enzymes were affected), suggesting oxidative stress. We also observed, for the first time in vivo and at environmental relevant concentrations, the disruption of several enzymes involved in the steroid and thyroid hormones biosynthesis. DNA and RNA synthesis was also impacted by changes in the purine and pyrimidine metabolisms. Overall, the multiomic workflow presented here increases the evidence on suspected effects of BP-4 exposure and identifies additional modes of action of the compounds that could have been overlooked by using single omic approaches.

Sub-lethal exposure to malaria vector control pesticides causes alterations in liver metabolomics and behaviour of the African clawed frog (Xenopus laevis)

In this study we explore the sub-lethal effects of two malaria vector control pesticides, deltamethrin and dichlorodiphenyltrichloroethane (DDT), on Xenopus laevis by incorporating different levels of biological organisation. Pesticide accumulation in frog tissue was measured alongside liver metabolomics and individual swimming behaviour to assess whether changes presented at these different levels, and if such changes could be linked between levels. Results showed evidence of concentration dependent accumulation of DDT and its metabolites, but no measurable accumulation of deltamethrin in adult X. laevis after 96 h of exposure. Both DDT and deltamethrin were shown to cause alterations in the liver metabolome of X. laevis. We also showed that some of these changes can be enhanced in exposure to a mixture of these two pesticides. Initial behavioural responses recorded directly after exposure were seen in the form of decreased activity, less alterations between mobility states, and less time spent at the water surface. This response persisted after 96 h of exposure to a mixture of the two pesticides. This study shows that sub-lethal exposure to pesticides can alter the biochemical homeostasis of frogs with the potential to cascade onto behavioural and ecological levels in mixture exposure scenarios.

Trends in the Application of "Omics" to Ecotoxicology and Stress Ecology

Our ability to predict and assess how environmental changes such as pollution and climate change affect components of the Earth's biome is of paramount importance. This need positioned the fields of ecotoxicology and stress ecology at the center of environmental monitoring efforts. Advances in these interdisciplinary fields depend not only on conceptual leaps but also on technological advances and data integration. High-throughput "omics" technologies enabled the measurement of molecular changes at virtually all levels of an organism's biological organization and thus continue to influence how the impacts of stressors are understood. This bibliometric review describes literature trends (2000-2020) that indicate that more different stressors than species are studied each year but that only a few stressors have been studied in more than two phyla. At the same time, the molecular responses of a diverse set of non-model species have been investigated, but cross-species comparisons are still rare. While transcriptomics studies dominated until 2016, a shift towards proteomics and multiomics studies is apparent. There is now a wealth of data at functional omics levels from many phylogenetically diverse species. This review, therefore, addresses the question of how to integrate omics information across species.

Applying a gene-suite approach to examine the physiological status of wild-caught walleye (Sander vitreus)

Molecular techniques have been increasingly used in a conservation physiology framework to provide valuable information regarding the mechanisms underlying responses of wild organisms to environmental and anthropogenic stressors. In the present study, we developed a reference gill transcriptome for walleye (Sander vitreus), allowing us to pair a gene-suite approach (i.e. multiple genes across multiple cellular processes) with multivariate statistics to examine the physiological status of wild-caught walleye. For molecular analyses of wild fish, the gill is a useful target for conservation studies, not only because of its importance as an indicator of the physiological status of fish but also because it can be biopsied non-lethally. Walleye were non-lethally sampled following short- (~1.5 months) and long-term (~3.5 months) confinement in the Delta Marsh, which is located south of Lake Manitoba in Manitoba, Canada. Large-bodied walleye are confined in the Delta Marsh from late April to early August by exclusion screens used to protect the marsh from invasive common carp (Cyprinus carpio), exposing fish to potentially stressful water quality conditions. Principal components analysis revealed patterns of transcript abundance consistent with exposure of fish to increasingly high temperature and low oxygen conditions with longer holding in the marsh. For example, longer-term confinement in the marsh was associated with increases in the mRNA levels of heat shock proteins and a shift in the mRNA abundance of aerobic to anaerobic metabolic genes. Overall, the results of the present study suggest that walleye confined in the Delta Marsh may be exhibiting sub-lethal responses to high temperature and low oxygen conditions. These results provide valuable information for managers invested in mediating impacts to a local species of conservation concern. More broadly, we highlight the usefulness of pairing transcriptomic techniques with multivariate statistics to address potential confounding factors that can affect measured physiological responses of wild-caught fish.

Are we forgetting the "proteomics" in multi-omics ecotoxicology?

Proteomics plays a significant role in discerning the effects of chemical exposures in animal taxa. Multi-omics applications have become more pervasive in toxicology, however questions remain about whether proteomics is being utilized by the community to its full potential - are we placing too much stock in transcriptomics and other omics approaches for developing adverse outcome pathways? Proteins are more relevant than transcripts because they are direct mediators of the resulting phenotype. There is also rarely perfect stoichiometry between transcript and protein abundance and transcript abundance may not accurately predict physiologic response. Proteins direct all levels of phenotype: structural proteins dictate physical form, enzymes catalyze biochemical reactions, and proteins act as signaling proteins, antibodies, transporters, ion pumps, and transcription factors to control gene expression. Molecular initiating events (MIEs) of AOPs predominantly occur at the level of the protein (e.g. ligand-receptor binding) and proteomics can elucidate novel MIEs and mapping KEs in AOPs. This critical review highlights the need for proteomics in multi-omics studies in environmental toxicology and outlines steps required for inclusion and wider acceptance in chemical risk assessment. We also present case studies of multi-omics approaches that utilize proteomics and discuss some of the challenges and opportunities for proteomics in comparative ecotoxicology. Our intention is not to minimize the importance of other omics technologies, as each has strengths and limitations, but rather to encourage researchers to consider proteomics-based methods in multi-omics studies and AOP development.

The effects of tissue fixation on sequencing and transcript abundance of nucleic acids from microdissected liver samples of smallmouth bass (Micropterus dolomieu)

There is an increasing emphasis on effects-based monitoring to document responses associated with exposure to complex mixtures of chemicals, climate change, pathogens, parasites and other environmental stressors in fish populations. For decades aquatic monitoring programs have included the collection of tissues preserved for microscopic pathology. Consequently, formalin-fixed, paraffin-embedded (FFPE) tissue can be an important reservoir of nucleic acids as technologies emerge that utilize molecular endpoints. Despite the cross-linking effects of formalin, its impact on nucleic acid quality and concentration, amplification, and sequencing are not well described. While fresh-frozen tissue is optimal for working with nucleic acids, FFPE samples have been shown to be conducive for molecular studies. Laser capture microdissection (LCM) is one technology which allows for collection of specific regions or cell populations from fresh or preserved specimens with pathological alterations, pathogens, or parasites. In this study, smallmouth bass (Micropterus dolomieu) liver was preserved in three different fixatives, including 10% neutral buffered formalin (NBF), Z-Fix® (ZF), and PAXgene® (PG) for four time periods (24 hr, 48 hr, seven days, and 14 days). Controls consisted of pieces of liver preserved in RNALater® or 95% ethanol. Smallmouth bass were chosen as they are an economically important sportfish and have been utilized as indicators of exposure to endocrine disruptors and other environmental stressors. Small liver sections were cut out with laser microdissection and DNA and RNA were purified and analyzed for nucleic acid concentration and quality. Sanger sequencing and the NanoString nCounter® technology were used to assess the suitability of these samples in downstream molecular techniques. The results revealed that of the formalin fixatives, NBF samples fixed for 24 and 48 hr were superior to ZF samples for both Sanger sequencing and the Nanostring nCounter®. The non-formalin PAXgene® samples were equally successful and they showed greater stability in nucleic acid quality and concentration over longer fixation times. This study demonstrated that small quantities of preserved tissue from smallmouth bass can be utilized in downstream molecular techniques; however, future studies will need to optimize the methods presented here for different tissue types, fish species, and pathological conditions.

Multifactorial Analysis of Environmental Metabolomic Data in Ecotoxicology: Wild Marine Mussel Exposed to WWTP Effluent as a Case Study

Environmental metabolomics is a powerful approach to investigate the response of organisms to contaminant exposure at a molecular scale. However, metabolomic responses to realistic environmental conditions can be hindered by factors intrinsic to the environment and the organism. Hence, a well-designed experimental exposure associated with adequate statistical analysis could be helpful to better characterize and relate the observed variability to its different origins. In the current study, we applied a multifactorial experiment combined to Analysis of variance Multiblock Orthogonal Partial Least Squares (AMOPLS), to assess the metabolic response of wild marine mussels, Mytilus galloprovincialis, exposed to a wastewater treatment plant effluent, considering gender as an experimental factor. First, the total observed variability was decomposed to highlight the contribution of each effect related to the experimental factors. Both the exposure and the interaction gender × exposure had a statistically significant impact on the observed metabolic alteration. Then, these metabolic patterns were further characterized by analyzing the individual variable contributions to each effect. A main change in glycerophospholipid levels was highlighted in both males and females as a common response, possibly caused by oxidative stress, which could lead to reproductive disorders, whereas metabolic alterations in some polar lipids and kynurenine pathway were rather gender-specific. This may indicate a disturbance in the energy metabolism and immune system only in males. Finally, AMOPLS is a useful tool facilitating the interpretation of complex metabolomic data and is expected to have a broad application in the field of ecotoxicology.

In-plate toxicometabolomics of single zebrafish embryos

Toxicometabolomic studies involving zebrafish embryos have become increasingly popular for linking apical endpoints to biochemical perturbations as part of adverse outcome pathway determination. These experiments involve pooling embryos to generate sufficient biomass for metabolomic measurement, which adds both time and cost. To address this limitation, we developed a high-throughput toxicometabolomic assay involving single zebrafish embryos. Incubation, microscopy, embryo extraction, and instrumental metabolomic analysis were all performed in the same 96-well plate, following acquisition of conventional toxicological endpoints. The total time for the assay (including testing of 6 doses/n = 12 embryos per dose plus positive and negative controls, assessing conventional endpoints, instrumental analysis, data processing and multivariate statistics) is <14 days. Metabolomic perturbations at low dose were linked statistically to those observed at high dose and in the presence of an adverse effect, thereby contextualizing omic data amongst apical endpoints. Overall, this assay enables collection of high resolution metabolomic data in a high throughput manner, suitable for mode of action hypothesis generation in the context of pharmaceutical or toxicological screening.

Complete transcriptome assembly and annotation of a critically important amphipod species in freshwater ecotoxicological risk assessment: Gammarus fossarum

Because of their crucial role in ecotoxicological risk assessment, amphipods (Crustacea) are commonly employed as model species in a wide range of studies. However, despite their ecological importance, their genome has not yet been completely annotated and molecular mechanisms underlying key pathways, such as the serotonin pathway, in development of ecotoxicological biomarkers of exposure to neuroactive pharmaceuticals are still poorly understood. Furthermore, genetic similarities and discrepancies with other model arthropods (e.g., Drosophila melanogaster) have not been completely clarified. In this report, we present a new transcriptome assembly of Gammarus fossarum, an important amphipod species, widespread in Central Europe. RNA-Seq with Illumina HiSeq technology was used to analyse samples extracted from total internal tissues. We used the Trinity and Trinotate software suites for transcriptome assembly and annotation, respectively. The quality of this assembly and the affiliated targeted homology searches greatly enrich the molecular knowledge on this species. Because of the lack of publicly available molecular information on the serotonin pathway, we also highlighted sequence homologies and divergences of the genes encoding the serotonin pathway components of the well-annotated arthropod D. melanogaster, and Crustacea with the corresponding genes of our assembly. An inferior number of hits was found when running a BLAST analysis of both D. melanogaster and Crustacea mRNA sequences encoding serotonin receptors available in GenBank against the total assembly, compared to other serotonin pathway components. A lack of information on important components for serotonin biosynthesis and vesicle endocytosis (i.e., tryptophan hydroxylase and vesicular monoamine transporter) in Crustacea was also brought to light. Our results will provide an extensive transcriptional resource for this important species in ecotoxicological risk assessment and highlight the need for a more detailed categorization of neuronal pathways components in invertebrates.

Metabolome response to anthropogenic contamination on microalgae: a review

BACKGROUND

Microalgae play a key role in ecosystems and are widely used in ecological status assessment. Research focusing on such organisms is then well developed and essential. Anyway, approaches for a better comprehension of their metabolome's response towards anthropogenic stressors are only emerging.

AIM OF REVIEW

This review presents the biochemical responses of various microalgae species towards several contaminants including metals and chemicals as pesticides or industrial compounds. We aim to provide a comprehensive and up-to-date overview of analytical approaches deciphering anthropogenic contaminants impact on microalgae metabolome dynamics, in order to bring out relevant biochemical markers that could be used for risk assessment.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Studies to date on ecotoxicological metabolomics on microalgae are highly heterogeneous in both analytical techniques and resulting metabolite identification. There is a real need for studies using complementary approaches to determine biomarkers usable for ecological risk assessment.

Comparative Proteomics on Deep-Sea Amphipods after in Situ Copper Exposure

The interest in deep-sea mining increased along with the environmental concerns of these activities to the deep-sea fauna. The discovery of optimal biomarkers of deep-sea mining activities in deep-sea species is a crucial step toward the supply of important ecological information for environmental impact assessment. In this study, an in situ copper exposure experiment was performed on deep-sea scavenging amphipods. Abyssorchomene distinctus individuals were selected among all the exposed amphipods for molecular characterization. Copper concentration within the gut was assessed, followed by a tandem mass tag-based coupled with two-dimensional liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) applied to identify and quantify the protein expression changes after 48 h of exposure. 2937 proteins were identified and annotated, and 1918 proteins among all identified proteins were assigned by at least two nonambiguous peptides. The screening process was performed based on the differences in protein abundance and the specific correlation between the proteins and copper in previous studies. These differentially produced proteins include Na⁺/K⁺ ATPase, cuticle, chitinase, and proteins with unknown function. Their abundances showed correlation with copper and had high sensitivity to indicate the copper level, being here proposed as biomarker candidates for deep-sea mining activities in the future. This is a key step in the development of environmental impact assessment of deep-sea mining activities integrating ecotoxicological data.

Mitochondrial dysfunction, apoptosis and transcriptomic alterations induced by four strobilurins in zebrafish (Danio rerio) early life stages

Though the toxicity of strobilurins on non-target aquatic organisms has been characterized, the associated toxic mechanisms have not been fully explored. The present study showed that the larval stage was the most sensitive developmental stage in zebrafish, and pyraclostrobin (PY) had the highest acute toxicity to embryos, larvae, juvenile and adult with 96 h-LC₅₀ at 0.048 mg/L, 0.029 mg/L, 0.039 mg/L, 0.031 mg/L respectively, when compared with the toxicity of trifloxystrobin (TR), kresoxim-methyl (KM) and azoxystrobin (AZ) at corresponding developmental stage. Then we investigated the transcriptomics and developmental toxicity of TR, KM, AZ and PY on zebrafish embryos after 72 h exposure. RNA-seq revealed that the pathways related to cell apoptosis and cancer, and cellular components organelle membrane and mitochondrion, were markedly affected after TR, KM, AZ and PY exposure during zebrafish early life stages. The results were further confirmed by the induction of antioxidant enzymes catalase (CAT) and superoxide dismutase (SOD) activities, the elevation of H₂O₂, malondialdehyde (MDA) and reactive oxygen species (ROS) level, as well as the reduction of intracellular calcium ions (Ca²⁺) and mitochondrial membrane potential (MMP), which indicated that strobilurins could cause mitochondrial dysfunction and cell apoptosis. The present study was performed a systematic analysis of strobilurins to zebrafish at multi-levels, which provided suggestions for further investigation of molecular mechanisms underlying the toxicity induced by strobilurins on aquatic organisms.

Comparative proteomics in the wild: Accounting for intrapopulation variability improves describing proteome response in a Gammarus pulex field population exposed to cadmium

High-throughput proteomics can be performed on animal sentinels for discovering key molecular biomarkers signing the physiological response and adaptation of organisms. Ecotoxicoproteomics is today amenable by means of proteogenomics to small arthropods such as Gammarids which are well known sentinels of aquatic environments. Here, we analysed two regional Gammarus pulex populations to characterize the potential proteome divergence induced in one site by natural bioavailable mono-metallic contamination (cadmium) compared to a non-contaminated site. Two RNAseq-derived protein sequence databases were established previously on male and female individuals sampled from the reference site. Here, individual proteomes were acquired on 10 male and 10 female paired organisms sampled from each site. Proteins involved in protein lipidation, carbohydrate metabolism, proteolysis, innate immunity, oxidative stress response and lipid transport were found more abundant in animals exposed to cadmium, while hemocyanins were found in lower abundance. The intrapopulation proteome variability of long-term exposed G. pulex was inflated relatively to the non-contaminated population. These results show that, while remaining a challenge for such organisms with not yet sequenced genomes, taking into account intrapopulation variability is important to better define the molecular players induced by toxic stress in a comparative field proteomics approach.

Amitriptyline at an Environmentally Relevant Concentration Alters the Profile of Metabolites Beyond Monoamines in Gilt-Head Bream

The antidepressant amitriptyline is a widely used selective serotonin reuptake inhibitor that is found in the aquatic environment. The present study investigates alterations in the brain and the liver metabolome of gilt-head bream (Sparus aurata) after exposure at an environmentally relevant concentration (0.2 µg/L) of amitriptyline for 7 d. Analysis of variance-simultaneous component analysis is used to identify metabolites that distinguish exposed from control animals. Overall, alterations in lipid metabolism suggest the occurrence of oxidative stress in both the brain and the liver-a common adverse effect of xenobiotics. However, alterations in the amino acid arginine are also observed. These are likely related to the nitric oxide system that is known to be associated with the mechanism of action of antidepressants. In addition, changes in asparagine and methionine levels in the brain and pantothenate, uric acid, and formylisoglutamine/N-formimino-L-glutamate levels in the liver could indicate variation of amino acid metabolism in both tissues; and the perturbation of glutamate in the liver implies that the energy metabolism is also affected. These results reveal that environmentally relevant concentrations of amitriptyline perturb a fraction of the metabolome that is not typically associated with antidepressant exposure in fish. Environ Toxicol Chem 2019;00:1-13. © 2019 SETAC.

Blood Transcriptomics Analysis of Fish Exposed to Perfluoro Alkyls Substances: Assessment of a Non-Lethal Sampling Technique for Advancing Aquatic Toxicology Research

In contrast to mammals, the blood from other vertebrates such as fish contains nucleated red cells. Using a fathead minnow ( Pimephales promelas) oligonucleotide microarray, we compared altered transcripts in the liver and whole blood after exposure to environmentally relevant concentrations of perfluorooctanesulfonic acid (PFOS) and a mixture of seven types of perfluoro alkyl substances (PFAS), including perfluorooctanoic acid (PFOA). We used quantitative polymerase chain reactions and cell-based assays to confirm the main effects and found that blood responded with a greater number of altered genes than the liver. The exposure to PFAS altered similar genes with central roles in a cellular pathway in both tissues, including estrogen receptor α and peroxisome proliferator activator β and γ, indicating that the genes previously associated with PFAS exposure are differentially expressed in blood and liver. The altered transcripts are involved with cholesterol metabolism and mitochondrial function. Our data confirmed that PFAS are weak xenoestrogens and exert effects on DNA integrity. Gene expression profiling from blood samples not related with the immune system, including very-low-density lipoprotein, vitellogenin, estrogen receptor, and thyroid hormone receptor, demonstrated that blood is a useful tissue for assessing endocrine disruption in non-mammalian vertebrates. We conclude that the use of blood for non-lethal sampling in genomics studies is informative and particularly useful for assessing the effects of pollution in endangered species. Further, using blood will reduce animal use and widen the experimental design options for studying the effects of contaminant exposure on wildlife.

Untargeted NMR-based metabolomics for field-scale monitoring: Temporal reproducibility and biomarker discovery in mosquitofish (Gambusia holbrooki) from a metal(loid)-contaminated wetland

There is considerable interest in applying omics techniques, which have proven extremely valuable for laboratory-based toxicology studies, towards field-scale ecotoxicology and environmental monitoring. Concerns that confounding factors in natural ecosystems may exacerbate variability in omics datasets must be addressed to validate the transition from laboratory to field. This study explores how temporal variability related to seasonal and climatic trends influence qualitative and quantitative metabolomics outcomes, in fish from reference and metal(loid)-polluted wetlands in Australia. Female mosquitofish (Gambusia holbrooki) were sampled on two separate occasions, from a rehabilitated tailings wetland at the site of historic antimony (Sb) processing and a reference wetland with comparable water quality. The first sampling coincided with greater monthly rainfall and colder water temperature, whereas the second sampling was drier and water was warmer. Despite temporal changes and associated differences in metal(loid) concentrations, site differences in metabolite profiles were qualitatively very similar between sampling events. However, quantitative differences were observed, with a greater number of significantly altered metabolites identified during the second sampling event, which coincided with greater metal(loid) concentrations in both water and fish. The majority of identified metabolites were elevated in fish from the contaminated wetland, but with notable decreases in several metabolites that are known to play a role in various aspects of metal(loid) binding, detoxification and excretion. Specifically, decreased aspartate, histidine, myo-inositol, taurine and choline were observed in fish from the contaminated wetland, and may therefore represent a metabolite suite that is broadly indicative of metal toxicity. Quantitative differences between sampling events are suggestive of a dose-response relationship observable at the cellular level which, if harnessed, may be useful for assigning levels of concern based on the degree of change in a multi-parameter set of metabolite biomarkers.

Development, characterization and comparisons of targeted and non-targeted metabolomics methods

The potential of a metabolomics method to detect statistically significant perturbations in the metabolome of an organism is enhanced by excellent analytical precision, unequivocal identification, and broad metabolomic coverage. While the former two metrics are usually associated with targeted metabolomics and the latter with non-targeted metabolomics, a systematic comparison of the performance of both approaches has not yet been carried out. The present work reports on the development and performance evaluation of separate targeted and non-targeted metabolomics methods. The targeted approach facilitated determination of 181 metabolites (quantitative analysis of 18 amino acids, 11 biogenic amines, 5 neurotransmitters, 5 nucleobases and semi-quantitative analysis of 50 carnitines, 83 phosphatidylcholines, and 9 sphingomyelins) using ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) and flow injection-tandem mass spectrometry (FI-MS/MS). Method accuracy and/or precision were assessed using replicate samples of NIST SRM1950 as well as fish liver and brain tissue from Gilthead Bream (Sparus aurata). The non-target approach involved UPLC-high resolution (Orbitrap) mass spectrometry (UPLC-HRMS). Testing of ionization mode and stationary phase revealed that a combination of positive electrospray ionization and HILIC chromatography produced the largest number of chromatographic features during non-target analysis. Furthermore, an evaluation of 4 different sequence drift correction algorithms, and combinations thereof, revealed that batchCorr produced the best precision in almost every test. However, even following correction of non-target data for signal drift, the precision of targeted data was better, confirming our existing assumptions about the strengths of targeted metabolomics. Finally, the accuracy of the online MS²-library mzCloud was evaluated using reference standards for 38 different metabolites. This is among the few studies that have systematically evaluated the performance of targeted and non-targeted metabolomics and provides new insight into the advantages and disadvantages of each approach.

Metabolite profiles of striped marsh frog (Limnodynastes peronii) larvae exposed to the anti-androgenic fungicides vinclozolin and propiconazole are consistent with altered steroidogenesis and oxidative stress

Amphibians use wetlands in urban and agricultural landscapes for breeding, growth and development. Fungicides and other pesticides used in these areas have therefore been identified as potential threats that could contribute towards amphibian population declines. However, relatively little is known about how such chemicals influence sensitive early life-stages or how short episodic exposures influence sub-lethal physiological and metabolic pathways. The present study applied untargeted metabolomics to evaluate effects in early post-hatch amphibian larvae exposed to the anti-androgenic fungicides vinclozolin and propiconazole. Recently hatched (Gosner developmental stage 25) striped marsh frog (Limnodynastes peronii) larvae were exposed for 96 h to vinclozolin at 17.5, 174.8 and 1748.6 nM and propiconazole at 5.8, 58.4 and 584.4 nM. Nuclear Magnetic Resonance (NMR) spectroscopy was performed on polar metabolites obtained from whole-body extracts. Both fungicides altered metabolite profiles compared to control animals at all concentrations tested, and there were notable differences between the two chemicals. Overall responses were consistent with altered steroidogenesis and/or cholesterol metabolism, with inconsistent responses between the two fungicides likely reflecting minor differences in the mechanisms of action of these chemicals. Broad down-regulation of the tricarboxylic acid (TCA) cycle was also observed and is indicative of oxidative stress. Interestingly, formic acid was significantly increased in larvae exposed to vinclozolin but not propiconazole, suggesting this metabolite may serve as a useful biomarker of exposure to androgen-receptor binding anti-androgenic contaminants. This study demonstrates the power of untargeted metabolomics for distinguishing between similarly acting, but distinct, pollutants and for unraveling non-endocrine responses resulting from exposure to known endocrine active contaminants.

1H NMR-based metabolomics reveals interactive effects between the carrier solvent methanol and a pharmaceutical mixture in an amphibian developmental bioassay with Limnodynastes peronii

Organic carrier solvents are used in aquatic toxicity testing to improve chemical solubility and facilitate the exploration of dose-response relationships. Both water- and solvent-control groups are normally included in these scenarios to ensure that the solvent itself has no effect on the test organism, but this fails to consider possible interactive effects between carrier solvents and contaminants of interest. We explored this topic by exposing Limnodynastes peronii tadpoles to a mixture of common water-soluble pharmaceuticals (diclofenac, metformin and valproic acid) in the presence and absence of the carrier solvent methanol, according to standard developmental bioassay methodology. Nuclear Magnetic Resonance (NMR) spectroscopy was applied as a platform for untargeted metabolomics, to compare broad sub-lethal hepatotoxicity in solvent- and solvent-free exposure scenarios. Considerable interactive effects were identified between the pharmaceutical mixture and a typical dose of methanol (0.003%). Specifically, pronounced differences were observed between the solvent- and solvent-free exposure groups for leucine, acetate, glutamine, citrate, glycogen, tyrosine, arginine, purine nucleotides and an unidentified metabolite at 6.53 ppm. Various other metabolites exhibited similar disparity related to the use of carrier solvent, but the interactions were non-significant. These results raise important questions about the use of carrier solvents for chemical exposures in aquatic ecotoxicology, and particularly for studies interested in sub-lethal mechanistic information and/or biomarker discovery.

Are we closer to the vision? A proposed framework for incorporating omics into environmental assessments

Environmental science has benefited a great deal from omics-based technologies. High-throughput toxicology has defined adverse outcome pathways (AOPs), prioritized chemicals of concern, and identified novel actions of environmental chemicals. While many of these approaches are conducted under rigorous laboratory conditions, a significant challenge has been the interpretation of omics data in "real-world" exposure scenarios. Clarity in the interpretation of these data limits their use in environmental monitoring programs. In recent years, one overarching objective of many has been to address fundamental questions concerning experimental design and the robustness of data collected under the broad umbrella of environmental genomics. These questions include: (1) the likelihood that molecular profiles return to a predefined baseline level following remediation efforts, (2) how reference site selection in an urban environment influences interpretation of omics data and (3) what is the most appropriate species to monitor in the environment from an omics point of view. In addition, inter-genomics studies have been conducted to assess transcriptome reproducibility in toxicology studies. One lesson learned from inter-genomics studies is that there are core molecular networks that can be identified by multiple laboratories using the same platform. This supports the idea that "omics-networks" defined a priori may be a viable approach moving forward for evaluating environmental impacts over time. Both spatial and temporal variability in ecosystem structure is expected to influence molecular responses to environmental stressors, and it is important to recognize how these variables, as well as individual factor (i.e. sex, age, maturation), may confound interpretation of network responses to chemicals. This mini-review synthesizes the progress made towards adopting these tools into environmental monitoring and identifies future challenges to be addressed, as we move into the next era of high throughput sequencing. A conceptual framework for validating and incorporating molecular networks into environmental monitoring programs is proposed. As AOPs become more defined and their potential in environmental monitoring assessments becomes more recognized, the AOP framework may prove to be the conduit between omics and penultimate ecological responses for environmental risk assessments.

The utility of transcriptomics in fish conservation

There is growing recognition of the need to understand the mechanisms underlying organismal resilience (i.e. tolerance, acclimatization) to environmental change to support the conservation management of sensitive and economically important species. Here, we discuss how functional genomics can be used in conservation biology to provide a cellular-level understanding of organismal responses to environmental conditions. In particular, the integration of transcriptomics with physiological and ecological research is increasingly playing an important role in identifying functional physiological thresholds predictive of compensatory responses and detrimental outcomes, transforming the way we can study issues in conservation biology. Notably, with technological advances in RNA sequencing, transcriptome-wide approaches can now be applied to species where no prior genomic sequence information is available to develop species-specific tools and investigate sublethal impacts that can contribute to population declines over generations and undermine prospects for long-term conservation success. Here, we examine the use of transcriptomics as a means of determining organismal responses to environmental stressors and use key study examples of conservation concern in fishes to highlight the added value of transcriptome-wide data to the identification of functional response pathways. Finally, we discuss the gaps between the core science and policy frameworks and how thresholds identified through transcriptomic evaluations provide evidence that can be more readily used by resource managers.

Exposure to wastewater effluent affects fish behaviour and tissue-specific uptake of pharmaceuticals

Pharmaceutical active compounds (PhACs) are increasingly being reported in wastewater effluents and surface waters around the world. The presence of these products, designed to modulate human physiology and behaviour, has created concern over whether PhACs similarly affect aquatic organisms. Though laboratory studies are beginning to address the effects of individual PhACs on fish behaviour, few studies have assessed the effects of exposure to complex, realistic wastewater effluents on fish behaviour. In this study, we exposed a wild, invasive fish species-the round goby (Neogobius melanostomus)-to treated wastewater effluent (0%, 50% or 100% effluent dilutions) for 28days. We then determined the impact of exposure on fish aggression, an important behaviour for territory acquisition and defense. We found that exposure to 100% wastewater effluent reduced the number of aggressive acts that round goby performed. We complimented our behavioural assay with measures of pharmaceutical uptake in fish tissues. We detected 11 of 93 pharmaceutical compounds that we tested for in round goby tissues, and we found that concentration was greatest in the brain followed by plasma, then gonads, then liver, and muscle. Fish exposed to 50% and 100% effluent had higher tissue concentrations of pharmaceuticals and concentrated a greater number of pharmaceutical compounds compare to control fish exposed to no (0%) effluent. Exposed fish also showed increased ethoxyresorufin-O-deethylase (EROD) activity in liver tissue, suggesting that fish were exposed to planar halogenated/polycyclic aromatic hydrocarbons (PHHs/PAHs) in the wastewater effluent. Our findings suggest that fish in effluent-dominated systems may have altered behaviours and greater tissue concentration of PhACs. Moreover, our results underscore the importance of characterizing exposure to multiple pollutants, and support using behaviour as a sensitive tool for assessing animal responses to complex contaminant mixtures, like wastewater effluent.

Altered expression of metabolites and proteins in wild and caged fish exposed to wastewater effluents in situ

Population growth has led to increased global discharges of wastewater. Contaminants that are not fully removed during wastewater treatment, such as pharmaceuticals and personal care products (PPCPs), may negatively affect aquatic ecosystems. PPCPs can bioaccumulate causing adverse health effects and behavioural changes in exposed fish. To assess the impact of PPCPs on wild fish, and to assess whether caged fish could be used as a surrogate for resident wild fish in future monitoring, we caged goldfish in a marsh affected by discharges of wastewater effluents (Cootes Paradise, Lake Ontario, Canada). We collected plasma from resident wild goldfish, and from goldfish that we caged in the marsh for three weeks. We analyzed the plasma proteome and metabolome of both wild and caged fish. We also compared proteomic and metabolic responses in caged and wild fish from the marsh to fish caged at a reference site (Jordan Harbour Conservation Area). We identified significant changes in expression of over 250 molecules that were related to liver necrosis, accumulation and synthesis of lipids, synthesis of cyclic AMP, and the quantity of intracellular calcium in fish from the wastewater affected marsh. Our results suggest that PPCPs could be affecting the health of wild fish populations.

Estrogen-responsive gene networks in the teleost liver: What are the key molecular indicators?

An overarching goal of environmental genomics is to leverage sensitive suites of markers that are robust and reliable to assess biological responses in a range of species inhabiting variable environments. The objective of this study was to identify core groups of transcripts and molecular signaling pathways that respond to 17alpha-ethylinestadiol (EE2), a ubiquitous estrogenic contaminant, using transcriptome datasets generated from six independent laboratories. We sought to determine which biomarkers and gene networks were those most robust and reliably detected in multiple laboratories. Six laboratories conducted microarray analysis in pieces of the same liver from male fathead minnows exposed to ∼15ng/L EE2 for 96h. There were common transcriptional networks identified in every dataset. These included down-regulation of gene networks associated with blood clotting, complement activation, triglyceride storage, and xenobiotic metabolism. Noteworthy was that more than ∼85% of the gene networks were suppressed by EE2. Leveraging both these data and those mined from the Comparative Toxicogenomics Database (CTD), we narrowed in on an EE2-responsive transcriptional network. All transcripts in this network responded ∼±5-fold or more to EE2, increasing reliability of detection. This network included estrogen receptor alpha, transferrin, myeloid cell leukemia 1, insulin like growth factor 1, insulin like growth factor binding protein 2, and methionine adenosyltransferase 2A. This estrogen-responsive interactome has the advantage over single markers (e.g. vitellogenin) in that these entities are directly connected to each other based upon evidence of expression regulation and protein binding. Thus, it represents an interacting functional suite of estrogenic markers. Vitellogenin, the gold standard for estrogenic exposures, can show high individual variability in its response to estrogens, and the use of a multi-gene approach for estrogenic chemicals is expected to improve sensitivity. In our case, the coefficient of variation was significantly lowered by the gene network (∼67%) compared to Vtg alone, supporting the use of this transcriptional network as a sensitive alternative for detecting estrogenic effluents and chemicals. We propose that screening chemicals for estrogenicity using interacting genes within a defined expression network will improve sensitivity, accuracy, and reduce the number of animals required for endocrine disruption assessments.

Ecotoxico-Proteomics for Aquatic Environmental Monitoring: First in Situ Application of a New Proteomics-Based Multibiomarker Assay Using Caged Amphipods

As a proof of principle, a selected reaction monitoring (SRM) mass spectrometry-based methodology was applied to the simultaneous quantification of dozens of protein biomarkers in caged amphipods (Gammarus fossarum). We evaluated the suitability of the methodology to assess complex field contaminations through its application in the framework of a regional river monitoring network. Thanks to the high throughput acquisition of biomarker levels in G. fossarum exposed in four reference and 13 contaminated sites, we analyzed the individual responses of 38 peptides reporting for 25 proteins of interest in 170 organisms. Responses obtained in contaminated sites included inductions of vitellogenin-like proteins in male organisms, inductions of Na⁺K⁺/ATPases, and strong inhibitions of molt-related proteins such as chitinase and JHE-carboxylesterase. Proteins from detoxification and immunity processes were also found modulated in abundance. Summarizing, the results presented here show that the SRM strategy developed for multibiomarker measurement paves a very promising way to define multiple indicators of the health status of sentinel organisms for environmental hazard assessment.

Reproductive Impacts of Endocrine-Disrupting Chemicals on Wildlife Species: Implications for Conservation of Endangered Species

Wildlife have proven valuable to our understanding of the potential effects of endocrine-disrupting chemicals (EDCs) on human health by contributing considerably to our understanding of the mechanisms and consequences of EDC exposure. But the threats EDCs present to populations of wildlife species themselves are significant, particularly for endangered species whose existence is vulnerable to any reproductive perturbation. However, few studies address the threats EDCs pose to endangered species owing to challenges associated with their study. Here, we highlight those barriers and review the available literature concerning EDC effects on endangered species. Drawing from other investigations into nonthreatened wildlife species, we highlight opportunities for new approaches to advance our understanding and potentially mitigate the effects of EDCs on endangered species to enhance their fertility.

Nonlethal laparoscopic detection of intersex (testicular oocytes) in largemouth bass (Micropterus salmoides) and smallmouth bass (Micropterus dolomieu)

Intersex in wild fish populations has received considerable attention in the scientific literature and public media. Conventional detection of testicular oocytes, the presence of immature oocytes within testis of male fish, commonly employs transverse sectioning of excised testis and is lethal. The present study used a nonlethal laparoscopic technique to collect biopsies of testis from black bass, entering the body cavity via the genital pore. Detection of testicular oocytes was compared between biopsy and conventional methods using 79 smallmouth bass (Micropterus dolomieu) from 8 sites and 68 largemouth bass (Micropterus salmoides) from 4 sites. The 2 methods performed similarly at sites where testicular oocyte severity was moderate or high (6 of 8 smallmouth bass sites), whereas transverse sectioning resulted in superior testicular oocyte detection at sites where severity was low (2 of 8 smallmouth bass sites and all 4 largemouth bass sites). In smallmouth bass, testicular oocyte prevalence by transverse and biopsy methods was strongly correlated across sites (r²  = 0.81) and severity reported by enumeration of testicular oocytes was moderately correlated across sites (r² = 0.59). Survival of a subset of largemouth bass (n = 20) to 28 d after laparoscopic surgery was 90%. The present research indicates that laparoscopy may be useful for monitoring the prevalence and severity of testicular oocytes in Micropterus species, particularly when lethal sampling is precluded. Environ Toxicol Chem 2017;36:1924-1933. © 2016 SETAC.

A multi-omic approach to elucidate low-dose effects of xenobiotics in zebrafish (Danio rerio) larvae

Regulatory-approved toxicity assays such as the OECD Fish Embryo Toxicity Assay (TG236) allow correlation of chemical exposure to adverse morphological phenotypes. However, these assays are ineffective in assessing sub-lethal (i.e. low-dose) effects, or differentiating between similar phenotypes induced by different chemicals. Inclusion of multi-omic analyses in studies investigating xenobiotic action provides improved characterization of biological response, thereby enhancing prediction of toxicological outcomes in whole animals in the absence of morphological effects. In the current study, we assessed perturbations in both the metabolome and transcriptome of zebrafish (Danio rerio; ZF) larvae exposed from 96 to 120h post fertilization to environmental concentrations of acetaminophen (APAP), diphenhydramine (DH), carbamazepine (CBZ), and fluoxetine (FLX); common pharmaceuticals with known mechanisms of action. Multi-omic responses were evaluated independently and integrated to identify molecular interactions and biological relevance of the responses. Results indicated chemical- and dose-specific changes suggesting differences in the time scale of transcript abundance and metabolite production. Increased impact on the metabolome relative to the transcriptome in FLX-treated animals suggests a stronger post-translational effect of the treatment. In contrast, the transcriptome showed higher sensitivity to perturbation in DH-exposed animals. Integration of 'omic' responses using multivariate approaches provided additional insights not obtained by independent 'omic' analyses and demonstrated that the most distinct overall response profiles were induced following low-dose exposure for all 4 pharmaceuticals. Importantly, changes in transcript abundance corroborated with predictions from metabolomic enrichment analyses and the identified perturbed biological pathways aligned with known xenobiotic mechanisms of action. This work demonstrates that a multi-omic toxicological approach, coupled with a sensitive animal model such as ZF larvae, can help characterize the toxicological relevance of acute low-dose chemical exposures.

Spotlight on environmental omics and toxicology: a long way in a short time

The applications for high throughput omics technologies in environmental science have increased dramatically in recent years. Transcriptomics, proteomics, and metabolomics have been used to study how chemicals in our environment affect both aquatic and terrestrial organisms, and the characterization of molecular initiating events is a significant goal in toxicology to better predict adverse responses to toxicants. This special journal edition demonstrates the scope of the science that leverages omics-based methods in both laboratory and wild populations within the context of environmental toxicology, ranging from fish to mammals. It is important to recognize that the environment comprises one axis of the One Health concept - the idea that human health is unequivocally intertwined to our environment and to the organisms that inhabit that environment. We have much to learn from a comparative approach, and studies that integrate the transcriptome, proteome, and the metabolome are expected to offer the most detailed mechanism-based adverse outcome pathways that are applicable for use in both environmental monitoring and risk assessment.

Xenobiotics Produce Distinct Metabolomic Responses in Zebrafish Larvae (Danio rerio)

Sensitive and quantitative protocols for characterizing low-dose effects are needed to meet the demands of 21st century chemical hazard assessment. To test the hypothesis that xenobiotic exposure at environmentally relevant concentrations produces specific biochemical fingerprints in organisms, metabolomic perturbations in zebrafish (Danio rerio) embryo/larvae were measured following 24 h exposures to 13 individual chemicals covering a wide range of contaminant classes. Measured metabolites (208 in total) included amino acids, biogenic amines, fatty acids, bile acids, sugars, and lipids. The 96-120 h post-fertilization developmental stage was the most appropriate model for detecting xenobiotic-induced metabolomic perturbations. Metabolomic fingerprints were largely chemical- and dose-specific and were reproducible in multiple exposures over a 16-month period. Furthermore, chemical-specific responses were detected in the presence of an effluent matrix; importantly, in the absence of morphological response. In addition to improving sensitivity for detecting biological responses to low-level xenobiotic exposures, these data can aid the classification of novel contaminants based on the similarity of metabolomic responses to well-characterized "model" compounds. This approach is clearly of use for rapid, sensitive, and specific analyses of chemical effect on organisms, and can supplement existing methods, such as the Zebrafish Embryo Toxicity assay (OECD TG236), with molecular-level information.

Plasma proteome profiles of White Sucker (Catostomus commersonii) from the Athabasca River within the oil sands deposit

There are questions about the potential for oil sands related chemicals to enter the Athabasca River, whether from tailing ponds, atmospheric deposition, precipitation, or transport of mining dust, at concentrations sufficient to negatively impact the health of biota. We applied shotgun proteomics to generate protein profiles of mature male and female White Sucker (Catostomus commersonii) that were collected from various sites along the main stem of the Athabasca River in 2011 and 2012. On average, 399±131 (standard deviation) proteins were identified in fish plasma from each location in both years. Ingenuity Pathway Analysis software was used to determine the proteins' core functions and to compare the datasets by location, year, and sex. Principal component analysis (PCA) was used to determine if variation in the number of proteins related to a core function among all male and female individuals from both sampling years was affected by location. The core biological functions of plasma proteins that were common to both sampling years for males and females from each location were also estimated separately (based on Ingenuity's Knowledge Base). PCA revealed site-specific differences in the functional characteristics of the plasma proteome from white sucker sampled from downstream of oil sands extraction facilities compared with fish from upstream. Plasma proteins that were unique to fish downstream of oil sands extraction were related to lipid metabolism, small molecule biochemistry, vitamin and mineral metabolism, endocrine system disorders, skeletal and muscular development and function, neoplasia, carcinomas, and gastrointestinal disease.

TPhP exposure disturbs carbohydrate metabolism, lipid metabolism, and the DNA damage repair system in zebrafish liver

Triphenyl phosphate is a high production volume organophosphate flame retardant that has been detected in multiple environmental media at increasing concentrations. The environmental and health risks of triphenyl phosphate have drawn attention because of the multiplex toxicity of this chemical compound. However, few studies have paid close attention to the impacts of triphenyl phosphate on liver metabolism. We investigated hepatic histopathological, metabolomic and transcriptomic responses of zebrafish after exposure to 0.050 mg/L and 0.300 mg/L triphenyl phosphate for 7 days. Metabolomic analysis revealed significant changes in the contents of glucose, UDP-glucose, lactate, succinate, fumarate, choline, acetylcarnitine, and several fatty acids. Transcriptomic analysis revealed that related pathways, such as the glycosphingolipid biosynthesis, PPAR signaling pathway and fatty acid elongation, were significantly affected. These results suggest that triphenyl phosphate exposure markedly disturbs hepatic carbohydrate and lipid metabolism in zebrafish. Moreover, DNA replication, the cell cycle, and non-homologous end-joining and base excision repair were strongly affected, thus indicating that triphenyl phosphate hinders the DNA damage repair system in zebrafish liver cells. The present study provides a systematic analysis of the triphenyl phosphate-induced toxic effects in zebrafish liver and demonstrates that low concentrations of triphenyl phosphate affect normal metabolism and cell cycle.