Defining the role of omics in assessing ecosystem health: Perspectives from the Canadian environmental monitoring program

Landscape transcriptomic analysis detects thermal stress responses and potential adaptive variation in wild brook trout (Salvelinus fontinalis) during successive heatwaves

Extreme weather events, such as heatwaves, are becoming more frequent and intense as a result of climate change. Critically, such extreme weather events can be more important drivers of extirpation and selection than changes in annual or seasonal averages and they pose a particularly large threat to poikilothermic organisms. In this study, we evaluated the transcriptomic response of a coldwater adapted fish species, the eastern brook trout (Salvelinus fontinalis), to two successive heatwaves during July and August 2022. We sampled brook trout at eight time points from four streams (N = 116 fish), sequenced mRNA from gill samples using TagSeq, and quantified expression levels of 32,670 unique transcripts. Multivariate analyses found that overall expression patterns in response to water temperature change were similar among streams. These analyses further detected groups of genes involved in immune response and oxygen carrier activity that were upregulated and downregulated respectively at higher water temperatures. We also detected 43 genes that were differentially expressed at different time points and followed the same expression pattern during the two heatwaves. Of these genes, 42 covaried with water temperature and most (27, 62.8 %) exhibited responses that varied by stream. Some of the differentially expressed genes, including heat shock proteins and cold-inducible RNA binding proteins, have been widely linked to temperature responses in experimental studies, whereas other genes we identified have functions that have not been well-studied in relationship to temperature or have unknown functions. This study shows the utility of landscape transcriptomic approaches to identify important biological processes governing wild organismal responses to short-term stressors. The results of this study can guide future investigations to identify phenotypic and genetic diversity that contribute to adaptive responses to heatwaves and improve predictions of how populations will respond to future climate change.

Effects of the organochlorine pesticide metabolite p,p'-DDE on the gastrointestinal lipidome in fish: A novel toxicity pathway for a legacy pollutant

Though phased out from use in the United States, environmental contamination by organochlorine pesticides (OCPs) remains a widespread issue, especially around intensive agricultural regions. OCPs, such as dichlorodiphenyltrichloroethane (DDT) and its primary metabolite, dichlorodiphenyldichloroethylene (DDE), have been detected in soils, sediments, surface waters, and biota decades after their discontinued use. As OCPs are persistent and can bioaccumulate in fats, these compounds can transfer and magnify across food webs. Freshwater predatory fish and birds can accumulate high OCP concentrations, leading to a myriad of deleterious impacts on organismal health. Studies have found evidence of reproductive disruption in predatory fish, such as the largemouth bass (LMB; Micropterus salmoides), associated with DDT and DDE exposure. DDE can act through estrogenic pathways and induce the expression of estrogenic signals in male animals; however, the molecular mechanism of disruption is unclear. Recently, metabolomics research has revealed corollary relationships between lipid signals and organic pollutant toxicity. Here, a two-month feeding experiment on LMB was conducted to assess the interactions of DDE (as p,p'-DDE) in food with gut and liver lipid signaling. Targeted lipidomic analysis revealed global alterations in the abundance of tissue lipids, especially cholesteryl esters and phospholipids, in LMB exposed to low levels of p,p'-DDE. Results from these studies indicate that p,p'-DDE may act through disruption of normal lipid homeostasis to cause toxicity in freshwater fish.

Environmental signatures and fish proteomics: A multidisciplinary study to identify the major stressors in estuaries located in French agricultural watersheds

Watersheds and estuaries are impacted by multiple anthropogenic stressors that affect their biodiversity and functioning. Assessing their ecological quality has consequently remained challenging for scientists and stakeholders. In this paper, we propose a multidisciplinary approach to identify the stressors in seven small French estuaries located in agricultural watersheds. We collected data from landscape (geography, hydrobiology) to estuary (pollutant chemistry) and fish individual scales (environmental signatures, proteomics). This integrative approach focused on the whole hydrosystems, from river basins to estuaries. To characterize each watershed, we attempted to determine the land use considering geographic indicators (agricultural and urbanised surfaces) and landscape patterns (hedges density and riparian vegetation). Juveniles of European flounder (Platichthys flesus) were captured in September, after an average residence of five summer months in the estuary. Analyses of water, sediments and biota allowed to determine the concentrations of dissolved inorganic nitrogen species, pesticides and trace elements in the systems. Environmental signatures were also measured in flounder tissues. These environmental parameters were used to establish a typology of the watersheds. Furthermore, data from proteomics on fish liver were combined with environmental signatures to determine the responses of fish to stressors in their environments. Differential protein abundances highlighted a dysregulation related to the detoxification of xenobiotics (mainly pesticides) in agricultural watersheds, characterized by intensive cereal and vegetable crops and high livestock. Omics also revealed a dysregulation of proteins associated with the response to hypoxia and heat stress in some estuaries. Furthermore, we highlighted a dysregulation of proteins involved in urea cycle, immunity and metabolism of fatty acids in several systems. Finally, the combination of environmental and molecular signatures appears to be a relevant method to identify the major stressors operating within hydrosystems.

Network-based integration of omics, physiological and environmental data in real-world Elbe estuarine Zander

Coastal and estuarine environments are under endogenic and exogenic pressures jeopardizing survival and diversity of inhabiting biota. Information of possible synergistic effects of multiple (a)biotic stressors and holobiont interaction are largely missing in estuaries like the Elbe but are of importance to estimate unforeseen effects on animals' physiology. Here, we seek to leverage host-transcriptional RNA-seq and gill mucus microbial 16S rRNA metabarcoding data coupled with physiological and abiotic measurements in a network analysis approach to decipher the impact of multiple stressors on the health of juvenile Sander lucioperca along one of the largest European estuaries. We find mesohaline areas characterized by gill tissue specific transcriptional responses matching osmosensing and tissue remodeling. Liver transcriptomes instead emphasized that zander from highly turbid areas were undergoing starvation which was supported by compromised body condition. Potential pathogenic bacteria, including Shewanella, Acinetobacter, Aeromonas and Chryseobacterium, dominated the gill microbiome along the freshwater transition and oxygen minimum zone. Their occurrence coincided with a strong adaptive and innate transcriptional immune response in host gill and enhanced energy demand in liver tissue supporting their potential pathogenicity. Taken together, we show physiological responses of a fish species and its microbiome to abiotic factors whose impact is expected to increase with consequences of climate change. We further present a method for the close-meshed detection of the main stressors and bacterial species with disease potential in a highly productive ecosystem.

Time Is Ripe for Targeting Per- and Polyfluoroalkyl Substances-Induced Hormesis: Global Aquatic Hotspots and Implications for Ecological Risk Assessment

Globally implemented ecological risk assessment (ERA) guidelines marginalize hormesis, a biphasic dose-response relationship characterized by low-dose stimulation and high-dose inhibition. The present study illuminated the promise of hormesis as a scientific dose-response model for ERA of per- and polyfluoroalkyl substances (PFAS) represented by perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS). A total of 266 hormetic dose-response relationships were recompiled from 1237 observations, covering 30 species from nine representative taxonomic groups. The standardized hormetic amplitudes followed the log-normal probability distribution, being subject to the limits of biological plasticity but independent of stress inducers. The SHapley Additive exPlanations algorithm revealed that the target endpoint was the most important variable explaining the hormetic amplitudes. Subsequently, quantitative frameworks were established to incorporate hormesis into the predicted no-effect concentration levels, with a lower induction dose and a zero-equivalent point but a broader hormetic zone for PFOS. Realistically, 10,117 observed concentrations of PFOA and PFOS were gathered worldwide, 4% of which fell within hormetic zones, highlighting the environmental relevance of hormesis. Additionally, the hormesis induction potential was identified in other legacy and emerging PFAS as well as their alternatives and mixtures. Collectively, it is time to incorporate the hormesis concept into PFAS studies to facilitate more realistic risk characterizations.

Metabolic disruptions and impaired reproductive fitness in wild-caught freshwater turtles (Emydura macquarii macquarii) exposed to elevated per- and polyfluoroalkyl substances (PFAS)

Per- and poly-fluoroalkyl substances (PFAS) pose a threat to organisms and ecosystems due to their persistent nature. Ecotoxicology endpoints used in regulatory guidelines may not reflect multiple, low-level but persistent stressors. This study examines the biological effects of PFAS on Eastern short-necked turtles in Queensland, Australia. In this study, blood samples were collected and analysed for PFAS, hormone levels, and functional omics endpoints. High levels of PFAS were found in turtles at the impacted site, with PFOS being the dominant constituent. The PFAS profiles of males and females differed, with males having higher PFAS concentrations. Hormone concentrations differed between impacted and reference sites in male turtles, with elevated testosterone and corticosterone indicative of stress. Further, energy utilisation, nucleotide synthesis, nitrogen metabolism, and amino acid synthesis were altered in both male and female turtles from PFAS-impacted sites. Both sexes show similar metabolic responses to environmental stressors from the PFAS-contaminated site, which may adversely affect their reproductive fitness. Purine metabolism, caffeine metabolism, and ferroptosis pathway changes in turtles can cause gout, cell death, and overall health problems. Further, the study showed that prolonged exposure to elevated PFAS levels in the wild could compromise turtle reproductive fitness by disrupting reproductive steroids and metabolic pathways.

Advancing environmental monitoring across the water continuum combining biomarker analysis in multiple sentinel species: A case study in the Seine-Normandie Basin (France)

Nowadays, biomarkers are recognized as valuable tools to complement chemical and ecological assessments in biomonitoring programs. They provide insights into the effects of contaminant exposures on individuals and establish connections between environmental pressure and biological response at higher levels. In the last decade, strong improvements in the design of experimental protocols and the result interpretation facilitated the use of biomarker across wide geographical areas, including aquatic continua. Notably, the statistical establishment of reference values and thresholds enabled the discrimination of contamination effects in environmental conditions, allowed interspecies comparisons, and eliminated the need of a reference site. The aim of this work was to study freshwater-estuarine-coastal water continua by applying biomarker measurements in multi-species caged organisms. During two campaigns, eight sentinel species, encompassing fish, mollusks, and crustaceans, were deployed to cover 25 sites from rivers to the sea. As much as possible, a common methodology was employed for biomarker measurements (DNA damage and phagocytosis efficiency) and data interpretation based on guidelines established using reference values and induction/inhibition thresholds (establishment of three effect levels). The methodology was successfully implemented and allowed us to assess the environmental quality. Employing multiple species per site enhances confidence in observed trends. The results highlight the feasibility of integrating biomarker-based environmental monitoring programs across a continuum scale. Biomarker results align with Water Framework Directive indicators in cases of poor site quality. Additionally, when discrepancies arise between chemical and ecological statuses, biomarker findings offer a comprehensive perspective to elucidate the disparities. Presented as a pilot project, this work contributes to gain insights into current biomonitoring needs, providing new questions and perspectives.

Review: Fish bile, a highly versatile biomarker for different environmental pollutants

Ecotoxicological assessments encompass a broad spectrum of biochemical endpoints and ecological factors, allowing for comprehensive assessments concerning pollutant exposure levels and their effects on both fish populations and surrounding ecosystems. While these evaluations offer invaluable insights into the overall health and dynamics of aquatic environments, they often provide an integrated perspective, making it challenging to pinpoint the precise sources and individual-level responses to environmental contaminants. In contrast, biliary pollutant excretion assessments represent a focused approach aimed at understanding how fish at the individual level respond to environmental stressors. In this sense, the analysis of pollutant profiles in fish bile not only serves as a valuable exposure indicator, but also provides critical information concerning the uptake, metabolism, and elimination of specific contaminants. Therefore, by investigating unique and dynamic fish responses to various pollutants, biliary assessments can contribute significantly to the refinement of ecotoxicological studies. This review aims to discuss the multifaceted utility of bile as a potent biomarker for various environmental pollutants in fish in targeted monitoring strategies, such as polycyclic aromatic hydrocarbons, metals, pesticides, pharmaceuticals, estrogenic compounds, resin acids, hepatotoxins and per- and polyfluorinated substances. The main caveats of this type of assessment are also discussed, as well as future directions of fish bile studies.

Reporting and reproducibility: Proteomics of fish models in environmental toxicology and ecotoxicology

Environmental toxicology and ecotoxicology research efforts are employing proteomics with fish models as New Approach Methodologies, along with in silico, in vitro and other omics techniques to elucidate hazards of toxicants and toxins. We performed a critical review of toxicology studies with fish models using proteomics and reported fundamental parameters across experimental design, sample preparation, mass spectrometry, and bioinformatics of fish, which represent alternative vertebrate models in environmental toxicology, and routinely studied animals in ecotoxicology. We observed inconsistencies in reporting and methodologies among experimental designs, sample preparations, data acquisitions and bioinformatics, which can affect reproducibility of experimental results. We identified a distinct need to develop reporting guidelines for proteomics use in environmental toxicology and ecotoxicology, increased QA/QC throughout studies, and method optimization with an emphasis on reducing inconsistencies among studies. Several recommendations are offered as logical steps to advance development and application of this emerging research area to understand chemical hazards to public health and the environment.

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.

Multi-omic responses of fish exposed to complex chemical mixtures in the Shenandoah River watershed

To evaluate relationships between different anthropogenic impacts, contaminant occurrence, and fish health, we conducted in situ fish exposures across the Shenandoah River watershed at five sites with different land use. Exposure water was analyzed for over 500 chemical constituents, and organismal, metabolomic, and transcriptomic endpoints were measured in fathead minnows. Adverse reproductive outcomes were observed in fish exposed in the upper watershed at both wastewater treatment plant (WWTP) effluent- and agriculture-impacted sites, including decreased gonadosomatic index and altered secondary sex characteristics. This was accompanied with increased mortality at the site most impacted by agricultural activities. Molecular biomarkers of estrogen exposure were unchanged and consistent with low or non-detectable concentrations of common estrogens, indicating that alternative mechanisms were involved in organismal adverse outcomes. Hepatic metabolomic and transcriptomic profiles were altered in a site-specific manner, consistent with variation in land use and contaminant profiles. Integrated biomarker response data were useful for evaluating mechanistic linkages between contaminants and adverse outcomes, suggesting that reproductive endocrine disruption, altered lipid processes, and immunosuppression may have been involved in these organismal impacts. This study demonstrated linkages between human-impact, contaminant occurrence, and exposure effects in the Shenandoah River watershed and showed increased risk of adverse outcomes in fathead minnows exposed to complex mixtures at sites impacted by municipal wastewater discharges and agricultural practices.

Expanding non-invasive approaches for fish-health monitoring: A survey of the epidermal mucous metabolomes of phylogenetically diverse freshwater fish species

There is a pressing need for more-holistic approaches to fisheries assessments along with growing demand to reduce the health impacts of sample collections. Metabolomic tools enable the use of sample matrices that can be collected with minimal impact on the organism (e.g., blood, urine, and mucus) and provide high-throughput, untargeted biochemical information without the requirement of a sequenced genome. These qualities make metabolomics ideal for monitoring a wide range of fish species, particularly those under protected status. In the current study, we surveyed the relative abundances of 120 endogenous metabolites in epidermal mucus across eight freshwater fish species belonging to seven phylogenetic orders. Principal component analysis was used to provide an overview of the data set, revealing strong interspecies relationships in the epidermal mucous metabolome. Normalized relative abundances of individual endogenous metabolites were then used to identify commonalities across multiple species, as well as those metabolites that showed notable species specificity. For example, taurine was measured in high relative abundance in the epidermal mucus of common carp (Cyprinus carpio), northern pike (Esox lucius), golden shiner (Notemigonus crysoleucas), rainbow trout (Oncorhynchus mykiss), and rainbow smelt (Osmerus mordax), whereas γ-amino butyric acid (GABA) exhibited a uniquely high relative abundance in flathead catfish (Pylodictis olivaris). Finally, hierarchical cluster analysis was used to evaluate species relatedness as characterized by both the epidermal mucous metabolome (phenotype) and genetic phylogeny (genotype). This comparison revealed species for which relatedness in the epidermal mucous metabolome composition closely aligns with phylogenetic relatedness (e.g., N. crysoleucas and C. carpio), as well as species for which these two measures are not well aligned (e.g., P. olivaris and Polyodon spathula). These, and other findings reported here, highlight novel areas for future research with fish, including development of epidermal mucous-based markers for non-invasive health monitoring, sex determination, and hypoxia tolerance.

Pesticide pestilence: Global scenario and recent advances in detection and degradation methods

Increased anthropogenic activities are confronted as the main cause for rising environmental and health concerns globally, presenting an indisputable threat to both environment and human well-being. Modern-day industrialization has given rise to a cascade of concurrent environmental and health challenges. The global human population is growing at an alarming rate, posing tremendous pressure on future food security, and healthy and environmentally sustainable diets for all. To feed all, the global food production needs to increase by 50% by 2050, but this increase has to occur from the limited arable land, and under the present-day climate variabilities. Pesticides have become an integral component of contemporary agricultural system, safeguarding crops from pests and diseases and their use must be reduce to fulfill the SDG (Sustainable Development Goals) agenda . However, their indiscriminate use, lengthy half-lives, and high persistence in soil and aquatic ecosystems have impacted global sustainability, overshot the planetary boundaries and damaged the pure sources of life with severe and negative impacts on environmental and human health. Here in this review, we have provided an overview of the background of pesticide use and pollution status and action strategies of top pesticide-using nations. Additionally, we have summarized biosensor-based methodologies for the rapid detection of pesticide residue. Finally, omics-based approaches and their role in pesticide mitigation and sustainable development have been discussed qualitatively. The main aim of this review is to provide the scientific facts for pesticide management and application and to provide a clean, green, and sustainable environment for future generations.

Integration of environmental signatures and omics-based approaches on the European flounder to assist with health assessment of estuarine ecosystems in Brittany, France

This study aimed to develop a multidisciplinary approach to assess the ecological status of six moderate-sized French estuaries. For each estuary, we gathered geographical information, hydrobiological data, chemistry of pollutants and fish biology, including integration of proteomics and transcriptomics data. This integrative study covered the entire hydrological system studied, from the watershed to the estuary, and considered all the anthropogenic factors that can impact this environment. To reach this goal, European flounder (Platichthys flesus) were collected from six estuaries in September, which ensures a minimum residence time of five months within an estuary. Geographical metrics are used to characterize land use in each watershed. The concentrations of nitrite, nitrate, organic pollutants, and trace elements were measured in water, sediments and biota. All of these environmental parameters allowed to set up a typology of estuaries. Classical fish biomarkers, coupled with molecular data from transcriptomics and shotgun proteomics, highlighted the flounder's responses to stressors in its environment. We analysed the protein abundances and gene expression levels in the liver of fish from the different estuaries. We showed clear positive deregulation of proteins associated with xenobiotic detoxification in a system characterized by a large population density and industrial activity, as well as in a predominantly agricultural catchment area (mostly cultures of vegetables and pig breeding) mainly impacted by pesticides. Fish from the latter estuary also displayed strong deregulation of the urea cycle, most probably related to high nitrogen load. Proteomic and transcriptomic data also revealed a deregulation of proteins and genes related to the response to hypoxia, and a probable endocrine disruption in some estuaries. Coupling these data allowed the precise identification of the main stressors interacting within each hydrosystem.

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.

Transcriptional Responses as Biomarkers of General Toxicity: A Systematic Review and Meta-analysis on Metal-Exposed Bivalves

Through a systematic review and a series of meta-analyses, we evaluated the general responsiveness of putative transcriptional biomarkers of general toxicity and chemical stress. We targeted metal exposures performed on bivalves under controlled laboratory conditions and selected six transcripts associated with general toxicity for evaluation: catalase, glutathione-S-transferase, heat shock proteins 70 and 90, metallothionein, and superoxide dismutase. Transcriptional responses (n = 396) were extracted from published scientific articles (k = 22) and converted to log response ratios (lnRRs). By estimating toxic units, we normalized different metal exposures to a common scale, as a proxy of concentration. Using Bayesian hierarchical random effect models, we then tested the effects of metal exposure on lnRR, both for metal exposure in general and in meta-regressions using toxic unit and exposure time as independent variables. Corresponding analyses were also repeated with transcript and tissue as additional moderators. Observed patterns were similar for general and for transcript- and tissue-specific responses. The expected overall response to arbitrary metal exposure was an lnRR of 0.50, corresponding to a 65% increase relative to a nonexposed control. However, when accounting for publication bias, the estimated "true" response showed no such effect. Furthermore, expected response magnitude increased slightly with exposure time, but there was little support for general monotonic concentration dependence with regard to toxic unit. Altogether, the present study reveals potential limitations that need consideration prior to applying the selected transcripts as biomarkers in environmental risk assessment. Environ Toxicol Chem 2023;42:628-641. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Unravelling the Role of Metabolites for Detecting Physiological State of Wild Animals: European Rabbit's (Oryctolagus cuniculus) Case

European wild rabbit (Oryctolagus cuniculus) has been defined as a keystone species in the Mediterranean ecosystem. Rabbits have been classed as "endangered" by the IUCN within their native range. In this sense, animal nutrition may play a fundamental and limiting role in the conservation of wild species. The overarching goal of ecological nutrition is to unravel the extensive web of nutritional links that direct animals in their interactions with their ecological environments. The main aim of this work was to evaluate the effect of different feed intake, geographic location, animal sex, and reproductive stage on glucose, non-esterified fatty acids (NEFA), and plasmatic urea nitrogen (PUN), albumin, glutamate, and total protein metabolites. Additionally, we examined the potential of these metabolites as biomarkers. Full stomach contents and blood samples were collected from European wild rabbits (n = 89) for the analysis of the metabolites described above. Our work shows that the levels of these metabolites are affected by the sex of the animals, as well as by their reproductive stage (glucose, NEFA and albumin). There were signs of better optimisation of resources by females than by other groups of animals. These data may be interesting in the study of nutritional components that could be affecting physiological state of this species.

Integration of Omics Approaches Enhances the Impact of Scientific Research in Environmental Applications

In the original article [...]

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.

Omics-based ecosurveillance for the assessment of ecosystem function, health, and resilience

Current environmental monitoring efforts often focus on known, regulated contaminants ignoring the potential effects of unmeasured compounds and/or environmental factors. These specific, targeted approaches lack broader environmental information and understanding, hindering effective environmental management and policy. Switching to comprehensive, untargeted monitoring of contaminants, organism health, and environmental factors, such as nutrients, temperature, and pH, would provide more effective monitoring with a likely concomitant increase in environmental health. However, even this method would not capture subtle biochemical changes in organisms induced by chronic toxicant exposure. Ecosurveillance is the systematic collection, analysis, and interpretation of ecosystem health-related data that can address this knowledge gap and provide much-needed additional lines of evidence to environmental monitoring programs. Its use would therefore be of great benefit to environmental management and assessment. Unfortunately, the science of 'ecosurveillance', especially omics-based ecosurveillance is not well known. Here, we give an overview of this emerging area and show how it has been beneficially applied in a range of systems. We anticipate this review to be a starting point for further efforts to improve environmental monitoring via the integration of comprehensive chemical assessments and molecular biology-based approaches. Bringing multiple levels of omics technology-based assessment together into a systems-wide ecosurveillance approach will bring a greater understanding of the environment, particularly the microbial communities upon which we ultimately rely to remediate perturbed ecosystems.

Metabolomics insight into the influence of environmental factors in responses of freshwater biofilms to the model herbicide diuron

Freshwater biofilms have been increasingly used during the last decade in ecotoxicology due to their ecological relevance to assess the effect(s) of environmental stress at the community level. Despite growing knowledge about the effect of various stressors on the structure and the function of these microbial communities, a strong research effort is still required to better understand their response to chemical stress and the influence of environmental stressors in this response. To tackle this challenge, untargeted metabolomics is an approach of choice because of its capacity to give an integrative picture of the exposure to multiple stress and associated effect as well as identifying the molecular pathways involved in these responses. In this context, the present study aimed to explore the use of an untargeted metabolomics approach to unravel at the molecular/biochemical level the response of the whole biofilm to chemical stress and the influence of various environmental factors in this response. To this end, archived high-resolution mass spectrometry data from previous experiments at our laboratory on the effect of the model photosynthesis inhibitor diuron on freshwater biofilm were investigated by using innovative solutions for OMICs data (e.g., DRomics) and more usual chemometric approaches (multivariate and univariate statistical analyses). The results showed a faster (1 min) and more sensitive response of the metabolome to diuron than usual functional descriptors, including photosynthesis. Also, the metabolomics response to diuron resulted from metabolites following various trends (increasing, decreasing, U/bell shape) along increasing concentration and time. This metabolomics response was influenced by the temperature, photoperiod, and flow. A focus on a plant-specific omega-3 (eicosapentaenoic acid) playing a key role in the trophic chain highlighted the potential relevance of metabolomics approach to establish the link between molecular alteration and ecosystem structure/functioning impairment but also how complex is the response and the influence of all the tested factors on this response at the metabolomics level. Altogether, our results underline that more fundamental researches are needed to decipher the metabolomics response of freshwater biofilm to chemical stress and its link with physiological, structural, and functional responses toward the unraveling of adverse outcome pathways (AOP) for key ecosystem functions (e.g., primary production).

Evaluation of transcriptional biomarkers using a high-resolution regression approach: Concentration-dependence of selected transcripts in copper-exposed freshwater mussels (Anodonta anatina)

We tested concentration-dependence of selected gene transcripts (cat, gst, hsp70, hsp90, mt and sod) for evaluation as biomarkers of chemical stress. Contrary to the common approach of factorial designs and few exposure concentrations, we used regression across a high-resolution concentration series. Specifically, freshwater mussels (Anodonta anatina) were acutely (96 h) exposed to Cu (13 nominal concentrations, measuring 0.13-1 600 µg/L), and transcripts were measured by RT-qPCR. In digestive glands, cat, hsp90 and mt decreased with water Cu (p < 0.05), but response magnitudes saturated at < 2-fold decreases. In gills, gst, hsp70, hsp90 and mt increased with water Cu (p < 0.05). While hsp70, hsp90 and mt exceeded 2-fold increases within the exposure range, high Cu concentrations were required (38-160 µg/L). Although gill responses were generally more robust compared to digestive glands, overall small response magnitudes and moderate sensitivity may set limit for potential application as general biomarkers of chemical stress.

Molecular biomarker responses in the freshwater mussel Anodonta anatina exposed to an industrial wastewater effluent

Using a selection of molecular biomarkers, we evaluated responses in freshwater mussels (Anodonta anatina) exposed to effluent from an industrial wastewater treatment facility. The aims of this work were to (1) assess biomarkers of general toxicity under sublethal exposure to an anthropogenic mixture of chemicals, represented by an arbitrary effluent, and (2) evaluate the potential of A. anatina as a bioindicator of pollution. Adult mussels (n = in total 32; 24 males and 8 females) were exposed (96 h) in the laboratory to a fixed dilution of effluent or to a control treatment of standardized freshwater. Metal concentrations were in general higher in the effluent, by an order of magnitude or more, compared to the control. Toxic unit estimates were used as proxies of chemical stress, and Cu, Ni, and Zn were identified as potential major contributors (Cu> Ni > Zn). Six transcriptional (cat, gst, hsp70, hsp90, mt, sod) and two biochemical (AChE, GST) biomarkers were analyzed in two tissues, gills, and digestive glands. Out of the 16 responses (eight biomarkers × two tissues), 14 effect sizes were small (within ± 28 % of control) and differences non-significant (p > 0.05). Results did however show that (1) AChE activity increased by 40% in gills of exposed mussels compared to control, (2) hsp90 expression was 100% higher in exposed female gills compared to control, and (3) three marker signals (AChE in both tissues, and hsp70 in gills) differed between sexes, independent of treatment. Results highlight a need for further investigation of molecular biomarker variability and robustness in A. anatina.

An optimized LC-HRMS untargeted metabolomics workflow for multi-matrices investigations in the three-spined stickleback

Environmental metabolomics has become a growing research field to understand biological and biochemical perturbations of organisms in response to various abiotic or biotic stresses. It focuses on the comprehensive and systematic analysis of a biologic system’s metabolome. This allows the recognition of biochemical pathways impacted by a stressor, and the identification of some metabolites as biomarkers of potential perturbations occurring in a body. In this work, we describe the development and optimization of a complete reliable methodology based on liquid chromatography coupled to high resolution mass spectrometry (LC-HRMS) for untargeted metabolomics studies within a fish model species, the three-spined stickleback (Gasterosteus aculeatus). We evaluated the differences and also the complementarities between four different matrices (brain, gills, liver and whole fish) to obtain metabolome information. To this end, we optimized and compared sample preparation and the analytical method, since the type and number of metabolites detected in any matrix are closely related to these latter. For the sample preparation, a solid-liquid extraction was performed on a low quantity of whole fish, liver, brain, or gills tissues using combinations of methanol/water/heptane. Based on the numbers of features observed in LC-HRMS and on the responses of analytical standards representative of different metabolites groups (amino acids, sugars…), we discuss the influence of the nature, volume, and ratio of extraction solvents, the sample weight, and the reconstitution solvent. Moreover, the analytical conditions (LC columns, pH and additive of mobile phases and ionization modes) were also optimized so as to ensure the maximum metabolome coverages. Thus, two complementary chromatographic procedures were combined in order to cover a broader range of metabolites: a reversed phase separation (RPLC) on a C18 column followed by detection with positive ionization mode (ESI+) and a hydrophilic interaction chromatography (HILIC) on a zwitterionic column followed by detection with negative ionization mode (ESI-). This work provides information on brain, gills, liver, vs the whole body contribution to the stickleback metabolome. These information would help to guide ecotoxicological and biomonitoring studies.

Human impacts on deep-sea sponge grounds: Applying environmental omics to monitoring

Sponges (Phylum Porifera) are the oldest extant Metazoans. In the deep sea, sponges can occur at high densities forming habitats known as sponge grounds. Sponge grounds can extend over large areas of up to hundreds of km² and are biodiversity hotspots. However, as human activities, including deep-water hydrocarbon extraction, continue to expand into areas harbouring sponge grounds, understanding how anthropogenic impacts affect sponges and the ecosystem services they provide at multiple biological scales (community, individual and (sub)cellular levels) is key for achieving sustainable management. This chapter (1) provides an update to the chapter of Advances in Marine Biology Volume 79 entitled "Potential Impacts of Offshore Oil and Gas Activities on Deep-Sea Sponges and the Habitats They Form" and (2) discusses the use of omics as a future tool for deep-sea ecosystem monitoring. While metagenomics and (meta)transcriptomics studies have contributed to improve our understanding of sponge biology in recent years, metabolomics analysis has mostly been used to identify natural products. The sponge metabolome, therefore, remains vastly unknown despite the fact that the metabolome is a key link between the genotype and phenotype, giving us a unique new insight to how key components of an ecosystem are functioning. As the fraction of the metabolome released into the seawater, the sponge exometabolome has only just started to be characterised in comparative environmental metabolomic studies. Yet, the sponge exometabolome constitute a unique opportunity for the identification of biomarkers of sponge health as compounds can be measured in seawater, bypassing the need for physical samples which can still be difficult to collect in the deep sea. Within sponge grounds, the characterisation of a shared sponge exometabolome could lead to the identification of biomarkers of ecosystem functioning and overall health. Challenges remain in establishing omics approaches in environmental monitoring but constant technological advances and reduction in costs means these techniques will become widely available in the future.

Modeling the metabolic profile of Mytilus edulis reveals molecular signatures linked to gonadal development, sex and environmental site

The monitoring of anthropogenic chemicals in the aquatic environment including their potential effects on aquatic organisms, is important for protecting life under water, a key sustainable development goal. In parallel with monitoring the concentrations of chemicals of concern, sentinel species are often used to investigate the biological effects of contaminants. Among these, bivalve molluscs such as mussels are filter-feeding and sessile, hence an excellent model system for measuring localized pollution. This study investigates the relationship between the metabolic state of the blue mussel (Mytilus edulis) and its physiology in different environments. We developed a computational model based on a reference site (relatively unpolluted) and integrated seasonal dynamics of metabolite relative concentrations with key physiological indicators and environmental parameters. The analysis of the model revealed that changes in metabolite levels during an annual cycle are influenced by water temperature and are linked to gonadal development. This work supports the importance of data-driven biology and its potential in environmental monitoring.

Transcriptome-Wide Prediction and Measurement of Combined Effects Induced by Chemical Mixture Exposure in Zebrafish Embryos

BACKGROUND

Humans and environmental organisms are constantly exposed to complex mixtures of chemicals. Extending our knowledge about the combined effects of chemicals is thus essential for assessing the potential consequences of these exposures. In this context, comprehensive molecular readouts as retrieved by omics techniques are advancing our understanding of the diversity of effects upon chemical exposure. This is especially true for effects induced by chemical concentrations that do not instantaneously lead to mortality, as is commonly the case for environmental exposures. However, omics profiles induced by chemical exposures have rarely been systematically considered in mixture contexts.

OBJECTIVES

In this study, we aimed to investigate the predictability of chemical mixture effects on the whole-transcriptome scale.

METHODS

We predicted and measured the toxicogenomic effects of a synthetic mixture on zebrafish embryos. The mixture contained the compounds diuron, diclofenac, and naproxen. To predict concentration- and time-resolved whole-transcriptome responses to the mixture exposure, we adopted the mixture concept of concentration addition. Predictions were based on the transcriptome profiles obtained for the individual mixture components in a previous study. Finally, concentration- and time-resolved mixture exposures and subsequent toxicogenomic measurements were performed and the results were compared with the predictions.

RESULTS

This comparison of the predictions with the observations showed that the concept of concentration addition provided reasonable estimates for the effects induced by the mixture exposure on the whole transcriptome. Although nonadditive effects were observed only occasionally, combined, that is, multicomponent-driven, effects were found for mixture components with anticipated similar, as well as dissimilar, modes of action.

DISCUSSION

Overall, this study demonstrates that using a concentration- and time-resolved approach, the occurrence and size of combined effects of chemicals may be predicted at the whole-transcriptome scale. This allows improving effect assessment of mixture exposures on the molecular scale that might not only be of relevance in terms of risk assessment but also for pharmacological applications. https://doi.org/10.1289/EHP7773.

Immunological Responses of Marine Bivalves to Contaminant Exposure: Contribution of the -Omics Approach

The increasing number of data studies on the biological impact of anthropogenic chemicals in the marine environment, together with the great development of invertebrate immunology, has identified marine bivalves as a key invertebrate group for studies on immunological responses to pollutant exposure. Available data on the effects of contaminants on bivalve immunity, evaluated with different functional and molecular endpoints, underline that individual functional parameters (cellular or humoral) and the expression of selected immune-related genes can distinctly react to different chemicals depending on the conditions of exposure. Therefore, the measurement of a suite of immune biomarkers in hemocytes and hemolymph is needed for the correct evaluation of the overall impact of contaminant exposure on the organism's immunocompetence. Recent advances in -omics technologies are revealing the complexity of the molecular players in the immune response of different bivalve species. Although different -omics represent extremely powerful tools in understanding the impact of pollutants on a key physiological function such as immune defense, the -omics approach has only been utilized in this area of investigation in the last few years. In this work, available information obtained from the application of -omics to evaluate the effects of pollutants on bivalve immunity is summarized. The data shows that the overall knowledge on this subject is still quite limited and that to understand the environmental relevance of any change in immune homeostasis induced by exposure to contaminants, a combination of both functional assays and cutting-edge technology (transcriptomics, proteomics, and metabolomics) is required. In addition, the utilization of metagenomics may explain how the complex interplay between the immune system of bivalves and its associated bacterial communities can be modulated by pollutants, and how this may in turn affect homeostatic processes of the host, host–pathogen interactions, and the increased susceptibility to disease. Integrating different approaches will contribute to knowledge on the mechanism responsible for immune dysfunction induced by pollutants in ecologically and economically relevant bivalve species and further explain their sensitivity to multiple stressors, thus resulting in health or disease.

20th Pollutant Responses in Marine Organisms (PRIMO 20): Global issues and fundamental mechanisms caused by pollutant stress in marine and freshwater organisms

The 20th Pollutant Responses in Marine Organisms (PRIMO 20) conference provided a forum for scientists from around the world to communicate novel toxicological research findings specifically focused on aquatic organisms, by combining applied and basic research at the intersection of environmental and mechanistic toxicology. The work highlighted in this special issue of Aquatic Toxicology, a special issue of Marine Environmental Research, and presented through posters and presentations, encompass important and emerging topics in freshwater and marine toxicology. This includes multiple types of emerging contaminants including microplastics and UV filtering chemicals. Other studies aimed to further our understanding of the effects of endocrine disrupting chemicals, pharmaceuticals, and personal care products. Further research presented in this virtual issue examined the interactive effects of chemicals and pathogens, while the final set of manuscripts demonstrates continuing efforts to combine traditional biomonitoring, data from -omic technologies, and modeling for use in risk assessment and management. An additional goal of PRIMO meetings is to address the link between environmental and human health. Several articles in this issue of Aquatic Toxicology describe the appropriateness of using aquatic organisms as models for human health, while the keynote speakers, as described in the editorial below, presented research that highlighted bioaccumulation of contaminants such as PFOS and mercury from fish to marine mammals and coastal human populations such as the Gullah/GeeChee near Charleston, South Carolina, USA.

Drivers of and Obstacles to the Adoption of Toxicogenomics for Chemical Risk Assessment: Insights from Social Science Perspectives

BACKGROUND

Some 20 y ago, scientific and regulatory communities identified the potential of omics sciences (genomics, transcriptomics, proteomics, metabolomics) to improve chemical risk assessment through development of toxicogenomics. Recognizing that regulators adopt new scientific methods cautiously given accountability to diverse stakeholders, the scope and pace of adoption of toxicogenomics tools and data have nonetheless not met the ambitious, early expectations of omics proponents.

OBJECTIVE

Our objective was, therefore, to inventory, investigate, and derive insights into drivers of and obstacles to adoption of toxicogenomics in chemical risk assessment. By invoking established social science frameworks conceptualizing innovation adoption, we also aimed to develop recommendations for proponents of toxicogenomics and other new approach methodologies (NAMs).

METHODS

We report findings from an analysis of 56 scientific and regulatory publications from 1998 through 2017 that address the adoption of toxicogenomics for chemical risk assessment. From this purposeful sample of toxicogenomics discourse, we identified major categories of drivers of and obstacles to adoption of toxicogenomics tools and data sets. We then mapped these categories onto social science frameworks for conceptualizing innovation adoption to generate actionable insights for proponents of toxicogenomics.

DISCUSSION

We identify the most salient drivers and obstacles. From 1998 through 2017, adoption of toxicogenomics was understood to be helped by drivers such as those we labeled Superior scientific understanding, New applications, and Reduced cost & increased efficiency but hindered by obstacles such as those we labeled Insufficient validation, Complexity of interpretation, and Lack of standardization. Leveraging social science frameworks, we find that arguments for adoption that draw on the most salient drivers, which emphasize superior and novel functionality of omics as rationales, overlook potential adopters' key concerns: simplicity of use and compatibility with existing practices. We also identify two perspectives-innovation-centric and adopter-centric-on omics adoption and explain how overreliance on the former may be undermining efforts to promote toxicogenomics. https://doi.org/10.1289/EHP6500.

Assessing legacy and endocrine disrupting pollutants in Boston Harbor with transcriptomic biomarkers

Within monitoring frameworks, biomarkers provide several benefits because they serve as intermediates between pollutant exposure and effects, and integrate the responses of contaminants that operate through the same mechanism of action. This study was designed to verify the use of transcriptomic biomarkers developed in our prior work (i.e., Coastal Biosensor of Endocrine Disruption; C-BED assay) on Mytilus edulis and identify additional biomarkers for legacy pollutants. M. edulis were collected from a reference site in Pemaquid, ME, USA and deployed by the Massachusetts Water Resources Authority (MWRA) at locations in and outside Boston Harbor, MA, USA: including (1) Boston Inner Harbor (IH), (2) the current outfall (OS), (3) 1 km away from the current outfall (LNB), and (4) Deer Island (DI), the site where untreated wastewater was formerly discharged into the bay. Differential gene expression was quantified with a high density microarray. Seven genes significantly correlated with whole tissue concentration of PAHs, and six genes significantly correlated with whole body concentrations of PCBs, two groups of legacy contaminants that were elevated at stations IH, OS, and DI. Enrichment analysis indicated that IH mussels had the highest induction of stress response genes, which correlated with the higher levels of contaminants measured at this site. Based on the C-BED assay gene analysis, stations IH and OS exhibited signs of endocrine disruption, which were further confirmed by incorporating the results for the C-BED assay within the Integrated Biomarker Response (IBR) approach. This study successfully demonstrated the potential use of transcriptomic biomarkers within a monitoring program to identify the presence and organismal responses to endocrine disrupting and legacy contaminant classes.

Transcriptional and biochemical biomarker responses in a freshwater mussel (Anodonta anatina) under environmentally relevant Cu exposure

Molecular biomarkers, like gene transcripts or enzyme activities, are potentially powerful tools for early warning assessment of pollution. However, a thorough understanding of response and baseline variation is required to distinguish actual effects from pollution. Here, we assess the freshwater mussel Anodonta anatina as a biomarker model species for freshwater ecosystems, by testing responses of six transcriptional (cat, gst, hsp70, hsp90, mt, and sod) and two biochemical (AChE and GST) biomarkers to environmentally relevant Cu water concentrations. Mussels (n = 20), collected from a stream free from point source pollution, were exposed in the laboratory, for 96 h, to Cu treatments (< 0.2 μg/L, 0.77 ± 0.87 μg/L, and 6.3 ± 5.4 μg/L). Gills and digestive glands were extracted and analyzed for transcriptional and biochemical responses. Biological and statistical effect sizes from Cu treatments were in general small (mean log₂ fold-change ≤ 0.80 and Cohen's f ≤ 0.69, respectively), and no significant treatment effects were observed. In contrast, four out of eight biomarkers (cat, gst, hsp70, and GST) showed a significant sex:tissue interaction, and additionally one (sod) showed significant overall effects from sex. Specifically, three markers in gills (cat, mt, GST) and one in digestive gland (AChE) displayed significant sex differences, independent of treatment. Results suggest that sex or tissue effects might obscure low-magnitude biomarker responses and potential early warnings. Thus, variation in biomarker baselines and response patterns needs to be further addressed for the future use of A. anatina as a biomarker model species.

Reproduction Stage Differentiates the Time-Course Regulation of Metabolites in Daphnia magna

Daphnia magna is a keystone indicator zooplankton used in environmental quality assessments. Comparative metabolomics, which contrasts small biomolecular regulations under different conditions, has emerged as a sensitive and informative investigation tool for the assessment of environmental stressors on D. magna. Baseline metabolomic variation is likely impacted by the asexual reproduction cycle of D. magna; however, the relationship between metabolite concentration regulation and reproduction cycle has not been investigated. This study investigated the time-course regulation of 51 metabolites during different reproduction stages to determine how the reproduction cycle controlled the metabolite profile of D. magna. Two-way analysis of variance (ANOVA) results reveal that most metabolites show significantly differentiated concentrations by individual or a combination of reproduction stages and sampling time. On the basis of the partial least-squares-discriminant analysis (PLS-DA) and ANOVA-simultaneous component analysis (ANOVA-SCA), stages 2 and 3 of reproduction show similarity in metabolite abundance regulation compared to stage 1. Metabolites were classified as being either dependent or independent of reproduction in the pairwise concentration correlation based on Pearson correlation coefficients. This study observed that the D. magna reproduction stage is an important consideration and potential variable and should be considered carefully when conducting metabolomic experiments using D. magna.

Metabolic profiling of Daphnia magna exposure to a mixture of hydrophobic organic contaminants in the presence of dissolved organic matter

The hydrophobic organic contaminants triclosan, triphenyl phosphate (TPhP) and diazinon sorb to dissolved organic matter (DOM) and this may alter their bioavailability and toxicity. ¹H nuclear magnetic resonance (NMR)-based metabolomics was used to investigate how DOM at 1 and 5 mg organic carbon/L may alter the metabolome of Daphnia magna from exposure to equitoxic mixtures of triclosan, TPhP and diazinon. These contaminants have different modes of action toward D. magna. The contaminant concentrations in each mixture were an equal percentage of their lethal concentration to 50% of the population (LC₅₀) values, which equates to 1250 μg/L TPhP, 330 μg/L triclosan and 0.9 μg/L diazinon. The ternary mixture exposure at 1% LC₅₀ values did not alter the D. magna metabolome. Contaminant mixture exposures at 5%, 10%, and 15% LC₅₀ values decreased glucose, serine and glycine concentrations and increased asparagine and threonine concentrations, suggesting disruptions in energy metabolism. The contaminant mixture had a unique mode of action in D. magna and DOM at 1 and 5 mg organic carbon/L did not change this mode of action. The estimated sorption of triclosan, TPhP or diazinon to DOM at 1 or 5 mg organic carbon/L in this experimental design was calculated to be <50% for each contaminant. This suggests that the mode of action of the contaminant mixture was not altered by DOM because the two environmentally relevant concentrations of DOM may have not substantially altered contaminant bioavailability. Our results indicate that DOM may not inevitably mitigate or alter the sub-lethal toxicity of a mixture of hydrophobic organic contaminants. This indicates the complexity of predicting the molecular-level toxicity of environmental mixtures. For adequate risk assessment of freshwater ecosystems, it is vital to account for the combined sub-lethal toxicity of an environmental mixture of contaminants.

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.

Map and model-moving from observation to prediction in toxicogenomics

BACKGROUND

Chemicals induce compound-specific changes in the transcriptome of an organism (toxicogenomic fingerprints). This provides potential insights about the cellular or physiological responses to chemical exposure and adverse effects, which is needed in assessment of chemical-related hazards or environmental health. In this regard, comparison or connection of different experiments becomes important when interpreting toxicogenomic experiments. Owing to lack of capturing response dynamics, comparability is often limited. In this study, we aim to overcome these constraints.

RESULTS

We developed an experimental design and bioinformatic analysis strategy to infer time- and concentration-resolved toxicogenomic fingerprints. We projected the fingerprints to a universal coordinate system (toxicogenomic universe) based on a self-organizing map of toxicogenomic data retrieved from public databases. Genes clustering together in regions of the map indicate functional relation due to co-expression under chemical exposure. To allow for quantitative description and extrapolation of the gene expression responses we developed a time- and concentration-dependent regression model. We applied the analysis strategy in a microarray case study exposing zebrafish embryos to 3 selected model compounds including 2 cyclooxygenase inhibitors. After identification of key responses in the transcriptome we could compare and characterize their association to developmental, toxicokinetic, and toxicodynamic processes using the parameter estimates for affected gene clusters. Furthermore, we discuss an association of toxicogenomic effects with measured internal concentrations.

CONCLUSIONS

The design and analysis pipeline described here could serve as a blueprint for creating comparable toxicogenomic fingerprints of chemicals. It integrates, aggregates, and models time- and concentration-resolved toxicogenomic data.

The use of gene expression to unravel the single and mixture toxicity of abamectin and difenoconazole on survival and reproduction of the springtail Folsomia candida

Pesticides risk assessments have traditionally focused on the effects on standard parameters, such as mortality, reproduction and development. However, one of the first signs of adverse effects that occur in organisms exposed to stress conditions is an alteration in their genomic expression, which is specific to the type of stress, sensitive to very low contaminant concentrations and responsive in a few hours. The aim of the present study was to evaluate the single and binary mixture toxicity of commercial products of abamectin (Kraft^® 36 EC) and difenoconazole (Score^® 250 EC) to Folsomia candida. Laboratory toxicity tests were conducted to access the effects of these pesticides on springtail survival, reproduction and gene expression. The reproduction assays gave EC₅₀ and EC₁₀ values, respectively, of 6.3 and 1.4 mg a.s./kg dry soil for abamectin; 1.0 and 0.12 mg a.s./kg dry soil for Kraft^® 36 EC; and 54 and 23 mg a.s./kg dry soil for Score^® 250 EC. Technical difenoconazole did not have any effect at the concentrations tested. No significant differences in gene expression were found between the abamectin concentrations tested (EC₁₀ and EC₅₀) and the solvent control. Exposure to Kraft^® 36 EC, however, significantly induced Cyp6 expression at the EC₅₀ level, while VgR was significantly downregulated at both the EC₁₀ and EC₅₀. Exposure to the simple pesticide mixture of Kraft^® 36 EC + Score^® 250 EC caused significant up regulation of ABC transporter, and significant down regulation of VgR relative to the controls. GABA receptor also showed significant down-regulation between the EC₁₀ and EC₅₀ mixture treatments. Results of the present study demonstrate that pesticide-induced gene expression effects precede and occur at lower concentrations than organism-level responses. Integrating "omic" endpoints in traditional bioassays may thus be a promising way forward in pesticide toxicity evaluations.

Profiling physicochemical and planktonic features from discretely/continuously sampled surface water

There is an increasing need for assessing aquatic ecosystems that are globally endangered. Since aquatic ecosystems are complex, integrated consideration of multiple factors utilizing omics technologies can help us better understand aquatic ecosystems. An integrated strategy linking three analytical (machine learning, factor mapping, and forecast-error-variance decomposition) approaches for extracting the features of surface water from datasets comprising ions, metabolites, and microorganisms is proposed herein. The three developed approaches can be employed for diverse datasets of sample sizes and experimentally analyzed factors. The three approaches are applied to explore the features of bay water surrounding Odaiba, Tokyo, Japan, as a case study. Firstly, the machine learning approach separated 681 surface water samples within Japan into three clusters, categorizing Odaiba water into seawater with relatively low inorganic ions, including Mg, Ba, and B. Secondly, the factor mapping approach illustrated Odaiba water samples from the summer as rich in multiple amino acids and some other metabolites and poor in inorganic ions relative to other seasons based on their seasonal dynamics. Finally, forecast-error-variance decomposition using vector autoregressive models indicated that a type of microalgae (Raphidophyceae) grows in close correlation with alanine, succinic acid, and valine on filters and with isobutyric acid and 4-hydroxybenzoic acid in filtrate, Ba, and average wind speed. Our integrated strategy can be used to examine many biological, chemical, and environmental physical factors to analyze surface water.

Reprint of: Environmental toxicology and omics: A question of sex

Molecular initiating events and downstream transcriptional/proteomic responses provide valuable information for adverse outcome pathways, which can be used predict the effects of chemicals on physiological systems. There has been a paucity of research that addresses sex-specific expression profiling in toxicology and due to cost, time, and logistic considerations, sex as a variable has not been widely considered. In response to this deficiency, federal agencies in the United States, Canada, and Europe have highlighted the importance of including sex as a variable in scientific investigations. Using case studies from both aquatic and mammalian toxicology, we report that there can be less than ~20-25% consensus in how the transcriptome and proteome of each sex responds to chemicals. Chemicals that have been shown to elicit sex-specific responses in the transcriptome or proteome include pharmaceuticals, anti-fouling agents, anticorrosive agents, and fungicides, among others. Sex-specific responses in the transcriptome and proteome are not isolated to whole animals, as investigations demonstrate that primary cell cultures isolated from each sex responds differently to toxicants. This signifies that sex is important, even in cell lines. Sex has significant implications for predictive toxicology, and both male and female data are required to improve robustness of adverse outcome pathways.

BIOLOGICAL SIGNIFICANCE

Clinical toxicology recognizes that sex is an important variable, as pharmacokinetics (ADME; absorption, distribution, metabolism, and excretion) can differ between females and males. However, few studies in toxicology have explored the implication of sex in relation to the transcriptome and proteome of whole organisms. High-throughput molecular approaches are becoming more frequently applied in toxicity screens (e.g. pre-clinical experiments, fish embryos, cell lines, synthetic tissues) and such data are expected to build upon reporter-based cell assays (e.g. receptor activation, enzyme inhibition) used in toxicant screening programs (i.e. Tox21, ToxCast, REACH). Thus, computational models can more accurately predict the diversity of adverse effects that can occur from chemical exposure within the biological system. Our studies and those synthesized from the literature suggest that the transcriptome and proteome of females and males respond quite differentially to chemicals. This has significant implications for predicting adverse effects in one sex when using molecular data generated in the other sex. While molecular initiating events are not expected to differ dramatically between females and males (i.e. an estrogen binds estrogen receptors in both sexes), it is important to acknowledge that the downstream transcriptomic and proteomic responses can differ based upon the presence/absence of co-regulators and inherent sex-specific variability in regulation of transcriptional and translational machinery. Transcriptomic and proteomic studies also reveal that cell processes affected by chemicals can differ due to sex, and this can undoubtedly lead to sex-specific physiological responses.

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.

Zinc Oxide Nanoparticle Caused Plasma Metabolomic Perturbations Correlate with Hepatic Steatosis

Zinc oxide nanoparticles (ZnO NPs), known for their chemical stability and strong adsorption, are used in everyday items such as cosmetics, sunscreens, and prophylactic drugs. However, they have also been found to adversely affect organisms; previously we found that ZnO NPs disrupt pubertal ovarian development, inhibit embryonic development by upsetting γ-H2AX and NF-κB pathways, and even disturb skin stem cells. Non-targeted metabolomic analysis of biological organisms has been suggested as an unbiased tool for the investigation of perturbations in response to NPs and their underlying mechanisms. Although metabolomics has been used in nanotoxicological studies, very few reports have used it to investigate the effects of ZnO NPs exposure. In the current investigation, through a metabolomics-based approach, we discovered that ZnO NPs caused changes in plasma metabolites involved in anti-oxidative mechanisms, energy metabolism, and lipid metabolism in hen livers. These results are in line with earlier findings that ZnO NPs perturb the tricarboxylic acid cycle and in turn result in the use of alternative energy sources. We also found that ZnO NPs disturbed lipid metabolism in the liver and consequently impacted blood lipid balance. Changes in plasma metabolomes were correlated with hepatic steatosis.

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.

How consistent are we? Interlaboratory comparison study in fathead minnows using the model estrogen 17α-ethinylestradiol to develop recommendations for environmental transcriptomics

Fundamental questions remain about the application of omics in environmental risk assessments, such as the consistency of data across laboratories. The objective of the present study was to determine the congruence of transcript data across 6 independent laboratories. Male fathead minnows were exposed to a measured concentration of 15.8 ng/L 17α-ethinylestradiol (EE2) for 96 h. Livers were divided equally and sent to the participating laboratories for transcriptomic analysis using the same fathead minnow microarray. Each laboratory was free to apply bioinformatics pipelines of its choice. There were 12 491 transcripts that were identified by one or more of the laboratories as responsive to EE2. Of these, 587 transcripts (4.7%) were detected by all laboratories. Mean overlap for differentially expressed genes among laboratories was approximately 50%, which improved to approximately 59.0% using a standardized analysis pipeline. The dynamic range of fold change estimates was variable between laboratories, but ranking transcripts by their relative fold difference resulted in a positive relationship for comparisons between any 2 laboratories (mean R²  > 0.9, p < 0.001). Ten estrogen-responsive genes encompassing a fold change range from dramatic (>20-fold; e.g., vitellogenin) to subtle (∼2-fold; i.e., block of proliferation 1) were identified as differentially expressed, suggesting that laboratories can consistently identify transcripts that are known a priori to be perturbed by a chemical stressor. Thus, attention should turn toward identifying core transcriptional networks using focused arrays for specific chemicals. In addition, agreed-on bioinformatics pipelines and the ranking of genes based on fold change (as opposed to p value) should be considered in environmental risk assessment. These recommendations are expected to improve comparisons across laboratories and advance the use of omics in regulations. Environ Toxicol Chem 2017;36:2593-2601. © 2017 SETAC.

Integrated assessment of toxic effects of maghemite (γ-Fe2O3) nanoparticles in zebrafish

The increasing use of nanotechnology in the last decade has raised concerns about the impact of nanoparticles in the environment. In particular, the potential harmful effects of iron oxide nanoparticles (IONPs) in aquatic organisms have been poorly addressed. We analyze here the toxic effects induced by IONPs in zebrafish using a combination of classical (genotoxicity, oxidative stress) and molecular (transcriptomic) methodologies. Adult animals were exposed for 96h to five sub-lethal IONP concentrations, ranging from 4.7 to 74.4mg/L. Comet and micronucleus assays revealed a significant number of DNA lesions induced by IONPs at all concentrations tested. Conversely, the thiobarbituric acid reactive substances (TBARS) test detected only a mild oxidative damage in liver cells (∼1.5-fold increase of malondialdehyde concentrations) and only at the two higher IONP concentrations tested. Microarray analysis of liver samples identified 953 transcripts (927 unique genes) differentially expressed between controls and IONP-exposed samples. Subsequent functional analysis identified genes related to cation/metal ion binding, membrane formation, and morphogenesis among the transcripts overrepresented upon IONP treatments, whereas mRNAs encompassing genes associated with RNA biogenesis, translation, ribosomes, and several metabolic processes became underrepresented in treated samples. Taken together, these results indicate considerable genotoxic effects of IONPs combined with general negative effect on cell growth and on the ability of the cell produce new proteins. On the contrary, IONPs showed only a limited capacity to induce oxidative stress. To our knowledge, this is the first study on IONPs toxicity using such an integrative approach in an aquatic organism.

Parental diuron-exposure alters offspring transcriptome and fitness in Pacific oyster Crassostrea gigas

One of the primary challenges in ecotoxicology is to contribute to the assessment of the ecological status of ecosystems. In this study, we used Pacific oyster Crassostrea gigas to explore the effects of a parental exposure to diuron, a herbicide frequently detected in marine coastal environments. The present toxicogenomic study provides evidence that exposure of oyster genitors to diuron during gametogenesis results in changes in offspring, namely, transcriptomic profile alterations, increased global DNA methylation levels and reduced growth and survival within the first year of life. Importantly, we highlighted the limitations to identify particular genes or gene expression signatures that could serve as biomarkers for parental herbicide-exposure and further for multigenerational and transgenerational effects of specific chemical stressors. By analyzing samples from two independent experiments, we demonstrated that, due to complex confounding effects with both tested solvent vehicles, diuron non-specifically affected the offspring transcriptome. These original results question the potential development of predictive genomic tools for detecting specific indirect impacts of contaminants in environmental risk assessments. However, our results indicate that chronic environmental exposure to diuron over several generations may have significant long term impacts on oyster populations with adverse health outcomes.