Ecotoxicogenomics: bridging the gap between genes and populations

Development of an effective QSAR-based hazard threshold prediction model for the ecological risk assessment of aromatic hydrocarbon compounds

The insufficient hazard thresholds of specific individual aromatic hydrocarbon compounds (AHCs) with diverse structures limit their ecological risk assessment. Thus, herein, quantitative structure-activity relationship (QSAR) models for estimating the hazard threshold of AHCs were developed based on the hazardous concentration for 5% of species (HC5) determined using the optimal species sensitivity distribution models and on the molecular descriptors calculated via the PADEL software and ORCA software. Results revealed that the optimal QSAR model, which involved eight descriptors, namely, Zagreb, GATS2m, VR3_Dzs, AATSC2s, GATS2c, ATSC2i, ω, and Vm, displayed excellent performance, as reflected by an optimal goodness of fit (R2adj = 0.918), robustness (Q2LOO = 0.869), and external prediction ability (Q2F1 = 0.760, Q2F2 = 0.782, and Q2F3 = 0.774). The hazard thresholds estimated using the optimal QSAR model were approximately close to the published water quality criteria developed by different countries and regions. The quantitative structure-toxicity relationship demonstrated that the molecular descriptors associated with electrophilicity and topological and electrotopological properties were important factors that affected the risks of AHCs. A new and reliable approach to estimate the hazard threshold of ecological risk assessment for various aromatic hydrocarbon pollutants was provided in this study, which can be widely popularised to similar contaminants with diverse structures.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

Genes-to-Pathways Species Conservation Analysis: Enabling the Exploration of Conservation of Biological Pathways and Processes Across Species

The last two decades have witnessed a strong momentum toward integration of cell-based and computational approaches in safety assessments. This is fueling a global regulatory paradigm shift toward reduction and replacement of the use of animals in toxicity tests while promoting the use of new approach methodologies. The understanding of conservation of molecular targets and pathways provides an opportunity to extrapolate effects across species and ultimately to determine the taxonomic applicability domain of assays and biological effects. Despite the wealth of genome-linked data available, there is a compelling need for improved accessibility, while ensuring that it reflects the underpinning biology. We present the novel pipeline Genes-to-Pathways Species Conservation Analysis (G2P-SCAN) to further support understanding on cross-species extrapolation of biological processes. This R package extracts, synthetizes, and structures the data available from different databases, that is, gene orthologs, protein families, entities, and reactions, linked to human genes and respective pathways across six relevant model species. The use of G2P-SCAN enables the overall analysis of orthology and functional families to substantiate the identification of conservation and susceptibility at the pathway level. In the present study we discuss five case studies, demonstrating the validity of the developed pipeline and its potential use as species extrapolation support. We foresee this pipeline will provide valuable biological insights and create space for the use of mechanistically based data to inform potential species susceptibility for research and safety decision purposes. Environ Toxicol Chem 2023;42:1152-1166. © 2023 UNILEVER GLOBAL IP LTD. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Gene expression changes in Epinephelus marginatus (Teleostei, Serranidae) liver reveals candidate molecular biomarker of iron ore contamination

Wastes from iron ore mining activities are potentially damaging to adjacent aquatic ecosystems. We aimed to determine biomarkers of environmental exposure to this xenobiotic in the dusky grouper Epinephelus marginatus by differential gene expression analysis. For this, fish were exposed to iron ore (15.2 mg/L) and gene expression in liver was assessed by RNA-Seq and compared to the control group. A total of 124 differentially expressed genes were identified, from which 52 were upregulated and 72 were downregulated in response to iron ore. From these, ferritin (medium subunit), cytochrome b reductase and epoxide hydrolase genes were selected for validation by RT-qPCR that confirmed the upregulation of epoxide hydrolase in fish exposed to iron ore.

Analysis of multi-omics data on the relationship between epigenetic changes and nervous system disorders caused by exposure to environmentally harmful substances

Environmentally hazardous substances and exposure to these can cause various diseases. Volatile organic compounds can easily evaporate into the atmosphere, thereby exerting toxic effects through either the skin or respiratory tract exposures. Toluene, a neurotoxin, has been widely used in various industries. However, it has a detrimental effect on the nervous system (such as hallucinations or memory impairment), while data on the mechanism underlaying its harmful effects remain limited. Therefore, this study investigates the effect of toluene on the nervous system via epigenetic and genetic changes of toluene-exposed individuals. We identified significant epigenetic changes and confirmed that the affected abnormally expressed genes negatively influenced the nervous system. In particular, we confirmed that the miR-15 family, upregulated by toluene, downregulated ABL2, which could affect the R as signaling pathway resulting in neuronal structural abnormalities. Our study suggests that miR-15a-5p, miR-15b-5p, miR-16-5p, miR-301a-3p, and lncRNA NEAT1 may represent effective epigenomic markers associated with neurodegenerative diseases caused by toluene.

Chironomidae (Midge) Sensitivities to Ammonia Using Multiple Endpoints in China and Australia for the Development of Water Quality Criteria for Freshwater River Systems in China

Deriving water quality criteria (WQC) for aquatic risk assessment requires sufficient toxicity data, which can determine the accuracy of WQC. Given that toxicity data vary between test species and endpoints, there is a great need to compare such data to generate the most suitable data set for WQC derivation. In the present study, a series of 11 ammonia exposure bioassays were conducted on Chironomidae species in either China or Australia, with test species and test endpoints varied (2 Chironomus sp., enzymatic up to lethal endpoints, and no-observed-effect concentration up to median lethal concentration [LC50] as endpoint metrics). There were no statistically significant differences between toxicity results generated from China compared to Australia using Chironomus sp., indicating that published data on native species generated in different countries could be appropriate for inclusion in the development of local Chinese WQC. In addition, the Chironomidae larvae laboratory-based toxicity value (LC50 = 384.6 mg/L) was lower than that of the in situ field-based toxicity value (LC50 ≥ 451.2 mg/L) where sensitive life stages are used, and, specifically for C. riparius, endpoints linked to biochemical and gene expression effects could be as sensitive as or more sensitive than chronic endpoints, both of which were more sensitive than acute endpoints. These findings help in the development of WQC by demonstrating the suitability of inclusion of toxicity data from a range of sources, as well as adding to the overall pool of knowledge regarding sensitivity to ammonia which can be used in aquatic risk assessment. Environ Toxicol Chem  2021;40:2899-2911. © 2021 SETAC.

Signaling pathway perturbation analysis for assessment of biological impact of cigarette smoke on lung cells

Exposure to cigarette smoke (CS) results in injury to the epithelial cells of the human respiratory tract and has been implicated as a causative factor in the development of chronic obstructive pulmonary disease and lung cancers. The application of omics-scale methodologies has improved the capacity to understand cellular signaling processes underlying response to CS exposure. We report here the development of an algorithm based on quantitative assessment of transcriptomic profiles and signaling pathway perturbation analysis (SPPA) of human bronchial epithelial cells (HBEC) exposed to the toxic components present in CS. HBEC were exposed to CS of different compositions and for different durations using an ISO3308 smoking regime and the impact of exposure was monitored in 2263 signaling pathways in the cell to generate a total effect score that reflects the quantitative degree of impact of external stimuli on the cells. These findings support the conclusion that the SPPA algorithm provides an objective, systematic, sensitive means to evaluate the biological impact of exposures to CS of different compositions making a powerful comparative tool for commercial product evaluation and potentially for other known or potentially toxic environmental smoke substances.

Multiple bioanalytical methods reveal a thyroid-disrupting mechanism related to the membrane receptor integrin αvβ3

Di(2-ethylhexyl) phthalate (DEHP), a manufactured chemical, is suitable for large-scale production and has extensive applications. Although restricted for use, DEHP is still ubiquitous in the environment and shows potential to disrupt the structure or function of the thyroid system. However, its toxic mechanism is complex and not clearly understood. In this study, a battery of methods was employed to investigate DEHP-induced thyroid-disrupting effects and their mechanism of action, focusing on a newly discovered membrane receptor-mediated mechanism. The results showed that DEHP promoted rat pituitary tumor (GH3) cell proliferation and c-fos gene expression at environment level concentrations (2 and 5 μmol/L) in a manner similar to that of the natural thyroid hormone 3,3',5-triiodo-L-thyronine (T3). The macromolecule DEHP-BSA cannot pass through the cell membrane to interact with nuclear receptors but upregulated the c-fos gene expression when administered at concentrations comparable to DEHP concentrations; molecular docking demonstrated that DEHP has affinity for the membrane receptor integrin α_vβ₃; DEHP at 2 μmol/L upregulated the β₃ gene expression in GH3 cells; after the addition of integrin α_vβ₃-inhibiting RGD peptide, DEHP-induced c-fos gene upregulation decreased. All of these findings support the supposition that DEHP-induced thyroid-disrupting effects might be mediated by the membrane receptor integrin α_vβ₃. Moreover, DEHP activated the downstream extracellular regulated protein kinase (ERK1/2) pathway, upregulating the gene expression of raf-1, MEK-1 and MAPK1 and increasing the protein levels of p-ERK; interestingly, ERK1/2 activation and c-fos upregulation induced by DEHP were attenuated by PD98059 (an ERK1/2 inhibitor). Taken together, the data suggest that the membrane receptor integrin α_vβ₃ and the downstream ERK1/2 pathway might be involved in DEHP-induced thyroid-disrupting effects. This study provides new insight into the thyroid-disrupting effect and the underlying mechanism and will advance the effort to construct adverse outcome pathways of DEHP and other thyroid hormone disrupting chemicals.

Ecotoxicology of persistent organic pollutants in birds

Considering the explosive growth of the list of persistent organic pollutants (POPs), the scientific community is combatting increasing challenges to protect humans and wildlife from the potentially negative consequences of POPs. Herein, we characterize the main aspects and progress in the ecotoxicology of POPs in avian species since 2000. The majority of previous efforts has revealed the global occurrence of high levels of various POPs in birds. Laboratory research and epidemiological studies imply that POPs exert a broad-spectrum of side-effects on birds by interfering with their endocrine, immune and neural system, reproduction, and development, and growth. However, inconsistent results suggest that the potential effects of POP exposure on the physiological parameters in birds are multifactorial, involving a multitude of biological processes, species-specific differences, gender, age and types of compounds. Great progress has been achieved in identifying the species-specific sensitivity to dioxin-like compounds, which is attributed to different amino acid residues in the ligand-binding domain of the aryl hydrocarbon receptor. Besides the conventional concentration additivity, several studies have suggested that different classes of POPs possibly act synergistically or antagonistically based on their concentration. However, ecotoxicology information is still recorded in a scattered and inadequate manner, including lack of enough avian species, limited number of POPs investigated, and insufficient geographical representation, and thus our understanding of the effects of POPs on birds remains rudimentary, although mechanistic understanding of their mode of action is progressing. Particularly, research on what happens to wild bird populations and their ecosystems under POP stress is still unavailable. Thus, our aim is to predict and trace the effects POPs at different biological organization levels, especially from the molecular, cellular and individual levels to the population, community and ecosystem levels because of the limited and scattered information, as mentioned above.

Cross-species extrapolation of chemical sensitivity

Ecosystems are usually populated by many species. Each of these species carries the potential to show a different sensitivity towards all of the numerous chemical compounds that can be present in their environment. Since experimentally testing all possible species-chemical combinations is impossible, the ecological risk assessment of chemicals largely depends on cross-species extrapolation approaches. This review overviews currently existing cross-species extrapolation methodologies, and discusses i) how species sensitivity could be described, ii) which predictors might be useful for explaining differences in species sensitivity, and iii) which statistical considerations are important. We argue that risk assessment can benefit most from modelling approaches when sensitivity is described based on ecologically relevant and robust effects. Additionally, specific attention should be paid to heterogeneity of the training data (e.g. exposure duration, pH, temperature), since this strongly influences the reliability of the resulting models. Regarding which predictors are useful for explaining differences in species sensitivity, we review interspecies-correlation, relatedness-based, traits-based, and genomic-based extrapolation methods, describing the amount of mechanistic information the predictors contain, the amount of input data the models require, and the extent to which the different methods provide protection for ecological entities. We develop a conceptual framework, incorporating the strengths of each of the methods described. Finally, the discussion of statistical considerations reveals that regardless of the method used, statistically significant models can be found, although the usefulness, applicability, and understanding of these models varies considerably. We therefore recommend publication of scientific code along with scientific studies to simultaneously clarify modelling choices and enable elaboration on existing work. In general, this review specifies the data requirements of different cross-species extrapolation methods, aiming to make regulators and publishers more aware that access to raw- and meta-data needs to be improved to make future cross-species extrapolation efforts successful, enabling their integration into the regulatory environment.

Visual system: An understudied target of aquatic toxicology

Visual system is increasingly recognized as a sensitive target of xenobiotics in aquatic ecosystems. Various environmental pollutants of distinct physicochemical properties are able to impair the retinal development and function of teleost fishes, including dioxin-like pollutants, flame retardants, pesticides, perfluoroalkyl acids, retinoic acids and metals. Considering the availability of developmental and functional database, zebrafish has been the most frequently used as the teleost model to study aquatic visual toxicology. A diversity of visual deficits has been displayed for fishes across multiple levels of biological organizations (e.g., molecule, cell, histology, physiology and behavior). Covering sensitive developmental windows of eyes during early embryogenesis, acute or chronic exposure to xenobiotics can disturb the expressions of visual gene and protein markers, which affect the retinal neurogenesis and induce degeneration of neurons. Morphological structures and physiological responses of retina and optic tectum are then disorganized, eventually compromising the performance of visually-mediated behaviors and recruitment of individuals. Environmental pollutants can cross the blood-retina barrier and accumulate in eyes, which might impact visual system directly. In addition, pollutants are very likely to interrupt retinal development and function indirectly by disturbing the signaling of retinoids and thyroid. However, exact mechanisms of visual toxicity are largely unknown currently. In this review, the development and structure of retina and available tools for studying visual science are described briefly. Advances in visual toxicology are summarized in detail and outlooks for future visual toxicity studies are discussed.

Multigene Biomarkers of Pyrethroid Exposure: Exploratory Experiments

We describe initial development of microarray-based assays for detecting 4 pyrethroid pesticides (bifenthrin, cypermethrin, esfenvalerate, and permethrin) in water. To facilitate comparison of transcriptional responses with gross apical responses, we estimated concentration-mortality curves for these pyrethroids using flow-through exposures of newly hatched Daphnia magna, Pimephales promelas adults, and 24 h posthatch P. promelas. Median lethal concentration (LC50) estimates were below most reported values, perhaps attributable to the use of flow-through exposures or of measured rather than nominal concentrations. Microarray analysis of whole P. promelas larvae and brains from exposed P. promelas adults showed that assays using either tissue type can detect these pyrethroids at concentrations below LC50 values reported for between 72 and 96% of aquatic species, depending on the pesticide. These estimates are conservative because they correspond to the lowest concentrations tested. This suggests that the simpler and less expensive whole-larval assay provides adequate sensitivity for screening contexts where acute aquatic lethality is observed, but the responsible agent is not known. Gene set analysis (GSA) highlighted several Gene Ontology (GO) terms consistent with known pyrethroid action, but the implications of other GO terms are less clear. Exploration of the sensitivity of results to changes in data processing suggests robustness of the detection assay results, but GSA results were sensitive to methodological variations. Environ Toxicol Chem 2019;38:2436-2446. Published 2019 Wiley Periodicals, Inc. on behalf of SETAC. This article is a US government work, and as such, is in the public domain in the United States of America.

Towards early-warning gene signature of Chlamydomonas reinhardtii exposed to Hg-containing complex media

The potential of using gene expression signature as a biomarker of toxicants exposure was explored in the microalga Chlamydomonas reinhardtii exposed 2 h to mercury (Hg) as inorganic mercury (IHg) and methyl mercury (MeHg) in presence of copper (Cu) and Suwannee River Humic Acid (SRHA). Total cellular Hg (THg = IHg + MeHg) decreased in presence of SRHA for 0.7 nM IHg and 0.4 nM MeHg, but increased for 70 nM IHg exposure. In mixtures of IHg + MeHg and (IHg or MeHg) + Cu, SRHA decreased THg uptake, except for 0.7 nM IHg + 0.4 nM MeHg which was unchanged (p-value>0.05). In the absence of SRHA, 0.5 μM Cu strongly decreased intracellular THg concentration for 70 nM IHg, while it had no effect for 0.7 nM IHg and 0.4 nM MeHg. The expression of single transcripts was not correlated with measured THg uptake, but a subset of 60 transcripts showed signatures specific to the exposed metal(s) and was congruent with exposure concentration. Notably, the range of fold change values of this subset correlated with THg bioaccumulation with a two-slope pattern in line with [THg]_intra/[THg]_med ratios. Gene expression signature seems a promising approach to complement chemical analyses to assess bioavailability of toxicants in presence of other metals and organic matter.

Differential modification of the C. elegans proteome in response to acute and chronic gamma radiation: Link with reproduction decline

Emission of ionizing radiation (IR) in the environment is a natural phenomenon which can be enhanced by human activities. Ecosystems are then chronically exposed to IR. But environmental risk assessment of chronic exposure suffers from a lack of knowledge. Extrapolation of data from acute to chronic exposure is not always relevant, and can lead to uncertainties as effects could be different between the two irradiation modes, especially regarding reproduction endpoint, which is an ecologically relevant parameter. In the present study, we decided to refine the understanding of the molecular mechanisms involved in response to acute and chronic γ-irradiation by a global proteome label free LC-MS/MS analysis. C. elegans were exposed to 3 common cumulated radiation doses for acute or chronic exposure condition and global modification of the proteome was studied. This analysis of protein expression has demonstrated the modulation of proteins involved in regulatory biological processes such as lipid transport, DNA replication, germ cell development, apoptosis, ion transport, cuticle development, and aging at lower doses than those for which individual effects on reproduction have been previously observed. Thus, these proteins could constitute early and sensitive markers of radio-induced reprotoxicity; more specifically HAT-1, RPS-19 in acute and VIT-3 for chronic conditions that are expressed in a dose-dependent manner. Finally, to focus on reproduction process, this analysis showed either repression or overexpression of 12 common proteins in organisms exposed to acute or chronic irradiation, respectively. These proteins include the vitellogenin cluster notably involved in lipid transport and oocyte maturation and proteins involved in cuticle development and molting i.e. COL-14, GLF-1, NOAH-1, NOAH-2, ACN-1. These results show that protein expression modulation is a sensitive and predictive marker of radio-induced reproductive effects, but also highlight limitation of data extrapolation from acute to chronic exposure for environmental risk assessment.

Variation in transcriptional responses to copper exposure across Daphnia pulex lineages

Copper pollution is pervasive in aquatic habitats and is particularly harmful to invertebrates sensitive to environmental changes such as Daphnia pulex. Mechanisms of toxicity and tolerance to copper are not well understood. We used RNA-sequencing to investigate these mechanisms in three genetically distinct D. pulex clonal lineages with different histories of copper exposure. Upregulated genes after copper exposure were enriched with Gene Ontology (GO) categories involved in digestion, molting and growth, whereas downregulated genes after copper exposure were enriched in the metal-regulatory system, immune response and epigenetic modifications. The three D. pulex clones in our study show largely similar transcriptional patterns in response to copper, with only a total of twenty genes differentially expressed in a single clonal lineages. We also detected lower relative expression of some genes known to be important for copper tolerance, metallothionein and glutathione-S-transferase, in a sensitive lineage sampled from an uncontaminated habitat. Daphnia-specific genes (without orthologs outside the genus) and Daphnia-specific duplications (genes duplicated in the Daphnia lineage) were overrepresented in differentially expressed genes, highlighting an important role for newly emerged genes in tolerating environmental stressors. The results indicate that the D. pulex lineages tested in this study generally respond to copper stress using the same major pathways, but that the more resistant clone with previous copper exposure might be better able to regulate key genes. This finding highlights the important nuances in gene expression among clones, shaped by historical exposure and influencing copper tolerance.

Life cycle environmental impacts of advanced wastewater treatment techniques for removal of pharmaceuticals and personal care products (PPCPs)

Pharmaceutical and personal care products (PPCPs) are of increasing interest because of their ecotoxicological properties and environmental impacts. Wastewater treatment plants (WWTPs) are the main pathway for their release into freshwaters due to the inefficiency of conventional WWTPs in removing many of these contaminants from effluents. Therefore, different advanced effluent treatment techniques have been proposed for their treatment. However, it is not known at present how effective these treatment methods are and whether on a life cycle basis they cause other environmental impacts which may outweigh the benefits of the treatment. In an effort to provide an insight into this question, this paper considers life cycle environmental impacts of the following advanced treatment techniques aimed at reducing freshwater ecotoxicity potential of PPCPs: granular activated carbon (GAC), nanofiltration (NF), solar photo-Fenton (SPF) and ozonation. The results suggest that on average NF has the lowest impacts for 13 out of 18 categories considered. GAC is the best alternative for five impacts, including metals and water depletion, but it has the highest marine eutrophication. SPF and ozonation are the least sustainable for eight impacts, including ecotoxicity and climate change. GAC and NF are also more efficient in treating heavy metals while avoiding generation of harmful by-products during the treatment, thus being more suitable for potable reuse of wastewater. However, releasing the effluent without advanced treatment to agricultural land achieves a much higher reduction of freshwater ecotoxicity than treating it by any of the advanced treatments and releasing to the environment. Therefore, the use of advanced effluent treatment for agricultural purposes is not recommended.

Fishing for contaminants: identification of three mechanism specific transcriptome signatures using Danio rerio embryos

In ecotoxicology, transcriptomics is an effective way to detect gene expression changes in response to environmental pollutants. Such changes can be used to identify contaminants or contaminant classes and can be applied as early warning signals for pollution. To do so, it is important to distinguish contaminant-specific transcriptomic changes from genetic alterations due to general stress. Here we present a first step in the identification of contaminant class-specific transcriptome signatures. Embryos of zebrafish (Danio rerio) were exposed to three substances (methylmercury, chlorpyrifos and Aroclor 1254, each from 24 to 48 hpf exposed) representing sediment typical contaminant classes. We analyzed the altered transcriptome to detect discriminative genes significantly regulated in reaction to the three applied contaminants. By comparison of the results of the three contaminants, we identified transcriptome signatures and biologically important pathways (using Cytoscape/ClueGO software) that react significantly to the contaminant classes. This approach increases the chance of finding genes that play an important role in contaminant class-specific pathways rather than more general processes.

Applying adverse outcome pathways and species sensitivity-weighted distribution to predicted-no-effect concentration derivation and quantitative ecological risk assessment for bisphenol A and 4-nonylphenol in aquatic environments: A case study on Tianjin City, China

Adverse outcome pathways (AOPs) are a novel concept that effectively considers the toxic modes of action and guides the ecological risk assessment of chemicals. To better use toxicity data including biochemical or molecular responses and mechanistic data, we further developed a species sensitivity-weighted distribution (SSWD) method for bisphenol A and 4-nonylphenol. Their aquatic predicted-no-effect concentrations (PNECs) were derived using the log-normal statistical extrapolation method. We calculated aquatic PNECs of bisphenol A and 4-nonylphenol with values of 4.01 and 0.721 µg/L, respectively. The ecological risk of each chemical in different aquatic environments near Tianjin, China, a coastal municipality along the Bohai Sea, was characterized by hazard quotient and probabilistic risk quotient assessment techniques. Hazard quotients of 7.02 and 5.99 at 2 municipal sewage sites using all of the endpoints were observed for 4-nonylphenol, which indicated high ecological risks posed by 4-nonylphenol to aquatic organisms, especially endocrine-disrupting effects. Moreover, a high ecological risk of 4-nonylphenol was indicated based on the probabilistic risk quotient method. The present results show that combining the SSWD method and the AOP concept could better protect aquatic organisms from adverse effects such as endocrine disruption and could decrease uncertainty in ecological risk assessment. Environ Toxicol Chem 2018;37:551-562. © 2017 SETAC.

A statistical framework for applying RNA profiling to chemical hazard detection

Use of 'omics technologies in environmental science is expanding. However, application is mostly restricted to characterizing molecular steps leading from toxicant interaction with molecular receptors to apical endpoints in laboratory species. Use in environmental decision-making is limited, due to difficulty in elucidating mechanisms in sufficient detail to make quantitative outcome predictions in any single species or in extending predictions to aquatic communities. Here we introduce a mechanism-agnostic statistical approach, supplementing mechanistic investigation by allowing probabilistic outcome prediction even when understanding of molecular pathways is limited, and facilitating extrapolation from results in laboratory test species to predictions about aquatic communities. We use concepts familiar to environmental managers, supplemented with techniques employed for clinical interpretation of 'omics-based biomedical tests. We describe the framework in step-wise fashion, beginning with single test replicates of a single RNA variant, then extending to multi-gene RNA profiling, collections of test replicates, and integration of complementary data. In order to simplify the presentation, we focus on using RNA profiling for distinguishing presence versus absence of chemical hazards, but the principles discussed can be extended to other types of 'omics measurements, multi-class problems, and regression. We include a supplemental file demonstrating many of the concepts using the open source R statistical package.

Development and application of the SSD approach in scientific case studies for ecological risk assessment

Species sensitivity distributions (SSDs) are used in ecological risk assessment for extrapolation of the results of toxicity tests with single species to a toxicity threshold considered protective of ecosystem structure and functioning. The attention to and importance of the SSD approach has increased in scientific and regulatory communities since the 1990s. Discussion and criticism have been triggered on the concept of the approach as well as its technical aspects (e.g., distribution type, number of toxicity endpoints). Various questions remain unanswered, especially with regard to different endpoints, statistical methods, and protectiveness of threshold levels, for example. In the present literature review (covering the period 2002-2013), case studies are explored in which the SSD approach was applied, as well as how endpoint types, species choice, and data availability affect SSDs. How statistical methods may be used to construct reliable SSDs and whether the lower 5th percentile hazard concentrations (HC5s) from a generic SSD can be protective for a specific local community are also investigated. It is shown that estimated protective concentrations were determined by taxonomic groups rather than the statistical method used to construct the distribution. Based on comparisons between semifield and laboratory-based SSDs, the output from a laboratory SSD was protective of semifield communities in the majority of studies. Environ Toxicol Chem 2016;35:2149-2161. © 2016 SETAC.

Global gene expression profiling using heterologous DNA microarrays to analyze alterations in the transcriptome of Mus spretus mice living in a heavily polluted environment

Microarray platforms are a good approach for assessing biological responses to pollution as they enable the simultaneous analyses of changes in the expression of thousands of genes. As an omic and non-targeted methodology, this technique is open to unforeseen responses under particular environmental conditions. In this study, we successfully apply a commercial oligonucleotide microarray containing Mus musculus whole-genome probes to compare and assess the biological effects of living in a heavily polluted settlement, the Domingo Rubio stream (DRS), at the Huelva Estuary (SW Spain), on inhabitant free-living Mus spretus mice. Our microarray results show that mice living in DRS suffer dramatic changes in gene and protein expression compared with reference specimens. DRS mice showed alteration in the oxidative status of hepatocytes, with activation of both the innate and the acquired immune responses and the induction of chronic inflammation, accompanied by metabolic alterations that imply the accumulation of lipids in the liver (hepatic steatosis). The identified deregulated genes may be useful as biomarkers of environmental pollution.

Fold-change threshold screening: a robust algorithm to unmask hidden gene expression patterns in noisy aggregated transcriptome data

Transcriptomics is often used to investigate changes in an organism's genetic response to environmental contamination. Data noise can mask the effects of contaminants making it difficult to detect responding genes. Because the number of genes which are found differentially expressed in transcriptome data is often very large, algorithms are needed to reduce the number down to a few robust discriminative genes. We present an algorithm for aggregated analysis of transcriptome data which uses multiple fold-change thresholds (threshold screening) and p values from Bayesian generalized linear model in order to assess the robustness of a gene as a potential indicator for the treatments tested. The algorithm provides a robustness indicator (ROBI) as well as a significance profile, which can be used to assess the statistical significance of a given gene for different fold-change thresholds. Using ROBI, eight discriminative genes were identified from an exemplary dataset (Danio rerio FET treated with chlorpyrifos, methylmercury, and PCB) which could be potential indicators for a given substance. Significance profiles uncovered genetic effects and revealed appropriate fold-change thresholds for single genes or gene clusters. Fold-change threshold screening is a powerful tool for dimensionality reduction and feature selection in transcriptome data, as it effectively reduces the number of detected genes suitable for environmental monitoring. In addition, it is able to unmask patterns in altered genetic expression hidden by data noise and reduces the chance of type II errors, e.g., in environmental screening.

Applications of biological tools or biomarkers in aquatic biota: A case study of the Tamar estuary, South West England

Biological systems are the ultimate recipients of pollutant-induced damage. Consequently, our traditional reliance on analytical tools is not enough to assess ecosystem health. Biological responses or biomarkers are therefore also considered to be important tools for environmental hazard and risk assessments. Due to historical mining, other anthropogenic activities, and its conservational importance (e.g. NATURA sites, SACs), the Tamar estuary in South West England is an ideal environment in which to examine applications of such biological tools. This review presents a thorough and critical evaluation of the different biological tools used in the Tamar estuary thus far, while also discussing future perspectives for biomarker studies from a global perspective. In particular, we focus on the challenges which hinder applications of biological tools from being more readily incorporated into regulatory frameworks, with the aim of enabling both policymakers and primary stakeholders to maximise the environmental relevance and regulatory usefulness of such tools.

Species sensitivity distributions for engineered nanomaterials

Engineered nanomaterials (ENMs) are a relatively new strain of materials for which little is understood about their impacts. A species sensitivity distribution (SSDs) is a cumulative probability distribution of a chemical's toxicity measurements obtained from single-species bioassays of various species that can be used to estimate the ecotoxicological impacts of a chemical. The recent increase in the availability of acute toxicity data for ENMs enabled the construction of 10 ENM-specific SSDs, with which we analyzed (1) the range of toxic concentrations, (2) whether ENMs cause greater hazard to an ecosystem than the ionic or bulk form, and (3) the key parameters that affect variability in toxicity. The resulting estimates for hazardous concentrations at which 5% of species will be harmed ranged from <1 ug/L for PVP-coated n-Ag to >3.5 mg/L for CNTs. The results indicated that size, formulation, and the presence of a coating can alter toxicity, and thereby corresponding SSDs. Few statistical differences were observed between SSDs of an ENM and its ionic counterpart. However, we did find a significant correlation between the solubility of ENMs and corresponding SSD. Uncertainty in SSD values can be reduced through greater consideration of ENM characteristics and physiochemical transformations in the environment.

Omics technologies and their applications to evaluate metal toxicity in mice M. spretus as a bioindicator

Metals are important components of living organisms since many biological functions critically depend on their interaction with some metal in the cell. However, human activities have increased toxic metal levels in the terrestrial and aquatic ecosystems affecting living organisms. The impact of metals on cellular metabolism and global homeostasis has been traditionally assessed in free-living organisms by using conventional biomarkers; however, to obtain a global vision of metal toxicity mechanisms and the responses that metals elicit in the organisms, new analytical methodologies are needed. We review the use of omics approaches to assess the response of living organisms under metal stress illustrating the possibilities of different methodologies on the basis of our previous results. Most of this research has been based on free-living mice Mus spretus, a conventional bioindicator used to monitor metal pollution in Doñana National Park (DNP) (SW Spain), which is an important European biological reserve for migrating birds affected by agricultural, mining and industrial activities. The benefits of using omic techniques such as heterologous microarrays, proteomics methodologies (2-DE, iTRAQ®), metallomics, ionomics or metabolomics has been remarked; however, the complexity of these areas requires the integration of omics to achieve a comprehensive assessment of their environmental status. This article is part of a Special Issue entitled: Environmental and structural proteomics.

BIOLOGICAL SIGNIFICANCE

This work presents new contributions in the study of environmental metal pollution in terrestrial ecosystems using Mus spretus mice as bioindicator in Doñana National Park (SW Spain) and surroundings. In addition, it has been demonstrated that the integration of omics multi-analytical approaches provides a very suitable approach for the study of the biological response and metal interactions in exposed and free-living mice (Mus musculus and Mus spretus, respectively) under metal pollution.

Heterologous microarray analysis of transcriptome alterations in Mus spretus mice living in an industrial settlement

This work demonstrates the successful application of a commercial oligonucleotide microarray containing Mus musculus whole-genome probes to assess the biological effects of an industrial settlement on inhabitant Mus spretus mice. The transcriptomes of animals in the industrial settlement contrasted with those of specimens collected from a nearby protected ecosystem. Proteins encoded by the differentially expressed genes were broadly categorized into six main functional classes. Immune-associated genes were mostly induced and related to innate and acquired immunity and inflammation. Genes sorted into the stress-response category were mainly related to oxidative-stress tolerance and biotransformation. Metabolism-associated genes were mostly repressed and related to lipid metabolic pathways; these included genes that encoded 11 of the 20 cholesterol biosynthetic pathway enzymes. Crosstalk between members of different functional categories was also revealed, including the repression of serine-protease genes and the induction of protease-inhibitor genes to control the inflammatory response. Absolute quantification of selected transcripts was performed via RT-PCR to verify the microarray results and assess interindividual variability. Microarray data were further validated by immunoblotting and by cholesterol and protein-thiol oxidation level determinations. Reported data provide a broad impression of the biological consequences of residing in an industrial area.

Transcriptomic profiling permits the identification of pollutant sources and effects in ambient water samples

Contaminant exposure is one possible contributor to population declines of endangered fish species in the Sacramento-San Joaquin Estuary, California, including the endangered delta smelt (Hypomesus transpacificus). Herein we investigated transcriptional responses in larval delta smelt resulting from exposure to water samples collected at the Department of Water Resources Field Station at Hood, a site of concern, situated upstream of known delta smelt habitat and spawning sites and downstream of the Sacramento Regional Wastewater Treatment Plant (SRWTP). Microarray assessments indicate impacts on energy metabolism, DNA repair mechanisms and RNA processing, the immune system, development and muscle function. Transcription responses of fish exposed to water samples from Hood were compared with exposures to 9% effluent samples from SRWTP, water from the Sacramento River at Garcia Bend (SRGB), upstream of the effluent discharge, and SRGB water spiked with 2mg/L total ammonium (9% effluent equivalent). Results indicate that transcriptomic profiles from Hood are similar to 9% SRWTP effluent and ammonium spiked SRGB water, but significantly different from SRGB. SRGB samples however were also significantly different from laboratory controls, suggesting that SRWTP effluent is not solely responsible for the responses determined at Hood, that ammonium exposure likely enhances the effect of multiple-contaminant exposures, and that the observed mortality at Hood is due to the combination of both effluent discharge and contaminants arising from upstream of the tested sites.

Development of aquatic life criteria in China: viewpoint on the challenge

While more developed countries have a well-established systems to develop water quality criteria (WQC), little research has been done on the adequacy of the current WQC to protect endemic species of China. In order to maintain the health of aquatic ecosystems in China, a series of projects to establish national WQC based on regional characteristics has recently been initiated. However, the establishment of a completely novel methodology would be costly and time consuming. Also, due to the similarities in physiologies and natural histories of classes of aquatic organisms, there is no reason to believe that WQC would not be sufficient to protect unique species in China. This review was undertaken to identify key outstanding issues regarding establishment of aquatic life criteria (ALC) to be applied in China, including prioritization of chemicals, test species, mode of action, field/semi-field data, and methods of aggregating the information and calculating the ALC. This was used to identify the principle issues that need to be addressed in order to better understand the methods for development of criteria for the protection of aquatic life and provide a reference to China and other developing countries committed to the establishment of their own WQC system.

Comparison of earthworm responses to petroleum hydrocarbon exposure in aged field contaminated soil using traditional ecotoxicity endpoints and 1H NMR-based metabolomics

(1)H NMR metabolomics and conventional ecotoxicity endpoints were used to examine the response of earthworms exposed to petroleum hydrocarbons (PHCs) in soil samples collected from a site that was contaminated with crude oil from a pipeline failure in the mid-1990s. The conventional ecotoxicity tests showed that the soils were not acutely toxic to earthworms (average survival ≥ 90%), but some soil samples impaired reproduction endpoints by >50% compared to the field control soil. Additionally, metabolomics revealed significant relationships between earthworm metabolic profiles (collected after 2 or 14 days of exposure) and soil properties including soil PHC concentration. Further comparisons by partial least squares regression revealed a significant relationship between the earthworm metabolomic data (collected after only 2 or 14 days) and the reproduction endpoints (measured after 63 days). Therefore, metabolomic responses measured after short exposure periods may be predictive of chronic, ecologically relevant toxicity endpoints for earthworms exposed to soil contaminants.

Gene expression responses for detecting sublethal effects of xenobiotics and whole effluents on a Xenopus laevis embryo assay

In the present study, the authors investigated the effects of bisphenol A, chlorpyrifos, methylparaben, and 2 effluent samples from wastewater treatment plants located in the province of Madrid, Spain, on the messenger RNA expression of specific genes involved in early development (ESR1, pax6, bmp4, and myf5) and a gene involved in the general stress response (hsp70) during Xenopus laevis embryo development. Gene expression was analyzed after 4 h, 24 h, and 96 h of exposure by semiquantitative reverse-transcriptase-polymerase chain reaction. Concentration ranges of the compounds and dilutions for the samples were selected to cause morphological alterations in embryos after 96 h of exposure. Transcript levels of ESR1, pax6, and hsp70 were differentially altered at early developmental stages with patterns specific to the contaminant and the exposure time. However, further studies are needed to establish transcript levels of specific genes as biomarkers of sublethal effects in an environmental risk-assessment framework. Besides, studies including more generic responses, such as genes encoding antioxidant enzymes, together with genes related to embryonic development have to be developed to look for a battery of mechanistic endpoints for the evaluation of chemical exposure at the molecular level in a first-tier assessment.

Proteomic and metabolomic responses of clam Ruditapes philippinarum to arsenic exposure under different salinities

Arsenic (As) contamination is a severe problem in the intertidal zones of the Bohai Sea (China) with wide salinity variation. In the present study, we combined proteomics and metabolomics to characterize the differential responses of arsenic in clam Ruditapes philippinarum under different salinities (31.1, 23.3 and 15.6 psu). Both proteomic and metabolomic responses indicated that varying salinities could significantly affect the toxicological responses of clams to As. Metabolic biomarkers revealed that the environmentally relevant arsenic (20 μg L(-1)) exposure induced disturbance in energy metabolism and/or osmotic regulation under different salinities, whereas protein biomarkers indicated oxidative stress, cellular injury and apoptosis and disturbance in energy metabolism. In addition, the up-regulated proteins including ATP synthase, succinyl-CoA synthetase and nucleoside diphosphate kinase were validated by related metabolites, succinate and ATP, which confirmed the disturbance in energy metabolism in clam gills at low salinity (15.6 psu). These findings provide important insights into toxicological effects of environmental contaminant at molecular levels using combined proteomics and metabolomics.

Response of Lumbriculus variegatus transcriptome and metabolites to model chemical contaminants

Assessment of the underlying molecular events leading to xenobiotic toxicity is challenging especially when techniques are applied in isolation. We examined transcriptional and metabolic changes in Lumbriculus variegatus exposed to benzo(a)pyrene (B(a)P), cadmium (Cd) or pentachlorophenol (PCP) by DNA microarrays (7422 ESTs) and gas chromatography-mass spectrometry (GC-MS), respectively. In addition, the DNA damage response of worms exposed to B(a)P was assessed by a capillary electrophoresis laser induced fluorescence (CE-LIF) immunoassay. We found elevated expression of oxidative stress responsive genes, which correlated positively with the changes in antioxidant vitamin precursors including alpha-tocopherol and cholecalciferol. Other genes with strong differential expressions were mostly involved in actin related processes and proteolysis, despite an apparent delayed Cd response. Phosphates, sugars and fatty acids were effectively reduced and suggested that chemical treatments may have interfered with energy metabolism. The increased amount of B(a)P diol-epoxide (BPDE)-DNA adducts in exposed worms appeared to correlate with the variability in uridine, inosine and xanthine, which are key components of nucleoside metabolism. This suggests that DNA damage was imminent or peaked within 6h. The results conformed to transcriptional changes in B(a)P exposed worms and compliment other approaches to elucidate underlying molecular changes.

Cytochrome P4501A mRNA and protein induction in striped bass (Morone saxatilis)

The striped bass (Morone saxatilis) supports a valuable recreational fishery and is among the most important piscivorous fish of the San Francisco Estuary. This species has suffered a significant decline in numbers over the past decades, and there is indication that contaminants are important contributors. Polycyclic aromatic hydrocarbons (PAHs) and polyhalogenated aromatic hydrocarbons (PHAHs) including PCBs and dioxins are widespread in the estuary, they typically bioaccumulate through trophic levels, reaching highest levels in top predators and are known to affect the fish health and development. The aim of this study was to investigate the dynamics of cytochrome P4501A (Cyp1a) induction simultaneously at different levels of biological organization (RNA transcription and protein synthesis) as a biomarker of exposure to PAHs and PHAHs. We utilized β-naphthoflavone (BNF) as a model PAH to induce Cyp1a responses in juvenile striped bass in both dose-response and time-response assessments and determined Cyp1a mRNA and protein levels. Significant responses were measured in both systems at 10 mg ΒΝF kg⁻¹, a concentration used for time-response studies. Messenger RNA levels peaked at 6 h post-injection, while protein levels increased progressively with time, significantly peaking at 96 h post-injection; both remaining elevated throughout the duration of the test (8 days). Our data suggest that rapid induction of gene transcription following exposure and subsequent cumulative protein synthesis could provide a useful means of identifying temporal variants in exposure to Cyp1a inducers in Morone saxatilis. The potential application of this combined Cyp1a gene and protein biomarker in this species for field studies is discussed.

Transcription profiling in environmental diagnostics: health assessments in Columbia River basin steelhead (Oncorhynchus mykiss)

The health condition of out-migrating juvenile salmonids can influence migration success. Physical damage, pathogenic infection, contaminant exposure, and immune system status can affect survival probability. The present study is part of a wider investigation of out-migration success in juvenile steelhead (Oncorhynchus mykiss) and focuses on the application of molecular profiling to assess sublethal effects of environmental stressors in field-collected fish. We used a suite of genes in O. mykiss to specifically assess responses that could be directly related to steelhead health condition during out-migration. These biomarkers were used on juvenile steelhead captured in the Snake River, a tributary of the Columbia River, in Washington, USA, and were applied on gill and anterior head kidney tissue to assess immune system responses, pathogen-defense (NRAMP, Mx, CXC), general stress (HSP70), metal-binding (metallothionein-A), and xenobiotic metabolism (Cyp1a1) utilizing quantitative polymerase chain reaction (PCR) technology. Upon capture, fish were ranked according to visual external physical conditions into good, fair, poor, and bad categories; gills and kidney tissues were then dissected and preserved for gene analyses. Transcription responses were tissue-specific for gill and anterior head kidney with less significant responses in gill tissue than in kidney. Significant differences between the condition ranks were attributed to NRAMP, MX, CXC, and Cyp1a1 responses. Gene profiling correlated gene expression with pathogen presence, and results indicated that gene profiling can be a useful tool for identifying specific pathogen types responsible for disease. Principal component analysis (PCA) further correlated these responses with specific health condition categories, strongly differentiating good, poor, and bad condition ranks. We conclude that molecular profiling is an informative and useful tool that could be applied to indicate and monitor numerous population-level parameters of management interest.

Soil ecotoxicology: state of the art and future directions

Developments in soil ecotoxicology started with observations on pesticide effects on soil invertebrates in the 1960s. To support the risk assessment of chemicals, in the 1980s and 1990s development of toxicity tests was the main issue, including single species tests and also more realistic test systems like model ecosystems and field tests focusing on structural and functional endpoints. In the mean time, awareness grew about issues like bioavailability and routes of exposure, while biochemical endpoints (biomarkers) were proposed as sensitive and potential early-warning tools. In recent years, interactions between different chemicals (mixture toxicity) and between chemical and other stressors attracted scientific interest. With the development of molecular biology, omics tools are gaining increasing interest, while the ecological relevance of exposure and effects is translating into concepts like (chemical) stress ecology, ecological vulnerability and trait-based approaches. This contribution addresses historical developments and focuses on current issues in soil ecotoxicology. It is concluded that soil ecotoxicological risk assessment would benefit from extending the available battery of toxicity tests by including e.g. isopods, by paying more attention to exposure, bioavailability and toxicokinetics, and by developing more insight into the ecology of soil organisms to support better understanding of exposure and long-term consequences of chemical exposure at the individual, population and community level. Ecotoxicogenomics tools may also be helpful in this, but will require considerable further research before they can be applied in the practice of soil ecotoxicological risk assessment.

Modeling physiological processes that relate toxicant exposure and bacterial population dynamics

Quantifying effects of toxicant exposure on metabolic processes is crucial to predicting microbial growth patterns in different environments. Mechanistic models, such as those based on Dynamic Energy Budget (DEB) theory, can link physiological processes to microbial growth.Here we expand the DEB framework to include explicit consideration of the role of reactive oxygen species (ROS). Extensions considered are: (i) additional terms in the equation for the "hazard rate" that quantifies mortality risk; (ii) a variable representing environmental degradation; (iii) a mechanistic description of toxic effects linked to increase in ROS production and aging acceleration, and to non-competitive inhibition of transport channels; (iv) a new representation of the "lag time" based on energy required for acclimation. We estimate model parameters using calibrated Pseudomonas aeruginosa optical density growth data for seven levels of cadmium exposure. The model reproduces growth patterns for all treatments with a single common parameter set, and bacterial growth for treatments of up to 150 mg(Cd)/L can be predicted reasonably well using parameters estimated from cadmium treatments of 20 mg(Cd)/L and lower. Our approach is an important step towards connecting levels of biological organization in ecotoxicology. The presented model reveals possible connections between processes that are not obvious from purely empirical considerations, enables validation and hypothesis testing by creating testable predictions, and identifies research required to further develop the theory.

Concentration-response concept in ecotoxicoproteomics: effects of different phenanthrene concentrations to the zebrafish (Danio rerio) embryo proteome

Concentration-response experiments, based on the testing of less replicates in favour of more exposure concentrations, represent the typical design of choice applied in toxicological and ecotoxicological effect assessment studies using traditional endpoints such as lethality. However, to our knowledge this concept has not found implementation in the increasingly applied OMICS techniques studying thousands of molecular endpoints at the same time. The present study is among the first applying the concentration-response concept for an ecotoxicoproteomics study. The effects of six different concentrations in the low effect range (<LC₂₀) of the PAH phenanthrene to the proteome of the ecotoxicological vertebrate model zebrafish (Danio rerio) embryo were investigated (two replicates per concentration) after 5 days exposure. Proteomics analyses were performed on organism extracts using 2-DE DIGE. Protein abundance profiles of around 713 protein spots were studied. About one-third of the protein signals could be detected to show robust reactions correlating with stressor concentration. Within this group, 65 protein signals showed significant changes compared to controls already at 1% lethal concentration (LC₀₁). Interestingly, 28 proteins significantly reacted at very low concentrations (<LC₀₁) and showed an exposure concentration dependent regulation status. Characteristic protein spots were identified by mass spectrometry. With the results of the present study the utility and several benefits using a concentration-response approach in proteomics studies could be shown. These included (i) knowledge about and the ability to model concentration dependent dynamics of molecular endpoints, (ii) to gain information about sensitivity of the molecular response in comparison to traditional endpoints and (iii) to help selecting the most promising protein spots for further investigations such as protein identification and biomarker studies. Using this experimental design based on testing of several exposure concentrations and less replicates might provide a step forward in getting increased output from toxicoproteomics studies.

Changes in gene transcription and whole organism responses in larval fathead minnow (Pimephales promelas) following short-term exposure to the synthetic pyrethroid bifenthrin

The combination of molecular and whole-organism endpoints in ecotoxicology provides valuable information about the ecological relevance of sublethal stressor effects in aquatic ecosystems such as those caused by the use of insecticides and translocation of their residues into surface waters. This study contributes knowledge about the sublethal effects of a common use insecticide, the synthetic pyrethroid bifenthrin, on larval fathead minnow (Pimephales promelas). Transcriptomic responses, assessed by quantitative real-time PCR, combined with individual effects on swimming performance were used to estimate the ecological relevance of insecticide impacts. Significant transcriptomic responses were observed at 0.07 μg L(-1) bifenthrin (lowest observed effect concentration, LOEC) but mostly followed a biphasic rather than a linear dose-response with increasing concentration. Transcript patterns for genes involved in detoxification, neuromuscular function and energy metabolism were linked to an impairment of swimming performance at ≥0.14 μg L(-1) bifenthrin. With increasing treatment concentration, a significant down-regulation was observed for genes coding for cyp3a, aspartoacylase, and creatine kinase, whereas metallothionein was up-regulated. Additionally, bifenthrin induced endocrine responses as evident from a significant up-regulation of vitellogenin and down-regulation of insuline-like growth factor transcripts. Recovery occurred after 6 days and was dependent on the magnitude of the initial stress. During the recovery period, down-regulation of vitellogenin was observed at lowest exposure concentrations. The data presented here emphasize that links can be made between gene transcription changes and behavioral responses which is of great value for the evaluation and interpretation of biomarker responses.

An integrated multi-disciplinary approach for studying multiple stressors in freshwater ecosystems: Daphnia as a model organism

The increased overexploitation of freshwater ecosystems and their extended watersheds often generates a cascade of anthropogenic stressors (e.g., acidification, eutrophication, metal contamination, Ca decline, changes in the physical environment, introduction of invasive species, over-harvesting of resources). The combined effect of these stressors is particularly difficult to study, requiring a coordinated multi-disciplinary effort and insights from various sub-disciplines of biology, including ecology, evolution, toxicology, and genetics. It also would benefit from a well-developed and broadly accepted model systems. The freshwater crustacean Daphnia is an excellent model organism for studying multiple stressors because it has been a chosen focus of study in all four of these fields. Daphnia is a widespread keystone species in most freshwater ecosystems, where it is routinely exposed to a multitude of anthropogenic and natural stressors. It has a fully sequenced genome, a well-understood life history and ecology, and a huge library of responses to toxicity. To make the case for its value as a model species, we consider the joint and separate effects of natural and three anthropogenic stressors-climatic change, calcium decline, and metal contaminants on daphniids. We propose that integrative approaches marrying various subfields of biology can advance our understanding of the combined effects of stressors. Such approaches can involve the measuring of multiple responses at several levels of biological organization from molecules to natural populations. For example, novel interdisciplinary approaches such as transcriptome profiling and mutation accumulation experiments can offer insights into how multiple stressors influence gene transcription and mutation rates across genomes, and, thus, help determine the causal mechanism between environmental stressors and population/community effects as well as long-term evolutionary patterns.

Omic approaches in environmental issues

Biomonitoring requires the application of batteries of different biomarkers, as environmental contaminants induce multiple responses in organisms that are not necessarily correlated. Omic technologies were proposed as an alternative to conventional biomarkers since these techniques quantitatively monitor many biological molecules in a high-throughput manner and thus provide a general appraisal of biological responses altered by exposure to contaminants. As the studies using omic technologies increase, it is becoming clear that any single omic approach may not be sufficient to characterize the complexity of ecosystems. This work aims to provide a preliminary working scheme for the use of combined transcriptomic and proteomic methodologies in environmental biomonitoring. There are difficulties in working with nonmodel organisms as bioindicators when combining several omic approaches. As a whole, our results with heterologous microarrays in M. spretus and suppressive subtractive hybridization (SSH) in P. clarkii indicated that animals sustaining a heavy pollution burden exhibited an enhanced immune response and/or cell apoptosis. The proteomic studies, although preliminary, provide a holistic insight regarding the manner by which pollution shifts protein intensity in two-dimensional gel electrophoresis (2-DE), completing the transcriptomic approach. In our study, the sediment element concentration was in agreement with the intensity of protein expression changes in C. maenas crabs. In conclusion, omics are useful technologies in addressing environmental issues and the determination of contamination threats.