Transcriptomic and Network Analyses Reveal Mechanistic-Based Biomarkers of Endocrine Disruption in the Marine Mussel, Mytilus edulis.

Transcriptomics, high-throughput assays, and adverse outcome pathways (AOP) are promising approaches applied to toxicity monitoring in the 21st century, but development of these methods is challenging for nonmodel organisms and emerging contaminants. For example, Endocrine Disrupting Compounds (EDCs) may cause reproductive impairments and feminization of male bivalves; however, the mechanism linked to this adverse outcome is unknown. To develop mechanism-based biomarkers that may be linked through an AOP, we exposed Mytilus edulis to 17-alpha-ethinylestradiol (5 and 50 ng/L) and 4-nonylphenol (1 and 100 µg/L) for 32 and 39 days. When mussels were exposed to these EDCs, we found elevated female specific transcripts and significant female-skewed sex ratios using a RT-qPCR assay. We performed gene expression analysis on digestive gland tissue using an M. edulis microarray and through network and targeted analyses identified the nongenomic estrogen signaling pathway and steroidogenesis pathway as the likely mechanisms of action for a putative AOP. We also identified several homologues to genes within the vertebrate steroidogenesis pathway including the cholesterol side chain cleavage complex. From this AOP, we designed the Coastal Biosensor for Endocrine Disruption (C-BED) assay which was confirmed in the laboratory and tested in the field.

Gene transcription, metabolite and lipid profiling in eco-indicator daphnia magna indicate diverse mechanisms of toxicity by legacy and emerging flame-retardants.

The use of chemical flame-retardants (FR) in consumer products has steadily increased over the last 30 years. Toxicity data exist for legacy FRs such as pentabromodiphenyl ether (pentaBDE), but less is known about effects of new formulations. To address this issue, the toxicity of seven FR chemicals and formulations was assessed on the freshwater crustacean Daphnia magna. Acute 48-h nominal LC50 values for penta- and octabromodiphenyl ether (pentaBDE, octaBDE), Firemaster 550 (FM550), Firemaster BZ-54 (BZ54), bis(2-ethylhexyl) tetrabromophthalate (BEH-TEBP), triphenyl phosphate (TPhP), and nonbrominated BEH-TEBP analog bis(2-ethylhexyl) phthalate (BEHP) ranged from 0.058 mg/L (pentaBDE) to 3.96 mg/L (octaBDE). mRNA expression, (1)H NMR-based metabolomic and lipidomic profiling at 1/10 LC50 revealed distinct patterns of molecular response for each exposure, suggesting pentaPBDE affects transcription and translation, octaBDE and BEH-TEBP affect glycosphingolipid biosynthesis and BZ54 affects Wnt and Hedgehog signal pathways as well as glycosaminoglycan degradation. Brominated components of FM550 (i.e., BZ54) were significantly higher in Daphnia after 48 h following 1/10 LC50 exposure. FM550 elicited significant mRNA changes at five concentrations across a range from 1/10(6) LC50 to 1/2 LC50. Analyses suggest FM550 impairs nutrient utilization or uptake in Daphnia.

Effects of strain and age on hepatic gene expression profiles in murine models of HFE-associated hereditary hemochromatosis.

Hereditary hemochromatosis is an iron overload disorder most commonly caused by a defect in the HFE gene. While the genetic defect is highly prevalent, the majority of individuals do not develop clinically significant iron overload, suggesting the importance of genetic modifiers. Murine hfe knockout models have demonstrated that strain background has a strong effect on the severity of iron loading. We noted that hepatic iron loading in hfe-/- mice occurs primarily over the first postnatal weeks (loading phase) followed by a timeframe of relatively static iron concentrations (plateau phase). We thus evaluated the effects of background strain and of age on hepatic gene expression in Hfe knockout mice (hfe-/-). Hepatic gene expression profiles were examined using cDNA microarrays in 4- and 8-week-old hfe-/- and wild-type mice on two different genetic backgrounds, C57BL/6J (C57) and AKR/J (AKR). Genes differentially regulated in all hfe-/- mice groups, compared with wild-type mice, including those involved in cell survival, stress and damage responses and lipid metabolism. AKR strain-specific changes in lipid metabolism genes and C57 strain-specific changes in cell adhesion and extracellular matrix protein genes were detected in hfe-/- mice. Mouse strain and age are each significantly associated with hepatic gene expression profiles in hfe-/- mice. These affects may underlie or reflect differences in iron loading in these mice.

Modulation of Ras signaling alters the toxicity of hydroquinone, a benzene metabolite and component of cigarette smoke.

BACKGROUND: Benzene is an established human leukemogen, with a ubiquitous environmental presence leading to significant population exposure. In a genome-wide functional screen in the yeast Saccharomyces cerevisiae, inactivation of IRA2, a yeast ortholog of the human tumor suppressor gene NF1 (Neurofibromin), enhanced sensitivity to hydroquinone, an important benzene metabolite. Increased Ras signaling is implicated as a causal factor in the increased pre-disposition to leukemia of individuals with mutations in NF1. METHODS: Growth inhibition of yeast by hydroquinone was assessed in mutant strains exhibiting varying levels of Ras activity. Subsequently, effects of hydroquinone on both genotoxicity (measured by micronucleus formation) and proliferation of WT and Nf1 null murine hematopoietic precursors were assessed. RESULTS: Here we show that the Ras status of both yeast and mammalian cells modulates hydroquinone toxicity, indicating potential synergy between Ras signaling and benzene toxicity. Specifically, enhanced Ras signaling increases both hydroquinone-mediated growth inhibition in yeast and genotoxicity in mammalian hematopoetic precursors as measured by an in vitro erythroid micronucleus assay. Hydroquinone also increases proliferation of CFU-GM progenitor cells in mice with Nf1 null bone marrow relative to WT, the same cell type associated with benzene-associated leukemia. CONCLUSIONS: Together our findings show that hydroquinone toxicity is modulated by Ras signaling. Individuals with abnormal Ras signaling could be more vulnerable to developing myeloid diseases after exposure to benzene. We note that hydroquinone is used cosmetically as a skin-bleaching agent, including by individuals with cafe-au-lait spots (which may be present in individuals with neurofibromatosis who have a mutation in NF1), which could be unadvisable given our findings.

Silver nanowire exposure results in internalization and toxicity to Daphnia magna.

Nanowires (NWs), high-aspect-ratio nanomaterials, are increasingly used in technological materials and consumer products and may have toxicological characteristics distinct from nanoparticles. We carried out a comprehensive evaluation of the physicochemical stability of four silver nanowires (AgNWs) of two sizes and coatings and their toxicity to Daphnia magna . Inorganic aluminum-doped silica coatings were less effective than organic poly(vinyl pyrrolidone) coatings at preventing silver oxidation or Ag(+) release and underwent a significant morphological transformation within 1 h following addition to low ionic strength Daphnia growth media. All AgNWs were highly toxic to D. magna but less toxic than ionic silver. Toxicity varied as a function of AgNW dimension, coating, and solution chemistry. Ag(+) release in the media could not account for observed AgNW toxicity. Single-particle inductively coupled plasma mass spectrometry distinguished and quantified dissolved and nanoparticulate silver in microliter-scale volumes of Daphnia magna hemolymph with a limit of detection of approximately 10 ppb. The silver levels within the hemolymph of Daphnia exposed to both Ag(+) and AgNW met or exceeded the initial concentration in the growth medium, indicating effective accumulation during filter feeding. Silver-rich particles were the predominant form of silver in hemolymph following exposure to both AgNWs and Ag(+). Scanning electron microscopy imaging of dried hemolymph found both AgNWs and silver precipitates that were not present in the AgNW stock or the growth medium. Both organic and inorganic coatings on the AgNW were transformed during ingestion or absorption. Pathway, gene ontology, and clustering analyses of gene expression response indicated effects of AgNWs distinct from ionic silver on Daphnia magna .

Toxicogenomic responses of nanotoxicity in Daphnia magna exposed to silver nitrate and coated silver nanoparticles.

Applications for silver nanomaterials in consumer products are rapidly expanding, creating an urgent need for toxicological examination of the exposure potential and ecological effects of silver nanoparticles (AgNPs). The integration of genomic techniques into environmental toxicology has presented new avenues to develop exposure biomarkers and investigate the mode of toxicity of novel chemicals. In the present study we used a 15k oligonucleotide microarray for Daphnia magna, a freshwater crustacean and common indicator species for toxicity, to differentiate between particle specific and ionic silver toxicity and to develop exposure biomarkers for citrate-coated and PVP-coated AgNPs. Gene expression profiles revealed that AgNO(3) and AgNPs have distinct expression profiles suggesting different modes of toxicity. Major biological processes disrupted by the AgNPs include protein metabolism and signal transduction. In contrast, AgNO(3) caused a downregulation of developmental processes, particularly in sensory development. Metal responsive and DNA damage repair genes were induced by the PVP AgNPs, but not the other treatments. In addition, two specific biomarkers were developed for the environmental detection of PVP AgNPs; although further verification under different environmental conditions is needed.

Genome-Wide Functional and Stress Response Profiling Reveals Toxic Mechanism and Genes Required for Tolerance to Benzo[a]pyrene in S. cerevisiae.

Benzo[a]pyrene (BaP) is a ubiquitous, potent, and complete carcinogen resulting from incomplete organic combustion. BaP can form DNA adducts but other mechanisms may play a role in toxicity. We used a functional toxicology approach in S. cerevisiae to assess the genetic requirements for cellular resistance to BaP. In addition, we examined translational activities of key genes involved in various stress response pathways. We identified multiple genes and processes involved in modulating BaP toxicity in yeast which support DNA damage as a primary mechanism of toxicity, but also identify other potential toxicity pathways. Gene ontology enrichment analysis indicated that DNA damage and repair as well as redox homeostasis and oxidative stress are key processes in cellular response to BaP suggesting a similar mode of action of BaP in yeast and mammals. Interestingly, toxicant export is also implicated as a potential novel modulator of cellular susceptibility. In particular, we identified several transporters with human orthologs (solute carrier family 22) which may play a role in mammalian systems.

Metabolomics of microliter hemolymph samples enables an improved understanding of the combined metabolic and transcriptional responses of Daphnia magna to cadmium.

Omic technologies offer unprecedented opportunities to better understand mode(s)-of-toxicity and downstream secondary effects by providing a holistic view of the molecular changes underlying physiological disruption. Crustacean hemolymph represents a largely untapped biochemical resource for such toxicity studies. We sought to characterize changes in the hemolymph metabolome and whole-body transcriptome to reveal early processes leading to chronic toxicity in the indicator species, Daphnia magna, after 24-h sublethal cadmium exposure (18 µg/L, corresponding to 1/10 LC(50)). We first confirmed that metabolites can be detected and identified in small volumes (∼3-6 µL) of D. magna hemolymph using Fourier transform ion cyclotron resonance mass spectrometry and NMR spectroscopy. Subsequently, mass spectrometry based metabolomics of hemolymph identified disruption to two major classes of metabolites: amino acids and fatty acids. These findings were compared to differentially expressed genes identified by a D. magna 44k oligonucleotide microarray, which included decreased levels of digestive enzymes and increased expression of cuticle proteins and oxidative stress response genes. The combination of metabolic and transcriptional changes revealed through KEGG pathway analysis and gene ontology, respectively, enabled a more complete understanding of how cadmium disrupts nutrient uptake and metabolism, ultimately resulting in decreased energy reserves and chronic toxicity.

Reverse engineering adverse outcome pathways.

The toxicological effects of many stressors are mediated through unknown, or incompletely characterized, mechanisms of action. The application of reverse engineering complex interaction networks from high dimensional omics data (gene, protein, metabolic, signaling) can be used to overcome these limitations. This approach was used to characterize adverse outcome pathways (AOPs) for chemicals that disrupt the hypothalamus-pituitary-gonadal endocrine axis in fathead minnows (FHM, Pimephales promelas). Gene expression changes in FHM ovaries in response to seven different chemicals, over different times, doses, and in vivo versus in vitro conditions, were captured in a large data set of 868 arrays. Potential AOPs of the antiandrogen flutamide were examined using two mutual information-based methods to infer gene regulatory networks and potential AOPs. Representative networks from these studies were used to predict network paths from stressor to adverse outcome as candidate AOPs. The relationship of individual chemicals to an adverse outcome can be determined by following perturbations through the network in response to chemical treatment, thus leading to the nodes associated with the adverse outcome. Identification of candidate pathways allows for formation of testable hypotheses about key biological processes, biomarkers, or alternative endpoints that can be used to monitor an AOP. Finally, the unique challenges facing the application of this approach in ecotoxicology were identified and a road map for the utilization of these tools presented.

Gene expression profiling in Daphnia magna part I: concentration-dependent profiles provide support for the No Observed Transcriptional Effect Level.

Ecotoxicogenomic approaches to environmental monitoring provide holistic information, offer insight into modes of action, and help to assess the causal agents and potential toxicity of effluents beyond the traditional end points of death and reproduction. Recent investigations of toxicant exposure indicate dose-dependent changes are a key issue in interpreting genomic studies. Additionally, there is interest in developing methods to integrate gene expression studies in environmental monitoring and regulation, and the No Observed Transcriptional Effect Level (NOTEL) has been proposed as a means for screening effluents and unknown chemicals fortoxicity. However, computational methods to determine the NOTEL have yet to be established. Therefore, we examined effects on gene expression in Daphnia magna following exposure to Cu, Cd, and Zn over a range of concentrations including a tolerated, a sublethal, and a nearly acutely toxic concentration. Each concentration produced a distinct gene expression profile. We observed differential expression of a very few genes at tolerated concentrations that were distinct from the expression profiles observed at concentrations associated with toxicity. These results suggest that gene expression analysis may offer a strategy for distinguishing toxic and nontoxic concentrations of metals in the environment and provide support for a NOTEL for metal exposure in D. magna. Mechanistic insights could be inferred from the concentration-dependent gene expression profiles including metal specific effects on disparate metabolic processes such as digestion, immune response, development and reproduction, and less specific stress responses at higher concentrations.

Gene expression profiling in Daphnia magna, part II: validation of a copper specific gene expression signature with effluent from two copper mines in California.

Genomic technologies show great potential for classifying disease states and toxicological impacts from exposure to chemicals into functional categories. In environmental monitoring, the ability to classify field samples and predict the pollutants present in these samples could contribute to monitoring efforts and the diagnosis of contaminated sites. Using gene expression analysis, we challenged our custom Daphnia magna cDNA microarray to determine the presence of a specific metal toxicant in blinded field samples collected from two copper mines in California. We compared the gene expression profiles from our field samples to previously established expression profiles for Cu, Cd, and Zn. The expression profiles from the Cu-containing field samples clustered with the laboratory-exposed Cu-specific gene expression profiles and included genes previously identified as copper biomarkers, verifying that gene expression analysis can predict environmental exposure to a specific pollutant. In addition, our study revealed that upstream field samples containing undetectable levels of Cu caused the differential expression of only a few genes, lending support for the concept of a no observed transcriptional effect level (NOTEL). If confirmed by further studies, the NOTEL may play an important role in discriminating polluted and nonpolluted sites in future monitoring efforts.

Daphnia magna ecotoxicogenomics provides mechanistic insights into metal toxicity.

Toxicogenomics has provided innovative approaches to chemical screening, risk assessment, and predictive toxicology. If applied to ecotoxicology, genomics tools could greatly enhance the ability to understand the modes of toxicity in environmentally relevant organisms. Daphnia magna, a small aquatic crustacean, is considered a "keystone" species in ecological food webs and is an indicator species for toxicant exposure. Our objective was to demonstrate the potential utility of gene expression profiling in ecotoxicology by identifying novel biomarkers and uncovering potential modes of action in D. magna. Using a custom D. magna cDNA microarray, we identified distinct expression profiles in response to sublethal copper, cadmium, and zinc exposures and discovered specific biomarkers of exposure including two probable metallothioneins, and a ferritin mRNA with a functional IRE. The gene expression patterns support known mechanisms of metal toxicity and reveal novel modes of action including zinc inhibition of chitinase activity. By integrating gene expression profiling into an environmentally important organism, this study provides experimental support for the utility of ecotoxicogenomics.