Towards a system level understanding of non-model organisms sampled from the environment: a network biology approach

Environmentally realistic concentrations of chlorinated, brominated, and fluorinated persistent organic pollutants induce the unfolded protein response as a shared stress pathway in the liver of Atlantic cod (Gadus morhua)

In the North Sea and North Atlantic coastal areas, fish experience relatively high background levels of persistent organic pollutants. This study aimed to compare the mode of action of environmentally relevant concentrations of mixtures of halogenated compounds in Atlantic cod. Juvenile male cod with mean weight of 840 g were exposed by gavage to dietary mixtures of chlorinated (PCBs, DDT analogs, chlordane, lindane, and toxaphene), brominated (PBDEs), and fluorinated (PFOS) compounds for 4 weeks. One group received a combined mixture of all three compound groups. The results showed that the accumulated levels of chemicals in cod liver after 4 weeks of exposure reflected concentrations found in wild fish in this region. Pathway analysis revealed that the treatment effects by each of the three groups of chemicals (chlorinated, brominated, and fluorinated) converged on activation of the unfolded protein response (UPR). Upstream regulator analysis predicted that almost all the key transcription factors (XBP1, ERN1, ATF4, EIF2AK3, and NFE2L2) regulating the UPR were significantly activated. No additive effect was observed in cod co-treated with all three compound groups. In conclusion, the genome-wide transcriptomic study suggests that the UPR pathway is a sensitive common target of halogenated organic environmental pollutants in fish.

Correlations between oxidative DNA damage and formation of hepatic tumours in two flatfish species from contaminated environments

Many species in aquatic environments face increased exposure to oncogenic pollution due to anthropogenic environmental change which can lead to higher cancer prevalence. The mechanistic relationship connecting environmental pollution and cancer is multi-factorial and poorly understood, and the specific mechanisms are so far still uncharacterized. One potential mediator between pollutant exposure and cancer is oxidative damage to DNA. We conducted a study in the field with two flatfish species, European flounder (Platichthys flesus L.) and common dab (Limanda limanda L.) with overlapping distribution and similar ecological niche, to investigate if the link between oncogenic pollutants and cancer described in ecotoxicological literature could be mediated by oxidative DNA damage. This was not the case for flounders as neither polycyclic aromatic hydrocarbon (PAH) bile metabolites nor metallic trace element concentrations were related to oxidative DNA damage measurements. However, dabs with higher PAH concentrations did exhibit increased oxidative damage. High oxidative DNA damage also did not predict neoplasm occurrence, rather, healthy individuals tended to have higher oxidative damage measurements compared to fishes with pre-neoplastic tumours. Our analyses showed that flounders had lower concentrations of PAH bile metabolites, suggesting that compared to dab this species is less exposed or better at eliminating these contaminants.

Three perspectives on the prediction of chemical effects in ecosystems

The increasing production, use and emission of synthetic chemicals into the environment represents a major driver of global change. The large number of synthetic chemicals, limited knowledge on exposure patterns and effects in organisms and their interaction with other global change drivers hamper the prediction of effects in ecosystems. However, recent advances in biomolecular and computational methods are promising to improve our capacity for prediction. We delineate three idealised perspectives for the prediction of chemical effects: the suborganismal, organismal and ecological perspective, which are currently largely separated. Each of the outlined perspectives includes essential and complementary theories and tools for prediction but captures only part of the phenomenon of chemical effects. Links between the perspectives may foster predictive modelling of chemical effects in ecosystems and extrapolation between species. A major challenge for the linkage is the lack of data sets simultaneously covering different levels of biological organisation (here referred to as biological levels) as well as varying temporal and spatial scales. Synthesising the three perspectives, some central aspects and associated types of data seem particularly necessary to improve prediction. First, suborganism- and organism-level responses to chemicals need to be recorded and tested for relationships with chemical groups and organism traits. Second, metrics that are measurable at many biological levels, such as energy, need to be scrutinised for their potential to integrate across levels. Third, experimental data on the simultaneous response over multiple biological levels and spatiotemporal scales are required. These could be collected in nested and interconnected micro- and mesocosm experiments. Lastly, prioritisation of processes involved in the prediction framework needs to find a balance between simplification and capturing the essential complexity of a system. For example, in some cases, eco-evolutionary dynamics and interactions may need stronger consideration. Prediction needs to move from a static to a real-world eco-evolutionary view.

Impacts of chemical stress, season, and climate change on the flounder population of the highly anthropised Seine estuary (France)

The main objective of this study was to improve our knowledge on the responses of fish populations to multistress (diffuse pollution and warming waters) in estuaries. Adult flounders were caught in two estuaries in the Eastern English Channel: the heavily polluted Seine estuary vs the moderately contaminated Canche estuary. Fish samplings were conducted in January just before the reproduction period, and in July when gonads were at rest. The overall rise in coastal winter water temperatures detected over the Channel impairs the flounder's phenology of reproduction in the two estuaries, inducing a delay of maturation process and probably also spawning. The higher liver histopathology index in Seine vs Canche could be the consequence of the fish exposition to a complex cocktail of contaminants in a strongly industrialized estuary. Higher levels of neurotoxicity, gill lipid peroxidation, and liver EROD activity were observed in Seine vs Canche. Furthermore, a possible impairment in mitochondrial metabolism was suggested in the Seine flounder population. We confirmed in this study the potential role of two membrane lipids (sphingomyelin and phosphatidylserine) in the resistance towards oxidative stress in Seine and Canche. Finally, we suggest that the Seine flounder population (and possibly the connected Eastern English Channel flounder populations over the French Coast) could be seriously impacted in the future by multistress: higher winter temperatures and chemical contamination.

Integrative proteomics and metabolomics approach to elucidate metabolic dysfunction induced by silica nanoparticles in hepatocytes

Silica nanoparticles (SiNPs) are derived from manufactured materials and the natural environment, and they cause detrimental effects on human health via various exposure routes. The liver is proven to be a key target organ for SiNP toxicity; however, the mechanisms causing toxicity remain largely uncertain. Here, we investigated the effects of SiNPs on the metabolic spectrum in hepatocytes via integrative analyses of proteomics and metabolomics. First, a proteomic analysis was used to screen for critical proteins (including RPL3, HSP90AA1, SOD, PGK1, GOT1, and PNP), indicating that abnormal protein synthesis, protein misfolding, oxidative stress, and metabolic dysfunction may contribute to SiNP-induced hepatotoxicity. Next, metabolomic data demonstrated that SiNPs caused metabolic dysfunction by altering vital metabolites (including glucose, alanine, GSH, CTP, and ATP). Finally, a systematic bioinformatic analysis of protein-metabolite interactions showed that SiNPs disturbed glucose metabolism (glycolysis and pentose phosphate pathways, amino acid metabolism (alanine, aspartate, and glutamate), and ribonucleotide metabolism (purine and pyrimidine). These metabolic dysfunctions could exacerbate oxidative stress and lead to liver injury. Moreover, SOD, TKT, PGM1, GOT1, PNP, and NME2 may be key proteins for SiNP-induced hepatotoxicity. This study revealed the metabolic mechanisms underlying SiNP-induced hepatotoxicity and illustrated that integrative omics analyses can be a powerful approach for toxicity evaluations and risk assessments of nanoparticles.

The Eco-Exposome Concept: Supporting an Integrated Assessment of Mixtures of Environmental Chemicals

Organisms are exposed to ever-changing complex mixtures of chemicals over the course of their lifetime. The need to more comprehensively describe this exposure and relate it to adverse health effects has led to formulation of the exposome concept in human toxicology. Whether this concept has utility in the context of environmental hazard and risk assessment has not been discussed in detail. In this Critical Perspective, we propose-by analogy to the human exposome-to define the eco-exposome as the totality of the internal exposure (anthropogenic and natural chemicals, their biotransformation products or adducts, and endogenous signaling molecules that may be sensitive to an anthropogenic chemical exposure) over the lifetime of an ecologically relevant organism. We describe how targeted and nontargeted chemical analyses and bioassays can be employed to characterize this exposure and discuss how the adverse outcome pathway concept could be used to link this exposure to adverse effects. Available methods, their limitations, and/or requirement for improvements for practical application of the eco-exposome concept are discussed. Even though analysis of the eco-exposome can be resource-intensive and challenging, new approaches and technologies make this assessment increasingly feasible. Furthermore, an improved understanding of mechanistic relationships between external chemical exposure(s), internal chemical exposure(s), and biological effects could result in the development of proxies, that is, relatively simple chemical and biological measurements that could be used to complement internal exposure assessment or infer the internal exposure when it is difficult to measure. Environ Toxicol Chem 2022;41:30-45. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Protein Abundance Prediction Through Machine Learning Methods

Proteins are responsible for most physiological processes, and their abundance provides crucial information for systems biology research. However, absolute protein quantification, as determined by mass spectrometry, still has limitations in capturing the protein pool. Protein abundance is impacted by translation kinetics, which rely on features of codons. In this study, we evaluated the effect of codon usage bias of genes on protein abundance. Notably, we observed differences regarding codon usage patterns between genes coding for highly abundant proteins and genes coding for less abundant proteins. Analysis of synonymous codon usage and evolutionary selection showed a clear split between the two groups. Our machine learning models predicted protein abundances from codon usage metrics with remarkable accuracy, achieving strong correlation with experimental data. Upon integration of the predicted protein abundance in enzyme-constrained genome-scale metabolic models, the simulated phenotypes closely matched experimental data, which demonstrates that our predictive models are valuable tools for systems metabolic engineering approaches.

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

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

Immunotoxicity of Xenobiotics in Fish: A Role for the Aryl Hydrocarbon Receptor (AhR)?

The impact of anthropogenic contaminants on the immune system of fishes is an issue of growing concern. An important xenobiotic receptor that mediates effects of chemicals, such as halogenated aromatic hydrocarbons (HAHs) and polyaromatic hydrocarbons (PAHs), is the aryl hydrocarbon receptor (AhR). Fish toxicological research has focused on the role of this receptor in xenobiotic biotransformation as well as in causing developmental, cardiac, and reproductive toxicity. However, biomedical research has unraveled an important physiological role of the AhR in the immune system, what suggests that this receptor could be involved in immunotoxic effects of environmental contaminants. The aims of the present review are to critically discuss the available knowledge on (i) the expression and possible function of the AhR in the immune systems of teleost fishes; and (ii) the impact of AhR-activating xenobiotics on the immune systems of fish at the levels of immune gene expression, immune cell proliferation and immune cell function, immune pathology, and resistance to infectious disease. The existing information indicates that the AhR is expressed in the fish immune system, but currently, we have little understanding of its physiological role. Exposure to AhR-activating contaminants results in the modulation of numerous immune structural and functional parameters of fish. Despite the diversity of fish species studied and the experimental conditions investigated, the published findings rather uniformly point to immunosuppressive actions of xenobiotic AhR ligands in fish. These effects are often associated with increased disease susceptibility. The fact that fish populations from HAH- and PAH-contaminated environments suffer immune disturbances and elevated disease susceptibility highlights that the immunotoxic effects of AhR-activating xenobiotics bear environmental relevance.

Transcriptomic, proteomic and metabolomic profiling unravel the mechanisms of hepatotoxicity pathway induced by triphenyl phosphate (TPP)

Triphenyl phosphate (TPP) has been found in various environmental media and in biota suggesting widespread human exposure. However, there is still insufficient information on the hepatotoxicity mechanisms of health risk exposed to TPP. In this study, TPP could induce human normal liver cell (L02) apoptosis, injury cell ultrastructure and elevate the levels of reactive oxygen species (ROS). The integrated multi-omic (transcriptomic, proteomic, and metabolomic) analysis was used to further investigate the mechanisms. Transcriptomic analysis revealed that TPP exposure markedly affected cell apoptosis, oncogene activation, REDOX homeostasis, DNA damage and repair. Additionally, proteomic analysis found that the related proteins associated with apoptosis, oxidative stress, metabolism and membrane structure were affected. And metabolomic analysis verified that the related metabolic pathways, such as glycolysis, citrate cycle, oxidative phosphorylation, lipid and protein metabolism, were also significantly disrupted. Based on the multi-omic results, a hypothesized network was constructed to discover the key molecular events in response to TPP and illustrate the mechanism of TPP-induced hepatotoxicity in L02 cells. Therefore, molecular responses could be elucidated at multiple biological levels, and multi-omic analysis could provide scientific tools for exploring potential mechanisms of toxicity and chemical risk assessment.

An Integrated Biomarker Approach Using Flounder to Improve Chemical Risk Assessments in the Heavily Polluted Seine Estuary

The objective of this study was to develop an integrative approach in ecotoxicology (from biomarkers to population genetics) to assess the ecological status of fish populations. Flounders (Platichthys flesus) collected after the spawning season in the heavily polluted Seine estuary were compared with the moderately polluted Bay of Douarnenez. The muscle energetic reserves were highly depleted in Seine vs. Douarnenez fish. The Seine fish displaying a reduced capacity to manage the oxidative stress and a higher energetic metabolism. An increase in the content of muscle membrane phospholipids (sphingomyelin, phosphatidylserine, free sterols) was detected in the Seine vs. Douarnenez fish. The data integration allowed to hypothesize relationships between membrane phospholipids, xenobiotic metabolism, bioenergetics, and antioxidant defence. The genetic diversity considering neutral markers was maintained in the heavily polluted Seine population compared with the Douarnenez population. Finally, we suggest that the high physiological cost of tolerance to toxicants in the Seine flounder population could compromise its capacity to respond in the future to an additional stressor like warming waters in shallow depth. Thus, this population could be submitted to an ecological risk.

EcoToxModules: Custom Gene Sets to Organize and Analyze Toxicogenomics Data from Ecological Species

Traditional results from toxicogenomics studies are complex lists of significantly impacted genes or gene sets, which are challenging to synthesize down to actionable results with a clear interpretation. Here, we defined two sets of 21 custom gene sets, called the functional and statistical EcoToxModules, in fathead minnow (Pimephales promelas) to (1) re-cast predefined molecular pathways into a toxicological framework and (2) provide a data-driven, unsupervised grouping of genes impacted by exposure to environmental contaminants. The functional EcoToxModules were identified by re-organizing KEGG pathways into biological processes that are more relevant to ecotoxicology based on the input from expert scientists and regulators. The statistical EcoToxModules were identified using co-expression analysis of publicly available microarray data (n = 303 profiles) measured in livers of fathead minnows after exposure to 38 different conditions. Potential applications of the EcoToxModules were demonstrated with two case studies that represent exposure to a pure chemical and to environmental wastewater samples. In comparisons to differential expression and gene set analysis, we found that EcoToxModule responses were consistent with these traditional results. Additionally, they were easier to visualize and quantitatively compare across different conditions, which facilitated drawing conclusions about the relative toxicity of the exposures within each case study.

Toxicological effects of As (V) in juvenile rockfish Sebastes schlegelii by a combined metabolomic and proteomic approach

Arsenic (As) is a metalloid element that is ubiquitous in the marine environment and its contamination has received worldwide attention due to its potential toxicity. Arsenic can induce multiple adverse effects, such as lipid metabolism disorder, immune system dysfunction, oxidative stress and carcinogenesis, in animals. Inorganic arsenic includes two chemical forms, arsenite (As (III)) and arsenate (As (V)), in natural environment. As (V) is the dominant form in natural waters. In the present study, metabolomic and proteomic alterations were investigated in juvenile rockfish Sebastes schlegelii exposed to environmentally relevant concentrations of As (V) for 14 d. The analysis of iTRAQ-based proteomics combined with untargeted NMR-based metabolomics indicated apparent toxicological effects induced by As (V) in juvenile rockfish. In details, the metabolites, including lactate, alanine, ATP, inosine and phosphocholine were significantly altered in As-treated groups. Proteomic responses suggested that As (V) could not only affected energy and primary metabolisms and signal transduction, but also influenced cytoskeleton structure in juvenile rockfish. This work suggested that the combined proteomic and metabolomic approach could shed light on the toxicological effects of pollutants in rockfish S. schlegelii.

Incorporating Suborganismal Processes into Dynamic Energy Budget Models for Ecological Risk Assessment

A working group at the National Institute for Mathematical and Biological Synthesis (NIMBioS) explored the feasibility of integrating 2 complementary approaches relevant to ecological risk assessment. Adverse outcome pathway (AOP) models provide "bottom-up" mechanisms to predict specific toxicological effects that could affect an individual's ability to grow, reproduce, and/or survive from a molecular initiating event. Dynamic energy budget (DEB) models offer a "top-down" approach that reverse engineers stressor effects on growth, reproduction, and/or survival into modular characterizations related to the acquisition and processing of energy resources. Thus, AOP models quantify linkages between measurable molecular, cellular, or organ-level events, but they do not offer an explicit route to integratively characterize stressor effects at higher levels of organization. While DEB models provide the inherent basis to link effects on individuals to those at the population and ecosystem levels, their use of abstract variables obscures mechanistic connections to suborganismal biology. To take advantage of both approaches, we developed a conceptual model to link DEB and AOP models by interpreting AOP key events as measures of damage-inducing processes affecting DEB variables and rates. We report on the type and structure of data that are generated for AOP models that may also be useful for DEB models. We also report on case studies under development that merge information collected for AOPs with DEB models and highlight some of the challenges. Finally, we discuss how the linkage of these 2 approaches can improve ecological risk assessment, with possibilities for progress in predicting population responses to toxicant exposures within realistic environments. Integr Environ Assess Manag 2018;14:615-624. © 2018 SETAC.

Tissue-Based Mapping of the Fathead Minnow (Pimephales promelas) Transcriptome and Proteome

Omics approaches are broadly used to explore endocrine and toxicity-related pathways and functions. Nevertheless, there is still a significant gap in knowledge in terms of understanding the endocrine system and its numerous connections and intricate feedback loops, especially in non-model organisms. The fathead minnow (Pimephales promelas) is a widely used small fish model for aquatic toxicology and regulatory testing, particularly in North America. A draft genome has been published, but the amount of available genomic or transcriptomic information is still far behind that of other more broadly studied species, such as the zebrafish. Here, we used a proteogenomics approach to survey the tissue-specific proteome and transcriptome profiles in adult male fathead minnow. To do so, we generated a draft transcriptome using short and long sequencing reads from liver, testis, brain, heart, gill, head kidney, trunk kidney, and gastrointestinal tract. We identified 30,378 different putative transcripts overall, with the assembled contigs ranging in size from 264 to over 9,720 nts. Over 17,000 transcripts were >1,000 nts, suggesting a robust transcriptome that can be used to interpret RNA sequencing data in the future. We also performed RNA sequencing and proteomics analysis on four tissues, including the telencephalon, hypothalamus, liver, and gastrointestinal tract of male fish. Transcripts ranged from 0 to 600,000 copies per gene and a large portion were expressed in a tissue-specific manner. Specifically, the telencephalon and hypothalamus shared the most expressed genes, while the gastrointestinal tract and the liver were quite distinct. Using protein profiling techniques, we identified a total of 4,045 proteins in the four tissues investigated, and their tissue-specific expression pattern correlated with the transcripts at the pathway level. Similarly to the findings with the transcriptomic data, the hypothalamus and telencephalon had the highest degree of similarity in the proteins detected. The main purpose of this analysis was to generate tissue-specific omics data in order to support future aquatic ecotoxicogenomic and endocrine-related studies as well as to improve our understanding of the fathead minnow as an ecological model.

Mitochondrial activity in fern spores of Cyathea costaricensis as an indicator of the impact of land use and water quality in rivers running through cloud forests

Early-warning biomarkers, such as mitochondrial activity, have become a key tool in ecosystem assessment. This study aims to evaluate the response of mitochondrial activity in spores of the autochthonous fern Cyathea costaricensis as a bioassessment tool concurrently with land use and physicochemical evaluation in 11 sites along Bobos River, Veracruz, Mexico, to assess river water quality. Bobos River is located in the Nautla basin, northeastern Veracruz (Mexico); the upper river runs through a protected natural area (Filobobos River and adjacent areas). The study involved three monitoring periods: February, June and September 2014. In each study site, physicochemical water quality parameters were recorded to calculate the Water Quality Index (WQI); also, study sites were characterized in terms of land use. Water samples were collected to perform bioassays where spores of C. costaricensis were exposed to samples to assess mitochondrial activity; a positive control exposure test was run under controlled conditions to maximize mitochondrial activity. A Principal Component Analysis was performed to correlate land-use attributes with environmental variables and mitochondrial activity. Three river sections were identified: the upper portion was characterized by the dominance of native vegetation, the highest WQI (in September), and the lowest mitochondrial activity (63.87%-77.47%), related to the geological nature of the basin and high hardness levels. Mitochondrial activity peaked in September (98.32% ± 9.01), likely resulting from nutrient enrichment in the rainy season, and was lowest in February (74.54% ± 1.60) (p < 0.05). Mitochondrial activity was found to be a good benchmark for the assessment of water quality, reflecting the effects of physicochemical characteristics. Mitochondrial activity showed changes along the river and between seasons, associated with environmental characteristics such as land use and the geological nature of the basin, as well as with those related to human impacts.

Toxicological effects of tris(2-chloropropyl) phosphate in human hepatic cells

Organophosphate flame retardants (OPFRs) are widely used as flame retardants which are ubiquitous in various environment media. As many of OPFRs are toxic and persistent, concerns have been raised in regards to their environmental impact. In this study, the toxicological effects of tris(2-chloropropyl) phosphate (TCPP) in human L02 cells was investigated by cell proliferation and apoptosis, oxidative stress, metabolomic and proteomic responses as well as gene expressions related to apoptosis. Results showed that TCPP did not significantly affect the L02 cell apoptosis, however, a significant increase of ROS production was observed in L02 cells with TCPP treatment compared with that in control group (p < 0.05). The expression levels of Bcl-2 family-encoding genes (Bax, Hrk and Bax/Bcl-2) were up-regulated significantly in 10^-4 M group (p < 0.05). Metabolomic and proteomic responses indicated that TCPP mainly caused disturbance in cell growth/division and gene expression, energy and material metabolism, signal transduction, defense and cytoskeleton, which was further confirmed by the western blot analysis.

Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia

Ecological genomics aims to understand the functional association between environmental gradients and the genes underlying adaptive traits. Many genes that are identified by genome-wide screening in ecologically relevant species lack functional annotations. Although gene functions can be inferred from sequence homology, such approaches have limited power. Here, we introduce ecological regulatory genomics by presenting an ontology-free gene prioritization method. Specifically, our method combines transcriptome profiling with high-throughput cis-regulatory sequence analysis in the water fleas Daphnia pulex and Daphnia magna. It screens coexpressed genes for overrepresented DNA motifs that serve as transcription factor binding sites, thereby providing insight into conserved transcription factors and gene regulatory networks shaping the expression profile. We first validated our method, called Daphnia-cisTarget, on a D. pulex heat shock data set, which revealed a network driven by the heat shock factor. Next, we performed RNA-Seq in D. magna exposed to the cyanobacterium Microcystis aeruginosa. Daphnia-cisTarget identified coregulated gene networks that associate with the moulting cycle and potentially regulate life history changes in growth rate and age at maturity. These networks are predicted to be regulated by evolutionary conserved transcription factors such as the homologues of Drosophila Shavenbaby and Grainyhead, nuclear receptors, and a GATA family member. In conclusion, our approach allows prioritising candidate genes in Daphnia without bias towards prior knowledge about functional gene annotation and represents an important step towards exploring the molecular mechanisms of ecological responses in organisms with poorly annotated genomes.

Proteomic and metabolomic responses in D-shape larval mussels Mytilus galloprovincialis exposed to cadmium and arsenic

Cadmium (Cd) and arsenic (As) are the main metal/metalloid contaminants in the coastal environments of the Bohai Sea, China. In this work, a combined proteomic and metabolomic approach was applied to investigate the biological effects of Cd and As (V) in the early life stage (D-shape larvae) of mussel Mytilus galloprovincialis. Results indicated that Cd was a potential immune toxicant to D-shape larval mussel because of the numerous proteomic responses related to immune system. Additionally, Cd induced oxidative stress, cellular injury and disturbance in nucleic acid metabolism in D-shape larval mussels. However, only two identified proteins were significantly altered in As (V)-treated group, suggesting that D-shape larval mussel was less sensitive to As (V) than to Cd at protein level. These two proteins in response to As (V) suggested that As (V) influenced anti-oxidative system and cell proliferation in D-shape larval mussels. Metabolic responses indicated that Cd and As (V) induced disturbances in osmotic regulation and energy metabolism in D-shape larval mussels via different metabolic pathways. In addition, Cd reduced lipid metabolism as well. This work demonstrated that a combination of proteomics and metabolomics could provide an insightful view in the biological effects of pollutants in mussel M. galloprovincialis at an early life stage.

Changing expression of vertebrate immunity genes in an anthropogenic environment: a controlled experiment

Background

The effect of anthropogenic environments on the function of the vertebrate immune system is a problem of general importance. For example, it relates to the increasing rates of immunologically-based disease in modern human populations and to the desirability of identifying optimal immune function in domesticated animals. Despite this importance, our present understanding is compromised by a deficit of experimental studies that make adequately matched comparisons between wild and captive vertebrates.

Results

We transferred post-larval fishes (three-spined sticklebacks), collected in the wild, to an anthropogenic (captive) environment. We then monitored, over 11 months, how the systemic expression of immunity genes changed in comparison to cohort-matched wild individuals in the originator population (total n = 299). We found that a range of innate (lyz, defbl2, il1r-like, tbk1) and adaptive (cd8a, igmh) immunity genes were up-regulated in captivity, accompanied by an increase in expression of the antioxidant enzyme, gpx4a. For some genes previously known to show seasonality in the wild, this appeared to be reduced in captive fishes. Captive fishes tended to express immunity genes, including igzh, foxp3b, lyz, defbl2, and il1r-like, more variably. Furthermore, although gene co-expression patterns (analyzed through gene-by-gene correlations and mutual information theory based networks) shared common structure in wild and captive fishes, there was also significant divergence. For one gene in particular, defbl2, high expression was associated with adverse health outcomes in captive fishes.

Conclusion

Taken together, these results demonstrate widespread regulatory changes in the immune system in captive populations, and that the expression of immunity genes is more constrained in the wild. An increase in constitutive systemic immune activity, such as we observed here, may alter the risk of immunopathology and contribute to variance in health in vertebrate populations exposed to anthropogenic environments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0751-8) contains supplementary material, which is available to authorized users.

Seasonal immunoregulation in a naturally-occurring vertebrate

Background

Fishes show seasonal patterns of immunity, but such phenomena are imperfectly understood in vertebrates generally, even in humans and mice. As these seasonal patterns may link to infectious disease risk and individual condition, the nature of their control has real practical implications. Here we characterize seasonal dynamics in the expression of conserved vertebrate immunity genes in a naturally-occurring piscine model, the three-spined stickleback.

Results

We made genome-wide measurements (RNAseq) of whole-fish mRNA pools (n = 36) at the end of summer and winter in contrasting habitats (riverine and lacustrine) and focussed on common trends to filter habitat-specific from overarching temporal responses. We corroborated this analysis with targeted year-round whole-fish gene expression (Q-PCR) studies in a different year (n = 478). We also considered seasonal tissue-specific expression (6 tissues) (n = 15) at a third contrasting (euryhaline) locality by Q-PCR, further validating the generality of the patterns seen in whole fish analyses. Extremes of season were the dominant predictor of immune expression (compared to sex, ontogeny or habitat). Signatures of adaptive immunity were elevated in late summer. In contrast, late winter was accompanied by signatures of innate immunity (including IL-1 signalling and non-classical complement activity) and modulated toll-like receptor signalling. Negative regulators of T-cell activity were prominent amongst winter-biased genes, suggesting that adaptive immunity is actively down-regulated during winter rather than passively tracking ambient temperature. Network analyses identified a small set of immune genes that might lie close to a regulatory axis. These genes acted as hubs linking summer-biased adaptive pathways, winter-biased innate pathways and other organismal processes, including growth, metabolic dynamics and responses to stress and temperature. Seasonal change was most pronounced in the gill, which contains a considerable concentration of T-cell activity in the stickleback.

Conclusions

Our results suggest major and predictable seasonal re-adjustments of immunity. Further consideration should be given to the effects of such responses in seasonally-occurring disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2701-7) contains supplementary material, which is available to authorized users.

The interactive effects of mercury and selenium on metabolic profiles, gene expression and antioxidant enzymes in halophyte Suaeda salsa

Suaeda salsa is the pioneer halophyte in the Yellow River Delta and was consumed as a popular vegetable. Mercury has become a highly risky contaminant in the sediment of intertidal zones of the Yellow River Delta. In this work, we investigated the interactive effects of mercury and selenium in S. salsa on the basis of metabolic profiling, antioxidant enzyme activities and gene expression quantification. Our results showed that mercury exposure (20 μg L(-1)) inhibited plant growth of S. salsa and induced significant metabolic responses and altered expression levels of INPS, CMO, and MDH in S. salsa samples, together with the increased activities of antioxidant enzymes including SOD and POD. Overall, these results indicated osmotic and oxidative stresses, disturbed protein degradation and energy metabolism change in S. salsa after mercury exposures. Additionally, the addition of selenium could induce both antagonistic and synergistic effects including alleviating protein degradation and aggravating osmotic stress caused by mercury.

Comparative investigations on the biological effects of As (III) and As (V) in clam Ruditapes philippinarum using multiple biomarkers

Inorganic arsenic is a known pollutant with two chemical forms, arsenite (As (III)) and arsenate (As (V)), in marine environment. Clam Ruditapes philippinarum is an important fishery species along the Bohai coast. In this study, the biological effects induced by the two arsenic chemical forms (arsenite and arsenate) were compared using multiple biochemical indices in the digestive glands of clam R. philippinarum. The production of reactive oxygen species, antioxidant enzyme activities and metabolic responses exhibited that both As (III) and As (V) induced immune, oxidative and osmotic stresses in clam digestive glands. The differential metabolic biomarkers, histidine and taurine, indicated the differential responsive mechanisms in osmotic regulation in clam digestive glands. In addition, both arsenic treatments enhanced the anaerobiosis metabolism in clam digestive glands. Overall, this work illustrated that arsenite and arsenate induced similar biological effects in clams, which might be accounted for the biological transformation of arsenate to arsenite in clams.

Systems Biology Approach Reveals a Calcium-Dependent Mechanism for Basal Toxicity in Daphnia magna

The expanding diversity and ever increasing amounts of man-made chemicals discharged to the environment pose largely unknown hazards to ecosystem and human health. The concept of adverse outcome pathways (AOPs) emerged as a comprehensive framework for risk assessment. However, the limited mechanistic information available for most chemicals and a lack of biological pathway annotation in many species represent significant challenges to effective implementation of this approach. Here, a systems level, multistep modeling strategy demonstrates how to integrate information on chemical structure with mechanistic insight from genomic studies, and phenotypic effects to define a putative adverse outcome pathway. Results indicated that transcriptional changes indicative of intracellular calcium mobilization were significantly overrepresented in Daphnia magna (DM) exposed to sublethal doses of presumed narcotic chemicals with log Kow ≥ 1.8. Treatment of DM with a calcium ATPase pump inhibitor substantially recapitulated the common transcriptional changes. We hypothesize that calcium mobilization is a potential key molecular initiating event in DM basal (narcosis) toxicity. Heart beat rate analysis and metabolome analysis indicated sublethal effects consistent with perturbations of calcium preceding overt acute toxicity. Together, the results indicate that altered calcium homeostasis may be a key early event in basal toxicity or narcosis induced by lipophilic compounds.

Immunology in wild nonmodel rodents: an ecological context for studies of health and disease

Transcriptomic methods are set to revolutionize the study of the immune system in naturally occurring nonmodel organisms. With this in mind, the present article focuses on ways in which the use of 'nonmodel' rodents (not the familiar laboratory species) can advance studies into the classical, but ever relevant, epidemiologic triad of immune defence, infectious disease and environment. For example, naturally occurring rodents are an interesting system in which to study the environmental stimuli that drive the development and homeostasis of the immune system and, by extension, to identify where these stimuli are altered in anthropogenic environments leading to the formation of immunopathological phenotypes. Measurement of immune expression may help define individual heterogeneity in infectious disease susceptibility and transmission and facilitate our understanding of infection dynamics and risk in the natural environment; furthermore, it may provide a means of surveillance that can filter individuals carrying previously unknown acute infections of potential ecological or zoonotic importance. Finally, the study of immunology in wild animals may reveal interactions within the immune system and between immunity and other organismal traits that are not observable under restricted laboratory conditions. Potentiating much of this is the possibility of combining gene expression profiles with analytical tools derived from ecology and systems biology to reverse engineer interaction networks between immune responses, other organismal traits and the environment (including symbiont exposures), revealing regulatory architecture. Such holistic studies promise to link ecology, epidemiology and immunology in natural systems in a unified approach that can illuminate important problems relevant to human health and animal welfare and production.

Assessment of the European flounder responses to chemical stress in the English Channel, considering biomarkers and life history traits

A multi-biomarker approach was developed to evaluate responses of European flounder (Platichthys flesus) in three contrasted estuaries over the English Channel: the Canche (pristine site), Tamar (heavy metals and PAHs contamination) and Seine (heavily pollution with a complex cocktail of contaminants). The condition factor and several biomarkers of the immune system, antioxidant enzymes, energetic metabolism and detoxification processes were investigated in young-of-the-year (0+) and one-year-old (1+) flounder. Results underlined the difference between the pristine site and the Seine estuary which showed a lower condition factor, a modulation of the immune system, a higher Cytochrome C oxidase activity, and an up-regulation of BHMT expression. The moderate biomarker responses in the Tamar fish could be linked to the specific contamination context of this estuary. Flounder life history traits were analyzed by otolith microchemistry, in order to depict how the fish use their habitat and thus respond to chemical stress in estuaries.

Does the chronic chemical contamination of a European flounder population decrease its thermal tolerance?

Juvenile flounders (Platichthys flesus), collected in two estuaries with similar temperature regimes (the heavily polluted Seine and the moderately contaminated Vilaine), were submitted to a common garden experiment. After an acclimation period, both populations were challenged by a thermal stress (9-24°C for 15days, then maintenance at 24°C for 19days). The condition factor of the Vilaine fish increased in both conditions, while it decreased for the heated Seine flounders after 34days. The expression of genes related to the energetic metabolism was measured in the liver. The expression levels for ATP-F0 and COII were significantly reduced for heated vs. standard fish from both estuaries, while a decrease of the 12S expression was detected only in heated vs. standard fish from the Seine estuary. Thus, it is suggested that highly contaminated fish from Seine could display a lower tolerance to thermal stress, compared to moderately contaminated fish from Vilaine.

Multiple biomarkers of biological effects induced by cadmium in clam Ruditapes philippinarum

Cadmium (Cd) is a known heavy metal pollutant in the Bohai Sea. Manila clam Ruditapes philippinarum is an important fishery species along the Bohai coast. In this study, the biological effects induced by two concentrations (20 and 200 μg/L) of Cd were characterized using multiple biochemical indices in the digestive glands of clam R. philippinarum. The total hemocyte counts, reactive oxygen species productions and antioxidant enzyme activities exhibited that Cd induced dose-dependent immune and oxidative stresses in clam digestive glands. Metabolic responses indicated that both Cd exposures caused immune stress marked by the elevated branched chain amino acids (valine, leucine and isoleucine), together with the disturbance in energy metabolism. The differential metabolic biomarkers related to osmotic stress, including homarine, betaine, tyrosine and phenylalanine, suggested the differential responsive mechanisms in clam digestive glands induced by Cd exposures. In addition, both Cd treatments enhanced the anaerobiosis metabolism in clam digestive glands via differential metabolic pathways.

Considerations for the use of transcriptomics in identifying the ‘genes that matter’ for environmental adaptation

Transcriptomics has emerged as a powerful approach for exploring physiological responses to the environment. However, like any other experimental approach, transcriptomics has its limitations. Transcriptomics has been criticized as an inappropriate method to identify genes with large impacts on adaptive responses to the environment because: (1) genes with large impacts on fitness are rare; (2) a large change in gene expression does not necessarily equate to a large effect on fitness; and (3) protein activity is most relevant to fitness, and mRNA abundance is an unreliable indicator of protein activity. In this review, these criticisms are re-evaluated in the context of recent systems-level experiments that provide new insight into the relationship between gene expression and fitness during environmental stress. In general, these criticisms remain valid today, and indicate that exclusively using transcriptomics to screen for genes that underlie environmental adaptation will overlook constitutively expressed regulatory genes that play major roles in setting tolerance limits. Standard practices in transcriptomic data analysis pipelines may also be limiting insight by prioritizing highly differentially expressed and conserved genes over those genes that undergo moderate fold-changes and cannot be annotated. While these data certainly do not undermine the continued and widespread use of transcriptomics within environmental physiology, they do highlight the types of research questions for which transcriptomics is best suited and the need for more gene functional analyses. Such information is pertinent at a time when transcriptomics has become increasingly tractable and many researchers may be contemplating integrating transcriptomics into their research programs.

Identification of differentially expressed proteins of brain tissue in response to methamidophos in flounder (Paralichthys olivaceus)

Methamidophos (MAP), an organophosphorus pesticide used around the world, has been associated with a wide spectrum of toxic effects on organisms in the environment. In this study, the flounder Paralichthys olivaceus was subjected to 10 mg/L MAP for 72 h and 144 h, and the morphological and proteomic changes in the brain were observed, analyzed and compared with those in the non-exposed control group. Under the light microscope and transmission electron microscope, MAP had evidently induced changes in or damage to the flounder tissues. Gas chromatography analysis demonstrated that the MAP residues were significantly accumulated in the flounder brain tissues. Proteomic changes in the brain tissue were revealed using two-dimensional gel electrophoresis and 27 protein spots were observed to be significantly changed by MAP exposure. The results indicated that the regulated proteins were involved in immune and stress responses, protein biosynthesis and modification, signal transduction, organismal development, and 50% of them are protease. qRT-PCR was used to further detect the corresponding change of transcription. These data may be beneficial to understand the molecular mechanism of MAP toxicity in flounder, be very useful for MAP-resistance screening in flounder culture. According to our results and analyzing, heat shock protein 90 (HSP90) and granzyme K (GzmK) had taken important part in immune response to MAP-stress and could be biomarkers for MAP-stress in flounder.

Multiple biomarkers of the cytotoxicity induced by BDE-47 in human embryonic kidney cells

Polybrominated diphenyl ethers (PBDEs) are widely used as brominated flame-retardants in a variety of industrial products. Among these PBDEs, 2,2',4,4'-tetra-bromodiphenyl ether (BDE-47) is one of the most predominant congeners inducing multiple toxicities, including hepatotoxicity, neurotoxicity, cytotoxicity, genotoxicity, carcinogenecity and immunotoxicity in human body. In this study, the cytotoxicity of BDE-47 in human embryonic kidney cells (HEK293) was investigated by a set of bioassays, including cell proliferation, apoptosis, oxidative stress and metabolic responses as well as gene expressions related to apoptosis. Results showed that BDE-47 induced an inverted U-shaped curve of cell proliferation in HEK293 cells from 10(-6) to 10(-4) M. Cell apoptosis and ROS overproduction were detected at 10(-5) M of BDE-47 (p<0.05). In addition, the expressions of Bcl-2 family-encoding genes (Bad, Hrk and Bcl-2) increased significantly in 10(-4)M group (p<0.05). Metabolic responses indicated that BDE-47 mainly caused disturbance in energy metabolism marked by differentially altered ethanol, glutathione, creatine, aspartate, UDP-glucose and NAD(+). The increased lactate/alanine ratios indicated the higher reductive state induced by BDE-47 in all exposures confirmed by the overproduction of ROS.

In Silico analysis of perturbed steroidogenesis and gonad growth in fathead minnows (P. promelas) exposed to 17α-ethynylestradiol

The multi-factorial nature of adverse reproductive effects mediated by endocrine disrupting compounds (or EDCs) makes understanding the mechanistic basis of reproductive dysfunction a highly pertinent area of research. As a consequence, a main motivator for continued research is to integrate 'multi-leveled' complexity (i.e., from genes to phenotype) using mathematical methods capable of encapsulating properties of physiological relevance. In this study, an in silico stoichiometric model of piscine steroidogenesis was augmented with a 'biomass' reaction associating the underlying stoichiometry of steroidogenesis with a reaction representative of gonad growth. The ability of the in silico model to predict perturbed steroidogenesis and subsequent effects on gonad growth was tested by exposing reproductively active male and female fathead minnows (Pimephales promelas) to 88 ng/L of the synthetic estrogen, 17α-ethynylestradiol (EE2). The in silico model was parameterized (or constrained) with experimentally quantified concentrations of selected steroid hormones (using mass spectrometry) and fold changes in gene expression (using RT-qPCR) for selected steroidogenic enzyme genes, in gonads of male and female fish. Once constrained, the optimization framework of flux balance analysis (FBA) was used to calculate an optimal flux through the biomass reaction (analogous to gonad growth) and associated steroidogenic flux distributions required to generate biomass. FBA successfully predicted effects of EE2 exposure on fathead minnow gonad growth (%gonadosomatic index or %GSI) and perturbed production of steroid hormones. Specifically, FBA accurately predicted no effects of exposure on male %GSI and a significant reduction for female %GSI. Furthermore, in silico simulations accurately identified disrupted reaction fluxes catalyzing productions of androgens (in male fish) and progestogens (in female fish), an observation which agreed with in vivo experimentation. The analyses presented is the first-ever to successfully associate underlying flux properties of the steroidogenic network with gonad growth in fish, an approach which can incorporate in silico predictions with toxicological risk assessments.

Future water quality monitoring--adapting tools to deal with mixtures of pollutants in water resource management

Environmental quality monitoring of water resources is challenged with providing the basis for safeguarding the environment against adverse biological effects of anthropogenic chemical contamination from diffuse and point sources. While current regulatory efforts focus on monitoring and assessing a few legacy chemicals, many more anthropogenic chemicals can be detected simultaneously in our aquatic resources. However, exposure to chemical mixtures does not necessarily translate into adverse biological effects nor clearly shows whether mitigation measures are needed. Thus, the question which mixtures are present and which have associated combined effects becomes central for defining adequate monitoring and assessment strategies. Here we describe the vision of the international, EU-funded project SOLUTIONS, where three routes are explored to link the occurrence of chemical mixtures at specific sites to the assessment of adverse biological combination effects. First of all, multi-residue target and non-target screening techniques covering a broader range of anticipated chemicals co-occurring in the environment are being developed. By improving sensitivity and detection limits for known bioactive compounds of concern, new analytical chemistry data for multiple components can be obtained and used to characterise priority mixtures. This information on chemical occurrence will be used to predict mixture toxicity and to derive combined effect estimates suitable for advancing environmental quality standards. Secondly, bioanalytical tools will be explored to provide aggregate bioactivity measures integrating all components that produce common (adverse) outcomes even for mixtures of varying compositions. The ambition is to provide comprehensive arrays of effect-based tools and trait-based field observations that link multiple chemical exposures to various environmental protection goals more directly and to provide improved in situ observations for impact assessment of mixtures. Thirdly, effect-directed analysis (EDA) will be applied to identify major drivers of mixture toxicity. Refinements of EDA include the use of statistical approaches with monitoring information for guidance of experimental EDA studies. These three approaches will be explored using case studies at the Danube and Rhine river basins as well as rivers of the Iberian Peninsula. The synthesis of findings will be organised to provide guidance for future solution-oriented environmental monitoring and explore more systematic ways to assess mixture exposures and combination effects in future water quality monitoring.

Toxicological evaluation of two pedigrees of clam Ruditapes philippinarum as bioindicators of heavy metal contaminants using metabolomics

Heavy metal pollution has been of great concern in the Bohai marine environment. Manila clam Ruditapes philippinarum has been used as a bioindicator in marine toxicology. In this study, NMR-based metabolomics was used to ascertain whether there were significant biological differences between two dominant pedigrees (White and Zebra) of clam and evaluate the suitability of two pedigrees for marine environmental toxicology, together with antioxidant enzymatic analysis. Our results indicated that there were significant biological differences between White and Zebra clams based on the metabolic profiles and antioxidant enzyme activities. In details, the metabolic profiles showed higher levels of amino acids and succinate in Zebra clam digestive glands and higher levels of ATP in White clam digestive glands, respectively. The superoxide dismutase activities in control White and Zebra clam samples were significantly different. Additionally, White clam was more sensitive to Cd based on the significant accumulation of Cd, antioxidant enzymatic alterations and sensitive metabolic changes. Overall, we concluded that White clam could be a preferable bioindicator for marine environmental toxicology.

Spatial patterns in markers of contaminant exposure, glucose and glycogen metabolism, and immunological response in juvenile winter flounder (Pseudoplueronectes americanus)

Inshore winter flounder (Pseudoplueronectes americanus) populations in NY, USA have reached record low numbers in recent years, and recruitment into the fishery appears to be limited by survival of post-settlement juvenile fish. In order to identify cellular pathways associated with site-specific variation in condition and mortality, we examined differential mRNA expression in juvenile winter flounder collected from six different bays across a gradient in human population density and sewage inputs. Illumina sequencing of pooled samples of flounder from contrasting degraded sites and less impacted sites was used to guide our choice of targets for qPCR analysis. 253 transcripts of >100bp were differentially expressed, with 60% showing strong homology to mostly teleost sequences within the NCBI database. Based on these data, transcripts representing nine genes of interest associated with contaminant exposure, immune response and glucose and glycogen metabolism were examined by qPCR in individual flounder from each site. Statistically significant site-specific differences were observed in expression of all but one gene, although patterns in expression were complex with only one (vitellogenin), demonstrating a west to east gradient consistent with known loadings of municipal sewage effluent. Principal components analysis (PCA) identified relationships among the genes evaluated. Our data indicate that juvenile winter flounder are responding to estrogenic chemicals in more urbanized coastal bays, and suggests potential mechanistic links between immune response, contaminant exposure and energy metabolism.

Molecular basis of carcinogenicity of tungsten alloy particles

The tungsten alloy of 91% tungsten, 6% nickel and 3% cobalt (WNC 91-6-3) induces rhabdomyosarcoma when implanted into a rat thigh muscle. To investigate whether this effect is species-specific human HSkMc primary muscle cells were exposed to WNC 91-6-3 particles and responses were compared with those from a rat skeletal muscle cell line (L6-C11). Toxicity was assessed by the adenylate kinase assay and microscopy, DNA damage by the Comet assay. Caspase 3 enzyme activity was measured and oligonucleotide microarrays were used for transcriptional profiling. WNC 91-6-3 particles caused toxicity in cells adjacent to the particles and also increased DNA strand breaks. Inhibition of caspase 3 by WNC 91-6-3 occurred in rat but not in human cells. In both rat and human cells, the transcriptional response to WNC 91-6-3 showed repression of transcripts encoding muscle-specific proteins with induction of glycolysis, hypoxia, stress responses and transcripts associated with DNA damage and cell death. In human cells, genes encoding metallothioneins were also induced, together with genes related to angiogenesis, dysregulation of apoptosis and proliferation consistent with pre-neoplastic changes. An alloy containing iron, WNF 97-2-1, which is non-carcinogenic in vivo in rats, did not show these transcriptional changes in vitro in either species while the corresponding cobalt-containing alloy, WNC 97-2-1 elicited similar responses to WNC 91-6-3. Tungsten alloys containing both nickel and cobalt therefore have the potential to be carcinogenic in man and in vitro assays coupled with transcriptomics can be used to identify alloys, which may lead to tumour formation, by dysregulation of biochemical processes.

The SOLUTIONS project: challenges and responses for present and future emerging pollutants in land and water resources management

SOLUTIONS (2013 to 2018) is a European Union Seventh Framework Programme Project (EU-FP7). The project aims to deliver a conceptual framework to support the evidence-based development of environmental policies with regard to water quality. SOLUTIONS will develop the tools for the identification, prioritisation and assessment of those water contaminants that may pose a risk to ecosystems and human health. To this end, a new generation of chemical and effect-based monitoring tools is developed and integrated with a full set of exposure, effect and risk assessment models. SOLUTIONS attempts to address legacy, present and future contamination by integrating monitoring and modelling based approaches with scenarios on future developments in society, economy and technology and thus in contamination. The project follows a solutions-oriented approach by addressing major problems of water and chemicals management and by assessing abatement options. SOLUTIONS takes advantage of the access to the infrastructure necessary to investigate the large basins of the Danube and Rhine as well as relevant Mediterranean basins as case studies, and puts major efforts on stakeholder dialogue and support. Particularly, the EU Water Framework Directive (WFD) Common Implementation Strategy (CIS) working groups, International River Commissions, and water works associations are directly supported with consistent guidance for the early detection, identification, prioritisation, and abatement of chemicals in the water cycle. SOLUTIONS will give a specific emphasis on concepts and tools for the impact and risk assessment of complex mixtures of emerging pollutants, their metabolites and transformation products. Analytical and effect-based screening tools will be applied together with ecological assessment tools for the identification of toxicants and their impacts. The SOLUTIONS approach is expected to provide transparent and evidence-based candidates or River Basin Specific Pollutants in the case study basins and to assist future review of priority pollutants under the WFD as well as potential abatement options.

A metabolomic investigation of the effects of metal pollution in oysters Crassostrea hongkongensis

Metal pollution has been of great concern in the estuaries in Southern China. In this study, metabolic differences between oysters Crassostrea hongkongensis from clean and metal-polluted sites were characterized using NMR-based metabolomics. We collected oyster samples from one clean (Jiuzhen) and two metal polluted sites (Baijiao and Fugong). The metal concentrations in oyster gills indicated that both the Baijiao and Fugong sites were severely polluted by several metals, including Cr, Ni, Cu, Zn, Ag, Cd and Pb. In particular, Cu and Zn were the major contaminants from the Baijiao and Fugong sites. Compared with those oysters from the clean site (JZ), metal pollution in BJ and FG induced disturbances in osmotic regulation and energy metabolism via different metabolic pathways, as indicated by different metabolic biomarkers. This study demonstrates that NMR-based metabolomics is a useful tool for characterizing metabolic responses induced by metal pollution.

A 2-DE-based proteomic study on the toxicological effects of cisplatin in L02 cells

Cisplatin is a chemotherapeutic agent for the treatment of various cancers. In this study, cisplatin-induced effects were characterized in vitro model of human liver cells (L02) using 2-DE-based proteomics. Results indicated that different cisplatin treatments primarily induced disturbances in protein synthesis and oxidative stress via differential mechanisms. Since the experimental concentrations of cisplatin described a hormesis effect in cell proliferation of L02 cells, it was expected to reveal the hormesis effects using proteomic markers. However, only confilin-1 was commonly up-regulated in three concentrations of cisplatin treatments showing a hormesis effects with a U-shape regulation. These results were highly consistent with many other toxico-proteomic studies, indicating that the toxico-proteomic responses based on dose-dependent protein responses were incongruent with the theoretically linear or hormetic concentration-effect relationship. Our findings suggested that a macroscopic hormesis phenomenon on the cell proliferation could not be reflected by proteomic responses induced by cisplatin treatments.

Cross-species toxicogenomic analyses and phenotypic anchoring in response to groundwater low-level pollution

Background

Comparison of toxicogenomic data facilitates the identification of deregulated gene patterns and maximizes health risk prediction in human.

Results

Here, we performed phenotypic anchoring on the effects of acute exposure to low-grade polluted groundwater using mouse and zebrafish. Also, we evaluated two windows of chronic exposure in mouse, starting in utero and at the end of lactation. Bioinformatic analysis of livers microarray data showed that the number of deregulated biofunctions and pathways is higher after acute exposure, compared to the chronic one. It also revealed specific profiles of altered gene expression in all treatments, pointing to stress response/mitochondrial pathways as major players of environmental toxicity. Of note, dysfunction of steroid hormones was also predicted by bioinformatic analysis and verified in both models by traditional approaches, serum estrogens measurement and vitellogenin mRNA determination in mice and zebrafish, respectively.

Conclusions

In our report, phenotypic anchoring in two vertebrate model organisms highlights the toxicity of low-grade pollution, with varying susceptibility based on exposure window. The overlay of zebrafish and mice deregulated pathways, more than single genes, is useful in risk identification from chemicals implicated in the observed effects.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1067) contains supplementary material, which is available to authorized users.

Spatial variation in transcript and protein abundance of Atlantic salmon during feeding migration in the Baltic Sea

The fitness and reproductive output of fishes can be affected by environmental disturbances. In this study, transcriptomics and label-free proteomics were combined to investigate Atlantic salmon (Salmo salar) sampled from three different field locations within the Baltic Sea (Baltic Main Basin (BMB), Gulf of Finland (GoF), and Bothnian Sea (BS)) during marine migration. The expression of several stress related mRNAs and proteins of xenobiotic metabolism, oxidative stress, DNA damage, and cell death were increased in salmon from GoF compared to salmon from BMB or BS. Respiratory electron chain and ATP synthesis related gene ontology-categories were upregulated in GoF salmon, whereas those associated with RNA processing and synthesis, translation, and protein folding decreased. Differences were seen also in metabolism and immune function related gene expression. Comparisons of the transcriptomic and proteomic profiles between salmon from GoF and salmon from BMB or BS suggest environmental stressors, especially exposure to contaminants, as a main explanation for differences. Salmon feeding in GoF are thus “disturbed by hazardous substances”. The results may also be applied in evaluating the conditions of pelagic ecosystems in the different parts of Baltic Sea.

Responses of the European flounder (Platichthys flesus) to a mixture of PAHs and PCBs in experimental conditions

A multibiomarker approach was developed to evaluate the juvenile European flounder responses to a complex mixture of 9 polycyclic aromatic hydrocarbons (PAHs) and 12 polychlorinated biphenyls (PCBs). Exposure was performed through contaminated food pellets displaying: (1) PAH and PCB levels similar to those detected in the heavily polluted Seine estuary, respectively in sediments and in flatfish and (2) ten times these concentrations. Several biomarkers of the immune system (e.g., lysozyme concentration and gene expression of complement component C3 and TNF-receptor), DNA damage (e.g., Comet assay), energetic metabolism (e.g., activity of cytochrome C oxidase), detoxification process (e.g., cytochrome P450 1A1 expression level: CYP1A1; betaine homocysteine methyl transferase expression level: BHMT) were investigated after 14 and 29 days of contamination, followed by a 14-days recovery period. After 29 days of contamination, the detoxification activity (CYP1A1 expression level) was positively correlated with DNA damages; the increase of the BHMT expression level could also be related to the detoxification process. Furthermore, after the recovery period, some biomarkers were still upregulated (i.e., CYP1A1 and BHMT expression levels). The immune system was significantly modulated by the chemical stress at the two concentration levels, and the lysozyme appeared to be the most sensitive marker of the mixture impact.

Gender-specific metabolic responses in hepatopancreas of mussel Mytilus galloprovincialis challenged by Vibrio harveyi

Mussel Mytilus galloprovincialis is a marine aquaculture shellfish and frequently studied in shellfish immunology. In this work, the gender-specific metabolic responses induced by Vibrio harveyi in hepatopancreas from M. galloprovincialis were characterized using NMR-based metabolomics. In details, V. harveyi challenge increased the levels of amino acids including (valine, leucine, isoleucine, threonine, alanine, arginine and tyrosine) and ATP, and decreased the level of glucose in male mussel hepatopancreas. In V. harveyi-challenged female mussel hepatopancreas, both threonine and AMP were significantly elevated, and choline, phoshphocholine, sn-glycero-3-phosphocholine, taurine, betaine and ATP were depleted. Obviously, only threonine was similarly altered to that in V. harveyi-challenged male mussel hepatopancreas. These findings confirmed the gender-specific metabolic responses in mussels challenged by V. harveyi. Overall, V. harveyi induced an enhanced energy demand through activated glycolysis and immune response indicated by increased BCAAs in male mussel hepatopancreas. In female mussel hepatopancreas, V. harveyi basically caused disturbances in both osmotic regulation and energy metabolism through the metabolic pathways of conversions of phosphocholine and ADP to choline and ATP, and sn-glycero-3-phosphocholine and H2O into choline and sn-glycerol 3-phosphate. The altered mRNA expression levels of related genes (Cu/Zn-SOD, HSP90, lysozyme and defensin) suggested that V. harveyi induced obvious oxidative and immune stresses in both male and female mussel hepatopancreas. This work demonstrated that V. harveyi could induce gender-specific metabolic responses in mussel M. galloprovincialis hepatopancreas using NMR-based metabolomics.

Bioremediation in marine ecosystems: a computational study combining ecological modeling and flux balance analysis

The pressure to search effective bioremediation methodologies for contaminated ecosystems has led to the large-scale identification of microbial species and metabolic degradation pathways. However, minor attention has been paid to the study of bioremediation in marine food webs and to the definition of integrated strategies for reducing bioaccumulation in species. We propose a novel computational framework for analysing the multiscale effects of bioremediation at the ecosystem level, based on coupling food web bioaccumulation models and metabolic models of degrading bacteria. The combination of techniques from synthetic biology and ecological network analysis allows the specification of arbitrary scenarios of contaminant removal and the evaluation of strategies based on natural or synthetic microbial strains. In this study, we derive a bioaccumulation model of polychlorinated biphenyls (PCBs) in the Adriatic food web, and we extend a metabolic reconstruction of Pseudomonas putida KT2440 (iJN746) with the aerobic pathway of PCBs degradation. We assess the effectiveness of different bioremediation scenarios in reducing PCBs concentration in species and we study indices of species centrality to measure their importance in the contaminant diffusion via feeding links. The analysis of the Adriatic sea case study suggests that our framework could represent a practical tool in the design of effective remediation strategies, providing at the same time insights into the ecological role of microbial communities within food webs.

Molecular responses of European flounder (Platichthys flesus) chronically exposed to contaminated estuarine sediments

Molecular responses to acute toxicant exposure can be effective biomarkers, however responses to chronic exposure are less well characterised. The aim of this study was to determine chronic molecular responses to environmental mixtures in a controlled laboratory setting, free from the additional variability encountered with environmental sampling of wild organisms. Flounder fish were exposed in mesocosms for seven months to a contaminated estuarine sediment made by mixing material from the Forth (high organics) and Tyne (high metals and tributyltin) estuaries (FT) or a reference sediment from the Ythan estuary (Y). Chemical analyses demonstrated that FT sediment contained significantly higher concentrations of key environmental pollutants (including polycyclic aromatic hydrocarbons (PAHs), chlorinated biphenyls and heavy metals) than Y sediment, but that chronically exposed flounder showed a lack of differential accumulation of contaminants, including heavy metals. Biliary 1-hydroxypyrene concentration and erythrocyte DNA damage increased in FT-exposed fish. Transcriptomic and (1)H NMR metabolomic analyses of liver tissues detected small but statistically significant alterations between fish exposed to different sediments. These highlighted perturbance of immune response and apoptotic pathways, but there was a lack of response from traditional biomarker genes. Gene-chemical association annotation enrichment analyses suggested that polycyclic aromatic hydrocarbons were a major class of toxicants affecting the molecular responses of the exposed fish. This demonstrated that molecular responses of sentinel organisms can be detected after chronic mixed toxicant exposure and that these can be informative of key components of the mixture.

Metabolic profiling of the tissue-specific responses in mussel Mytilus galloprovincialis towards Vibrio harveyi challenge

Mussel Mytilus galloprovincialis is a marine aquaculture shellfish distributing widely along the coast in north China. In this work, we studied the differential metabolic responses induced by Vibrio harveyi in digestive gland and gill tissues from M. galloprovincialis using NMR-based metabolomics. The differential metabolic responses in the two tissue types were detected, except the similarly altered taurine and betaine. These metabolic responses suggested that V. harveyi mainly induced osmotic disruption and reduced energy demand via the metabolic pathways of glucose synthesis and ATP/AMP conversion in mussel digestive gland. In mussel gill tissues, V. harveyi basically caused osmotic stress and possible reduced energy demand as shown by the elevated phosphocholine that is involved in one of the metabolic pathways of ATP synthesis from ADP and phosphocholine. The altered mRNA expression levels of related genes (superoxide dismutase with copper and zinc, heat shock protein 90, defensin and lysozyme) suggested that V. harveyi induced clear oxidative and immune stresses in both digestive gland and gill tissues. However, the mRNA expression levels of both lysozyme and defensin in digestive gland were more significantly up-regulated than those in gill from V. harveyi-challenged mussel M. galloprovincialis, meaning that the immune organ, digestive gland, was more sensitive than gill. Overall, our results indicated that V. harveyi could induce tissue-specific metabolic responses in mussel M. galloprovincialis.

Transcriptome-metabolome wide association study (TMWAS) of maneb and paraquat neurotoxicity reveals network level interactions in toxicologic mechanism

A combination of the herbicide paraquat (PQ) and fungicide maneb (MB) has been linked to Parkinson's disease. Previous studies show that this involves an additive toxicity with at least two different mechanisms. However, detailed understanding of mixtures is often difficult to elucidate because of the multiple ways by which toxic agents can interact. In the present study, we used a combination of transcriptomics and metabolomics to investigate mechanisms of toxicity of PQ and MB in a neuroblastoma cell line. Conditions were studied with concentrations of PQ and MB that each individually caused 20% cell death and together caused 50% cell death. Transcriptomic and metabolomic samples were collected at time points prior to significant cell death. Statistical and bioinformatic methods were applied to the resulting 30,869 transcripts and 1358 metabolites. Results showed that MB significantly changed more transcripts and metabolites than PQ, and combined PQ + MB impacted more than MB alone. Transcriptome-metabolome-wide association study (TMWAS) showed that significantly changed transcripts and metabolites mapped to two network substructures, one associating with significant effects of MB and the other included features significantly associated with PQ + MB. The latter contained 4 clusters of genes and associated metabolites, with one containing genes for two cation transporters and a cation transporter regulatory protein also recognized as a pro-apoptotic protein. Other clusters included stress response genes and transporters linked to cytoprotective mechanisms. MB also had a significant network structure linked to cell proliferation. Together, the results show that the toxicologic mechanism of the combined neurotoxicity of PQ and MB involves network level interactions and that TMWAS provides an effective approach to investigate such complex mechanisms.

Endocrinology: advances through omics and related technologies

The rapid development of new omics technologies to measure changes at genetic, transcriptomic, proteomic, and metabolomics levels together with the evolution of methods to analyze and integrate the data at a systems level are revolutionizing the study of biological processes. Here we discuss how new approaches using omics technologies have expanded our knowledge especially in nontraditional models. Our increasing knowledge of these interactions and evolutionary pathway conservation facilitates the use of nontraditional species, both invertebrate and vertebrate, as new model species for biological and endocrinology research. The increasing availability of technology to create organisms overexpressing key genes in endocrine function allows manipulation of complex regulatory networks such as growth hormone (GH) in transgenic fish where disregulation of GH production to produce larger fish has also permitted exploration of the role that GH plays in testis development, suggesting that it does so through interactions with insulin-like growth factors. The availability of omics tools to monitor changes at nearly any level in any organism, manipulate gene expression and behavior, and integrate data across biological levels, provides novel opportunities to explore endocrine function across many species and understand the complex roles that key genes play in different aspects of the endocrine function.

A design automation framework for computational bioenergetics in biological networks

The bioenergetic activity of mitochondria can be thoroughly investigated by using computational methods. In particular, in our work we focus on ATP and NADH, namely the metabolites representing the production of energy in the cell. We develop a computational framework to perform an exhaustive investigation at the level of species, reactions, genes and metabolic pathways. The framework integrates several methods implementing the state-of-the-art algorithms for many-objective optimization, sensitivity, and identifiability analysis applied to biological systems. We use this computational framework to analyze three case studies related to the human mitochondria and the algal metabolism of Chlamydomonas reinhardtii, formally described with algebraic differential equations or flux balance analysis. Integrating the results of our framework applied to interacting organelles would provide a general-purpose method for assessing the production of energy in a biological network.