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

Case Study in 21st-Century Ecotoxicology: Using In Vitro Aromatase Inhibition Data to Predict Reproductive Outcomes in Fish In Vivo

To reduce the use of intact animals for chemical safety testing, while ensuring protection of ecosystems and human health, there is a demand for new approach methodologies (NAMs) that provide relevant scientific information at a quality equivalent to or better than traditional approaches. The present case study examined whether bioactivity and associated potency measured in an in vitro screening assay for aromatase inhibition could be used together with an adverse outcome pathway (AOP) and mechanistically based computational models to predict previously uncharacterized in vivo effects. Model simulations were used to inform designs of 60-h and 10-21-day in vivo exposures of adult fathead minnows (Pimephales promelas) to three or four test concentrations of the in vitro aromatase inhibitor imazalil ranging from 0.12 to 260 µg/L water. Consistent with an AOP linking aromatase inhibition to reproductive impairment in fish, exposure to the fungicide resulted in significant reductions in ex vivo production of 17β-estradiol (E2) by ovary tissue (≥165 µg imazalil/L), plasma E2 concentrations (≥74 µg imazalil/L), vitellogenin (Vtg) messenger RNA expression (≥165 µg imazalil/L), Vtg plasma concentrations (≥74 µg imazalil/L), uptake of Vtg into oocytes (≥260 µg imazalil/L), and overall reproductive output in terms of cumulative fecundity, number of spawning events, and eggs per spawning event (≥24 µg imazalil/L). Despite many potential sources of uncertainty in potency and efficacy estimates based on model simulations, observed magnitudes of apical effects were quite consistent with model predictions, and in vivo potency was within an order of magnitude of that predicted based on in vitro relative potency. Overall, our study suggests that NAMs and AOP-based approaches can support meaningful reduction and refinement of animal testing. Environ Toxicol Chem 2023;42:100-116. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Effects of nanoplastic on toxicity of azole fungicides (ketoconazole and fluconazole) in zebrafish embryos

The ubiquity of nanoplastics (NPs) raises concerns about their interactions and combined toxicity with other common contaminants. Although azoles are present throughout the natural environment, their interactions with NP are not well known. We investigated the effects of polystyrene (PS) NP on the toxicity of ketoconazole (KCZ) and fluconazole (FCZ) in zebrafish embryos using the developmental toxicity, oxidative-stress-related biochemical parameters, and expression of genes related to neurotoxicity (ache), cardiotoxicity (gata4, bmp4), inflammation (il1b), oxidative stress (sod1, sod2, cyp1a), and apoptosis (bax, bcl2). Co-exposure to NP (1 mg/L) and KCZ/FCZ (1 mg/L) for 96 h reduced the hatching rate, survival rate, and heart rate and increased the malformation rate and catalase activity. The bax/bcl2 ratio, an apoptosis indicator, was higher after NP, KCZ, or FCZ treatment. However, the bax/bcl2 ratio after exposure to NP + KCZ or NP + FCZ was much higher than that after single exposure. Overall, the results indicated that NP aggravated the toxicity of azole by significantly increasing the reactive oxygen species, lipid peroxidation and altering the expression of oxidative-stress- and apoptosis-related genes. The interactive toxicity of PS NP with KCZ/FCZ reported in this study emphasises the need for caution in the release of azole fungicides in the environment.

A comprehensive review on environmental toxicity of azole compounds to fish

BACKGROUND

Azoles are considered as one of the most efficient fungicides for the treatment of humans, animals, and plant fungal pathogens. They are of significant clinical importance as antifungal drugs and are widely used in personal care products, ultraviolet stabilizers, and in aircraft for its anti-corrosive properties. The prevalence of azole compounds in the natural environment and its accumulation in fish raises questions about its impact on aquatic organisms.

OBJECTIVES

The objective of this paper is to review the scientific studies on the effects of azole compounds in fish and to discuss future opportunities for the risk evaluation.

METHODS

A systematic literature search was conducted on Web of Science, PubMed, and ScienceDirect to locate peer-reviewed scientific articles on occurrence, environmental fate, and toxicological impact of azole fungicides on fish.

RESULTS

Studies included in this review provide ample evidence that azole compounds are not only commonly detected in the natural environment but also cause several detrimental effects on fish. Future studies with environmentally relevant concentrations of azole alone or in combination with other commonly occurring contaminants in a multigenerational study could provide a better understanding.

CONCLUSION

Based on current knowledge and studies reporting adverse biological effects of azole on fish, considerable attention is required for better management and effective ecological risk assessment of these emerging contaminants.

Building and Applying Quantitative Adverse Outcome Pathway Models for Chemical Hazard and Risk Assessment

An important goal in toxicology is the development of new ways to increase the speed, accuracy, and applicability of chemical hazard and risk assessment approaches. A promising route is the integration of in vitro assays with biological pathway information. We examined how the adverse outcome pathway (AOP) framework can be used to develop pathway-based quantitative models useful for regulatory chemical safety assessment. By using AOPs as initial conceptual models and the AOP knowledge base as a source of data on key event relationships, different methods can be applied to develop computational quantitative AOP models (qAOPs) relevant for decision making. A qAOP model may not necessarily have the same structure as the AOP it is based on. Useful AOP modeling methods range from statistical, Bayesian networks, regression, and ordinary differential equations to individual-based models and should be chosen according to the questions being asked and the data available. We discuss the need for toxicokinetic models to provide linkages between exposure and qAOPs, to extrapolate from in vitro to in vivo, and to extrapolate across species. Finally, we identify best practices for modeling and model building and the necessity for transparent and comprehensive documentation to gain confidence in the use of qAOP models and ultimately their use in regulatory applications. Environ Toxicol Chem 2019;38:1850-1865. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

In-silico approach for drug induced liver injury prediction: Recent advances

Drug induced liver injury (DILI) is the prime cause of liver disfunction which may lead to mild non-specific symptoms to more severe signs like hepatitis, cholestasis, cirrhosis and jaundice. Not only the prescription medications, but the consumption of herbs and health supplements have also been reported to cause these adverse reactions resulting into high mortality rates and post marketing withdrawal of drugs. Due to the continuously increasing DILI incidences in recent years, robust prediction methods with high accuracy, specificity and sensitivity are of priority. Bioinformatics is the emerging field of science that has been used in the past few years to explore the mechanisms of DILI. The major emphasis of this review is the recent advances of in silico tools for the diagnostic and therapeutic interventions of DILI. These tools have been developed and widely used in the past few years for the prediction of pathways induced from both hepatotoxic as well as hepatoprotective Chinese drugs and for the identification of DILI specific biomarkers for prognostic purpose. In addition to this, advanced machine learning models have been developed for the classification of drugs into DILI causing and non-DILI causing. Moreover, development of 3 class models over 2 class offers better understanding of multi-class DILI risks and at the same time providing authentic prediction of toxicity during drug designing before clinical trials.

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.

Population-relevant endpoints in the evaluation of endocrine-active substances (EAS) for ecotoxicological hazard and risk assessment

For ecotoxicological risk assessment, endocrine disruptors require the establishment of an endocrine mode of action (MoA) with a plausible link to a population-relevant adverse effect. Current ecotoxicity test methods incorporate mostly apical endpoints although some also include mechanistic endpoints, subcellular-through-organ level, which can help establish an endocrine MoA. However, the link between these endpoints and adverse population-level effects is often unclear. The case studies of endocrine-active substances (EAS) (tributyltin, ethinylestradiol, perchlorate, trenbolone, propiconazole, and vinclozolin) evaluated from the Society of Environmental Toxicology and Chemistry (SETAC) Pellston Workshop^® "Ecotoxicological Hazard and Risk Assessment Approaches for Endocrine-Active Substances (EHRA)" were used to evaluate the population relevance of toxicity endpoints in various taxa according to regulatory endocrine-disruptor frameworks such as the Organisation for Economic Co-operation and Development (OECD) Conceptual Framework for Testing and Assessment of Endocrine Disruptors. A wide variety of potentially endocrine-relevant endpoints were identified for mollusks, fish, amphibians, birds, and mammals, although the strength of the relationship between test endpoints and population-level effects was often uncertain. Furthermore, testing alone is insufficient for assessing potential adaptation and recovery processes in exposed populations. For this purpose, models that link effects observed in laboratory tests to the dynamics of wildlife populations appear to be necessary, and their development requires reliable and robust data. As our understanding of endocrine perturbations and key event relationships improves, adverse population-level effects will be more easily and accurately predicted. Integr Environ Assess Manag 2017;13:317-330. © 2017 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals, Inc. on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Recommended approaches to the scientific evaluation of ecotoxicological hazards and risks of endocrine-active substances

A SETAC Pellston Workshop® "Environmental Hazard and Risk Assessment Approaches for Endocrine-Active Substances (EHRA)" was held in February 2016 in Pensacola, Florida, USA. The primary objective of the workshop was to provide advice, based on current scientific understanding, to regulators and policy makers; the aim being to make considered, informed decisions on whether to select an ecotoxicological hazard- or a risk-based approach for regulating a given endocrine-disrupting substance (EDS) under review. The workshop additionally considered recent developments in the identification of EDS. Case studies were undertaken on 6 endocrine-active substances (EAS-not necessarily proven EDS, but substances known to interact directly with the endocrine system) that are representative of a range of perturbations of the endocrine system and considered to be data rich in relevant information at multiple biological levels of organization for 1 or more ecologically relevant taxa. The substances selected were 17α-ethinylestradiol, perchlorate, propiconazole, 17β-trenbolone, tributyltin, and vinclozolin. The 6 case studies were not comprehensive safety evaluations but provided foundations for clarifying key issues and procedures that should be considered when assessing the ecotoxicological hazards and risks of EAS and EDS. The workshop also highlighted areas of scientific uncertainty, and made specific recommendations for research and methods-development to resolve some of the identified issues. The present paper provides broad guidance for scientists in regulatory authorities, industry, and academia on issues likely to arise during the ecotoxicological hazard and risk assessment of EAS and EDS. The primary conclusion of this paper, and of the SETAC Pellston Workshop on which it is based, is that if data on environmental exposure, effects on sensitive species and life-stages, delayed effects, and effects at low concentrations are robust, initiating environmental risk assessment of EDS is scientifically sound and sufficiently reliable and protective of the environment. In the absence of such data, assessment on the basis of hazard is scientifically justified until such time as relevant new information is available. Integr Environ Assess Manag 2017;13:267-279. © 2017 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals, Inc. on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

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.

Editor's Highlight: Computational Modeling of Plasma Vitellogenin Alterations in Response to Aromatase Inhibition in Fathead Minnows

In vertebrates, conversion of testosterone into 17β-estradiol (E2) is catalyzed by cytochrome P450 (CYP) 19A aromatase. An important role of E2 in oviparous vertebrates such as fish is stimulation of hepatic synthesis of the glycolipoprotein vitellogenin (VTG), an egg yolk precursor essential to oocyte development and larval survival. In fathead minnows (FHMs) (Pimephales promelas) exposed to the aromatase inhibitor fadrozole, plasma VTG levels do not change in concert with plasma E2 levels. Specifically, while plasma VTG and E2 levels both drop quickly when aromatase is first inhibited, the recovery of plasma VTG upon cessation of aromatase inhibition is substantially delayed relative to the recovery of plasma E2. We modified an existing computational model of the FHM hypothalamic-pituitary-gonadal axis to evaluate alternative hypotheses that might explain this delay. In the first hypothesis, a feedback loop involving active transport of VTG from the blood into the ovary is used. The activity of the transporter is negatively regulated by ovarian VTG. In the second hypothesis, a type 1 coherent feed-forward loop is implemented in the liver. This loop has 2 arms, both requiring E2 activation. The first arm describes direct, canonical E2-driven transcriptional induction of VTG, and the second describes an E2-driven intermediate transcriptional regulator that is also required for VTG synthesis. Both hypotheses accurately described the observed VTG dynamics. This result could be used to guide design of laboratory experiments intended to determine if either of the motifs, or perhaps even both of them, actually do control VTG dynamics in FHMs exposed to aromatase inhibitors.