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.

Responsive genes for fish embryos have been identified so far for some endocrine pathways but not for androgens. Using transcriptome analysis and multiple concentration-response modeling, we identified putative androgen-responsive genes in zebrafish embryos exposed to 0.05-5000 nM 11-ketotestosterone for 24 h. Four selected genes with sigmoidal concentration-dependent expression profiles (EC50 = 6.5-30.0 nM) were characterized in detail. The expression of cyp2k22 and slco1f4 was demonstrated in the pronephros; lipca was detected in the liver, and sult2st3 was found in the olfactory organs and choroid plexus. Their expression domains, the function of human orthologs, and a pathway analysis suggested a role of these genes in the metabolism of hormones. Hence, it was hypothesized that they were induced to compensate for elevated hormone levels. The induction of sult2st3 and cyp2k22 by 11-ketotestosterone was repressed by co-exposure to the androgen receptor antagonist nilutamide supporting a potential androgen receptor mediated regulation. Sensitivity (expressed as EC50 values) of sult2st3 and cyp2k22 gene expression induction after exposure to other steroidal hormones (11-ketotestosterone ∼ testosterone > progesterone > cortisol > ethinylestradiol) correlated with their known binding affinities to zebrafish androgen receptor. Hence, these genes might represent potential markers for screening of androgenic compounds in the zebrafish embryo.

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.