Skip to main content

Gene expression responses in male fathead minnows exposed to binary mixtures of an estrogen and antiestrogen



Aquatic organisms are continuously exposed to complex mixtures of chemicals, many of which can interfere with their endocrine system, resulting in impaired reproduction, development or survival, among others. In order to analyze the effects and mechanisms of action of estrogen/anti-estrogen mixtures, we exposed male fathead minnows (Pimephales promelas) for 48 hours via the water to 2, 5, 10, and 50 ng 17α-ethinylestradiol (EE2)/L, 100 ng ZM 189,154/L (a potent antiestrogen known to block activity of estrogen receptors) or mixtures of 5 or 50 ng EE2/L with 100 ng ZM 189,154/L. We analyzed gene expression changes in the gonad, as well as hormone and vitellogenin plasma levels.


Steroidogenesis was down-regulated by EE2 as reflected by the reduced plasma levels of testosterone in the exposed fish and down-regulation of genes in the steroidogenic pathway. Microarray analysis of testis of fathead minnows treated with 5 ng EE2/L or with the mixture of 5 ng EE2/L and 100 ng ZM 189,154/L indicated that some of the genes whose expression was changed by EE2 were blocked by ZM 189,154, while others were either not blocked or enhanced by the mixture, generating two distinct expression patterns. Gene ontology and pathway analysis programs were used to determine categories of genes for each expression pattern.


Our results suggest that response to estrogens occurs via multiple mechanisms, including canonical binding to soluble estrogen receptors, membrane estrogen receptors, and other mechanisms that are not blocked by pure antiestrogens.


Worldwide, aquatic organisms are exposed to mixtures of chemicals (e.g., pharmaceuticals, pesticides, and industrial chemicals), which enter the environment through wastewaters as well as other sources. Many of these chemicals are capable of interfering with endocrine signaling via a complex array of biomolecules (e.g., hormones) to regulate processes such as reproduction and metabolism. These endocrine disrupting chemicals (EDCs) alter signaling through a variety of mechanisms including binding to soluble sex hormone receptors or membrane receptors and acting as agonists or antagonists, or by inhibiting/inducing enzymes and proteins, which produce naturally occurring steroid hormones. Compared to other chemical pollutants, EDCs are likely to have effects at relatively low concentrations [1].

Of the EDCs, xenoestrogens have been the most studied because estrogenic effects have been observed in field studies of fish and wildlife populations [24]. In oviparous animals such as fish, a sensitive and robust biomarker (i.e. vitellogenin, VTG) exists for evaluating exposure to xenoestrogens. Early studies of sewage treatment effluents attributed the feminization of fish to exposure to mixtures of natural (e.g., estrone and 17β-estradiol, E2) and synthetic (e.g., 17α-ethinylestradiol, EE2) estrogens [1, 5]. One of the most potent estrogens known is EE2, a pharmaceutical that is one of the active ingredients in contraceptives. Indeed, EE2 has been shown to be up to 27 times more potent than E2 [6]. In the United States, EE2 use is estimated at 170 kg/yr [7]; and in the United Kingdom, its use is roughly 26 kg/yr [8]. Measured EE2 surface water concentrations in the United States, United Kingdom, The Netherlands, and Germany range from 0.5 to 15 ng/L [7], and it has been frequently measured in United States streams [9].

In laboratory studies, exposures of fish to environmentally relevant EE2 concentrations cause a variety of effects that include testis-ova (the appearance of both sperm and egg follicles in the testis), increased plasma VTG concentrations, reduced gonad size, and altered sex ratios. Studies have used exposure durations of various lengths, including short (≤ 7 days of exposure), intermediate (7 to 28 days exposure), and long (> 28 days) term. In female fish, environmentally relevant EE2 exposures can increase plasma VTG concentrations [1012] and decrease egg production [13] in long-term studies, but seem to have little or no effect on fecundity for intermediate length exposures [10, 12]. In some studies, long-term exposure to EE2 completely inhibits spawning in fish [11, 14].

Long-term EE2 exposure of embryos has been shown to disrupt sexual differentiation of male fish. Fathead minnow (FHM, Pimephales promelas) embryos continuously exposed to EE2 concentrations as low as 4 ng/L did not clearly sexually differentiate at 176 days post-fertilization [12]. Similarly, continuous exposure of zebrafish (Danio rerio) embryos to EE2 concentrations as low as 3 ng/L resulted in all fish having ovaries [11]. EE2 also reduced gonad size and circulating testosterone (T) levels [15], increased VTG [11, 12, 16], and arrested the developmental transition of the gonads of genetically male zebrafish [11]. The steroid also can cause hepatotoxicity, nephrotoxicity and gonadotoxicity [17]. Overall, studies to date suggest that exposure to EE2 elicits adverse effects on fish reproduction primarily through the feminization of male fish, and in females through cessation of spawning. These findings have alerted scientists and environmental regulators to the potential for severe adverse effects on aquatic populations [18], and, potentially, aquatic ecosystems [19]. The current research was conducted to provide a better understanding of the mechanistic basis for effects of estrogenic chemicals in fish.

Effects on gene expression have been investigated with short- and intermediate-term exposures to EE2 [2022] in order to discover gene expression profiles indicative of potential adverse effects. In addition to affecting gene expression through soluble nuclear hormone receptors, it is now clear that sex hormones can also bind directly to membrane receptors and enact immediate changes in signaling via non-genomic pathways [23, 24]. Specific sex hormone receptors in membranes have been identified in fish testis and ovaries for E2 [25, 26], T [27] and progestins [28]. It is difficult to distinguish gene transcription regulation through classical receptor-dependent mechanisms, where estrogen receptor homo- and heterodimers bind to estrogen receptor elements in promoters, from action due to binding of estrogen receptors (ERs) to other transcription factors that activate through Sp1 (stimulatory protein 1) or AP-1 (activating protein 1) binding sites or that activate signaling cascades that start at the membrane. ZM189,154 (ZM) was produced by Astra-Zeneca (Alderly Park, Cheshire, UK) and there are reports that it functions as a "pure" antiestrogen in mammals [29] and in fish [30, 31], meaning that it will bind to and inhibit activation of the ERs in all tissues. But even pure antiestrogens appear to fail in this regard with some genes that are regulated by E2 [32, 33]. ICI 182,780, the most studied pure antiestrogen, can bind to membrane receptors of GnRH-producing GT1-7 cells and displace binding of E2 coupled to bovine serum albumin [34], suggesting that its binding to membrane receptors is inhibited, but it is not clear if this influences all E2 membrane activity [32]. The Atlantic croaker G protein-coupled receptor 30 has been shown to function as a membrane-bound estrogen receptor and its function is agonized by ICI 182,780 [35]. Other E2 activated pathways may not be inhibited by ICI 182,780, as has been shown for E2-stimulated gene regulation through an SP1 site [33]. ZM interactions with membrane receptors have not been studied.

Unlike mammalian species, as many as three to four different ERs have been identified in teleost fish [31, 3638] making evaluation of gene regulation by different ER isotypes even more challenging to understand than in mammalian systems. Using in vitro transfection experiments for largemouth bass (Micropterus salmoides) ERs, we have determined that ZM is equally effective at antagonizing the three soluble receptors [31].

A few studies have investigated the effects of estrogenic mixtures on fish [20, 39, 40] and the binary mixture of E2 with tamoxifen and letrizole, two antiestrogens [41]. However, no studies in fish have investigated the effects of a mixture of EE2 with the potent anti-estrogen, ZM. In this study, the objective was to determine changes in steroidogenesis and in gene expression profiles associated with different exposures by exposing adult male FHM to aqueous doses of EE2, (2, 5, 10 and 50 ng/L); to the pure antiestrogen, ZM (100 ng/L); and to mixtures of EE2 and ZM. The hypothesis we tested was that ZM in the mixture would block the action of EE2 on soluble ERs in the FHM gonad and effectively block gene expression changes observed with EE2 alone.


Water Chemistry

Two distinct experiments were performed. In Exp 1, FHM were treated with three concentrations of EE2 (2, 10 and 50 ng EE2/L), 100 ng ZM/L or a mixture of 50 ng EE2 with 100 ng ZM/L. In Exp 2, FHM were treated with vehicle, 5 ng EE2/L or with a mixture of 5 ng EE2 with 100 ng ZM/L. Water concentrations of EE2 alone and in the mixture were close to target values but decreased after 24 h when they were again renewed to target concentrations (Table 1). Actual concentrations of ZM were not measured.

Table 1 Chemical analysis of water exposures.

Biological responses

There were no mortalities in any of the treatments. Changes in plasma T and VTG were assessed only for a subset of the exposures for Exp 1 (10 and 50 ng EE2/L, 100 ng ZM/L and the mixture of 50 ng EE2/L and 100 ng ZM/L) and only plasma VTG was assessed for exposures for Exp 2. Within 48 h, plasma T levels in males were dramatically reduced in all treatments that were measured for Exp 1 (Figure 1A). In the same time frame there was a significant increase in plasma VTG for the two EE2 concentrations tested, and for the mixture of 50 ng EE2/L and 100 ng ZM/L (Figure 1B). Exposure to 100 ng ZM/L alone did not induce VTG. In the second experiment plasma VTG was significantly up-regulated for the 5 ng EE2/L and for the mixture of 5 ng EE2 with 100 ng ZM/L (Figure 1C).

Figure 1
figure 1

Phenotypic anchoring measurements for male fathead minnows. (A and B) Experiment 1. (A) Plasma T concentrations (ng/mL), (B) Plasma vitellogenin concentrations (μg/mL) in the same fish samples. (Cntrl) triethylene glycol control, (EE2-10) 10 ng EE2/L, (EE2-50) 50 ng EE2/L, (ZM-100) 100 ng ZM/L, (Mix 1:2) 50 ng EE2/L and 100 ng ZM/L. (C) Experiment 2. Plasma VTG concentrations (μg/mL). (EE2-5) 5 ng EE2/L; (Mix 1:20) 5 ng EE2/L and 100 ng ZM/L. Significance ** P ≤ 0.001 and * P ≤ 0.05.

Microarray Results

As described in the Methods section, two microarray experiments were performed, one using testis from FHM exposed to 50 ng EE2/L, 100 ng ZM/L and a combination of both and another using testis from FHM exposed to 5 ng EE2/L or to a combination of 5 ng EE2 with 100 ng ZM/L. The rest of the samples from other EE2 doses were reserved for the quantitative real time PCR (Q-PCR) experiments described below. Samples from the first experiment were analyzed using a 2,000 gene oligonucleotide microarray, and the results are shown in the two sided hierarchical cluster in Figure 2A. The heat map represents genes differentially expressed (p < 0.01) between testis of vehicle control and treated fish. We analyzed four biological samples for each of the exposures; each column in Figure 2A represents one of the samples. As expected, control fish clustered together, whereas fish treated with EE2 alone or with a combination of EE2 and ZM formed a different cluster. Exposure to ZM alone showed the least difference compared to solvent controls; however, even in this comparison there were some differences, suggesting that ZM can influence up- and down-regulation of gene expression in males. There were minor differences between the non-solvent and solvent controls (data not shown).

Figure 2
figure 2

Bi-directional hierarchical cluster analysis of gene expression changes. Green indicates down-regulation relative to control and red indicates up-regulation relative to control. Fathead minnows were exposed to (A) 50 ng EE2/L (EE2), 100 ng ZM/L (ZM), a mixture of 50 ng EE2/L and 100 ng ZM/L (Mix), or TEG control (Cntrl). Array analysis was on the 2 K array. (B) 5 ng EE2/L (EE2), a mixture of 5 ng EE2/L with 100 ng ZM/L (EE2/ZM) or TEG control (Cntrl). Array analysis was on the 22 K array. Top, clustering was performed by treatment; side, clustering was performed by gene. Each column represents a different array.

Exposure to 50 ng EE2/L caused many differences in gene expression. The mixture of 100 ng ZM/L and 50 ng EE2/L reversed the change for several genes affected by EE2 alone, but at this 2:1 ratio the antiestrogen concentration seemed insufficient to totally block the effects of EE2.

Based on these initial results, we conducted a second study, this time using 5 ng EE2/L and a mixture of 100 ng ZM/L with 5 ng EE2/L (a ratio of 20:1; Figure 2B). For this analysis, we used a newer 22,000 gene array that had subsequently become available. Exposure to 5 ng EE2/L increased plasma VTG (Figure 1C), while the 20-fold excess of ZM in the mixture did not affect this increase. A group of 173 genes was altered (p < 0.01) after exposure to either 5 ng EE2/L or to the mixture of 5 ng EE2/L and 100 ng ZM/L (Figure 3). These changes are plotted in order of their degree of expression change for EE2 (Figure 3A), with 83 genes up-regulated and 90 genes down-regulated. Keeping the same order of genes, their fold-expression is plotted for the mixture (Figure 3B). It is clear from this graph that while ZM blocks the EE2 effects for some genes, it does not do so for all. There also appears to be a few genes in the middle of this distribution that are significantly altered only by the mixture and not by EE2 alone.

Figure 3
figure 3

Comparison of overall gene regulation. (A) 5 ng EE2/L and (B) a mixture of 5 ng EE2/L and 100 ng ZM/L as determined by the 22 K array. The genes were ordered according to their expression level in the EE2 treatment (p < 0.01) and represent median expression values of the four arrays for each condition.

Of the 173 regulated genes, 71 genes were modulated by EE2 and blocked by ZM (i.e. reduced expression relative to EE2 alone) in the mixture treatment (Figure 4A and 4B). These genes are likely directly regulated by one or more of the soluble ERs and include "cellular processes involved in calcium-dependent cell-cell adhesion," "sugar transporters," "gonadal mesoderm development," "protein repair," and "proteolysis and gas transport" (see Additional file 1). Expression of the remaining 102 genes modulated by EE2 was either not affected or was enhanced in either direction by the addition of ZM (Figure 4C and 4D). Many of these genes appear to be involved in signaling cascades, as well as other functions such as "peptide crosslinking," "amino acid biosynthesis and metabolism," "regulation of the immune response," "lipid modification," or "response to stress and to radiation" (see Additional file 2).

Figure 4
figure 4

Competitive and non-competitive blocking of gene expression by ZM. (A) Genes differentially changed by 5 ng EE2/L and (B) blocked by the treatment with the mixture of 5 ng EE2/L and 100 ng ZM/L (p < 0.01). (C) Genes differentially changed by 5 ng EE2/L but (D) whose expression was either not changed or enhanced by the treatment with the mixture of 5 ng EE2/L and 100 ng ZM/L. Genes are plotted in order of their change in expression with EE2 (p < 0.01).

Quantitative real-time reverse transcriptase PCR (Q-PCR)

Genes that were tested by Q-PCR (Figure 5) were used to both validate the arrays and to focus on genes whose protein products are involved in steroidogenesis and were expected to be affected by EE2 [42]. Of the genes tested, steroidogenic acute regulatory protein (StAR), cholesterol side-chain cleavage enzyme (P450scc), cytochrome P450 17α hydroxylase, 17,20 lyase (CYP17) and inhibin were significantly down-regulated by 2 to 50 ng EE2/L. Genes for hydroxysteroid dehydrogenases (HSDs) 3β-HSD and 11β-HSD and cytochrome P450 aromatase A-isoform (CYP19A) were not significantly altered, but 11β-HSD and CYP19A showed a downward tendency.

Figure 5
figure 5

Q-PCR analysis of mRNAs for proteins involved in steroidogenesis. Q-PCR results are expressed as fold difference compared to control. Panel A, StAR; Panel B, P450scc; Panel C, CYP17; Panel D, Inhibin; Panel E, CYP19; Panel F, 3β-HSD; Panel G, 11β-HSD; Panel H, model for steroidogenesis. Green boxes refer to mRNAs that are significantly decreased by the treatment in accordance with the Q-PCR graphs illustrated within the panels. Yellow boxes refer to mRNAs that are not significantly changed by the treatments. FHM were treated with TEG, 2, 5, 10 or 50 ng EE2/L. StAR, steroidogenic acute regulatory protein, P450scc, Cytochrome P450 side chain cleavage enzyme, CYP17, Cytochrome P450 17, CYP 19, gonadal aromatase, inhibin, hydroxysteroid dehydrogenases including 3β-HSD and 11β-HSD.

Functional Analysis

While it is interesting to identify individual genes regulated by EE2, most biological processes occur through functional pathways. To assess this, we first assigned as many of the FHM genes as possible to GO categories and to human homologs and then used this information to visualize pathways via Pathway Studio®, software from Ariadne Genomics (Rockville, MD, USA). Of the 1,048 genes regulated by any treatment (p < 0.05), we were able to assign GO categories to 684 genes (65%). Of these we were able to assign human homologs to 536 genes (51% of the original group).

Because of its environmental significance, we focused on the 5 ng EE2/L data for GO analyses. The data set was reduced by statistically determining GO categories for biological processes that were over-represented among the regulated genes which are arranged by increasing p-value (up to 0.05) in Table 2. Since GO categories are listed in a hierarchical format, we removed higher order categories if a lower category was present. We found 39 GO biological process categories up-regulated and 51 categories down-regulated. The most significantly up-regulated GO category was "development," while "DNA replication," "response to radiation," "mutagenesis," "DNA repair," "response to light stimulus," "response to DNA damage stimulus," "DNA metabolism," and "response to endogenous stimulus" were the most significantly down-regulated categories.

Table 2 GO biological processes that are regulated by 5 ng/L EE2



We chose three test concentrations of EE2 (2, 5, and 10 ng EE2/L) with known environmental relevance, and one concentration (50 ng EE2/L) higher than normally seen in the environment [43, 44]. In our experiments, 10 and 50 ng EE2/L decreased plasma T levels, while 5 to 50 ng EE2/L increased plasma VTG concentrations in male fish. To our best knowledge, ZM is not present in the environment, although it represents a potentially important mechanism of action, ER antagonism [45]. The concentration we used, 100 ng ZM/L, and the time of exposure, 48 h, are lower and shorter, respectively than in most other studies [46, 47] where ZM has been shown to have effects in fish. We chose 100 ng ZM/L to attempt to discern intermediate effects on sensitive genes.

In our study, ZM treatment alone or in the mixture with EE2 decreased plasma T levels after 48 h but alone it did not induce plasma VTG concentrations nor did it inhibit the increase in VTG induced by EE2 in the mixture in males. In a study by Panter et al [47] ZM significantly decreased VTG after 4 d in E2-treated juvenile FHM, but only at a concentration of 76 μg/L, a concentration almost 100-times greater than tested in our experiment. In mummichog (Fundulus heteroclitus), there was decrease of plasma T levels in males exposed for 7 days to 250 ng ZM/L but not when treated with 100 ng ZM/L; in those studies there were no effects on VTG levels in males or females with as much as 1,000 ng ZM/L [46].

The Q-PCR data on mRNAs for specific enzymes involved in the biosynthesis of T suggest that the depression of plasma T levels may have occurred directly at the level of steroidogenesis, possibly by direct ER-mediated control of promoters. In the case of CYP17, its down-regulation was blocked by ZM (microarray data), suggesting that it may be regulated via ERs.

Pathway Analysis

Pathway Studio® [48] was used to visualize changes in gene expression from exposure to 5 ng EE2/L, or to the mixture of 5 ng EE2/L and a 100 ng ZM/L. This software can be used effectively to compare expression changes with the much larger database of human protein interactions, but only if gene identities are converted to their human homologs. Important caveats for this type of analysis are that there may be many fish genes for which there are no human homologs (e.g., VTG), and some genes in fish belonging to gene families conserved in mammals may actually function differently in fish due to chromosomal duplications. Given these caveats, this type of analysis can help visualize interactions among gene products and their localization in cellular compartments and assist in the formulation of hypotheses that can be tested in future research.

An interactome is defined as a set of genes whose protein products are functionally linked together either by direct binding, regulation of activity, regulation of expression, promoter binding, protein modification or molecular transport [48]. Using the databases available in PubMed at NCBI (National Center for Biotechnology Information,, we have identified interactomes for both the genes that were regulated by EE2 and then blocked by the combination of EE2 and ZM (Figure 6A – called "competitive interactome" in the discussion below), and those that were regulated by EE2 and either not affected by ZM or enhanced (in either direction) by the combination (Figure 6B – called "non-competitive interactome"). We analyzed these separately in order to determine the types of genes that were included in each.

Figure 6
figure 6

Pathway Studio analysis. (A) Competitive interactome consisting of genes that are regulated by 5 ng EE2/L and blocked by the mixture of 5 ng EE2/L and 100 ng ZM/L. (B) Non-competitive interactome consisting of genes that are regulated by 5 ng EE2/L but either not changed or enhanced by the mixture of 5 ng EE2/L and 100 ng ZM/L. Red indicates up-regulation, green indicates down-regulation.

To reduce the complexity of the data, we only examined genes whose products had direct interactions with other entities. We were only able to assign human homologs to about half of the regulated genes, thus our data set underestimates the genes that are directly linked. Entities that showed no linkages to other entities were removed from the figures, but all entities for which we have human homologs are listed in additional files 1 and 2. While it is possible to allow missing entities in the figures in an effort to link all of the entities, this was not attempted because we wanted to exemplify direct interactomes for which there were expression data. Pathway Studio® assigns gene products to cellular compartments depending on their cellular GO terms.

Many of the genes that are found in the "competitive interactome" are known to be regulated by E2 and antagonized by estrogen antagonists such as ICI 182,780 in mammalian systems. These genes fit a classical pattern of regulation via soluble ERs. For some of these genes, there is evidence that they contain EREs in their promoters in mammalian systems. For example, angiotensinogen (AGT) is highly prominent in the example shown (Figure 6A) and serves as a node for the "competitive interactome". AGT is normally secreted from the liver into the blood, but there are reports indicating secretion by other organs as well [49]. Our data suggest that it is also produced by the testis of FHM. AGT helps to control blood pressure and, as illustrated by the large number of interactions in Figure 6A, can interact with other proteins in a complex way. In mammals, AGT contains an ERE in its promoter which is up-regulated by both natural and synthetic estrogenic steroids [50]. Thus, the observation that the EE2-enhanced expression of this gene is blocked by the EE2-ZM mixture lends support to the assumption that the competitive interactome includes genes directly regulated by soluble ERs. The other genes in this interactome have also been implicated in E2 signalling in mammalian systems or in cell culture, but there are insufficient data in the literature to determine whether they are all regulated directly by soluble ERs. Furthermore, depending on the tissue the direction of regulation may differ from what we observed in the testis of FHM.

Endothelin receptor type A (EDNRA) is up-regulated by E2 during the proliferative phase of the endometrial glandular epithelium [51]. There are no reports on whether or not this up-regulation is through an ERE in the promoter. Different from our results, clusterin (CLU) is down-regulated by E2 in rat endometrium, but this regulation is reversed by tamoxifen, another ER antagonist [52]. CLU is known to be regulated by TGF beta and c-fos through an AP-1 site [53]. CLU is a glycoprotein also known as testosterone-repressed prostate message-2 [54]. This gene is expressed in mammalian testis and apparently has many roles including involvement in apoptosis of the seminal vesicle [55]. Complement component factor H (Mud1) contains an imperfect palindrome motif in its promoter in L cells that is present in EREs [56], suggesting that this gene could be directly regulated by soluble ERs. Inhibitor of DNA binding 2 (ID2) has been shown to be down regulated by E2 in MCF-7 cells [57] and is directly related to the down-regulation of inhibin alpha, which in turn has a role inhibiting the secretion of FSH from pituitary gonadotrophs [58].

Several genes interconnected with this set were down-regulated by EE2 and blocked by ZM in our experiment. For example, receptor (calcitonin) activity modifying protein 1 (RAMP1) has been shown to be down-regulated by E2 in rat placenta [59]. Different from our study, some of the down-regulated genes in the FHM testis have been shown to be up-regulated by E2 in various mammalian tissues. Pleiotrophin (PTN) is up-regulated in human endometrial epithelial cells [60], and annexin A1 (ANXA1) is up-regulated in a lymphoblastic leukemia cell line [61]. Caspase 8 (CASP8) is regulated by activation of human ERβ but not by ERα [62]. Tissue inhibitor of metalloproteinase 3 (TIMP3) is increased in breast cancer cell growth [63]. It is not clear why the direction of regulation is different in FHM testis, but this may be a tissue specific effect.

In the case of the "non-competitive interactome" (Figure 6B), the genes were differentially expressed in response to EE2 exposure and either were not affected by ZM or further amplified by ZM in the mixture. We did not expect to see many genes in this category. This type of effect could be due to activation by E2 on non-canonical response elements, as recently demonstrated in transgenic mice expressing a reporter construct containing SP1 sites [33] or by secondary effects that may have occurred in the 48 h timeframe. In Figure 6B, we have accentuated one pathway to illustrate this effect. Prostaglandin-endoperoxide synthase 2 (PTGS2) is a central node in this figure and is induced 1.5 fold by EE2 alone but not changed appreciably by the mixture (1.8 fold change). PTGS2 is involved in the synthesis of prostaglandins from arachidonic acid and is influenced by E2 in mammalian tissues [64]. Important in this set of genes is the gene for "signal transducer and activator of transcription 1" (STAT1) which is a critical transcription factor involved in the JAK-STAT signalling pathway central for innate immunity [65] and apoptosis [66], among other functions. This transcription factor is also activated via the retinoic acid receptor signalling pathway [67], thus bridging both the E2 and retinoic acid pathways. Also important is histone deacetylase 9 (HDAC9), a gene product that is involved in chromatin remodelling, allowing access of transcription factors to regions in DNA. No information exists regarding the influence of estrogen on HDAC9 in mammalian tissues but inactivation of other histone deacetlyases is an important step for ER activation in cell lines no longer responsive to E2 [68]. Tnf receptor-associated factor 6 (TRAF6) is a protein known to be involved in signal transduction through membrane receptors [69], and RAS p21 protein activator 1 appears to play a role in Ras GTPase mediated signal transduction [70].

Among the down-regulated genes were several transcription factors, including transcription factor 12 (TCF12) involved in control of immunoglobulin transcription [71], transcription factor 8 (TCF8), a negative regulator of cadherin [72], myeloblastosis oncogene (MYB), known to be involved in estrogen signalling in some breast cancer cells [73] and AF4/FMR2 family, member 4 (AF5Q31), which functions as a transcriptional regulator in testicular somatic cells, essential for male germ cell differentiation and survival [74]. Karyopherin (importin) alpha 4 (KPNA4) helps import proteins into the nucleus during spermatogenesis [75]. RAN binding protein 3 isoform b (RANBP3) links the RAS and PI3-kinase signaling pathways with nuclear transport [76] and Eph receptor A2 (EPHA2), is a protein in the tyrosine kinase family [77]. The roles of other genes in this interactome are listed in additional file 2.

It is not clear at this time why so many genes have escaped antagonism by ZM in the mixture exposure. It is possible that they are activated via non-canonical ER interactions with other transcription factors [33] or are the results of activation through G protein-coupled receptor 30 (GPR30), a membrane-bound estrogen receptor [35]. A large percentage of genes in the non-competitive interactome function in non-genomic signaling pathways, raising the possibility that these genes are all regulated via membrane receptors which escape antagonism. Further work will be required to sort out exactly how each of these genes is regulated.


We used genomics to try to elucidate the mechanisms of action of estrogenic and anti-estrogenic compounds and their potential effects on aquatic organisms. Our data provides some insight into the estrogen-regulated effects, suggesting that response to estrogens occurs via different mechanisms. The use of an estrogen/antiestrogen mixture provides a distinction among different modes of action of estrogenic compounds: through canonical binding to soluble ERs; membrane ERs; or some other potential mechanisms that may not be blocked by pure antiestrogens.


Fish Exposure and Tissue Collection

Reproductively-mature, pond-reared FHM were purchased from Andersen Minnow Farm, AR, 4 days prior to starting the experiment. Upon arrival, the fish were treated for parasites and bacteria by a prophylactic salt-water dip (3%, 1 min). Males were separated from the population the following day, and acclimated in the treatment aquaria for 48 h. The water used for this study was carbon-filtered, dechlorinated tap water.

The exposure system consisted of 40 L glass aquaria. Each exposure was conducted in quadruplicate and each aquarium contained eight male FHM in 25 L of treatment water. Test chemicals for each treatment group (100 L for 4 aquaria) were prepared in separate (by treatment) 250 L fiberglass tanks the day of exposure. Aquaria were equilibrated with test chemicals for 24 h prior to the introduction of fish. Test solutions were renewed to 90% of the 25 L exposure volume after 24 h and the exposure was ended at 48 h. The positions of the treatment tanks were randomized and test initiation times were staggered to ensure an exposure/sampling interval of 48 h. The fish were not fed the day before and during the experiment. Temperature was maintained at 25°C with a photoperiod of 16 h light: 8 h dark.

Exposure Solutions

EE2 was purchased from Sigma Chemical Company (St. Louis, MO). ZM189,154 was a generous gift from AstraZeneca. Working solutions for each test chemical consisted of 1 mg/ml test compound in 70% triethylene glycol (TEG) and 28.5% ethanol. This working solution was further diluted to make stock solutions for each treatment (nominal concentrations of 2, 5, 10 and 50 ng EE2/L and 100 ng ZM/L and mixtures containing 5 ng EE2/L or 50 ng EE2/L and 100 ng ZM/L), so as to maintain a concentration of 50 μl TEG/L of test water. EE2 concentrations spanned the environmentally relevant levels (2–10 ng EE2/L) to a concentration higher than would typically occur in the environment. The antiestrogen ZM concentration was chosen to be higher than EE2 so that it could effectively block action on the ERs. The concentration of ZM in the two mixtures was 100 ng ZM/L.

Water samples were collected at the start of the exposures (0 hr), and after 24 h (both before and after change of tank water) and after 48 h. A sample of the test solution (1 L) was collected in an amber glass bottle with a teflon cap and stored at 4°C. The water was passed through an AccuBond II ODS-C18 solid phase extraction column (Agilent, Palo Alto, CA) and the EE2 was eluted with 5 ml dichloromethane. After drying, the EE2residue was reconstituted in 1 ml of buffer and analyzed using an enzyme-linked immunosorbent assay (ELISA) kit (Abraxis, Los Angeles, CA), following the manufacturer's instructions. The detection limit for this assay is 50 ng EE2/L in the reconstituted solution. We were not able to determine actual concentrations of ZM, so we report only nominal concentrations.

All procedures involving live fish were reviewed and approved by the University of Florida Institutional Animal Care and Use Committee (IACUC). At the conclusion of the exposures, fish were anesthetized (MS-222), weighed to the nearest 0.1 g and blood samples were collected from the caudal vasculature for analysis of VTG and T concentrations, as described below. The testes were removed and cut into small pieces. Dissected tissues were flash frozen using liquid nitrogen and stored at -80°C until needed.

Vitellogenin Assay

Plasma concentrations of VTG were determined by ELISA using a monoclonal antibody, 2D3, previously validated for the FHM [78]. The limit of detection for the FHM VTG ELISA in plasma was 0.5 μg/mL. All assays were performed in triplicate and reported as the mean of the three measurements. The coefficient of variation was < 10% for all samples analyzed. Inter and intra-assay variability was routinely measured by analyzing positive controls on several plates and found to be < 10% and < 5%, respectively.

Testosterone Radioimmunoassay

Plasma concentrations of T were measured using a radioimmunoassay (RIA) validated for the FHM based on a slight modification of a previously published protocol [79]. The antibody against T, 20-TR05, was purchased from Fitzgerald Industries International, Concord, MA. Tritiated label ([1,2,6,7-3H] T) was from GE Healthcare (Piscataway, NJ). The T standard (Sigma T-1500) was obtained from Sigma Chemical Company (St. Louis, MO). Plasma samples (12 μL each) were extracted with 2 mL of ethyl ether, as described previously [79]. The extraction efficiency was 93%. Samples were analyzed in duplicate. The intraassay coefficients of variance were generally ≤ 5% and all samples were run in one assay to prevent interassay variability.

Data analysis

Plasma concentrations of T were analyzed by one-way ANOVA, followed by Fisher Protected Least Significant Difference (PLSD) test for post-hoc analysis. All analyses were carried out using StatView 5.0 (SAS Institute, Inc., Cary, NC). Homoscedasticity was assessed using F – tests, and, where necessary (p < 0.05), data were log transformed [80]. All data are reported as nontransformed values, as mean ± SEM, and significance was determined at p-value < 0.05. Plasma VTG concentrations were analyzed by Dunnett's pairwise multiple comparisons on log transformed data.

RNA Extraction

Total RNA was isolated from 30–50 mg FHM gonadal tissue with the RNA Stat-60 reagent (Tel-test, Friendswood, TX), as previously described [81]. Total RNA was treated with DNase and the quality assessed with an Agilent 2100 BioAnalyzer (Agilent, Palo Alto, CA), and the quantity determined on a NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, DE). RNA was stored at -80°C until further use.


Fathead minnow microarrays manufactured by Agilent (Palo Alto, CA) were purchased from EcoArray (Alachua, FL). For the first experiment we used a targeted 2,000 gene array (GPL6516) while for the other we employed a 22,000 gene array (4 × 44 K format, GPL7282). Array hybridizations were performed using a reference design. The reference material, which was used for all studies, consisted of equal amounts of RNA from both female and male tissues (liver, brain and gonad). Four replicates consisting of four different individuals were analyzed for each of the treatments (solvent (TEG) control, non-solvent control, EE2, ZM, EE2/ZM). The cDNA synthesis, cRNA labeling and hybridization were performed following the manufacturer's kits and protocols (Agilent Low RNA Input Fluorescent Linear Amplification Kit and Agilent 60-mer oligo microarray processing protocol; Agilent, Palo Alto, CA). The gonad samples were labeled with Cy5 while the reference sample was labeled with Cy3. Once the labeling was complete, samples were hybridized to the microarray using conditions recommended by the manufacturer. After hybridizing for 17 h, microarrays were washed and then scanned with a laser-based detection system (Agilent, Palo Alto, CA). Text versions of the Agilent raw data have been deposited at the Gene Expression Omnibus website (GEO:; Accession series record number GSE14235).


Microarray image processing and data pre-processing were performed using Agilent's Feature Extraction software v 9.5 (Agilent, 2007). The intensity of each spot was summarized by the median pixel intensity. A log2 transformed signal ratio between the experimental channel and the reference channel was calculated for each spot, followed by within-array LOWESS transformation and between array scale normalization on median intensities [82].

Two-way ANOVA was performed on log2 transformed signal ratios of each probe individually, followed by Tukey-HSD pair-wise comparisons to determine genes whose expression was significantly regulated by the treatments. A p-value ≤ 0.05 was used as the cutoff. Genes whose fold expression changes were less than 1.5 fold were eliminated from further analyses irrespective of statistical significance.

GeneOntology (GO) annotations were provided by EcoArray Inc. based largely on homologies between FHM genes and human genes. Overrepresentation of differentially expressed genes in the biological process GO category was determined by Fisher Exact Test with a p-value ≤ 0.05 as a cutoff, and the false discovery rate was determined [83]. PathwayStudio® software [48] from Ariadne Genomics (Rockville, MD) was used to determine the list of common regulators among the genes that were differentially expressed in the experiments.

Real-time Polymerase Chain Reaction (Q-PCR)

Total RNA (1 μg) was reverse transcribed into cDNA using 3 μl random primers (0.1 μg/μl), 0.8 μl dNTP mix (25 mM each dNTP), 2 μl transcription buffer (10×), 1 μl StrataScript RT (50 U/μl), and 0.5 μl RNAse Block (40 U/μl) in a final volume of 20 μl (all reverse transcription reagents were from Stratagene, La Jolla, CA). The resulting cDNA was used as a template for Q-PCR. Specific primers for selected genes were designed to perform Q-PCR (Table 3).

Table 3 Real time PCR primers

Each Q-PCR reaction consisted of 1× iQ SYBR Green Supermix (Bio-Rad, Hercules, CA), 0.4 μM primers and 1 μl of cDNA in a 25 μl reaction. The Q-PCR conditions were 95°C for 3 min and 40 cycles at 95°C for 15 sec and 60°C for 1 min in an iCycler Thermal Cycler (Bio-Rad, Hercules, CA). The Q-PCR results were normalized to 18S rRNA (Applied Biosystems, Foster City, CA) and analyzed using the ΔΔCt method, compared to the vehicle controls. We measured the following mRNAs: cytochrome P450 17α hydroxylase, 17,20 lyase (CYP17) [84], steroidogenic acute regulatory protein (StAR) [85], cholesterol side-chain cleavage enzyme (P450scc) [85], hydroxysteroid dehydrogenases (HSDs) 3β-HSD [85] and 11β-HSD [86], inhibin (INHB) and cytochrome P450 aromatase A-isoform (CYP19A) [87].



AF4/FMR2 family, member 4




annexin A1


activating protein 1


3 beta hydroxysteroid dehydrogenase


11 beta hydroxysteroid dehydrogenase


caspase 8


CBL E3 ubiquitin protein ligase





CYP 17:

cytochrome P450, family 17, subfamily a, polypeptide 1


cytochrome P450 aromatase A-isoform


dual specific phosphatase 6


Endocrine disrupting chemicals


17β estradiol


17α ethinylestradiol


enzyme-linked immunosorbent assay


endothelin receptor type A


Eph receptor A2


estrogen receptor


fathead minnow


growth arrest and DNA-damage-inducible 45 beta


gene ontology


histone deacetylase 9


inhibitor of DNA binding 2


inhibin alpha


karyopherin (importin) alpha 4


myeloblastosis oncogene

Mud 1:

complement component factor H


nuclear receptor subfamily 4, group A, member 3


quantitative real time reverse transcriptase polymerase chain reaction


Cytochrome P450 side chain cleavage enzyme


prostaglandin-endoperoxide synthase 2




purine rich element binding protein A


receptor (calcitonin) activity modifying protein 1


RAN binding protein 3 isoform


RAS p21 protein activator 1




SH2 domain containing 3C


stimulatory protein 1


steroidogenic acute regulatory protein


signal transducer and activator of transcription 1




transcription factor 8


transcription factor 12


tissue inhibitor of metalloproteinase 3


Tnf receptor-associated factor 6




zinc finger protein 36, C3H type-like 1


ZM 189,154.


  1. Sumpter JP: Endocrine disrupters in the aquatic environment: An overview. Acta Hydrochimica et Hydrobiologica. 2005, 33: 9-16. 10.1002/aheh.200400555.

    Article  CAS  Google Scholar 

  2. Damstra T, Barlow S, Bergman A, Kavlock RJ, Kraak Van der G: Global Assessment of the State-of-the-Science of Endocrine Disruptors. International Programme on Chemical Safety, World Health Organization: 2002; Geneva, Switzerland. 2002, 180-

    Google Scholar 

  3. Pickering AD, Sumpter JP: Comprehending endocrine disrupters in aquatic environments. Environmental Science & Technology. 2003, 37 (17): 331A-336A. 10.1021/es032570f.

    Article  CAS  Google Scholar 

  4. Jobling S, Nolan M, Tyler CR, Brighty G, Sumpter JP: Widespread sexual disruption in wild fish. Environmental Science & Technology. 1998, 32 (17): 2498-2506. 10.1021/es9710870.

    Article  CAS  Google Scholar 

  5. Folmar LC, Denslow ND, Rao V, Chow M, Crain DA, Enblom J, Marcino J, Guillette LJ: Vitellogenin induction and reduced serum testosterone concentrations in feral male carp (Cyprinus carpio) captured near a major metropolitan sewage treatment plant. Environ Health Perspect. 1996, 104 (10): 1096-1101. 10.2307/3433123.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Thorpe KL, Cummings RI, Hutchinson TH, Scholze M, Brighty G, Sumpter JP, Tyler CR: Relative potencies and combination effects of steroidal estrogens in fish. Environ Sci Technol. 2003, 37 (6): 1142-1149. 10.1021/es0201348.

    Article  CAS  PubMed  Google Scholar 

  7. Ericson JF, Laenge R, Sullivan DE: Comment on "Pharmaceuticals, Hormones, and Other Organic Wastewater Contaminants in U.S. Streams, 1999–2000: A National Reconnaissance. Environmental Science and Technology. 2002, 36 (18): 4005-4006. 10.1021/es0200903.

    Article  CAS  PubMed  Google Scholar 

  8. Ankley GT, Brooks BW, Huggett DB, Sumpter JP: Repeating history: Pharmaceuticals in the environment. Environmental Science & Technology. 2007, 41 (24): 8211-8217. 10.1021/es072658j.

    Article  CAS  Google Scholar 

  9. Kolpin DW, Furlong ET, Meyer MT, Thurman EM, Zaugg SD, Barber LB, Buxton HT: Pharmaceuticals, hormones, and other organic wastewater contaminants in U.S. streams, 1999–2000: A national reconnaissance. Environ Sci Technol. 2002, 36 (6): 1202-1211. 10.1021/es011055j.

    Article  CAS  PubMed  Google Scholar 

  10. Pawlowski S, van Aerle R, Tyler CR, Braunbeck T: Effects of 17alpha-ethinylestradiol in a fathead minnow (Pimephales promelas) gonadal recrudescence assay. Ecotoxicol Environ Saf. 2004, 57 (3): 330-345. 10.1016/j.ecoenv.2003.07.019.

    Article  CAS  PubMed  Google Scholar 

  11. Fenske M, Maack G, Schafers C, Segner H: An environmentally relevant concentration of estrogen induces arrest of male gonad development in zebrafish, Danio rerio. Environ Toxicol Chem. 2005, 24 (5): 1088-1098. 10.1897/04-096R1.1.

    Article  CAS  PubMed  Google Scholar 

  12. Lange R, Hutchinson TH, Croudace CP, Siegmund F, Schweinfurth H, Hampe P, Panter GH, Sumpter JP: Effects of the synthetic estrogen 17 alpha-ethinylestradiol on the life-cycle of the fathead minnow (Pimephales promelas). Environ Toxicol Chem. 2001, 20 (6): 1216-1227. 10.1897/1551-5028(2001)020<1216:EOTSEE>2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  13. Santos EM, Paull GC, Van Look KJW, Workman VL, Holt WV, Van Aerle R, Kille P, Tyler CR: Gonadal transcriptome responses and physiological consequences of exposure to oestrogen in breeding zebrafish (Danio rerio). Aquatic Toxicology. 2007, 83 (2): 134-142. 10.1016/j.aquatox.2007.03.019.

    Article  CAS  PubMed  Google Scholar 

  14. Parrott JL, Blunt BR: Life-cycle exposure of fathead minnows (Pimephales promelas) to an ethinylestradiol concentration below 1 ng/L reduces egg fertilization success and demasculinizes males. Environ Toxicol. 2005, 20 (2): 131-141. 10.1002/tox.20087.

    Article  CAS  PubMed  Google Scholar 

  15. Martyniuk CJ, Xiong H, Crump K, Chiu S, Sardana R, Nadler A, Gerrie ER, Xia X, Trudeau VL: Gene expression profiling in the neuroendocrine brain of male goldfish (Carassius auratus) exposed to 17alpha-ethinylestradiol. Physiol Genomics. 2006, 27 (3): 328-336. 10.1152/physiolgenomics.00090.2006.

    Article  CAS  PubMed  Google Scholar 

  16. Andersen L, Holbech H, Gessbo A, Norrgren L, Petersen GI: Effects of exposure to 17alpha-ethinylestradiol during early development on sexual differentiation and induction of vitellogenin in zebrafish (Danio rerio). Comp Biochem Physiol C Toxicol Pharmacol. 2003, 134 (3): 365-374. 10.1016/S1532-0456(03)00006-1.

    Article  PubMed  Google Scholar 

  17. Weber LP, Balch GC, Metcalfe CD, Janz DM: Increased kidney, liver, and testicular cell death after chronic exposure to 17alpha-ethinylestradiol in medaka (Oryzias latipes). Environ Toxicol Chem. 2004, 23 (3): 792-797. 10.1897/02-570.

    Article  CAS  PubMed  Google Scholar 

  18. Miller DH, Jensen KM, Villeneuve DL, Kahl MD, Makynen EA, Durhan EJ, Ankley GT: Linkage of biochemical responses to population-level effects: a case study with vitellogenin in the fathead minnow (Pimephales promelas). Environ Toxicol Chem. 2007, 26 (3): 521-527. 10.1897/06-318R.1.

    Article  CAS  PubMed  Google Scholar 

  19. Kidd KA, Blanchfield PJ, Mills KH, Palace VP, Evans RE, Lazorchak JM, Flick RW: Collapse of a fish population after exposure to a synthetic estrogen. Proc Natl Acad Sci USA. 2007, 104 (21): 8897-8901. 10.1073/pnas.0609568104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Filby AL, Neuparth T, Thorpe KL, Owen R, Galloway TS, Tyler CR: Health impacts of estrogens in the environment, considering complex mixture effects. Environ Health Perspect. 2007, 115 (12): 1704-1710.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Filby AL, Thorpe KL, Maack G, Tyler CR: Gene expression profiles revealing the mechanisms of anti-androgen- and estrogen-induced feminization in fish. Aquat Toxicol. 2007, 81 (2): 219-231. 10.1016/j.aquatox.2006.12.003.

    Article  CAS  PubMed  Google Scholar 

  22. Martyniuk CJ, Gerrie ER, Popesku JT, Ekker M, Trudeau VL: Microarray analysis in the zebrafish (Danio rerio) liver and telencephalon after exposure to low concentration of 17alpha-ethinylestradiol. Aquatic Toxicology. 2007, 84 (1): 38-49. 10.1016/j.aquatox.2007.05.012.

    Article  CAS  PubMed  Google Scholar 

  23. Thomas P, Dressing G, Pang Y, Berg H, Tubbs C, Benninghoff A, Doughty K: Progestin, estrogen and androgen G-protein coupled receptors in fish gonads. Steroids. 2006, 71 (4): 310-316. 10.1016/j.steroids.2005.09.015.

    Article  CAS  PubMed  Google Scholar 

  24. Braun AM, Thomas P: Androgens inhibit estradiol-17beta synthesis in Atlantic croaker (Micropogonias undulatus) ovaries by a nongenomic mechanism initiated at the cell surface. Biol Reprod. 2003, 69 (5): 1642-1650. 10.1095/biolreprod.103.015479.

    Article  CAS  PubMed  Google Scholar 

  25. Loomis AK, Thomas P: Effects of estrogens and xenoestrogens on androgen production by Atlantic croaker testes in vitro: evidence for a nongenomic action mediated by an estrogen membrane receptor. Biol Reprod. 2000, 62 (4): 995-1004. 10.1095/biolreprod62.4.995.

    Article  CAS  PubMed  Google Scholar 

  26. Thomas P, Doughty K: Disruption of rapid, nongenomic steroid actions by environmental chemicals: interference with progestin stimulation of sperm motility in Atlantic croaker. Environ Sci Technol. 2004, 38 (23): 6328-6332. 10.1021/es0403662.

    Article  CAS  PubMed  Google Scholar 

  27. Braun AM, Thomas P: Biochemical characterization of a membrane androgen receptor in the ovary of the Atlantic croaker (Micropogonias undulatus). Biol Reprod. 2004, 71 (1): 146-155. 10.1095/biolreprod.103.025825.

    Article  CAS  PubMed  Google Scholar 

  28. Patino R, Thomas P: Characterization of membrane receptor activity for 17 alpha, 20 beta, 21-trihydroxy-4-pregnen-3-one in ovaries of spotted seatrout (Cynoscion nebulosus). Gen Comp Endocrinol. 1990, 78 (2): 204-217. 10.1016/0016-6480(90)90007-9.

    Article  CAS  PubMed  Google Scholar 

  29. Dukes M, Chester R, Yarwood L, Wakeling AE: Effects of a non-steroidal pure antioestrogen, ZM 189,154, on oestrogen target organs of the rat including bones. J Endocrinol. 1994, 141 (2): 335-341. 10.1677/joe.0.1410335.

    Article  CAS  PubMed  Google Scholar 

  30. Marlatt VL, Hewitt LM, Kraak Van Der G: Utility of in vitro test methods to assess the activity of xenoestrogens in fish. Environ Toxicol Chem. 2006, 25 (12): 3204-3212. 10.1897/05-686R.1.

    Article  CAS  PubMed  Google Scholar 

  31. Sabo-Attwood T, Blum JL, Kroll KJ, Patel V, Birkholz D, Szabo NJ, Fisher SZ, McKenna R, Campbell-Thompson M, Denslow ND: Distinct expression and activity profiles of largemouth bass (Micropterus salmoides) estrogen receptors in response to estradiol and nonylphenol. J Mol Endocrinol. 2007, 39 (4): 223-237. 10.1677/JME-07-0038.

    Article  CAS  PubMed  Google Scholar 

  32. Sismondi P, Biglia N, Ponzone R, Fuso L, Scafoglio C, Cicatiello L, Ravo M, Weisz A, Cimino D, Altobelli G, et al: Influence of estrogens and antiestrogens on the expression of selected hormone-responsive genes. Maturitas. 2007, 57 (1): 50-55. 10.1016/j.maturitas.2007.02.019.

    Article  CAS  PubMed  Google Scholar 

  33. Wu F, Xu R, Martin J, Safe S: In vivo profiling of estrogen receptor/specificity protein-dependent transactivation. Endocrinology. 2008, 149 (11): 5696-5705. 10.1210/en.2008-0720.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Morales A, Gonzalez M, Marin R, Diaz M, Alonso R: Estrogen inhibition of norepinephrine responsiveness is initiated at the plasma membrane of GnRH-producing GT1-7 cells. J Endocrinol. 2007, 194 (1): 193-200. 10.1677/JOE-06-0001.

    Article  CAS  PubMed  Google Scholar 

  35. Pang Y, Dong J, Thomas P: Estrogen signaling characteristics of Atlantic croaker G protein-coupled receptor 30 (GPR30) and evidence it is involved in maintenance of oocyte meiotic arrest. Endocrinology. 2008, 149 (7): 3410-3426. 10.1210/en.2007-1663.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Hawkins MB, Thornton JW, Crews D, Skipper JK, Dotte A, Thomas P: Identification of a third distinct estrogen receptor and reclassification of estrogen receptors in teleosts. Proc Natl Acad Sci USA. 2000, 97 (20): 10751-10756. 10.1073/pnas.97.20.10751.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Nagler JJ, Cavileer T, Sullivan J, Cyr DG, Rexroad C: The complete nuclear estrogen receptor family in the rainbow trout: discovery of the novel ERalpha2 and both ERbeta isoforms. Gene. 2007, 392 (1–2): 164-173. 10.1016/j.gene.2006.12.030.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Sabo-Attwood T, Kroll KJ, Denslow ND: Differential expression of largemouth bass (Micropterus salmoides) estrogen receptor isotypes alpha, beta, and gamma by estradiol. Mol Cell Endocrinol. 2004, 218 (1–2): 107-118. 10.1016/j.mce.2003.12.007.

    Article  CAS  PubMed  Google Scholar 

  39. Jukosky JA, Watzin MC, Leiter JC: The effects of environmentally relevant mixtures of estrogens on Japanese medaka (Oryzias latipes) reproduction. Aquat Toxicol. 2008, 86 (2): 323-331. 10.1016/j.aquatox.2007.11.012.

    Article  CAS  PubMed  Google Scholar 

  40. Lin LL, Janz DM: Effects of binary mixtures of xenoestrogens on gonadal development and reproduction in zebrafish. Aquat Toxicol. 2006, 80 (4): 382-395. 10.1016/j.aquatox.2006.10.004.

    Article  CAS  PubMed  Google Scholar 

  41. Sun L, Zha J, Wang Z: Effects of binary mixtures of estrogen and antiestrogens on Japanese medaka (Oryzias latipes). Aquat Toxicol. 2009, 93 (1): 83-89. 10.1016/j.aquatox.2009.03.010.

    Article  CAS  PubMed  Google Scholar 

  42. Villeneuve DL, Larkin P, Knoebl I, Miracle AL, Kahl MD, Jensen KM, Makynen EA, Durhan EJ, Carter BJ, Denslow ND, et al: A graphical systems model to facilitate hypothesis-driven ecotoxicogenomics research on the teleost brain-pituitary-gonadal axis. Environ Sci Technol. 2007, 41 (1): 321-330. 10.1021/es061739x.

    Article  CAS  PubMed  Google Scholar 

  43. Desbrow C, Routledge EJ, Brighty GC, Sumpter JP, Waldock M: Identification of estrogenic chemicals in STW effluent. 1. Chemical fractionation and in vitro biological screening. Environmental Science & Technology. 1998, 32 (11): 1549-1558. 10.1021/es9707973.

    Article  CAS  Google Scholar 

  44. Kashiwada S, Ishikawa H, Miyamoto N, Ohnishi Y, Magara Y: Fish test for endocrine-disruption and estimation of water quality of Japanese rivers. Water Res. 2002, 36 (8): 2161-2166. 10.1016/S0043-1354(01)00406-7.

    Article  CAS  PubMed  Google Scholar 

  45. Dukes M, Chester R, Yarwood L, Wakeling AE: Effects of a nonsteroidal pure antiestrogen, ZM-189,154, on estrogen target organs of the rat including bones. J Endocrinol. 1994, 141 (2): 335-341. 10.1677/joe.0.1410335.

    Article  CAS  PubMed  Google Scholar 

  46. MacLatchy DL, Courtenay SC, Rice CD, Kraak Van der GJ: Development of a short-term reproductive endocrine bioassay using steroid hormone and vitellogenin end points in the estuarine mummichog (Fundulus heteroclitus). Environ Toxicol Chem. 2003, 22 (5): 996-1008. 10.1897/1551-5028(2003)022<0996:DOASRE>2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  47. Panter G, Hutchinson T, Länge R, Lye C, JP S, Zerulla M, Tyler C: Utility of a juvenile fathead minnow screening assay for detecting (anti-) estrogenic substances. Environmental Toxicology and Chemistry. 2002, 21 (2): 319-326. 10.1897/1551-5028(2002)021<0319:UOAJFM>2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  48. Nikitin A, Egorov S, Daraselia N, Mazo I: Pathway studio – the analysis and navigation of molecular networks. Bioinformatics. 2003, 19 (16): 2155-2157. 10.1093/bioinformatics/btg290.

    Article  CAS  PubMed  Google Scholar 

  49. Wong MK, Ge W, Woo NY: Positive feedback of hepatic angiotensinogen expression in silver sea bream (Sparus sarba). Mol Cell Endocrinol. 2007, 263 (1–2): 103-111. 10.1016/j.mce.2006.09.001.

    Article  CAS  PubMed  Google Scholar 

  50. Krattenmacher R, Knauthe R, Parczyk K, Walker A, Hilgenfeldt U, Fritzemeier KH: Estrogen action on hepatic synthesis of angiotensinogen and IGF-I: direct and indirect estrogen effects. J Steroid Biochem Mol Biol. 1994, 48 (2–3): 207-214. 10.1016/0960-0760(94)90146-5.

    Article  CAS  PubMed  Google Scholar 

  51. Samsioe G: The endometrium: effects of estrogen and estrogen-progestogen replacement therapy. Int J Fertil Menopausal Stud. 1994, 39 (Suppl 2): 84-92.

    PubMed  Google Scholar 

  52. Heikaus S, Winterhager E, Traub O, Grummer R: Responsiveness of endometrial genes Connexin26, Connexin43, C3 and clusterin to primary estrogen, selective estrogen receptor modulators, phyto- and xenoestrogens. J Mol Endocrinol. 2002, 29 (2): 239-249. 10.1677/jme.0.0290239.

    Article  CAS  PubMed  Google Scholar 

  53. Jin G, Howe PH: Transforming growth factor beta regulates clusterin gene expression via modulation of transcription factor c-Fos. Eur J Biochem. 1999, 263 (2): 534-542. 10.1046/j.1432-1327.1999.00533.x.

    Article  CAS  PubMed  Google Scholar 

  54. Londou A, Mikrou A, Zarkadis IK: Cloning and characterization of two clusterin isoforms in rainbow trout. Mol Immunol. 2008, 45 (2): 470-478. 10.1016/j.molimm.2007.05.027.

    Article  CAS  PubMed  Google Scholar 

  55. Wong P, Pineault J, Lakins J, Taillefer D, Leger J, Wang C, Tenniswood M: Genomic organization and expression of the rat TRPM-2 (clusterin) gene, a gene implicated in apoptosis. J Biol Chem. 1993, 268 (7): 5021-5031.

    CAS  PubMed  Google Scholar 

  56. Munoz-Canoves P, Vik DP, Tack BF: Mapping of a retinoic acid-responsive element in the promoter region of the complement factor H gene. J Biol Chem. 1990, 265 (33): 20065-20068.

    CAS  PubMed  Google Scholar 

  57. Wang DY, McKague B, Liss SN, Edwards EA: Gene expression profiles for detecting and distinguishing potential endocrine-disrupting compounds in environmental samples. Environ Sci Technol. 2004, 38 (23): 6396-6406. 10.1021/es049235r.

    Article  CAS  PubMed  Google Scholar 

  58. Chaudhary J, Sadler-Riggleman I, Ague JM, Skinner MK: The helix-loop-helix inhibitor of differentiation (ID) proteins induce post-mitotic terminally differentiated Sertoli cells to re-enter the cell cycle and proliferate. Biol Reprod. 2005, 72 (5): 1205-1217. 10.1095/biolreprod.104.035717.

    Article  CAS  PubMed  Google Scholar 

  59. Dong YL, Vegiraju S, Chauhan M, Yallampalli C: Expression of calcitonin gene-related peptide receptor components, calcitonin receptor-like receptor and receptor activity modifying protein 1, in the rat placenta during pregnancy and their cellular localization. Mol Hum Reprod. 2003, 9 (8): 481-490. 10.1093/molehr/gag058.

    Article  CAS  PubMed  Google Scholar 

  60. Zhang L, Rees MC, Bicknell R: The isolation and long-term culture of normal human endometrial epithelium and stroma. Expression of mRNAs for angiogenic polypeptides basally and on oestrogen and progesterone challenges. J Cell Sci. 1995, 108 (Pt 1): 323-331.

    CAS  PubMed  Google Scholar 

  61. Castro-Caldas M, Duarte CB, Carvalho AR, Lopes MC: 17beta-estradiol promotes the synthesis and the secretion of annexin I in the CCRF-CEM human cell line. Mediators Inflamm. 2001, 10 (5): 245-251. 10.1080/09629350120093713.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Huang EJ, Wu CC, Lee SD, Chen JH, Liu JY, Ko JL, Lin JA, Lu MC, Chen LM, Huang CY, et al: Opposing action of estrogen receptors alpha and beta on tumor necrosis factor-alpha gene expression and caspase-8-mediated apoptotic effects in HA22T cells. Mol Cell Biochem. 2006, 287 (1–2): 137-145. 10.1007/s11010-005-9092-4.

    Article  CAS  PubMed  Google Scholar 

  63. Stevens TA, Meech R: BARX2 and estrogen receptor-alpha (ESR1) co-ordinately regulate the production of alternatively spliced ESR1 isoforms and control breast cancer cell growth and invasion. Oncogene. 2006, 25 (39): 5426-5435. 10.1038/sj.onc.1209529.

    Article  CAS  PubMed  Google Scholar 

  64. Shah BH: Estrogen stimulation of COX-2-derived PGI2 confers atheroprotection. Trends Endocrinol Metab. 2005, 16 (5): 199-201. 10.1016/j.tem.2005.05.008.

    Article  CAS  PubMed  Google Scholar 

  65. Durbin JE, Hackenmiller R, Simon MC, Levy DE: Targeted disruption of the mouse Stat1 gene results in compromised innate immunity to viral disease. Cell. 1996, 84 (3): 443-450. 10.1016/S0092-8674(00)81289-1.

    Article  CAS  PubMed  Google Scholar 

  66. Kim HS, Lee MS: STAT1 as a key modulator of cell death. Cell Signal. 2007, 19 (3): 454-465. 10.1016/j.cellsig.2006.09.003.

    Article  CAS  PubMed  Google Scholar 

  67. Shang Y, Baumrucker CR, Green MH: The induction and activation of STAT1 by all-trans-retinoic acid are mediated by RAR beta signaling pathways in breast cancer cells. Oncogene. 1999, 18 (48): 6725-6732. 10.1038/sj.onc.1203084.

    Article  CAS  PubMed  Google Scholar 

  68. Zhou Q, Atadja P, Davidson NE: Histone deacetylase inhibitor LBH589 reactivates silenced estrogen receptor alpha (ER) gene expression without loss of DNA hypermethylation. Cancer Biol Ther. 2007, 6 (1): 64-69.

    Article  CAS  PubMed  Google Scholar 

  69. Blair HC, Robinson LJ, Zaidi M: Osteoclast signalling pathways. Biochem Biophys Res Commun. 2005, 328 (3): 728-738. 10.1016/j.bbrc.2004.11.077.

    Article  CAS  PubMed  Google Scholar 

  70. Boon LM, Mulliken JB, Vikkula M: RASA1: variable phenotype with capillary and arteriovenous malformations. Curr Opin Genet Dev. 2005, 15 (3): 265-269. 10.1016/j.gde.2005.03.004.

    Article  CAS  PubMed  Google Scholar 

  71. Hikima J, Lennard ML, Wilson MR, Miller NW, Clem LW, Warr GW: Evolution of vertebrate E protein transcription factors: comparative analysis of the E protein gene family in Takifugu rubripes and humans. Physiol Genomics. 2005, 21 (2): 144-151. 10.1152/physiolgenomics.00312.2004.

    Article  CAS  PubMed  Google Scholar 

  72. Hurteau GJ, Carlson JA, Spivack SD, Brock GJ: Overexpression of the microRNA hsa-miR-200c leads to reduced expression of transcription factor 8 and increased expression of E-cadherin. Cancer Res. 2007, 67 (17): 7972-7976. 10.1158/0008-5472.CAN-07-1058.

    Article  CAS  PubMed  Google Scholar 

  73. Drabsch Y, Hugo H, Zhang R, Dowhan DH, Miao YR, Gewirtz AM, Barry SC, Ramsay RG, Gonda TJ: Mechanism of and requirement for estrogen-regulated MYB expression in estrogen-receptor-positive breast cancer cells. Proc Natl Acad Sci USA. 2007, 104 (34): 13762-13767. 10.1073/pnas.0700104104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Urano A, Endoh M, Wada T, Morikawa Y, Itoh M, Kataoka Y, Taki T, Akazawa H, Nakajima H, Komuro I, et al: Infertility with defective spermiogenesis in mice lacking AF5q31, the target of chromosomal translocation in human infant leukemia. Mol Cell Biol. 2005, 25 (15): 6834-6845. 10.1128/MCB.25.15.6834-6845.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  75. Hogarth CA, Calanni S, Jans DA, Loveland KL: Importin alpha mRNAs have distinct expression profiles during spermatogenesis. Dev Dyn. 2006, 235 (1): 253-262. 10.1002/dvdy.20569.

    Article  CAS  PubMed  Google Scholar 

  76. Yoon SO, Shin S, Liu Y, Ballif BA, Woo MS, Gygi SP, Blenis J: Ran-binding protein 3 phosphorylation links the Ras and PI3-kinase pathways to nucleocytoplasmic transport. Mol Cell. 2008, 29 (3): 362-375. 10.1016/j.molcel.2007.12.024.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  77. Fang WB, Brantley-Sieders DM, Hwang Y, Ham A, Chen J: Identification and functional analysis of phosphorylated tyrosine residues within EphA2 receptor tyrosine kinase. J Biol Chem. 2008, 283 (23): 16017-16026. 10.1074/jbc.M709934200.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  78. Hemming JM, Waller WT, Chow MC, Denslow ND, Venables B: Assessment of the estrogenicity and toxicity of a domestic wastewater effluent flowing through a constructed wetland system using biomarkers in male fathead minnows (Pimephales promelas rafinesque, 1820). Environ Toxicol Chem. 2001, 20 (10): 2268-2275. 10.1897/1551-5028(2001)020<2268:AOTEAT>2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  79. Jensen KM, Korte JJ, Kahl MD, Pasha MS, Ankley GT: Aspects of basic reproductive biology and endocrinology in the fathead minnow (Pimephales promelas). Comp Biochem Physiol C Toxicol Pharmacol. 2001, 128 (1): 127-141. 10.1016/S1532-0456(00)00185-X.

    Article  CAS  PubMed  Google Scholar 

  80. Sokal R, Rohlf F: Biometry. 1995, New York: Freeman and Company

    Google Scholar 

  81. Garcia-Reyero N, Barber DS, Gross TS, Johnson KG, Sepulveda MS, Szabo NJ, Denslow ND: Dietary exposure of largemouth bass to OCPs changes expression of genes important for reproduction. Aquat Toxicol. 2006, 78 (4): 358-369. 10.1016/j.aquatox.2006.05.003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  82. Zahurak M, Parmigiani G, Yu W, Scharpf RB, Berman D, Schaeffer E, Shabbeer S, Cope L: Pre-processing Agilent microarray data. BMC Bioinformatics. 2007, 8: 142-10.1186/1471-2105-8-142.

    Article  PubMed Central  PubMed  Google Scholar 

  83. Benjamini Y, Hochberg Y: Controlling the false discovery rate – a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995, 57 (1): 289-300.

    Google Scholar 

  84. Ankley GT, Jensen KM, Kahl MD, Makynen EA, Blake LS, Greene KJ, Johnson RD, Villeneuve DL: Ketoconazole in the fathead minnow (Pimephales promelas): reproductive toxicity and biological compensation. Environ Toxicol Chem. 2007, 26 (6): 1214-1223. 10.1897/06-428R.1.

    Article  CAS  PubMed  Google Scholar 

  85. Villeneuve DL, Blake LS, Brodin JD, Greene KJ, Knoebl I, Miracle AL, Martinovic D, Ankley GT: Transcription of key genes regulating gonadal steroidogenesis in control and ketoconazole- or vinclozolin-exposed fathead minnows. Toxicol Sci. 2007, 98 (2): 395-407. 10.1093/toxsci/kfm124.

    Article  CAS  PubMed  Google Scholar 

  86. Martinovic D, Blake LS, Durhan EJ, Greene KJ, Kahl MD, Jensen KM, Makynen EA, Villeneuve DL, Ankley GT: Reproductive toxicity of vinclozolin in the fathead minnow: confirming an anti-androgenic mode of action. Environ Toxicol Chem. 2008, 27 (2): 478-488. 10.1897/07-206R.1.

    Article  CAS  PubMed  Google Scholar 

  87. Villeneuve DL, Knoebl I, Kahl MD, Jensen KM, Hammermeister DE, Greene KJ, Blake LS, Ankley GT: Relationship between brain and ovary aromatase activity and isoform-specific aromatase mRNA expression in the fathead minnow (Pimephales promelas). Aquat Toxicol. 2006, 76 (3–4): 353-368. 10.1016/j.aquatox.2005.10.016.

    Article  CAS  PubMed  Google Scholar 

Download references


This work was supported by the Environmental Protection Agency STAR grant, (R831848) to ND, MS, EO and KW, and by a fellowship from the Spanish Ministry of Sciences and Technology (EX-2004-0986) co-funded by the European Union to NGR. We wish to thank Dr. Thomas Hutchinson, (formerly of AstraZeneca) for the generous gift of ZM189,154. The research described in this article does not necessarily reflect the views of the EPA and no official endorsement should be inferred. NDD holds equity in EcoArray, Inc., a company commercializing the microarray technology used in this study.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Nancy D Denslow.

Additional information

Authors' contributions

NGR, NDD, KHW, MSS and EFO, conceived of the study and helped draft the manuscript; NGR, DLV, EFO, and KJK, performed the experimental work; LL, performed bioinformatics analysis; EJP, GTA, revised the manuscript critically for important intellectual content. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Garcia-Reyero, N., Kroll, K.J., Liu, L. et al. Gene expression responses in male fathead minnows exposed to binary mixtures of an estrogen and antiestrogen. BMC Genomics 10, 308 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: