Ecotoxicogenomics: the challenge of integrating genomics into aquatic and terrestrial ecotoxicology

Collembase: a repository for springtail genomics and soil quality assessment

BACKGROUND

Environmental quality assessment is traditionally based on responses of reproduction and survival of indicator organisms. For soil assessment the springtail Folsomia candida (Collembola) is an accepted standard test organism. We argue that environmental quality assessment using gene expression profiles of indicator organisms exposed to test substrates is more sensitive, more toxicant specific and significantly faster than current risk assessment methods. To apply this species as a genomic model for soil quality testing we conducted an EST sequencing project and developed an online database.

DESCRIPTION

Collembase is a web-accessible database comprising springtail (F. candida) genomic data. Presently, the database contains information on 8686 ESTs that are assembled into 5952 unique gene objects. Of those gene objects approximately 40% showed homology to other protein sequences available in GenBank (blastx analysis; non-redundant (nr) database; expect-value < 10-5). Software was applied to infer protein sequences. The putative peptides, which had an average length of 115 amino-acids (ranging between 23 and 440) were annotated with Gene Ontology (GO) terms. In total 1025 peptides (approximately 17% of the gene objects) were assigned at least one GO term (expect-value < 10-25). Within Collembase searches can be conducted based on BLAST and GO annotation, cluster name or using a BLAST server. The system furthermore enables easy sequence retrieval for functional genomic and Quantitative-PCR experiments. Sequences are submitted to GenBank (Accession numbers: EV473060 - EV481745).

CONCLUSION

Collembase http://www.collembase.org is a resource of sequence data on the springtail F. candida. The information within the database will be linked to a custom made microarray, based on the Agilent platform, which can be applied for soil quality testing. In addition, Collembase supplies information that is valuable for related scientific disciplines such as molecular ecology, ecogenomics, molecular evolution and phylogenetics.

Proteomic analyses indicate induction of hepatic carbonyl reductase/20beta-hydroxysteroid dehydrogenase B in rainbow trout exposed to sewage effluent

Proteomic analyses were performed to identify regulated liver proteins in rainbow trout (Oncorhynchus mykiss) caged upstream and downstream from a sewage treatment works (STW). Two-dimensional gel electrophoresis, image analysis and FT-ICR mass-spectrometry revealed four regulated protein spots. The three down-regulated spots contained betaine aldehyde dehydrogenase, lactate dehydrogenase and an unidentified protein respectively. The only up-regulated spot consisted of both mitochondrial ATP synthase alpha-subunit and carbonyl reductase/20beta-hydroxysteroid dehydrogenase (CR/20beta-HSD). Further studies using quantitative PCR revealed a 13.5-fold induction of CR/20beta-HSD B mRNA following STW effluent exposure. The CR/20beta-HSD B gene was not regulated by 17alpha-ethinylestradiol, suggesting that its induction downstream from the STW is due to other factors than exposure to estrogens. Image analysis was initially performed on four gels from each group. These analyses suggested 15 regulated spots. However, validation of the 15 spots by increasing the number of replicates confirmed only four regulated spots. Hence, the present study also demonstrates the need for sufficient biological/technical replication in the interpretation of proteomic data.

A simple and rapid method for preserving RNA of aquatic invertebrates for ecotoxicogenomics

Here we describe a novel, inexpensive and simple method for preserving RNA that reduces handling stress in aquatic invertebrates following ecotoxicogenomic experimentation. The application of the method is based on transcriptomic experiments conducted on Daphnia magna, but may easily be applied on a range of other aquatic organisms of a particular size with e.g. amphipod Gammarus pulex representing an upper size limit. We explain in detail how to apply this new method, named the "Cylindrical Sieve (CS) system", and highlight its advantages and disadvantages.

Toxicogenomics and ecotoxicogenomics for studying endocrine disruption and basic biology

Chemicals released into the environment have the potential to disrupt the endocrine system in wild animals, mouse, and humans. To understand molecular mechanisms of chemical toxicity in various species, toxicogenomics/ecotoxicogenomics, describing the integration of genomics (trascriptomics, proteomics and metabolomics) into toxicology/ecotoxicology, needs to be established as a powerful tool for research. Ecotoxicogenomics is defined as the study of gene and protein expression in non-target organisms that is important in responses to environmental toxicant exposures. Estrogen-responsive genes and estrogen response element(s) in genes have been identified in the mouse reproductive tract by application of cDNA microarray technology. Additionally, functional mechanisms of tributyltin action via nuclear receptors such as retinoid X receptor alpha and peroxisome proliferator activated receptor gamma also have been identified using cDNA microarray. A microarray system has been established for Daphnia magna. Toxicogenomics/ecotoxicogenomics provide powerful tools to help us understand not only molecular mechanisms of chemical toxicity but also the basic biology of various animal species.

Estimating the efficiency of fish cross-species cDNA microarray hybridization

Using an available cross-species cDNA microarray is advantageous for examining multigene expression patterns in non-model organisms, saving the need for construction of species-specific arrays. The aim of the present study was to estimate relative efficiency of cross-species hybridizations across bony fishes, using bioinformatics tools. The methodology may serve also as a model for similar evaluations in other taxa. The theoretical evaluation was done by substituting comparative whole-transcriptome sequence similarity information into the thermodynamic hybridization equation. Complementary DNA sequence assemblages of nine fish species belonging to common families or suborders and distributed across the bony fish taxonomic branch were selected for transcriptome-wise comparisons. Actual cross-species hybridizations among fish of different taxonomic distances were used to validate and eventually to calibrate the theoretically computed relative efficiencies.

Metabolic profile biomarkers of metal contamination in a sentinel terrestrial species are applicable across multiple sites

In this study, we addressed the question of whether an omic approach could genuinely be useful for biomarker profile analysis across different field sites with different physicochemical characteristics. We collected earthworms (Lumbricus rubellus) from seven sites with very different levels of metal contamination and prevailing soil type and analyzed tissue extracts by 1H nuclear magnetic resonance spectroscopy. Pattern recognition analysis of the data showed that both site- and contaminant-specific effects on the metabolic profiles could be discerned. Zinc was identified as the probable major contaminant causing a metabolic change in the earthworms. Individual sites could be resolved on the basis of NMR spectral profiles by principal component analysis; these site differences may also have been caused by additional abiotic factors such as soil pH. Despite an inevitable degree of confounding between site and contaminant concentrations, it was possible to identify metabolites which were correlated with zinc across all different sites. This study therefore acts as a proof of principle for the use of NMR-based metabolic profiling as a diagnostic tool for ecotoxicological research in polluted field soils.

Sensitive and robust gene expression changes in fish exposed to estrogen--a microarray approach

Background

Vitellogenin is a well established biomarker for estrogenic exposure in fish. However, effects on gonadal differentiation at concentrations of estrogen not sufficient to give rise to a measurable vitellogenin response suggest that more sensitive biomarkers would be useful. Induction of zona pellucida genes may be more sensitive but their specificities are not as clear. The objective of this study was to find additional sensitive and robust candidate biomarkers of estrogenic exposure.

Results

Hepatic mRNA expression profiles were characterized in juvenile rainbow trout exposed to a measured concentration of 0.87 and 10 ng ethinylestradiol/L using a salmonid cDNA microarray. The higher concentration was used to guide the subsequent identification of generally more subtle responses at the low concentration not sufficient to induce vitellogenin. A meta-analysis was performed with data from the present study and three similar microarray studies using different fish species and platforms. Within the generated list of presumably robust responses, several well-known estrogen-regulated genes were identified. Two genes, confirmed by quantitative RT-PCR (qPCR), fulfilled both the criteria of high sensitivity and robustness; the induction of the genes encoding zona pellucida protein 3 and a nucleoside diphosphate kinase (nm23).

Conclusion

The cross-species, cross-platform meta-analysis correctly identified several robust responses. This adds confidence to our approach used for identifying candidate biomarkers. Specifically, we propose that analyses of an nm23 gene together with zona pellucida genes may increase the possibilities to detect an exposure to low levels of estrogenic compounds in fish.

Using mummichog (Fundulus heteroclitus) arrays to monitor the effectiveness of remediation at a superfund site in Charleston, South Carolina, U.S.A

We previously developed a cDNA array for mummichogs (Fundulus heteroclitus), an estuarine minnow, that is targeted for identifying differentially expressed genes from exposure to polycyclic aromatic hydrocarbons and several metals, including chromium. A chromium-contaminated Superfund site at Shipyard Creek in Charleston, South Carolina, USA, is undergoing remediation, providing us a unique opportunity to study the utility of arrays for monitoring the effectiveness of site remediation. Mummichogs were captured in Shipyard Creek in Charleston prior to remediation (2000) and after remediation began (2003 and 2005). Simultaneously, mummichogs were collected from a reference site at the Winyah Bay National Estuarine Research Reserve (NERR) in Georgetown, South Carolina, USA. The hepatic gene expression pattern of fish captured at Shipyard Creek in 2000 showed wide differences from the fish captured at NERR in 2000. Interestingly, as remediation progressed the gene expression pattern of mummichogs captured at Shipyard Creek became increasingly similar to those captured at NERR. The arrays acted as multidimensional biomarkers as the number of differentially expressed genes dropped from 22 in 2000 to four in 2003, and the magnitude of differential expression dropped from 3.2-fold in 2000 to no gene demonstrating a difference over 1.5-fold in 2003. Furthermore, the arrays indicated changes in the bioavailability of chromium caused by hydraulic dredging in the summer of 2005. This research is, to our knowledge, the first report using arrays as biomarkers for a weight-of-evidence hazard assessment and demonstrates that arrays can be used as multidimensional biomarkers to monitor site mitigation because the gene expression profile is associated with chromium bioavailability and body burden.

Analytical chemistry and ecotoxicology--tasks, needs and trends

Thousands of synthetic chemicals are continuously released into the environment, where they have the potential to produce adverse effects on ecosystems, even at low concentrations. To avoid adverse effects induced by environmental pollutants, science and society have developed methods and tools to (1) measure the fate and distribution of the pollutants (analytical chemistry), (2) analyze effects of pollutants on biota in standardized bioassays (ecotoxicology), and (3) combine the data for risk assessment procedures. Though such procedures are well established, new issues in environmental chemistry and ecotoxicology arise and continue to challenge scientists. Furthermore, analytical chemistry has undergone enormous technological progress and new techniques from the life sciences area have been added to the ecotoxicology toolbox. Here our view on the tasks and needs in environmental chemistry and ecotoxicology is presented and the current trends under development in the field are illustrated.

Microevolution and ecotoxicology of metals in invertebrates

Risk assessment of metal-contaminated habitats based on responses in the field is complicated by the evolution of local, metal-resistant ecotypes. The unpredictability of occurrence of genetically determined adaptive traits, in terms of site-specific geochemistry, a population's inferred exposure history, and in the physiology of resistance, exacerbates the problem. Micro-evolutionary events warrant the attention of ecotoxicologists because they undermine the application of the bedrock of toxicology, the dose-response curve, to in situ field assessments. Here we survey the evidence for the existence of genetically differentiated, metal-resistant, invertebrate populations; we also describe some of the molecular mechanisms underpinning the adaptations. Quantitative changes in tissue-metal partitioning, and in the molecular and cellular responses (biomarkers)to alterations in internal bioreactive metal pools, are widely accepted as indicators of toxicity and/or exposure in free-living organisms. Both can be modulated by resistance. The understanding that all genomes are intrinsicallyflexible, with subtle sequence changes in promoter regions or epigenetic adjustments conferring significant phenotypic consequences, is deemed highly relevant. Equally relevant is the systems biology insight that genes and proteins are woven into networks. We advocate that biomarker studies should work toward assimilating and exploiting these biological realities through monitoring the activities of suites of genes (transcriptomics) and their expressed products (proteomics), as well as profiling the metabolite signatures of individuals (metabolomics) and by using neutral genetic markers to genotype populations. Ecotoxicology requires robust tools that recognize the imprint of evolution on the constitution of field populations, as well as sufficient mechanistic understanding of the molecular-genetic observations to interpret them in meaningful environmental diagnostic ways.

Daphnia magna and ecotoxicogenomics: gene expression profiles of the anti-ecdysteroidal fungicide fenarimol using energy-, molting- and life stage-related cDNA libraries

In the present study, the existing life stage-specific cDNA library was extended with energy- and molting-related genes using Suppression Subtractive Hybridization PCR and a microarray for the aquatic test organism Daphnia magna was created. A gene set of 2455 fragments was produced belonging to different pathways such as carbohydrate and lipid metabolism, O2 transport and heme metabolism, immune response, embryo development, cuticula metabolism and visual perception pathways. Using this custom microarray, gene expression profiles were generated from neonates exposed to three concentrations of the anti-ecdysteroidal fungicide fenarimol (0.5, 0.75, 1 microg/ml) during 48 h and 96 h. In total, 59 non-redundant genes were differentially expressed, of which more genes were down- than up-regulated. The gene expression data indicated a main effect on molting specific pathways. At the highest concentration, a set of proteolytic enzymes - including different serine proteases and carboxypeptidases - were induced whereas different cuticula proteins were down-regulated (48 h). Moreover, effects on embryo development were demonstrated at the gene expression as well as at the organismal level. The embryo development related gene vitellogenin was differentially expressed after 96 h of exposure together with a significant increase in embryo abnormalities in the offspring. This study suggests that this Daphnia magna microarray is of great further value for the elucidation of molecular mechanisms of toxicity and for the future development of specific biomarkers for hazard characterization.

Mysid crustaceans as standard models for the screening and testing of endocrine-disrupting chemicals

Investigative efforts into the potential endocrine-disrupting effects of chemicals have mainly concentrated on vertebrates, with significantly less attention paid to understanding potential endocrine disruption in the invertebrates. Given that invertebrates account for at least 95% of all known animal species and are critical to ecosystem structure and function, it remains essential to close this gap in knowledge and research. The lack of progress regarding endocrine disruption in invertebrates is largely due to: (1) our ignorance of mode-of-action, physiological control, and hormone structure and function in invertebrates; (2) lack of a standardized invertebrate assay; (3) the irrelevance to most invertebrates of the proposed activity-based biological indicators for endocrine disruptor (ED) exposure (androgen, estrogen, and thyroid); (4) limited field studies. Past and ongoing research efforts using the standard invertebrate toxicity test model, the mysid shrimp, have aimed at addressing some of these issues. The present review serves as an update to a previous publication on the use of mysids for the evaluation of EDs (Verslycke et al. 2004a). It summarizes recent investigative efforts that have significantly advanced our understanding of invertebrate-specific endocrine toxicity, population modeling, field studies, and transgeneration standard test development using the mysid model.

Temporal changes in gene expression in rainbow trout exposed to ethynyl estradiol

We examined changes in the genomic response during continuous exposure to the xenoestrogen ethynyl estradiol. Isogenic rainbow trout Oncorhynchus mykiss were exposed to nominal concentrations of 100 ng/L ethynyl estradiol (EE2) for a period of 3 weeks. At fixed time points within the exposure, fish were euthanized, livers harvested and RNA extracted. Fluorescently labeled cDNA were generated and hybridized against a commercially available Salmonid array (GRASP project, University of Victoria, Canada) spotted with 16,000 cDNAs. The slides were scanned to measure abundance of a given transcript in each sample relative to controls. Data were analyzed via Genespring (Silicon Genetics) to identify a list of up and down regulated genes, and to determine gene clustering patterns that can be used as "expression signatures". Gene ontology was determined using the annotation available from the GRASP website. Our analysis indicates each exposure time period generated specific gene expression profiles. Changes in gene expression were best understood by grouping genes by their gene expression profiles rather than examining fold change at a particular time point. Many of the genes commonly used as biomarkers of exposure to xenoestrogens were not induced initially and did not have gene expression profiles typical of the majority of genes with altered expression.

A graphical systems model to facilitate hypothesis-driven ecotoxicogenomics research on the teleost brain-pituitary-gonadal axis

Graphical systems models are powerful tools that can help facilitate hypothesis-driven ecotoxicogenomic research and aid in mechanistic interpretation of results. This paper describes a novel graphical model of the teleost brain-pituitary-gonadal (BPG) axis designed for ecotoxicogenomics research on endocrine-disrupting chemicals using small fish models. The model incorporates six compartments representing the major organs involved in the fish reproductive axis and depicts the interactions of over 105 proteins and 40 simple molecules, transcriptional regulation of 25 genes, and over 300 different reactions/ processes. Application of the model is illustrated in the context of a study examining effects of the competitive aromatase inhibitor, fadrozole, on gene expression in gonad, brain, and liver tissue of fathead minnows. Changes in mRNA transcript abundance were measured using a fathead minnow oligonucleotide microarray and quantitative real-time polymerase chain reaction. Gene expression changes observed in the ovaries of females exposed to 6.3 microg fadrozole/L for7 d were functionally consistent with fadrozole's mechanism of action, and expected compensatory responses of the BPG axis to fadrozole's effects. Furthermore, microarray results helped identify additional elements (genes/ proteins) that could be included in the model to potentially increase its predictive capacity. With proper recognition of

Development of the GENIPOL European flounder (Platichthys flesus) microarray and determination of temporal transcriptional responses to cadmium at low dose

We have constructed a high density, 13 270-clone cDNA array for the sentinel fish species European flounder (Platichthys flesus), combining clones from suppressive subtractive hybridization and a liver cDNA library; DNA sequences of 5211 clones were determined. Fish were treated by single intraperitoneal injection with 50 micrograms cadmium chloride per kilogram body weight, a dose relevant to environmental exposures, and hepatic gene expression changes were determined at 1, 2, 4, 8, and 16 days postinjection in comparison to saline-treated controls. Gene expression responses were confirmed by real-time reverse transcription polymerase chain reaction (RT-PCR). Blast2GO gene ontology analysis highlighted a general induction of the unfolded protein response, response to oxidative stress, protein synthesis, transport, and degradation pathways, while apoptosis, cell cycle, cytoskeleton, and cytokine genes were also affected. Transcript levels of cytochrome P450 1A (CYP1A) were repressed and vitellogenin altered, real-time PCR showed induction of metallothionein. We thus describe the establishment of a useful resource for ecotoxicogenomics and the determination of the temporal molecular responses to cadmium, a prototypical heavy metal pollutant.

The pitfalls of proteomics experiments without the correct use of bioinformatics tools

The elucidation of the entire genomic sequence of various organisms, from viruses to complex metazoans, most recently man, is undoubtedly the greatest triumph of molecular biology since the discovery of the DNA double helix. Over the past two decades, the focus of molecular biology has gradually moved from genomes to proteomes, the intention being to discover the functions of the genes themselves. The postgenomic era stimulated the development of new techniques (e.g. 2-DE and MS) and bioinformatics tools to identify the functions, reactions, interactions and location of the gene products in tissues and/or cells of living organisms. Both 2-DE and MS have been very successfully employed to identify proteins involved in biological phenomena (e.g. immunity, cancer, host-parasite interactions, etc.), although recently, several papers have emphasised the pitfalls of 2-DE experiments, especially in relation to experimental design, poor statistical treatment and the high rate of 'false positive' results with regard to protein identification. In the light of these perceived problems, we review the advantages and misuses of bioinformatics tools - from realisation of 2-DE gels to the identification of candidate protein spots - and suggest some useful avenues to improve the quality of 2-DE experiments. In addition, we present key steps which, in our view, need to be to taken into consideration during such analyses. Lastly, we present novel biological entities named 'interactomes', and the bioinformatics tools developed to analyse the large protein-protein interaction networks they form, along with several new perspectives of the field.

Genotoxic risk identification of soil contamination at a major industrialized city in northeast China by a combination of in vitro and in vivo bioassays

The present study evaluated the genotoxicity of field soils in the Tianjin area, one of the most industrialized contaminated areas in northeast China. The genotoxicity of organic extracts of 41 soils was assayed by an in vitro SOS/ umu bioassay with Salmonella typhimurium TA 1535/pSK 1002. From the 41 soil samples, 11 samples were selected to confirm the genotoxic effect by in vivo single-cell gel electrophoresis (comet assay) using earthworms (Eisenia fetida). The results obtained demonstrated that, in the in vitro assay, genotoxicity expressed as induction ratios (IR) ranged from 1.00 to 4.60, and in the in vivo assay, the genotoxicity expressed as tail moment (TM) varied from 14.6 to 57.8 microm. All samples with high genotoxicity assessed by the SOS/umu bioassay possessed significantly high genotoxic effects in the comet assay, and there was a correlation (R2 = 0.736, p < 0.05) between IR and TM in both bioassays. It is concluded that soils in the Tianjin area were seriously contaminated by organic genotoxicants and higher levels of genotoxic effects existed in soils in the urban area of Tianjin as well as in areas near the coastal towns in the northeast part of the city. It can be concluded that a combination of in vivo and in vitro bioassays as a powerful and efficient genotoxicity-assessing tool could facilitate the assessment of genotoxic risk at a regional scale.

Hepatic gene expression profiling using Genechips in zebrafish exposed to 17alpha-ethynylestradiol

Genomic, proteomic, and metabolomic technologies continue to receive increasing interest from environmental toxicologists. This interest is due to the great potential of these technologies to identify detailed modes of action and to provide assistance in the evaluation of a contaminant's risk to aquatic organisms. Our experimental model is the zebrafish (Danio rerio) exposed to reference endocrine disrupting compounds in order to investigate compound-induced changes in gene transcript profiles. Adult, female zebrafish were exposed to 0, 15, 40, and 100ng/L of 17alpha-ethynylestradiol (EE2) and concentration and time-dependent changes in hepatic gene expression were examined using Affymetrix GeneChip Zebrafish Genome Microarrays. At 24, 48, and 168h, fish were sacrificed and liver mRNA was extracted for gene expression analysis (24 and 168h only). In an effort to link gene expression changes to effects on higher levels of biological organization, body and ovary weights were measured and blood was collected for measurement of plasma steroid hormones (17beta-estradiol (E2), testosterone (T)) and vitellogenin (VTG) using ELISA. EE2 exposure significantly affected gene expression, GSI, E2, T, and VTG. We observed 1622 genes that were significantly affected (p< or =0.001) in a concentration-dependent manner by EE2 exposure at either 24 or 168h. Gene ontology (GO) analysis revealed that EE2 exposure affected genes involved in hormone metabolism, vitamin A metabolism, steroid binding, sterol metabolism, and cell growth. Plasma VTG was significantly increased at 24, 48, and 168h (p< or =0.05) at 40 and 100ng/L and at 15ng/L at 168h. E2 and T were significantly reduced following EE2 exposure at 48 and 168h. GSI was decreased in a concentration-dependent manner at 168h. In this study, we identified genes involved in a variety of biological processes that have the potential to be used as markers of exposure to estrogenic substances. Future work will evaluate the use of these genes in zebrafish exposed to weak estrogens to determine if these genes are indicative of exposure to estrogens with varying potencies.

Gene Expression Profiles in Fathead Minnow Exposed to 2,4-DNT: Correlation with Toxicity in Mammals

Toxicogenomics, the genome-wide analysis of gene expression to study the effect of toxicants, has great potential for use in environmental toxicology. Applied to standard test organisms, it has possible applications in aquatic toxicology as a sensitive monitoring tool to detect the presence of contaminants while providing information on the mechanisms of action of these pollutants. We describe the use of a complementary DNA (cDNA) microarray of the fathead minnow (Pimephales promelas) a standard sentinel organism in aquatic toxicology, to better understand the mechanisms of toxicity of 2,4-dinitrotoluene (2,4-DNT) which is released in the environment through military and industrial use. We have constructed a fathead minnow microarray containing 5000 randomly picked anonymous cDNAs from a whole fish cDNA library. Expression profiles were analyzed in fish exposed to 2,4-DNT for 10 days at three concentrations (11, 22, and 44 microM, respectively) below the measured median lethal concentration (58 microM). Sequence analysis of cDNAs corresponding to differentially expressed genes affected by exposure revealed that lipid metabolism and oxygen transport genes were prominently affected in a dose-specific manner. We measured liver lipids and demonstrate that lipid metabolism is indeed perturbed following exposure. These observations correlate well with available toxicological data on 2,4-DNT. We present possible modes of action of 2,4-DNT toxicity and suggest that fathead minnow cDNA microarrays can be useful to identify mechanisms of toxicity in fish and as a predictive tool for toxicity in mammals.

The fathead minnow in aquatic toxicology: past, present and future

This paper reviews the roles of the fathead minnow (Pimephales promelas) as a small fish model in the field of aquatic toxicology. The species has been (and is) extensively used both for regulatory testing and research, especially in North America. For example, tests with the fathead minnow, ranging from 48-h lethality through partial and full life-cycle assays, are routinely used for regulatory programs aimed at assessing potential risks of new chemicals such as high-production volume materials and pesticides, as well as impacts of complex mixtures like effluents. The species also has been used for a wide variety of research applications focused on topics like the development of quantitative structure-activity relationship models, mixture toxicity, extrapolation of the effects of chemicals across species, and understanding the results of laboratory assays relative to impacts in the field. Attributes of the fathead minnow also make it an excellent model for addressing new challenges in aquatic toxicology, including identification of sensitive life-stages/endpoints for chemicals with differing modes/mechanisms of action, predicting population-level effects based on data collected from lower levels of biological organization, and exploring/understanding the emerging role of genomics in research and regulation.

Expression analysis of sex-specific and 17beta-estradiol-responsive genes in the Japanese medaka, Oryzias latipes, using oligonucleotide microarrays

Gene profiling of Japanese medaka (Oryzias latipes) was performed using an oligonucleotide DNA microarray representing 22,587 TIGR O. latipes gene indices (OLGIs). The average correlation coefficients for gene expression between individual mature fish were high (>0.95) for both female and male, indicating that the physiological status of medaka is highly reproducible under prescribed growth conditions. Of the 22,587 OLGIs, 2575 showed significant differences in expression between female and male. Exposure to 17beta-estradiol (E2) revealed 381 E2-responsive OLGIs in male medaka. Feminization and male-dysfunction factors of the E2-treated males calculated using the combination of Pearson correlation coefficient and Euclidean distances indicate that E2 treatment "weakly feminized" male medaka, while male physiological functions were not significantly disrupted. This study demonstrates the possibility of using medaka microarrays to estimate the overall effects of hormonally active chemicals.

Soil screening for identifying ecological risk stressors using a battery of in vitro cell bioassays

Soil screening could be a process of identifying and defining areas, contaminants, and condition at the sites that warrant further attention for developing ecological risk assessments. In present work, a total of 41 surface soil samples from Tianjin, China were sampled and the soil organic extracts were evaluated using a battery of in vitro cell bioassays. The battery included ethoxyresorfin O-deethylase (EROD) with H4IIE rat hepatoma cells bioassay for Aryl hydrocarbon (Ah) receptor (Ah-agonists) effects, the SOS/umu bioassay for genotoxic effects, and human estrogen receptor recombinant yeast bioassay for estrogenic effects. The results have showed that total estrogenic effects in these soil samples was measured to be between 0.1 and 14.2ng EEQkg(-1) soil (d.w.); Ah-agonists effects assayed by EROD bioassay varied from 2.8ng TEQkg(-1) soil (d.w.) to 42.6ng TEQkg(-1) soil (d.w.), and the amounts of soil weight required for the extracts to lead positive result (IR 2.0) in the SOS/umu bioassay were between 3.9 and 31.3mg (d.w.) per well. In addition, the geographic distributions of Ah-agonists effects and genotoxic effects in Tianjin area exhibited a strong positive correlation with each other. However, the distribution of estrogenic effects with high levels in northwest Tianjin was markedly different from that of Ah-agonists effects, where the high levels were distributed in the urban of Tianjin, as well as coastal towns. It has been concluded that the toxicity assessment of surface soil using a battery of in vitro cell bioassays could provided meaningful information regarding characterization procedure in ecological risk assessment.

Gene expression patterns in rainbow trout, Oncorhynchus mykiss, exposed to a suite of model toxicants

The increased availability and use of DNA microarrays has allowed the characterization of gene expression patterns associated with exposure to different toxicants. An important question is whether toxicant induced changes in gene expression in fish are sufficiently diverse to allow for identification of specific modes of action and/or specific contaminants. In theory, each class of toxicant may generate a gene expression profile unique to its mode of toxic action. In this study, isogenic (cloned) rainbow trout Oncorhynchus mykiss were exposed to sublethal levels of a series of model toxicants with varying modes of action, including ethynylestradiol (xeno-estrogen), 2,2,4,4'-tetrabromodiphenyl ether (BDE-47, thyroid active), diquat (oxidant stressor), chromium VI, and benzo[a]pyrene (BaP) for a period of 1-3 weeks. An additional experiment measured trenbolone (anabolic steroid; model androgen) induced gene expression changes in sexually mature female trout. Following exposure, fish were euthanized, livers removed and RNA extracted. Fluorescently labeled cDNA were generated and hybridized against a commercially available Atlantic Salmon/Trout array (GRASP project, University of Victoria) spotted with 16,000 cDNA's. The slides were scanned to measure abundance of a given transcript in each sample relative to controls. Data were analyzed via Genespring (Silicon Genetics) to identify a list of up- and downregulated genes, as well as to determine gene clustering patterns that can be used as "expression signatures". The results indicate each toxicant exposure caused between 64 and 222 genes to be significantly altered in expression. Most genes exhibiting altered expression responded to only one of the toxicants and relatively few were co-expressed in multiple treatments. For example, BaP and Diquat, both of which exert toxicity via oxidative stress, upregulated 28 of the same genes, of over 100 genes altered by either treatment. Other genes associated with steroidogenesis, p450 and estrogen responsive genes appear to be useful for selectively identifying toxicant mode of action in fish, suggesting a link between gene expression profile and mode of toxicity. Our array results showed good agreement with quantitative real time polymerase chain reaction (qRT PCR), which demonstrates that the arrays are an accurate measure of gene expression. The specificity of the gene expression profile in response to a model toxicant, the link between genes with altered expression and mode of toxic action, and the consistency between array and qRT PCR results all suggest that cDNA microarrays have the potential to screen environmental contaminants for biomarkers and mode of toxic action.

Derivation of species-specific hybridization-like knowledge out of cross-species hybridization results

BACKGROUND

One of the approaches for conducting genomics research in organisms without extant microarray platforms is to profile their expression patterns by using Cross-Species Hybridization (CSH). Several different studies using spotted microarray and CSH produced contradicting conclusions in the ability of CSH to reflect biological processes described by species-specific hybridization (SSH).

RESULTS

We used a tomato-spotted cDNA microarray to examine the ability of CSH to reflect SSH data. Potato RNA was hybridized to spotted cDNA tomato and potato microarrays to generate CSH and SSH data, respectively. Difficulties arose in obtaining transcriptomic data from CSH that reflected those obtained from SSH. Nevertheless, once the data was filtered for those corresponding to matching probe sets, by restricting proper cutoffs of probe homology, the CSH transcriptome data showed improved reflection of those of the SSH.

CONCLUSIONS

This study evaluated the relative performance of CSH compared to SSH, and proposes methods to ensure that CSH closely reflects the biological process analyzed by SSH.

Application of ecotoxicogenomics for studying endocrine disruption in vertebrates and invertebrates

Chemicals released into the environment potentially disrupt the endocrine system in wild animals and humans. Developing organisms are particularly sensitive to estrogenic chemicals. Exposure to estrogens or estrogenic chemicals during critical periods of development induces persistent changes in both reproductive and nonreproductive organs, including persistent molecular alterations. Estrogen-responsive genes and critical developmental windows of various animal species, therefore, need to be identified for investigators to understand the molecular basis of estrogenic activity during embryonic development. For investigators to understand molecular mechanisms of toxicity in various species, toxicogenomics/ecotoxicogenomics, defined as the integration of genomics (transcriptomics, proteomics, metabolomics) into toxicology and ecotoxicology, need to be established as powerful tools for research. As the initial step toward using genomics to examine endocrine-disrupting chemicals, estrogen receptors and other steroid hormone receptors have been cloned in various species, including reptiles, amphibians, and fish, and alterations in the expression of these genes in response to chemicals were investigated. We are identifying estrogen-responsive genes in mouse reproductive tracts using cDNA microarrays and trying to establish microarray systems in the American alligator, roach, medaka, and water fleas (Daphnia magna). It is too early to define common estrogen-responsive genes in various animal species; however, toxicogenomics and ectotoxicogenomics provide powerful tools to help us understand the molecular mechanism of chemical toxicities in various animal species.

Ten challenges for improved ecotoxicological testing in environmental risk assessment

New regulations, in particular the new European chemicals legislation (REACH), will increase the demands on environmental risk assessment (ERA). The requirements on efficient ecotoxicological testing systems are summarized, and 10 major issues for the improvement of ERA practices are discussed, namely: (1) the choice of representative test species, (2) the development of test systems that are relevant for ecosystems in different parts of the world, (3) the inclusion of sensitive life stages in test systems, (4) the inclusion of endpoints on genetic variation in populations, (5) using mechanistic understanding of toxic effects to develop more informative and efficient test systems, (6) studying disruption in invertebrate endocrine mechanisms, that may differ radically from those we know from vertebrates, (7) developing standardized methodologies for testing of poorly water-soluble substances, (8) taking ethical considerations into account, in particular by reducing the use of vertebrates in ecotoxicological tests, (9) using a systematic (statistical) approach in combination with mechanistic knowledge to combine tests efficiently into testing systems, and (10) developing ERA so that it provides the information needed for precautionary decision-making.

The potential of ecotoxicoproteomics in environmental monitoring: biomarker profiling in mussel plasma using ProteinChip array technology

New global technologies, allowing simultaneous analysis of thousands of genes, proteins, and metabolites (so-called "omics" technologies), are being adopted rapidly by industry, academia, and regulatory agencies. This study evaluated the potential of proteomics in ecotoxicological research (i.e., ecotoxicoproteomics). Filter-feeding mussels (Mytilus edulis) were exposed continuously for 3 wk to oil, or oil spiked with alkylphenols and extra polycyclic aromatic hydrocarbons. The influence of chronic exposure on mussel plasma protein expression was investigated utilizing ProteinChip array technology in combination with surface-enhanced laser desorption/ionization time-of-flight mass spectrometry (SELDI TOF MS). Results indicated that exposure to spiked oil had a more significant effect on protein expression in mussels than oil alone. In total, 83 mass peaks (intact or modified proteins/peptides) were significantly altered by spiked oil, while 49 were altered by oil. In exposed organisms, the majority of peaks were upregulated compared to controls (i.e., 69% in oil and 71% in spiked oil). Some peaks (32 in total) were affected by both treatments; however, the degree of response was higher in the spiked oil group for 25 of the 32 commonly affected features. Additionally, certain peaks revealed exposure- or gender-specific responses. Multivariate analysis with regression tree-based methods detected protein patterns associated with exposure that correctly classified masked samples with 90-95% accuracy. Similarly, 92% of females and 85% of males were correctly classified (independent of exposure). Results indicate that proteomics have the potential to make a valuable contribution to environmental monitoring and risk assessment.

Molecular impact of propiconazole on Daphnia magna using a reproduction-related cDNA array

We have developed a first version cDNA microarray of the cladoceran Daphnia magna. Through Suppression Subtractive Hybridisation PCR (SSH-PCR) 855 life stage-specific cDNAs were collected and used to document the toxicological mode of action of the pesticide propiconazole. DNA sequencing analysis revealed gene fragments related to important functional classes such as embryo development, energy metabolism, molting and cell cycle. Major changes in transcription were observed in organisms exposed for 4 and 8 days to 1 microg/mL. After 4 days a 3-fold down-regulation of the gene encoding the yolk protein, vitellogenin, was observed indicating impaired oocyte maturation. Moreover, genes such as a larval-specific gene and chaperonin were repressed, whereas the heat shock 90 protein and ATP synthase were induced. Organismal effects clearly confirmed the major molecular findings: at the highest concentration (1 microg/mL) adult growth was significantly (p < 0.05) impaired and increased developmental effects in the offspring could be noted. We have demonstrated the potential of microarray analysis in toxicity screening with D. magna. The use of vitellogenin mRNA as a rapid biomarker of reproductive effects in chronic toxicity studies with cladocerans is suggested.

Recent applications of DNA microarray technology to toxicology and ecotoxicology

Gene expression is a unique way of characterizing how cells and organisms adapt to changes in the external environment. The measurements of gene expression levels upon exposure to a chemical can be used both to provide information about the mechanism of action of the toxicant and to form a sort of "genetic signature" for the identification of toxic products. The development of high-quality, commercially available gene arrays has allowed this technology to become a standard tool in molecular toxicology. Several national and international initiatives have provided the proof-of-principle tests for the application of gene expression for the study of the toxicity of new and existing chemical compounds. In the last few years the field has progressed from evaluating the potential of the technology to illustrating the practical use of gene expression profiling in toxicology. The application of gene expression profiling to ecotoxicology is at an earlier stage, mainly because of the the many variables involved in analyzing the status of natural populations. Nevertheless, significant studies have been carried out on the response to environmental stressors both in model and in nonmodel organisms. It can be easily predicted that the development of stressor-specific signatures in gene expression profiling in ecotoxicology will have a major impact on the ecotoxicology field in the near future. International collaborations could play an important role in accelerating the application of genomic approaches in ecotoxicology.

Analysis of expressed sequence tags of the water flea Daphnia magna

To study gene expression in the water flea Daphnia magna we constructed a cDNA library and characterized the expressed sequence tags (ESTs) of 7210 clones. The EST sequences clustered into 2958 nonredundant groups. BLAST analyses of both protein and DNA databases showed that 1218 (41%) of the unique sequences shared significant similarities to known nucleotide or amino acid sequences, whereas the remaining 1740 (59%) showed no significant similarities to other genes. Clustering analysis revealed particularly high expression of genes related to ATP synthesis, structural proteins, and proteases. The cDNA clones and EST sequence information should be useful for future functional analysis of daphnid biology and investigation of the links between ecology and genomics.

Fish as models for environmental genomics

Fish offer important advantages for defining the organism-environment interface and responses to natural or anthropogenic stressors. Genomic approaches using fish promise increased investigative power, and have already provided insights into the mechanisms that underlie short-term and long-term environmental adaptations. The range of fish species for which genomic resources are available is increasing, but will require significant further expansion for the optimal application of fish environmental genomics.

Xenopus as a model organism in developmental chemical genetic screens

Chemical genetics is a potentially powerful tool for studying developmental processes in vertebrate systems. We present data showing Xenopus laevis as a model organism in which systematic chemical genetic screens can be carried out. Previous forward chemical genetic screens, including those with developing zebrafish embryos, have demonstrated the nature and value of biological information gained with this approach. We show how amenable Xenopus is to chemical genetics by investigating a series of compounds either with known biochemical effects, or previously identified to give developmental phenotypes, on a range of biological functions, including the development of pigmentation, the heart and the central nervous system in zebrafish. We have found that the compounds give comparable phenotypes when applied to developing Xenopus embryos. We have also studied the penetrance and expressivity of these chemical genetic phenotypes in relation to genetic variation and the developmental window during which the compound is present. Finally, we assess the feasibility and the potential throughput of a screen in this vertebrate species.

Resolving mechanisms of toxicity while pursuing ecotoxicological relevance?

In this age of modern biology, aquatic toxicological research has pursued mechanisms of action of toxicants. This has provided potential tools for ecotoxicologic investigations. However, problems of biocomplexity and issues at higher levels of biological organization remain a challenge. In the 1980s and 1990s and continuing to a lesser extent today, organisms residing in highly contaminated field sites or exposed in the laboratory to calibrated concentrations of individual compounds were carefully analyzed for their responses to priority pollutants. Correlation of biochemical and structural analyses in cultured cells and tissues, as well as the in vivo exposures led to the production and application of biomarkers of exposure and effect and to our awareness of genotoxicity and its chronic manifestations, such as neoplasms, in wild fishes. To gain acceptance of these findings in the greater environmental toxicology community, "validation of the model" versus other, better-established often rodent models, was necessary and became a major focus. Resultant biomarkers were applied to heavily contaminated and reference field sites as part of effects assessment and with investigations following large-scale disasters such as oil spills or industrial accidents. Over the past 15 years, in the laboratory, small aquarium fish models such as medaka (Oryzias latipes), zebrafish (Danio rerio), platyfish (Xiphophorus species), fathead minnow (Pimephales promelas), and sheepshead minnow (Cyprinodon variegatus) were increasingly used establishing mechanisms of toxicants. Today, the same organisms provide reliable information at higher levels of biological organization relevant to ecotoxicology. We review studies resolving mechanisms of toxicity and discuss ways to address biocomplexity, mixtures of contaminants, and the need to relate individual level responses to populations and communities.

An integrated biomarker-based strategy for ecotoxicological evaluation of risk in environmental management

Environmental impacts by both natural events and man-made interventions are a fact of life; and developing the capacity to minimise these impacts and their harmful consequences for biological resources, ecosystems and human health is a daunting task for environmental legislators and regulators. A major challenge in impact and risk assessment, as part of integrated environmental management (IEM), is to link harmful effects of pollution (including toxic chemicals) in individual sentinel animals to their ecological consequences. This obstacle has resulted in a knowledge-gap for those seeking to develop effective policies for sustainable use of resources and environmental protection. Part of the solution to this problem may lie with the use of diagnostic clinical-type laboratory-based ecotoxicological tests or biomarkers, utilising sentinel animals as integrators of pollution, coupled with direct immunochemical tests for contaminants. These rapid and cost-effective ecotoxicological tools can provide information on the health status of individuals and populations based on relatively small samples of individuals. In the context of ecosystem status or health of the environment, biomarkers are also being used to link processes of molecular and cellular damage through to higher levels (i.e., prognostic capability), where they can result in pathology with reduced physiological performance and reproductive success. Complex issues are involved in evaluating environmental risk, such as the effects of the physico-chemical environment on the speciation and uptake of pollutant chemicals and inherent inter-individual and inter-species differences in vulnerability to toxicity; and the toxicity of complex mixtures. Effectively linking the impact of pollutants through the various hierarchical levels of biological organisation to ecosystem and human health requires a pragmatic integrated approach based on existing information that either links or correlates processes of pollutant uptake, detoxication and pathology with each other and higher level effects. It is further proposed here that this process will be facilitated by pursuing a holistic or whole systems approach with the development of computational simulation models of cells, organs and animals in tandem with empirical data (i.e., the middle-out approach). In conclusion, an effective integrated environmental management strategy to secure resource sustainability requires an integrated capability for risk assessment and prediction. Furthermore, if such a strategy is to influence and help in the formulation of environmental policy decisions, then it is crucial to demonstrate scientific robustness of predictions concerning the long-term consequences of pollution to politicians, industrialists and environmental managers; and also increase stakeholder awareness of environmental problems.