A transcriptomics-based biological framework for studying mechanisms of endocrine disruption in small fish species

This study sought to construct a transcriptomics-based framework of signal transduction pathways, transcriptional regulatory networks, and the hypothalamic-pituitary gonadal (HPG) axis in zebrafish (Danio rerio) to facilitate formulation of specific, testable hypotheses regarding the mechanisms of endocrine disruption in fish. For the analyses involved, we used data from a total of more than 300 microarrays representing 58 conditions, which encompassed 4 tissue types from zebrafish of both genders exposed for 1 of 3 durations to 10 different test chemicals (17α-ethynyl estradiol, fadrozole, 17β-trenbolone, fipronil, prochloraz, flutamide, muscimol, ketoconazole, trilostane, and vinclozolin). Differentially expressed genes were identified by one class t-tests for each condition, and those with false discovery rates of less than 40% and treatment/control ratios ≥1.3-fold were mapped to orthologous human, mouse, and rat pathways by Ingenuity Pathway Analysis to look for overrepresentation of known biological pathways. To complement the analysis of known biological pathways, the genes regulated by approximately 1800 transcription factors were inferred using the ARACNE mutual information-based algorithm. The resulting gene sets for all transcriptional factors, along with a group of compiled HPG-axis genes and approximately 130 publicly available biological pathways, were analyzed for their responses to the 58 treatment conditions by Gene Set Enrichment Analysis (GSEA) and its variant, Extended-GSEA. The biological pathways and transcription factors associated with multiple distinct treatments showed substantial interactions among the HPG-axis, TGF-beta, p53, and several of their cross-talking partners. These candidate networks/pathways have a variety of profound impacts on such cellular functions as stress response, cell cycle, and apoptosis.