Alternative developmental toxicity models for assessing the in vivo embryotoxicity of azoles

The implementation of regulations for protecting both humans and the environment from potential chemical health hazards, as well as the increase of global pressure for reducing, refining and replacing animal experiments promote the development and application of alternatives to in vivo developmental toxicity studies. Due to the complexity of the reproductive cycle, combined in vitro approaches, focusing on morphological, molecular and toxicokinetic parameters, could better define the developmental toxicity of chemicals. In this thesis, azoles, which are a group of chemicals with antifungal activity, are under investigation. These compounds show marked differences in developmental toxicity potency and similarities with retinoic acid (RA)- related teratogenicity. Chapter 1 of this thesis introduced information regarding the background of reproductive and developmental toxicology, including scientific concerns and the impact of past teratogenic outcomes on the society. For screening developmental teratogens, in vitro approaches have been proposed and successfully applied. Their combination may better mimic the in vivo embryo and, therefore, increase the accuracy in predicting possible developmental toxicants. Additional co-implementation of molecular approaches may give an insight in the mode of action underlying the observed effects. Azoles were selected in the present thesis due to evidence for possibly increasing developmental toxicity through dysregulating the balance of the RA pathways in the mammalian system. The chapter also described the objectives and outline of the research. In chapter 2, we examined the time- dependent developmental effects in rat embryos exposed in vitro to flusilazole (FLU), and their link to RA mediated pathways. To this end, we assessed the effects of 4-hour exposure of whole embryo culture (WEC) embryos to 300μM FLU during four developmental time windows (0-4, 4-8, 24-28 and 44-48 h), evaluated morphological parameters, as well as expression and localization of five genes directly or indirectly linked with the RA pathway. A stage- specific gene expression response of cultured rat embryos exposed to FLU was detected, which preceded the development of morphologically observable malformations. During all the tested time windows, the most pronounced effect was observed in the regulation of RA-related genes. Therefore, it was concluded that such biomarkers can be employed as useful tools for early detection of possible teratogenic properties of compounds that belong to the triazole- group or of compounds with a similar teratogenic mode of action. Chapter 3 provides mechanistic insight into the embryotoxicity of six azoles tested in the rat WEC. Here, we evaluated dose-dependent embryotoxicity of azoles in the rat WEC, calculating the concentration at which the total morphological score (TMS) is 10% decreased (ID10). For the azoles tested we compared the in vitro ID10 for embryotoxicity to the in vivo effective doses, while we also performed a comparative analysis for understanding the toxicological and pharmacological mode of action of azoles in the rat WEC at the level of the transcriptome. Functional analysis of differential gene expression after 4 hours exposure at the ID10 revealed regulation of the sterol biosynthesis pathway and embryonic development genes, dominated by genes in the RA pathway, albeit in a differential way. FLU, ketoconazole and triadimefon were the most potent compounds affecting the RA pathway, while in terms of regulation of sterol function, difenoconazole and ketoconazole showed the most pronounced effects. A similar analysis at the 24-hour time point indicated an additional time-dependent difference in the aforementioned pathways regulated by FLU. Strong in vivo embryotoxic azoles showed also an increased regulation of the RA pathway when tested in vitro. On the other hand, weak or non- embryotoxic azoles showed a non-significant effect on genes that belong in the RA pathway. These observations led us to the conclusion that the toxicological mode of action of azoles was mediated through the RA pathway. In summary, the rat WEC assay in combination with transcriptomics could add mechanistic insight into the embryotoxic potency ranking and functional efficacy of the tested compounds, showing Cyp26a1 and Cyp51 as leader biomarkers of the off- and on- target effects, respectively. Similarly to the previous chapter, in chapter 4, the potency ranking of the majority of the twelve tested azoles obtained based on the TMS in the WEC assay was in agreement with the in vivo potency ranking. Additionally, our expanded transcriptomics data, including gene specific responses of twelve azoles tested at their ID10 in the rat WEC for 4 hours, confirmed the observations of chapter 3 with another set of azoles. Potent embryotoxicants in both in vivo and in vitro assays caused more pronounced effects on the dysregulation of RA- mediated genes. Furthermore, azoles with more pronounced effects on the sterol biosynthesis mediated pathway were tested at a higher concentration, but with the same level of effect (ID10). Due to the increased concentration needed for reaching the same level of morphological effects and the absence of RA-mediated pathway regulation, these azoles were considered as more favourable candidates for clinical and agricultural use. Focusing on monitoring the fungicidal activity of azoles, we also detected an increased sensitivity of the expression of Msmo1, which is an enzyme participating in converting lanosterol for synthesizing cholesterol in the mammalian sterol biosynthesis, together with Cyp51and Nsdhl. This observation led us to the conclusion that Msmo1 could be a better biomarker of effect on the sterol biosynthesis pathway compared to the classical biomarker of this pathway, Cyp51, and this may be of use for further improvement of the assessment of fungicidal activity of azoles or chemicals with similar mode of action. Chapter 5 shows the value of combining toxico-dynamic and -kinetic in vitro approaches for embryotoxicity testing of azoles. We also report on the alterations in gene expression induced by azoles. Both the WEC assay and the embryonic stem cells test (EST) predicted the in vivo potency ranking of the twelve tested azoles with moderate accuracy. Combining these results with relative placental transfer rates (Papp values) as determined in the BeWo cell culture model, increased the predictability of both WEC and EST, with R2 values increasing from 0.51 to 0.87 and from 0.35 to 0.60, respectively. The comparison of these in vitro systems correlated well (R2 = 0.67), correctly identifying the strong and weak embryotoxicants. Evaluating specific gene responses related with the toxicological and fungicidal mode of action of the tested azoles in WEC and EST, we observed that the differential regulation of Dhrs3 and Msmo1 reached higher magnitudes in both systems compared to Cyp26a1 and Cyp51. Establishing sensitive biomarkers across all the in vitro systems for studying the underlying mechanism of action of chemicals, such as azoles, is valuable for comparing alternative in vitro models and for improving insight in the mechanism of developmental toxicity of chemicals. Chapter 6 of this thesis presented the general discussion and future perspectives on different topics raised based on the results obtained in the previously described experimental chapters. The results suggested that the combination of in vitro assays for screening the developmental toxicity of azoles may lead to predictions that are more accurate and in agreement with the in vivo observations. The addition of toxico-kinetics, which the BeWo placental transfer model offered, notably improved the correlations of in vivo and in vitro data. Furthermore, the co- implementation of transcriptomics and the identification of gene biomarkers revealed that despite the tested azoles were classified in the same chemical group, they might have a different mode of toxicological action. In conclusion, future combination of in vitro and in silico alternative approaches appear to be of advantage for screening and prioritizing chemical testing, in the process of assessing the consequences of chemical exposure for human health and the environment.</p

A transcriptomic approach for evaluating the relative potency and mechanism of action of azoles in the rat Whole Embryo Culture.

We evaluated the effect of six azoles on embryonic development in the rat whole embryo culture (WEC). Using the total morphological scoring system (TMS), we calculated the ID10concentration (effective dose for 10% decrease in TMS). For evaluating gene specific responses, we combined previously and newly collected transcriptomics data of rat WEC exposed to a total of twelve azoles at their ID10for 4h. Results revealed shared expressions responses in genes involved in the retinoic acid (RA) and sterol biosynthesis pathways, which are respectively representatives of developmental toxicity and targeted fungicidal action of the azoles. Azoles with more pronounced effects on the regulation of RA-associated genes were generally characterized as more potent embryotoxicants. Overall, compounds with strong sterol biosynthesis related responses and low RA related responses were considered as more favourable candidates, as they specifically regulated genes related to a desired target response. Among the identified sterol associated genes, we detected that methylsterol monooxygenase 1 (Msmo1) was more sensitively induced compared to Cyp51, a classical biomarker of this pathway. Therefore, we suggest that Msmo1 could be a better biomarker for screening the fungicidal value of azoles. In summary, we conclude that the embryonic regulation of RA and sterol metabolic pathways could be indicators for ranking azoles as embryotoxicants and determining their drug efficacy

Flusilazole induces spatio-temporal expression patterns of retinoic acid-, differentiation- and sterol biosynthesis-related genes in the rat Whole Embryo Culture

Embryotoxic responses are critically dependent on the timing of exposure during embryo development. Here, we examined the time- dependent developmental effects in rat embryos exposed to flusilazole (FLU), and their link to retinoic acid (RA) mediated pathways. To this end, we assessed the effects of 4. h exposure of rat embryos in vitro to 300. μM FLU during four developmental time windows (0-4, 4-8, 24-28 and 44-48. h), evaluating morphological parameters, expression and localization of five genes directly or indirectly linked with the RA pathway. These were RA- (Cyp26a1 and Dhrs3), differentiation- (Gbx2 and Cdx1) and sterol biosynthesis- (Cyp51) related genes. Extended exposure for 48. h to 300. μM FLU resulted in morphological changes, typical for triazoles and RA, while the 4. h exposure times did not. Time dependent significant upregulation of the five selected genes was observed. These results corroborate that the embryotoxic responses to FLU are correlated with the regulation of the RA pathway. Thus, these gene expression markers can be considered early biomarkers of FLU-induced potential developmental toxicity later in the development.</p

Impact of Application of Abscisic Acid, Benzothiadiazole and Chitosan on Berry Quality Characteristics and Plant Associated Microbial Communities of Vitis vinifera L var. Mouhtaro Plants

The phenolic profile of the grape berries is a key quality factor for the red grapevine varieties and several techniques have been applied to improve it. An innovative technique is the application of resistance elicitors and phytohormones. In the present study, leaves and berries of a Greek red indigenous variety (Mouhtaro) sprayed with two elicitors, benzothiadiazole and chitosan and a plant hormone abscisic acid, during veraison. Physicochemical and phenolic characteristics of the berries and microbial communities of rhizosphere, phyllosphere and carposphere were analyzed at harvest. Differences in the microbial communities on different plant compartments were observed after the application of the plant activators. Chitosan treatment increased the abundance of the beneficial lactic acid bacteria, while the abscisic acid treatment decreased the presence of spoilage fungi on the carposphere. Treatments differentiate total phenolics, anthocyanins and in the chemical characteristics of grape must with chitosan treated grapes had increased anthocyanins and skin-derived phenolics that correlated positively with the microbial taxa that was discriminant by LefSe analysis. This research provides an overview of the effect of plant activators on the microbial ecology and grape quality of the Greek variety Mouhtaro and presents the potential of new and innovative approaches in the field of sustainable viticulture

Embryotoxic and pharmacologic potency ranking of six azoles in the rat whole embryo culture by morphological and transcriptomic analysis

Differential gene expression analysis in the rat whole embryo culture (WEC) assay provides mechanistic insight into the embryotoxicity of test compounds. In our study, we hypothesized that comparative analysis of the transcriptomes of rat embryos exposed to six azoles (flusilazole, triadimefon, ketoconazole, miconazole, difenoconazole and prothioconazole) could lead to a better mechanism-based understanding of their embryotoxicity and pharmacological action. For evaluating embryotoxicity, we applied the total morphological scoring system (TMS) in embryos exposed for 48 h. The compounds tested showed embryotoxicity in a dose-response fashion. Functional analysis of differential gene expression after 4 h exposure at the ID10 (effective dose for 10% decreased TMS), revealed the sterol biosynthesis pathway and embryonic development genes, dominated by genes in the retinoic acid (RA) pathway, albeit in a differential way. Flusilazole, ketoconazole and triadimefon were the most potent compounds affecting the RA pathway, while in terms of regulation of sterol function, difenoconazole and ketoconazole showed the most pronounced effects. Dose-dependent analysis of the effects of flusilazole revealed that the RA pathway related genes were already differentially expressed at low dose levels while the sterol pathway showed strong regulation at higher embryotoxic doses, suggesting that this pathway is less predictive for the observed embryotoxicity. A similar analysis at the 24-hour time point indicated an additional time-dependent difference in the aforementioned pathways regulated by flusilazole. In summary, the rat WEC assay in combination with transcriptomics could add a mechanistic insight into the embryotoxic potency ranking and pharmacological mode of action of the tested compounds

Corrigendum to ‘Embryotoxic and pharmacologic potency ranking of six azoles in the rat whole embryo culture by morphological and transcriptomic analysis’

The authors regret that in the Results section - 3.7 - the data regarding the inhibition of Cyp26a1 provided by BASF SE were not IC50 values but rather Kd (dissociation constant) values. The authors would like to apologise for any inconvenience caused.</p

Embryotoxic and pharmacologic potency ranking of six azoles in the rat whole embryo culture by morphological and transcriptomic analysis

Differential gene expression analysis in the rat whole embryo culture (WEC) assay provides mechanistic insight into the embryotoxicity of test compounds. In our study, we hypothesized that comparative analysis of the transcriptomes of rat embryos exposed to six azoles (flusilazole, triadimefon, ketoconazole, miconazole, difenoconazole and prothioconazole) could lead to a better mechanism-based understanding of their embryotoxicity and pharmacological action. For evaluating embryotoxicity, we applied the total morphological scoring system (TMS) in embryos exposed for 48h. The compounds tested showed embryotoxicity in a dose-response fashion. Functional analysis of differential gene expression after 4h exposure at the ID10 (effective dose for 10% decreased TMS), revealed the sterol biosynthesis pathway and embryonic development genes, dominated by genes in the retinoic acid (RA) pathway, albeit in a differential way. Flusilazole, ketoconazole and triadimefon were the most potent compounds affecting the RA pathway, while in terms of regulation of sterol function, difenoconazole and ketoconazole showed the most pronounced effects. Dose-dependent analysis of the effects of flusilazole revealed that the RA pathway related genes were already differentially expressed at low dose levels while the sterol pathway showed strong regulation at higher embryotoxic doses, suggesting that this pathway is less predictive for the observed embryotoxicity. A similar analysis at the 24-hour time point indicated an additional time-dependent difference in the aforementioned pathways regulated by flusilazole. In summary, the rat WEC assay in combination with transcriptomics could add a mechanistic insight into the embryotoxic potency ranking and pharmacological mode of action of the tested compounds

A comparison of the embryonic stem cell test and whole embryo culture assay combined with the BeWo placental passage model for predicting the embryotoxicity of azoles

In the present study, we show the value of combining toxico-dynamic and -kinetic in vitro approaches for embryotoxicity testing of azoles. Both the whole embryo culture (WEC) and the embryonic stem cells test (EST) predicted the in vivo potency ranking of twelve tested azoles with moderate accuracy. Combining these results with relative placental transfer rates (Papp values) as determined in the BeWo cell culture model, increased the predictability of both WEC and EST, with R2 values increasing from 0.51 to 0.87 and from 0.35 to 0.60, respectively. The comparison of these in vitro systems correlated well (R2 = 0.67), correctly identifying the in vivo strong and weak embryotoxicants. Evaluating also specific gene responses related with the retinoic acid and sterol biosynthesis pathways, which represent the toxicological and fungicidal mode of action of azoles respectively in the WEC and EST, we observed that the differential regulation of Dhrs3 and Msmo1 reached higher magnitudes in both systems compared to Cyp26a1 and Cyp51. Establishing sensitive biomarkers across the in vitro systems for studying the underlying mechanism of action of chemicals, such as azoles, is valuable for comparing alternative in vitro models and for improving insight in the mechanism of developmental toxicity of chemicals

A comparison of the embryonic stem cell test and whole embryo culture assay combined with the BeWo placental passage model for predicting the embryotoxicity of azoles.

In the present study, we show the value of combining toxico-dynamic and -kinetic in vitro approaches for embryotoxicity testing of azoles. Both the whole embryo culture (WEC) and the embryonic stem cells test (EST) predicted the in vivo potency ranking of twelve tested azoles with moderate accuracy. Combining these results with relative placental transfer rates (Papp values) as determined in the BeWo cell culture model, increased the predictability of both WEC and EST, with R2 values increasing from 0.51 to 0.87 and from 0.35 to 0.60, respectively. The comparison of these in vitro systems correlated well (R2 = 0.67), correctly identifying the in vivo strong and weak embryotoxicants. Evaluating also specific gene responses related with the retinoic acid and sterol biosynthesis pathways, which represent the toxicological and fungicidal mode of action of azoles respectively in the WEC and EST, we observed that the differential regulation of Dhrs3 and Msmo1 reached higher magnitudes in both systems compared to Cyp26a1 and Cyp51. Establishing sensitive biomarkers across the in vitro systems for studying the underlying mechanism of action of chemicals, such as azoles, is valuable for comparing alternative in vitro models and for improving insight in the mechanism of developmental toxicity of chemicals