FutureTox IV Workshop Summary: Predictive Toxicology for Healthy Children.

FutureTox IV, a Society of Toxicology Contemporary Concepts in Toxicology workshop, was held in November 2018. Building upon FutureTox I, II, and III, this conference focused on the latest science and technology for in vitro profiling and in silico modeling as it relates to predictive developmental and reproductive toxicity (DART). Publicly available high-throughput screening data sets are now available for broad in vitro profiling of bioactivities across large inventories of chemicals. Coupling this vast amount of mechanistic data with a deeper understanding of molecular embryology and post-natal development lays the groundwork for using new approach methodologies (NAMs) to evaluate chemical toxicity, drug efficacy, and safety assessment for embryo-fetal development. NAM is a term recently adopted in reference to any technology, methodology, approach, or combination thereof that can be used to provide information on chemical hazard and risk assessment to avoid the use of intact animals (U.S. Environmental Protection Agency [EPA], Strategic plan to promote the development and implementation of alternative test methods within the tsca program, 2018, https://www.epa.gov/sites/production/files/2018-06/documents/epa_alt_strat_plan_6-20-18_clean_final.pdf). There are challenges to implementing NAMs to evaluate chemicals for developmental toxicity compared with adult toxicity. This forum article reviews the 2018 workshop activities, highlighting challenges and opportunities for applying NAMs for adverse pregnancy outcomes (eg, preterm labor, malformations, low birth weight) as well as disorders manifesting postnatally (eg, neurodevelopmental impairment, breast cancer, cardiovascular disease, fertility). DART is an important concern for different regulatory statutes and test guidelines. Leveraging advancements in such approaches and the accompanying efficiencies to detecting potential hazards to human development are the unifying concepts toward implementing NAMs in DART testing. Although use of NAMs for higher level regulatory decision making is still on the horizon, the conference highlighted novel testing platforms and computational models that cover multiple levels of biological organization, with the unique temporal dynamics of embryonic development, and novel approaches for estimating toxicokinetic parameters essential in supporting in vitro to in vivo extrapolation.

Chronic Toxicity of Cobalt to Marine Organisms: Application of a Species Sensitivity Distribution Approach to Develop International Water Quality Standards.

Water quality standards for cobalt (Co) have yet to be developed for the European Union or the United States. The primary objective of the present study was to produce a data set comprising marine Co toxicity data that could be used by both the European Union and the United States to determine a predicted-no-effect concentration and ambient water quality criteria, respectively. Ten marine species, ranging from algae to fish, were subjected to chronic Co toxicity tests that were designed to meet international water quality testing standards. Chronic 10% effect concentration values ranged from a low of 1.23 µg dissolved Co/L for red algae (Champia parvula) to a high of 31 800 µg dissolved Co/L for sheepshead minnow (Cyprinodon variegatus). The species sensitivity ranking for chronic marine Co exposure was as follows (from most to least sensitive): C. parvula > Neanthes arenaceodentata (polychaete) > Americamysis bahia (mysid shrimp) > Skeletonema costatum (marine diatom) > Dendraster excentricus (sand dollar) > Mytilus galloprovincialis (mussel) > Strongylocentrotus purpuratus (purple sea urchin) > Crassostrea gigas (oyster) > Dunaliella tertiolecta (marine flagellate) > C. variegatus. Chronic test results indicated that invertebrate and plant species were substantially more sensitive to Co exposure than fish. The chronic toxicity data were used to calculate a species sensitivity distribution, from which a hazardous concentration 5th percentile of 7.09 µg dissolved Co/L (95% CI 0.025-47.3 µg Co/L) was derived. Environ Toxicol Chem 2021;40:1405-1418. © 2021 SETAC.

A cross-platform approach to characterize and screen potential neurovascular unit toxicants.

Development of the neurovascular unit (NVU) is a complex, multistage process that requires orchestrated cell signaling mechanisms across several cell types and ultimately results in formation of the blood-brain barrier. Typical high-throughput screening (HTS) assays investigate single biochemical or single cell responses following chemical insult. As the NVU comprises multiple cell types interacting at various stages of development, a methodology combining high-throughput results across pertinent cell-based assays is needed to investigate potential chemical-induced disruption to the development of this complex cell system. To this end, we implemented a novel method for screening putative NVU disruptors across diverse assay platforms to predict chemical perturbation of the developing NVU. HTS assay results measuring chemical-induced perturbations to cellular key events across angiogenic and neurogenic outcomes in vitro were combined to create a cell-based prioritization of NVU hazard. Chemicals were grouped according to similar modes of action to train a logistic regression literature model on a training set of 38 chemicals. This model utilizes the chemical-specific pairwise mutual information score for PubMed MeSH annotations to represent a quantitative measure of previously published results. Taken together, this study presents a methodology to investigate NVU developmental hazard using cell-based HTS assays and literature evidence to prioritize screening of putative NVU disruptors towards a knowledge-driven characterization of neurovascular developmental toxicity. The results from these screening efforts demonstrate that chemicals representing a range of putative vascular disrupting compound (pVDC) scores can also produce effects on neurogenic outcomes and characterizes possible modes of action for disrupting the developing NVU.

Profiling the ToxCast Library With a Pluripotent Human (H9) Stem Cell Line-Based Biomarker Assay for Developmental Toxicity.

The Stemina devTOX quickPredict platform is a human pluripotent stem cell-based assay that predicts the developmental toxicity potential based on changes in cellular metabolism following chemical exposure [Palmer, J. A., Smith, A. M., Egnash, L. A., Conard, K. R., West, P. R., Burrier, R. E., Donley, E. L. R., and Kirchner, F. R. (2013). Establishment and assessment of a new human embryonic stem cell-based biomarker assay for developmental toxicity screening. Birth Defects Res. B Dev. Reprod. Toxicol. 98, 343-363]. Using this assay, we screened 1065 ToxCast phase I and II chemicals in single-concentration or concentration-response for the targeted biomarker (ratio of ornithine to cystine secreted or consumed from the media). The dataset from the Stemina (STM) assay is annotated in the ToxCast portfolio as STM. Major findings from the analysis of ToxCast_STM dataset include (1) 19% of 1065 chemicals yielded a prediction of developmental toxicity, (2) assay performance reached 79%-82% accuracy with high specificity (> 84%) but modest sensitivity (< 67%) when compared with in vivo animal models of human prenatal developmental toxicity, (3) sensitivity improved as more stringent weights of evidence requirements were applied to the animal studies, and (4) statistical analysis of the most potent chemical hits on specific biochemical targets in ToxCast revealed positive and negative associations with the STM response, providing insights into the mechanistic underpinnings of the targeted endpoint and its biological domain. The results of this study will be useful to improving our ability to predict in vivo developmental toxicants based on in vitro data and in silico models.

Molecular characterization of a toxicological tipping point during human stem cell differentiation.

Endoderm gives rise to the gut tube in the early embryo. We differentiated human induced pluripotent stem cells (hiPSCs) to embryonic endoderm to identify a "tipping point" at which the developing system did not recover from perturbations caused by exposure to all-trans retinoic acid (ATRA). Differentiating hiPSC-derived endoderm exposed to five concentrations of ATRA between 0.001 and 10 µM at 6 h, 96 h, or 192 h was assessed for forkhead box A2 (FOXA2) protein expression and global gene transcript expression. A tipping point of 17 ±â€¯11 nM was identified where patterns of differentially expressed genes supported a developmental trajectory shift indicating a potential teratogenic outcome. This concentration is between women's endogenous ATRA blood plasma levels and teratogenic levels of circulating isotretinoin, an ATRA isomer used to treat acne. Taken together, these data suggest that this approach is a sensitive method to identify a point of departure for chemicals that impact early embryonic development.

Systems Modeling of Developmental Vascular Toxicity.

The more than 80,000 chemicals in commerce present a challenge for hazard assessments that toxicity testing in the 21st century strives to address through high-throughput screening (HTS) assays. Assessing chemical effects on human development adds an additional layer of complexity to the screening, with a need to capture complex and dynamic events essential for proper embryo-fetal development. HTS data from ToxCast/Tox21 informs systems toxicology models, which incorporate molecular targets and biological pathways into mechanistic models describing the effects of chemicals on human cells, 3D organotypic culture models, and small model organisms. Adverse Outcome Pathways (AOPs) provide a useful framework for integrating the evidence derived from these in silico and in vitro systems to inform chemical hazard characterization. To illustrate this formulation, we have built an AOP for developmental toxicity through a mode of action linked to embryonic vascular disruption (Aop43). Here, we review the model for quantitative prediction of developmental vascular toxicity from ToxCast HTS data and compare the HTS results to functional vascular development assays in complex cell systems, virtual tissues, and small model organisms. ToxCast HTS predictions from several published and unpublished assays covering different aspects of the angiogenic cycle were generated for a test set of 38 chemicals representing a range of putative vascular disrupting compounds (pVDCs). Results boost confidence in the capacity to predict adverse developmental outcomes from HTS in vitro data and model computational dynamics for in silico reconstruction of developmental systems biology. Finally, we demonstrate the integration of the AOP and developmental systems toxicology to investigate the unique modes of action of two angiogenesis inhibitors.

Blood-brain barrier development: Systems modeling and predictive toxicology.

The blood-brain barrier (BBB) serves as a gateway for passage of drugs, chemicals, nutrients, metabolites, and hormones between vascular and neural compartments in the brain. Here, we review BBB development with regard to the microphysiology of the neurovascular unit (NVU) and the impact of BBB disruption on brain development. Our focus is on modeling these complex systems. Extant in silico models are available as tools to predict the probability of drug/chemical passage across the BBB; in vitro platforms for high-throughput screening and high-content imaging provide novel data streams for profiling chemical-biological interactions; and engineered human cell-based microphysiological systems provide empirical models with which to investigate the dynamics of NVU function. Computational models are needed that bring together kinetic and dynamic aspects of NVU function across gestation and under various physiological and toxicological scenarios. This integration will inform adverse outcome pathways to reduce uncertainty in translating in vitro data and in silico models for use in risk assessments that aim to protect neurodevelopmental health.

Developmental bisphenol A exposure impairs sperm function and reproduction in zebrafish.

The developmental and reproductive toxicity of bisphenol A (BPA) has been demonstrated in a variety of model systems. Zebrafish (Danio rerio) were waterborne-exposed to BPA during three different developmental stages: embryonic period:6 h post fertilization (hpf) to 5 months post fertilization (mpf); larval period: 6 days post fertilization (dpf) to 5 mpf; and sexually mature period: 3 mpf to 5 mpf. Evaluations included F0 adult growth, reproduction parameters, and F1 offspring development. BPA exposure did not affect zebrafish growth in any of exposure groups. Testis weight was decreased only following the 6 hpf to 5 mpf 0.001 µM BPA exposure. The lowest effect level indicated by a reduction in sperm volume, density, motility, and velocity across a range of exposure durations was 0.001 µM, with all but sperm density significant for the longest exposure duration, which was also the only significant endpoint for the lowest exposure concentration in the 3-5 mpf exposure group. Nonmonotonic concentration-response curves were noted for all F0 reproductive endpoints for at least one of the two longest exposure durations. For the F1 offspring of fish exposed from 6 hpf to 5 mpf, malformations and mortality were increased following 0.001 µM BPA exposure, while egg production and fertilization were reduced in higher concentration treatment groups. Overall, BPA exposure during three different developmental periods impaired zebrafish reproductive development, with most significance changes found in the lowest concentration treatment groups. Genetic impacts on gamete development may underlie the secondary effects of reduced fertilization rate, embryonic mortality, and malformations.

Global gene expression analysis reveals pathway differences between teratogenic and non-teratogenic exposure concentrations of bisphenol A and 17ß-estradiol in embryonic zebrafish.

Transient developmental exposure to 0.1µM bisphenol A (BPA) results in larval zebrafish hyperactivity and learning impairments in the adult, while exposure to 80µM BPA results in teratogenic responses, including craniofacial abnormalities and edema. The mode of action underlying these effects is unclear. We used global gene expression analysis to identify candidate genes and signaling pathways that mediate BPA's developmental toxicity in zebrafish. Exposure concentrations were selected and anchored to the positive control, 17ß-estradiol (E2), based on previously determined behavioral or teratogenic phenotypes. Functional analysis of differentially expressed genes revealed distinct expression profiles at 24h post fertilization for 0.1µM versus 80µM BPA and 0.1µM versus 15µM E2 exposure, identification of prothrombin activation as a top canonical pathway impacted by both 0.1µM BPA and 0.1µM E2 exposure, and suppressed expression of several genes involved in nervous system development and function following 0.1µM BPA exposure.

Experimental evidence for evolved tolerance to avian malaria in a wild population of low elevation Hawai'i 'Amakihi (Hemignathus virens).

Introduced vector-borne diseases, particularly avian malaria (Plasmodium relictum) and avian pox virus (Avipoxvirus spp.), continue to play significant roles in the decline and extinction of native forest birds in the Hawaiian Islands. Hawaiian honeycreepers are particularly susceptible to avian malaria and have survived into this century largely because of persistence of high elevation refugia on Kaua'i, Maui, and Hawai'i Islands, where transmission is limited by cool temperatures. The long term stability of these refugia is increasingly threatened by warming trends associated with global climate change. Since cost effective and practical methods of vector control in many of these remote, rugged areas are lacking, adaptation through processes of natural selection may be the best long-term hope for recovery of many of these species. We document emergence of tolerance rather than resistance to avian malaria in a recent, rapidly expanding low elevation population of Hawai'i 'Amakihi (Hemignathus virens) on the island of Hawai'i. Experimentally infected low elevation birds had lower mortality, lower reticulocyte counts during recovery from acute infection, lower weight loss, and no declines in food consumption relative to experimentally infected high elevation Hawai'i 'Amakihi in spite of similar intensities of infection. Emergence of this population provides an exceptional opportunity for determining physiological mechanisms and genetic markers associated with malaria tolerance that can be used to evaluate whether other, more threatened species have the capacity to adapt to this disease.

Persistent adult zebrafish behavioral deficits results from acute embryonic exposure to gold nanoparticles.

As the number of products containing nanomaterials increase, human exposure to nanoparticles (NPs) is unavoidable. Presently, few studies focus on the potential long-term consequences of developmental NP exposure. In this study, zebrafish embryos were acutely exposed to three gold NPs that possess functional groups with differing surface charge. Embryos were exposed to 50 µg/mL of 1.5 nm gold nanoparticles (AuNPs) possessing negatively charged 2-mercaptoethanesulfonic acid (MES) or neutral 2-(2-(2-mercaptoethoxy)ethoxy)ethanol (MEEE) ligands or 10 µg/mL of the AuNPs possessing positively charged trimethylammoniumethanethiol (TMAT). Both MES- and TMAT-AuNP exposed embryos exhibited hypo-locomotor activity, while those exposed to MEEE-AuNPs did not. A subset of embryos that were exposed to 1.5 nm MES- and TMAT-AuNPs during development from 6 to 120 h post fertilization was raised to adulthood. Behavioral abnormalities and the number of survivors into adulthood were evaluated at 122 days post fertilization. We found that both treatments induced abnormal startle behavior following a tap stimulus. However, the MES-AuNPs exposed group also exhibited abnormal adult behavior in the light and had a lower survivorship into adulthood. This study demonstrates that acute, developmental exposure to 1.5 nm MES- and TMAT-AuNPs, two NPs differing only in the functional group, affects larval behavior, with behavioral effects persisting into adulthood.

Neurodevelopmental low-dose bisphenol A exposure leads to early life-stage hyperactivity and learning deficits in adult zebrafish.

Developmental bisphenol A (BPA) exposure has been implicated in adverse behavior and learning deficits. The mode of action underlying these effects is unclear. The objectives of this study were to identify whether low-dose, developmental BPA exposure affects larval zebrafish locomotor behavior and whether learning deficits occur in adults exposed during development. Two control compounds, 17ß-estradiol (an estrogen receptor ligand) and GSK4716 (a synthetic estrogen-related receptor gamma ligand), were included. Larval toxicity assays were used to determine appropriate BPA, 17ß-estradiol, and GSK4716 concentrations for behavior testing. BPA tissue uptake was analyzed using HPLC and lower doses were extrapolated using a linear regression analysis. Larval behavior tests were conducted using a ViewPoint Zebrabox. Adult learning tests were conducted using a custom-built T-maze. BPA exposure to <30µM was non-teratogenic. Neurodevelopmental BPA exposure to 0.01, 0.1, or 1µM led to larval hyperactivity or learning deficits in adult zebrafish. Exposure to 0.1µM 17ß-estradiol or GSK4716 also led to larval hyperactivity. This study demonstrates the efficacy of using the zebrafish model for studying the neurobehavioral effects of low-dose developmental BPA exposure.

The anti-inflammatory drug leflunomide is an agonist of the aryl hydrocarbon receptor.

BACKGROUND: The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that mediates the toxicity and biological activity of dioxins and related chemicals. The AhR influences a variety of processes involved in cellular growth and differentiation, and recent studies have suggested that the AhR is a potential target for immune-mediated diseases. METHODOLOGY/PRINCIPAL FINDINGS: During a screen for molecules that activate the AhR, leflunomide, an immunomodulatory drug presently used in the clinic for the treatment of rheumatoid arthritis, was identified as an AhR agonist. We aimed to determine whether any biological activity of leflunomide could be attributed to a previously unappreciated interaction with the AhR. The currently established mechanism of action of leflunomide involves its metabolism to A771726, possibly by cytochrome P450 enzymes, followed by inhibition of de novo pyrimidine biosynthesis by A771726. Our results demonstrate that leflunomide, but not its metabolite A771726, caused nuclear translocation of AhR into the nucleus and increased expression of AhR-responsive reporter genes and endogenous AhR target genes in an AhR-dependent manner. In silico Molecular Docking studies employing AhR ligand binding domain revealed favorable binding energy for leflunomide, but not for A771726. Further, leflunomide, but not A771726, inhibited in vivo epimorphic regeneration in a zebrafish model of tissue regeneration in an AhR-dependent manner. However, suppression of lymphocyte proliferation by leflunomide or A771726 was not dependent on AhR. CONCLUSIONS: These data reveal that leflunomide, an anti-inflammatory drug, is an agonist of the AhR. Our findings link AhR activation by leflunomide to inhibition of fin regeneration in zebrafish. Identification of alternative AhR agonists is a critical step in evaluating the AhR as a therapeutic target for the treatment of immune disorders.