A predictive data-driven framework for endocrine prioritization: a triazole fungicide case study.

The US Environmental Protection Agency Endocrine Disruptor Screening Program (EDSP) is a tiered screening approach to determine the potential for a chemical to interact with estrogen, androgen, or thyroid hormone systems and/or perturb steroidogenesis. Use of high-throughput screening (HTS) to predict hazard and exposure is shifting the EDSP approach to (1) prioritization of chemicals for further screening; and (2) targeted use of EDSP Tier 1 assays to inform specific data needs. In this work, toxicology data for three triazole fungicides (triadimefon, propiconazole, and myclobutanil) were evaluated, including HTS results, EDSP Tier 1 screening (and other scientifically relevant information), and EPA guideline mammalian toxicology study data. The endocrine-related bioactivity predictions from HTS and information that satisfied the EDSP Tier 1 requirements were qualitatively concordant. Current limitations in the available HTS battery for thyroid and steroidogenesis pathways were mitigated by inclusion of guideline toxicology studies in this analysis. Similar margins (3-5 orders of magnitude) were observed between HTS-predicted human bioactivity and exposure values and between in vivo mammalian bioactivity and EPA chronic human exposure estimates for these products' registered uses. Combined HTS hazard and human exposure predictions suggest low priority for higher-tiered endocrine testing of these triazoles. Comparison with the mammalian toxicology database indicated that this HTS-based prioritization would have been protective for any potential in vivo effects that form the basis of current risk assessment for these chemicals. This example demonstrates an effective, human health protective roadmap for EDSP evaluation of pesticide active ingredients via prioritization using HTS and guideline toxicology information.

Phenobarbital and propiconazole toxicogenomic profiles in mice show major similarities consistent with the key role that constitutive androstane receptor (CAR) activation plays in their mode of action.

Toxicogenomics (TGx) is employed frequently to investigate underlying molecular mechanisms of the compound of interest and, thus, has become an aid to mode of action determination. However, the results and interpretation of a TGx dataset are influenced by the experimental design and methods of analysis employed. This article describes an evaluation and reanalysis, by two independent laboratories, of previously published TGx mouse liver microarray data for a triazole fungicide, propiconazole (PPZ), and the anticonvulsant drug phenobarbital (PB). Propiconazole produced an increase incidence of liver tumors in male CD-1 mice only at a dose that exceeded the maximum tolerated dose (2500 ppm). Firstly, we illustrate how experimental design differences between two in vivo studies with PPZ and PB may impact the comparisons of TGx results. Secondly, we demonstrate that different researchers using different pathway analysis tools can come to different conclusions on specific mechanistic pathways, even when using the same datasets. Finally, despite these differences the results across three different analyses also show a striking degree of similarity observed for PPZ and PB treated livers when the expression data are viewed as major signaling pathways and cell processes affected. Additional studies described here show that the postulated key event of hepatocellular proliferation was observed in CD-1 mice for both PPZ and PB, and that PPZ is also a potent activator of the mouse CAR nuclear receptor. Thus, with regard to the events which are hallmarks of CAR-induced effects that are key events in the mode of action (MOA) of mouse liver carcinogenesis with PB, PPZ-induced tumors can be viewed as being promoted by a similar PB-like CAR-dependent MOA.

Current and future use of genomics data in toxicology: opportunities and challenges for regulatory applications.

Toxicogenomics is the application of toxicology, genetics, molecular biology and environmental health to describe the response of organisms to environmental stimuli. The field of toxicogenomics has developed over the past 15 years mainly due to advances in toxicology, molecular genetics and cell biology. Its prospective use to resolve crucial data gaps and data inconsistencies could improve risk assessment by providing additional data to increase the understanding of mechanisms and modes of action (MOA) and enhance the reliability of dose-response extrapolation. Thus, toxicogenomics holds promise for advancing the scientific basis of risk assessments. However, one of the current issues is how genomic/transcriptional data is being used to further describe a MOA for oncogenicity and, in turn, its potential uses in cancer risk assessment. This commentary identifies how toxicogenomics could be used on a case by case basis to add information to a MOA addressing both the opportunities and challenges this technology holds. In addition, some pitfalls to avoid in the generation and interpretation of toxicogenomic data and validation issues that need to be addressed before toxicogenomics can be used in the risk assessment process and regulatory decisions are discussed.

Inhibition of rat and human steroidogenesis by triazole antifungals.

Environmental chemicals that alter steroid production could interfere with male reproductive development and function. Three agricultural antifungal triazoles that are known to modulate expression of cytochrome P450 (CYP) genes and enzymatic activities were tested for effects on steroidogenesis using rat in vivo (triadimefon), rat in vitro (myclobutanil and triadimefon), and human in vitro (myclobutanil, propiconazole, and triadimefon) model systems. Hormone production was measured in testis organ cultures from untreated adult and neonatal rats, following in vitro exposure to 1, 10, or 100 muM of myclobutanil or triadimefon. Myclobutanil and triadimefon reduced media levels of testosterone by 40-68% in the adult and neonatal testis culture, and altered steroid production in a manner that indicated CYP17-hydroxylase/17,20 lyase (CYP17A1) inhibition at the highest concentration tested. Rat to human comparison was explored using the H295R (human adrenal adenocarcinoma) cell line. Following 48 h exposure to myclobutanil, propiconazole, or triadimefon at 1, 3, 10, 30, or 100 muM, there was an overall decrease in estradiol, progesterone, and testosterone by all three triazoles. These data indicate that myclobutanil, propiconazole, and triadimefon are weak inhibitors of testosterone production in vitro. However, in vivo exposure of rats to triazoles resulted in increased serum and intra-testicular testosterone levels. This discordance could be due to higher concentrations of triazoles tested in vitro, and differences within an in vitro model system lacking hepatic metabolism and neuroendocrine control.

Mode of action for reproductive and hepatic toxicity inferred from a genomic study of triazole antifungals.

The mode of action for the reproductive toxicity of some triazole antifungals has been characterized as an increase in serum testosterone and hepatic response, and reduced insemination and fertility indices. In order to refine our mechanistic understanding of these potential modes of action, gene expression profiling was conducted on liver and testis from male Wistar Han IGS rats exposed to myclobutanil (500, 2000 ppm), propiconazole (500, 2500 ppm), or triadimefon (500, 1800 ppm) from gestation day six to postnatal day 92. Gene expression profiles indicated that all three triazoles significantly perturbed the fatty acid, steroid, and xenobiotic metabolism pathways in the male rat liver. In addition, triadimefon modulated expression of genes in the liver from the sterol biosynthesis pathway. Although expression of individual genes were affected, there were no common pathways modulated by all three triazoles in the testis. The pathways identified in the liver included numerous genes involved in phase I-III metabolism (Aldh1a1, Cyp1a1, Cyp2b2, Cyp3a1, Cyp3a2, Slco1a4, Udpgtr2), fatty acid metabolism (Cyp4a10, Pcx, Ppap2b), and steroid metabolism (Ugt1a1, Ugt2a1) for which expression was altered by the triazoles. These differentially expressed genes form part of a network involving lipid, sterol, and steroid homeostatic pathways regulated by the constitutive androstane (CAR), pregnane X (PXR), peroxisome proliferator-activated alpha, and other nuclear receptors in liver. These relatively high dose and long-term exposures to triazole antifungals appeared to perturb fatty acid and steroid metabolism in the male rat liver predominantly through the CAR and PXR signaling pathways. These toxicogenomic effects describe a plausible series of key events contributing to the disruption in steroid homeostasis and reproductive toxicity of select triazole antifungals.

Toxicogenomic effects common to triazole antifungals and conserved between rats and humans.

The triazole antifungals myclobutanil, propiconazole and triadimefon cause varying degrees of hepatic toxicity and disrupt steroid hormone homeostasis in rodent in vivo models. To identify biological pathways consistently modulated across multiple timepoints and various study designs, gene expression profiling was conducted on rat livers from three separate studies with triazole treatment groups ranging from 6 h after a single oral gavage exposure, to prenatal to adult exposures via feed. To explore conservation of responses across species, gene expression from the rat liver studies were compared to in vitro data from rat and human primary hepatocytes exposed to the triazoles. Toxicogenomic data on triazoles from 33 different treatment groups and 135 samples (microarrays) identified thousands of probe sets and dozens of pathways differentially expressed across time, dose, and species--many of these were common to all three triazoles, or conserved between rodents and humans. Common and conserved pathways included androgen and estrogen metabolism, xenobiotic metabolism signaling through CAR and PXR, and CYP mediated metabolism. Differentially expressed genes included the Phase I xenobiotic, fatty acid, sterol and steroid metabolism genes Cyp2b2 and CYP2B6, Cyp3a1 and CYP3A4, and Cyp4a22 and CYP4A11; Phase II conjugation enzyme genes Ugt1a1 and UGT1A1; and Phase III ABC transporter genes Abcb1 and ABCB1. Gene expression changes caused by all three triazoles in liver and hepatocytes were concentrated in biological pathways regulating lipid, sterol and steroid homeostasis, identifying a potential common mode of action conserved between rodents and humans. Modulation of hepatic sterol and steroid metabolism is a plausible mode of action for changes in serum testosterone and adverse reproductive outcomes observed in rat studies, and may be relevant to human risk assessment.

Disruption of testosterone homeostasis as a mode of action for the reproductive toxicity of triazole fungicides in the male rat.

Triazole fungicides associated with a range of reported male reproductive effects in experimental animals were selected to assess potential toxic modes of action. Wistar Han rats were fed myclobutanil (M: 100, 500, or 2000 ppm), propiconazole (P: 100, 500, or 2500 ppm), or triadimefon (T: 100, 500, or 1800 ppm) from gestation day 6 to postnatal day (PND) 120. One male per litter was necropsied on PND1, 22, 50, or 92. Measurements included anogenital distance (AGD) at PND0, body and organ weights, serum hormone levels, age at preputial separation (PPS), sperm morphology and motility, and fertility and fecundity. AGD was increased by the high dose of all three triazoles, indicating hypervirilization. Triadimefon delayed PPS, consistent with delayed puberty, at 1800 ppm. Relative liver weights were increased at PND1, 50, and 92 by all three triazoles. Hepatocellular hypertrophy was present at PND50 from propiconazole and triadimefon and at PND92 from all three high-dose triazole treatments. Relative pituitary weights were decreased at PND92 by middle- and high-dose myclobutanil treatment. Absolute testis weights were increased at PND1 by myclobutanil, at PND22 by myclobutanil and triadimefon, and at PND50 by propiconazole and triadimefon treatment. Relative ventral prostate weights were increased at PND92 by myclobutanil and triadimefon treatment. Serum testosterone was increased at PND50 by triadimefon and at PND92/99 by all three triazole treatments. Insemination and fertility were impaired by myclobutanil and triadimefon treatment. In addition to the reproductive system effects, total serum thyroxine levels were decreased at PND92 by high-dose triadimefon. These reproductive effects are consistent with the disruption of testosterone homeostasis as a key event in the mode of action for triazole-induced reproductive toxicity.

Effect of conazole fungicides on reproductive development in the female rat.

Three triazole fungicides were evaluated for effects on female rat reproductive development. Rats were exposed via feed to propiconazole (P) (100, 500, or 2500 ppm), myclobutanil (M) (100, 500, or 2000 ppm), or triadimefon (T) (100, 500, or 1800 ppm) from gestation day 6 to postnatal day (PND) 98. Body weight (BW) and anogenital distance (AGD) at PND 0, age and BW at vaginal opening (VO), estrous cyclicity, and body and organ weight at necropsy were measured. BW at PND 0 was unaffected by treatment. AGD was increased by M2000. VO was delayed by M2000 and T1800. Estrous cyclicity was initially disrupted by P500, P2500 and T1800, but later normalized. At PND 99 there was a decrease in BW by T1800, an increase in liver weight by P2500 and T1800, and an increase in ovarian weight by M2000 and T1800. It is concluded that exposure to P, M and T adversely impacted female rodent reproductive development.

Gene expression profiling in the liver of CD-1 mice to characterize the hepatotoxicity of triazole fungicides.

Four triazole fungicides used in agricultural or pharmaceutical applications were examined for hepatotoxic effects in mouse liver. Besides organ weight, histopathology, and cytochrome P450 (CYP) enzyme induction, DNA microarrays were used to generate gene expression profiles and hypotheses on potential mechanisms of action for this class of chemicals. Adult male CD-1 mice were exposed daily for 14 days to fluconazole, myclobutanil, propiconazole, or triadimefon at three dose levels by oral gavage. Doses were based on previous studies that resulted in liver hypertrophy or hepatotoxicity. All four triazoles caused hepatocyte hypertrophy, and all except triadimefon increased relative liver/body weight ratios at the middle and high dose levels. CYP enzyme activities were also induced by all four triazoles at the middle and high doses as measured by the dealkylations of four alkoxyresorufins, although some differences in substrate specificity were observed. Consistent with this common histopathology and biochemistry, several CYP and xenobiotic metabolizing enzyme (XME) genes were differentially expressed in response to all four (Cyp2d26 and Cyp3a11), or three of the four (Cyp2c40, Cyp2c55, Ces2, Slco1a4) triazoles. Differential expression of numerous other CYP and XME genes discriminated between the various triazoles, consistent with differences in CYP enzyme activities, and indicative of possible differences in mechanisms of hepatotoxicity or dose response. Multiple isoforms of Cyp1a, 2b, 2c, 3a, and other CYP and XME genes regulated by the nuclear receptors constitutive androstane receptor (CAR) and pregnane X receptor (PXR) were differentially expressed following triazole exposure. Based on these results, we expanded on our original hypothesis that triazole hepatotoxicity was mediated by CYP induction, to include additional XME genes, many of which are modulated by CAR and PXR.

Gene expression profiling in liver and testis of rats to characterize the toxicity of triazole fungicides.

Four triazole fungicides were studied using toxicogenomic techniques to identify potential mechanisms of action. Adult male Sprague-Dawley rats were dosed for 14 days by gavage with fluconazole, myclobutanil, propiconazole, or triadimefon. Following exposure, serum was collected for hormone measurements, and liver and testes were collected for histology, enzyme biochemistry, or gene expression profiling. Body and testis weights were unaffected, but liver weights were significantly increased by all four triazoles, and hepatocytes exhibited centrilobular hypertrophy. Myclobutanil exposure increased serum testosterone and decreased sperm motility, but no treatment-related testis histopathology was observed. We hypothesized that gene expression profiles would identify potential mechanisms of toxicity and used DNA microarrays and quantitative real-time PCR (qPCR) to generate profiles. Triazole fungicides are designed to inhibit fungal cytochrome P450 (CYP) 51 enzyme but can also modulate the expression and function of mammalian CYP genes and enzymes. Triazoles affected the expression of numerous CYP genes in rat liver and testis, including multiple Cyp2c and Cyp3a isoforms as well as other xenobiotic metabolizing enzyme (XME) and transporter genes. For some genes, such as Ces2 and Udpgtr2, all four triazoles had similar effects on expression, suggesting possible common mechanisms of action. Many of these CYP, XME and transporter genes are regulated by xeno-sensing nuclear receptors, and hierarchical clustering of CAR/PXR-regulated genes demonstrated the similarities of toxicogenomic responses in liver between all four triazoles and in testis between myclobutanil and triadimefon. Triazoles also affected expression of multiple genes involved in steroid hormone metabolism in the two tissues. Thus, gene expression profiles helped identify possible toxicological mechanisms of the triazole fungicides.

Fluconazole-induced hepatic cytochrome P450 gene expression and enzymatic activities in rats and mice.

This study was undertaken to examine the effects of the triazole antifungal agent fluconazole on the expression of hepatic cytochrome P450 (Cyp) genes and the activities of Cyp enzymes in male Sprague-Dawley rats and male CD-1 mice. Alkoxyresorufin O-dealkylation (AROD) methods were used as measures of Cyp enzyme activities. Western analyses identified specific Cyp isoforms. Quantitative real-time reverse-transcription polymerase chain reaction (quantitative real time-RT-PCR) assays were used to quantitate the mRNA expression of specific Cyp genes induced by this conazole. Rats and mice were administered fluconazole 2, 25, or 50 mg/kg bw/d by gavage daily for 14 days. In rats, fluconazole treatment (50 mg/kg bw/d) significantly induced pentoxyresorufin O-dealkylation (PROD), benzyloxyresorufin O-dealkylation (BROD), and ethoxyresorufin O-dealkylation (EROD) hepatic microsomal activities. Fluconazole treatment significantly increased rat hepatic mRNA expression of CYP2B1 and CYP3A23/3A1 with dose-related responses. The highest dose of fluconazole gave a 128-fold induction of CYP2B1 and a 4.6-fold induction of CYP3A23/3A1 mRNA. CYP3A2 mRNA levels were also overexpressed 5.6-7.2-fold depending on dose. Western immunoblots of rat hepatic microsomal proteins identified Cyp isoforms: CYP1A1, CYP1A2, CYP2B1/2, CYP3A23/3A1, and Cyp3A2 with increased levels of CYP2B1/2 and CYP3A23/3A1 proteins. In mice, fluconazole induced BROD, PROD, EROD, and methoxyresorufin O-dealkylation hepatic microsomal activities after treatment with 25 and 50 mg/kg bw/d. Fluconazole increased mouse hepatic mRNA expression of Cyp2b10 (1.9-fold) and Cyp3a11 (2.6-fold) in the 50 mg/kg bw/d treatment group. In summary, these results indicated that fluconazole, a triazole-containing conazole, clearly induced CYP2B and CYP3A families of isoforms in rat liver and Cyp2b and Cyp3a families of isoforms in mouse liver.

A database for tracking toxicogenomic samples and procedures.

Reproductive toxicogenomic studies generate large amounts of toxicological and genomic data. On the toxicology side, a substantial quantity of data accumulates from conventional endpoints such as histology, reproductive physiology and biochemistry. The largest source of genomics data is DNA microarrays, which generate enormous amounts of information in the course of profiling gene expression. Thus, data storage and management become essential and require a more sophisticated system than lab notebooks and electronic spreadsheets. We developed a database for tracking toxicogenomic samples and procedures (TSP 1.0) for our reproductive studies based on the MIAME-Tox guidelines and relational database theory. This database stores the various types of data from both toxicological and genomic assays in a hierarchical fashion. The user-friendly interface provides easy procedures for researchers to add, edit, save, delete, and navigate different records. Finally, TSP facilitates exporting microarray data into public databases.

Propiconazole-induced cytochrome P450 gene expression and enzymatic activities in rat and mouse liver.

Propiconazole is a N-substituted triazole used as a fungicide on fruits, grains, seeds, hardwoods, and conifers. In the present study, propiconazole was examined for its effects on the expression of hepatic cytochrome P450 genes and on the activities of P450 enzymes in male Sprague-Dawley rats and male CD-1 mice. Rats and mice were administered propiconazole by gavage daily for 14 days at doses of 10, 75, and 150 mg/kg body weight/day. Quantitative real time RT-PCR assays of rat hepatic RNA samples from animals treated at the 150 mg/kg body weight/day dose revealed significant mRNA overexpression of the following genes compared to control: CYP1A2 (1.62-fold), CYP2B1 (10.8-fold), CYP3A1/CYP3A23 (2.78-fold), and CYP3A2 (1.84-fold). In mouse liver, propiconazole produced mRNA overexpression of Cyp2b10 (2.39-fold) and Cyp3a11 (5.19-fold). mRNA expression of CYP1A1 was not detected in liver tissues from treated or controls animals from either species. Propiconazole significantly induced both pentoxyresorufin O-dealkylation (PROD) and methoxyresorufin O-dealkylation (MROD) activities in both rat and mouse liver at the 150 mg/kg body weight/day and 75 mg/kg body weight/day doses. In summary, these results indicated that propiconazole induced CYP1A2 in rat liver and CYP2B and CYP3A families of isoforms in rat and mouse liver.