Molecular aspects of azoles-induced teratogenesis

Identifying chemicals based on receptor binding/bioactivation/mechanistic explanation associated with potential to elicit hepatotoxicity and to support structure activity relationship-based read-across

The liver is the most common target organ in toxicology studies. The development of chemical structural alerts for identifying hepatotoxicity will play an important role in in silico model prediction and help strengthen the identification of analogs used in structure activity relationship (SAR)- based read-across. The aim of the current study is development of an SAR-based expert-system decision tree for screening of hepatotoxicants across a wide range of chemistry space and proposed modes of action for clustering of chemicals using defined core chemical categories based on receptor-binding or bioactivation. The decision tree is based on ∼ 1180 different chemicals that were reviewed for hepatotoxicity information. Knowledge of chemical receptor binding, metabolism and mechanistic information were used to group these chemicals into 16 different categories and 102 subcategories: four categories describe binders to 9 different receptors, 11 categories are associated with possible reactive metabolites (RMs) and there is one miscellaneous category. Each chemical subcategory has been associated with possible modes of action (MOAs) or similar key structural features. This decision tree can help to screen potential liver toxicants associated with core structural alerts of receptor binding and/or RMs and be used as a component of weight of evidence decisions based on SAR read-across, and to fill data gaps.

Triazole fungicides induce adipogenesis and repress osteoblastogenesis in zebrafish

Triazoles are a major group of azole fungicides commonly used in agriculture, and veterinary and human medicine. Maternal exposure to certain triazole antifungal medication causes congenital malformations, including skeletal malformations. We hypothesized that triazoles used as pesticides in agriculture also pose a risk of causing skeletal malformations in developing embryos. In this study, teratogenic effects of three commonly used triazoles, cyproconazole, paclobutrazol, and triadimenol, were investigated in zebrafish, Danio rerio. Exposure to the triazole fungicides caused bone and cartilage malformations in developing zebrafish larvae. Data from whole-embryo transcriptomics with cyproconazole suggested that exposure to this compound induces adipogenesis while repressing skeletal development. Confirming this finding, the expression of selected bone and cartilage marker genes were significantly downregulated with triazoles exposure as determined by quantitative PCR. The expression of selected adipogenic genes was upregulated by the triazoles. Furthermore, exposure to each of the three triazoles induced adipogenesis and lipid droplet formation in vitro in 3T3-L1 pre-adipocyte cells. In vivo in zebrafish larvae, cyproconazole exposure caused lipid accumulation. These results suggest that exposure to triazoles promotes adipogenesis at the expense of skeletal development, and thus they expand the chemical group of bona fide bone to fat switchers.

Triazole fungicide tebuconazole induces apoptosis through ROS-mediated endoplasmic reticulum stress pathway

Tebuconazole (TEB) is a common triazole fungicide that has been widely applied in the treatment of fungal diseases. It is reported that TEB could exert harmful effects on mammals' health. However, the molecular mechanism involved in TEB toxicity remain undefined. Our study aimed to investigate the mechanisms of TEB-induced toxicity in intestinal cells. We found that TEB stimulates apoptosis through the mitochondrial pathway. Additionally, TEB triggers endoplasmic reticulum (ER) stress as demonstrated by the activation of the three arms of unfolded protein response (UPR). The incubation with the chemical chaperone 4-phenylbutyrate (4-PBA) alleviated ER stress and reduced TEB-induced apoptosis, suggesting that ER stress plays an important role in mediating TEB-induced toxicity. Furthermore, inhibition of ROS by N-acetylcysteine (NAC) inhibited TEB-induced ER stress and apoptosis. Taken together, these findings suggest that TEB exerts its toxic effects in HCT116 cells by inducing apoptosis through ROS-mediated ER stress and mitochondrial apoptotic pathway.

A microRNA or messenger RNA point of departure estimates an apical endpoint point of departure in a rat developmental toxicity model

Traditional developmental toxicity testing practice examines fetal apical endpoints to identify a point of departure (POD) for risk assessment. A potential new testing paradigm involves deriving a POD from a comprehensive analysis of molecular-level change. Here, the rat ketoconazole endocrine-mediated developmental toxicity model was used to test the hypothesis that maternal epigenomic (miRNA) and transcriptomic (mRNA) PODs are similar to fetal apical endpoint PODs. Sprague-Dawley rats were exposed from gestation day (GD) 6-21 to 0, 0.063, 0.2, 0.63, 2, 6.3, 20, or 40 mg/kg/day ketoconazole. Dam systemic, liver, and placenta PODs, along with GD 21 fetal resorption, body weight, and skeletal apical PODs were derived using BMDS software. GD 21 dam liver and placenta miRNA and mRNA PODs were obtained using three methods: a novel individual molecule POD accumulation method, a first mode method, and a gene set method. Dam apical POD values ranged from 2.0 to 38.6 mg/kg/day; the lowest value was for placenta histopathology. Fetal apical POD values were 10.9-20.3 mg/kg/day; the lowest value was for fetal resorption. Dam liver miRNA and mRNA POD values were 0.34-0.69 mg/kg/day, and placenta miRNA and mRNA POD values were 2.53-6.83 mg/kg/day. Epigenomic and transcriptomic POD values were similar across liver and placenta. Deriving a molecular POD from dam liver or placenta was protective of a fetal apical POD. These data support the conclusion that a molecular POD can be used to estimate, or be protective of, a developmental toxicity apical POD.

Grouping of chemicals into mode of action classes by automated effect pattern analysis using the zebrafish embryo toxicity test

A central element of high throughput screens for chemical effect assessment using zebrafish is the assessment and quantification of phenotypic changes. By application of an automated and more unbiased analysis of these changes using image analysis, patterns of phenotypes may be associated with the mode of action (MoA) of the exposure chemical. The aim of our study was to explore to what extent compounds can be grouped according to their anticipated toxicological or pharmacological mode of action using an automated quantitative multi-endpoint zebrafish test. Chemical-response signatures for 30 endpoints, covering phenotypic and functional features, were generated for 25 chemicals assigned to 8 broad MoA classes. Unsupervised clustering of the profiling data demonstrated that chemicals were partially grouped by their main MoA. Analysis with a supervised clustering technique such as a partial least squares discriminant analysis (PLS-DA) allowed to identify markers with a strong potential to discriminate between MoAs such as mandibular arch malformation observed for compounds interfering with retinoic acid signaling. The capacity for discriminating MoAs was also benchmarked to an available battery of in vitro toxicity data obtained from ToxCast library indicating a partially similar performance. Further, we discussed to which extent the collected dataset indicated indeed differences for compounds with presumably similar MoA or whether other factors such as toxicokinetic differences could have an important impact on the determined response patterns.

Transcriptomic changes upon epoxiconazole exposure in a human stem cell-based model of developmental toxicity

Conazole fungicides such as epoxiconazole are mostly used on cereals of crops to inhibit fungal growth through direct inhibition of sterol 14α-demethylase (CYP51A1). However, this enzyme is highly conserved and in humans it is part of the steroid hormone biosynthesis pathway. Endocrine disrupting effects of epoxiconazole have been shown in rodents and have been substantiated by in vitro data, however, the underlying molecular mechanisms are not clear. We took advantage of a human stem cell based in vitro model for developmental toxicity to study the molecular effects of epoxiconazole. This model is based on 3D cultures of embryoid bodies and differentiation into cardiomyocytes, which mimics the early stages of embryonic development. We have previously shown that epoxiconazole impairs differentiation of these embryoid bodies and therefore has the potential to affect human embryonic development. We employed global transcriptome analysis using RNA sequencing and found that the steroid biosynthesis pathway including CYP51A1, the human sterol 14α-demethylase, was highly deregulated by epoxiconazole in our model. We confirmed that most genes of the steroid biosynthesis pathway were upregulated, including CYP51A1, suggesting a compensatory mechanism at the gene expression level. Our data suggest that epoxiconazole acts mainly by decreasing cholesterol biosynthesis in the cells. We conclude that epoxiconazole bears the potential to harm human embryonic development through inhibition of the steroid biosynthesis pathway. As this may be a common feature of compounds that target sterol 14α-demethylase, we add evidence to the assumption that conazole fungicides may be human developmental toxicants.

An adverse outcome pathway on the disruption of retinoic acid metabolism leading to developmental craniofacial defects

Adverse outcome pathway (AOP) is a conceptual framework that links a molecular initiating event (MIE) via intermediate key events (KEs) with adverse effects (adverse outcomes, AO) relevant for risk assessment, through defined KE relationships (KERs). The aim of the present work is to describe a linear AOP, supported by experimental data, for skeletal craniofacial defects as the AO. This AO was selected in view of its relative high incidence in humans and the suspected relation to chemical exposure. We focused on inhibition of CYP26, a retinoic acid (RA) metabolizing enzyme, as MIE, based on robust previously published data. Conazoles were selected as representative stressors. Intermediate KEs are RA disbalance, aberrant HOX gene expression, disrupted specification, migration, and differentiation of neural crest cells, and branchial arch dysmorphology. We described the biological basis of the postulated events and conducted weight of evidence (WoE) assessments. The biological plausibility and the overall empirical evidence were assessed as high and moderate, respectively, the latter taking into consideration the moderate evidence for concordance of dose-response and temporal relationships. Finally, the essentiality assessment of the KEs, considered as high, supported the robustness of the presented AOP. This AOP, which appears of relevance to humans, thus contributes to mechanistic underpinning of selected test methods, thereby supporting their application in integrated new approach test methodologies and strategies and application in a regulatory context.

Multi-level analysis of exposure to triazole fungicides through treated seed ingestion in the red-legged partridge

Triazole fungicides are the most widely used products to treat cereal seeds. Granivorous birds, such as red-legged partridges (Alectoris rufa), which consume seeds left on the surface of fields after sowing, have a high risk of exposure. As triazole fungicides can affect sterol synthesis, we tested the hypothesis that treated seed consumption could alter the synthesis of sex hormones and reduce the reproductive capacity of partridges. We exposed adult partridges to seeds treated with four different formulations containing triazoles as active ingredients (flutriafol, prothioconazole, tebuconazole, and a mixture of the latter two) simulating a field exposure during the late autumn sowing season. All treatments produced biochemical changes and an overexpression of genes encoding for enzymes involved in the biosynthesis of sterols and steroid hormones, such as PMVK, ABCA1, MVD, PSCK9, DHCR7 and HSD17B7. Plasma levels of oestradiol were reduced in partridges exposed to tebuconazole. We also monitored reproduction 3 months after exposure (laying date, egg fertilization and hatching rates). We observed a 14-day delay in the laying onset of partridges that had been exposed to flutriafol as compared to controls. These results show that the consumption of seeds treated with triazole fungicides has the potential to affect granivorous bird reproduction. We recommend the evaluation of lagged reproductive effects as part of the protocols of environmental risk assessment of pesticides in wild birds in light of the effects resulting from the exposure to triazole-treated seeds.

Transcriptome changes of Takifugu obscurus liver after acute exposure to the oxygenated-PAH 9,10-phenanthrenequione

Contamination with polycyclic aromatic hydrocarbons (PAHs) causes noticeable ecological problems in aquatic ecosystems. 9,10-Phenanthrenequione (9,10-PQ) is an oxidized PAH and is highly toxic to aquatic animals. However, the effects of 9,10-PQ on the molecular metabolism of fish remain largely unknown. In this study, Takifugu obscurus juveniles were acutely exposed to 44.30 µg/L 9,10-PQ for 3 days. The transcriptome profile changes in their livers were compared between the 9,10-PQ treatment group and the control using T. rubripes as the reference genome. The results identified 22,414 genes in our transcriptome. Among them, 767 genes were differentially expressed after exposure to 9,10-PQ, which enriched 16 KEGG pathways. Among them, the glycolysis, phagosome, and FOXO signaling pathways were significantly activated in 9,10-PQ treatment compared with the control. These data indicate that 9,10-PQ increased the glycolysis capacity to produce more energy for resistance and harmed immune function. Moreover, several genes related to tumorigenesis were significantly upregulated in response to 9,10-PQ, displaying the carcinogenic toxicity of 9,10-PQ to T. obscurus. Genes in steroid biosynthesis pathways were downregulated in the 9,10-PQ treatment group, suggesting interference with the endocrine system. Overall, these findings provide information to help evaluate the environmental risks that oxygenated-PAHs present to T. obscurus.

The Ascidian Embryo Teratogenicity assay in Ciona intestinalis as a new teratological screening to test the mixture effect of the co-exposure to ethanol and fluconazole

The aim of this work was to evaluate the Ascidian Embryo Teratogenicity assay (AET) as new alternative invertebrate model to test the developmental effects of the co-exposure to ethanol and fluconazole. Ciona intestinalis embryos were exposed to the azolic fungicide fluconazole, (FLUCO, 7.8-250μM), to ethanol (Eth, 0.01-0.5%) and to their mixture (0.01% Eth+FLUCO 7.8-250μM) from neurula to larval stage. At the end of the exposure period, larvae were morphologically evaluated and benchmark analysis performed by using the PROAST modelling software. Both compounds were teratogenic in a concentration-related manner, particularly affecting the pigmented organs. The co-exposure to Eth enhanced the effects of FLUCO, the additive hypothesis was not rejected by the modelling. The results demonstrated that AET could be considered a good vertebrate-free alternative model for toxicological investigation in embryos.

Zebrafish as an Alternative Vertebrate Model for Investigating Developmental Toxicity-The Triadimefon Example

Triadimefon is a widely used triazole fungicide known to cause severe developmental defects in several model organisms and in humans. The present study evaluated in detail the developmental effects seen in zebrafish embryos exposed to triadimefon, confirmed and expanded upon previous phenotypic findings and compared them to those observed in other traditional animal models. In order to do this, we exposed embryos to 2 and 4 µg/mL triadimefon and evaluated growth until 120 h post-fertilization (hpf) through gross morphology examination. Our analysis revealed significant developmental defects at the highest tested concentration including somite deformities, severe craniofacial defects, a cleft phenotype along the three primary neural divisions, a rigorously hypoplastic or even absent mandible and a hypoplastic morphology of the pharyngeal arches. Interestingly, massive pericardial edemas, abnormal shaped hearts, brachycardia and inhibited or absent blood circulation were also observed. Our results revealed that the presented zebrafish phenotypes are comparable to those seen in other organism models and those derived from human observations as a result of triadimefon exposure. We therefore demonstrated that zebrafish provide an excellent system for study of compounds with toxic significance and can be used as an alternative model for developmental toxicity studies to predict effects in mammals.

Fetal safety profile of aromatase inhibitors: Animal data

Aromatase inhibitors (AIs) are a class of drugs that act by blocking the production of estrogens from androgens. The current review concentrates on the prenatal developmental toxicity of AIs in experimental models. Available data indicate that AIs may affect pregnancy at human therapeutic or lower doses. The window of vulnerability to AIs is not limited to organogenesis, but also includes the preimplantation stage and fetal periods. Decreased embryo/fetal survival was the prominent treatment-related effect. Morphological anomalies noted in fetuses exposed to AIs included skeletal anomalies, abnormal head morphology, increased ano-genital distance in female fetuses, and minor urinary tract system anomalies. Placental enlargement was consistently reported in rats and non-human primates after maternal treatment with several AIs. In conclusion, data from basic scientific research suggest that low intensity exposure to AIs applied during a wide gestational window can profoundly affect prenatal development.

Developmental exposures to an azole fungicide triadimenol at environmentally relevant concentrations cause reproductive dysfunction in females of medaka fish

Triadimenol is an effective meatabolite derived from the triazole fungicide triadimenfon. It is an agriculturally important reagent of environmentally emerging concern because of its broad use, persistent occurrence in the environment and greater fungicidal or toxic potency than the parent compound. However, the ecotoxicological impact of triadimenol on fish populations remains unclear. In this study, we investigated developmental toxicity and endocrine disruption effects in medaka fish (Oryzias latipes) exposed at an early life stage to triadimenol. First, mortality, gross development and oxidative stress responses were assessed with triadimenol exposure (3-3000 μg/L) during the embryonic stage. Then, medaka at a sensitive stage of early sexual development underwent 35-day continuous chronic exposure to triadimenol, and the endocrine disruption effects were assessed in adulthood and the next generation. Embryonic exposure to triadimenol did not induce significant teratogenic effects or oxidative stress in embryos or hatchlings. However, early-life exposure to triadimenol under environmentally relevant concentrations (3-30 μg/L) and 300 μg/L persistently altered ovary development and reproduction in female adults and skewed the sex ratio in progeny. As well, triadimenol exposure interrupted the hormone balance, as seen by the expression of genes responsible for estrogen metabolism and egg reproduction. Environmentally relevant triadimenol exposure in medaka fish at early life stages may have ecotoxicological impact in aquatic environments. Along with previous studies, we suggest that conazoles share similar modes of action in disrupting hormone homeostasis and reproduction in fish and mammals.

Extended evaluation on the ES-D3 cell differentiation assay combined with the BeWo transport model, to predict relative developmental toxicity of triazole compounds

The mouse embryonic stem D3 (ES-D3) cell differentiation assay is based on the morphometric measurement of cardiomyocyte differentiation and is a promising tool to detect developmental toxicity of compounds. The BeWo transport model, consisting of BeWo b30 cells grown on transwell inserts and mimicking the placental barrier, is useful to determine relative placental transport velocities of compounds. We have previously demonstrated the usefulness of the ES-D3 cell differentiation assay in combination with the in vitro BeWo transport model to predict the relative in vivo developmental toxicity potencies of a set of reference azole compounds. To further evaluate this combined in vitro toxicokinetic and toxicodynamic approach, we combined ES-D3 cell differentiation data of six novel triazoles with relative transport rates obtained from the BeWo model and compared the obtained ranking to the developmental toxicity ranking as derived from in vivo data. The data show that the combined in vitro approach provided a correct prediction for in vivo developmental toxicity, whereas the ES-D3 cell differentiation assay as stand-alone did not. In conclusion, we have validated the combined in vitro approach for developmental toxicity, which we have previously developed with a set of reference azoles, for a set of six novel triazoles. We suggest that this combined model, which takes both toxicodynamic and toxicokinetic aspects into account, should be further validated for other chemical classes of developmental toxicants.

Biomonitoring of the Genotoxic and Hepatotoxic Effects and Oxidative Stress Potentials of Itraconazole in Pregnant Rats

BACKGROUND

Pregnant women are more susceptible to both vaginal colonization and infection by yeast. One hundred million fungal infected patients have been treated worldwide with itraconazole (Caputo, 2003.

METHOD

Itraconazole was administrated orally to pregnant rats at doses of 75, 100, or 150 mg/kg during gestational days (GD) 1 to 7 or GD 8 to 14 or GD 14 to 20. The genotoxicity and hepatotoxicity of the antifungal drug itraconazole were assessed during different periods of pregnancy using different methods.

RESULTS

It was found that itraconazole was a genotoxic drug for both mothers and fetuses. This finding was observed via significant elevation in the estimated comet assay parameters (percentage of fragmented DNA, tail moment, and olive moment), percentage of fragmented DNA measured by diphenylamine assay and mixed smearing and laddering of DNA fragments of liver samples. In addition, itraconazole caused significant elevation in the level of hepatic malondialdehyde and depletion in the catalase activity and glutathione level. Furthermore, itraconazole induced histopathological alterations in the hepatic tissues of both mothers and fetuses.

CONCLUSION

These findings indicate that itraconazole administration at doses of 75, 100, or 150 mg/kg during pregnancy induced maternal and fetal toxicity that could be induced by the genotoxicity and the oxidative damage.

Triazole-induced toxicity in developing rare minnow (Gobiocypris rarus) embryos

Using rare minnow (Gobiocypris rarus) at early-life stages as experimental models, the developmental toxicity of five widely used triazole fungicides (myclobutanil, fluconazole, flusilazole, triflumizole, and epoxiconazole) were investigated following exposure to 1-15 mg/L for 72 h. Meanwhile, morphological parameters (body length, body weight, and heart rate), enzyme activities (superoxide dismutase (SOD), glutathione S-transferase (GST), adenosine triphosphatase (ATPase), and acetyl cholinesterase (AChE)), and mRNA levels (hsp70, mstn, mt, apaf1, vezf1, and cyp1a) were also recorded following exposure to 0.2, 1.0, and 5.0 mg/L for 72 h. Results indicated that increased malformation and mortality, decreased body length, body weight, and heart rate provide a concentration-dependent pattern; values of 72 h LC50 (median lethal concentration) and EC50 (median effective concentration) ranged from 3 to 12 mg/L. Most importantly, the results of the present study suggest that even at the lowest concentration, 0.2 mg/L, five triazole fungicides also caused notable changes in enzyme activities and mRNA levels. Overall, the present study points out that those five triazole fungicides are highly toxic to the early development of G. rarus embryos. The information presented in this study will be helpful in better understanding the toxicity induced by triazole fungicides in fish embryos.

Lipid peroxidation, detoxification capacity, and genome damage in mice after transplacental exposure to pharmaceutical drugs

Data on genome damage, lipid peroxidation, and levels of glutathione peroxidase (GPX) in newborns after transplacental exposure to xenobiotics are rare and insufficient for risk assessment. The aim of the current study was to analyze, in an animal model, transplacental genotoxicity, lipid peroxidation, and detoxification disturbances caused by the following drugs commonly prescribed to pregnant women: paracetamol, fluconazole, 5-nitrofurantoin, and sodium valproate. Genome damage in dams and their newborn pups transplacentally exposed to these drugs was investigated using the in vivo micronucleus (MN) assay. The drugs were administered to dams intraperitoneally in three consecutive daily doses between days 12 and 14 of pregnancy. The results were correlated, with detoxification capacity of the newborn pups measured by the levels of GPX in blood and lipid peroxidation in liver measured by malondialdehyde (HPLC-MDA) levels. Sodium valproate and 5-nitrofurantoin significantly increased MN frequency in pregnant dams. A significant increase in the MN frequency of newborn pups was detected for all drugs tested. This paper also provides reference levels of MDA in newborn pups, according to which all drugs tested significantly lowered MDA levels of newborn pups, while blood GPX activity dropped significantly only after exposure to paracetamol. The GPX reduction reflected systemic oxidative stress, which is known to occur with paracetamol treatment. The reduction of MDA in the liver is suggested to be an unspecific metabolic reaction to the drugs that express cytotoxic, in particular hepatotoxic, effects associated with oxidative stress and lipid peroxidation.

Use of oral fluconazole during pregnancy and the risk of birth defects

BACKGROUND

Case reports suggest that long-term, high-dose fluconazole treatment for severe fungal infections during pregnancy causes a pattern of birth defects. It is unclear whether commonly used lower doses increase the risk of specific birth defects.

METHODS

In a registry-based cohort of liveborn infants in Denmark, we evaluated first-trimester oral fluconazole exposure and the risk of birth defects overall and of birth defects previously linked to azole antifungal agents.

RESULTS

The majority of fluconazole-exposed pregnancies were in women who received common therapeutic doses of 150 mg (56% of pregnancies) or 300 mg (31%). Oral fluconazole exposure was not associated with an increased risk of birth defects overall (210 birth defects among 7352 fluconazole-exposed pregnancies [prevalence, 2.86%] and 25,159 birth defects among 968,236 unexposed pregnancies [prevalence, 2.60%]; adjusted prevalence odds ratio, 1.06; 95% confidence interval [CI], 0.92 to 1.21). In addition, oral fluconazole exposure was not associated with a significantly increased risk of 14 of 15 types of birth defects previously linked to azole antifungal agents: craniosynostosis, other craniofacial defects, middle-ear defects, cleft palate, cleft lip, limb defects, limb-reduction defects, polydactyly, syndactyly, diaphragmatic hernia, heart defects overall, pulmonary-artery hypoplasia, ventricular septal defects, and hypoplastic left heart. A significantly increased risk of tetralogy of Fallot was observed (7 cases in fluconazole-exposed pregnancies [prevalence, 0.10%] as compared with 287 cases in unexposed pregnancies [prevalence, 0.03%]; adjusted prevalence odds ratio, 3.16; 95% CI, 1.49 to 6.71).

CONCLUSIONS

Oral fluconazole was not associated with a significantly increased risk of birth defects overall or of 14 of the 15 specific birth defects of previous concern. Fluconazole exposure may confer an increased risk of tetralogy of Fallot. (Funded by the Danish Medical Research Council.).

An optimized gene set for transcriptomics based neurodevelopmental toxicity prediction in the neural embryonic stem cell test

The murine neural embryonic stem cell test (ESTn) is an in vitro model for neurodevelopmental toxicity testing. Recent studies have shown that application of transcriptomics analyses in the ESTn is useful for obtaining more accurate predictions as well as mechanistic insights. Gene expression responses due to stem cell neural differentiation versus toxicant exposure could be distinguished using the Principal Component Analysis based differentiation track algorithm. In this study, we performed a de novo analysis on combined raw data (10 compounds, 19 exposures) from three previous transcriptomics studies to identify an optimized gene set for neurodevelopmental toxicity prediction in the ESTn. By evaluating predictions of 200,000 randomly selected gene sets, we identified genes which significantly contributed to the prediction reliability. A set of 100 genes was obtained, predominantly involved in (neural) development. Further stringency restrictions resulted in a set of 29 genes that allowed for 84% prediction accuracy (area under the curve 94%). We anticipate these gene sets will contribute to further improve ESTn transcriptomics studies aimed at compound risk assessment.

Compound-specific effects of diverse neurodevelopmental toxicants on global gene expression in the neural embryonic stem cell test (ESTn)

Alternative assays for developmental toxicity testing are needed to reduce animal use in regulatory toxicology. The in vitro murine neural embryonic stem cell test (ESTn) was designed as an alternative for neurodevelopmental toxicity testing. The integration of toxicogenomic-based approaches may further increase predictivity as well as provide insight into underlying mechanisms of developmental toxicity. In the present study, we investigated concentration-dependent effects of six mechanistically diverse compounds, acetaldehyde (ACE), carbamazepine (CBZ), flusilazole (FLU), monoethylhexyl phthalate (MEHP), penicillin G (PENG) and phenytoin (PHE), on the transcriptome and neural differentiation in the ESTn. All compounds with the exception of PENG altered ESTn morphology (cytotoxicity and neural differentiation) in a concentration-dependent manner. Compound induced gene expression changes and corresponding enriched gene ontology biological processes (GO-BP) were identified after 24h exposure at equipotent differentiation-inhibiting concentrations of the compounds. Both compound-specific and common gene expression changes were observed between subsets of tested compounds, in terms of significance, magnitude of regulation and functionality. For example, ACE, CBZ and FLU induced robust changes in number of significantly altered genes (≥ 687 genes) as well as a variety of GO-BP, as compared to MEHP, PHE and PENG (≤ 55 genes with no significant changes in GO-BP observed). Genes associated with developmentally related processes (embryonic morphogenesis, neuron differentiation, and Wnt signaling) showed diverse regulation after exposure to ACE, CBZ and FLU. In addition, gene expression and GO-BP enrichment showed concentration dependence, allowing discrimination of non-toxic versus toxic concentrations on the basis of transcriptomics. This information may be used to define adaptive versus toxic responses at the transcriptome level.

Comparison of the mouse Embryonic Stem cell Test, the rat Whole Embryo Culture and the Zebrafish Embryotoxicity Test as alternative methods for developmental toxicity testing of six 1,2,4-triazoles

The relatively high experimental animal use in developmental toxicity testing has stimulated the search for alternatives that are less animal intensive. Three widely studied alternative assays are the mouse Embryonic Stem cell Test (EST), the Zebrafish Embryotoxicity Test (ZET) and the rat postimplantation Whole Embryo Culture (WEC). The goal of this study was to determine their efficacy in assessing the relative developmental toxicity of six 1,2,4-triazole compounds,(1) flusilazole, hexaconazole, cyproconazole, triadimefon, myclobutanil and triticonazole. For this purpose, we analyzed effects and relative potencies of the compounds in and among the alternative assays and compared the findings to their known in vivo developmental toxicity. Triazoles are antifungal agents used in agriculture and medicine, some of which are known to induce craniofacial and limb abnormalities in rodents. The WEC showed a general pattern of teratogenic effects, typical of exposure to triazoles, mainly consisting of reduction and fusion of the first and second branchial arches, which are in accordance with the craniofacial malformations reported after in vivo exposure. In the EST all triazole compounds inhibited cardiomyocyte differentiation concentration-dependently. Overall, the ZET gave the best correlation with the relative in vivo developmental toxicities of the tested compounds, closely followed by the EST. The relative potencies observed in the WEC showed the lowest correlation with the in vivo developmental toxicity data. These differences in the efficacy between the test systems might be due to differences in compound kinetics, in developmental stages represented and in the relative complexity of the alternative assays.

Gene set assembly for quantitative prediction of developmental toxicity in the embryonic stem cell test

The embryonic stem cell test (EST) is an in vitro method for predicting developmental toxicity based on compound-induced inhibition of embryonic stem cell (ESC) differentiation. We previously described how gene expression analysis in the EST can be used to describe normal ESC differentiation as well as identify compound developmental toxicity, by means of our differentiation track algorithm. In this study, we combined raw data from our three previous studies in a new integrated analysis, to identify a gene set that allows for improved prediction. By evaluating predictions of 100,000 randomly selected gene sets, we identified which genes contribute significantly to the prediction reliability. By additional cross-validation, we identified a set of 52 genes that allows for improved prediction of toxicity. The correlation between the predictions using this gene set and the magnitude of the EST endpoint was 0.85, and the gene set predicted developmental toxicity with 83% accuracy (area under the curve 89%). If compounds with ineffective data because of a too low tested concentration or too much variation between samples were excluded, even 100% accuracy could be reached based on 15 compounds. This novel gene set consists mainly of genes involved in the stem cell differentiation or other developmental processes. We expect that this set can be of use in future studies aimed at improving the EST for risk assessment, thus making a next step towards regulatory implementation of this method.

Discriminating classes of developmental toxicants using gene expression profiling in the embryonic stem cell test

The embryonic stem cell test (EST) has been shown to be a promising in vitro method for the prediction of developmental toxicity. In our previous studies, we demonstrated that the implementation of gene expression analysis in the EST may further improve the identification of developmental toxicants. In the present study, we investigated if gene expression profiling could be used to discriminate compound classes with distinct modes of action (MoA) using the EST protocol. Gene expression data of our previous study were used and were analyzed of embryonic stem cell (ESC) differentiation cultures exposed to six compounds belonging to two classes with distinct MoA, namely phthalates and triazoles. We used three approaches to study class-characteristic gene regulation that may be useful for discrimination of compound classes. First, at the individual gene level, gene signatures characteristic for each class were identified that successfully discriminated both classes using principal component analysis. Second, at the functional level, enriched gene ontology (GO) biological processes showed their usefulness for class discrimination via hierarchical clustering. Third, two previously identified gene sets, which we designed to predict developmental toxicity, appeared successful in separating phthalate from triazole compounds. In summary, we established the possibility to discriminate between compound classes in the EST system using three different specific transcriptomics-based approaches. Differential gene expression information specific for the class of compound under study may be employed to optimize prioritization of compounds within that class for further testing.