Distinctive metabolite profiles in in-migrating Sockeye salmon suggest sex-linked endocrine perturbation

Contaminants of Concern and Spatiotemporal Metabolomic Changes in Quagga Mussels (Dreissena bugensis rostriformis) from the Milwaukee Estuary (Wisconsin, USA)

Environmental metabolomics has emerged as a promising technique in the field of biomonitoring and as an indicator of aquatic ecosystem health. In the Milwaukee Estuary (Wisconsin, USA), previous studies have used a nontargeted metabolomic approach to distinguish between zebra mussels (Dreissena polymorpha) collected from sites of varying contamination. To further elucidate the potential effects of contaminants on bivalve health in the Milwaukee Estuary, the present study adopted a caging approach to study the metabolome of quagga mussels (Dreissena bugensis rostriformis) deployed in six sites of varying contamination for 2, 5, or 55 days. Caged mussels were co-deployed with two types of passive sampler (polar organic chemical integrative samplers and semipermeable membrane devices) and data loggers. In conjunction, in situ quagga mussels were collected from the four sites studied previously and analyzed for residues of contaminants and metabolomics using a targeted approach. For the caging study, temporal differences in the metabolomic response were observed with few significant changes observed after 2 and 5 days, but larger differences (up to 97 significantly different metabolites) to the metabolome in all sites after 55 days. A suite of metabolic pathways were altered, including biosynthesis and metabolism of amino acids, and upmodulation of phospholipids at all sites, suggesting a potential biological influence such as gametogenesis. In the caging study, average temperatures appeared to have a greater effect on the metabolome than contaminants, despite a large concentration gradient in polycyclic aromatic hydrocarbons residues measured in passive samplers and mussel tissue. Conversely, significant differences between the metabolome of mussels collected in situ from all three contaminated sites and the offshore reference site were observed. Overall, these findings highlight the importance of contextualizing the effects of environmental conditions and reproductive processes on the metabolome of model organisms to facilitate the wider use of this technique for biomonitoring and environmental health assessments. Environ Toxicol Chem 2024;43:307-323. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Multi-omics responses in tree swallow (Tachycineta bicolor) nestlings from the Maumee Area of Concern, Maumee River, Ohio

A multi-omics approach was utilized to identify altered biological responses and functions, and to prioritize contaminants to assess the risks of chemical mixtures in the Maumee Area of Concern (AOC), Maumee River, OH, USA. The Maumee AOC is designated by the United States Environmental Protection Agency as having significant beneficial use impairments, including degradation of fish and wildlife populations, bird or animal deformities or reproduction problems, and loss of fish and wildlife habitat. Tree swallow (Tachycineta bicolor) nestlings were collected at five sites along the Maumee River, which included wastewater treatment plants (WWTPs) and industrial land-use sites. Polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), polycyclic aromatic hydrocarbons (PAHs), polychlorinated dibenzo p dioxins and furans (PCDD/Fs), and chlorinated pesticide concentrations were elevated in Maumee tree swallows, relative to a remote reference site, Star Lake, WI, USA. Liver tissue was utilized for non-targeted transcriptome and targeted metabolome evaluation. A significantly differentially expressed gene cluster related to a downregulation in cell growth and cell cycle regulation was identified when comparing all Maumee River sites with the reference site. There was an upregulation of lipogenesis genes, such as PPAR signaling (HMGCS2, SLC22A5), biosynthesis of unsaturated fatty acids (FASN, SCD, ELOVL2, and FADS2), and higher lipogenesis related metabolites, such as docosapentaenoic acid (DPA), docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA), and arachidonic acid (AA) at two industrial land-use sites, Ironhead and Maumee, relative to WWTP sites (Perrysburg and SideCut), and the reference site. Toledo Water, in the vicinity of the other two industrial sites and also adjacent to a WWTP, showed a mix of signals between industrial land-use and WWTP land-use. PAHs, oxychlordane, and PBDEs were determined to be the most likely causes of the differentiation in biological responses, including de novo lipogenesis and biosynthesis of unsaturated fatty acids.

Predicting bioconcentration factor and estrogen receptor bioactivity of bisphenol a and its analogues in adult zebrafish by directed message passing neural networks

The bioconcentration factor (BCF) is a key parameter for bioavailability assessment of environmental pollutants in regulatory frameworks. The comparative toxicology and mechanism of action of congeners are also of concern. However, there are limitations to acquire them by conducting field and laboratory experiments while machinelearning is emerging as a promising predictive tool to fill the gap. In this study, the Direct Message Passing Neural Network (DMPNN) was applied to predict logBCFs of bisphenol A (BPA) and its four analogues (bisphenol AF (BPAF), bisphenol B (BPB), bisphenol F (BPF) and bisphenol S (BPS)). For the test set, the Pearson correlation coefficient (PCC) and mean square error (MSE) were 0.85 and 0.52 respectively, suggesting a good predictive performance. The predicted logBCFs values by the DMPNN ranging from 0.35 (BPS) to 2.14 (BPAF) coincided well with those by the classical EPI Suite (BCFBAF model). Besides, estrogen receptor α (ERα) bioactivity of these bisphenols was also predicted well by the DMPNN, with a probability of 97.0 % (BPB) to 99.7 % (BPAF), which was validated by the extent of vitellogenin (VTG) induction in male zebrafish as a biomarker except BPS. Thus, with little need for expert knowledge, DMPNN is confirmed to be a useful tool to accurately predict logBCF and screen for estrogenic activity from molecular structures. Moreover, a gender difference was noted in the changes of three endpoints (logBCF, ER binding affinity and VTG levels), the rank order of which was BPAF > BPB > BPA > BPF > BPS consistently, and abnormal amino acid metabolism is featured as an omics signature of abnormal hormone protein expression.

A Multi-Matrix Metabolomic Approach in Ringed Seals and Beluga Whales to Evaluate Contaminant and Climate-Related Stressors

As a high trophic-level species, ringed seals (Pusa hispida) and beluga whales (Delphinapterus leucas) are particularly vulnerable to elevated concentrations of biomagnifying contaminants, such as polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs) and mercury (Hg). These species also face climate-change-related impacts which are leading to alterations in their diet and associated contaminant exposure. The metabolomic profile of marine mammal tissues and how it changes to environmental stressors is poorly understood. This study characterizes the profiles of 235 metabolites across plasma, liver, and inner and outer blubber in adult ringed seals and beluga whales and assesses how these profiles change as a consequence of contaminants and dietary changes. In both species, inner and outer blubber were characterized by a greater proportion of lipid classes, whereas the dominant metabolites in liver and plasma were amino acids, carbohydrates, biogenic amines and lysophosphatidylcholines. Several metabolite profiles in ringed seal plasma correlated with δ13C, while metabolite profiles in blubber were affected by hexabromobenzene in ringed seals and PBDEs and Hg in belugas. This study provides insight into inter-matrix similarities and differences across tissues and suggests that plasma and liver are more suitable for studying changes in diet, whereas liver and blubber are more suitable for studying the impacts of contaminants.

The plasma metabolome of Atlantic salmon as studied by 1H NMR spectroscopy using standard operating procedures: effect of aquaculture location and growth stage

INTRODUCTION

Metabolomics applications to the aquaculture research are increasing steadily. The use of standardized proton nuclear magnetic resonance (¹H NMR) spectroscopy can provide the aquaculture industry with an unbiased, reproducible, and high-throughput screening tool, which can help to diagnose nutritional and disease-related metabolic disorders in farmed fish.

OBJECTIVE

Standard operating procedures developed for analysing (human) plasma by ¹H NMR were applied to fingerprint the metabolome in plasma samples collected from Atlantic salmon. The aim was to explore the metabolome of salmon plasma in relation to growth stage and sampling site.

METHODS

A total of 72 salmon were collected from three aquaculture sites in Norway (Lat. 65, 67, and 70 °N) and over two sampling events (December 2017 and November 2018). Plasma drawn from each salmon was measured by ¹H NMR and metabolites were quantified using the SigMa software. The NMR data was analysed by principal component analysis (PCA) and ANOVA-simultaneous component analysis (ASCA).

RESULTS

Important metabolic differences were evidenced, with adult salmon having a much higher content of very low-density lipoproteins and cholesterol in their plasma, while smolts displayed significantly higher levels of propylene glycol. Overall, 24% of the metabolite variation was due to the growth stage, whereas 12% of the metabolite variation was related to the aquaculture site and practice (p < 0.001).

CONCLUSION

This study provides a baseline investigation of the plasma metabolome of the Atlantic salmon and demonstrates how ¹H NMR metabolomics can be used in future investigations for comparing aquaculture practices and their influence on the fish metabolome.

Lipid metabolites as indicators of body condition in highly contaminant-exposed belugas from the endangered St. Lawrence Estuary population (Canada)

The endangered St. Lawrence Estuary (SLE) beluga population is declining and has shown no sign of recovery over the past decades despite several protective measures. Changes in the availability of food resources and exposure to organohalogen contaminants have been suggested as potential factors limiting the recovery of this population. Studies on SLE belugas have suggested that contaminant exposure may perturb energy metabolism, however, whether this translates into changes in energy reserves (lipid composition) and body condition is unknown. The objective of this study was to investigate the relationships between body condition and concentrations of organohalogens (polychlorinated biphenyls, organochlorine pesticides, and flame retardants) and a range of lipid metabolites (fatty acids, acylcarnitines, lysophosphatidylcholines, phosphatidylcholines, and sphingomyelins) in blubber samples collected from 51 SLE beluga carcasses recovered between 1998 and 2016 for which the cause of mortality was documented. Blubber Σ₉fatty acid concentrations in SLE belugas significantly decreased between 1998 and 2016, suggesting a decline in energy reserves over the past two decades. Concentrations of several phosphatidylcholine analogues were greater in blubber of beluga males and/or females that were in poor body condition compared to those in good body condition. Moreover, concentrations of phosphatidylcholine acyl-alkyl C32:2 were greater in females that died from primary starvation (poor body condition). Greater concentrations of Σ₁₂emerging flame retardants were also found in blubber of SLE beluga females that were in poorer body condition. This study suggests that the use of membrane lipids including phosphatidylcholine concentrations may be a good indicator of body condition and energy reserve status in blubber of marine mammals.

Distribution and potential risk assessment of antibiotic pollution in the main drinking water sources of Nanjing, China

The distribution character of 41 antibiotics belonging to 6 groups, sulfonamides (SAs), quinolones (QUs), tetracyclines (TCs), macrolides (MLs), penicillins (PLs), and chloramphenicol (CHL), was investigated in drinking water sources of Nanjing during 2017-2019. MLs (42.98%) were the most abundant category, followed by SAs (25.94%) and QUs (22.52%). The dominant antibiotic was ofloxacin (OFX) in Dec. 2017 (average concentration, 3.14 ng/L; range, ND-35.20 ng/L) and Nov. 2018 (2.16 ng/L, ND-12.26 ng/L), and sulfadiazine (SDZ) in Mar. 2019 (16.37 ng/L, ND-25.90 ng/L). For Dec. 2017, the total concentrations in Zhongshan Waterworks (S15) and Jinniushan Reservoir (S16) were significantly higher than the other sampling sites, which may be attributed to point source pollution. The ecological and human risk of the main antibiotics was assessed by risk quotients (RQs) and target hazard quotient (THQ), respectively. Most of the RQ values were below 0.1, except enrofloxacin (ERX, 0.11) and enoxacin (ENX, 0.62) in Dec. 2017, lomefloxacin (LOM, 0.14) in Nov. 2018, and LOM (0.28) and ERX (0.10) in Mar. 2019. This indicated that the risk of the target antibiotics to aquatic organisms in the 3 years was moderate or low level. Meanwhile, results of the THQ values showed that antibiotic exposure caused no risk to human health. This research provides scientific information for antibiotic pollution control and enriches environmental monitoring data in the drinking water sources.

Metabolomic profiles of the endangered St. Lawrence Estuary beluga population and associations with organohalogen contaminants

The endangered beluga (Delphinapterus leucas) population residing in the St. Lawrence Estuary (SLE; Eastern Canada) is declining. The elevated tissue concentrations of a wide range of organohalogen contaminants might play a role in the non-recovery of this whale population. Organohalogens have been reported to impair the regulation of several metabolic products from cellular reactions in mammals such as amino acids and fatty acids. The objective of this study was to investigate a suite of organohalogens including polychlorinated biphenyls, organochlorine pesticides, short-chain chlorinated paraffins (SCCPs), polybrominated diphenyl ethers, and selected emerging flame retardants in blubber (biopsy) collected from 40 SLE male belugas, and their relationships to skin concentrations of targeted metabolites (i.e., 21 amino acids, 22 biogenic amines, 18 fatty acids, and 17 energy metabolites). A cluster analysis based on metabolomic profiles distinguished two main subgroups of belugas in the upper and lower sector of their summer habitat in the SLE. These results indicate that ecological factors such as local prey availability and diet composition played a role in shaping the metabolite profiles of belugas. Moreover, SCCP concentrations in SLE male belugas correlated negatively with those of four unsaturated fatty acids (C16:1ω7, C22:5ω3c1, C22:5ω3c2, and C22:6ω3), and positively with those of acetylornithine (biogenic amine). These findings suggest that biological functions such as lipid metabolism represent potential targets for organohalogens in this population, and further our understanding on potential health risks associated with elevated organohalogen exposure in cetaceans. Our results also underscore the necessity of considering ecological factors (e.g., diet and habitat use) in metabolomic studies.

Environmental exposure of northern pike to a primary wastewater effluent: Impact on the lipidomic profile and lipid metabolism

Lipids play important roles in growth, reproduction, locomotion, and migration of fish. Municipal effluents, which are complex mixtures of biological and chemical compounds including flame retardants, have been shown to alter lipid metabletabolism in environmentally and experimentally exposed fish. Down-regulation of several genes coding for fatty acid metabolism enzymes has previously been reported in male northern pike (Esox lucius) collected in the St. Lawrence River (QC, Canada) downstream of a major primary wastewater treatment plant (WWTP) point of discharge. The main objective of this study was to investigate the effects of exposure to the Montreal's WWTP effluent on the lipidomic profile (i.e., fatty acids, acylcarnitines, and phospholipids) as well as the transcription of genes related to lipid metabolism in the liver of northern pike collected upstream and downstream of this WWTP effluent. Halogenated flame retardant concentrations were also determined in pike liver and used as markers of exposure to this effluent. Greater concentrations of saturated and monounsaturated lysophosphatidylcholines (LPCs) and lower concentrations of polyunsaturated LPCs were determined in the liver of pike collected downstream of the WWTP compared to those collected upstream. Lower mRNA levels of peroxisome proliferator-activated receptor alpha (pparα), a major regulator of lipid metabolism, were also measured in pike exposed to Montreal's WWTP effluent. In addition, the relative contributions (%) of LPC 18:2 and LPC14:0, compounds used as markers of obesity and inflammation, were significantly correlated with halogenated flame retardant concentrations and fish girth. Results of the present study suggest that chronic environmental exposure to a primary WWTP effluent can modulate the transcription of genes related to lipid metabolism, and hence affect the hepatic phospholipid composition of pike from the St. Lawrence River.

Responses in the crucian carp (Carassius auratus) exposed to environmentally relevant concentration of 17α-Ethinylestradiol based on metabolomics

17α-ethynylestradiol (EE2), a ubiquitous synthetic endocrine disrupting chemical, was the principal component of contraceptive drugs and one of common hormone medications. The detrimental impact of EE2 on the reproduction of organisms was widely recognized. However, the underlying mechanisms of physiological and metabolome effects of EE2 on freshwater fish are still unclear. This study investigated the toxic effects and related mechanisms of EE2 on freshwater fish crucian carp (Carassius auratus) based on metabolomics. Crucian carp were exposed to EE2 at environmentally relevant concentration for 9 days, 18 days, and 27 days, and the biological responses were explored through analysis of the physiological endpoints, steroid hormones, and metabolome. The physiological endpoints of crucian carp had no distinct change after EE2 exposure. However, metabolomics analysis probed significant deviation based on chemometrics, indicating that the metabolomics approach was more sensitive to the effects of EE2 at environmentally relevant concentration to freshwater fish than the traditional endpoints. The alterations of 24 metabolites in gonad and 16 metabolites in kidney were induced by treatment with EE2, respectively, which suggesting the perturbations in amino acid metabolism, lipid metabolism, energy metabolism, and oxidative stress. Moreover, EE2 exposure could induce the disruption of lipid metabolism and then broke the homeostasis of endogenous steroid hormones. Metabolomics provided a new strategy for the studies on contaminant exposure at a low dose in the short term and gave important information for the toxicology and mechanism of EE2.

Concentrations of legacy and new contaminants are related to metabolite profiles in Hudson Bay polar bears

There are very few metabolomics assessments based on field accumulated, uncontrolled contaminant exposures in wildlife, particularly in the Arctic. In the present study, targeted metabolomics and contaminant data were analyzed together to assess potential influences of contaminant exposure on the hepatic metabolome of male polar bears (n = 29) from the southern and western Hudson Bay (SHB and WHB respectively), Canada. The 29 metabolites identified as important in the differentiation of the two subpopulations after partial least squares discriminant analysis (PLS-DA) included phosphatidylcholines (PCs), acylcarnitines (ACs; involved in β-oxidation of fatty acids), and the fatty acid (FA) arachidonic acid (ARA). Perfluorinated alkyl substances, polybrominated diphenyl ethers, dichlorodiphenyldichloroethylene (p,p'-DDE) and some highly chlorinated ortho-polychlorinated biphenyl congeners were greater in the SHB bears and were consistently inversely correlated with discriminating ACs and PCs between the subpopulations. The concentrations of discriminatory, legacy organochlorine pesticides along with one tetrachlorobiphenyl were greater in the WHB and were directly correlated with the VIP-identified ACs and PCs. ARA, glycerophospholipid and several amino acid metabolic pathways were identified as different between subpopulations and/or were impacted. ARA is an important, conditionally essential, dietary n-6 FA and is also part of the inflammation response, and elevated concentrations in the SHB could be related to differences in chronic contaminant exposure and/or differences in diet and/or season, among a number of possible explanations. Dietary tracers (stable isotopes of carbon and nitrogen) were correlated with some discriminatory metabolites, supporting the hypothesis that dietary variation was also an important factor in the differentiation of the subpopulations. The results suggest linkages between contaminant exposure in Hudson Bay polar bears and elements of the hepatic metabolome, particularly those related to lipid metabolism.

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.

Omics in Zebrafish Teratogenesis

The genome revolution represents a complete change on our view of biological systems. The quantitative determination of changes in all major molecular components of the living cells, the "omics" approach, opened whole new fields for all health sciences. Genomics, transcriptomics, proteomics, metabolomics, and others, together with appropriate prediction and modeling tools, will mark the future of developmental toxicity assessment both for wildlife and humans. This is especially true for disciplines, like teratology, which rely on studies in model organisms, as studies at lower levels of organization are difficult to implement. Rodents and frogs have been the favorite models for studying human reproductive and developmental disorders for decades. Recently, the study of the development of zebrafish embryos (ZE) is becoming a major alternative tool to adult animal testing. ZE intrinsic characteristics makes this model a unique system to analyze in vivo developmental alterations that only can be studied applying in toto approaches. Moreover, under actual legislations, ZE is considered as a replacement model (and therefore, excluded from animal welfare regulations) during the first 5 days after fertilization. Here we review the most important components of the zebrafish toolbox available for analyzing early stages of embryotoxic events that could eventually lead to teratogenesis.

Reduced anxiety is associated with the accumulation of six serotonin reuptake inhibitors in wastewater treatment effluent exposed goldfish Carassius auratus

Pharmaceuticals and personal care products (PPCPs) have been found in wastewater treatment plant (WWTP) effluents and their recipient watersheds. To assess the potential of WWTP effluents to alter fish behaviour, we caged male goldfish (Carassius auratus) for 21-days at three sites along a contamination gradient downstream from a WWTP which discharges into Cootes Paradise Marsh, on the western tip of Lake Ontario. We also included a fourth caging site as an external reference site within Lake Ontario at the Jordan Harbour Conservation Area. We then measured concentrations of PPCPs and monoamine neurotransmitters in caged goldfish plasma, and conducted behavioural assays measuring activity, startle response, and feeding. We detected fifteen different PPCPs in goldfish plasma including six serotonin reuptake inhibitors (amitriptyline, citalopram, fluoxetine/norfluoxetine, sertraline, venlafaxine, and diphenhydramine). Plasma concentrations of serotonin were significantly greater in plasma of fish caged closer to the WWTP effluent outfall site. The fish caged near and downstream of the WWTP effluent were bolder, more exploratory, and more active overall than fish caged at the reference site. Taken together, our results suggest that fish downstream of WWTPs are accumulating PPCPs at levels sufficient to alter neurotransmitter concentrations and to also impair ecologically-relevant behaviours.

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.

Metabolic disruption of zebrafish (Danio rerio) embryos by bisphenol A. An integrated metabolomic and transcriptomic approach

Although bisphenol A (BPA) is commonly recognized as an endocrine disruptor, the metabolic consequences of its exposure are still poorly understood. In this study, we present a non-targeted LC-MS based metabolomic analysis in combination with a full-genome, high-throughput RNA sequencing (RNA-Seq) to reveal the metabolic effects and the subjacent regulatory pathways of exposing zebrafish embryos to BPA during the first 120 hours post-fertilization. We applied multivariate data analysis methods to extract biochemical information from the LC-MS and RNA-Seq complex datasets and to perform testable predictions of the phenotypic adverse effects. Metabolomic and transcriptomic data revealed a similar subset of altered pathways, despite the large difference in the number of identified biomarkers (around 50 metabolites and more than 1000 genes). These results suggest that even a moderate coverage of zebrafish metabolome may be representative of the global metabolic changes. These multi-omic responses indicate a specific metabolic disruption by BPA affecting different signaling pathways, such as retinoid and prostaglandin metabolism. The combination of transcriptomic and metabolomic data allowed a dynamic interpretation of the results that could not be drawn from either single dataset. These results illustrate the utility of -omic integrative analyses for characterizing the physiological effects of toxicants beyond the mere indication of the affected pathways.

Bisphenol analogues in surface water and sediment from the shallow Chinese freshwater lakes: Occurrence, distribution, source apportionment, and ecological and human health risk

Compared to Bisphenol A (BPA), current knowledge on the spatial distribution, potential sources and environmental risk assessment of other bisphenol analogues (BPs) remains limited. The occurrence, distribution and sources of seven BPs were investigated in the surface water and sediment from Taihu Lake and Luoma Lake, which are the Chinese shallow freshwater lakes. Because there are many industries and living areas around Taihu Lake, the total concentrations of ∑BPs were much higher than that in Luoma Lake, which is away from the industry-intensive areas. For the two lakes, BPA was still the dominant BPs in both surface water and sediment, followed by BPF and BPS. The spatial distribution and principal component analysis showed that BPs in Luoma Lake was relatively homogeneous and the potential sources were relatively simple than that in Taihu Lake. The spatial distribution of BPs in sediment of Taihu Lake indicated that ∑BPs positively correlated with the TOC content. For both Taihu Lake and Luoma Lake, the risk assessment at the sampling sites showed that no high risk in surface water and sediment (RQ_t < 1.0, and EEQ_t < 1.0 ng E₂/L).

Differential metabolic profiles associated to movement behaviour of stream-resident brown trout (Salmo trutta)

The mechanisms that can contribute in the fish movement strategies and the associated behaviour can be complex and related to the physiology, genetic and ecology of each species. In the case of the brown trout (Salmo trutta), in recent research works, individual differences in mobility have been observed in a population living in a high mountain river reach (Pyrenees, NE Spain). The population is mostly sedentary but a small percentage of individuals exhibit a mobile behavior, mainly upstream movements. Metabolomics can reflect changes in the physiological process and can determine different profiles depending on behaviour. Here, a non-targeted metabolomics approach was used to find possible changes in the blood metabolomic profile of S. trutta related to its movement behaviour, using a minimally invasive sampling. Results showed a differentiation in the metabolomic profiles of the trouts and different level concentrations of some metabolites (e.g. cortisol) according to the home range classification (pattern of movements: sedentary or mobile). The change in metabolomic profiles can generally occur during the upstream movement and probably reflects the changes in metabolite profile from the non-mobile season to mobile season. This study reveals the contribution of the metabolomic analyses to better understand the behaviour of organisms.

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.

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.

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

There are multiple sources of biological and technical variation in a typical ecotoxicology study that may not be revealed by traditional endpoints but that become apparent in an omics dataset. As researchers increasingly apply omics technologies to environmental studies, it will be necessary to understand and control the main source(s) of variability to facilitate meaningful interpretation of such data. For instance, can variability in omics studies be addressed by changing the approach to study design and data analysis? Are there statistical methods that can be employed to correctly interpret omics data and make use of unattributed, inherent variability? The present study presents a review of experimental design and statistical considerations applicable to the use of omics methods in systems toxicology studies. In addition to highlighting potential sources that contribute to experimental variability, this review suggests strategies with which to reduce and/or control such variability so as to improve reliability, reproducibility, and ultimately the application of omics data for systems toxicology.