Despite the abundant literature on the adverse effects of Bisphenol A (BPA) as endocrine disruptor, its toxicity mechanisms are still poorly understood. We present here a study of its effects on the zebrafish eleutheroembryo transcriptome at concentrations ranging from 0.1 to 4â¯mgâ¯L-1, this latter representing the lowest observed effect concentration (LOEC) found in our study at three different macroscopical endpoints (survival, hatching and swim bladder inflation). Multivariate data analysis methods identified both monotonic and bi-phasic patterns of dose-dependent responses. Functional analyses of genes affected by BPA exposure suggest an interaction of BPA with different signaling pathways, being the estrogenic and retinoid receptors two likely targets. In addition, we identified an apparently unrelated inhibitory effect on, among others, visual function genes. We interpret our data as the result of a sum of underlying, independent molecular mechanisms occurring simultaneously at the exposed animals, well below the macroscopic LOEC, but related to at least some of the observed morphological alterations, particularly in eye size and yolk sac resorption. Our data supports the idea that the physiological effects of BPA cannot be only explained by its rather weak interaction with the estrogen receptor, and that multivariate analyses are required to analyze the effects of toxicants like BPA, which interact with different cellular targets producing complex phenotypes.
KEYWORDS: ANOVA-PLS, BPA, Differentially expressed genes (DEGs), Hormone-response, Obesogens, RNA-Seq