Flusilazole induced developmental toxicity, neurotoxicity, and cardiovascular toxicity via apoptosis and oxidative stress in zebrafish

Flusilazole is a well-known triazole fungicide applied to various crops and fruits worldwide. Flusilazole residues are frequently detected in the environment, and many researchers have reported the hazardous effects of flusilazole on non-target organisms; however, the developmental toxicity of flusilazole has not been fully elucidated. In this study, we investigated flusilazole-induced developmental defects in zebrafish, which are used in toxicology studies to assess the toxic effects of chemicals on aquatic species or vertebrates. We confirmed that flusilazole exposure affected the viability and hatching rate of zebrafish larvae, and resulted in morphological defects, reduced body length, diminished eye and head sizes, and inflated pericardial edema. Apoptosis, oxidative stress, and inflammation were also observed. These factors interrupted the normal organ formation during early developmental stages, and transgenic models were used to identify organ defects. We confirmed the effects of flusilazole on the nervous system using olig2:dsRed transgenic zebrafish, and on the cardiovascular system using cmlc2:dsRed and fli1:eGFP transgenic zebrafish. Our results demonstrate the developmental toxicity of flusilazole and its mechanisms in zebrafish as well as the detrimental effects of flusilazole.

Network analysis of toxic endpoints of fungicides in zebrafish

Zebrafish being the best animal model to study, every attempt has been made to decipher the toxic mechanism of every fungicide of usage and interest. It is important to understand the multiple targets of a toxicant to estimate the toxic potential in its totality. A total of 22 fungicides of different classes like amisulbrom, azoxystrobin, carbendazim, carboxin, chlorothalonil, difenoconazole, etridiazole, flusilazole, fluxapyroxad, hexaconazole, kresoxim methyl, mancozeb, myclobutanil, prochloraz, propiconazole, propineb, pyraclostrobin, tebuconazole, thiophanate-methyl, thiram, trifloxystrobin and ziram were reviewed and analyzed for their multiple explored targets in zebrafish. Toxic end points in zebrafish are highly informative when it comes to network analysis. They provide a window into the molecular and cellular pathways that are affected by a certain toxin. This can then be used to gain insights into the underlying mechanisms of toxicity and to draw conclusions on the potential of a particular compound to induce toxicity. This knowledge can then be used to inform decisions about drug development, environmental regulation, and other areas of research. In addition, the use of zebrafish toxic end points can also be used to better understand the effects of environmental pollutants on ecosystems. By understanding the pathways affected by a given toxin, researchers can determine how pollutants may interact with the environment and how this could lead to health or environmental impacts.

Transgenerational effects of triazole fungicides on gene expression and egg compounds in non-exposed offspring: A case study using Red-Legged Partridges (Alectoris rufa)

Triazole fungicides are widely used to treat cereal seeds before sowing. Granivorous birds like the Red-legged Partridge (Alectoris rufa) have high exposure risk because they ingest treated seeds that remain on the field surface. As triazole fungicides can act as endocrine disruptors, affecting sterol synthesis and reproduction in birds several months after exposure, we hypothesized that these effects could also impact subsequent generations of exposed birds. To test this hypothesis, we exposed adult partridges (F0) to seeds treated at commercial doses with four different formulations containing triazoles as active ingredients (flutriafol, prothioconazole, tebuconazole, and a mixture of the latter two), simulating field exposure during late autumn sowing. During the subsequent reproductive season, two to four months after exposure, we examined compound allocation of steroid hormones, cholesterol, vitamins, and carotenoids in eggs laid by exposed birds (F1), as well as the expression of genes encoding enzymes involved in sterol biosynthesis in one-day-old chicks of this F1. One year later, F1 animals were paired again to investigate the expression of the same genes in the F2 chicks. We found changes in the expression of some genes for all treatments and both generations. Additionally, we observed an increase in estrone levels in eggs from partridges treated with flutriafol compared to controls, a decrease in tocopherol levels in partridges exposed to the mixture of tebuconazole and prothioconazole, and an increase in retinol levels in partridges exposed to prothioconazole. Despite sample size limitations, this study provides novel insights into the mechanisms of action of the previously observed effects of triazole fungicide-treated seeds on avian reproduction with evidence that the effects can persist beyond the exposure windows, affecting unexposed offspring of partridges fed with treated seeds. The results highlight the importance of considering long-term chronic effects when assessing pesticide risks to wild birds.

A realistic mixture of ubiquitous persistent organic pollutants affects bone and cartilage development in zebrafish by interaction with nuclear receptor signaling

"Persistent organic pollutants (POPs)" have a plethora of deleterious effects on humans and the environment due to their bioaccumulative, persistent, and mimicking properties. Individually, each of these chemicals has been tested and its effects measured, however they are rather found as parts of complex mixtures of which we do not fully grasp the extent of their potential consequences. Here we studied the effects of realistic, environmentally relevant mixtures of 29 POPs on cartilage and bone development using zebrafish as a model species. We observed developmental issues in cartilage, in the form of diverse malformations such as micrognathia, reduced size of the Meckel's and other structures. Also, mineralized bone formation was disrupted, hence impacting the overall development of the larvae at later life stages. Assessment of the transcriptome revealed disruption of nuclear receptor pathways, such as androgen, vitamin D, and retinoic acid, that may explain the mechanisms of action of the compounds within the tested mixtures. In addition, clustering of the compounds using their chemical signatures revealed structural similarities with the model chemicals vitamin D and retinoic acid that can explain the effects and/or enhancing the phenotypes we witnessed. Further mechanistic studies will be required to fully understand this kind of molecular interactions and their repercussions in organisms. Our results contribute to the already existing catalogue of deleterious effects caused by exposure to POPs and help to understand the potential consequences in at risk populations.

Dynamic regulation of gene expression and morphogenesis in the zebrafish embryo test after exposure to all-trans retinoic acid

The zebrafish embryotoxicity test (ZET) is widely used in developmental toxicology. The analysis of gene expression regulation in ZET after chemical exposure provides mechanistic information about the effects of chemicals on morphogenesis in the test. The gene expression response magnitude has been shown to change with exposure duration. The objective of this work is to study the effect of the exposure duration on the magnitude of gene expression changes in the all-trans retinoic acid (ATRA) signaling pathway in the ZET. Retinoic acid regulation is a key driver of morphogenesis and is therefore employed here as an indicator for the regulation of developmental genes. A teratogenic concentration of 7.5 nM of ATRA was given at 3 hrs post fertilization (hpf) for a range of exposure durations until 120 hrs of development. The expression of a selection of genes related to ATRA signaling and downstream developmental genes was determined. The highest magnitudes of gene expression regulation were observed after 2-24 hrs exposure with an optimal response after 4 hrs. Longer exposures showed a decrease in the gene expression response, although continued exposure to 120 hpf caused malformations and lethality. This study shows that assessment of gene expression regulation at early time points after the onset of exposure in the ZET may be optimal for the prediction of developmental toxicity. We believe these results could help optimize sensitivity in future studies with ZET.

Differential expression of genes in C. elegans reveals transcriptional responses to indirect-acting xenobiotic compounds and insensitivity to 2,3,7,8-tetrachlorodibenzodioxin

Caenorhabditis elegans is a well-established model organism for toxicity testing of chemical substances. We recently demonstrated its potential for bioanalysis of the toxic potency of chemical contaminants in water. While many detoxification genes are homologues to those in mammalians, C. elegans is reported to be deficient in cytochrome CYP1-like P450 metabolism and that its aryl hydrocarbon receptor (AhR) homolog encoded by ahr-1 purportedly does not interact with dioxins or any other known xenobiotic ligand. This suggests that C. elegans is insensitive for compounds that require bioactivation (indirectly acting compounds) and for dioxins or dioxin-like compounds. This study analysed genome-wide gene expression of the nematode in response to 30 μM of aflatoxin B1 (AFB1), benzo(a)pyrene (B(a)P), Aroclor 1254 (PCB1254), and 10 μM of 2,3,7,8-tetrachlorodibenzodioxin (TCDD). After 24 h of exposure in the early L4 larval stage, microarray analysis revealed 182, 86, and 321 differentially expressed genes in the nematodes treated with 30 μM of AFB1, B(a)P, and PCB1254, respectively. Among these genes, many encode xenobiotic-metabolizing enzymes, and their transcription levels were among the highest-ranked fold-changed genes. Interestingly, only one gene (F59B1.8) was upregulated in the nematodes exposed to 10 μM TCDD. Genes related to metabolic processes and catalytic activity were the most induced by exposure to 30 μM of AFB1, B(a)P, and PCB1254. Despite the genotoxic nature of AFB1 and B(a)P, no differential expression was found in the genes encoding DNA repair and cell cycle checkpoint proteins. Analysis of concentration-response curves was performed to determine the Lowest Observed Transcriptomic Effect Levels (LOTEL) of AFB1, B(a)P, and PCB1254. The obtained LOTEL values showed that gene expression changes in C. elegans are more sensitive to toxicants than reproductive effects. Overall, transcriptional responses of metabolic enzymes suggest that the nematode does metabolize AFB1, B(a)P, and PCB1254. Our findings also support the assumption that the transcription factor AhR homolog in C. elegans does not bind typical xenobiotic ligands, rendering the nematode transcriptionally insensitive to TCDD effects.

Effect of fluensulfone on different functional genes of root-knot nematode Meloidogyne incognita

Meloidogyne incognita is an obligate plant-parasitic nematode causing serious damage to agricultural crops. Major constraints in nematode management arose due to the limited availability of non-fumigant nematicides in conjunction with the considerable ill effects of fumigants on human and non-target organisms. Recently, fluensulfone has been reported to be an effective non-fumigant nematicide against plant-parasitic nematodes and the model nematode Caenorhabditis elegans. The nematicidal efficacy varies according to its concentration at the time of application, exposure timing, nematode species variability, and even across subpopulations within the same species. It interferes with the key physiological processes of nematodes, like motility, behavior, chemosensation, stylet thrusting, infectivity, metabolism, lipid consumption, tissue integrity, oviposition, egg hatching, and survival. However, the molecular basis of these multivariate physiological anomalies is still largely unknown. Quantitative real-time PCR was carried out to understand the acute transcriptional perturbation of 30 functional genes associated with key physiological and life processes in a M. incognita population, following exposure of 10, 50, and 100 ppm of fluensulfone for 5 and 10 hr. The chemical treatment resulted in significant downregulation of all the neuropeptidergic genes, with concomitant repression of majority of genes related to chemosensation, esophageal gland secretion, parasitism, fatty acid metabolism, and G-protein coupled receptors. Collectively, the parasitism genes were found to be perturbed at highest magnitude, followed by the GPCRs and neuropeptidergic genes. These results establish the wide ranging effect of fluensulfone on various metabolic and physiological pathways of nematode.

Preclinical Toxicity and Safety of MM-129-First-in-Class BTK/PD-L1 Inhibitor as a Potential Candidate against Colon Cancer

MM-129 is a novel inhibitor targeting BTK/PI3K/AKT/mTOR and PD-L1, as it possesses antitumor activity against colon cancer. To evaluate the safety profile of MM-129, we conducted a toxicity study using the zebrafish and rodent model. MM-129 was also assessed for pharmacokinetics features through an in vivo study on Wistar rats. The results revealed that MM-129 exhibited favorable pharmacokinetics with quick absorption and 68.6% of bioavailability after intraperitoneal administration. No serious adverse events were reported for the use of MM-129, confirming a favorable safety profile for this compound. It was not fatal and toxic to mice at an anticancer effective dose of 10 μmol/kg. At the end of 14 days of administering hematological and biochemical parameters, liver and renal functions were all at normal levels. No sublethal effects were either detected in zebrafish embryos treated with a concentration of 10 μM. MM-129 has the potential as a safe and well-tolerated anticancer formulation for future treatment of patients with colon cancer.

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.

Role of oxidative stress in clofazimine-induced cardiac dysfunction in a zebrafish model

BACKGROUND

Clofazimine (CFZ), a riminophenazine, is now commonly used in the treatment of multidrug-resistant tuberculosis. However, its use may be potentially associated with cardiac dysfunction in some individuals. In this study, the zebrafish heart, by merit of its developmental and genetic characteristics being in homology with that of human, was chosen as an animal model for evaluation of such dysfunction.

METHODS

Morphological and physiological parameters were used to assess cardiac dysfunction. Transcriptome analysis was performed, followed by validation with real-time quantitative PCR, for delineation of the relevant genomics.

RESULTS

Exposure of 2 dpf zebrafish to 4 mg/L CFZ for 2 days, adversely affected cardiac functions including significant decreases in HR, SV, CO, and FS, with observable pathophysiological developments of pericardial effusion and blood accumulation in the heart, in comparison with the control group. In addition, genes which respond to xenobiotic stimulus, related to oxygen transport, glutathione metabolism and extracellular matrix -receptor interactions, were significantly enriched among the differentially up-regulated genes. Antioxidant response element motif was enriched in the 5000 base pair upstream regions of the differentially expressed genes. Co-administration of N-acetylcysteine was shown to protect zebrafish against the development of CFZ-induced cardiac dysfunction.

CONCLUSIONS

This study suggests an important role of oxidative stress as a major pathogenetic mechanism of riminophenazine-induced cardiac dysfunction.

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.

Molecular fingerprints of conazoles via functional genomic profiling of Saccharomyces cerevisiae

Conazoles were designed to inhibit ergosterol biosynthesis. Conazoles have been widely used as agricultural fungicides and are frequently detected in the environment. Although conazoles have been reported to have adverse effects, such as potential carcinogenic effects, the underlying molecular mechanisms of toxicity remain unclear. Here, the molecular fingerprints of five conazoles (propiconazole (Pro), penconazole (Pen), tebuconazole (Teb), flusilazole (Flu) and epoxiconazole (Epo)) were assessed in Saccharomyces cerevisiae (yeast) via functional genome-wide knockout mutant profiling. A total of 169 (4.49%), 176 (4.67%), 198 (5.26%), 218 (5.79%) and 173 (4.59%) responsive genes were identified at three concentrations (IC₅₀, IC₂₀ and IC₁₀) of Pro, Pen, Teb, Flu and Epo, respectively. The five conazoles tended to have similar gene mutant fingerprints and toxicity mechanisms. "Ribosome" (sce03010) and "cytoplasmic translation" (GO: 0002181) were the common KEGG pathway and GO biological process term by gene set enrichment analysis of the responsive genes, which suggested that conazoles influenced protein synthesis. Conazoles also affected fatty acids synthesis because "biosynthesis of unsaturated fatty acids" pathway was among the top-ranked KEGG pathways. Moreover, two genes, YGR037C (acyl-CoA-binding protein) and YCR034W (fatty acid elongase), were key fingerprints of conazoles because they played vital roles in conazole-induced toxicity. Overall, the fingerprints derived from the yeast functional genomic screening provide an alternative approach to elucidate the molecular mechanisms of environmental pollutant conazoles.

Underlying mechanisms of reproductive toxicity caused by multigenerational exposure of 2, bromo-4, 6-dinitroaniline (BDNA) to Zebrafish (Danio rerio) at environmental relevant levels

2-bromo-4, 6-dinitroaniline (BDNA) is a mutagenic aromatic amine involved in the production and degradation of Disperse blue 79, one of the most extensively used brominated azo dyes. In our previous study, a multigenerational exposure of BDNA (0.5, 5, 50 and 500 μg/L) to zebrafish from F₀ adult to F₂ larvae including a recovery group in F₂ larvae was conducted. The effects on apical points observed in individuals and the long-term effects predicted on population were all related to reproduction. In this study, we performed molecular analysis to elucidate the underlying mechanisms of the reproductive toxicity of BDNA. In F₁ generation, measurement of vitellogenin and transcription levels of genes associated with hypothalamus-pituitary-gland (HPG) axis, estrogen receptor (ER) and androgen receptor (AR) were conducted. There was a decrease in VTG level in the blood of F₁ female fish and transcription of genes related to ER was more affected than that of genes related to AR. These results were consistent with adverse effects that sexual differentiation was biased towards males and fecundity was impaired in a concentration-dependent manner in adults of F₁ generation after 150 days exposure. In F₂ generation, global gene transcriptions of F₂ larvae were investigated. It was uncovered that processes related to apoptosis, development and DNA damage were strongly affected. Alterations to these biological pathways accounted for the irreversible parental influence on a significant decrease in hatchability and increase in abnormality of F₂ larvae. All evidence suggested that the multigenerational exposure of BDNA posed lasting effects transmitted from parents to offspring that persisted after exposure ceased.

Map and model-moving from observation to prediction in toxicogenomics

BACKGROUND

Chemicals induce compound-specific changes in the transcriptome of an organism (toxicogenomic fingerprints). This provides potential insights about the cellular or physiological responses to chemical exposure and adverse effects, which is needed in assessment of chemical-related hazards or environmental health. In this regard, comparison or connection of different experiments becomes important when interpreting toxicogenomic experiments. Owing to lack of capturing response dynamics, comparability is often limited. In this study, we aim to overcome these constraints.

RESULTS

We developed an experimental design and bioinformatic analysis strategy to infer time- and concentration-resolved toxicogenomic fingerprints. We projected the fingerprints to a universal coordinate system (toxicogenomic universe) based on a self-organizing map of toxicogenomic data retrieved from public databases. Genes clustering together in regions of the map indicate functional relation due to co-expression under chemical exposure. To allow for quantitative description and extrapolation of the gene expression responses we developed a time- and concentration-dependent regression model. We applied the analysis strategy in a microarray case study exposing zebrafish embryos to 3 selected model compounds including 2 cyclooxygenase inhibitors. After identification of key responses in the transcriptome we could compare and characterize their association to developmental, toxicokinetic, and toxicodynamic processes using the parameter estimates for affected gene clusters. Furthermore, we discuss an association of toxicogenomic effects with measured internal concentrations.

CONCLUSIONS

The design and analysis pipeline described here could serve as a blueprint for creating comparable toxicogenomic fingerprints of chemicals. It integrates, aggregates, and models time- and concentration-resolved toxicogenomic data.

Is fish embryo test (FET) according to OECD 236 sensible enough for delivering quality data for effluent risk assessment?

Over the past few years, the fish embryo test (FET) has become widely accepted as an animal-friendly protocol for ecotoxicological research. As Organisation for Economic Co-operation and Development (OECD) 236, the FET has been widely applied for simple mixture exposures under the Registration, Evaluation, Authorisation, and Restriction of Chemicals regulation of the European Union; and now its use is spreading worldwide as a supposedly reliable whole-effluent test (i.e., the testing of complex mixture exposures). However, comparative peer-reviewed data regarding the FET's efficiency for whole-effluent tests are virtually nonexistent. The primary objective of the present study was to make the first comparative test between the FET according to OECD 236 and other standard and slightly modified standard fish protocols used worldwide for whole-effluent tests. For that, we used an untreated hospital effluent considered to be highly toxic but disposed of in municipal sewerage. The base methods were OECD 203 (juvenile), US Environmental Protection Agency method 2000.0 (larvae), and OECD 236 (embryo). We also evaluated the addition of 3 virtually costless sublethal metrics (immobility, nonhatching, and pericardial edema) that could enhance the sensitivity of OECD 236. We observed acute toxicity in all 8 methodologies tested, with a clear escalation in sensitivity (larvae > juvenile ≥ embryo). Larvae were the most sensitive life stage for whole-effluent tests. The addition of sublethal metrics to OECD 236 enhanced its previous sensitivity in over 30%. Thus we conclude that OECD 236 acts below its potential and that the embryonic stage (as used in the FET) may not be the most sensitive life stage for whole-effluent tests. Environ Toxicol Chem 2018;37:2925-2932. © 2018 SETAC.

A Qualitative Modeling Approach for Whole Genome Prediction Using High-Throughput Toxicogenomics Data and Pathway-Based Validation

Efficient high-throughput transcriptomics (HTT) tools promise inexpensive, rapid assessment of possible biological consequences of human and environmental exposures to tens of thousands of chemicals in commerce. HTT systems have used relatively small sets of gene expression measurements coupled with mathematical prediction methods to estimate genome-wide gene expression and are often trained and validated using pharmaceutical compounds. It is unclear whether these training sets are suitable for general toxicity testing applications and the more diverse chemical space represented by commercial chemicals and environmental contaminants. In this work, we built predictive computational models that inferred whole genome transcriptional profiles from a smaller sample of surrogate genes. The model was trained and validated using a large scale toxicogenomics database with gene expression data from exposure to heterogeneous chemicals from a wide range of classes (the Open TG-GATEs data base). The method of predictor selection was designed to allow high fidelity gene prediction from any pre-existing gene expression data set, regardless of animal species or data measurement platform. Predictive qualitative models were developed with this TG-GATES data that contained gene expression data of human primary hepatocytes with over 941 samples covering 158 compounds. A sequential forward search-based greedy algorithm, combining different fitting approaches and machine learning techniques, was used to find an optimal set of surrogate genes that predicted differential expression changes of the remaining genome. We then used pathway enrichment of up-regulated and down-regulated genes to assess the ability of a limited gene set to determine relevant patterns of tissue response. In addition, we compared prediction performance using the surrogate genes found from our greedy algorithm (referred to as the SV2000) with the landmark genes provided by existing technologies such as L1000 (Genometry) and S1500 (Tox21), finding better predictive performance for the SV2000. The ability of these predictive algorithms to predict pathway level responses is a positive step toward incorporating mode of action (MOA) analysis into the high throughput prioritization and testing of the large number of chemicals in need of safety evaluation.

The Transcriptome of the Zebrafish Embryo After Chemical Exposure: A Meta-Analysis

Numerous studies have been published in the past years investigating the transcriptome of the zebrafish embryo (ZFE) upon being subjected to chemical stress. Aiming at a more mechanistic understanding of the results of such studies, knowledge about commonalities of transcript regulation in response to chemical stress is needed. Thus, our goal in this study was to identify and interpret genes and gene sets constituting a general response to chemical exposure. Therefore, we aggregated and reanalyzed published toxicogenomics data obtained with the ZFE. We found that overlap of differentially transcribed genes in response to chemical stress across independent studies is generally low and the most commonly differentially transcribed genes appear in less than 50% of all treatments across studies. However, effect size analysis revealed several genes showing a common trend of differential expression, among which genes related to calcium homeostasis emerged as key, especially in exposure settings up to 24 h post-fertilization. Additionally, we found that these and other downregulated genes are often linked to anatomical regions developing during the respective exposure period. Genes showing a trend of increased expression were, among others, linked to signaling pathways (e.g., Wnt, Fgf) as well as lysosomal structures and apoptosis. The findings of this study help to increase the understanding of chemical stress responses in the developing zebrafish embryo and provide a starting point to improve experimental designs for this model system. In future, improved time- and concentration-resolved experiments should offer better understanding of stress response patterns and access to mechanistic information.

Editor's Highlight: Transgenic Zebrafish Reporter Lines as Alternative In Vivo Organ Toxicity Models

Organ toxicity, particularly liver toxicity, remains one of the major reasons for the termination of drug candidates in the development pipeline as well as withdrawal or restrictions of marketed drugs. A screening-amenable alternative in vivo model such as zebrafish would, therefore, find immediate application in the early prediction of unacceptable organ toxicity. To identify highly upregulated genes as biomarkers of toxic responses in the zebrafish model, a set of well-characterized reference drugs that cause drug-induced liver injury (DILI) in the clinic were applied to zebrafish larvae and adults. Transcriptome microarray analysis was performed on whole larvae or dissected adult livers. Integration of data sets from different drug treatments at different stages identified common upregulated detoxification pathways. Within these were candidate biomarkers which recurred in multiple treatments. We prioritized 4 highly upregulated genes encoding enzymes acting in distinct phases of the drug metabolism pathway. Through promoter isolation and fosmid recombineering, eGFP reporter transgenic zebrafish lines were generated and evaluated for their response to DILI drugs. Three of the 4 generated reporter lines showed a dose and time-dependent induction in endodermal organs to reference drugs and an expanded drug set. In conclusion, through integrated transcriptomics and transgenic approaches, we have developed parallel independent zebrafish in vivo screening platforms able to predict organ toxicities of preclinical drugs.

Retinoic acid in developmental toxicology: Teratogen, morphogen and biomarker

This review explores the usefulness retinoic acid (RA) related physiological factors as possible biomarkers of embryotoxicity. RA is involved in the morphogenesis of the early embryo as well as in the development and maturation of a wide variety of organ anlagen. The region-specific homeostasis of RA in the embryo is in many ways the driving force determining developmental cell proliferation versus differentiation. As a consequence, RA concentrations are carefully controlled in time and space in the developing embryo. RA deficiency and overdosing both result in characteristic patterns of malformations that may involve many different organ systems. The central role of RA in embryo development provides us with a set of sensitive biomarkers that may be employed in developmental toxicity testing. This includes the synthesizing and metabolizing enzymes of RA, but also a myriad of related morphogenetic factors and their genes, of which the expression may be affected by changes in RA balance. Several examples of embryotoxicants interfering with the homeostasis of RA and related parameters have been described. A preliminary adverse outcome pathway framework for RA mediated malformations has been published. Expansion of this framework and its application in developmental toxicity testing may allow the detection of a large variety of embryotoxicants with diverse modes of action. RA homeostasis therefore provides a promising set of molecular tools that may be employed in the advancement of mode of action driven animal-free developmental toxicity testing.

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.

Gene Expression Profiling in Fish Toxicology: A Review

In this review, we present an overview of transcriptomic responses to chemical exposures in a variety of fish species. We have discussed the use of several molecular approaches such as northern blotting, differential display reverse transcription-polymerase chain reaction (DDRT-PCR), suppression subtractive hybridization (SSH), real time quantitative PCR (RT-qPCR), microarrays, and next-generation sequencing (NGS) for measuring gene expression. These techniques have been mainly used to measure the toxic effects of single compounds or simple mixtures in laboratory conditions. In addition, only few studies have been conducted to examine the biological significance of differentially expressed gene sets following chemical exposure. Therefore, future studies should focus more under field conditions using a multidisciplinary approach (genomics, proteomics and metabolomics) to understand the synergetic effects of multiple environmental stressors and to determine the functional significance of differentially expressed genes. Nevertheless, recent developments in NGS technologies and decreasing costs of sequencing holds the promise to uncover the complexity of anthropogenic impacts and biological effects in wild fish populations.

Using Transcriptomics to Evaluate Thresholds in Genotoxicity Dose–Response

Several genotoxic chemicals have been reported to produce threshold-shaped dose–response curves for mutation and genotoxicity assays, both in vivo and in vitro, challenging the current default practice for risk assessment of genotoxic chemicals, which assumes a linear dose–response below the lowest tested dose. Statistical methods cannot determine whether a biological threshold exists with sufficient confidence to overturn this assumption of linearity. Indeed, to truly define the shape of the dose–response curves, we must look to the underlying biology and develop targeted experiments to identify and measure the key processes governing the response of the cell to DNA damage. This chapter describes a series of studies aimed at defining the key transcriptional responses. Two approaches were taken to evaluate transcriptional responses preventing micronucleus induction: (1) comparison of gene signatures for several prototype compounds at a single chemical dose that led to a similar activation of the p53-DNA damage pathway (i.e. 1.5-fold increase in total p53); and (2) evaluation of a subset of chemicals with in-depth dose–response studies. The goal of these efforts was to determine the transcriptional pathways responsible for maintaining homeostasis at low levels of DNA damage, i.e., the biological underpinning of threshold-shaped dose–response curves for mutagenicity.

Fluoride caused thyroid endocrine disruption in male zebrafish (Danio rerio)

Excessive fluoride in natural water ecosystem has the potential to detrimentally affect thyroid endocrine system, but little is known of such effects or underlying mechanisms in fish. In the present study, we evaluated the effects of fluoride on growth performance, thyroid histopathology, thyroid hormone levels, and gene expressions in the HPT axis in male zebrafish (Danio rerio) exposed to different determined concentrations of 0.1, 0.9, 2.0 and 4.1 M of fluoride to investigate the effects of fluoride on thyroid endocrine system and the potential toxic mechanisms caused by fluoride. The results indicated that the growth of the male zebrafish used in the experiments was significantly inhibited, the thyroid microtrastructure was changed, and the levels of T3 and T4 were disturbed in fluoride-exposed male fish. In addition, the expressional profiles of genes in HPT axis displayed alteration. The expressions of all studied genes were significantly increased in all fluoride-exposed male fish after exposure for 45 days. The transcriptional levels of corticotrophin-releasing hormone (CRH), thyroid-stimulating hormone (TSH), thyroglobulin (TG), sodium iodide symporter (NIS), iodothyronine I (DIO1), and thyroid hormone receptor alpha (TRα) were also elevated in all fluoride-exposed male fish after 90 days of exposure, while the inconsistent expressions were found in the mRNA of iodothyronineⅡ (DIO2), UDP glucuronosyltransferase 1 family a, b (UGT1ab), transthyretin (TTR), and thyroid hormone receptor beta (TRβ). These results demonstrated that fluoride could notably inhibit the growth of zebrafish, and significantly affect thyroid endocrine system by changing the microtrastructure of thyroid, altering thyroid hormone levels and endocrine-related gene expressions in male zebrafish. All above indicated that fluoride could pose a great threat to thyroid endocrine system, thus detrimentally affected the normal function of thyroid of male zebrafish.

High-throughput concentration-response analysis for omics datasets

Omics-based methods are increasingly used in current ecotoxicology. Therefore, a large number of observations for various toxic substances and organisms are available and may be used for identifying modes of action, adverse outcome pathways, or novel biomarkers. For these purposes, good statistical analysis of toxicogenomic data is vital. In contrast to established ecotoxicological techniques, concentration-response modeling is rarely used for large datasets. Instead, statistical hypothesis testing is prevalent, which provides only a limited scope for inference. The present study therefore applied automated concentration-response modeling for 3 different ecotoxicotranscriptomic and ecotoxicometabolomic datasets. The modeling process was performed by simultaneously applying 9 different regression models, representing distinct mechanistic, toxicological, and statistical ideas that result in different curve shapes. The best-fitting models were selected by using Akaike's information criterion. The linear and exponential models represented the best data description for more than 50% of responses. Models generating U-shaped curves were frequently selected for transcriptomic signals (30%), and sigmoid models were identified as best fit for many metabolomic signals (21%). Thus, selecting the models from an array of different types seems appropriate, because concentration-response functions may vary because of the observed response type, and they also depend on the compound, the organism, and the investigated concentration and exposure duration range. The application of concentration-response models can help to further tap the potential of omics data and is a necessary step for quantitative mixture effect assessment at the molecular response level.

The fish embryo toxicity test as a replacement for the larval growth and survival test: A comparison of test sensitivity and identification of alternative endpoints in zebrafish and fathead minnows

The fish embryo toxicity (FET) test has been proposed as an alternative to the larval growth and survival (LGS) test. The objectives of the present study were to evaluate the sensitivity of the FET and LGS tests in fathead minnows (Pimephales promelas) and zebrafish (Danio rerio) and to determine if the inclusion of sublethal metrics as test endpoints could enhance test utility. In both species, LGS and FET tests were conducted using 2 simulated effluents. A comparison of median lethal concentrations determined via each test revealed significant differences between test types; however, it could not be determined which test was the least and/or most sensitive. At the conclusion of each test, developmental abnormalities and the expression of genes related to growth and toxicity were evaluated. Fathead minnows and zebrafish exposed to mock municipal wastewater-treatment plant effluent in a FET test experienced an increased incidence of pericardial edema and significant alterations in the expression of genes including insulin-like growth factors 1 and 2, heat shock protein 70, and cytochrome P4501A, suggesting that the inclusion of these endpoints could enhance test utility. The results not only show the utility of the fathead minnow FET test as a replacement for the LGS test but also provide evidence that inclusion of additional endpoints could improve the predictive power of the FET test.

Antibiotic toxicity and absorption in zebrafish using liquid chromatography-tandem mass spectrometry

Evaluation of drug toxicity is necessary for drug safety, but in vivo drug absorption is varied; therefore, a rapid, sensitive and reliable method for measuring drugs is needed. Zebrafish are acceptable drug toxicity screening models; we used these animals with a liquid chromatography-tandem mass spectrometry (LC-MS/MS) method in a multiple reaction monitoring mode to quantify drug uptake in zebrafish to better estimate drug toxicity. Analytes were recovered from zebrafish homogenate by collecting supernatant. Measurements were confirmed for drugs in the range of 10-1,000 ng/mL. Four antibiotics with different polarities were tested to explore any correlation of drug polarity, absorption, and toxicity. Zebrafish at 3 days post-fertilization (dpf) absorbed more drug than those at 6 h post-fertilization (hpf), and different developmental periods appeared to be differentially sensitive to the same compound. By observing abnormal embryos and LD50 values, zebrafish embryos at 6 hpf were considered to be suitable for evaluating embryotoxicity. Also, larvae at 3 dpf were adapted to measure acute drug toxicity in adult mammals. Thus, we can exploit zebrafish to study drug toxicity and can reliably quantify drug uptake with LC-MS/MS. This approach will be helpful for future studies of toxicology in zebrafish.

Egg-in-cube: design and fabrication of a novel artificial eggshell with functionalized surface

An eggshell is a porous microstructure that regulates the passage of gases to allow respiration. The chick embryo and its circulatory system enclosed by the eggshell has become an important model for biomedical research such as the control of angiogenesis, cancer therapy, and drug delivery test, because the use of embryo is ethically acceptable and it is inexpensive and small. However, chick embryo and extra-embryonic blood vessels cannot be accessed freely and has poor observability because the eggshell is tough and cannot be seen through, which limits its application. In this study, a novel artificial eggshell with functionalized surface is proposed, which allows the total amount of oxygen to pass into the egg for the chick embryo culturing and has high observability and accessibility for embryo manipulation. First, a 40-mm enclosed cubic-shaped eggshell consisting of a membrane structure and a rigid frame structure is designed, and then the threshold of the membrane thickness suitable for the embryo survival is figured out according to the oxygen-permeability of the membrane structure. The designed artificial eggshell was actually fabricated by using polydimethylsiloxane (PDMS) and polycarbonate (PC) in the current study. Using the fabricated eggshell, chick embryo and extra-embryonic blood vessels can be observed from multiple directions. To test the effectiveness of the design, the cubic eggshells were used to culture chick embryos and survivability was confirmed when PDMS membranes with adequate oxygen permeability were used. Since the surface of the eggshell is transparent, chick embryo tissue development could be observed during the culture period. Additionally, the chick embryo tissues could be accessed and manipulated from outside the cubic eggshell, by using mechanical tools without breakage of the eggshell. The proposed "Egg-in-Cube" with functionalized surface has great potential to serve as a promising platform for biomedical research.

Toxicogenomic analysis in the combined effect of tributyltin and benzo[a]pyrene on the development of zebrafish embryos

There is a growing recognition that the toxic effects of chemical mixtures are been an important issue in toxicological sciences. Tributyltin (TBT) and benzo[a]pyrene (BaP) are widespread pollutants that occur simultaneously in the aquatic environments. This study was designed to examine comprehensively the combined effects of TBT and BaP on zebrafish (Danio rerio) embryos using toxicogenomic approach combined with biochemical detection and morphological analysis, and tried to gain insight into the mechanisms underlying the combined effects of TBT and BaP. The results of toxicogenomic data indicated that: (1) TBT cotreatment rescued the embryos from decreased hatching ratio caused by BaP alone, while the alteration of gene expression (in this article the phrase gene expression is used as a synonym to gene transcription, although in is acknowledged that gene expression can also be regulated by, e.g., translation and mRNA or protein stability) relative to zebrafish hatching in the BaP groups was resumed by the cotreatment with TBT; (2) BaP cotreatment decreased TBT-mediated dorsal curvature, and alleviated the perturbation of Notch pathway caused by TBT alone; (3) cotreatment with TBT decreased BaP-mediated bradycardia, which might be due to that TBT cotreatment alleviated the perturbation in expression of genes related to cardiac muscle cell development and calcium handling caused by BaP alone; 4) TBT cotreatment brought an antagonistic effect on the BaP-mediated oxidative stress and DNA damage. These results suggested that toxicogenomic approach was available for analyzing combined toxicity with high sensitivity and accuracy, which might improve our understanding and predictability for the combined effects of chemicals.

An adverse outcome pathway framework for neural tube and axial defects mediated by modulation of retinoic acid homeostasis

Developmental toxicity can be caused through a multitude of mechanisms and can therefore not be captured through a single simple mechanistic paradigm. However, it may be possible to define a selected group of overarching mechanisms that might allow detection of the vast majority of developmental toxicants. Against this background, we have explored the usefulness of retinoic acid mediated regulation of neural tube and axial patterning as a general mechanism that, when perturbed, may result in manifestations of developmental toxicity that may cover a large part of malformations known to occur in experimental animals and in man. Through a literature survey, we have identified key genes in the regulation of retinoic acid homeostasis, as well as marker genes of neural tube and axial patterning, that may be used to detect developmental toxicants in in vitro systems. A retinoic acid-neural tube/axial patterning adverse outcome pathway (RA-NTA AOP) framework was designed. The framework was tested against existing data of flusilazole exposure in the rat whole embryo culture, the zebrafish embryotoxicity test, and the embryonic stem cell test. Flusilazole is known to interact with retinoic acid homeostasis, and induced common and unique NTA marker gene changes in the three test systems. Flusilazole-induced changes were similar in directionality to gene expression responses after retinoic acid exposure. It is suggested that the RA-NTA framework may provide a general tool to define mechanistic pathways and biomarkers of developmental toxicity that may be used in alternative in vitro assays for the detection of embryotoxic compounds.

Use of toxicogenomics to predict the potential toxic effect of Benzo(a)pyrene on zebrafish embryos: ocular developmental toxicity

Benzo(a)pyrene (BaP) is a representative polycyclic aromatic hydrocarbon (PAH), which is ubiquitous in the environment. The toxic effects of BaP on fish embryos have been described in detail, but some potentially toxic effects of BaP might have been neglected owing to the limitations of traditional techniques. In the present research, global transcriptional patterns were used to study the potentially toxic effects of BaP, as well as its underlying toxicological mechanisms. The expression levels of multiple genes were significantly changed by BaP exposure. The results of ontology assignments and cluster analysis showed that BaP could affect the processes of photoreceptor maintenance and phototransduction. We also conducted an experiment on phototactic response and found that larvae exposed to BaP displayed a decreasing response to light. In addition, BaP exposure decreased the cellular density of the ganglion cell layer (GCL) significantly. These results suggested that BaP exposure induced visual system developmental defects and dysfunction by perturbation of photoreceptor development related genes. Our results were helpful for an understanding of the toxicity of BaP. This study also indicated that microarray analysis was effective for predicting the potential toxicity of chemicals with high sensitivity and accuracy.

Teratological effects of a panel of sixty water-soluble toxicants on zebrafish development

The zebrafish larva is a promising whole-animal model for safety pharmacology, environmental risk assessment, and developmental toxicity. This model has been used for the high-throughput toxicity screening of various compounds. Our aim here is to identify possible phenotypic markers of teratogenicity in zebrafish embryos that could be used for the assaying compounds for reproductive toxicity. We have screened a panel of 60 water-soluble toxicants to examine their effects on zebrafish development. A total of 22,080 wild-type zebrafish larvae were raised in 250 μL defined buffer in 96-well plates at a plating density of one embryo per well. They were exposed for a 96-h period starting at 24 h post-fertilization. A logarithmic concentration series was used for range-finding, followed by a narrower geometric series for developmental toxicity assessment. A total of 9017 survivors were analyzed at 5 days post-fertilization for nine phenotypes, namely, (1) normal, (2) pericardial oedema, (3) yolk sac oedema, (4) melanophores dispersed, (5) bent tail tip, (6) bent body axis, (7) abnormal Meckel's cartilage, (8) abnormal branchial arches, and (9) uninflated swim bladder. For each toxicant, the EC50 (concentration required to produce one or more of these abnormalities in 50% of embryos) was also calculated. For the majority of toxicants (55/60) there was, at the population level, a statistically significant, concentration-dependent increase in the incidence of abnormal phenotypes among survivors. The commonest abnormalities were pericardial oedema, yolk sac oedema, dispersed melanophores, and uninflated swim bladder. It is possible therefore that these could prove to be general indicators of reproductive toxicity in the zebrafish embryo assay.

Current perspectives on the use of alternative species in human health and ecological hazard assessments

BACKGROUND

Traditional animal toxicity tests can be time and resource intensive, thereby limiting the number of chemicals that can be comprehensively tested for potential hazards to humans and/or to the environment.

OBJECTIVE

We compared several types of data to demonstrate how alternative models can be used to inform both human and ecological risk assessment.

METHODS

We reviewed and compared data derived from high throughput in vitro assays to fish reproductive tests for seven chemicals. We investigated whether human-focused assays can be predictive of chemical hazards in the environment. We examined how conserved pathways enable the use of nonmammalian models, such as fathead minnow, zebrafish, and Xenopus laevis, to understand modes of action and to screen for chemical risks to humans.

RESULTS

We examined how dose-dependent responses of zebrafish embryos exposed to flusilazole can be extrapolated, using pathway point of departure data and reverse toxicokinetics, to obtain human oral dose hazard values that are similar to published mammalian chronic toxicity values for the chemical. We also examined how development/safety data for human health can be used to help assess potential risks of pharmaceuticals to nontarget species in the environment.

DISCUSSION

Using several examples, we demonstrate that pathway-based analysis of chemical effects provides new opportunities to use alternative models (nonmammalian species, in vitro tests) to support decision making while reducing animal use and associated costs.

CONCLUSIONS

These analyses and examples demonstrate how alternative models can be used to reduce cost and animal use while being protective of both human and ecological health.

Gene expression responses for detecting sublethal effects of xenobiotics and whole effluents on a Xenopus laevis embryo assay

In the present study, the authors investigated the effects of bisphenol A, chlorpyrifos, methylparaben, and 2 effluent samples from wastewater treatment plants located in the province of Madrid, Spain, on the messenger RNA expression of specific genes involved in early development (ESR1, pax6, bmp4, and myf5) and a gene involved in the general stress response (hsp70) during Xenopus laevis embryo development. Gene expression was analyzed after 4 h, 24 h, and 96 h of exposure by semiquantitative reverse-transcriptase-polymerase chain reaction. Concentration ranges of the compounds and dilutions for the samples were selected to cause morphological alterations in embryos after 96 h of exposure. Transcript levels of ESR1, pax6, and hsp70 were differentially altered at early developmental stages with patterns specific to the contaminant and the exposure time. However, further studies are needed to establish transcript levels of specific genes as biomarkers of sublethal effects in an environmental risk-assessment framework. Besides, studies including more generic responses, such as genes encoding antioxidant enzymes, together with genes related to embryonic development have to be developed to look for a battery of mechanistic endpoints for the evaluation of chemical exposure at the molecular level in a first-tier assessment.

Thyroid endocrine disruption in zebrafish larvae following exposure to hexaconazole and tebuconazole

The widely used triazole fungicides have the potential to disrupt endocrine system, but little is known of such effects or underlying mechanisms of hexaconazole (HEX) and tebuconazole (TEB) in fish. In the present study, zebrafish (Danio rerio) embryos were exposed to various concentrations of HEX (0.625, 1.25 and 2.5 mg/L) and TEB (1, 2 and 4 mg/L) from fertilization to 120 h post-fertilization (hpf). The whole body content of thyroid hormone and transcription of genes in the hypothalamic-pituitary-thyroid (HPT) axis were analyzed. The results showed that thyroxine (T4) levels were significantly decreased, while triiodothyronine (T3) concentrations were significantly increased after exposure to HEX and TEB, indicating thyroid endocrine disruption. Exposure to HEX significantly induced the transcription of all the measured genes (i.e., corticotrophin-releasing hormone (CRH), thyroid-stimulating hormone (TSHβ), sodium/iodide symporter (NIS), transthyretin (TTR), uridine diphosphate glucuronosyltransferase (UGT1ab), thyronine deiodinase (Dio1 and Dio2), thyroid hormone receptors (TRα and TRβ) in the HPT axis, but did not affect the transcription of thyroglobulin (TG). However, TEB exposure resulted in the upregulation of all the measured genes, excepting that TG, Dio1and TRα had not changed significantly. The overall results indicated that exposure to HEX and TEB could alter thyroid hormone levels as well as gene transcription in the HPT axis in zebrafish larvae.