Dynamics in zebrafish development define transcriptomic specificity after angiogenesis inhibitor exposure

Testing for developmental toxicity is an integral part of chemical regulations. The applied tests are laborious and costly and require a large number of vertebrate test animals. To reduce animal numbers and associated costs, the zebrafish embryo was proposed as an alternative model. In this study, we investigated the potential of transcriptome analysis in the zebrafish embryo model to support the identification of potential biomarkers for key events in developmental toxicity, using the inhibition of angiogenesis as a proof of principle. Therefore, the effects on the zebrafish transcriptome after exposure to the tyrosine kinase inhibitors, sorafenib (1.3 µM and 2.4 µM) and SU4312 (1 µM, 2 µM, and 5 µM), and the putative vascular disruptor compound rotenone (25 nM and 50 nM) were analyzed. An early (2 hpf-hours post fertilization) and a late (24 hpf) exposure start with a time resolved transcriptome analysis was performed to compare the specificity and sensitivity of the responses with respect to anti-angiogenesis. We also showed that toxicodynamic responses were related to the course of the internal concentrations. To identify differentially expressed genes (DEGs) the time series data were compared by applying generalized additive models (GAMs). We observed mainly unspecific developmental toxicity in the early exposure scenario, while a specific repression of vascular related genes was only partially observed. In contrast, differential expression of vascular-related genes could be identified clearly in the late exposure scenario. Rotenone did not show angiogenesis-specific response on a transcriptomic level, indicating that the observed mild phenotype of angiogenesis inhibition may represent a secondary effect.