RNA-Seq analysis of transcriptome responses in Atlantic cod (Gadus morhua) precision-cut liver slices exposed to benzo[a]pyrene and 17α-ethynylestradiol

Polycyclic aromatic hydrocarbons such as benzo[a]pyrene (BaP) that activate the aryl hydrocarbon receptor (Ahr) pathway, and endocrine disruptors acting through the estrogen receptor pathway are among environmental pollutants of major concern. In this work, we exposed Atlantic cod (Gadus morhua) precision-cut liver slices (PCLS) to BaP (10 nM and 1000 nM), ethynylestradiol (EE2) (10 nM and 1000 nM), and equimolar mixtures of BaP and EE2 (10 nM and 1000 nM) for 48 h, and performed RNA-Seq based transcriptome mapping followed by systematic bioinformatics analyses. Our gene expression analysis showed that several genes were differentially expressed in response to BaP and EE2 treatments in PCLS. Strong up-regulation of genes coding for the cytochrome P450 1a (Cyp1a) enzyme and the Ahr repressor (Ahrrb) was observed in BaP treated PCLS. EE2 treatment of liver slices strongly up-regulated genes coding for precursors of vitellogenin (Vtg) and eggshell zona pellucida (Zp) proteins. As expected, pathway enrichment and network analysis showed that the Ahr and estrogen receptor pathways are among the top affected by BaP and EE2 treatments, respectively. Interestingly, two genes coding for fibroblast growth factor 3 (Fgf3) and fibroblast growth factor 4 (Fgf4) were up-regulated by EE2 in this study. To our knowledge, the fgf3 and fgf4 genes have not previously been described in relation to estrogen signaling in fish liver, and these results suggest the modulation of the FGF signaling pathway by estrogens in fish. The signature expression profiles of top differentially expressed genes in response to the single compound (BaP or EE2) treatment were generally maintained in the expression responses to the equimolar binary mixtures. However, in the mixture-treated groups, BaP appeared to have anti-estrogenic effects as observed by lower number of differentially expressed putative EE2 responsive genes. Our in-depth quantitative analysis of changes in liver transcriptome in response to BaP and EE2, using PCLS tissue culture provides further mechanistic insights into effects of the compounds. Moreover, the analyses demonstrate the usefulness of PCLS in cod for omics experiments.