- Research article
- Open Access
Changes in the gene expression profiles of the brains of male European eels (Anguilla anguilla) during sexual maturation
BMC Genomicsvolume 15, Article number: 799 (2014)
The vertebrate brain plays a critical role in the regulation of sexual maturation and reproduction by integrating environmental information with developmental and endocrine status. The European eel Anguilla anguilla is an important species in which to better understand the neuroendocrine factors that control reproduction because it is an endangered species, has a complex life cycle that includes two extreme long distance migrations with both freshwater and seawater stages and because it occupies a key position within the teleost phylogeny. At present, mature eels have never been caught in the wild and little is known about most aspects of reproduction in A. anguilla. The goal of this study was to identify genes that may be involved in sexual maturation in experimentally matured eels. For this, we used microarrays to compare the gene expression profiles of sexually mature to immature males.
Using a false discovery rate of 0.05, a total of 1,497 differentially expressed genes were identified. Of this set, 991 were expressed at higher levels in brains (forebrain and midbrain) of mature males while 506 were expressed at lower levels relative to brains of immature males. The set of up-regulated genes includes genes involved in neuroendocrine processes, cell-cell signaling, neurogenesis and development. Interestingly, while genes involved in immune system function were down-regulated in the brains of mature males, changes in the expression levels of several receptors and channels were observed suggesting that some rewiring is occurring in the brain at sexual maturity.
This study shows that the brains of eels undergo major changes at the molecular level at sexual maturity that may include re-organization at the cellular level. Here, we have defined a set of genes that help to understand the molecular mechanisms controlling reproduction in eels. Some of these genes have previously described functions while many others have roles that have yet to be characterized in a reproductive context. Since most of the genes examined here have orthologs in other vertebrates, the results of this study will contribute to the body of knowledge concerning reproduction in vertebrates as well as to an improved understanding of eel biology.
Animal reproduction depends on the precise timing and coordination of behavioral and physiological processes. For successful reproduction to occur, animals need to reach the appropriate developmental stage, locate and attract mates and respond appropriately to reproductively relevant signals. The accuracy of each aspect of reproduction is essential to ensure that genetic material is passed on to the next generation.
A major function of the brain is to establish the appropriate developmental and endocrine status for reproduction, and then to co-ordinate this with aspects of behavior and sexual maturation. In vertebrates, this coordination is established through the hypothalamus and pituitary gland . The pituitary gland secretes follicle-stimulating hormone (FSH) and luteinizing hormone (LH) which are gonadotropins that act on the gonads to induce sex steroid synthesis and promote gametogenesis. The pituitary gonadotropes are in the vicinity of neuron endings that release neuropeptides and neurotransmitters, such as gonadotropin releasing-hormone (GnRH) and dopamine, which stimulate and inhibit gonadotropin secretion respectively. Following gonadotropin release, sex steroids produced by gonad feedback to the brain and pituitary provide information about reproductive status [2, 3].
Brain and pituitary gene expression changes that are associated with reproduction have been observed in several teleosts. These include changes in the pattern of GnRH expression at maturation in masu salmon Oncorhynchus masou[4, 5], European sea bass Dicentrarchus labrax, red seabream Pagrus major, goldfish  and European eel Anguilla anguilla. Changes in the expression levels of LHβ and FSHβ gene subunits have also been observed in the pituitary of several species, including European sea bass and European eel [9, 10]. Additionally, changes in the levels of kisspeptin system genes have also been observed in several teleosts [11, 12], as have changes in aromatase expression and activity [13, 14]. Many of these gene expression changes occur in areas of the brain that are associated with the detection and processing of sensory information.
Due to its complex life cycle and extreme long distance migrations, A. anguilla is an intriguing model in which to investigate reproduction related changes in the brain. The European eel is believed to spawn in a single location in the Sargasso Sea . After spawning, larvae are transported by surface currents in the Gulf Stream to the shores of Europe and North Africa. When larvae reach the continental shelf, they undergo metamorphosis into glass eels that migrate into fresh, brackish and coastal waters as yellow eels. After a highly variable feeding period of 5–30 years, yellow eels metamorphose into silver eels that migrate from their widely dispersed foraging areas to the common spawning ground in the Sargasso Sea. This migration is approximately 5,000-6,000 km long. The changes that occur during metamorphosis into silver eels include an increase in eye size, a change in body coloration and a regression of the digestive system. Eels do not feed during the transatlantic adult migration and rely exclusively on fat content for energy. After reaching the spawning grounds in the Sargasso Sea, eels reproduce once then die shortly thereafter . A. anguilla is also an important species to better understand the molecular aspects of reproduction because it is an endangered species, is of interest for aquaculture purposes and occupies a key position within the teleost phylogeny.
While much progress has been made over the last century, many aspects of eel reproduction remain poorly understood, including the physiological and molecular mechanisms responsible for sexual maturation and the coordination of reproductive events . European eels are sexually immature when they start their migration across the Atlantic to the spawning grounds; it is not clear if they mature sexually during migration or after they arrive at Sargasso Sea. Maturation experiments through hormonal induction have shown dramatic changes at the physiological level (e.g. an increase gonad size, an increase in the production of sex steroids and other sex hormones, and the progression of gonadogenesis from the immature stage to the gamete producing stage) as well as changes in secondary sexual characteristics (e.g. darkened skin, increased ocular indices, and soft swollen abdomens). Final maturation in eels does not appear to be triggered by swimming alone and is likely influenced by exposure to environmental cues as the eels swim across the Atlantic . In addition, reproductive pheromones are believed to have a role in the regulation of final maturation and reproduction in eels [17, 18]. If sexual maturation in A. anguilla is triggered by environmental signals, then signal detection is expected to occur through sensory pathways and signal processing is expected to occur in the brain before the neuroendocrine response is activated.
At present, mature adult eels have never been caught in the wild and eels do not mature in captivity without treatment with gonadotropic hormones. The inability to reach sexual maturation in captivity is due to a lack in sufficient pituitary gonadotropin production . In eels, maturation is under double neuroendocrine control which occurs in some teleosts [2, 3] and which can be reversed by gonadotropic hormone treatments [9, 20, 21]. The dopaminergic inhibitory process leads to insufficient GnRH stimulation, and to the inhibition of pituitary gonadotropes FSH and LH, which exert direct and indirect control on gonad development by acting on steroid metabolism [3, 22].
The goal of the present gene profiling study was to identify changes in the brain that are involved in the coordination of reproductive events in the European eel, including those involved in neuroendocrine regulation of gametogenesis. More specifically, genes involved in sexual maturation and reproduction are expected to be expressed at higher levels in the brain of sexually mature males compared to immature males. Additionally, since the coordination of reproductive events is influenced by environmental cues that must be integrated with the physiological and developmental status of the animal, it is expected that changes at the molecular level occur in the areas associated with the processing of sensory signals (e.g. telencephalon, diencephalon and olfactory bulb).
Using a false discovery rate of 0.05, a total of 1,497 differentially expressed genes were identified. Of this set, 991 were expressed at higher levels in brains (forebrain and midbrain) of sexually mature males while 506 were expressed at lower levels relative to brains of sexually immature males. When searched against the RefSeq proteins using BLASTX , these transcripts retrieved 732 and 389 unique sequences, respectively.
Functional annotation using DAVID
To obtain a list of gene identifiers that could be used for the functional analysis of differentially expressed genes, the A. anguilla transcriptome assembly contigs that were used to design the microarray were used in a BLASTX search of SwissProt proteins. The result of this search revealed that 11,130 of the 13,523 transcripts had alignments with e-values of 0.01 or less in SwissProt. When the contigs without hits in SwissProt were searched against RefSeq proteins and then the nr database, alignments for additional 1,149 and 110 sequences were retrieved, respectively. Using this approach, sequence identifiers for a total of 12,389 contigs were retrieved. For the analyses with DAVID , only differentially expressed transcripts with BLAST alignment e-values of 0.001 or less were considered.
Functional annotation analysis with DAVID uncovered a number of significantly enriched terms within the list of up- and down-regulated genes (Additional file 1: Table S1). Figure 1 details the proportion of over-represented terms belonging to the following GO categories: Biological process, Cellular component and Molecular Function. For each of the three main GO categories, terms were assigned to subcategories.
Cluster analysis with DAVID identified eight significantly enriched clusters of up-regulated genes and five significantly enriched clusters of down-regulated genes (Additional file 2: Table S2 and Additional file 3: Table S3). For this analysis, only clusters with EASE scores greater than or equal to 1.3 were considered. Clusters of up-regulated genes include many that are involved in DNA binding, calcium singling and nervous system development and differentiation. Clusters of down-regulated genes include genes for RNA processing, immune system function and nucleotide binding.
Differential expression: keyword analysis
To identify specific sets of differentially expressed genes, the protein descriptions obtained from the BLASTX search of the eukaryotic proteins in RefSeq were used. In total, 12,269 out of 13,523 assembly contigs showed alignments with an e-value score of 0.01 or better in the RefSeq proteins. The complete list of differentially expressed genes can be found in Additional file 4: Table S1. Among this set, we identified 64 neuroendocrine genes that are differentially expressed in the brains of sexually mature male eels (Table 1).
The results of the DAVID analysis suggest that some re-organization of signaling pathways in the brain occurs at sexual maturation. Therefore, it seemed likely that there would also be changes in the expression levels for genes that are involved directly in cell-cell communication such as receptors and channels. We therefore searched for up and down regulated channels and G protein-coupled receptors (GPCRs). Seven GPCRs were differentially expressed. Four of these belong to the Glutamate family of GPCRs , one belongs to the Rhodopsin-like family and two belong to the Secretin family (Table 2). One of the Glutamate family receptors, an extracellular calcium sensing receptor-like gene, was expressed at higher levels in the brains of sexually mature males. The other three were expressed at lower levels in the brains of sexually mature males compared with sexually immature males. The receptor from the Rhodopsin-like GPCR family was the serotonin receptor (5-hydroxytryptamine receptor 1B) which was expressed at higher levels in the brains of sexually mature males than in immature males (Table 2). The GPCR 64-like gene and the GPCR 125-like gene belong to the Secretin family and were expressed at higher levels in the brains of sexually mature males. Additionally, there were eleven differentially expressed channels (Table 2) seven of which were expressed at higher levels in the brains of sexually mature males in comparison to the brains of immature males.
Quantitative polymerase chain reaction
Eight differentially expressed genes were chosen for quantitative polymerase chain (qPCR) using the same RNA samples that were used for the microarray (Figure 2). This list includes: galanin peptide (GALP), glutamate decarboxylase 1 (GAD1), corticotropin releasing factor binding protein (CRFBP), C-type natriuretic peptide (CNP), gamma-aminobutyric acid receptor subunit δ (GABRD), neuropeptide Y (NPY), mammalian ependymin-related protein (MERP) and the iroquois-class homeodomain protein 2 (IRX2). The TATA box binding protein (TBP) and elongation factor 1-alpha 1 (EEF1A1) were used as controls.
To verify that the fold changes observed in the microarray experiments correlated with the fold changes observed in the qPCR experiments, Spearman’s rank order correlation was used. The results of this test indicate a positive relationship exists between the fold changes observed in the microarray and qPCR experiments (n = 24, r = 0.511, p = 0.011). The qPCR expression levels between immature and mature samples were further evaluated using Student's t-tests. GALP, CRFPB, CNP, MERP and IRX2 were expressed at significantly higher levels in the brains of mature males compared with mature (P < 0.05) while the differences in expression of NPY, GAD1 and GABDR were below the statistical significance level.
This study shows that major changes in the gene expression profiles occur in the brains of male European eels at sexual maturity. This shift includes an increase in the expression levels of genes involved in cell signaling, including receptor and channels, as well as genes involved in development and differentiation. Several neuroendocrine genes were also expressed at higher levels in the brains of mature males in comparison to the brains of immature males. These genes are likely involved in gonadal maturation and development and act through the brain-pituitary-gonad axis. Genes involved in immune system function were down-regulated in the brains of sexually mature males, suggesting that a re-allocation of resources occurs at this stage. Interestingly, the expression levels of several receptors and channels change at sexual maturity suggesting that some neuronal rewiring occurs at this stage.
More transcripts for neuropeptide Y (NPY) were observed in the brains of mature male eels in comparison to immature males. NPY is a short peptide whose expression pattern has been described in the brains of several teleosts. For example, in goldfish, NPY expression has been detected in the forebrain including the telencephalon, preoptic area, the olfactory bulb and thalamic regions . In the Japanese eel Anguilla japonica, NPY is highly expressed in the telencephalon, optic tectum thalamus and hypothalamus . NPY stimulates LH release in fish , has a role in growth hormone (GH) and gonadotropin hormone (GTH) secretion [29–31] and may coordinate growth, feeding and reproduction in fish . In goldfish, fasting and exposure to sex steroids increases NPY mRNA levels in the brain [32, 33]. Furthermore, the NPY response in goldfish depends on reproductive status . Food deprivation in A. japonica also causes increases of NPY mRNA levels in the brain . Changes in feeding behavior and regulation of energy expenditure are expected to accompany the maturation process in fish, particularly in species such as A. anguilla that allocate a lot of resources to reproduction . Since the European eel does not feed at sexual maturity, it is possible that increased levels of NPY are related to nutritional and developmental status.
Genes involved in gamma-aminobutyric acid (GABA) signaling were also expressed at higher levels in the brains of mature males. For example, glutamate decarboxylase 1 (GAD1), an enzyme involved in the production of GABA from L-glutamic acid, was more abundant in the brains of sexually mature than immature males. GABA is a neurotransmitter involved in the control of pituitary hormone secretion in teleosts . GABA has a stimulatory effect on LH by enhancing GnRH release and inhibiting dopamine turnover [2, 37]. The GABA A receptor subunit δ (GABRD) was expressed at a higher level in the brains of sexually mature male eels than in immature males. GABA A receptors are hetero-pentamer ligand gated channels that are gated by GABA binding. It appears that different combinations of receptor subunits may confer functional specificity . In goldfish, four GABA A subunits (α 1, β 4, γ 1 and 2) were expressed at higher levels in the brains of sexually mature females . Similarly, the GABA A receptor γ 2 subunit was expressed at higher levels in the brains of goldfish exposed to prostaglandin-F2α (PGF2α) than in controls . In the current study, the δ subunit was up-regulated in sexually mature males; this is intriguing given that the δ subunit appears to be involved in neurosteroid modulation of the GABA A receptor .
In the brains of sexually mature males, the 5-hydroxytryptamine (serotonin) receptor 1B was expressed at a higher level than in immature males. The brain serotonergic system is one of the mechanisms used to cope with the dramatic changes during maturation. This system plays an important role in a complex neuroendocrine loop that functions in homeostasis and acclimation during physiological or environmental challenges [42, 43].
Corticotropin releasing factor binding protein (CRFBP) was also expressed at a higher level in the brains of mature males than in immature males. CRFBP has an inhibitory role on CRF which is a regulator of stress responses . In fish, acute stress results in an increase in CRFBP expression in the hypothalamus . For European eels, reproduction is a stressful event because the spawning migration journey is very long (5,000-6,000 km). In addition, it is presumably challenging for both sexes to coordinate the temporal and spatial aspects of spawning in the vast Sargasso Sea which is 1,000 km wide and 3,000 km long. At the mature stage, European eels are likely to reduce their threshold for assessing risks and threats since they have a single opportunity to reproduce. Therefore the observed increases in CRFBP may be part of the stress response and of changes in behavior that are associated with reproduction.
An increase in the levels of mammalian ependymin-related protein (MERP) in the brains of mature males was observed. MERPs belong to a phylogenetically diverse clade that is related to ependymins  which are glycoproteins that are found in abundance in the cerebrospinal fluid of fish. Ependymins are involved in synaptic changes during memory consolidation and regeneration of the optic nerve [47–49]. Changes in the levels of ependymin-II expression have been observed in the hypothalamus in goldfish at breeding season  and in the telencephalon following exposure to a female pheromone .
Several genes involved in calmodulin (CALM) signaling were differentially expressed in the brains of sexually mature males (Table 1). For example, CALM is involved in calcium signaling and nervous system function including memory consolidation  and was up-regulated in the brains of sexually mature male eels. In goldfish, CALM levels also change with reproductive status. During breeding season, CALM and calcium/calmodulin-dependant protein kinase 2 (CaMK2) levels are highest in the hypothalamus in female goldfish  and CALM levels are highest in the telencephalon of males following exposure to PGF2α. Sexual maturation alters the eel’s internal state and changes in the way the animal relates to its environment. In the European eel, the memory circuitry is expected to be important for reproduction due to the long migration across the Atlantic to the Sargasso Sea. Though CaM signaling genes are likely involved in many processes in the brains of mature eels, it is possible that the changes observed here are also related to memory.
The microarray results showed a high positive correlation with those obtained in selected genes by qPCR thus providing a global validation of the microarray results. NPY, GAD1 and GABDR did not seem to follow the same pattern; this could be related to the low expression levels of these genes (and the associated higher error rate of measurement) or to the existence of alternative transcripts that are differentially detected by the two methods. However, these possibilities were not further investigated.
Recent studies in fish show that specific cells in the brain respond to light and exert neuroendocrine control of testis development [51, 52]. In A. anguilla, behavioral evidence indicates that photoreceptor sites that are independent of the pineal exist in the brain . Eels become more sexually mature after months of migration in deep waters while only surfacing at night. In A. anguilla, spawning is believed to occur in the top 200 meters  where there is more light than in deeper water. The results of the current study show that several genes with roles in photosignal transduction and visual system development were up-regulated in the brains of mature males. This list includes retinal rod rhodopsin-sensitive cGMP 3' 5'-cyclic phosphodiesterase subunit γ and cyclic nucleotide-gated cation channel α-3 which is involved in olfactory and visual signal transduction. It also includes the calcium binding protein 4, which has an important role in synaptic function in photoreceptor cells  and the iroquois-class homeodomain protein (IRX5). IRX5 is a transcription factor that is involved in cone cell differentiation and retinal, craniofacial and gonad development [55, 56]. The transcription factor HES4 was also up-regulated in the brains of sexually mature males. HES4 belongs to a family of genes that have roles in cell differentiation and has a role in retina development . While it is possible that these genes have other roles in the brains of mature males, an intriguing avenue for future research would be to determine whether there is a link between neurohormone systems, gonad development and photoreceptors in the brains of A. anguilla like there is in masu salmon .
Several genes involved in neurogenesis, cellular differentiation and development were up-regulated in the brain of sexually mature males (Additional file 1: Table S1). This set includes genes from the GO categories such as GO:0050767 (regulation of neurogenesis), GO:0007411 (axon guidance), GO:0045664 (regulation of neuron differentiation) and GO:0007423 (sensory organ development). In addition, many differentially expressed ion channels and GPCRs were differentially expressed (Table 2) which suggests that some re-organization and/or re-growth occurs at maturation. The calcium-activated potassium channel would tend to hyperpolarize (inhibit) neurons in response to neuromodulators/neurotransmitters. Down-regulation of these channels suggests a less responsive, more ‘hard-wired’ brain, which may be advantageous for certain aspects of reproduction.
GnRH, aromatase, kisspeptin system genes and sex steroid receptors are also known to have roles in reproduction. The microarray included a probe for one isoform of GnRH1 and no differences in expression were observed between immature and mature males. The array did not contain probes for other GnRH genes or isoforms. Although GnRH1 controls gonadotrophin secretion in the pituitary , the fact that the fish were artificially induced to mature using human chorionic gonadotrophin suggests that no positive or negative feedback mechanisms towards GnRH1 were occurring at the time of sampling. The microarray also contained probes for transcripts that are similar to the estrogen-related receptor beta like 1 and the androgen receptor beta and no significant differences in expression were detected. Kisspeptin system genes and aromatase were not analyzed in this survey. At present, mature eels have never been caught in the wild and the only way to obtain sexually mature animals is through artificial induction of maturation. It is therefore not clear to what extent changes that occur at the molecular level in the brain are the result of specific gonadal feedback mechanisms.
The eyes and the olfactory epithelium change with sexual maturation in A. anguilla and the American eel Anguilla rostrata[59, 60] suggesting that there is a change in the use of sensory modalities with reproduction. As mentioned previously, the data from this study suggest that the forebrains of male eels undergo some sort of rewiring at sexual maturation. As the areas of the brain tested here are associated with the processing and conversion of sensory input to physiological responses, it is possible that changes in the brain at sexual maturity are influenced by changes in sensory system requirements that occur at reproduction.
In this study, we used nearly 15,000 microarray probes to identify gene expression changes in the brains of eels that occur at sexual maturation. Of the set of differentially expressed genes, 991 were expressed at higher levels in sexually mature males while 506 were expressed at lower levels than in immature males. This finding shows that the brains of eels undergo major changes at the molecular level during reproduction. The set of up-regulated genes includes genes that are involved in neuroendocrine processes which likely act through the brain-pituitary-gonad axis. It also includes genes involved in cell-cell signaling, neurogenesis and development indicating that the brain is undergoing developmental changes that are associated with maturation and reproduction. A major advantage to the approach employed here is that it allows many genes to be tested at the same time. It also allows for the identification of genes not previously known to be involved in particular processes such as reproduction. Consequently, some of these differentially expressed genes identified here have roles that have been described in other studies while others have roles that have yet to be characterized. Furthermore, since the majority of the genes tested here have orthologs in other vertebrates, the results from this study can be used to understand reproduction across a diversity of vertebrate species. Most importantly however, is that we have defined a set of genes that appear to be important for sexual maturation in A. anguilla. These findings therefore contribute to the body of knowledge concerning eel biology and will, in the future, contribute to aquaculture and conservation efforts.
Approximately four hundred sexually immature European eels (80–90 g/fish) were obtained from a local distributor (Pescafial, Seville, Spain) and kept for acclimation in a 2000 L tank with freshwater (0 PSU salinity) and under natural temperature and photoperiod. After 14 days, approximately 120 fish were transferred to two independent temperature controlled recirculation water units connected to two 500 L cylindro-conical tanks (30 fish per tank) where the fish were acclimated to seawater (35 PSU salinity) over two weeks. Recirculation units were kept at a 20-22°C at 12 L/12D photoperiod and the fish were fed with commercial pellets (DIBAQ No. 3, Proaqua, Spain) at 3% of body weight per day. Males in one of the seawater units were injected intramuscularly every week for 140 days with 2000 IU hCG/kg (human chorionic gonadotropin, Sigma–Aldrich Chemical) in 0.9% saline to induce sexual maturation. Males in the other recirculation unit were injected weekly over the same period with 0.9% NaCl (vehicle).
Finally, eels were anesthetized in a solution of 100 mg/L tricaine methanesulfonate (MS-222, Sigma Aldrich) then decapitated. Tissues (brain, including the olfactory bulbs, telencephalon, diencephalon and mesencephalon) were collected from sexually immature (n = 12) and sexually mature males (n = 12) and frozen immediately in liquid nitrogen until the tissue collections were complete. They were then placed in RNAlater® (Sigma R0901) and kept overnight at 4°C then transferred to -20°C. Each tissue sample was pooled with the tissue from two other samples. Therefore each replicate (n = 4 mature, n = 4 immature) contains tissue from three individuals. Sexual maturation was characterized as follows : seawater immature fish were designated stage III (presence of spermatogonia and spermatocytes), whereas sexually mature males were designated stage VI (presence of only spermatozoa and few spermatogonia which is called the functional maturation stage). The average gonadosomal index (GSI) values for sexually mature and immature males were 8.41 ± 0.56% and 0.084 ± 0.01% respectively. The mean eye diameters for mature males were 8.19 ± 0.16 mm and 6.51 ± 0.13 mm for immature males. Brain samples were shipped on dry ice. Fish care and experimentation complied with national legislation for the use of laboratory animals under a group-1 license issued by the Directorate-General for Veterinary, Ministry of Agriculture, Rural Development and Fisheries of Portugal.
RNA was isolated from the brain using the RNAeasy Mini Kit (Qiagen) following the manufacter’s protocol. RNA was extracted separately for olfactory bulbs, telencephalons and diencephalon/mesencephalons. Afterwards, RNA concentration was measured using a Nanodrop ND-100 spectrophotometer (NanoDrop Technologies) and RNA integrity and quality were estimated using an Agilent 2100 Bioanalyzer (Agilent Technologies). A minimum RNA integrity number (RIN) of 7.5 was considered. Because RNA concentrations of olfactory bulb, telencephalon and diencephalon were below the amount needed for microarray analysis, a separate experiment for each part of the brain was not possible. Therefore, we combined equal concentrations of olfactory bulb, telencephalon and diencephalon/mesencephalon as a single brain sample for each replicate. Each replicate therefore contains RNA from the olfactory bulb, telencephalon and diencephalon/mesencephalon from three individuals.
Microarray analysis was conducted using an European eel-specific array consisting of a total of 14,913 probes based on a large collection of high-throughput transcriptomic sequences . Probe sequences and further details on the microarray platform can be found on the GEO database under accession number GPL15124.
Sample labelling and hybridization were conducted following the details in . Hybridized slides were scanned at 5 μm resolution using an Agilent DNA microarray scanner. Slides were scanned at two different sensitivity levels (XDR Hi 100% and XDR Lo 10%) to increase the power to detect both lowly and highly expressed genes. The two linked images generated were analyzed together. Data were extracted and background subtracted using the standard procedure in Agilent Feature Extraction (FE) software v. 9.5.1. Hybridization success was evaluated using flag values, excluding those intensities not equal to 1. Data was normalized using a quantile normalization procedure using R (http://www.r-project.org). Normalized fluorescence data from the arrays have been deposited in the GEO database (http://www.ncbi.nlm.nih.gov/geo) under accession number GSE55858.
Differentially transcribed genes across samples were identified using the program SAM (Significance Analysis of Microarrays) version 4.0 . Groups (mature vs. immature males) were compared using the two-class unpaired test and up and down regulated genes were identified. Only genes with a minimum fold change of 1.5 were considered.
Functional annotation using DAVID
To obtain a list of gene identifiers that could be used for the functional analysis of differentially expressed genes, BLASTX (2.2.25+ suite)  was used to search SwissProt. For this analysis, the SwissProt database was chosen to maximize the number of sequence identifiers that DAVID version 6.7  would recognize. In this search, the European eel transcriptome assembly contigs (n = 13,523) that were used to design the array  were used as queries. Assembly contigs that did not have a BLASTX hit in SwissProt meeting the e-value cutoff of 0.01 were then searched against RefSeq proteins. Those without hits in SwissProt or RefSeq were then searched against the National Center for Biotechnology Information (NCBI) nr database using the same e-value cutoff.
For functional annotation analysis with DAVID, only BLAST alignments with an e-value less than or equal to 0.001 were considered. To identify differentially expressed genes, a false discovery rate (q-value) of 0.05 was used. For the DAVID analysis, the standard default settings that include an EASE score of 0.1 were used. For the functional annotation analysis, the annotation categories selected were: GO terms, Interpro, KEGG pathways and SP_PIR keywords. For the cluster analysis with DAVID, only clusters with EASE scores of greater than or equal to 1.3 were considered. To obtain a more descriptive view of up and down-regulated processes, GO terms belonging to the three GO categories (Biological Processes, Cellular Cycle and Molecular Function) were further divided into subcategories.
Differential expression: keyword analysis
To identify specific subsets of genes that are differentially expressed in the brains of sexually mature male eels, the set of annotations from a BLASTX search of the RefSeq eukaryotic proteins (n = 3,360,895) was used. This approach was chosen because RefSeq, like SwissProt, is a curated dataset yet contains greater number and more diverse set of sequences from fish than SwissProt. In this search, the European eel assembly contigs were used as queries and an e-value cutoff of 0.01. Protein descriptions were retrieved for each blast hit and the results were loaded into a MySQL table which was used for a keyword search to identify genes and sets of genes that are differentially expressed and are involved in neuroendocrine and reproductive processes.
Quantitative polymerase chain reaction
The expression levels for eight differentially expressed genes were determined by quantitative polymerase chain reaction (qPCR) with the same RNA pools that were used for the microarray. This set was chosen because it includes genes with small and large fold changes (fold changes of 1.63-9.78). It includes GALP, GAD1, CNP, CRFBP, GABRD, NPY, MERP and IRX2. TBP and EEF1A1 were used as controls. For qPCR, three replicates for each condition (immature and immature) were used instead of four because there was not enough material left in the other two samples. cDNA was synthesized from RNA in 20 μl reactions containing 500 ng of DNase-treated RNA, 200 ng of random hexamers (Jena Biosciences, Germany), 100 U of RevertAid (Fermentas, Thermo Fisher Scientific, USA) reverse transcriptase (RT) and 8 U of RiboLockRNase Inhibitor (Fermentas). Reactions were incubated for 10 minutes (min) at 25°C and 60 min at 42°C followed by enzyme inactivation for 10 min at 70°C. Quantifications were carried out using the relative standard curve method and the EvaGreen chemistry in duplicate 15 μl reactions containing: 2 μl of cDNA (diluted 1:4), 300 nM each specific primer (Additional file 5: Table S5) and 1x Sso Fast EvaGreen Supermix (Bio-Rad Laboratories, USA). Reactions were carried out using a Bio-Rad iClycler iQ5 qPCR thermocycler with the following cycling conditions: 30 seconds (s) at 95°C, 45 cycles of 5 s at 95°C, 10 s at 60°C and final melting curves between 65 and 98°C. Target specificity for each qPCR reaction was confirmed by verifying that the amplicon was the expected size by gel electrophoresis and by the presence of a single peak in each melting curve. A cDNA synthesis reaction without reverse transcriptase was used as a negative control to ensure all traces of genomic DNA had been removed and no amplification was observed for all primer pairs. Standard curves relating amplification cycle with initial template quantity were generated using serial dilutions of specific RT-PCR products for each gene (obtained from brain cDNA using the same primers) and included in all quantification plates, with qPCR efficiency ranging between 85–97% with R2 > 0.98. The geometric mean of reference genes TBP and EEF1A1 was used to normalize expression data. Relative expression levels were calculated by dividing the expression of each target gene by the geometric mean of the two reference genes.
To verify that the fold changes observed in the microarray experiments correlated with the fold changes observed in the qPCR experiments, Spearman’s rank order correlation was used. Fold changes for each biological replicate were calculated using the mean expression values from the immature brain samples. Student's t-tests were used to determine if expression levels measured by qPCR were significantly different in the brains of immature and mature males. The normality and equal variance assumptions were verified using Shapiro-Wilk’s test and equal variance tests.
Calcium/calmodulin-dependant protein kinase 2
C-type natriuretic peptide
Corticotropin releasing factor binding protein
Elongation factor 1-alpha 1
γ-aminobutyric acid receptor subunit δ
Glutamate decarboxylase 1
G protein-coupled receptors
Human chorionic gonadotropin
Transcription factor HES-4
Iroquois-class homeodomain protein 2
Iroquois-class homeodomain protein 5
Mammalian ependymin-related protein
Quantitative polymerase chain reaction
TATA binding protein.
Weltzien FA, Andersson E, Andersen Ø, Shalchian-Tabrizi K, Norberg B: The brain–pituitary–gonad axis in male teleosts, with special emphasis on flatfish (Pleuronectiformes). Comp Biochem Physiol A Mol Integr Physiol. 2004, 137 (3): 447-477.
Zohar Y, AntonioMuñoz-Cueto J, Elizur A, Kah O: Neuroendocrinology of reproduction in teleost fish. Gen Comp Endocrinol. 2010, 165 (3): 438-455.
Dufour S, Sebert ME, Weltzien FA, Rousseau K, Pasqualini C: Neuroendocrine control by dopamine of teleost reproduction. J Fish Biol. 2010, 76 (1): 129-160.
Amano M, Kitamura S, Ikuta K, Suzuki Y, Aida K: Activation of salmon GnRH mRNA expression prior to differentiation of precocious males in masu salmon. Gen Comp Endocrinol. 1997, 105 (3): 365-371.
Ando H, Sasaki Y, Okada H, Urano A: Prepubertal increases in the levels of two salmon gonadotropin-releasing hormone mRNAs in the ventral telencephalon and preoptic area of masu salmon. Neurosci Lett. 2001, 307 (2): 93-96.
Rodriguez L, Carrillo M, Sorbera LA, Soubrier MA, Mañanós E, Holland MCH, Zohar Y, Zanuy S: Pituitary levels of three forms of GnRH in the male European sea bass (Dicentrarchus labrax, L.) during sex differentiation and first spawning season. Gen Comp Endocrinol. 2000, 120 (1): 67-74.
Okuzawa K, Gen K, Bruysters M, Bogerd J, Gothilf Y, Zohar Y, Kagawa H: Seasonal variation of the three native gonadotropin-releasing hormone messenger ribonucleic acids levels in the brain of female red seabream. Gen Comp Endocrinol. 2003, 130 (3): 324-332.
Canosa LF, Stacey N, Peter RE: Changes in brain mRNA levels of gonadotropin-releasing hormone, pituitary adenylate cyclase activating polypeptide, and somatostatin during ovulatory luteinizing hormone and growth hormone surges in goldfish. Am J Physiol Regul Integr Comp Physiol. 2008, 295 (6): R1815-R1821.
Peñaranda DS, Pérez L, Gallego V, Jover M, Tveiten H, Baloche S, Dufour S, Asturiano JF: Molecular and physiological study of the artificial maturation process in European eel males: From brain to testis. Gen Comp Endocrinol. 2010, 166 (1): 160-171.
Schmitz M, Aroua S, Vidal B, Le Belle N, Elie P, Dufour S: Differential regulation of luteinizing hormone and follicle-stimulating hormone expression during ovarian development and under sexual steroid feedback in the European eel. Neuroendocrinology. 2005, 81 (2): 107-119.
Mechaly AS, Viñas J, Piferrer F: Sex-specific changes in the expression of kisspeptin, kisspeptin receptor, gonadotropins and gonadotropin receptors in the Senegalese sole (Solea senegalensis) during a full reproductive cycle. Comp Biochem Physiol A Mol Integr Physiol. 2012, 162 (4): 364-371.
Alvarado MV, Carrillo M, Felip A: Expression of kisspeptins and their receptors, gnrh-/gnrhr-II-1a and gonadotropin genes in the brain of adult male and female European sea bass during different gonadal stages. Gen Comp Endocrinol. 2013, 187: 104-116.
Jeng SR, Dufour S, Chang CF: Differential expression of neural and gonadal aromatase enzymatic activities in relation to gonadal development in Japanese eel, Anguilla japonica. J Exp Zool Part A. 2005, 303A (9): 802-812.
Jeng SR, Pasquier J, Yueh WS, Chen GR, Lee YH, Dufour S, Chang CF: Differential regulation of the expression of cytochrome P450 aromatase, estrogen and androgen receptor subtypes in the brain-pituitary-ovarian axis of the Japanese eel (Anguilla japonica) reveals steroid dependent and independent mechanisms. Gen Comp Endocrinol. 2012, 175 (1): 163-172.
Van den Thillart G, Rankin J, Dufour S: Spawning Migration of the European eel: Reproduction Index, a Useful Tool for Conservation Management. 2009, The Netherlands: Springer
van Ginneken VJT, Maes GE: The European eel (Anguilla anguilla, Linnaeus), its lifecycle, evolution and reproduction: a literature review. Rev Fish Biol Fisheries. 2005, 15 (4): 367-398.
Huertas M, Canário AVM, Hubbard PC: Chemical communication in the genus Anguilla: a minireview. Behaviour. 2008, 145: 1389-1407.
Sorensen PW, Winn HE: The induction of maturation and ovulation in American eels, Anguilla rostrata (LeSueur), and the relevance of chemical and visual cues to male spawning behaviour. J Fish Biol. 1984, 25 (3): 261-268.
Dufour S, Delerue-Le Belle N, Fontaine YA: Effects of steroid hormones on pituitary immunoreactive gonadotropin in European freshwater eel, Anguilla anguilla L. Gen Comp Endocrinol. 1983, 52 (2): 190-197.
Dufour S, Lopez E, Lemenn F, Lebelle N, Baloche S, Fontaine YA: Stimulation of gonadotropin-release and of ovarian development, by the administration of a gonadoliberin agonist and dopamine antagonists, in female silver eel pretreated with estradiol. Gen Comp Endocrinol. 1988, 70 (1): 20-30.
Vidal B, Pasqualini C, Le Belle N, Holland MCH, Sbaihi M, Vernier P, Zohar Y, Dufour S: Dopamine inhibits luteinizing hormone synthesis and release in the juvenile European eel: A neuroendocrine lock for the onset of puberty. Biol Reprod. 2004, 71 (5): 1491-1500.
Aroua S, Rousseau K, Schmitz M, Chang C-F, Dufour S: The gonadoliberin(S)-gonadotropin(S) axis in the eel: expression and regulation under induced maturation and sex steroid feedbacks. Spawning Migration of the European eel: Reproduction Index, a Useful Tool for Conservation Management, Volume 30. Edited by: van den Thillart G, Dufour S, Rankin JC. 2009, The Netherlands: Springer, 253-278.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402.
Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008, 4 (1): 44-57.
Schiöth HB, Fredriksson R: The GRAFS classification system of G-protein coupled receptors in comparative perspective. Gen Comp Endocrinol. 2005, 142: 94-101.
Kah O, Pontet A, Danger JM, Dubourg P, Pelletier G, Vaudry H, Calas A: Characterization, cerebral distribution and gonadotropin release activity of neuropeptide Y (NPY) in the goldfish. Fish Physiol Biochem. 1989, 7 (1–4): 69-76.
Li S, Zhao L, Xiao L, Liu Q, Zhou W, Qi X, Chen H, Yang H, Liu X, Zhang Y, Lin H: Structural and functional characterization of neuropeptide Y in a primitive teleost, the Japanese eel (Anguilla japonica). Gen Comp Endocrinol. 2012, 179 (1): 99-106.
Levavi-Sivan B, Bogerd J, Mañanós EL, Gómez A, Lareyre JJ: Perspectives on fish gonadotropins and their receptors. Gen Comp Endocrinol. 2010, 165 (3): 412-437.
Peng C, Humphries S, Peter RE, Rivier JE, Blomqvist AG, Larhammar D: Actions of goldfish neuropeptide Y on the secretion of growth hormone and gonadotropin-II in female goldfish. Gen Comp Endocrinol. 1993, 90 (3): 306-317.
Peng C, Huang YP, Peter RE: Neuropeptide Y stimulates growth hormone and gonadotropin release from the goldfish pituitary in vitro. Neuroendocrinology. 1990, 52 (1): 28-34.
Breton B, Mikolajczyk T, Popek W, Bieniarz K, Epler P: Neuropeptide Y stimulates in vivo gonadotropin secretion in teleost fish. Gen Comp Endocrinol. 1991, 84 (2): 277-283.
Peng C, Gallin W, Peter RE, Blomqvist AG, Larhammar D: Neuropeptide-Y gene expression in the goldfish brain: distribution and regulation by ovarian steroids. Endocrinology. 1994, 134 (3): 1095-1103.
Narnaware YK, Peter RE: Effects of food deprivation and refeeding on neuropeptide Y (NPY) mRNA levels in goldfish. Comp Biochem Physiol B Biochem Mol Biol. 2001, 129 (2–3): 633-637.
Peng C, Trudeau VL, Peter RE: Seasonal variation of neuropeptide Y actions on growth hormone and gonadotropin-II secretion in the goldfish: effects of sex steroids. J Neuroendocrinol. 1993, 5 (3): 273-280.
Tuomi J, Hakala T, Haukioja E: Alternative concepts of reproductive effort, costs of reproduction, selection in life-history evolution. Am Zool. 1983, 23: 25-34.
Trudeau VL, Spanswick D, Fraser EJ, Lariviere K, Crump D, Chiu S, MacMillan M, Schulz RW: The role of amino acid neurotransmitters in the regulation of pituitary gonadotropin release in fish. Biochem Cell Biol. 2000, 78 (3): 241-259.
Popesku JT, Martyniuk CJ, Mennigen J, Xiong H, Zhang D, Xia X, Cossins AR, Trudeau VL: The goldfish (Carassius auratus) as a model for neuroendocrine signaling. Mol Cell Endocrinol. 2008, 293 (1–2): 43-56.
Wang M: Neurosteroids and GABA-A receptor function. Front Endocrinol. 2011, 2: 44-
Zhang D, Xiong H, Mennigen JA, Popesku JT, Marlatt VL, Martyniuk CJ, Crump K, Cossins AR, Xia X, Trudeau VL: Defining global neuroendocrine gene expression patterns associated with reproductive seasonality in fish. PLoS One. 2009, 4 (6): e5816-
Lado WE, Zhang D, Mennigen JA, Zamora JM, Popesku JT, Trudeau VL: Rapid modulation of gene expression profiles in the telencephalon of male goldfish following exposure to waterborne sex pheromones. Gen Comp Endocrinol. 2013, 192: 204-213.
Stell BM, Brickley SG, Tang CY, Farrant M, Mody I: Neuroactive steroids reduce neuronal excitability by selectively enhancing tonic inhibition mediated by δ subunit-containing GABAA receptors. Proc Natl Acad Sci U S A. 2003, 100 (24): 14439-14444.
Chaouloff F: Physiopharmacological interactions between stress hormones and central serotonergic systems. Brain Res Rev. 1993, 18 (1): 1-32.
Winberg S, Winberg Y, Fernald RD: Effect of social rank on brain monoaminergic activity in a cichlid fish. Brain Behav Evol. 1997, 49 (4): 230-236.
Wendelaar Bonga SE: The stress response in fish. Physiol Rev. 1997, 77 (3): 591-625.
Huising MO, Metz JR, van Schooten C, Taverne-Thiele AJ, Hermsen T, Verburg-van Kemenade BML, Flik G: Structural characterisation of a cyprinid (Cyprinus carpio L.) CRH, CRH-BP and CRH-R1, and the role of these proteins in the acute stress response. J Mol Endocrinol. 2004, 32 (3): 627-648.
Suárez-Castillo EC, Garcia-Arrarás JE: Molecular evolution of the ependymin protein family: a necessary update. BMC Evol Biol. 2007, 7: 23-
Shashoua VE: Ependymin, a brain extracellular glycoprotein, and CNS plasticity. Ann N Y Acad Sci. 1991, 627 (1): 94-114.
Benowitz LI, Shashoua VE, Yoon MG: Specific changes in rapidly transported proteins during regeneration of the goldfish optic nerve. J Neurosci. 1981, 1 (3): 300-307.
Pradel G, Schachner M, Schmidt R: Inhibition of memory consolidation by antibodies against cell adhesion molecules after active avoidance conditioning in zebrafish. J Neurobiol. 1999, 39 (2): 197-206.
Xia ZG, Storm DR: The role of calmodulin as a signal integrator for synaptic plasticity. Nat Rev Neurosci. 2005, 6 (4): 267-276.
Fischer RM, Fontinha BM, Kirchmaier S, Steger J, Bloch S, Inoue D, Panda S, Rumpel S, Tessmar-Raible K: Co-Expression of VAL- and TMT-opsins uncovers ancient photosensory interneurons and motorneurons in the vertebrate brain. PLoS Biol. 2013, 11 (6): e1001585-
Nakane Y, Ikegami K, Iigo M, Ono H, Takeda K, Takahashi D, Uesaka M, Kimijima M, Hashimoto R, Arai N, Suga T, Kosuge K, Abe T, Maeda R, Senga T, Amiya N, Azuma T, Amano M, Abe H, Yamamoto N, Yoshimura T: The saccus vasculosus of fish is a sensor of seasonal changes in day length. Nat Commun. 2013, 4: 2108-
Vanveen T, Hartwig HG, Müller K: Light-dependent motor-activity and photonegative behavior in eel (Anguilla anguilla L.) Evidence for extraretinal extrapineal photoreception. J Comp Physiol. 1976, 111 (2): 209-219.
Haeseleer F, Imanishi Y, Maeda T, Possin DE, Maeda A, Lee A, Rieke F, Palczewski K: Essential role of Ca2+-binding protein 4, a Ca(v)1.4 channel regulator, in photoreceptor synaptic function. Nat Neurosci. 2004, 7 (10): 1079-1087.
Bonnard C, Strobl AC, Shboul M, Lee H, Merriman B, Nelson SF, Ababneh OH, Uz E, Güran T, Kayserili H, Hamamy H, Reversade B: Mutations in IRX5 impair craniofacial development and germ cell migration via SDF1. Nat Genet. 2012, 44 (6): 709-713.
Cheng CW, Chow RL, Lebel M, Sakuma R, Cheung HO-L, Thanabalasingham V, Zhang X, Bruneau BG, Birch DG, Hui C, McInnes RR, Cheng SH: The Iroquois homeobox gene, Irx5, is required for retinal cone bipolar cell development. Dev Biol. 2005, 287 (1): 48-60.
El Yakoubi W, Borday C, Hamdache J, Parain K, Tran HT, Vleminckx K, Perron M, Locker M: Hes4 controls proliferative properties of neural stem cells during retinal ontogenesis. Stem Cells. 2012, 30 (12): 2784-2795.
Guilgur LG, Moncaut NP, Canário AVM, Somoza GM: Evolution of GnRH ligands and receptors in gnathostoma. Comp Biochem Physiol A Mol Integr Physiol. 2006, 144 (3): 272-283.
Pankhurst NW, Lythgoe JN: Changes in vision and olfaction during sexual maturation in the European eel Anguilla anguilla (L.). J Fish Biol. 1983, 23 (2): 229-240.
Sorensen PW, Pankhurst NW: Histological changes in the gonad, skin, intestine and olfactory epithelium of artificially-matured male American eels, Anguilla rostrata (LeSueur). J Fish Biol. 1988, 32 (2): 297-307.
Huertas M, Scott AP, Hubbard PC, Canário AVM, Cerda J: Sexually mature European eels (Anguilla anguilla L.) stimulate gonadal development of neighbouring males: Possible involvement of chemical communication. Gen Comp Endocrinol. 2006, 147 (3): 304-313.
Pujolar JM, Marino IAM, Milan M, Coppe A, Maes GE, Capoccioni F, Ciccotti E, Bervoets L, Covaci A, Belpaire C, Cramb G, Patarnello T, Bargelloni L, Bortoluzzi S, Zane L: Surviving in a toxic world: transcriptomics and gene expression profiling in response to environmental pollution in the critically endangered European eel. BMC Genomics. 2012, 13 (1): 507-
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121.
The research that led to these results was funded by the Foundation for Science and Technology of Portugal (FCT), through project PTDC/MAR/113608/2009 and fellowships to AMC (SFRH/ BPD/85408/2012), JLS (SFRH/BPD/67008/2009), MH (SFRH/BPD/26339/2005), IAMM (CPDR123580/12) and RSTM (SFRH/BPD/66742/2009). The authors are grateful to John S. Taylor for his comments and suggestions during manuscript preparation and to Soraia Santos for her help with sample collection.
The authors declare that they have no competing interests.
AMC, MH and AVMC conceived and designed the study. JMP, IAMM and LZ were responsible for the microarray design. MH was responsible for injecting and maintaining the animals. Tissue sampling was done by AMC, PCH and MH. MM, IAMM, LB and TP conducted the microarray experiments and contributed to the acquisition of data. Data analysis and interpretation was done by AMC, PCH, JLS, JMP and MM. RSTM conducted the qPCR experiments and AMC prepared the manuscript with contributions from PCH, AVMC, JLS, JMP, MM, MH and LZ. All authors read and approved the final manuscript.