Natàlia Garcia-Reyero

Systems biology: Leading the revolution in ecotoxicology

The rapid development of new technologies such as transcriptomics, proteomics, and metabolomics (Omics) are changing the way ecotoxicology is practiced. The data deluge has begun with genomes of over 65 different aquatic species that are currently being sequenced, and many times that number with at least some level of transcriptome sequencing. Integrating these top-down methodologies is an essential task in the field of systems biology. Systems biology is a biology-based interdisciplinary field that focuses on complex interactions in biological systems, with the intent to model and discover emergent properties of the system. Recent studies demonstrate that Omics technologies provide valuable insight into ecotoxicity, both in laboratory exposures with model organisms and with animals exposed in the field. However, these approaches require a context of the whole animal and population to be relevant. Powerful approaches using reverse engineering to determine interacting networks of genes, proteins, or biochemical reactions are uncovering unique responses to toxicants. Modeling efforts in aquatic animals are evolving to interrelate the interacting networks of a system and the flow of information linking these elements. Just as is happening in medicine, systems biology approaches that allow the integration of many different scales of interaction and information are already driving a revolution in understanding the impacts of pollutants on aquatic systems.

Effects of a short-term exposure to the fungicide prochloraz on endocrine function and gene expression in female fathead minnows (Pimephales promelas).

Prochloraz is a fungicide known to cause endocrine disruption through effects on the hypothalamic-pituitary-gonadal (HPG) axis. To determine the short-term impacts of prochloraz on gene expression and steroid production, adult female fathead minnows (Pimephales promelas) were exposed to the chemical (0 or 300 μg/L) for a time-course of 6, 12 and 24 h. Consistent with inhibition of cytochrome P450 17α-hydroxylase/17,20-lyase (CYP17) and aromatase (CYP19), known molecular targets of prochloraz, plasma 17β-estradiol (E2) was reduced within 6 h. Ex vivo E2 production was significantly reduced at all time-points, while ex vivo testosterone (T) production remained unchanged. Consistent with the decrease in E2 levels, plasma concentrations of the estrogen-responsive protein vitellogenin were significantly reduced at 24 h. Genes coding for CYP19, CYP17, and steroidogenic acute regulatory protein were up-regulated in a compensatory manner in ovaries of the prochloraz-treated fish. In addition to targeted quantitative real-time polymerase chain reaction analyses, a 15k feature fathead minnow microarray was used to determine gene expression profiles in ovaries. From time-point to time-point, the microarray results showed a relatively rapid change in the differentially expressed gene (DEG) profiles associated with the chemical exposure. Functional analysis of the DEGs indicated changes in expression of genes associated with cofactor and coenzyme binding (GO:0048037 and 0050662), fatty acid binding (GO:0005504) and organelle organization and biogenesis (GO:0006996). Overall, the results from this study are consistent with compensation of the fish HPG axis to inhibition of steroidogenesis by prochloraz, and provide further insights into relatively rapid, system-wide, effects of a model chemical stressor on fish.

Influence of ovarian stage on transcript profiles in fathead minnow (Pimephales promelas) ovary tissue

Interpretation of toxicogenomic experiments conducted with ovary tissue from asynchronous-spawning small fish species is complicated by background variation in the relative abundance and proportion of follicles at different stages within the ovary tissue sample. This study employed both real-time quantitative polymerase chain reaction and a 15,000 gene oligonucleotide microarray to examine variation in the fathead minnow (Pimephales promelas) ovarian transcriptional profile as a function of quantitative and qualitative differences in ovarian histology. The objectives were to provide data that could potentially aid interpretation of future toxicogenomics experiments, identify putative stage-related transcriptional markers, and generate insights into basic biological regulation of asynchronous oocyte development. Multiple lines of evidence from the present study indicate that variation in the transcriptional profile is primarily dependent on the relative abundance of previtellogenic versus vitellogenic follicles in the ovary. Due to the relatively small proportions of mature ovulated follicles or atretic follicles in the overall follicle population, few potential transcriptional markers of maturation, ovulation, or atresia could be identified. However, among the 460 differentially expressed genes identified in the present study, several targets, including HtrA serine peptidase 3 (htra3), tissue inhibitor of metalloproteinase 3 (timp3), aquaporin 8 (aqp8), transgelin 2 like (tagln2), Nedd4 family interacting protein 2 (ndfip2), chemokine ligand 12a (cxcl12a), midkine-related growth factor (mdka), and jagged 1b (jag 1b) exhibited responses and functional properties that support them as candidate molecular markers of significant shift in gross ovarian stage. Genes associated with a diversity of functions including cellular development, morphogenesis, coated vesicle transport, sexual reproduction, and neuron development, among others, were statistically enriched within the list of 460 genes differentially expressed among different ovarian classes. Overall, results of this study provide insights into background variation in ovary transcript profiles that should aid and enhance the interpretation of toxicogenomic data generated in experiments conducted with small, asynchronous-spawning fish species.

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.

Conserved toxic responses across divergent phylogenetic lineages: a meta-analysis of the neurotoxic effects of RDX among multiple species using toxicogenomics

At military training sites, a variety of pollutants such as hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), may contaminate the area originating from used munitions. Studies investigating the mechanism of toxicity of RDX have shown that it affects the central nervous system causing seizures in humans and animals. Environmental pollutants such as RDX have the potential to affect many different species, therefore it is important to establish how phylogenetically distant species may respond to these types of emerging pollutants. In this paper, we have used a transcriptional network approach to compare and contrast the neurotoxic effects of RDX among five phylogenetically disparate species: rat (Sprague-Dawley), Northern bobwhite quail (Colinus virginianus), fathead minnow (Pimephales promelas), earthworm (Eisenia fetida), and coral (Acropora formosa). Pathway enrichment analysis indicated a conservation of RDX impacts on pathways related to neuronal function in rat, Northern bobwhite quail, fathead minnows and earthworm, but not in coral. As evolutionary distance increased common responses decreased with impacts on energy and metabolism dominating effects in coral. A neurotransmission related transcriptional network based on whole rat brain responses to RDX exposure was used to identify functionally related modules of genes, components of which were conserved across species depending upon evolutionary distance. Overall, the meta-analysis using genomic data of the effects of RDX on several species suggested a common and conserved mode of action of the chemical throughout phylogenetically remote organisms.

I. Effects of a dopamine receptor antagonist on fathead minnow, Pimephales promelas, reproduction

Neurotransmitters such as dopamine play an important role in regulating fish reproduction. However, the potential for neuroendocrine active chemicals to disrupt fish reproduction has not been well studied, despite emerging evidence of their discharge into aquatic environments. This study is the first to apply the fathead minnow 21 d reproduction assay developed for the US Endocrine Disruptor Screening Program to evaluate the reproductive toxicity of a model neuroendocrine active chemical, the dopamine 2 receptor antagonist, haloperidol. Continuous exposure to up to 20 imcrog haloperidol/L had no significant effects on fathead minnow fecundity, secondary sex characteristics, gonad histology, or plasma steroid and vitellogenin concentrations. The only significant effect observed was an increase in gonadotropin-releasing hormone (cGnRH) transcripts in the male brain. Results suggest that non-lethal concentrations of haloperidol do not directly impair fish reproduction. Potential effects of haloperidol on reproductive behaviors and gene expression were examined in a companion study.

Fathead minnow steroidogenesis: in silico analyses reveals tradeoffs between nominal target efficacy and robustness to cross-talk

BackgroundInterpreting proteomic and genomic data is a major challenge in predictive ecotoxicology that can be addressed by a systems biology approach. Mathematical modeling provides an organizational platform to consolidate protein dynamics with possible genomic regulation. Here, a model of ovarian steroidogenesis in the fathead minnow, Pimephales promelas, (FHM) is developed to evaluate possible transcriptional regulation of steroid production observed in microarray studies.ResultsThe model was developed from literature sources, integrating key signaling components (G-protein and PKA activation) with their ensuing effect on steroid production. The model properly predicted trajectory behavior of estradiol and testosterone when fish were exposed to fadrozole, a specific aromatase inhibitor, but failed to predict the steroid hormone behavior occurring one week post-exposure as well as the increase in steroid levels when the stressor was removed. In vivo microarray data implicated three modes of regulation which may account for over-production of steroids during a depuration phase (when the stressor is removed): P450 enzyme up-regulation, inhibin down-regulation, and luteinizing hormone receptor up-regulation. Simulation studies and sensitivity analysis were used to evaluate each case as possible source of compensation to endocrine stress.ConclusionsSimulation studies of the testosterone and estradiol response to regulation observed in microarray data supported the hypothesis that the FHM steroidogenesis network compensated for endocrine stress by modulating the sensitivity of the ovarian network to global cues coming from the hypothalamus and pituitary. Model predictions of luteinizing hormone receptor regulation were consistent with depuration and in vitro data. These results challenge the traditional approach to network elucidation in systems biology. Generally, the most sensitive interactions in a network are targeted for further elucidation but microarray evidence shows that homeostatic regulation of the steroidogenic network is likely maintained by a mildly sensitive interaction. We hypothesize that effective network elucidation must consider both the sensitivity of the target as well as the targets robustness to biological noise (in this case, to cross-talk) when identifying possible points of regulation.

II: Effects of a dopamine receptor antagonist on fathead minnow dominance behavior and ovarian gene expression in the fathead minnow and zebrafish.

Neurotransmitters such as dopamine play an important role in reproductive behaviors and signaling. Neuroendocrine-active chemicals in the environment have potential to interfere with and/or alter these processes. A companion study with the dopamine 2 receptor antagonist, haloperidol, found no evidence of a direct effect of the chemical on fish reproduction. This study considered haloperidols potential effects on behavior and ovarian gene expression. Male fathead minnows exposed to 50 microg haloperidol/L for 96 h were found to be significantly more dominant than control males. In terms of molecular signaling, investigated using oligonucleotide microarrays, there was little similarity in the identity and functions of genes differentially expressed in the ovaries of fathead minnows (Pimephales promelas) versus zebrafish (Danio rerio) exposed under the same conditions. Results suggest that non-lethal concentrations of haloperidol do not induce ovarian molecular responses that could serve as biomarkers of exposure to D2R antagonists, but may impact behavior.

Use of Genetic Manipulation for Evaluating and Understanding Adverse Outcome Pathways

Innovations in biology have brought forth a new era of genetic manipulation ranging from the creation of molecular scissors for targeted single-nucleotide alterations, to a simultaneous inactivation of 62 genes in pig embryos to “humanize” transplant tissue. Genetic engineering advances allow for novel testing paradigms to understand chemical interactions and information flow in biological systems. Emerging platforms may provide mechanistic knowledge of chemical stressor interactions in biological systems to facilitate the development of alternative testing methods, as well as prioritize higher tier toxicity testing for risk assessment. This chapter will discuss recent advances in genetic manipulation and describe how these techniques improve our understanding of toxicity across multiple biological scales. These efforts will ultimately aid in validation of Adverse Outcome Pathway (AOP) key event relationships for ecological risk assessment.