- Research article
- Open Access
Revealing genes associated with vitellogenesis in the liver of the zebrafish (Danio rerio) by transcriptome profiling
© Levi et al; licensee BioMed Central Ltd. 2009
- Received: 05 September 2008
- Accepted: 31 March 2009
- Published: 31 March 2009
In oviparous vertebrates, including fish, vitellogenesis consists of highly regulated pathways involving 17β-estradiol (E2). Previous studies focused on a relatively small number of hepatic expressed genes during vitellogenesis. This study aims to identify hepatic genes involved in vitellogenesis and regulated by E2, by using zebrafish microarray gene expression profiling, and to provide information on functional distinctive genes expressed in the liver of a vitellogenic female, using zebrafish as a model fish.
Genes associated with vitellogenesis were revealed by the following paired t-tests (SAM) comparisons: a) two-month old vitellogenic (Vit2) females were compared with non-vitellogenic (NV) females, showing 825 differentially expressed transcripts during early stages of vitellogenesis, b) four-month old vitellogenic (Vit4) females were compared with NV females, showing 1,046 differentially expressed transcripts during vitellogenesis and c) E2-treated males were compared with control males, showing 1,828 differentially expressed transcripts regulated by E2. A Venn diagram revealed 822 common transcripts in the three groups, indicating that these transcripts were involved in vitellogenesis and putatively regulated by E2. In addition, 431 transcripts were differentially expressed in Vit2 and Vit4 females but not in E2-treated males, indicating that they were putatively not up-regulated by E2. Correspondence analysis showed high similarity in expression profiles of Vit2 with Vit4 and of NV females with control males. The E2-treated males differed from the other groups. The repertoire of genes putatively regulated by E2 in vitellogenic females included genes associated with protein synthesis and reproduction. Genes associated with the immune system processes and biological adhesion, were among the genes that were putatively not regulated by E2. E2-treated males expressed a large array of transcripts that were not associated with vitellogenesis.
The study revealed several genes that were not reported before as being regulated by E2. Also, the hepatic expression of several genes was reported here for the first time.
Gene expression profiling of liver samples revealed 1,046 differentially expressed transcripts during vitellogenesis of which at least ~64% were regulated by E2. The results raise the question on the regulation pattern and temporal pleiotropic expression of hepatic genes in vitellogenic females.
- Control Male
- Immune System Process
- Vitellogenic Stage
- Vitellogenic Female
- Yolk Precursor Protein
The accumulation of yolk in oocytes of oviparous animals during oocyte development is essential for proper embryonic development after fertilization and is therefore, a key process in successful reproduction. In fish, the egg yolk protein precursors (vitellogenins) are synthesized in the liver, secreted to the plasma and transported to the oocytes for uptake in a process known as vitellogenesis. Several metabolic changes occur during vitellogenesis in the maturing female fish as reflected in the pronounced increase in liver weight, RNA content, lipid deposition, glycogen depletion, plasma proteins, calcium and magnesium and phosphoprotein content [1, 3].
The most dominant trigger of vitellogenin (vtg) expression is the ovarian steroid hormone 17β-estradiol (E2) that is synthesized under the regulation of the hypothalamic–pituitary–gonad axis [reviewed in ]. Most data to-date supports the premise that the action of estrogens is mediated principally through specific nuclear Estrogen receptors (ERs). In the "classical" or "genomic" mechanism of E2 action, estrogens diffuse into the cell and bind to ERs, which are located in the cytosol or the nucleus of target cells. After ligand binding, the ERs form homo- or hetero dimers that bind to specific palindromic estrogen response elements (ERE) sequences  in the promoter region of estrogen-responsive genes, resulting in recruitment of coactivators or corepressors to the promoter. Subsequently this leads to increased or decreased mRNA levels and associated protein synthesis, resulting in the physiological response . Two main ERs (ERa and ERb) were characterized in mammals, birds and fish. Three ER subtypes were described so far for fish and include the Estrogen receptor 1, Estrogen receptor 2b and Estrogen receptor 2a [with the gene names of estrogen receptor 1 (esr1), estrogen receptor 2b (esr2b) and estrogen receptor 2a (esr2a), respectively) [7, 9]. Some of the effects of estrogens are so rapid that they cannot depend on RNA and protein synthesis and are known as non-genomic actions. They involve activating protein-kinase cascades, leading eventually to regulation of gene expression through phosphorylation and activation of transcription factors (TFs) within the nucleus [10, 12].
Hepatic expression of vtg is tightly coupled to E2-dependent up-regulation of esr1 expression [13, 15]. Vtg is specific to maturing females and therefore assessment of vtg expression or Vtg plasma levels is considered a useful approach in evaluating female maturity related with peripheral gonadal steroid changes . This protein is normally not detected in males or juveniles, but yolk precursor proteins can be detected in males or juveniles exposed to estrogens. Hepatocytes synthesize yolk precursor proteins when stimulated with exogenous estrogens or substances that mimic estrogens. Several changes in hepatic morphology such as proliferation of the endoplasmic reticulum and the Golgi apparatus also accompany estrogen stimulation. These aspects were investigated in several oviparous species [4, 17–23].
Vtg is a large (MW; 250–600 kDa) and complex calcium-binding phospholipoglycoprotein and in order to reach the end product found in the plasma, substantial post-translational modification must occur within liver cells. First, the protein backbone of Vtg is synthesized on membrane-bound ribosomes and subsequently the Vtg molecule is lipidated, glycosylated and phosphorylated [Reviewed in ]. In addition, Vtg may carry additional compounds such as retinal that are also transported to the developing oocytes . The genes involved in these processes have not been fully elucidated. In zebrafish, seven vtg genes were previously identified  but recent proteome profiling data from maturing ovarian follicles indicates the occurrence of eight vtg' s . The proteins fall into three main families represented by Vitellogenin 1 or VtgAo1 (with five corresponding genes, vtg 1, 4, 5, 6 and 7), Vitellogenin 2 or VtgAo2 (with two vtg2 genes) and Vitellogenin 3 or Vtg C (encoded by vtg3). Many more genes appear in the genome of zebrafish and fourteen of these genes were tightly linked to chromosome 22, while the phosvitinless gene (vtg3) was located on chromosome 11 .
In recent years, the study of hepatic expressed genes involved in fish vitellogenesis focused on few genes such as vtg, esr1, insulin-like growth factor 1 (igf1), zona pellucida glycoproteins (zp's), choriogenin H, cytochrome p450, family 1, subfamily a (cyp1a; also known as cyp1a1) and peroxisome proliferator-activated receptors (ppar's) [3, 15, 28–32] that are known to be regulated by estrogen. It is also well known that teleost apolipoproteins such as Apolipoprotein A-I (Apoa1), Apolipoprotein A-II (Apoa2) and Apolipoprotein E (Apoe) are regulated by E2 and presumably contribute to changes in the lipoprotein classes during vitellogenesis in fish [33, 35]. A high-throughput expression genomics approach would provide complementary information to the single-gene approaches used so far. Large-scale microarrays, available for model fish species including zebrafish, provide the opportunity to simultaneously monitor the expression of thousands of genes in different physiological stages during vitellogenesis. This approach has already been used with success to elucidate the zebrafish embryonic transcriptome [36, 37], to understand the molecular pathways defining gender specificity in zebrafish [38, 39] and to explore hepatic gene expression after exposure to different estrogens [40, 44].
The present study aims to identify genes involved in vitellogenesis and putatively regulated by E2 in the liver of zebrafish as a model fish, by using zebrafish oligonucleotide microarrays. Comparison of the hepatic expression profiles of vitellogenic and non-vitellogenic females provides information on the genes associated with vitellogenesis. In order to reveal E2-regulated genes, E2 treatment was administered to males for 48 hr at levels detected in the plasma of vitellogenic female. Genes suggested to be regulated by E2 were revealed by comparing the gene expression profiles of E2-treated males with those of vitellogenic females. The results also provide comparative information on the hepatic transcriptome profiles of 2- and 4- month old vitellogenic and non-vitellogenic females, E2-treated and control males and of the resemblance in gene expression profiles of these five groups. Distinctive putative pathways for the liver of vitellogenic females were found by analyses of the molecular functions and biological processes of the different treatment groups.
Oocyte developmental stages in ovaries of vitellogenic and non-vitellogenic females
The association of E2 plasma levels with hepatic transcript levels of esr1 and vtg3 in the experimental groups
Transcriptome analysis revealing gene expression patterns during vitellogenesis
A list of 20 most regulated annotated genes during vitellogenesis and after E2 treatment of males.
Function according to http://www.ncbi.nlm.nih.gov/sites/entrez gene
Vitellogenin 3, phosvitinless
Egg yolk precursor, phosvitinless.
Similar to reticulon 1
Associates with the endoplasmic reticulum.
Nothepsin, cathepsin e
Eukaryotic aspartyl (acid) protease.
Similar to low density lipoprotein receptor
Plays a central role in cholesterol metabolism.
Moderately similar to spectrin repeat containing, nuclear envelope 1
Involved in cytoskeletal structure.
Similar to ectonucleoside triphosphate diphosphohydrolase 4
Cleaving nucleotide tri- and diphosphates.
Similar to follistatin-like 1
Play an important role in tissue specific regulation.
Similar to human Prefoldin subunit 6
Binds and stabilizes newly synthesized polypeptides.
Lag1 homolog, ceramide synthase 2
May play a role in the regulation of cell growth.
Ankyrin repeat domain 6
Ankyrin repeats mediate protein-protein interactions.
Similar to nitric oxide synthase interacting protein
Promotes translocation of eNOS from the plasma membrane to intracellular sites.
Similar to heterogeneous nuclear ribonucleoprotein h1
RNA recognition/binding motif.
Similar to family with sequence similarity 46 c
The function of this gene is unknown.
Egg yolk precursor.
Similar to human guanylate binding protein1
Specifically bind guanine nucleotides (GMP, GDP, and GTP).
Similar to human protein fam20c precursor?
Has a crucial role in normal bone development.
Similar to glutamate receptor, ionotropic, kainate 1
Ligand-gated ion channel.
Similar to lectin, mannose-binding, 1 precursor
Mannose-specific lectin, a member of a Mannose-specific lectin, a member of a the secretory pathway of animal cells.
Similar to rna polymerase ii associated protein 1
The function of this gene is unknown.
Estrogen receptor 1
A ligand-activated transcription factor.
Genes suggested to be regulated by E2 during vitellogenesis
Genes putatively not directly regulated by E2 during vitellogenesis
Paired t-test analysis (SAM) comparing Vit2 or Vit4 females with NV females revealed 263 or 374, respectively transcripts that were not regulated by E2 (Fig 6B; see Additional file 2). The proportion of the down-regulated transcripts was higher (60.0% and 56.4% for Vit2 and Vit4 females, respectively) than the up-regulated transcripts indicating a slightly reverse trend from the E2 regulated genes.
Comparison between E2-treated males, control males and females
Putative processes associated with vitellogenesis using Gene Ontology terms
The number of genes regulated or not-regulated by E2 in selected GO functions.
Not regulated by E2
Lipid metabolic process
Immune system process
Transcription factor activity
Validation of microarray results
Validation of the microarray results was performed by testing the relative expression of 16 genes (see Additional file 4) in the same RNA samples that were used for the chip hybridization, by real-time PCR. The 16 tested transcripts were: 1) genes highly expressed in vitellogenic females and known to be induced by E2 treatment [esr1, vtg1, vtg3, nots and cytochrome p450 2k6 (cyp2k6)], 2) genes that were highly expressed in males and known to be down-regulated by E2 (cyp1a1 and igf1), and 3) genes that showed significantly different expression levels in the tested groups [retinoic acid receptor alpha a (raraa), alcohol dehydrogenase 5 (adh5), alcohol dehydrogenase 8b (adh8b), igf1 and retinol dehydrogenase 10 (rdh10), retinol dehydrogenase 14 (rdh14), dehydrogenase/reductase (SDR family) member 10 (dhrs10), stearoyl-desaturase (sCd), fatty acid desaturase 2 (fads2) and steroidogenic acute regulatory protein (star)]. A very high correlation was found between the microarray and real-time PCR results, with regression coefficients (R2 values) ranging from 0.9102 to 0.9340 (see Additional file 5).
This study aims at identifying and characterizing genes associated with vitellogenesis and defining the role of E2 in their regulation by using five physiological groups: non-vitellogenic females, vitellogenic females (2- and 4-month old), E2-treated males and control males. The expression levels of esr1 and vtg3 corresponding with the E2 plasma levels of females and E2 treated males, confirmed the efficacy of the E2 treatment. They also depicted the physiological state of the vitellogenic females, indicating that Vit2 females were at an interim stage between NV and Vit4 females.
The results provided a novel insight into the number and scope of hepatically regulated genes during vitellogenesis and indicated that most genes were regulated during the early stages of the process as young vitellogenic females (Vit2) differed by only 33 transcripts from older females (Vit4). Moreover, only ~64% of the transcripts regulated during vitellogenesis were suggested to be also regulated by E2. The resemblance in the gene expression pattern between non-vitellogenic females and males, stresses the specific change in pattern taking place during vitellogenesis in females. This change cannot be simply attributed to E2 as E2 treatment of males at physiological concentrations, resulted in a 1.8 fold higher number of genes than those regulated during vitellogenesis. These results also emphasize that the wide effects of xenobiotics with estrogen activity [45, 47] are not confined to genes associated with oocyte development. The following discussion sections highlight selected specific genes and putative pathways that were regulated during vitellogenesis and E2 treatment of males. Some of these putative pathways, were previously shown to be regulated in E2- treated males .
A list of 20 most differentially expressed hepatic transcripts reveals genes regulated by E2 and novel hepatic transcripts
The list of the 20 most differentially expressed hepatic transcripts includes genes known to be regulated by E2, genes that were not recorded previously as regulated by E2 or genes that were not previously reported to be expressed in the liver. Eight of the 20 most differentially expressed genes (Table 1), were previously reported to be regulated by E2, including: vtg1 and vtg3 [25, 48–50], nots [42, 51], syne1 , fst1 , nosip , grik1  and esr1 [15, 49]. Transcripts that were not associated previously with E2 regulation include reticulon1 (rtn1), entpd4 and lman1. A few transcripts (syne1, fstl1 and the fam20c) were not reported previously in hepatic cells. Proteins associated with cytoskeleton formation (syne1 and ank) showed higher up-regulation in vitellogenic females than after E2 treatment of males, supporting a putative role in the growth of the liver  and of the dramatic increase in protein synthesis and secretion by the endoplasmic reticulum  during vitellogenesis.
Zebrafish display eight different variant gene sequences [25, 27] for vtg genes but the array used here included probes only for vtg3 and vtg1. Since vtg1 shows high sequence similarity with vtg4, vtg5, vtg6 and vtg7 (all coding for VtgAo1), transcript levels for vtg1 may also reflect the expression of these genes. The higher expression levels of vtg3 compared with vtg1, may suggest their differential regulation by E2 stems from differences in the estrogen response elements (ERE's) in the promoter regions .
Putative molecular functions and biological processes regulated during vitellogenesis and by E2 treatment
Numerous prominent and putatively regulated functions were revealed in this study to take place during vitellogenesis and include lipid metabolism and lipid binding, hormone and cholesterol binding, transcription factor activity, immune response, immune system processes and cell adhesion (Table 2 and see Additional file 3). The task of allocating processes and functional significance to genes that were putatively regulated by E2 and also to those that were putatively not regulated by E2, was faced with a general difficulty for zebrafish. Linking specific pathways or modes of function with genes by the general descriptors provided by Gene Ontology (GO) for zebrafish is problematic due to the incomplete annotation of the zebrafish genome and a deficiency in functional studies. Consequently, several GO terms rely on homology of putative functions described in higher vertebrate species.
Lipid metabolism associated with vitellogenesis and E2-treatment of males
Several genes with a role in lipid metabolism were reported to be regulated by E2 [57, 58]. Genes associated with lipid metabolic processes putatively regulated in vitellogenic females or putatively regulated by E2, were also identified in the present study. Transcripts indicating a change in plasma lipoproteins were identified here, supporting previous published results on higher levels of plasma lipoproteins during vitellogenesis in fish [34, 35]. Plasma lipoproteins associated with transport of lipids are mainly synthesized in the liver and intestine . The protein components of lipoproteins, the apolipoproteins, form distinct complexes and two gene clusters, one consisting of apoa1, apoa4 and apoc3 and the other of apoe (apoeb for zebrafish), apoc, apoc2 and apoa4. These apolipoproteins are known from mammals and some were also characterized in fish . In the current study, there were no significant changes in apoa1 gene expression, in contrast with previous studies reporting on the downregulation of apoa1 gene expression after E2 and EE2 treatment [34, 42, 43]. Interestingly, Apoa1 was also found to serve as an antimicrobial protein and to be associated with the immune response system in fish . However, the expression of apoa4 and apoeb was down-regulated in vitellogenic females in comparison to NV females or after E2 treatment of males (see Additional file 1 and Additional file 3).
Several genes in the lipid metabolic processes are associated with the PPARs signaling pathway [62, 63] and some of them were found to be regulated in this study. The Ppars are ligand-inducible TFs belonging to the nuclear hormone receptor superfamily and are important regulators of lipid and energy homeostasis. Three isotypes have been identified in mammals, birds and amphibians, termed Ppar alpha (Ppara), Ppar beta (Pparb) (also known as Ppar delta) and Ppar gamma (Pparg) and each isotype is a product of one gene and shows distinct tissue distribution . Ppara functions in regulating reversible induction of β-oxidation in specific tissues but mainly the liver . Unexpectedly, there was no significant change in the expression of ppara in the liver of vitellogenic females or E2-treated males. The Pparb presumably functions in global control of lipid homeostasis and cellular proliferation and differentiation in mammals, is expressed in the liver and is moderately activated by a range of unsaturated fatty acids. Multiple/isoforms genes were found for pparb in several fish species, with two genes for pparb in zebrafish [64, 67]. The pparb1 was up-regulated in Vit4 females and control males while pparb2 was down-regulated in Vit4 and Vit2 females and by E2 treatment. In mammals, Pparg is associated with fat accumulation, particularly in adipocytes and in lipid accumulation in macrophages . Here, the pparg coactivator related 1 (ppargcr1) was up-regulated in Vit4 females and after E2 treatment (see Additional file 1).
A large array of transcripts putatively associated with lipid metabolism, was found to be regulated during vitellogenesis and by E2 treatment (see Additional file 3). Among them were genes coding for enzymes associated with fatty acid elongation and metabolism, and bile acid metabolism. Also, several transcripts of putative members of the cytochrome 450 superfamily were found to be regulated during vitellogenesis. This group of monooxygenases catalyzes several reactions involved in xenobiotic and drugs metabolism, synthesis of cholesterol, steroids and lipids.
Biosynthesis and catabolism of steroids
One of the explanations for the higher transcript levels of star, hsd3b7 and hsd17b 3 in NV females, could be related to the variable levels of steroid hormones in the plasma of females during onset of vitellogenesis [77, 78].
Hormone binding and TFs
The process of vitellogenesis is known to be regulated by E2 through the induction of esr1, as mentioned before. A similar putative expression pattern for genes coding for two additional hormone binding receptors was suggested to take place during vitellogenesis and regulated by E2 (Table 2, see Additional file 3). These include the progesterone receptor membrane component 2 (pgmrc) and the kappa opioid receptor 1 (oprk1) as shown by the regulation of their respective ESTs; an EST similar to pgmrc2 (AW153364) and an EST similar to oprk1 (BG883146). Transcripts of pgmrc genes were found in livers of mammals and oviparous vertebrates including zebrafish but their functions remain unknown [79, 80]. The gene pgmrc2 is known to be involved in progestin signaling in several vertebrate reproductive tissues and in the brain (reviewed in ) and Pgmrc1 (closely related to Pgmrc2) was suggested to have a role in the regulation of oocyte maturation in trout ovarian follicles . Since PGMRC1 was suggested to be regulated by testosterone in porcine hepatocytes , the high transcript levels of pgmrc2 found in vitellogenic females could be linked with the elevated testosterone levels associated with higher E2 plasma levels [77, 78]. The opioid receptors have multiple effects on reproductive, endocrine and immune functions. Transcripts of oprk were found to be widely expressed in rat tissues including in the liver, yet their function in the liver remained unknown .
Immune system processes and immune response
Surprisingly, all transcripts related to immune system processes in this study, were down-regulated in vitellogenic females and putatively not regulated by E2 (Fig 8). Three of these transcripts are also related to the immune response GO term (Table 2, see Additional file 3). The other transcripts related to the immune response term were all up-regulated by E2 treatment and some were also up-regulated during vitellogenesis (Table 2, see Additional file 3). This is in contrast to previous studies where E2 treatment repressed the expression of immune system and immune response related transcripts [92, 93]. Down-regulation in transcript levels in vitellogenic females (but not regulated by E2), was observed for genes coding for a class of intracellular molecules that play a role in coupling T-cell antigen receptor stimulation to the activation of integrins, major histocompatibility complex class genes (mhc's), proteases mediating programmed cell death or apoptosis [caspase 8 (casp8)] and a gene coding for an adherens junction protein [catenin beta 1 (ctnnb1)], functioning in communication and adhesion between cells, and anchoring the actin cytoskeleton. The genes showing high transcript levels in vitellogenic females and E2 treated males include an interferon activated gene, a TF that regulates mhc class II genes and the oprk1 previously described to be affected by steroid hormones, including E2 [94, 95].
Microarray profiling of liver samples revealed expression patterns characteristic of vitellogenic females of which only ~64% were found to be putatively regulated by E2. The repertoire of regulated genes implicates a wide range of functions especially those associated with protein synthesis, lipid metabolism, steroid biosynthesis, hormone binding and TF. Genes associated with the immune system and biological adhesion were among the genes that were up-regulated in vitellogenic females but not in E2-treated males, indicating that they were putatively not regulated by E2. E2-treated males expressed a large array of genes that were not associated with vitellogenesis. The study revealed several genes (rtn1, entpd4, lman1) that were not reported before as being regulated by E2. Also, the hepatic expression of several genes (e.g., syne1, fstl1, fam20c, star) is reported here for the first time for fish liver. In general, these results raise the question on the identity of the factors that regulate the pleiotropic expression of hepatic genes in vitellogenic females, in addition to E2.
Zebrafish were purchased from a local fish supplier (A&H Holdings, Israel LTD). All fish were maintained in 5-liter aquaria with UV treated, recycled and dechlorinated water and at ambient temperature of 25 ± 2°C with a light/dark cycle of 14/10 hr. The fish were fed twice a day, with shrimp nauplii (PGT, Eilat, Israel) in the morning and fish eggs in the afternoon. Non-vitellogenic females were kept under a light/dark cycle of 6/18 hr and fed twice a day with dry pellet food to avoid access to steroid compounds that maybe found in live food. For histological analysis the fish body or ovaries were fixed in Bouin's. Fixation, sectioning, and histological examination were performed according to . Paraffin sections of 4–7 µm were stained with hematoxylin and eosin. The terminology used by  for the zebrafish was adopted for this study. All fish were anaesthetized with Tricaine (Sigma-Aldrich, USA) before experimental procedures  and treatment of fish adhered with institutional regulations.
Samples collection and 17β-estradiol treatment
Two experiments were performed, one for evaluating the effect of E2 exposure in males and a second experiment for comparing gene expression in the liver of non-vitellogenic females with that of vitellogenic females. The experiments designs were as follows: 1) Four months old zebrafish were divided into three groups consisting of 8 fish in each group: i) adult spermeating males weighing 2.45 ± 0.207 g (N = 8) were exposed to E2 (Sigma-Aldrich, USA) by immersion for 48 hr (group E2-treated males). The concentration used was 5 μg/L (18 nM), as 3–4 ng/ml was determined to be the E2 natural concentration in the plasma of adult vitellogenic female ZF . A period of 48 hr of exposure was chosen as the highest expression level of esr1 was reached after 12 hr  and the expression of vtg stabilized after 48 hr of exposure to E2 and lasted for 17 days . The hepatosomatic index (HIS) of E2-treated males was 5.6 ± 0.6 after treatment. ii) control males weighing 2.26 ± 0.246 g (N = 8;) and a HIS of 3.2 ± 0.4; iii) vitellogenic females weighing 3.95 ± 0.303 g (N = 8) and showing a HIS of 5.9 ± 0.4 (group Vit4). Four replicate samples were prepared for each group and each replicate consisted of a pooled sample from livers of two fish. The samples were pooled after RNA extraction (see below). The fish in this experiment were kept with a light/dark cycle of 14/10 hr. The fish were fed twice a day, with shrimp nauplii (PGT, Eilat, Israel) in the morning and fish eggs in the afternoon. 2) In order to reveal the differences between vitellogenic and non-vitellogenic females, a second experiment was designed. One month old zebrafish were divided into two groups consisting of 32 fish in each group. Due to the small size of the liver, pooled samples from eight fish were prepared for each of the four replicates in the expression studies. The pooling of samples was done after total RNA extraction (see below). The groups consisted of: i) fish that were kept under a light/dark cycle of 14/10 hr for 5 weeks weighing 2.56 ± 0.426 g (N = 32) and a HIS of 5.6 ± 0.4, with ovaries in vitellogenic stage (group Vit2). Fish were kept at 25 ± 2°C and fed twice a day, with shrimp nauplii (PGT, Eilat, Israel) in the morning and with fish eggs in the afternoon. ii) fish that were kept under a light/dark cycle of 6/18 hr for the same time period and weighing 1.44 ± 0.391 g (N = 32) and showing a HIS of 4.6 ± 0.3, with non-vitellogenic ovaries (group NV). This group was fed only with dry pelleted food. In order to show the reversibility of this condition, some of the non-vitellogenic females were placed in tanks under the regular 10 hr light/14 hr dark photoperiod cycle. Fish were kept at 25 ± 2°C and were fed twice a day, with shrimp nauplii (PGT, Eilat, Israel) in the morning and with fish eggs in the afternoon. Ovaries collected after three weeks from these females were in the vitellogenic stage.
Determination of sex and of developmental stages of ovaries, were done by microscopic examination of the gonads. Oocyte stages were determined according to . The livers collected from all five groups were frozen instantly in liquid nitrogen and stored in -80°C until further use.
Measurements of E2 concentration in plasma was performed in similar designed separate experiments. Blood samples were collected using Micro-Hematocrit Tubes with Heparin (VWR, USA) from all tested groups and stored at -80°C. E2 concentrations were measured using Estradiol EIA Kit (Cayman, USA) according to the manufacture's protocol.
Total RNA was extracted from whole livers of zebrafish using Trizol reagent (GIBCO, USA) according to the manufacture's protocol followed by a clean-up and DNase treatement using RNeasy MiniElut Kit (Qiagene, Germany). After clean-up, 3 μl of the RNA samples were separated on 1.2% agarose gel to evaluate their quality and concentration. According to the gel picture, pooled samples from the livers of two fish were prepared for each group in Experiment 1 and of livers from eight fish for each group in the second experiment. Each RNA pool was quantified (A260) and assessed for purity (A260:A280 ratio) using Gene Quant (Amersham, UK) and by visual inspection of the 3 μg RNA separated on a denaturing gel.
Zebrafish Oligonucleotide Microarray
Zebrafish microarrays were prepared by the Kimmel Cancer Center at Thomas Jefferson University (TJU), Philadelphia, USA. The microarray is a single color system based on zebrafish oligonucleotide library from Compugen/Sigma Genosys and consists of 16,399 oligonucleotide probes (65 nt), representing 16,288 unique gene clusters. In order to minimize non-specific binding, CodeLink slides with a special coating were used. Also β-actin internal controls were used to monitor the labeling and hybridization quality.
Processed chips were scanned by using a Perkin Elmer ScanArray® XL 5000 Scanner, software version 3.1. Images were quantified by PerkinElmer (USA) Quant Array® Software 3.0. Quantization used the fixed circle method and outputs were total intensities. Microarray data were normalized across all arrays using quantile normalization of data in log base 2 scale . This method corrects background noise and non-specific hybridization.
Statistical analysis was performed using Significant Analysis of Microarray (SAM) software . For the multiclass analysis a false discovery rate (FDR) of 0.1% was used. For paired t-test comparisons between the different groups a FDR of 0.03% was used. Correspondence Analysis (COA) was performed using MultiExperiment Viewer (MeV) version 4.1. Clustering analysis was performed using Coupled Two Way Clustering (CTWC) algorithm [103, 105].
cDNA synthesis and Real-Time PCR
Real-time PCR was performed using the same RNA samples used for microarray hybridization (n = 16). For cDNA synthesis, 4 μg of total RNA were mixed with 0.1 µg of oligo-d(T) (Promega, USA), 4 µl of Bio-RT 5× buffer, 2 µM dNTP mix (Promega, USA), 200 U of Bio-RT (Bio-lab, Israel) and H2O to reach a final volume of 20µl. After an incubation of 1 h at 37°C, 80 µl of H2O were added to the reaction. The PCR mixture consisted of 0.5 ul of cDNA sample, 70 nM of each primer (see Additional file 4) and 12.5 μl of SYBR Green master mix (ABgene, UK), in a final volume of 25 μl.
Amplification was carried out in a GenAmp 5700 thermocycler (PE Applied Biosystems, USA) and according to the manufacture's protocol. Amplification was performed in triplicates and the results were analyzed with REST-384 version 2 . The relative expression of the 16 tested genes (see Additional file 4) was calculated using zebrafish elongation factor 1 alpha (ef1a) as a reference gene. The gene ef1a was recently shown to be a suitable reference gene for tissue analysis, developmental and E2 exposure studies of zebrafish [41, 107].
The microarray annotations were updated using BlastX program against nr database of the GeneBank. Blast2GO software  was used for achieving GO annotations for the 2,523 differentially expressed genes found by SAM analysis (see Additional file 1 and Additional file 3). Combinations of gene lists were performed using Gene List Venn Diagrams software  (Fig. 5).
This study was supported by the Israel Science Foundation Grant 1184/04. We would like to thank Ms. Hana Bernard for assistance in the graphics. We would also like to thank the reviewers for useful comments.
- Wiegand M: Composition, accumulation and utilization of yolk lipids in teleost fish. Rev Fish Biol Fish. 1996, 6: 259-286.Google Scholar
- Bjornsson BT, Haux C, Forlin L, Deftos LJ: The involvement of calcitonin in the reproductive physiology of the rainbow trout. J Endocrinol. 1986, 108: 17-23.PubMedGoogle Scholar
- Arukwe A, Goksoyr A: Eggshell and egg yolk proteins in fish: hepatic proteins for the next generation: oogenetic, population, and evolutionary implications of endocrine disruption. Comp Hepatol. 2003, 2: 4-PubMed CentralPubMedGoogle Scholar
- Polzonetti-Magni AM, Mosconi G, Soverchia L, Kikuyama S, Carnevali O: Multihormonal control of vitellogenesis in lower vertebrates. Int Rev Cytol. 2004, 239: 1-46.PubMedGoogle Scholar
- Gruber CJ, Gruber DM, Gruber IM, Wieser F, Huber JC: Anatomy of the estrogen response element. Trends Endocrinol Metab. 2004, 15: 73-78.PubMedGoogle Scholar
- Klinge CM, Jernigan SC, Mattingly KA, Risinger KE, Zhang J: Estrogen response element-dependent regulation of transcriptional activation of estrogen receptors alpha and beta by coactivators and corepressors. J Mol Endocrinol. 2004, 33: 387-410.PubMedGoogle Scholar
- 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: 10751-10756.PubMed CentralPubMedGoogle Scholar
- Ma CH, Dong KW, Yu KL: cDNA cloning and expression of a novel estrogen receptor beta-subtype in goldfish (Carassius auratus). Biochim Biophys Acta. 2000, 1490: 145-152.PubMedGoogle Scholar
- Menuet A, Pellegrini E, Anglade I, Blaise O, Laudet V, Kah O, Pakdel F: Molecular characterization of three estrogen receptor forms in zebrafish: binding characteristics, transactivation properties, and tissue distributions. Biol Reprod. 2002, 66: 1881-1892.PubMedGoogle Scholar
- Bjornstrom L, Sjoberg M: Mechanisms of estrogen receptor signaling: convergence of genomic and nongenomic actions on target genes. Mol Endocrinol. 2005, 19: 833-842.PubMedGoogle Scholar
- Hall JM, Couse JF, Korach KS: The multifaceted mechanisms of estradiol and estrogen receptor signaling. J Biol Chem. 2001, 276: 36869-36872.PubMedGoogle Scholar
- Klinge CM: Estrogen receptor interaction with estrogen response elements. Nucleic Acids Res. 2001, 29: 2905-2919.PubMed CentralPubMedGoogle Scholar
- Pakdel F, Feon S, Le Gac F, Le Menn F, Valotaire Y: In vivo estrogen induction of hepatic estrogen receptor mRNA and correlation with vitellogenin mRNA in rainbow trout. Mol Cell Endocrinol. 1991, 75: 205-212.PubMedGoogle Scholar
- Flouriot G, Pakdel F, Valotaire Y: Transcriptional and post-transcriptional regulation of rainbow trout estrogen receptor and vitellogenin gene expression. Mol Cell Endocrinol. 1996, 124: 173-183.PubMedGoogle Scholar
- Menuet A, Le Page Y, Torres O, Kern L, Kah O, Pakdel F: Analysis of the estrogen regulation of the zebrafish estrogen receptor (ER) reveals distinct effects of ERalpha, ERbeta1 and ERbeta2. J Mol Endocrinol. 2004, 32: 975-986.PubMedGoogle Scholar
- Heppell SA, Denslow ND, Folmar LC, Sullivan CV: Universal assay of vitellogenin as a biomarker for environmental estrogens. Environ Health Perspect. 1995, 103 (Suppl 7): 9-15.PubMed CentralPubMedGoogle Scholar
- Carnevali O, Mosconi G: In vitro induction of vitellogenin synthesis in Rana esculenta: role of the pituitary. Gen Comp Endocrinol. 1992, 86: 352-358.PubMedGoogle Scholar
- Carnevali O, Mosconi G, Yamamoto K, Kobayashi T, Kikuyama S, Polzonetti-Magni AM: Hormonal control of in vitro vitellogenin synthesis in Rana esculenta liver: effects of mammalian and amphibian growth hormone. Gen Comp Endocrinol. 1992, 88: 406-414.PubMedGoogle Scholar
- Carnevali O, Sabbieti MG, Mosconi G, Polzonetti-Magni AM: Multihormonal control of vitellogenin mRNA expression in the liver of frog, Rana esculenta. Mol Cell Endocrinol. 1995, 114: 19-25.PubMedGoogle Scholar
- Gobbetti A, Polzonetti-Magni A, Zerani M, Carnevali O, Botte V: Vitellogenin hormonal control in the green frog, Rana esculenta. Interplay between estradiol and pituitary hormones. Comp Biochem Physiol A. 1985, 82: 855-858.PubMedGoogle Scholar
- Polzonetti-Magni AM, Mosconi G, Carnevali O, Yamamoto K, Hanaoka Y, Kikuyama S: Gonadotropins and reproductive function in the anuran amphibian, Rana esculenta. Biol Reprod. 1998, 58: 88-93.PubMedGoogle Scholar
- Mosconi G, Carnevali O, Franzoni MF, Cottone E, Lutz I, Kloas W, Yamamoto K, Kikuyama S, Polzonetti-Magni AM: Environmental estrogens and reproductive biology in amphibians. Gen Comp Endocrinol. 2002, 126: 125-129.PubMedGoogle Scholar
- Tyler CR, Sumpter JP, Bromage NR: In-vivo ovarian uptake and processing of vitellogenin in the Rainbow trout, Salmo Gairdneri. The Journal of Experimental Zoology. 1988, 246: 171-179.Google Scholar
- Irie T, Seki T: Retinoid composition and retinal localization in the eggs of teleost fishes. Comp Biochem Physiol B Biochem Mol Biol. 2002, 131: 209-219.PubMedGoogle Scholar
- Wang H, Tan JT, Emelyanov A, Korzh V, Gong Z: Hepatic and extrahepatic expression of vitellogenin genes in the zebrafish, Danio rerio. Gene. 2005, 356: 91-100.PubMedGoogle Scholar
- Ziv T, Gattegno T, Chapovetsky V, Wolf H, Barnea E, Lubzens E, Admon A: Comparative proteomics of the developing fish (zebrafish and gilthead seabream) oocytes. Comparative Biochemistry and Physiology. 2008, Part D: 12-35.Google Scholar
- Finn RN, Kristoffersen BA: Vertebrate vitellogenin gene duplication in relation to the "3R hypothesis": correlation to the pelagic egg and the oceanic radiation of teleosts. PLoS ONE. 2007, 2: e169-PubMed CentralPubMedGoogle Scholar
- Davis LK, Pierce AL, Hiramatsu N, Sullivan CV, Hirano T, Grau EG: Gender-specific expression of multiple estrogen receptors, growth hormone receptors, insulin-like growth factors and vitellogenins, and effects of 17 beta-estradiol in the male tilapia (Oreochromis mossambicus). Gen Comp Endocrinol. 2008, 156: 544-551.PubMedGoogle Scholar
- Modig C, Modesto T, Canario A, Cerda J, von Hofsten J, Olsson PE: Molecular characterization and expression pattern of zona pellucida proteins in gilthead seabream (Sparus aurata). Biol Reprod. 2006, 75: 717-725.PubMedGoogle Scholar
- Murata K, Sugiyama H, Yasumasu S, Iuchi I, Yasumasu I, Yamagami K: Cloning of cDNA and estrogen-induced hepatic gene expression for choriogenin H, a precursor protein of the fish egg envelope (chorion). Proc Natl Acad Sci USA. 1997, 94: 2050-2055.PubMed CentralPubMedGoogle Scholar
- Bemanian V, Male R, Goksoyr A: The aryl hydrocarbon receptor-mediated disruption of vitellogenin synthesis in the fish liver: Cross-talk between AHR- and ERalpha-signalling pathways. Comp Hepatol. 2004, 3: 2-PubMed CentralPubMedGoogle Scholar
- Ibabe A, Herrero A, Cajaraville MP: Modulation of peroxisome proliferator-activated receptors (PPARs) by PPAR(alpha)- and PPAR(gamma)-specific ligands and by 17beta-estradiol in isolated zebrafish hepatocytes. Toxicol In Vitro. 2005, 19: 725-735.PubMedGoogle Scholar
- Bon E, Barbe U, Nunez Rodriguez J, Cuisset B, Pelissero C, Sumpter JP, Le Menn F: Plasma vitellogenin levels during the annual reproductive cycle of the female rainbow trout (Oncorhynchus mykiss): establishment and validation of an ELISA. Comp Biochem Physiol B Biochem Mol Biol. 1997, 117: 75-84.PubMedGoogle Scholar
- Wallaert C, Babin PJ: Age-related, sex-related, and seasonal changes of plasma lipoprotein concentrations in trout. J Lipid Res. 1994, 35: 1619-1633.PubMedGoogle Scholar
- Jensen B, Taylor M: Lipid transport in female Fundulus heteroclitus during the reproductive season. Fish Physiology and Biochemistry. 2002, 25: 141-151.Google Scholar
- Mathavan S, Lee SG, Mak A, Miller LD, Murthy KR, Govindarajan KR, Tong Y, Wu YL, Lam SH, Yang H, et al: Transcriptome analysis of zebrafish embryogenesis using microarrays. PLoS Genet. 2005, 1: 260-276.PubMedGoogle Scholar
- Ouyang M, Garnett AT, Han TM, Hama K, Lee A, Deng Y, Lee N, Liu HY, Amacher SL, Farber SA, et al: A web based resource characterizing the zebrafish developmental profile of over 16, 000 transcripts. Gene Expr Patterns. 2008, 8: 171-180.PubMed CentralPubMedGoogle Scholar
- Santos EM, Workman VL, Paull GC, Filby AL, Van Look KJ, Kille P, Tyler CR: Molecular basis of sex and reproductive status in breeding zebrafish. Physiol Genomics. 2007, 30: 111-122.PubMedGoogle Scholar
- Sreenivasan R, Cai M, Bartfai R, Wang X, Christoffels A, Orban L: Transcriptomic analyses reveal novel genes with sexually dimorphic expression in the zebrafish gonad and brain. PLoS ONE. 2008, 3: e1791-PubMed CentralPubMedGoogle Scholar
- Yang L, Kemadjou JR, Zinsmeister C, Bauer M, Legradi J, Muller F, Pankratz M, Jakel J, Strahle U: Transcriptional profiling reveals barcode-like toxicogenomic responses in the zebrafish embryo. Genome Biol. 2007, 8: R227-PubMed CentralPubMedGoogle Scholar
- Kausch U, Alberti M, Haindl S, Budczies J, Hock B: Biomarkers for exposure to estrogenic compounds: gene expression analysis in zebrafish (Danio rerio). Environ Toxicol. 2008, 23: 15-24.PubMedGoogle Scholar
- 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. Aquat Toxicol. 2007, 84: 38-49.PubMedGoogle Scholar
- Hoffmann JL, Torontali SP, Thomason RG, Lee DM, Brill JL, Price BB, Carr GJ, Versteeg DJ: Hepatic gene expression profiling using Genechips in zebrafish exposed to 17alpha-ethynylestradiol. Aquat Toxicol. 2006, 79: 233-246.PubMedGoogle Scholar
- Lam SH, Mathavan S, Tong Y, Li H, Karuturi RK, Wu Y, Vega VB, Liu ET, Gong Z: Zebrafish whole-adult-organism chemogenomics for large-scale predictive and discovery chemical biology. PLoS Genet. 2008, 4: e1000121-PubMed CentralPubMedGoogle Scholar
- Arukwe A, Celius T, Walther BT, Goksoyr A: Effects of xenoestrogen treatment on zona radiata protein and vitellogenin expression in Atlantic salmon (Salmo salar). Aquat Toxicol. 2000, 49: 159-170.PubMedGoogle Scholar
- Navas JM, Segner H: Antiestrogenicity of beta-naphthoflavone and PAHs in cultured rainbow trout hepatocytes: evidence for a role of the arylhydrocarbon receptor. Aquat Toxicol. 2000, 51: 79-92.PubMedGoogle Scholar
- Mortensen AS, Arukwe A: Targeted salmon gene array (SalArray): a toxicogenomic tool for gene expression profiling of interactions between estrogen and aryl hydrocarbon receptor signalling pathways. Chem Res Toxicol. 2007, 20: 474-488.PubMedGoogle Scholar
- Miracle A, Ankley G, Lattier D: Expression of two vitellogenin genes (vg1 and vg3) in fathead minnow (Pimephales promelas) liver in response to exposure to steroidal estrogens and androgens. Ecotoxicol Environ Saf. 2006, 63: 337-342.PubMedGoogle Scholar
- Islinger M, Willimski D, Volkl A, Braunbeck T: Effects of 17a-ethinylestradiol on the expression of three estrogen-responsive genes and cellular ultrastructure of liver and testes in male zebrafish. Aquat Toxicol. 2003, 62: 85-103.PubMedGoogle Scholar
- Tong Y, Shan T, Poh YK, Yan T, Wang H, Lam SH, Gong Z: Molecular cloning of zebrafish and medaka vitellogenin genes and comparison of their expression in response to 17beta-estradiol. Gene. 2004, 328: 25-36.PubMedGoogle Scholar
- Riggio M, Scudiero R, Filosa S, Parisi E: Sex- and tissue-specific expression of aspartic proteinases in Danio rerio (zebrafish). Gene. 2000, 260: 67-75.PubMedGoogle Scholar
- Rau SW, Dubal DB, Bottner M, Gerhold LM, Wise PM: Estradiol attenuates programmed cell death after stroke-like injury. J Neurosci. 2003, 23: 11420-11426.PubMedGoogle Scholar
- Cheng GF, Yuen CW, Ge W: Evidence for the existence of a local activin follistatin negative feedback loop in the goldfish pituitary and its regulation by activin and gonadal steroids. J Endocrinol. 2007, 195: 373-384.PubMedGoogle Scholar
- Harada H, Bharwani S, Pavlick KP, Korach KS, Grisham MB: Estrogen receptor-alpha, sexual dimorphism and reduced-size liver ischemia and reperfusion injury in mice. Pediatr Res. 2004, 55: 450-456.PubMedGoogle Scholar
- Matagne V, Lebrethon MC, Gerard A, Bourguignon JP: Kainate/estrogen receptor involvement in rapid estradiol effects in vitro and intracellular signaling pathways. Endocrinology. 2005, 146: 2313-2323.PubMedGoogle Scholar
- Olsson PE, Zafarullah M, Gedamu L: A role of metallothionein in zinc regulation after oestradiol induction of vitellogenin synthesis in rainbow trout, Salmo gairdneri. Biochem J. 1989, 257: 555-559.PubMed CentralPubMedGoogle Scholar
- Moens LN, Ven van der K, Van Remortel P, Del-Favero J, De Coen WM: Gene expression analysis of estrogenic compounds in the liver of common carp (Cyprinus carpio) using a custom cDNA microarray. J Biochem Mol Toxicol. 2007, 21: 299-311.PubMedGoogle Scholar
- Moens LN, Ven van der K, Van Remortel P, Del-Favero J, De Coen WM: Expression profiling of endocrine-disrupting compounds using a customized Cyprinus carpio cDNA microarray. Toxicol Sci. 2006, 93: 298-310.PubMedGoogle Scholar
- Babin PJ, Vernier JM: Plasma lipoproteins in fish. J Lipid Res. 1989, 30: 467-489.PubMedGoogle Scholar
- Kondo H, Morinaga K, Misaki R, Nakaya M, Watabe S: Characterization of the pufferfish Takifugu rubripes apolipoprotein multigene family. Gene. 2005, 346: 257-266.PubMedGoogle Scholar
- Villarroel F, Bastias A, Casado A, Amthauer R, Concha MI: Apolipoprotein A-I, an antimicrobial protein in Oncorhynchus mykiss: evaluation of its expression in primary defence barriers and plasma levels in sick and healthy fish. Fish Shellfish Immunol. 2007, 23: 197-209.PubMedGoogle Scholar
- Keller JM, Collet P, Bianchi A, Huin C, Bouillaud-Kremarik P, Becuwe P, Schohn H, Domenjoud L, Dauca M: Implications of peroxisome proliferator-activated receptors (PPARS) in development, cell life status and disease. Int J Dev Biol. 2000, 44: 429-442.PubMedGoogle Scholar
- Michalik L, Auwerx J, Berger JP, Chatterjee VK, Glass CK, Gonzalez FJ, Grimaldi PA, Kadowaki T, Lazar MA, O'Rahilly S, et al: International Union of Pharmacology. LXI. Peroxisome proliferator-activated receptors. Pharmacol Rev. 2006, 58: 726-741.PubMedGoogle Scholar
- Leaver MJ, Ezaz MT, Fontagne S, Tocher DR, Boukouvala E, Krey G: Multiple peroxisome proliferator-activated receptor beta subtypes from Atlantic salmon (Salmo salar). J Mol Endocrinol. 2007, 38: 391-400.PubMedGoogle Scholar
- Leaver MJ, Boukouvala E, Antonopoulou E, Diez A, Favre-Krey L, Ezaz MT, Bautista JM, Tocher DR, Krey G: Three peroxisome proliferator-activated receptor isotypes from each of two species of marine fish. Endocrinology. 2005, 146: 3150-3162.PubMedGoogle Scholar
- Ibabe A, Grabenbauer M, Baumgart E, Fahimi HD, Cajaraville MP: Expression of peroxisome proliferator-activated receptors in zebrafish (Danio rerio). Histochem Cell Biol. 2002, 118: 231-239.PubMedGoogle Scholar
- Ibabe A, Grabenbauer M, Baumgart E, Volkl A, Fahimi HD, Cajaraville MP: Expression of peroxisome proliferator-activated receptors in the liver of gray mullet (Mugil cephalus). Acta Histochem. 2004, 106: 11-19.PubMedGoogle Scholar
- Strauss JF, Kallen CB, Christenson LK, Watari H, Devoto L, Arakane F, Kiriakidou M, Sugawara T: The steroidogenic acute regulatory protein (StAR): a window into the complexities of intracellular cholesterol trafficking. Recent Prog Horm Res. 1999, 54: 369-394.PubMedGoogle Scholar
- Christenson LK, Strauss JF: Steroidogenic acute regulatory protein: an update on its regulation and mechanism of action. Arch Med Res. 2001, 32: 576-576.PubMedGoogle Scholar
- Stocco DM: StAR protein and the regulation of steroid hormone biosynthesis. Annu Rev Physiol. 2001, 63: 193-213.PubMedGoogle Scholar
- Borg B: Androgens in teleost fishes. Camp Biochem Physiol. 1994, 109C: 219-245.Google Scholar
- Hall EA, Ren S, Hylemon PB, Rodriguez-Agudo D, Redford K, Marques D, Kang D, Gil G, Pandak WM: Detection of the steroidogenic acute regulatory protein, StAR, in human liver cells. Biochim Biophys Acta. 2005, 1733: 111-119.PubMedGoogle Scholar
- Blum JL, Nyagode BA, James MO, Denslow ND: Effects of the pesticide methoxychlor on gene expression in the liver and testes of the male largemouth bass (Micropterus salmoides). Aquat Toxicol. 2008, 86: 459-469.PubMed CentralPubMedGoogle Scholar
- Filby AL, Thorpe KL, Tyler CR: Multiple molecular effect pathways of an environmental oestrogen in fish. J Mol Endocrinol. 2006, 37: 121-134.PubMedGoogle Scholar
- Mindnich R, Haller F, Halbach F, Moeller G, Hrabe de Angelis M, Adamski J: Androgen metabolism via 17beta-hydroxysteroid dehydrogenase type 3 in mammalian and non-mammalian vertebrates: comparison of the human and the zebrafish enzyme. J Mol Endocrinol. 2005, 35: 305-316.PubMedGoogle Scholar
- Kusakabe M, Nakamura I, Young G: 11beta-hydroxysteroid dehydrogenase complementary deoxyribonucleic acid in rainbow trout: cloning, sites of expression, and seasonal changes in gonads. Endocrinology. 2003, 144: 2534-2545.PubMedGoogle Scholar
- Modesto T, Canario AV: Morphometric changes and sex steroid levels during the annual reproductive cycle of the Lusitanian toadfish, Halobatrachus didactylus. Gen Comp Endocrinol. 2003, 131: 220-231.PubMedGoogle Scholar
- Laidley CW, Thomas P: Changes in plasma sex steroid-binding protein levels associated with ovarian recrudescence in the spotted seatrout (Cynoscion nebulosus). Biol Reprod. 1997, 56: 931-937.PubMedGoogle Scholar
- Thomas P: Characteristics of membrane progestin receptor alpha (mPRalpha) and progesterone membrane receptor component 1 (PGMRC1) and their roles in mediating rapid progestin actions. Front Neuroendocrinol. 2008, 29: 292-312.PubMed CentralPubMedGoogle Scholar
- Mourot B, Nguyen T, Fostier A, Bobe J: Two unrelated putative membrane-bound progestin receptors, progesterone membrane receptor component 1 (PGMRC1) and membrane progestin receptor (mPR) beta, are expressed in the rainbow trout oocyte and exhibit similar ovarian expression patterns. Reprod Biol Endocrinol. 2006, 4: 6-PubMed CentralPubMedGoogle Scholar
- Meyer C, Schmieding K, Falkenstein E, Wehling M: Are high-affinity progesterone binding site(s) from porcine liver microsomes members of the sigma receptor family?. Eur J Pharmacol. 1998, 347: 293-299.PubMedGoogle Scholar
- Wittert G, Hope P, Pyle D: Tissue distribution of opioid receptor gene expression in the rat. Biochem Biophys Res Commun. 1996, 218: 877-881.PubMedGoogle Scholar
- Rowlands JC, Gustafsson JA: Aryl hydrocarbon receptor-mediated signal transduction. Crit Rev Toxicol. 1997, 27: 109-134.PubMedGoogle Scholar
- Matthews J, Gustafsson J: Estrogen receptor and aryl hydrocarbon receptor signaling pathways. Nuclear Receptor Signaling. 2006, 4: e016-PubMed CentralPubMedGoogle Scholar
- Klinge CM, Bodenner DL, Desai D, Niles RM, Traish AM: Binding of type II nuclear receptors and estrogen receptor to full and half-site estrogen response elements in vitro. Nucleic Acids Res. 1997, 25: 1903-1912.PubMed CentralPubMedGoogle Scholar
- Klinge CM, Kaur K, Swanson HI: The aryl hydrocarbon receptor interacts with estrogen receptor alpha and orphan receptors COUP-TFI and ERRalpha1. Arch Biochem Biophys. 2000, 373: 163-174.PubMedGoogle Scholar
- Anderson MJ, Olsen H, Matsumura F, Hinton DE: In vivo modulation of 17 beta-estradiol-induced vitellogenin synthesis and estrogen receptor in rainbow trout (Oncorhynchus mykiss) liver cells by beta-naphthoflavone. Toxicol Appl Pharmacol. 1996, 137: 210-218.PubMedGoogle Scholar
- Arukwe A, Nordbo B: Hepatic biotransformation responses in Atlantic salmon exposed to retinoic acids and 3, 3', 4, 4'-tetrachlorobiphenyl (PCB congener 77). Comp Biochem Physiol C Toxicol Pharmacol. 2008, 147: 470-482.PubMedGoogle Scholar
- Arukwe A, Nordtug T, Kortner TM, Mortensen AS, Brakstad OG: Modulation of steroidogenesis and xenobiotic biotransformation responses in zebrafish (Danio rerio) exposed to water-soluble fraction of crude oil. Environ Res. 2008, 107: 362-370.PubMedGoogle Scholar
- Mortensen AS, Arukwe A: Activation of estrogen receptor signaling by the dioxin-like aryl hydrocarbon receptor agonist, 3, 3', 4, 4', 5-pentachlorobiphenyl (PCB126) in salmon in vitro system. Toxicol Appl Pharmacol. 2008, 227: 313-324.PubMedGoogle Scholar
- Mortensen AS, Arukwe A: Estrogenic effect of dioxin-like aryl hydrocarbon receptor (AhR) agonist (PCB congener 126) in salmon hepatocytes. Mar Environ Res. 2008, 66: 119-120.PubMedGoogle Scholar
- Williams TD, Diab AM, George SG, Sabine V, Chipman JK: Gene expression responses of European flounder (Platichthys flesus) to 17-beta estradiol. Toxicol Lett. 2007, 168: 236-248.PubMedGoogle Scholar
- Tilton SC, Givan SA, Pereira CB, Bailey GS, Williams DE: Toxicogenomic profiling of the hepatic tumor promoters indole-3-carbinol, 17beta-estradiol and beta-naphthoflavone in rainbow trout. Toxicol Sci. 2006, 90: 61-72.PubMedGoogle Scholar
- Piva F, Limonta P, Dondi D, Pimpinelli F, Martini L, Maggi R: Effects of steroids on the brain opioid system. J Steroid Biochem Mol Biol. 1995, 53: 343-348.PubMedGoogle Scholar
- Maggi R, Ma ZQ, Pimpinelli F, Maggi A, Martini L: Decrease of the number of opioid receptors and of the responsiveness to morphineduring neuronal differentiation induced by 17beta-estradiol in estrogen receptor-transfected neuroblastoma cells (SK-ER3). Neuroendocrinology. 1999, 69: 54-62.PubMedGoogle Scholar
- Biran J, Ben-Dor S, Levavi-Sivan B: Molecular identification and functional characterization of the kisspeptin/kisspeptin receptor system in lower vertebrates. Biol Reprod. 2008, 79: 776-786.PubMedGoogle Scholar
- Selman K, Wallace RA, Sarka A, Qi X: Stages of Oocyte Development in the Zebrafish, Brachydanio rerio. Journal of Morphology. 1993, 218: 203-224.Google Scholar
- Westerfield M: The Zebrafish Book. 1995, Oregon: University of Oregon Press, 3Google Scholar
- Heiden TK, Carvan MJ, Hutz RJ: Inhibition of follicular development, vitellogenesis, and serum 17beta-estradiol concentrations in zebrafish following chronic, sublethal dietary exposure to 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin. Toxicol Sci. 2006, 90: 490-499.PubMedGoogle Scholar
- Cohen A, Shmoish M, Levi L, Cheruti U, Levavi-Sivan B, Lubzens E: Alterations in micro-ribonucleic acid expression profiles reveal a novel pathway for estrogen regulation. Endocrinology. 2008, 149: 1687-1696.PubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193.PubMedGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121.PubMed CentralPubMedGoogle Scholar
- Coupled Two Way Clustering (CTWC). [http://ctwc.weizmann.ac.il/]
- Blatt M, Wiseman S, Domany E: Superparamagnetic clustering of data. Phys Rev Lett. 1996, 76: 3251-3254.PubMedGoogle Scholar
- Getz G, Levine E, Domany E: Coupled two-way clustering analysis of gene microarray data. Proc Natl Acad Sci USA. 2000, 97: 12079-12084.PubMed CentralPubMedGoogle Scholar
- Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: e36-PubMed CentralPubMedGoogle Scholar
- Tang R, Dodd A, Lai D, McNabb WC, Love DR: Validation of zebrafish (Danio rerio) reference genes for quantitative real-time RT-PCR normalization. Acta Biochim Biophys Sin (Shanghai). 2007, 39: 384-390.Google Scholar
- Blast2GO. [http://www.blast2go.de/]
- Gene List Venn Diagrams (GeneVenn). [http://mcbc.usm.edu/genevenn/genevenn.htm]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.