A computational model of the hypothalamic: pituitary: gonadal axis in female fathead minnows (Pimephales promelas) exposed to 17α-ethynylestradiol and 17β-trenbolone

PBK-TD modelling of the gonadotropic axis: Case study with two azole fungicides in female zebrafish

Endocrine disruptors (EDs) can disrupt the gonadotropic axis, which consists of the Hypothalamus-Pituitary-Gonads (HPG), notably by altering aromatase (cyp19a), a key enzyme regulating the endocrine system and reproductive function in fish. The effects of EDs can be predicted by integrating both toxicokinetic (TK) and toxicodynamic (TD) processes in order to relate adverse outcomes to external exposures. In this study, we developed a physiologically based kinetic-toxicodynamic model to simulate the disruption of the HPG axis (PBK-TD, hereafter named PBK-HPG) in female zebrafish exposed to either of two aromatase inhibitors, imazalil or prochloraz. The model was calibrated using Bayesian methods and supported by novel experimental data, including measurements of vitellogenin, 17β-estradiol, and 11-ketotestosterone levels, along with in vivo monitoring of the cyp19a1a gene in transgenic cyp19a1a-GFP ebrafish. Seamless integration of a PBK model within a TD model of the HPG-axis, provided the link between external exposure and internal levels of imazalil and prochloraz in key organs, allowing for mechanistic predictions of their inhibitory effects on gonadal aromatase. Our PBK-HPG model accurately predicted both baseline homeostasis and the effects of aromatase inhibition, with all endocrine endpoints including reproductive disruption, i.e., decreased egg production, falling within a twofold range of both experimental and literature data. Therefore, our PBK-HPG model could further support the development of a mechanistic qAOP with TK considerations. The model offers significant potential for improving environmental risk assessments of EDs and possibly other stressors across species.

Putative adverse outcome pathway development based on physiological responses of female fathead minnows to model estrogen versus androgen receptor agonists

Several adverse outcome pathways (AOPs) have linked molecular initiating events like aromatase inhibition, androgen receptor (AR) agonism, and estrogen receptor (ER) antagonism to reproductive impairment in adult fish. Estrogen receptor agonists can also cause adverse reproductive effects, however, the early key events (KEs) in an AOP leading to this are mostly unknown. The primary aim of this study was to develop hypotheses regarding the potential mechanisms through which exposure to ER agonists might lead to reproductive impairment in female fish. Mature fathead minnows were exposed to 1 or 10 ng 17α-ethynylestradiol (EE2)/L or 10 or 100 µg bisphenol A (BPA)/L for 14 d. The response to EE2 and BPA was contrasted with the effects of 500 ng/L of 17β-trenbolone (TRB), an AR agonist, as well as TRB combined with the low and high concentrations of EE2 or BPA tested individually. Exposure to 10 ng EE2/L, 100 µg BPA/L, TRB, or the various mixtures with TRB caused significant decreases in plasma concentrations of 17β-estradiol. Exposure to TRB alone caused a significant reduction in plasma vitellogenin (VTG), but VTG was unaffected or even increased in females exposed to EE2 or BPA alone or, in most cases, in mixtures with TRB. Over the course of the 14-d exposure, the only treatments that clearly did not affect egg production were 1 ng EE2/L and 10 µg BPA/L. Based on these results and knowledge of hypothalamic-pituitary-gonadal axis function, we hypothesize an AOP whereby decreased production of maturation-inducing steroid leading to impaired oocyte maturation and ovulation, possibly due to negative feedback or direct inhibitory effects of membrane ER activation, could be responsible for causing adverse reproductive impacts in female fish exposed to ER agonists.

Transcriptomic evidence of luteinizing hormone-releasing hormone agonist (LHRH-A) regulation on lipid metabolism in grass carp (Ctenopharyngodon idella)

In this study, RNA sequencing was used to identify the hepatic gene expression profile in grass carp associated with luteinizing hormone-releasing hormone agonist (LHRH-A) treatment. A total of 93,912,172 reads were generated by HiSeq 4000 sequencing platform. After filtering, 83,450,860 clean reads were mapped to the reference genome. By calculating the FPKM of genes, 1475 differentially expressed genes were identified. PPAR signaling pathway was enriched with upregulated genes in LHRH-A injection group showing the regulation of the lipid metabolism by LHRH-A. The expression of eight key genes in PPAR signaling pathway was confirmed by qPCR and the results suggested that ACSL4A, ACSL4B, ANGPTL4, LPL, RXRBA and SLC27A1B were significantly stimulated by LHRH-A injection. This investigation provides the evidence that LHRH-A could play a role in lipid metabolism.

Incorporating Suborganismal Processes into Dynamic Energy Budget Models for Ecological Risk Assessment

A working group at the National Institute for Mathematical and Biological Synthesis (NIMBioS) explored the feasibility of integrating 2 complementary approaches relevant to ecological risk assessment. Adverse outcome pathway (AOP) models provide "bottom-up" mechanisms to predict specific toxicological effects that could affect an individual's ability to grow, reproduce, and/or survive from a molecular initiating event. Dynamic energy budget (DEB) models offer a "top-down" approach that reverse engineers stressor effects on growth, reproduction, and/or survival into modular characterizations related to the acquisition and processing of energy resources. Thus, AOP models quantify linkages between measurable molecular, cellular, or organ-level events, but they do not offer an explicit route to integratively characterize stressor effects at higher levels of organization. While DEB models provide the inherent basis to link effects on individuals to those at the population and ecosystem levels, their use of abstract variables obscures mechanistic connections to suborganismal biology. To take advantage of both approaches, we developed a conceptual model to link DEB and AOP models by interpreting AOP key events as measures of damage-inducing processes affecting DEB variables and rates. We report on the type and structure of data that are generated for AOP models that may also be useful for DEB models. We also report on case studies under development that merge information collected for AOPs with DEB models and highlight some of the challenges. Finally, we discuss how the linkage of these 2 approaches can improve ecological risk assessment, with possibilities for progress in predicting population responses to toxicant exposures within realistic environments. Integr Environ Assess Manag 2018;14:615-624. © 2018 SETAC.

Statistical relationship between metabolic decomposition and chemical uptake predicts bioconcentration factor data for diverse chemical exposures

BACKGROUND

A challenge of in vitro to in vivo extrapolation (IVIVE) is to predict the physical state of organisms exposed to chemicals in the environment from in vitro exposure assay data. Although toxicokinetic modeling approaches promise to bridge in vitro screening data with in vivo effects, they are often encumbered by a need for redesign or re-parameterization when applied to different tissues or chemicals.

RESULTS

We demonstrate a parameterization of reverse toxicokinetic (rTK) models developed for the adult zebrafish (Danio rerio) based upon particle swarm optimizations (PSO) of the chemical uptake and degradation rates that predict bioconcentration factors (BCF) for a broad range of chemicals. PSO reveals a relationship between chemical uptake and decomposition parameter values that predicts chemical-specific BCF values with moderate statistical agreement to a limited yet diverse chemical dataset, and all without a need to retrain the model to new data.

CONCLUSIONS

The presented model requires only the octanol-water partitioning ratio to predict BCFs to a fidelity consistent with existing QSAR models. This success begs re-evaluation of the modeling assumptions; specifically, it suggests that chemical uptake into arterial blood may be limited by transport across gill membranes (diffusion) rather than by counter-current flow between gill lamellae (convection). Therefore, more detailed molecular modeling of aquatic respiration may further improve predictive accuracy of the rTK approach.

A critical review of the environmental occurrence and potential effects in aquatic vertebrates of the potent androgen receptor agonist 17β-trenbolone

Trenbolone acetate is widely used in some parts of the world for its desirable anabolic effects on livestock. Several metabolites of the acetate, including 17β-trenbolone, have been detected at low nanograms per liter concentrations in surface waters associated with animal feedlots. The 17β-trenbolone isomer can affect androgen receptor signaling pathways in various vertebrate species at comparatively low concentrations/doses. The present article provides a comprehensive review and synthesis of the existing literature concerning exposure to and biological effects of 17β-trenbolone, with an emphasis on potential risks to aquatic animals. In vitro studies indicate that, although 17β-trenbolone can activate several nuclear hormone receptors, its highest affinity is for the androgen receptor in all vertebrate taxa examined, including fish. Exposure of fish to nanograms per liter water concentrations of 17β-trenbolone can cause changes in endocrine function in the short term, and adverse apical effects in longer exposures during development and reproduction. Impacts on endocrine function typically are indicative of inappropriate androgen receptor signaling, such as changes in sex steroid metabolism, impacts on gonadal stage, and masculinization of females. Exposure of fish to 17β-trenbolone during sexual differentiation in early development can greatly skew sex ratios, whereas adult exposures can adversely impact fertility and fecundity. To fully assess ecosystem-level risks, additional research is warranted to address uncertainties as to the degree/breadth of environmental exposures and potential population-level effects of 17β-trenbolone in sensitive species. Environ Toxicol Chem 2018;37:2064-2078. Published 2018 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

In Silico Computational Transcriptomics Reveals Novel Endocrine Disruptors in Largemouth Bass ( Micropterus salmoides)

In recent years, decreases in fish populations have been attributed, in part, to the effect of environmental chemicals on ovarian development. To understand the underlying molecular events we developed a dynamic model of ovary development linking gene transcription to key physiological end points, such as gonadosomatic index (GSI), plasma levels of estradiol (E2) and vitellogenin (VTG), in largemouth bass ( Micropterus salmoides). We were able to identify specific clusters of genes, which are affected at different stages of ovarian development. A subnetwork was identified that closely linked gene expression and physiological end points and by interrogating the Comparative Toxicogenomic Database (CTD), quercetin and tretinoin (ATRA) were identified as two potential candidates that may perturb this system. Predictions were validated by investigation of reproductive associated transcripts using qPCR in ovary and in the liver of both male and female largemouth bass treated after a single injection of quercetin and tretinoin (10 and 100 μg/kg). Both compounds were found to significantly alter the expression of some of these genes. Our findings support the use of omics and online repositories for identification of novel, yet untested, compounds. This is the first study of a dynamic model that links gene expression patterns across stages of ovarian development.

Physiological modes of action across species and toxicants: the key to predictive ecotoxicology

As ecotoxicologists we strive for a better understanding of how chemicals affect our environment. Humanity needs tools to identify those combinations of man-made chemicals and organisms most likely to cause problems. In other words: which of the millions of species are at risk from pollution? And which of the tens of thousands of chemicals contribute most to the risk? We identified our poor knowledge on physiological modes of action (how a chemical affects the energy allocation in an organism), and how they vary across species and toxicants, as a major knowledge gap. We also find that the key to predictive ecotoxicology is the systematic, rigorous characterization of physiological modes of action because that will enable more powerful in vitro to in vivo toxicity extrapolation and in silico ecotoxicology. In the near future, we expect a step change in our ability to study physiological modes of action by improved, and partially automated, experimental methods. Once we have populated the matrix of species and toxicants with sufficient physiological mode of action data we can look for patterns, and from those patterns infer general rules, theory and models.

Quantitative Adverse Outcome Pathways and Their Application to Predictive Toxicology

A quantitative adverse outcome pathway (qAOP) consists of one or more biologically based, computational models describing key event relationships linking a molecular initiating event (MIE) to an adverse outcome. A qAOP provides quantitative, dose-response, and time-course predictions that can support regulatory decision-making. Herein we describe several facets of qAOPs, including (a) motivation for development, (b) technical considerations, (c) evaluation of confidence, and (d) potential applications. The qAOP used as an illustrative example for these points describes the linkage between inhibition of cytochrome P450 19A aromatase (the MIE) and population-level decreases in the fathead minnow (FHM; Pimephales promelas). The qAOP consists of three linked computational models for the following: (a) the hypothalamic-pitutitary-gonadal axis in female FHMs, where aromatase inhibition decreases the conversion of testosterone to 17β-estradiol (E2), thereby reducing E2-dependent vitellogenin (VTG; egg yolk protein precursor) synthesis, (b) VTG-dependent egg development and spawning (fecundity), and (c) fecundity-dependent population trajectory. While development of the example qAOP was based on experiments with FHMs exposed to the aromatase inhibitor fadrozole, we also show how a toxic equivalence (TEQ) calculation allows use of the qAOP to predict effects of another, untested aromatase inhibitor, iprodione. While qAOP development can be resource-intensive, the quantitative predictions obtained, and TEQ-based application to multiple chemicals, may be sufficient to justify the cost for some applications in regulatory decision-making.

Physiological fidelity or model parsimony? The relative performance of reverse-toxicokinetic modeling approaches

Background

Physiologically-based toxicokinetic (PBTK) models are often developed to facilitate in vitro to in vivo extrapolation (IVIVE) using a top-down, compartmental approach, favoring architectural simplicity over physiological fidelity despite the lack of general guidelines relating model design to dynamical predictions. Here we explore the impact of design choice (high vs. low fidelity) on chemical distribution throughout an animal’s organ system.

Results

We contrast transient dynamics and steady states of three previously proposed PBTK models of varying complexity in response to chemical exposure. The steady states for each model were determined analytically to predict exposure conditions from tissue measurements. Steady state whole-body concentrations differ between models, despite identical environmental conditions, which originates from varying levels of physiological fidelity captured by the models. These differences affect the relative predictive accuracy of the inverted models used in exposure reconstruction to link effects-based exposure data with whole-organism response thresholds obtained from in vitro assay measurements.

Conclusions

Our results demonstrate how disregarding physiological fideltiy in favor of simpler models affects the internal dynamics and steady state estimates for chemical accumulation within tissues, which, in turn, poses significant challenges for the exposure reconstruction efforts that underlie many IVIVE methods. Developing standardized systems-level models for ecological organisms would not only ensure predictive consistency among future modeling studies, but also ensure pragmatic extrapolation of in vivo effects from in vitro data or modeling exposure-response relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0407-3) contains supplementary material, which is available to authorized users.

Toxicokinetic models and related tools in environmental risk assessment of chemicals

Environmental risk assessment of chemicals for the protection of ecosystems integrity is a key regulatory and scientific research field which is undergoing constant development in modelling approaches and harmonisation with human risk assessment. This review focuses on state-of-the-art toxicokinetic tools and models that have been applied to terrestrial and aquatic species relevant to environmental risk assessment of chemicals. Both empirical and mechanistic toxicokinetic models are discussed using the results of extensive literature searches together with tools and software for their calibration and an overview of applications in environmental risk assessment. These include simple tools such as one-compartment models, multi-compartment models to physiologically-based toxicokinetic (PBTK) models, mostly available for aquatic species such as fish species and a number of chemical classes including plant protection products, metals, persistent organic pollutants, nanoparticles. Data gaps and further research needs are highlighted.

How Adverse Outcome Pathways Can Aid the Development and Use of Computational Prediction Models for Regulatory Toxicology

Efforts are underway to transform regulatory toxicology and chemical safety assessment from a largely empirical science based on direct observation of apical toxicity outcomes in whole organism toxicity tests to a predictive one in which outcomes and risk are inferred from accumulated mechanistic understanding. The adverse outcome pathway (AOP) framework provides a systematic approach for organizing knowledge that may support such inference. Likewise, computational models of biological systems at various scales provide another means and platform to integrate current biological understanding to facilitate inference and extrapolation. We argue that the systematic organization of knowledge into AOP frameworks can inform and help direct the design and development of computational prediction models that can further enhance the utility of mechanistic and in silico data for chemical safety assessment. This concept was explored as part of a workshop on AOP-Informed Predictive Modeling Approaches for Regulatory Toxicology held September 24-25, 2015. Examples of AOP-informed model development and its application to the assessment of chemicals for skin sensitization and multiple modes of endocrine disruption are provided. The role of problem formulation, not only as a critical phase of risk assessment, but also as guide for both AOP and complementary model development is described. Finally, a proposal for actively engaging the modeling community in AOP-informed computational model development is made. The contents serve as a vision for how AOPs can be leveraged to facilitate development of computational prediction models needed to support the next generation of chemical safety assessment.

Perspectives for integrating human and environmental exposure assessments

Integrated Risk Assessment (IRA) has been defined by the EU FP7 HEROIC Coordination action as "the mutual exploitation of Environmental Risk Assessment for Human Health Risk Assessment and vice versa in order to coherently and more efficiently characterize an overall risk to humans and the environment for better informing the risk analysis process" (Wilks et al., 2015). Since exposure assessment and hazard characterization are the pillars of risk assessment, integrating Environmental Exposure assessment (EEA) and Human Exposure assessment (HEA) is a major component of an IRA framework. EEA and HEA typically pursue different targets, protection goals and timeframe. However, human and wildlife species also share the same environment and they similarly inhale air and ingest water and food through often similar overlapping pathways of exposure. Fate models used in EEA and HEA to predict the chemicals distribution among physical and biological media are essentially based on common properties of chemicals, and internal concentration estimations are largely based on inter-species (i.e. biota-to-human) extrapolations. Also, both EEA and HEA are challenged by increasing scientific complexity and resources constraints. Altogether, these points create the need for a better exploitation of all currently existing data, experimental approaches and modeling tools and it is assumed that a more integrated approach of both EEA and HEA may be part of the solution. Based on the outcome of an Expert Workshop on Extrapolations in Integrated Exposure Assessment organized by the HEROIC project in January 2014, this paper identifies perspectives and recommendations to better harmonize and extrapolate exposure assessment data, models and methods between Human Health and Environmental Risk Assessments to support the further development and promotion of the concept of IRA. Ultimately, these recommendations may feed into guidance showing when and how to apply IRA in the regulatory decision-making process for chemicals.

The pros and cons of ecological risk assessment based on data from different levels of biological organization

Ecological risk assessment (ERA) is the process used to evaluate the safety of manufactured chemicals to the environment. Here we review the pros and cons of ERA across levels of biological organization, including suborganismal (e.g., biomarkers), individual, population, community, ecosystem and landscapes levels. Our review revealed that level of biological organization is often related negatively with ease at assessing cause-effect relationships, ease of high-throughput screening of large numbers of chemicals (it is especially easier for suborganismal endpoints), and uncertainty of the ERA because low levels of biological organization tend to have a large distance between their measurement (what is quantified) and assessment endpoints (what is to be protected). In contrast, level of biological organization is often related positively with sensitivity to important negative and positive feedbacks and context dependencies within biological systems, and ease at capturing recovery from adverse contaminant effects. Some endpoints did not show obvious trends across levels of biological organization, such as the use of vertebrate animals in chemical testing and ease at screening large numbers of species, and other factors lacked sufficient data across levels of biological organization, such as repeatability, variability, cost per study and cost per species of effects assessment, the latter of which might be a more defensible way to compare costs of ERAs than cost per study. To compensate for weaknesses of ERA at any particular level of biological organization, we also review mathematical modeling approaches commonly used to extrapolate effects across levels of organization. Finally, we provide recommendations for next generation ERA, submitting that if there is an ideal level of biological organization to conduct ERA, it will only emerge if ERA is approached simultaneously from the bottom of biological organization up as well as from the top down, all while employing mathematical modeling approaches where possible to enhance ERA. Because top-down ERA is unconventional, we also offer some suggestions for how it might be implemented efficaciously. We hope this review helps researchers in the field of ERA fill key information gaps and helps risk assessors identify the best levels of biological organization to conduct ERAs with differing goals.

Limitations of toxicity characterization in life cycle assessment: Can adverse outcome pathways provide a new foundation?

Life cycle assessment (LCA) has considerable merit for holistic evaluation of product planning, development, production, and disposal, with the inherent benefit of providing a forecast of potential health and environmental impacts. However, a technical review of current life cycle impact assessment (LCIA) methods revealed limitations within the biological effects assessment protocols, including: simplistic assessment approaches and models; an inability to integrate emerging types of toxicity data; a reliance on linear impact assessment models; a lack of methods to mitigate uncertainty; and no explicit consideration of effects in species of concern. The purpose of the current study is to demonstrate that a new concept in toxicological and regulatory assessment, the adverse outcome pathway (AOP), has many useful attributes of potential use to ameliorate many of these problems, to expand data utility and model robustness, and to enable more accurate and defensible biological effects assessments within LCIA. Background, context, and examples have been provided to demonstrate these potential benefits. We additionally propose that these benefits can be most effectively realized through development of quantitative AOPs (qAOPs) crafted to meet the needs of the LCIA framework. As a means to stimulate qAOP research and development in support of LCIA, we propose 3 conceptual classes of qAOP, each with unique inherent attributes for supporting LCIA: 1) mechanistic, including computational toxicology models; 2) probabilistic, including Bayesian networks and supervised machine learning models; and 3) weight of evidence, including models built using decision-analytic methods. Overall, we have highlighted a number of potential applications of qAOPs that can refine and add value to LCIA. As the AOP concept and support framework matures, we see the potential for qAOPs to serve a foundational role for next-generation effects characterization within LCIA. Integr Environ Assess Manag 2016;12:580-590. Published 2015. This article is a US Government work and is in the public domain in the USA.

Next-generation ecological risk assessment: Predicting risk from molecular initiation to ecosystem service delivery

Ecological risk assessment is the process of evaluating how likely it is that the environment may be impacted as the result of exposure to one or more chemicals and/or other stressors. It is not playing as large a role in environmental management decisions as it should be. A core challenge is that risk assessments often do not relate directly or transparently to protection goals. There have been exciting developments in in vitro testing and high-throughput systems that measure responses to chemicals at molecular and biochemical levels of organization, but the linkage between such responses and impacts of regulatory significance - whole organisms, populations, communities, and ecosystems - are not easily predictable. This article describes some recent developments that are directed at bridging this gap and providing more predictive models that can make robust links between what we typically measure in risk assessments and what we aim to protect.

A Computational Model of the Rainbow Trout Hypothalamus-Pituitary-Ovary-Liver Axis

Reproduction in fishes and other vertebrates represents the timely coordination of many endocrine factors that culminate in the production of mature, viable gametes. In recent years there has been rapid growth in understanding fish reproductive biology, which has been motivated in part by recognition of the potential effects that climate change, habitat destruction and contaminant exposure can have on natural and cultured fish populations. New approaches to understanding the impacts of these stressors are being developed that require a systems biology approach with more biologically accurate and detailed mathematical models. We have developed a multi-scale mathematical model of the female rainbow trout hypothalamus-pituitary-ovary-liver axis to use as a tool to help understand the functioning of the system and for extrapolation of laboratory findings of stressor impacts on specific components of the axis. The model describes the essential endocrine components of the female rainbow trout reproductive axis. The model also describes the stage specific growth of maturing oocytes within the ovary and permits the presence of sub-populations of oocytes at different stages of development. Model formulation and parametrization was largely based on previously published in vivo and in vitro data in rainbow trout and new data on the synthesis of gonadotropins in the pituitary. Model predictions were validated against several previously published data sets for annual changes in gonadotropins and estradiol in rainbow trout. Estimates of select model parameters can be obtained from in vitro assays using either quantitative (direct estimation of rate constants) or qualitative (relative change from control values) approaches. This is an important aspect of mathematical models as in vitro, cell-based assays are expected to provide the bulk of experimental data for future risk assessments and will require quantitative physiological models to extrapolate across biological scales.

Reproduction in fishes and other vertebrates represents the timely coordination of many endocrine factors that culminate in the production of mature, viable gametes. Improving the ability to estimate reproductive performance in fish is important, due to the growth of the aquaculture industry and the need to maintain adequate broodstock and concerns over the effects of anthropogenic stressors on feral fish populations. We present here a quantitative, mathematical model of the female rainbow trout reproductive cycle. We show how the model is able to accurately describe experimentally measured data associated with pituitary, ovarian and liver reproductive performance. We also use the model to describe similar data sets collected in rainbow trout by other researchers. An important value of quantitative biological models is the ability to simulate various physiological conditions, real or hypothetical. We demonstrate this by predicting the effects of exposure to an endocrine disruptor on oocyte growth. The need to limit cost and animal usage will encourage future experimental studies to use in vitro methods. The model presented here can assist with the extrapolation of in vitro effects to the whole fish.

Internal exposure dynamics drive the Adverse Outcome Pathways of synthetic glucocorticoids in fish

The Adverse Outcome Pathway (AOP) framework represents a valuable conceptual tool to systematically integrate existing toxicological knowledge from a mechanistic perspective to facilitate predictions of chemical-induced effects across species. However, its application for decision-making requires the transition from qualitative to quantitative AOP (qAOP). Here we used a fish model and the synthetic glucocorticoid beclomethasone dipropionate (BDP) to investigate the role of chemical-specific properties, pharmacokinetics, and internal exposure dynamics in the development of qAOPs. We generated a qAOP network based on drug plasma concentrations and focused on immunodepression, skin androgenisation, disruption of gluconeogenesis and reproductive performance. We showed that internal exposure dynamics and chemical-specific properties influence the development of qAOPs and their predictive power. Comparing the effects of two different glucocorticoids, we highlight how relatively similar in vitro hazard-based indicators can lead to different in vivo risk. This discrepancy can be predicted by their different uptake potential, pharmacokinetic (PK) and pharmacodynamic (PD) profiles. We recommend that the development phase of qAOPs should include the application of species-specific uptake and physiologically-based PK/PD models. This integration will significantly enhance the predictive power, enabling a more accurate assessment of the risk and the reliable transferability of qAOPs across chemicals.

Reproductive effects in fathead minnows (Pimphales promelas) following a 21 d exposure to 17α-ethinylestradiol

17α-ethinylestradiol (EE2) is a synthetic estrogen that is an active ingredient in oral contraception and hormone replacement therapy. Surveys of wastewater treatment plant effluents and surface waters throughout the world have reported EE2 concentrations in the ng/L range, and these low levels can cause significant reproductive effects in fish. This study tested the effects of three environmentally relevant EE2 concentrations: 0.47, 1.54 and 3.92 ng/L using a 21 d short-term reproductive assay to investigate the effects of EE2 on fathead minnow (Pimephales promelas) reproduction. The two highest EE2 concentrations tested in this study caused significant liver gene expression and induction of vitellogenin plasma protein in male fathead minnows. Exposure to 3.92 ng EE2/L increased the production of plasma vitellogenin in the females. Plasma estradiol concentrations were significantly reduced in females exposed to 1.54 and 3.92 ng EE2/L. All three tested concentrations significantly reduced fathead minnow egg production after a 21 d exposure to EE2. The results of this study indicate that the previously reported no observed adverse effect concentration (NOAEC) for EE2 on fathead minnow egg production (1.0 ng/L) may be too high. Because all three treatments resulted in significantly reduced egg production, the lowest observed adverse effect concentration (LOAEC) for EE2 on fathead minnow egg production is 0.47 ng EE2/L. This research estimates a NOAEC for fathead minnow reproduction at 0.24 ng EE2/L following a 21 d exposure. Additionally, induction of vitellogenin is a sensitive indicator of estrogen exposure but does not appear to be predictive of fathead minnow egg production.

Predicting Fecundity of Fathead Minnows (Pimephales promelas) Exposed to Endocrine-Disrupting Chemicals Using a MATLAB®-Based Model of Oocyte Growth Dynamics

Fish spawning is often used as an integrated measure of reproductive toxicity, and an indicator of aquatic ecosystem health in the context of forecasting potential population-level effects considered important for ecological risk assessment. Consequently, there is a need for flexible, widely-applicable, biologically-based models that can predict changes in fecundity in response to chemical exposures, based on readily measured biochemical endpoints, such as plasma vitellogenin (VTG) concentrations, as input parameters. Herein we describe a MATLAB® version of an oocyte growth dynamics model for fathead minnows (Pimephales promelas) with a graphical user interface based upon a previously published model developed with MCSim software and evaluated with data from fathead minnows exposed to an androgenic chemical, 17β-trenbolone. We extended the evaluation of our new model to include six chemicals that inhibit enzymes involved in steroid biosynthesis: fadrozole, ketoconazole, propiconazole, prochloraz, fenarimol, and trilostane. In addition, for unexposed fathead minnows from group spawning design studies, and those exposed to the six chemicals, we evaluated whether the model is capable of predicting the average number of eggs per spawn and the average number of spawns per female, which was not evaluated previously. The new model is significantly improved in terms of ease of use, platform independence, and utility for providing output in a format that can be used as input into a population dynamics model. Model-predicted minimum and maximum cumulative fecundity over time encompassed the observed data for fadrozole and most propiconazole, prochloraz, fenarimol and trilostane treatments, but did not consistently replicate results from ketoconazole treatments. For average fecundity (eggs•female(-1)•day(-1)), eggs per spawn, and the number of spawns per female, the range of model-predicted values generally encompassed the experimentally observed values. Overall, we found that the model predicts reproduction metrics robustly and its predictions capture the variability in the experimentally observed data.

Development and application of the adverse outcome pathway framework for understanding and predicting chronic toxicity: I. Challenges and research needs in ecotoxicology

To elucidate the effects of chemicals on populations of different species in the environment, efficient testing and modeling approaches are needed that consider multiple stressors and allow reliable extrapolation of responses across species. An adverse outcome pathway (AOP) is a concept that provides a framework for organizing knowledge about the progression of toxicity events across scales of biological organization that lead to adverse outcomes relevant for risk assessment. In this paper, we focus on exploring how the AOP concept can be used to guide research aimed at improving both our understanding of chronic toxicity, including delayed toxicity as well as epigenetic and transgenerational effects of chemicals, and our ability to predict adverse outcomes. A better understanding of the influence of subtle toxicity on individual and population fitness would support a broader integration of sublethal endpoints into risk assessment frameworks. Detailed mechanistic knowledge would facilitate the development of alternative testing methods as well as help prioritize higher tier toxicity testing. We argue that targeted development of AOPs supports both of these aspects by promoting the elucidation of molecular mechanisms and their contribution to relevant toxicity outcomes across biological scales. We further discuss information requirements and challenges in application of AOPs for chemical- and site-specific risk assessment and for extrapolation across species. We provide recommendations for potential extension of the AOP framework to incorporate information on exposure, toxicokinetics and situation-specific ecological contexts, and discuss common interfaces that can be employed to couple AOPs with computational modeling approaches and with evolutionary life history theory. The extended AOP framework can serve as a venue for integration of knowledge derived from various sources, including empirical data as well as molecular, quantitative and evolutionary-based models describing species responses to toxicants. This will allow a more efficient application of AOP knowledge for quantitative chemical- and site-specific risk assessment as well as for extrapolation across species in the future.

Determining the effects of a mixture of an endocrine disrupting compound, 17α-ethinylestradiol, and ammonia on fathead minnow (Pimephales promelas) reproduction

Aquatic organisms are exposed to a multitude of contaminants and to fully understand the impact of multiple stressors on fish populations, we must first understand the mechanism of action for each toxicant and how the combined effects manifest at the level of the individual. 17α-ethinylestradiol (EE2) has been known to cause adverse reproductive effects including reduced fecundity and fertility, intersex and skewed sex ratios in fish by mimicking naturally produced estrogen at low concentrations. Ammonia can cause adverse reproductive and mortality effects in individual fish through effects or damage to the central nervous system. Both EE2 and ammonia are found in most municipal effluents in various concentrations. A flow-through diluter system was used to test the individual effects of these two contaminants at their respective no observable adverse effect concentration (NOAEC) as well as their combined effects on fathead minnow, (Pimephales promelas) reproduction in a mixture exposure. While neither contaminant nor their mixture altered reproduction in terms of fecundity, their mixture resulted in significant fathead minnow mortality during a 21 d exposure. This study demonstrated the need to consider mixture effects when assessing risk for toxicity testing with multiple stressors.

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.

A physiologically based toxicokinetic model for the zebrafish Danio rerio

Zebrafish (Danio rerio) is a widely used model for toxicological studies, in particular those related to investigations on endocrine disruption. The development and regulatory use of in vivo and in vitro tests based on this species can be enhanced by toxicokinetic modeling. For this reason, we propose a physiologically based toxicokinetic (PBTK) model for zebrafish describing the uptake and disposition of organic chemicals. The model is based on literature data on zebrafish, other cyprinidae and other fish families, new experimental physiological information (volumes, lipids and water contents) obtained from zebrafish, and chemical-specific parameters predicted by generic models. The relevance of available models predicting the latter parameters was evaluated with respect to gill uptake and partition coefficients in zebrafish. This evaluation benefited from the fact that the influence of confounding factors such as body weight and temperature on ventilation rate was included in our model. The predictions for six chemicals (65 data points) yielded by our PBTK model were compared to available toxicokinetics data for zebrafish and 88% of them were within a factor of 5 of the corresponding experimental values. Sensitivity analysis highlighted that the 1-octanol/water partition coefficient, the metabolism rate, and all the parameters that enable the prediction of assimilation efficiency and partitioning of chemicals need to be precisely determined in order to allow an effective toxicokinetic modeling.

High-throughput sequencing and pathway analysis reveal alteration of the pituitary transcriptome by 17α-ethynylestradiol (EE2) in female coho salmon, Oncorhynchus kisutch

Considerable research has been done on the effects of endocrine disrupting chemicals (EDCs) on reproduction and gene expression in the brain, liver and gonads of teleost fish, but information on impacts to the pituitary gland are still limited despite its central role in regulating reproduction. The aim of this study was to further our understanding of the potential effects of natural and synthetic estrogens on the brain-pituitary-gonad axis in fish by determining the effects of 17α-ethynylestradiol (EE2) on the pituitary transcriptome. We exposed sub-adult coho salmon (Oncorhynchus kisutch) to 0 or 12 ng EE2/L for up to 6 weeks and effects on the pituitary transcriptome of females were assessed using high-throughput Illumina(®) sequencing, RNA-Seq and pathway analysis. After 1 or 6 weeks, 218 and 670 contiguous sequences (contigs) respectively, were differentially expressed in pituitaries of EE2-exposed fish relative to control. Two of the most highly up- and down-regulated contigs were luteinizing hormone β subunit (241-fold and 395-fold at 1 and 6 weeks, respectively) and follicle-stimulating hormone β subunit (-3.4-fold at 6 weeks). Additional contigs related to gonadotropin synthesis and release were differentially expressed in EE2-exposed fish relative to controls. These included contigs involved in gonadotropin releasing hormone (GNRH) and transforming growth factor-β signaling. There was an over-representation of significantly affected contigs in 33 and 18 canonical pathways at 1 and 6 weeks, respectively, including circadian rhythm signaling, calcium signaling, peroxisome proliferator-activated receptor (PPAR) signaling, PPARα/retinoid x receptor α activation, and netrin signaling. Network analysis identified potential interactions between genes involved in circadian rhythm and GNRH signaling, suggesting possible effects of EE2 on timing of reproductive events.

Developing predictive approaches to characterize adaptive responses of the reproductive endocrine axis to aromatase inhibition: I. Data generation in a small fish model

Adaptive or compensatory responses to chemical exposure can significantly influence in vivo concentration-duration-response relationships. This study provided data to support development of a computational dynamic model of the hypothalamic-pituitary-gonadal axis of a model vertebrate and its response to aromatase inhibitors as a class of endocrine active chemicals. Fathead minnows (Pimephales promelas) were either exposed to the aromatase inhibitor fadrozole (0.5 or 30 μg/l) continuously for 1, 8, 12, 16, 20, 24, or 28 days or exposed for 8 days and then held in control water (no fadrozole) for an additional 4, 8, 12, 16, or 20 days. The time course of effects on ovarian steroid production, circulating 17β-estradiol (E2) and vitellogenin (VTG) concentrations, and expression of steroidogenesis-related genes in the ovary was measured. Exposure to 30 μg fadrozole/l significantly reduced plasma E2 and VTG concentrations after just 1 day and those effects persisted throughout 28 days of exposure. In contrast, ex vivo E2 production was similar to that of controls on day 8-28 of exposure, whereas transcripts coding for aromatase and follicle-stimulating hormone receptor were elevated, suggesting a compensatory response. Following cessation of fadrozole exposure, ex vivo E2 and plasma E2 concentrations exceeded and then recovered to control levels, but plasma VTG concentrations did not, even after 20 days of depuration. Collectively these data provide several new insights into the nature and time course of adaptive responses to an aromatase inhibitor that support development of a computational model (see companion article).

A computational model to predict rat ovarian steroid secretion from in vitro experiments with endocrine disruptors

A finely tuned balance between estrogens and androgens controls reproductive functions, and the last step of steroidogenesis plays a key role in maintaining that balance. Environmental toxicants are a serious health concern, and numerous studies have been devoted to studying the effects of endocrine disrupting chemicals (EDCs). The effects of EDCs on steroidogenic enzymes may influence steroid secretion and thus lead to reproductive toxicity. To predict hormonal balance disruption on the basis of data on aromatase activity and mRNA level modulation obtained in vitro on granulosa cells, we developed a mathematical model for the last gonadal steps of the sex steroid synthesis pathway. The model can simulate the ovarian synthesis and secretion of estrone, estradiol, androstenedione, and testosterone, and their response to endocrine disruption. The model is able to predict ovarian sex steroid concentrations under normal estrous cycle in female rat, and ovarian estradiol concentrations in adult female rats exposed to atrazine, bisphenol A, metabolites of methoxychlor or vinclozolin, and letrozole.

Constraints-based stoichiometric analysis of hypoxic stress on steroidogenesis in fathead minnows, Pimephales promelas

In this study, an in silico genome-scale metabolic model of steroidogenesis was used to investigate the effects of hypoxic stress on steroid hormone productions in fish. Adult female fathead minnows (Pimephales promelas) were exposed to hypoxia for 7 days with fish sub-sampled on days 1, 3 and 7 of exposure. At each time point, selected steroid enzyme gene expressions and steroid hormone productions were quantified in ovaries. Fold changes in steroid enzyme gene expressions were used to qualitatively scale transcript enzyme reaction constraints (akin to the range of an enzyme’s catalytic activity) in the in silico model. Subsequently, in silico predicted steroid hormone productions were qualitatively compared with experimental results. Key findings were as follows. (1) In silico gene deletion analysis identified highly conserved ‘essential’ genes required for steroid hormone productions. These agreed well (75%) with literature-published genes downregulated in vertebrates (fish and mammal) exposed to hypoxia. (2) Quantification of steroid hormones produced ex vivo from ovaries showed a significant reduction for 17β-estradiol and 17α,20β-dihydroxypregnenone production after 24 h (day 1) of exposure. This lowered 17β-estradiol production was concomitant with downregulation of cyp19a1a gene expression in ovaries. In silico predictions showed agreement with experimentation by predicting effects on estrogen (17β-estradiol and estrone) production. (3) Stochastic sampling of in silico reactions indicated that cholesterol uptake and catalysis to pregnenolone along with estrogen methyltransferase and glucuronidation reactions were also impacted by hypoxia. Taken together, this in silico analysis introduces a powerful model for pathway analysis that can lend insights on the effects of various stressor scenarios on metabolic functions.