The Expression Pattern of Melatonin Receptor 1a Gene during Early Life Stages in the Nile tilapia (Oreochromis niloticus)

Ye Hwa Jin, Jin Woo Park, Jung-Hyun Kim, Joon Yeong Kwon†
Author Information & Copylight
Corresponding author: Joon Yeong Kwon, Department of Aquatic Life Medical Sciences, Sunmoon University, Asan 336-708, Korea. Phone: +82-41-530-2284, Fax: +82-41-530-2917, E-mail:

Copyright © 2013 © Korean Society of Developmental Biology. All Rights Reserved. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Feb 17, 2013 ; Revised: Mar 1, 2013 ; Accepted: Mar 11, 2013


The action of melatonin within the body of animals is known to be mediated by melatonin receptors. Three different types of melatonin receptors have been identified so far in fish. However, which of these are specifically involved in puberty onset is not known in fish. We cloned and analyzed the sequence of melatonin receptor 1a (mel 1a) gene in Nile tilapia Oreochromis niloticus. In addition, we examined the tissue distribution of gene expressions for three types of receptors, mel 1a, 1b and lc and investigated which of them is involved in the onset of puberty by comparing their expression with that of gonadotropin-releasing hormone receptor I (GnRHr I) gene using quantitative real-time PCR from 1 week post hatch (wph) to 24 wph. The mel 1a gene of Nile tilapia consisted of two exons and one bulky intron between them. Mel 1a gene was found to be highly conserved gene showing high homology with the corresponding genes from different teleost. All three types of melatonin receptor genes were expressed in the brain, eyes and ovary in common. Expression of mel 1a gene was the most abundant and ubiquitous among 3 receptors in the brain, liver, gill, ovary, muscle, eye, heart, intestine, spleen and kidney. Mel 1b and mel 1c genes were, however, expressed in fewer tissues at low level. During the development post hatch, expressions of both mel 1a and GnRHr I genes significantly increased at 13 wph which was close to the putative timing of puberty onset in this species. These results suggest that among three types of receptors mel 1a is most likely associated with the action of melatonin in the onset of puberty in Nile tilapia.


Melatonin mediates environmental cues such as light and photoperiod and consequently influences important physiological processes including reproduction in mammals (Malpaux et al., 2001; Pevet, 2003) and in teleosts (Mayer et al., 1997; Falcón et al., 2007). However, the role of melatonin in the process of teleost reproduction seems less clear due to the descrepancies observed in the previous studies. In European eel, Anguilla anguilla, the effect of melatonin treatment on reproduction was inhibitory (Sébert et al., 2008) while the same treatments were stimulatory in carp, Cyprinus catla (Bhattacharya et al., 2007), seabass, Dicentrarchus labrax (Bayarri et al., 2004) and zebrafish, Danio rerio (Carnevali et al., 2011). Even within the same species, the effect of melatonin administration was dependent on gonadal stage or sex of subjects (Bhattacharya et al., 2007). This indicates that the action of melatonin on fish reproduction is far more complicated than we thought.

The effects of melatonin within the body of animals are known to be mediated through melatonin receptors that belong to the G protein-coupled receptor (GPCR) superfamily (Dubocovich, 1995; Reppert et al., 1996). Supporting this, in one of teleost fish, seabass, the effects of melatonin rely not only on its rhythmic production, but also on the rhythmic expression of its receptors (Bayarri et al., 2010). There are three different subtypes of melatonin receptors, mel 1a, mel 1b and mel 1c. While the mel 1a and mel 1b subtypes are found in all vertebrates investigated, the mel 1c subtype has been found only in non-mammalian so far (Ebisawa et al., 1994; Reppert et al., 1995; Wiechmann et al., 1999). However, which of these are specifically involved in fish reproduction including puberty is not known.

Puberty is the developmental period when an individual matures sexually for the first time, and activation of brain-pituitary-gonadal axis (BPG axis) is prerequisite for puberty. This activation of BPG axis is initiated by an increase of gonadotropin-releasing hormone (GnRH). GnRH released from hypothalamus promotes the maturation of gonad by inducing release of luteinizing hormone (LH) and follicle-stimulating hormone (FSH) via its cognate receptor located in pituitary gonadotropic cells. GnRHs are key regulators in sexual development and reproduction in vertebrates and there are two or three piscine GnRH genes encoding GnRH peptides of different structures (Parhar, 2002). Among three types of GnRH identified in Nile tilapia Oreochromis niloticus (Parhar et al., 2003), GnRH I is located in the preoptic area (POA) and pituitary showing a close relationship with the activation of BPG axis (White et al., 1995; Parhar et al., 1996; Gothilf et al., 1996). Thus, the expression of GnRH receptor I (GnRHr I) gene could be indicative of the activation of BPG axis in this species.

In this study, we cloned and analyzed the sequence of melatonin receptor 1a gene in Nile tilapia and examined the tissue distribution of three receptors, mel 1a, mel 1b and mel lc genes. Finally, we investigated which receptor is involved in the onset of puberty by comparing their gene expression with GnRHr I gene expression using quantitative real-time PCR.


1. Rapid amplification of 5’/3’ cDNA ends (RACE) of mel 1a gene and sequence analysis

Mature male Nile tilapia (O. niloticus) was anaesthetized by benzocaine (50 ppm), and the whole brain was removed (body length: 16.2 cm, body weight: 70.5 g). The brain was homogenized with TRI REAGENT® (Molecular Research Center, Inc., Cincinnati, USA) to extract total RNA. Before RACE PCR, partial sequence of mel 1a was extended by RT-PCR using gene specific primers (Table 1) based on the mel 1a sequence of a close species: O. mossambicus (GenBank accession No. FJ531797.1). cDNA was synthesized from 2 μg of the total RNA with M-MLV reverse transcriptase (Promega, Madison, USA) and oligo d(T)15 primer (Promega, Madison, USA). PCR was carried out in 25 μl reaction mixture containing 1 μl cDNA, 400 nM forward and reverse primers (Table 1), 12.5 μl Go Taq Green Master Mix (Promega, Madison, USA) and nuclease-free water up to 25 μl. RT-PCR product was sequenced, and consequently the partial sequence (492 bp) was expanded to 858 bp. Based on the expanded mel 1a sequence, two forward primers and three reverse primers were designed (Table 1). RACE PCR was conducted using nTaq DNA Polymerase (Enzynomics, Daejeon, Korea) and 5'/3' RACE Kit, 2nd Generation (Roche, Mannheim, Germany) according to the protocol. RACE PCR products were cloned using TOPcloner TA kit and Chemically Competent E. coli (Enzynomics, Daejeon, Korea). Next, plasmid DNA was extracted and sequenced to find out the full sequence of mel 1a including 5' and 3'-UTR. The obtained mel 1a mRNA sequence was translated by ExPASy Translate Tool ( and the amino acid sequence was used to conserved domains analysis (http://www. Mel 1a amino acid sequences from various species were obtained from GenBank ( and aligned by ClustalW Multiple Alignment program and trimmed. A phylogenetic tree was formed using MEGA5 (Tamura et al., 2011) and Neighbor-Joining method (Saitou and Nei, 1987).

Table 1. Sequences of primers used to extend mel 1a partial sequence and RACE PCR
Extension of mel 1a Forward 5'-ACG CAG GGA ACA TCT TTG TG-3'
3' RACE Forward 1 5'-GCT GAC GCC ACA TGA CGT TA-3'
Forward 2 5'-TTC ATC GGA CTG GCT GTA GC-3'
5' RACE Reverse 1 5'-ATC GTG CCC AAC CTG TTT GT-3'
Reverse 2 5'-TTC CTC ATG GGT GTC AGC GT-3'
Reverse 3 5'-GCA GTA GCA GAC CTT GTG GT-3'
Download Excel Table
2. Tissue distribution of melatonin receptor genes (mel 1a, mel lb and mel lc)

To investigate the tissue specific expressions of mel 1a, mel 1b and mel 1c, total RNA was isolated from whole brain, liver, gill, ovary, muscle, eye, heart, intestine, spleen and kidney of an adult female tilapia (body length: 15.0 cm, body weight: 57.9 g) and subjected to RT-PCR using primers for mel 1a, mel 1b and mel 1c genes (Table 2). PCR step was a cycle of preheating at 94°C for 2 min and 35 cycles of amplification (for mel 1a, mel 1b and mel 1c genes) or 25 cycles (for β-actin gene) including denaturation at 94°C for 30 sec, annealing (at 53°C for mel 1a, 57°C for mel 1b, 53°C for mel 1c or 55°C for β-actin gene) for 1 min and extension at 72°C for 1 min, followed by a cycle of final extension at 72°C for 5 min. The PCR products were electrophoresed on 1% agarose gels and stained with ethidium bromide. The products were, then, visualized under UV light and analyzed by Image Analysis System (Kodak DC290).

Table 2. Sequences of primers used for expression studies
Sequence Product size (bp)
Tissue distribution (RT-PCR) mel 1a F 5'-ACG CAG GGA ACA TCT TTG TG-3' 318
mel 1b F 5'-TCG CCC AGA ATG TCA GCA GT-3' 221
mel 1c F 5'-TAC CCC TAC CCA CTG GTC CT-3' 234
β-actin F 5'-AAT CGT GCG TGA CAT CAA GG-3' 392
β-actin R 5'-AGT ATT TAC GCT CAG GTG GG-3'
Quantitative real-time PCR qmel 1a F 5'-GGT TGT GAT TCC CCT CAT CCC-3' 105
qβ-actin F 5'-AGC ATC CCG TCC TGC TCA CA-3' 121
qβ-actin R 5'-AGC ACA GCC TGG ATG GCA AC-3'

GnRHr I gene primer from Martinez-Chavez et al. (2008). Primers for mel 1a, mel 1b (GenBank accession No. XM003446809.1), mel 1c (Kim et al., unpublished) and β-actin (GenBank accession No. XM003443127.1) genes were newly designed.

Download Excel Table
3. Expression pattern of mel 1a gene during early development

Experimental fish were produced from Nile tilapia obtained from Chungju, Chungcheongbuk-do in Korea. They were reared in a closed recirculating aquaculture system at 27 ± 1°C under natural photoperiod (sampling period: July to December) and fed twice a day. Fish from a batch were sampled 1 to 4 week intervals at noon (12 pm) from 1 week post-hatch (wph) to 24 wph. Depending on the growth of the fish, their heads (1—7 wph) or brains (10—24 wph) were sampled (n=4). Protocols of total RNA extraction and cDNA synthesis are the same as above. We examined the expression of GnRHr I and mel 1a genes together with the expression of β-actin gene for normalization. Quantitative real-time PCR reactions were conducted using QuantiMix SYBR (PhileKorea, Seoul, Korea) and Eco Real-Time PCR System (illumina®, San Diego, USA) in a 20 μl reaction mixture containing 1 μl cDNA, 10 μl QuantiMix SYBR, 250 nM forward and reverse primers (Table 2) and nuclease-free water up to 20 μl. Melting curve analysis was also included at the end of the PCR. Each expression value of the mRNAs was presented as Ct (threshold cycles) values. The levels of the mRNAs were normalized against the amount of β-actin mRNA and the relative level was determined by the comparative threshold cycle method, 2-ΔCt.

4. Statistical analyses

Data were presented as mean ± S.E.M. For statistical analysis of mel 1a and GnRHr I mRNA expression, expression level observed at each wph was compared with the level observed immediately previous wph within this study by independent t-test (p<0.05) using SPSS version 18.0 for windows.


1. Mel 1a gene sequence

The sequence of mel 1a mRNA cloned from the brain of Nile tilapia consisted of an ORF (1053 nt), a 82 nt 5'-UTR and a 49 nt 3'-UTR (Fig. 1). Deduced amino acid sequence of mel 1a contained a signal peptide region, predicted transmembrane domains and conserved N-R-Y motif. By comparison of mel 1a mRNA sequence with that of genome sequence in this species (GenBank accession No. NT167456), mel 1a gene of O. niloticus was found to be comprised of two relatively short exons (265 and 936 bp, respectively) and one bulky intron of 35,231 bp between the two exons (Fig. 2).

Fig. 1. Nucleotide and deduced amino acid sequence of O. niloticus mel 1a. The length of ORF is 1053 nt and encodes a protein of 350 amino acids, leaving a 82 nt 5'-UTR and a 49 nt 3'-UTR. Start and stop (*) codons are boxed. Predicted trans-membrane (TM1-7) domains are underlined. The gray background shows signal peptide. The dotted box is the conserved NRY motif.
Download Origianl Figure
Fig. 2. Mel 1a gene structure of O. niloticus. It is comprised of 2 exons of 265 bp and 936 bp, respectively. Between exon1 and exon2, there is one bulky intron of 35,231 bp.
Download Origianl Figure

Multiple alignment of Nile tilapia mel 1a amino acid sequence with those of other fish species and subsequent phylogenetic analysis (Fig. 3) showed high similarity of this gene with that of other fish species such as Mozambique tilapia (O. mossambicus, 99%), European seabass (D. labrax, 97%), rabbitfish (Siganus guttatus, 96%) and grass puffer (Takifugu niphobles, 96%).

Fig. 3. Phylogenetic analysis of mel 1a genes from various vertebrates. This analysis was performed with multiple alignments from amino acid sequences using the ClustalW program and neighbor joining method. GenBank accession Numbers for the respective animals: human (Homo sapiens) NP005949.1; mouse (Mus musculus) NP032665.1; chicken (Gallus gallus) AAA92498.1; zebrafish (Danio rerio) NP571468.1; rabbitfish (Siganus guttatus) ABG77572.1; Mozambique tilapia (O. mossambicus) ACL78767.1; European seabass (Dicentrarchus labrax) ACB13280.1; rainbow trout (Oncorhynchus mykiss) AAF00191.1; sweetfish (Plecoglossus altivelis) BAI65858.1
Download Origianl Figure
2. Tissue distribution of melatonin receptor genes (mel 1a, lb and lc)

Tissue specific expression of melatonin receptors analyzed by RT-PCR was shown in Fig. 4. Mel 1a gene was expressed abundantly and ubiquitously in all tissues investigated including the brain and various peripheral tissues: liver, gill, ovary, muscle, eye, heart, intestine, spleen and kidney. In contrast, mel 1b and mel 1c genes were faintly detected in only several tissues. Mel 1b transcript existed in the brain, gill, ovary, eye and kidney, and mel 1c transcript was distributed in the brain, gill, ovary, eye, heart and intestine. Three different melatonin receptors were commonly detected in the brain, gill, ovary and eye.

Fig. 4. Expression of mel 1a, mel 1b and mel 1c gene from various tissues (B: brain, L: liver, G: gill, O: ovary, M: muscle, E: eye, H: heart, I: intestine, S: spleen and K: kidney) of Nile tilapia. The expressions of melatonin receptor genes and β-actin gene (as a control) were examined by using RT-PCR. PCR products were run on 1% agarose gels.
Download Origianl Figure
3. Expression pattern of mel 1a gene during early development

The relative expression level of mel 1a gene significantly increased at 4 wph and decreased at 7 wph compared with the respective previous investigation. Then it remained low until 12 wph and increased significantly at 13 wph and 14 wph compared with the previous investigation (Fig. 5A). The relative expression level constantly increased until 15 wph and reached the highest level though not significantly different from the level at 14 wph. Then, it declined rapidly at 16 wph and 20 wph. On the whole, the level of mel 1a transcript exhibited a prominent increase between 13—16 wph.

Fig. 5. Relative expression pattern of mel 1a mRNA (A) and GnRHr I mRNA (B) in the brain of Nile tilapia from 1 wph to 24 wph. Differences in the mRNA levels were estimated by quantitative real-time PCR. The relative abundance of the mRNAs was normalized to the amount of β-actin mRNA by the comparative threshold cycle method. Each value represents mean ± SEM (n = 4). Asterisks (*) indicate significant difference from the immediately previous observation within this study (P<0.05).
Download Origianl Figure

The first significant increase in the level of relative GnRHr I gene expression was observed at 13 wph (Fig. 5B). Then, the expression level maintained up to 15 wph though it tended to be decreased. Overall, expression profile of mel 1a gene that exhibited higher expression between 13—16 wph approximately corresponded to the profile of GnRHr I gene.


1. Mel 1a gene sequence

Mel 1a gene cloned from Nile tilapia was a rhodopsin like GPCR and had 7 transmembrane domains. While most of rhodopsin-like GPCRs have a consensus sequence of D/E-R-Y at the end of third transmembrane domain, melatonin receptors have a conserved sequence of N-R-Y at the same location (Kokkola et al., 2005). Anterior two residues of N-R-Y motif are known to be involved directly in the activation of rhodopsin and β-adrenergic receptors (Ballesteros et al., 2001; Gether et al., 2002). The mel 1a gene structure of O. niloticus that had two exons and one bulky intron (35 kilobases) was similar to that of mouse (Roca et al., 1996). In addition, the mel 1a gene in tilapia showed high similarity to the same gene from different teleost species. These data indicates that mel 1a is highly conserved gene.

The coding sequence (CDS) of mel 1a gene identified in this study was 100% identical to CDS of predicted O. niloticus melatonin receptor subtype 1A-A-like (GenBank accession No. XM003449697) predicted from the genomic sequence (NT167456) using GNOMON method. The presence of 1A-like sequence (GenBank accession No. XM003447612.1) distinguished from 1A-A-like suggests the possible existence of two subtypes of mel 1a gene in this species. It is necessary to confirm this possibility by further studies. So far, teleost species which have two subtypes of mel 1a are zebrafish (Reppert et al., 1995), rainbow trout (Mazurais et al., 1999) and goldfish (Ikegami et al., 2009).

2. Tissue distribution of melatonin receptor genes (mel 1a, lb and lc)

All three types of Nile tilapia melatonin receptor genes were expressed in the brain, eye and ovary in common. Especially, they were abundant in the brain and eye. Even though these receptors were expressed abundantly in the brain and eye in common, mel 1a mRNA was the most abundant suggesting that this receptor mediate most of melatonin actions at the level of central nervous system. Mel 1a gene was also prominently expressed in all peripheral tissues investigated in this study, including liver, gill, ovary, muscle, heart, intestine, spleen and kidney while mel 1b and mel 1c genes were expressed in fewer tissues. The functional significance of peripheral melatonin receptors mediating the melatonin effects needs further investigation.

3. Expression pattern of mel 1a gene during early development

In this study, expression of GnRHr I gene was used to verify the timing of puberty onset in Nile tilapia, and it was compared with the expression pattern of mel 1a gene since the increase of GnRHr I gene expression during puberty had been frequently reported in fish (Mohamed et al., 2007; Martinez-Chavez et al., 2008; Park et al., 2012). Martinez-Chavez et al. (2008) reported 9—13 wph as the timing of puberty onset in Nile tilapia by putting together the results of gonadal histological observation and the significant expression of GnRHr I. In the present study, the first significant rise of GnRHr I gene expression occurred at 13 wph, implying the timing of puberty onset. Difference in the timing of puberty onset between these two studies may have resulted from different rearing conditions, strains and growth rates between two experimental systems (Park et al., 2012).

Relative expression of both mel 1a and GnRHr I gene increased significantly at 13 wph, the putative timing of puberty. This coinciding elevation of both gene expression indicates melatonin/mel 1a system may have a role in puberty onset through GnRH system. In male masu salmon Oncorhynchus masou pituitary, GnRH and LH levels were decreased by melatonin treatment, resulting the decrease of gonado-somatic index (GSI) and plasma T level (Amano et al., 2004). Similarly, melatonin administration in A. anguilla partly inhibited gonadal development by stimulating dopaminergic inhibitory actions on pituitary-gonad axis (Sébert et al., 2008). In zebrafish, however, melatonin administration via water led to a significant increase of GnRH III gene expression in the brain and LH gene expression in the pituitary resulting the increase of GSI (Carnevali et al., 2011). In addition, numerous melatonin-binding sites were distributed in POA where GnRH and dopamine neurons were located in various teleost species (reviewed by Dufour et al., 2010) and in tilapia hybrid (O. aureus × O. niloticus) (Levavi-Sivan et al., 2004). Despite of the discrepancies between studies, all of these data commonly indicate that melatonin would have negative or positive effects on reproduction via GnRH/dopamine system in teleost.

The concomitant increase of both GnRHr I and mel 1a gene expression suggests the possible involvement of melatonin in the onset of puberty via mel 1a. In support of this, in male seabass, the plasma melatonin level increased during the first sexual maturation (Bayarri et al., 2010). In Senegalese sole Solea senegalensis exhibiting annual reproductive rhythms, a seasonal variation of mel 1a gene expression in pituitary was observed, showing an increase after main spawning season and decrease in second spawning period (Confente et al., 2010). Changes of plasma melatonin level positively corresponded to the seasonal variation in the same species (Vera et al., 2007).

In summary, we have cloned melatonin receptor 1a gene from the brain of Nile tilapia. This gene was highly conserved among species based on the phylogenetic and sequence analysis. Mel 1a mRNA was most abundant in all tissues investigated in this study, suggesting that the considerable melatonin actions are mediated by mel 1a in these peripheral tissues. Gene expression profiles of GnRHr I revealed a putative timing of puberty onset at 13 wph, and expression of mel 1a gene showed a close relationship with puberty onset in this species. Further studies should clarify the exact function of mel 1a gene that was concomitantly increased with GnRHr I gene during the putative puberty onset.


This work was supported by Korea Research Foundation Grant funded by the Korean Government (KRF-2008-521-F00038).



Amano M, Iigo M, Ikuta K, Kitamura S, Okuzawa K, Yamada H, Yamamori K. Disturbance of plasma melatonin profile by high dose melatonin administration inhibits testicular maturation of precocious male masu salmon. Zool Sci. 2004; 21 p. 79-85.


Ballesteros JA, Jensen AD, Liapakis G, Rasmussen SG, Shi L, Gether U, Javitch JA. Activation of the beta 2-adrenergic receptor involves disruption of an ionic lock between the cytoplasmic ends of transmembrane segments 3 and 6. J Biol Chem. 2001; 276 p. 29171-29177.


Bayarri MJ, Falcón J, Zanuy S, Carrillo M. Continuous light and melatonin daily and seasonal variations of brain binding sites and plasma concentration during the first reproductive cycle of sea bass. Gen Comp Endocrinol. 2010; 169 p. 58-64.


Bayarri MJ, Rodriguez L, Zanuy S, Madrid JA, Sanchez- Vazquez FJ, Kagawa H, Okuzawa K, Carrillo M. Effect of photoperiod manipulation on the daily rhythms of melatonin and reproductive hormones in caged European sea bass (Dicentrarchus labrax). Gen Comp Endocrinol. 2004; 136 p. 72-81.


Bhattacharya S, Chattoraj A, Maitra SK. Melatonin in the regulation of annual testicular events in carp Catla catla evidence from the studies on the effects of exogenous melatonin, continuous light, and continuous darkness. Chronobiol Int. 2007; 24 p. 629-650.


Carnevali O, Gioacchini G, Maradonna F, Olivotto I, Migliarini B. Melatonin Induces follicle maturation in Danio rerio. PLoS One. 2011; 6 p. e19978.


Confente F, Rendón MC, Besseau L, Falcón J, Muñoz-Cueto JA. Melatonin receptors in a pleuronectiform species, Solea senegalensis: Cloning, tissue expression, day-night and seasonal variations. Gen Comp Endocrinol. 2010; 167 p. 202-214.


Dubocovich ML. Melatonin receptors: are there multiple subtypes?. Trends Pharmacol Sci. 1995; 16 p. 50-56.


Dufour S, Sébert ME, Weltzien FA, Rousseau K, Pasqualini C. Neuroendocrine control by dopamine of teleost reproduction. J Fish Biol. 2010; 76 p. 129-160.


Ebisawa T, Karne S, Lerner MR. Expression cloning of a high-affinity melatonin receptor from Xenopus dermal melanophores. Proc Natl Acad Sci USA. 1994; 91 p. 6133-6137.


Falcón J, Besseau L, Sauzet S, Boeuf G. Melatonin effects on the hypothalamo-pituitary axis in fish. Trends Endocrinol Metab. 2007; 18 p. 81-88.


Gether U, Asmar F, Meinild AK, Rasmussen SGF. Structural basis for activation of G-protein-coupled receptors. Pharmacol Toxicol. 2002; 91 p. 304-312.


Gothilf Y, Munoz-Cueto JA, Sagrillo CA, Selmanoff M, Chen TT, Kah O, Elizur A, Zohar Y. Three forms of gonadotropin-releasing hormone in a perciform fish (Sparus aurata): complementary deoxyribonucleic acid characterization and brain localization. Biol Reprod. 1996; 55 p. 636-645.


Ikegami T, Azuma K, Nakamura M, Suzuki N, Hattori A, Ando H. Diurnal expressions of four subtypes of melatonin receptor genes in the optic tectum and retina of goldfish. Comp Biochem Physiol A. 2009; 152 p. 219-224.


Kokkola T, Salo OMH, Poso A, Laitinen JT. The functional role of cysteines adjacent to the NRY motif of the human MT1 melatonin receptor. J Pineal Res. 2005; 39 p. 1-11.


Levavi-Sivan B, Safarian H, Rosenfeld H, Elizur A, Avitan A. Regulation of gonadotropin-releasing hormone (GnRH)-receptor gene expression in tilapia: effect of GnRH and dopamine. Biol Reprod. 2004; 70 p. 1545-1551.


Malpaux B, Migaud M, Tricoire H, Chemineau P. Biology of mammalian photoperiodism and the critical role of the pineal gland and melatonin. J Biol Rhythms. 2001; 16 p. 336-347.


Martinez-Chavez CC, Minghetti M, Migaud H. GPR54 and rGnRH I gene expression during the onset of puberty in Nile tilapia. Gen Comp Endocrinol. 2008; 156 p. 224-233.


Mayer I, Bornestaf C, Borg B. Melatonin in nonmammalian vertebrates physiological role in reproduction?. Comp Biochem Physiol A. 1997; 118 p. 515-531.


Mazurais D, Brierley I, Anglade I, Drew J, Randall C, Bromage N, Michel D, Kah O, Williams LM. Central melatonin receptors in the rainbow trout comparative distribution of ligand binding and gene expression. J Comp Neurol. 1999; 409 p. 313-324.


Mohamed JS, Benninghoff AD, Holt GJ, Khan IA. Developmental expression of the G protein-coupled receptor 54 and three GnRH mRNAs in the teleost fish cobia. J Mol Endocrinol. 2007; 38 p. 235-244.


Parhar IS, Pfaff DW, Schwanzel-Fukuda M. Gonadotropinreleasing hormone gene expression in teleosts. Mol Brain Res. 1996; 41 p. 216-227.


Parhar IS. Cell migration and evolutionary significance of GnRH subtypes. Prog Brain Res. 2002; 141 p. 3-17.


Parhar IS, Ogawa S, Hamada T, Sakuma Y. Singlecell real-time quantitative polymerase chain reaction of immunofluorescently identified neurons of gonadotropinreleasing hormone subtypes in cichlid fish. Endocrinology. 2003; 144 p. 3297-3300.


Park JW, Kim JH, Jin YH, Kwon JY. Expression Profiles of Kiss2, GPR54 and GnRH Receptor? mRNAs in the Early Life Stage of Nile Tilapia, Oreochromis niloticus. Dev Reprod. 2012; 16 p. 31-38.


Pevet P. Melatonin from seasonal to circadian signal. J Neuroendocrinol. 2003; 15 p. 422-426.


Reppert SM, Weaver DR, Cassone VM, Godson C, Kolakowski LF. Melatonin receptors are for the birds molecular analysis of two receptor subtypes differentially expressed in chick brain. Neuron. 1995; 15 p. 1003-1015.


Reppert SM, Weaver DR, Godson C. Melatonin receptors step into the light cloning and classification of subtypes. Trends Pharmacol Sci. 1996; 17 p. 100-102.


Roca AL, Godson C, Weaver DR, Reppert SM. Structure, characterization, and expression of the gene encoding the mouse mel 1a melatonin receptor. Endocrinology. 1996; 137 p. 3469-3477.


Saitou N, Nei M. The neighbor-joining method a new method: for reconstructing phylogenetic trees. Mol Biol Evol. 1987; 4 p. 406-425.


Sébert ME, Legros C, Weltzien FA, Malpaux B, Chemineau P, Dufour S. Melatonin activates brain dopaminergic systems in the eel with an inhibitory impact on reproductive function. J Neuroendocrinol. 2008; 20 p. 917-929.


Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011; 28 p. 2731-2739.


Vera LM, De Oliveira C, López-Olmeda JF, Ramos J, Mañanós E, Madrid JA, Sánchez-Vázquez FJ. Seasonal and daily plasma melatonin rhythms and reproduction in Senegal sole kept under natural photoperiod and natural or controlled water temperature. J Pineal Res. 2007; 43 p. 50-55.


White SA, Kasten TL, Bond CT, Adelman JP, Fernald RD. Three gonadotropin-releasing hormone genes in one organism suggest novel roles for an ancient peptide. Proc Natl Acad Sci USA. 1995; 92 p. 8363-8367.


Wiechmann AF, Campbell LD, Defoe DM. Melatonin receptor RNA expression in Xenopus retina. Mol Brain Res. 1999; 297 p. 297-303.