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

2,4-Dinitrotoluene (2,4-DNT) exposure induces liver developmental damage and perturbs lipid metabolism and oxygen transport gene expression in zebrafish (Danio rerio)

2,4-Dinitrotoluene (2, 4-DNT) is a common environmental pollutant. The toxic effect on mammals of 2,4-DNT has been well studied, but its toxicity on aquatic organisms is little known. In this study, 126 healthy female zebrafish (Danio rerio) were exposed to different concentrations of 2,4-DNT (0, 2, 4, 8, 12 and 16 mg/L) to determine 96-h semi-lethal concentrations (LC₅₀). And then, 90 female zebrafish were exposed to 0, 2, 4 and 8 mg/L 2,4-DNT for 5 days to study liver toxicity. Exposed zebrafish developed hypoxia features, such as floating head and breathing rapidly, and then died. 96-h LC₅₀ of 2,4-DNT in zebrafish was 9.36 mg/L. Histological data revealed that 2,4-DNT severely damaged the liver tissues, following with the round nucleus, dense interstitial tissue, dense arranged hepatocyte cords and more inflammatory cells. Additionally, the further result showed that the lower levels of lipid transport and metabolism (apoα2, mtp, ppar-α and acox) were noticed. But, exposed to 2,4-DNT for 5 days significantly upregulated the expression levels of genes involved in respiration (hif1a, tfa and ho1, p < 0.05). These results indicated that 2,4-DNT exposure disturbed lipid transport and metabolism and oxygen supply in zebrafish, which could contribute to severe damage in liver and death.

Example of Adverse Outcome Pathway Concept Enabling Genome-to-Phenome Discovery in Toxicology

The following article represents a mini-review of an intensive 10-year progression of genome-to-phenome (G2P) discovery guided by the adverse outcome pathway (AOP) concept. This example is presented as a means to stimulate crossover of this toxicological concept to enhance G2P discovery within the broader biological sciences community. The case study demonstrates the benefits of the AOP approach for establishing causal linkages across multiple levels of biological organization ultimately linking molecular initiation (often at the genomic scale) to organism-level phenotypes of interest. The case study summarizes a US military effort to identify the mechanism(s) underlying toxicological phenotypes of lethargy and weight loss in response to nitroaromatic munitions exposures, such as 2,4,6-trinitrotoluene. Initial key discoveries are described including the toxicogenomic results that nitrotoluene exposures inhibited expression within the peroxisome proliferator activated receptor α (PPARα) pathway. We channeled the AOP concept to test the hypothesis that inhibition of PPARα signaling in nitrotoluene exposures impacted lipid metabolic processes, thus affecting systemic energy budgets, ultimately resulting in body weight loss. Results from a series of transcriptomic, proteomic, lipidomic, in vitro PPARα nuclear signaling, and PPARα knock-out investigations ultimately supported various facets of this hypothesis. Given these results, we next proceeded to develop a formalized AOP description of PPARα antagonism leading to body weight loss. This AOP was refined through intensive literature review and polished through multiple rounds of peer-review leading to final international acceptance as an Organisation for Economic Cooperation and Development-approved AOP. Briefly, that AOP identifies PPARα antagonist binding as the molecular initiating event (MIE) leading to a series of key events including inhibition of nuclear transactivation for genes controlling lipid metabolism and ketogenesis, inhibition of fatty acid beta-oxidation and ketogenesis dynamics, negative energy budget, and ultimately the adverse outcome (AO) of body-weight loss. Given that the PPARα antagonism MIE represented a reliable indicator of AO progression within the pathway, a phylogenetic analysis was conducted which indicated that PPARα amino acid relatedness generally tracked species relatedness. Additionally, PPARα amino acid relatedness analysis using the Sequence Alignment to Predict Across Species Susceptibility predicted susceptibility to the MIE across vertebrates providing context for AOP extrapolation across species. Overall, we hope this illustrative example of how the AOP concept has benefited toxicology sows a seed within the broader biological sciences community to repurpose the concept to facilitate enhanced G2P discovery in biology.

2,4-Dinitrotoluene (DNT) Perturbs Yolk Absorption, Liver Development and Lipid Metabolism/Oxygen Transport Gene Expression in Zebrafish Embryos and Larvae

2,4-dinitrotoluene (2,4-DNT) is a common environmental pollutant, and was classified as a group 2B human carcinogenic compound by the International Agency for Research on Cancer. This study determined the toxic effects of 2,4-DNT exposure on zebrafish at the embryo-larvae stage, in terms of organ morphogenesis and the expression pattern of selected target genes related to lipid metabolism and oxygen transportation. The results showed that the 120-h post-fertilization LC50 of 2,4-DNT was 9.59 mg/L with a 95% confidence interval of 8.89-10.44 mg/L. The larvae treated with 2,4-DNT showed toxic symptoms including smaller body, less skin pigment production, yolk malabsorption, and disordered liver development. Further studies on the expression of genes related to lipid transport and metabolism, and respiration indicated that they were significantly affected by 2,4-DNT. It is concluded that 2,4-DNT exposure perturbed liver development and yolk absorption in early-life zebrafish, and disturbed the lipid metabolism /oxygen transport gene expression.

Multiple environmental stressors induce complex transcriptomic responses indicative of phenotypic outcomes in Western fence lizard

BACKGROUND

The health and resilience of species in natural environments is increasingly challenged by complex anthropogenic stressor combinations including climate change, habitat encroachment, and chemical contamination. To better understand impacts of these stressors we examined the individual- and combined-stressor impacts of malaria infection, food limitation, and 2,4,6-trinitrotoluene (TNT) exposures on gene expression in livers of Western fence lizards (WFL, Sceloporus occidentalis) using custom WFL transcriptome-based microarrays.

RESULTS

Computational analysis including annotation enrichment and correlation analysis identified putative functional mechanisms linking transcript expression and toxicological phenotypes. TNT exposure increased transcript expression for genes involved in erythropoiesis, potentially in response to TNT-induced anemia and/or methemoglobinemia and caused dose-specific effects on genes involved in lipid and overall energy metabolism consistent with a hormesis response of growth stimulation at low doses and adverse decreases in lizard growth at high doses. Functional enrichment results were indicative of inhibited potential for lipid mobilization and catabolism in TNT exposures which corresponded with increased inguinal fat weights and was suggestive of a decreased overall energy budget. Malaria infection elicited enriched expression of multiple immune-related functions likely corresponding to increased white blood cell (WBC) counts. Food limitation alone enriched functions related to cellular energy production and decreased expression of immune responses consistent with a decrease in WBC levels.

CONCLUSIONS

Despite these findings, the lizards demonstrated immune resilience to malaria infection under food limitation with transcriptional results indicating a fully competent immune response to malaria, even under bio-energetic constraints. Interestingly, both TNT and malaria individually increased transcriptional expression of immune-related genes and increased overall WBC concentrations in blood; responses that were retained in the TNT x malaria combined exposure. The results demonstrate complex and sometimes unexpected responses to multiple stressors where the lizards displayed remarkable resiliency to the stressor combinations investigated.

Evaluation of ecotoxicological effects of benzophenone UV filters: Luminescent bacteria toxicity, genotoxicity and hormonal activity

The widespread use of organic ultraviolet (UV) filters in personal care products raises concerns about their potentially hazardous effects on human and ecosystem health. In this study, the toxicities of four commonly used benzophenones (BPs) UV filters including benzophenone (BP), 2-Hydroxybenzophenone (2HB), 2-Hydroxy-4-methoxybenzophenone (BP3), and 2-Hydroxy-4-methoxybenzophenone-5-sulfonicacid (BP4) in water were assayed in vitro using Vibrio fischeri, SOS/umu assay, and yeast estrogen screen (YES) assay, as well as in vivo using zebrafish larvae. The results showed that the luminescent bacteria toxicity, expressed as logEC₅₀, increased with the lipophilicity (logKow) of BPs UV filters. Especially, since 2HB, BP3 and BP4 had different substituent groups, namely -OH, -OCH₃ and -SO₃H, respectively, these substituent functional groups had a major contribution to the lipophilicity and acute toxicity of these BPs. Similar tendency was observed for the genotoxicity, expressed as the value of induction ratio=1.5. Moreover, all the target BPs UV filters showed estrogenic activity, but no significant influences of lipophilicity on the estrogenicity were observed, with BP3 having the weakest estrogenic efficiency in vitro. Although BP3 displayed no noticeable adverse effects in any in vitro assays, multiple hormonal activities were observed in zebrafish larvae including estrogenicity, anti-estrogenicity and anti-androgenicity by regulating the expression of target genes. The results indicated potential hazardous effects of BPs UV filters and the importance of the combination of toxicological evaluation methods including in vitro and in vivo assays.

Gene Expression Profiling in Fish Toxicology: A Review

In this review, we present an overview of transcriptomic responses to chemical exposures in a variety of fish species. We have discussed the use of several molecular approaches such as northern blotting, differential display reverse transcription-polymerase chain reaction (DDRT-PCR), suppression subtractive hybridization (SSH), real time quantitative PCR (RT-qPCR), microarrays, and next-generation sequencing (NGS) for measuring gene expression. These techniques have been mainly used to measure the toxic effects of single compounds or simple mixtures in laboratory conditions. In addition, only few studies have been conducted to examine the biological significance of differentially expressed gene sets following chemical exposure. Therefore, future studies should focus more under field conditions using a multidisciplinary approach (genomics, proteomics and metabolomics) to understand the synergetic effects of multiple environmental stressors and to determine the functional significance of differentially expressed genes. Nevertheless, recent developments in NGS technologies and decreasing costs of sequencing holds the promise to uncover the complexity of anthropogenic impacts and biological effects in wild fish populations.

Assessment of multiple hormone activities of a UV-filter (octocrylene) in zebrafish (Danio rerio)

In this study, zebrafish (Danio rerio) were exposed to a UV-filter-octocrylene (OCT) with elevated concentrations for 28 d. The total body accumulation of OCT in zebrafish was found to reach 2321.01 ("L" level), 31,234.80 ("M" level), and 70,593.38 ng g(-1) ("H" level) when the average OCT exposure concentration was controlled at 28.61, 505.62, and 1248.70 μg L(-1), respectively. Gross and histological observations as well as RT-qPCR analysis were conducted to determine the effects of OCT accumulation on zebrafish. After exposure, the gonad-somatic index and percentage of vitellogenic oocytes were found to increase significantly in the ovaries of female zebrafish at the H accumulation level. Significant up-regulation of esr1 and cyp19b were observed in the gonads, as well as vtg1 in the livers for both female and male zebrafish. At M and H accumulation levels, apparent down-regulation of ar was observed in the ovaries and testis of the female and male zebrafish, respectively. Although the extent of the effects on zebrafish differed at different accumulation levels, the induction of vtg1 and histological changes in the ovaries are indications of estrogenic activity and the inhibition of esr1 and ar showed antiestrogenic and antiandrogenic activity, respectively. Thus, as OCT could easily accumulate in aquatic life such as zebrafish, one of its most of concern hazards would be the disturbance of the histological development and its multiple hormonal activities.

Systems toxicology identifies mechanistic impacts of 2-amino-4,6-dinitrotoluene (2A-DNT) exposure in Northern Bobwhite

Background

A systems toxicology investigation comparing and integrating transcriptomic and proteomic results was conducted to develop holistic effects characterizations for the wildlife bird model, Northern bobwhite (Colinus virginianus) dosed with the explosives degradation product 2-amino-4,6-dinitrotoluene (2A-DNT). A subchronic 60d toxicology bioassay was leveraged where both sexes were dosed via daily gavage with 0, 3, 14, or 30 mg/kg-d 2A-DNT. Effects on global transcript expression were investigated in liver and kidney tissue using custom microarrays for C. virginianus in both sexes at all doses, while effects on proteome expression were investigated in liver for both sexes and kidney in males, at 30 mg/kg-d.

Results

As expected, transcript expression was not directly indicative of protein expression in response to 2A-DNT. However, a high degree of correspondence was observed among gene and protein expression when investigating higher-order functional responses including statistically enriched gene networks and canonical pathways, especially when connected to toxicological outcomes of 2A-DNT exposure. Analysis of networks statistically enriched for both transcripts and proteins demonstrated common responses including inhibition of programmed cell death and arrest of cell cycle in liver tissues at 2A-DNT doses that caused liver necrosis and death in females. Additionally, both transcript and protein expression in liver tissue was indicative of induced phase I and II xenobiotic metabolism potentially as a mechanism to detoxify and excrete 2A-DNT. Nuclear signaling assays, transcript expression and protein expression each implicated peroxisome proliferator-activated receptor (PPAR) nuclear signaling as a primary molecular target in the 2A-DNT exposure with significant downstream enrichment of PPAR-regulated pathways including lipid metabolic pathways and gluconeogenesis suggesting impaired bioenergetic potential.

Conclusion

Although the differential expression of transcripts and proteins was largely unique, the consensus of functional pathways and gene networks enriched among transcriptomic and proteomic datasets provided the identification of many critical metabolic functions underlying 2A-DNT toxicity as well as impaired PPAR signaling, a key molecular initiating event known to be affected in di- and trinitrotoluene exposures.

Electronic supplementary material

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

A Bayesian model for the identification of differentially expressed genes in Daphnia magna exposed to munition pollutants

In this article we propose a Bayesian hierarchical model for the identification of differentially expressed genes in Daphnia magna organisms exposed to chemical compounds, specifically munition pollutants in water. The model we propose constitutes one of the very first attempts at a rigorous modeling of the biological effects of water purification. We have data acquired from a purification system that comprises four consecutive purification stages, which we refer to as "ponds," of progressively more contaminated water. We model the expected expression of a gene in a pond as the sum of the mean of the same gene in the previous pond plus a gene-pond specific difference. We incorporate a variable selection mechanism for the identification of the differential expressions, with a prior distribution on the probability of a change that accounts for the available information on the concentration of chemical compounds present in the water. We carry out posterior inference via MCMC stochastic search techniques. In the application, we reduce the complexity of the data by grouping genes according to their functional characteristics, based on the KEGG pathway database. This also increases the biological interpretability of the results. Our model successfully identifies a number of pathways that show differential expression between consecutive purification stages. We also find that changes in the transcriptional response are more strongly associated to the presence of certain compounds, with the remaining contributing to a lesser extent. We discuss the sensitivity of these results to the model parameters that measure the influence of the prior information on the posterior inference.

Endocrinology: advances through omics and related technologies

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

Validation of a genomics-based hypothetical adverse outcome pathway: 2,4-dinitrotoluene perturbs PPAR signaling thus impairing energy metabolism and exercise endurance

2,4-dinitrotoluene (2,4-DNT) is a nitroaromatic used in industrial dyes and explosives manufacturing processes that is found as a contaminant in the environment. Previous studies have implicated antagonism of PPARα signaling as a principal process affected by 2,4-DNT. Here, we test the hypothesis that 2,4-DNT-induced perturbations in PPARα signaling and resultant downstream deficits in energy metabolism, especially from lipids, cause organism-level impacts on exercise endurance. PPAR nuclear activation bioassays demonstrated inhibition of PPARα signaling by 2,4-DNT whereas PPARγ signaling increased. PPARα (-/-) and wild-type (WT) female mice were exposed for 14 days to vehicle or 2,4-DNT (134 mg/kg/day) and performed a forced swim to exhaustion 1 day after the last dose. 2,4-DNT significantly decreased body weights and swim times in WTs, but effects were significantly mitigated in PPARα (-/-) mice. 2,4-DNT decreased transcript expression for genes downstream in the PPARα signaling pathway, principally genes involved in fatty acid transport. Results indicate that PPARγ signaling increased resulting in enhanced cycling of lipid and carbohydrate substrates into glycolytic/gluconeogenic pathways favoring energy production versus storage in 2,4-DNT-exposed WT and PPARα (-/-) mice. PPARα (-/-) mice appear to have compensated for the loss of PPARα by shifting energy metabolism to PPARα-independent pathways resulting in lower sensitivity to 2,4-DNT when compared with WT mice. Our results validate 2,4-DNT-induced perturbation of PPARα signaling as the molecular initiating event for impaired energy metabolism, weight loss, and decreased exercise performance.

Gene expression of fathead minnows (Pimephales promelas) exposed to two types of treated municipal wastewater effluents

Contaminants of emerging concern (CECs) in treated municipal effluents have the potential to adversely impact exposed organisms prompting elevated public concern. Using transcriptomic tools, we investigated changes in gene expression and cellular pathways in the liver of male fathead minnows (Pimephales promelas) exposed to 5% concentrations of full secondary-treated (HTP) or advanced primary-treated (PL) municipal wastewater effluents containing CECs. Gene expression changes were associated with apical end points (plasma vitellogenin and changes in secondary sexual characteristics). Of 32 effluent CECs analyzed, 28 were detected including pharmaceuticals, personal care products, hormones, and industrial compounds. Exposure to both effluents produced significantly higher levels of plasma VTG and changes in secondary sexual characteristics (e.g., ovipositor development). Transcript patterns differed between effluents, with <10% agreement in the detected response (e.g., altered production of transcripts involved in xenobiotic detoxification, oxidative stress, and apoptosis were observed following exposure to both effluents). Exposure to PL effluent caused changes in transcription of genes involved in metabolic pathways (e.g., lipid transport and steroid metabolism). Exposure to HTP effluent affected transcripts involved in signaling pathways (e.g., focal adhesion assembly and extracellular matrix). The results suggest a potential association between some transcriptomic changes and physiological responses following effluent exposure. This study identified responses in pathways not previously implicated in exposure to complex chemical mixtures containing CECs, which are consistent with effluent exposure (e.g., oxidative stress) in addition to other pathway responses specific to the effluent type.

Assessment of chemical mixtures and groundwater effects on Daphnia magna transcriptomics

Small organisms can be used as biomonitoring tools to assess chemicals in the environment. Chemical stressors are especially hard to assess and monitor when present as complex mixtures. Here, fifteen polymerase chain reaction assays targeting Daphnia magna genes were calibrated to responses elicited in D. magna exposed for 24 h to five different doses each of the munitions constituents 2,4,6-trinitrotoluene, 2,4-dinitrotoluene, 2,6-dinitrotoluene, trinitrobenzene, dinitrobenzene, or 1,3,5-trinitro-1,3,5-triazacyclohexane. A piecewise-linear model for log-fold expression changes in gene assays was used to predict response to munitions mixtures and contaminated groundwater under the assumption that chemical effects were additive. The correlations of model predictions with actual expression changes ranged from 0.12 to 0.78 with an average of 0.5. To better understand possible mixture effects, gene expression changes from all treatments were compared using high-density microarrays. Whereas mixtures and groundwater exposures had genes and gene functions in common with single chemical exposures, unique functions were also affected, which was consistent with the nonadditivity of chemical effects in these mixtures. These results suggest that, while gene behavior in response to chemical exposure can be partially predicted based on chemical exposure, estimation of the composition of mixtures from chemical responses is difficult without further understanding of gene behavior in mixtures. Future work will need to examine additive and nonadditive mixture effects using a much greater range of different chemical classes in order to clarify the behavior and predictability of complex mixtures.

Effects of the UV-filter 2-ethyl-hexyl-4-trimethoxycinnamate (EHMC) on expression of genes involved in hormonal pathways in fathead minnows (Pimephales promelas) and link to vitellogenin induction and histology

UV-filters are increasingly used in cosmetics and in the protection of materials against UV-irradiation, and ultimately they reach aquatic systems. The lipophilic UV-filter 2-ethyl-hexyl-4-trimethoxycinnamate (EHMC) belongs to one of the most frequently used UV-filters and accumulates in aquatic animals. Despite its ubiquitous presence in water and biota, very little is known about its potential hormonal effects on aquatic organisms. In our study, we evaluated the effects of measured water concentration of 5.4, 37.5, 244.5 and 394 μg/L EHMC on the expression of genes involved in hormonal pathways in the liver, testis and brain of male and female fathead minnows (Pimephales promelas). We compare the transcription profile with the plasma vitellogenin (VTG) content, secondary sex characteristics, and gonad histology. Transcripts of the androgen receptor (ar) were significantly down-regulated in the liver of females at 37.5, 244.5 μg/L and 394 μg/L EHMC. Additionally, the 3β-hydroxysteroid dehydrogenase (3β-HSD) transcript was significantly decreased in the liver of males at 37.5, 244.5 and 394 μg/L EHMC, and at 244.5 and 394 μg/L EHMC in females. The expressional changes were tissue-specific in most cases, being most significant in the liver. Vitellogenin plasma concentration was significantly increased at 244.5 μg/L EHMC in males. EHMC induced significant histological changes in testes and ovaries at 394 μg/L. Testes displayed a decrease in spermatocytes, and ovaries a decrease in previtellogenic oocytes. The induction of VTG plasma concentration and the histological changes in gonads suggest an estrogenic and/or antiandrogenic activity of EHMC. On the other hand, the gene expression profile shows an antiestrogenic (e.g.: down-regulation of esr1) activity of EHMC. In conclusion, our data demonstrate that EHMC displays low but multiple hormonal activities in fish.

Analysis of common and specific mechanisms of liver function affected by nitrotoluene compounds

Background

Nitrotoluenes are widely used chemical manufacturing and munitions applications. This group of chemicals has been shown to cause a range of effects from anemia and hypercholesterolemia to testicular atrophy. We have examined the molecular and functional effects of five different, but structurally related, nitrotoluenes on using an integrative systems biology approach to gain insight into common and disparate mechanisms underlying effects caused by these chemicals.

Methodology/Principal Findings

Sprague-Dawley female rats were exposed via gavage to one of five concentrations of one of five nitrotoluenes [2,4,6-trinitrotoluene (TNT), 2-amino-4,6-dinitrotoluene (2ADNT) 4-amino-2,6-dinitrotoulene (4ADNT), 2,4-dinitrotoluene (2,4DNT) and 2,6-dinitrotoluene (2,6DNT)] with necropsy and tissue collection at 24 or 48 h. Gene expression profile results correlated well with clinical data and liver histopathology that lead to the concept that hematotoxicity was followed by hepatotoxicity. Overall, 2,4DNT, 2,6DNT and TNT had stronger effects than 2ADNT and 4ADNT. Common functional terms, gene expression patterns, pathways and networks were regulated across all nitrotoluenes. These pathways included NRF2-mediated oxidative stress response, aryl hydrocarbon receptor signaling, LPS/IL-1 mediated inhibition of RXR function, xenobiotic metabolism signaling and metabolism of xenobiotics by cytochrome P450. One biological process common to all compounds, lipid metabolism, was found to be impacted both at the transcriptional and lipid production level.

Conclusions/Significance

A systems biology strategy was used to identify biochemical pathways affected by five nitroaromatic compounds and to integrate data that tie biochemical alterations to pathological changes. An integrative graphical network model was constructed by combining genomic, gene pathway, lipidomic, and physiological endpoint results to better understand mechanisms of liver toxicity and physiological endpoints affected by these compounds.

Investigations of transcript expression in fathead minnow (Pimephales promelas) brain tissue reveal toxicological impacts of RDX exposure

Production, usage and disposal of the munitions constituent (MC) cyclotrimethylenetrinitramine (RDX) has led to environmental releases on military facilities. The chemical attributes of RDX are conducive for leaching to surface water which may put aquatic organisms at risk of exposure. Because RDX has been observed to cause aberrant neuromuscular effects across a wide range of animal phyla, we assessed the effects of RDX on central nervous system (CNS) functions in the representative aquatic ecotoxicological model species, fathead minnow (Pimephales promelas). We developed a fathead minnow brain-tissue cDNA library enriched for transcripts differentially expressed in response to RDX and trinitrotoluene (TNT) exposure. All 4,128 cDNAs were sequenced, quality filtered and assembled yielding 2230 unique sequences and 945 significant blastx matches (E ≤10(-5)). The cDNA library was leveraged to create custom-spotted microarrays for use in transcript expression assays. The impact of RDX on transcript expression in brain tissue was examined in fathead minnows exposed to RDX at 0.625, 2.5, 5, 10mg/L or an acetone-spike control for 10 days. Overt toxicity of RDX in fathead minnow occurred only at the highest exposure concentration resulting in 50% mortality and weight loss. Conversely, Bayesian analysis of microarray data indicated significant changes in transcript expression at concentrations as low as 0.625 mg/L. In total, 154 cDNAs representing 44 unique transcripts were differentially expressed in RDX exposures, the majority of which were validated by reverse transcriptase-quantitative PCR (RT-qPCR). Investigation of molecular pathways, gene ontology (GO) and individual gene functions affected by RDX exposures indicated changes in metabolic processes involved in: oxygen transport, neurological function, calcium binding/signaling, energy metabolism, cell growth/division, oxidative stress and ubiquitination. In total, our study indicated that RDX exposure affected molecular processes critical to CNS function in fathead minnow.

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.

Regulation of iodothyronine deiodinases and sodium iodide symporter mRNA expression by perchlorate in larvae and adult Chinese rare minnow (Gobiocypris rarus)

Perchlorate is a widespread contaminant in the aquatic environment. In the present work, the expressions of deiodinase enzymes (d1, d2, and d3) and sodium iodide symporter (nis) genes were determined after larval and adult rare minnow (Gobiocypris rarus) exposed to 5 and 50 μg/L perchlorate for 21 days. The results showed that deflation of swim bladder development was observed in larvae at 50 μg/L perchlorate treatment. An up-regulation of the d2 and nis mRNA levels were observed in the larve and in brain of adults. Meanwhile the expressions of d3 mRNA levels were significantly down-regulated in the liver. These results indicate the changes in d2, nis, and d3 mRNA expression brings about increased outer-ring deiodination, idodine uptake, and a further decrease of inner-ring deiodination, respectively reflecting auto-regulation of hypothalamic-pituitary-thyroid (HPT) axis in adult after perchlorate exposure. The larval fish development could be affected by perchlorate at environmentally relevant concentrations.

Regulation of thyroid hormone related genes mRNA expression by exogenous T₃ in larvae and adult Chinese rare minnow (Gobiocypris rarus)

In this study, the expression time and profiles of thyroid hormone receptor alpha (trα), type I and II deiodinase enzymes (d1 and d2), transthyretin (ttr), sodium iodide symporter (nis), and thyrotropin-releasing hormone receptor (trhr) genes in Chinese rare minnow (Gobiocypris rarus) were determined using real-time PCR. Meanwhile, the changes of these genes were investigated by exogenous T(3) (3.8 nM) in larvae and adult fish. The retardation of swim bladder development and growth inhibition were observed for larvae, and the transcription of trα, d1, d2, nis, and trhr was significantly down-regulated at the end of exposure (21 d). In adults, a down-regulation of trα, d1, nis, and trhr mRNA levels occurred at 7th or 14th day of exposure, but returned back to their normal levels similar to control at the end of exposure. The down-regulation of gene mRNA expression could serve as a compensatory mechanism for the activation of thyroid system.

A new approach to construct pathway connected networks and its application in dose responsive gene expression profiles of rat liver regulated by 2,4DNT

Abstract

Background

Military and industrial activities have lead to reported release of 2,4-dinitrotoluene (2,4DNT) into soil, groundwater or surface water. It has been reported that 2,4DNT can induce toxic effects on humans and other organisms. However the mechanism of 2,4DNT induced toxicity is still unclear. Although a series of methods for gene network construction have been developed, few instances of applying such technology to generate pathway connected networks have been reported.

Results

Microarray analyses were conducted using liver tissue of rats collected 24h after exposure to a single oral gavage with one of five concentrations of 2,4DNT. We observed a strong dose response of differentially expressed genes after 2,4DNT treatment. The most affected pathways included: long term depression, breast cancer regulation by stathmin1, WNT Signaling; and PI3K signaling pathways. In addition, we propose a new approach to construct pathway connected networks regulated by 2,4DNT. We also observed clear dose response pathway networks regulated by 2,4DNT.

Conclusions

We developed a new method for constructing pathway connected networks. This new method was successfully applied to microarray data from liver tissue of 2,4DNT exposed animals and resulted in the identification of unique dose responsive biomarkers in regards to affected pathways.

Promise and progress in environmental genomics: a status report on the applications of gene expression-based microarray studies in ecologically relevant fish species

The advent of any new technology is typically met with great excitement. So it was a few years ago, when the combination of advances in sequencing technology and the development of microarray technology made measurements of global gene expression in ecologically relevant species possible. Many of the review papers published around that time promised that these new technologies would revolutionize environmental biology as they had revolutionized medicine and related fields. A few years have passed since these technological advancements have been made, and the use of microarray studies in non-model fish species has been adopted in many laboratories internationally. Has the relatively widespread adoption of this technology really revolutionized the fields of environmental biology, including ecotoxicology, aquaculture and ecology, as promised? Or have these studies merely become a novelty and a potential distraction for scientists addressing environmentally relevant questions? In this review, the promises made in early review papers, in particular about the advances that the use of microarrays would enable, are summarized; these claims are compared to the results of recent studies to determine whether the forecasted changes have materialized. Some applications, as discussed in the paper, have been realized and have led to advances in their field, others are still under development.

In vitro gene regulatory networks predict in vivo function of liver

BACKGROUND

Evolution of toxicity testing is predicated upon using in vitro cell based systems to rapidly screen and predict how a chemical might cause toxicity to an organ in vivo. However, the degree to which we can extend in vitro results to in vivo activity and possible mechanisms of action remains to be fully addressed.

RESULTS

Here we use the nitroaromatic 2,4,6-trinitrotoluene (TNT) as a model chemical to compare and determine how we might extrapolate from in vitro data to in vivo effects. We found 341 transcripts differentially expressed in common among in vitro and in vivo assays in response to TNT. The major functional term corresponding to these transcripts was cell cycle. Similarly modulated common pathways were identified between in vitro and in vivo. Furthermore, we uncovered the conserved common transcriptional gene regulatory networks between in vitro and in vivo cellular liver systems that responded to TNT exposure, which mainly contain 2 subnetwork modules: PTTG1 and PIR centered networks. Interestingly, all 7 genes in the PTTG1 module were involved in cell cycle and downregulated by TNT both in vitro and in vivo.

CONCLUSIONS

The results of our investigation of TNT effects on gene expression in liver suggest that gene regulatory networks obtained from an in vitro system can predict in vivo function and mechanisms. Inhibiting PTTG1 and its targeted cell cycle related genes could be key mechanism for TNT induced liver toxicity.

Seasonal reorganization of the xylem transcriptome at different tree ages reveals novel insights into wood formation in Pinus radiata

*Seasonal wood development produces earlywood (EW) and latewood (LW) with distinct properties. The molecular mechanisms controlling EW and LW formation at different tree ages are poorly understood. *Seasonal reorganization of the xylem transcriptome was investigated in Pinus radiata at four tree ages using cDNA microarrays. Transcriptome profiles were compared with seasonal wood variation measured by SilviScan (CSIRO, Clayton, Australia). *The xylem transcriptome was considerably reorganized during seasonal change, and this reorganization showed a maturation-related pattern. The greater reorganization occurred at the transition (30%) and juvenile (21%) stages, but it declined with tree maturity (11-13%). However, this pattern does not correlate well with maturation-related patterns of seasonal wood variation. In total, 319 seasonal-responsive xylem candidate genes were identified. Many transcripts involved in primary and secondary wall biosynthesis were preferentially accumulated in EW and LW, respectively. A large proportion (45-81%) of the candidate genes are preferentially regulated at a single age and their transcript abundance may influence maturation-related patterns of seasonal wood variation. *Seasonal reorganization of the xylem transcriptome is significantly affected by tree age. Physiological changes at the transition stage may contribute to its greater seasonal transcriptome reorganization. Identified stage-preferential xylem transcripts could influence seasonal wood variation at different tree ages.

From raw materials to validated system: the construction of a genomic library and microarray to interpret systemic perturbations in Northern bobwhite

The limited availability of genomic tools and data for nonmodel species impedes computational and systems biology approaches in nonmodel organisms. Here we describe the development, functional annotation, and utilization of genomic tools for the avian wildlife species Northern bobwhite (Colinus virginianus) to determine the molecular impacts of exposure to 2,6-dinitrotoluene (2,6-DNT), a field contaminant of military concern. Massively parallel pyrosequencing of a normalized multitissue library of Northern bobwhite cDNAs yielded 71,384 unique transcripts that were annotated with gene ontology (GO), pathway information, and protein domain analysis. Comparative genome analyses with model organisms revealed functional homologies in 8,825 unique Northern bobwhite genes that are orthologous to 48% of Gallus gallus protein-coding genes. Pathway analysis and GO enrichment of genes differentially expressed in livers of birds exposed for 60 days (d) to 10 and 60 mg/kg/d 2,6-DNT revealed several impacts validated by RT-qPCR including: prostaglandin pathway-mediated inflammation, increased expression of a heme synthesis pathway in response to anemia, and a shift in energy metabolism toward protein catabolism via inhibition of control points for glucose and lipid metabolic pathways, PCK1 and PPARGC1, respectively. This research effort provides the first comprehensive annotated gene library for Northern bobwhite. Transcript expression analysis provided insights into the metabolic perturbations underlying several observed toxicological phenotypes in a 2,6-DNT exposure case study. Furthermore, the systemic impact of dinitrotoluenes on liver function appears conserved across species as PPAR signaling is similarly affected in fathead minnow liver tissue after exposure to 2,4-DNT.

Metatranscriptomic analysis of the response of river biofilms to pharmaceutical products, using anonymous DNA microarrays

Pharmaceutical products are released at low concentrations into aquatic environments following domestic wastewater treatment. Such low concentrations have been shown to induce transcriptional responses in microorganisms, which could have consequences on aquatic ecosystem dynamics. In order to test if these transcriptional responses could also be observed in complex river microbial communities, biofilm reactors were inoculated with water from two rivers of differing trophic statuses and subsequently treated with environmentally relevant doses (ng/liter to microg/liter range) of four pharmaceuticals (erythromycin [ER], gemfibrozil [GM], sulfamethazine [SN], and sulfamethoxazole [SL]). To monitor functional gene expression, we constructed a 9,600-feature anonymous DNA microarray platform onto which cDNA from the biofilms was hybridized. Pharmaceutical treatments induced both positive and negative transcriptional responses from biofilm microorganisms. For instance, ER induced the transcription of several stress, transcription, and replication genes, while GM, a lipid regulator, induced transcriptional responses from several genes involved in lipid metabolism. SN caused shifts in genes involved in energy production and conversion, and SL induced responses from a range of cell membrane and outer envelope genes, which in turn could affect biofilm formation. The results presented here demonstrate for the first time that low concentrations of small molecules can induce transcriptional changes in a complex microbial community. The relevance of these results also demonstrates the usefulness of anonymous DNA microarrays for large-scale metatranscriptomic studies of communities from differing aquatic ecosystems.

Biomarker discovery and transcriptomic responses in Daphnia magna exposed to munitions constituents

Ecotoxicogenomic approaches are emerging as alternative methods in environmental monitoring because they allow insight into pollutant modes of action and help assess the causal agents and potential toxicity beyond the traditional end points of death, growth, and reproduction. Gene expression analysis has shown particular promise for identifying gene expression biomarkers of chemical exposure that can be further used to monitor specific chemical exposures in the environment. We focused on the development of gene expression markers to detect and discriminate between chemical exposures. Using a custom cDNA microarray for Daphnia magna, we identified distinct expression fingerprints in response to exposure at sublethal concentrations of Cu, Zn, Pb, and munitions constituents. Using the results obtained from microarray analysis, we chose a suite of potential biomarkers for each of the specific exposures. The selected potential biomarkers were tested in independent chemical exposures for specificity using quantitative reverse transcription polymerase chain reaction. Six genes were confirmed as differentially regulated bythe selected chemical exposures. Furthermore, each exposure was identified by response of a unique combination (suite) of individual gene expression biomarkers. These results demonstrate the potential for discovery and validation of novel biomarkers of chemical exposures using gene expression analysis, which could have broad applicability in environmental monitoring.

Identification of copper-responsive genes in an early life stage of the fathead minnow Pimephales promelas

While physiological changes associated with copper toxicity have been studied in adult fathead minnow, Pimephales promelas, little is known about the effect of copper on newly hatched larvae. As a result we initiated an investigation on the mechanism of copper toxicity in 24 h post-hatch larvae using gene expression changes to identify responsive genes. Fish were exposed to copper concentrations of 0, 50, 125 and 200 mug/L in a 48 h toxicity test. Total RNA from survivors was used in a differential display assay to screen for differentially expressed gene products. Altogether, 654 copper-responsive differentially expressed bands were collected. Database searches found homology for 261 sequences. One hundred and sixty-one bands were homologous to NCBI genes of known function, of which 69 were individual genes. The most abundant categories of functional genes responding to copper were involved in protein synthesis/translational machinery and contractile proteins. Twenty-one dose-responsive genes were measured for expression changes using real-time quantitative PCR. Differential gene expression was validated for 11 of 13 genes, when a 1.2 times qPCR difference between the copper and control samples was observed. Transcripts identified as titin, cytochrome b, fast muscle specific heavy myosin chain 4, fast muscle troponin I, proteasome 26S subunit and troponin T3a were induced over twofold. Differential display bands identified as 60S ribosomal proteins L27 and L12 were repressed approximately threefold. We conclude that copper exposure affects several cellular pathways in larval fathead minnows with protein synthesis, ribosome structure, and muscle contractile proteins being the most sensitive to this stress.

The outflow tract of the heart in fishes: anatomy, genes and evolution

A large number of congenital heart defects associated with mortality in humans are those that affect the cardiac outflow tract, and this provides a strong imperative to understand its development during embryogenesis. While there is wide phylogenetic variation in adult vertebrate heart morphology, recent work has demonstrated evolutionary conservation in the early processes of cardiogenesis, including that of the outflow tract. This, along with the utility and high reproductive potential of fish species such as Danio rerio, Oryzias latipes etc., suggests that fishes may provide ideal comparative biological models to facilitate a better understanding of this poorly understood region of the heart. In this review, the authors present the current understanding of both phylogeny and ontogeny of the cardiac outflow tract in fishes and examine how new molecular studies are informing the phylogenetic relationships and evolutionary trajectories that have been proposed. The authors also attempt to address some of the issues of nomenclature that confuse this area of research.

The effects of selenium on oxidative stress biomarkers in the freshwater characid fish matrinxã, Brycon cephalus (Günther, 1869) exposed to organophosphate insecticide Folisuper 600 BR (methyl parathion)

Methyl parathion (MP), an organophosphate widely applied in agriculture and aquaculture, induces oxidative stress due to free radical generation and changes in the antioxidant defense system. The antioxidant roles of selenium (Se) were evaluated in Brycon cephalus exposed to 2 mg L(-1) of Folisuper 600 BR (MP commercial formulation - MPc, 600 g L(-1)) for 96 h. Catalase (CAT), glutathione peroxidase (GPx), superoxide dismutase (SOD), glutathione S-transferase (GST), reduced glutathione (GSH) and lipid peroxidation (LPO) levels in the gills, white muscle and liver were evaluated in fish fed on diets containing 0 or 1.5 mg Se kg(-1) for 8 weeks. In fish treated with a Se-free diet, the MPc exposure increased SOD and CAT activities in all tissues. However, the GPx activity decreased in white muscle and gills whereas no alterations were observed in the liver. MPc also increased GST activity in all tissues with a concurrent decrease in GSH levels. LPO values increased in white muscle and gills and did not change in liver after MPc exposure. A Se-supplemented diet reversed these findings, preventing increases in LPO levels and concurrent decreases in GPx activity in gills and white muscle. Similarly, GSH levels were maintained in all tissue after MPc exposure. These results suggest that dietary Se supplementation protects cells against MPc-induced oxidative stress.

Toxicogenomics of water chemistry influence on chronic lead exposure to the fathead minnow (Pimephales promelas)

Establishment of water quality criteria (WQC), intended to protect aquatic life, continues to rely principally on water hardness (i.e. Ca(2+)) for lead (Pb) despite growing evidence that other chemical parameters also strongly influence toxicity. To more clearly define the water chemistry parameters mediating Pb toxicity, we evaluated the effects of hardness as CaSO(4) and dissolved organic carbon (DOC) as humic acid during chronic (150 days) exposures to the fathead minnow. Measured Pb concentrations ranged from 157+/-5 nM (33+/-1 microg/L) Pb in base water to 177+/-7 (37+/-1 microg/L) and 187+/-7 nM (39+/-1 microg/L) Pb in CaSO(4)- or HA-supplemented water, respectively. Fish were collected at 2, 4, 10, 30, 63, 90 and 150 days of exposure. Traditional toxicological endpoints were examined alongside gene transcription analyses to help clarify the underlying mechanisms of Pb toxicity and to identify candidate molecular markers that might ultimately serve as robust indicators of exposure and effect. Addition of CaSO(4) did not prevent whole body Pb accumulation whereas DOC afforded strong protection (about half the amount accumulated by fish in base water) suggesting that current, hardness-based WQC are likely inaccurate for predicting chronic Pb effects in aquatic systems. Custom-made microarrays were co-hybridized with base water samples+/-Pb up to the 30 days time point. Quantitative PCR was employed to verify gene transcription responses and to extend analysis to the CaSO(4) and HA treatments and the 150 days time point. Identification of four genes by microarray analysis revealed clear Pb-induced responses over time: glucose-6-phosphate dehydrogenase, glutathione-S-transferase, ferritin and beta-globin. Results obtained by qPCR were in strong agreement with microarray data by regression analysis (r(2)=0.82, slope=1.28). The associated pathways implicated herein for these genes provide further evidence supporting roles for anemia and neurological disorders in chronic Pb toxicity. Effects of water chemistry on Pb accumulation and gene transcription responses were in close parallel, though alterations in ionoregulatory and morphological endpoints were not observed. Whereas DOC was protective against Pb accumulation and mRNA expression changes, Ca(2+) was not. Additionally, several hypothesis-driven genes (ECaC, DMT-1, and ALA-D) were examined by qPCR but revealed either no change or small Pb-induced responses lacking any clear influence attributable to water chemistry. These findings should help pave the way toward development of a new chronic Pb BLM and a Pb-responsive gene transcript profile for fathead minnows, both of which would greatly aid future environmental monitoring and regulatory strategies for Pb.

Gene expression profiles in rainbow trout, Onchorynchus mykiss, exposed to a simple chemical mixture

Among proposed uses for microarrays in environmental toxiciology is the identification of key contributors to toxicity within a mixture. However, it remains uncertain whether the transcriptomic profiles resulting from exposure to a mixture have patterns of altered gene expression that contain identifiable contributions from each toxicant component. We exposed isogenic rainbow trout Onchorynchus mykiss, to sublethal levels of ethynylestradiol, 2,2,4,4-tetrabromodiphenyl ether, and chromium VI or to a mixture of all three toxicants Fluorescently labeled complementary DNA (cDNA) were generated and hybridized against a commercially available Salmonid array spotted with 16,000 cDNAs. Data were analyzed using analysis of variance (p<0.05) with a Benjamani-Hochberg multiple test correction (Genespring [Agilent] software package) to identify up and downregulated genes. Gene clustering patterns that can be used as "expression signatures" were determined using hierarchical cluster analysis. The gene ontology terms associated with significantly altered genes were also used to identify functional groups that were associated with toxicant exposure. Cross-ontological analytics approach was used to assign functional annotations to genes with "unknown" function. Our analysis indicates that transcriptomic profiles resulting from the mixture exposure resemble those of the individual contaminant exposures, but are not a simple additive list. However, patterns of altered genes representative of each component of the mixture are clearly discernible, and the functional classes of genes altered represent the individual components of the mixture. These findings indicate that the use of microarrays to identify transcriptomic profiles may aid in the identification of key stressors within a chemical mixture, ultimately improving environmental assessment.

Toxicogenomic analysis provides new insights into molecular mechanisms of the sublethal toxicity of 2,4,6-trinitrotoluene in Eisenia fetida

Xenobiotics such as explosives and pesticides released into the environment can have lethal and sublethal impacts on soil organisms such as earthworms with potential subsequent impacts at highertrophic levels. To better understand the molecular toxicological mechanisms of 2,4,6-trinitrotoluene (TNT), a commonly used explosive, in Eisenia fetida, earthworms were exposed to a gradient of TNT-spiked soils for 28 days and impacts on gene expression were examined using a 4032 cDNA microarray. Reproduction was increased at low doses of TNT, whereas high doses of TNT reduced juvenile production. On the basis of reproduction responses to TNT, four treatments, that is, control, 2, 10.6, and 38.7 mg/kg, were selected for gene expression studies in a balanced interwoven loop design microarray experiment in which the expression of 311 transcripts was significantly affected. Reverse-transcription quantitative polymerase chain reaction (RT-QPCR) data on 68 selected differentially and nondifferentially expressed transcripts showed a significant correlation with microarray results. The expression of genes involved in multiple biological processes was altered, including muscle contraction, neuronal signaling and growth, ubiquitinylation, fibrinolysis and coagulation, iron and calcium homeostasis, oxygen transport, and immunity. Chitinase activity assays confirmed down-regulation of chitinase genes as indicated by array and RT-QPCR data. An acute toxicity test provided evidence that dermal contact with TNT can cause bleeding, inflammation, and constriction, which may be explained by gene expression results. Sublethal doses of TNT affected the nervous system, caused blood disorders similar to methemoglobinemia, and weakened immunity in E. fetida.

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.