help button home button Endocrine Society Endocrinology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Endocrinology, doi:10.1210/en.2006-1400
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Flanagan, C. A.
Right arrow Articles by Illing, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Flanagan, C. A.
Right arrow Articles by Illing, N.
Endocrinology Vol. 148, No. 10 5060-5071
Copyright © 2007 by The Endocrine Society

Expression, Structure, Function, and Evolution of Gonadotropin-Releasing Hormone (GnRH) Receptors GnRH-R1SHS and GnRH-R2PEY in the Teleost, Astatotilapia burtoni

Colleen A. Flanagan, Chun-Chun Chen, Marla Coetsee, Sipho Mamputha, Kathleen E. Whitlock, Nicholas Bredenkamp, Logan Grosenick, Russell D. Fernald and Nicola Illing

University of Cape Town/Medical Research Council Research Group for Receptor Biology, Institute of Infectious Diseases and Molecular Medicine and Division of Medical Biochemistry (C.A.F., M.C., S.M.), and Department of Medicine (C.A.F.), University of Cape Town Faculty of Health Sciences, Observatory 7925, South Africa; Department of Biological Sciences and Program in Neuroscience (C.-C.C., L.G., R.D.F.), Stanford University, Stanford, California 94305; Department of Molecular Biology and Genetics (K.E.W.), Cornell University, Ithaca, New York 14853; and Department of Molecular and Cell Biology (N.B., N.I.), University of Cape Town, Rondebosch 7700, South Africa

Address all correspondence and requests for reprints to: R. D. Fernald, Department of Biological Sciences and Program in Neuroscience, Stanford University, Stanford, California 94305-2130. E-mail: rfernald{at}stanford.edu.


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Multiple GnRH receptors are known to exist in nonmammalian species, but it is uncertain which receptor type regulates reproduction via the hypothalamic-pituitary-gonadal axis. The teleost fish, Astatotilapia burtoni, is useful for identifying the GnRH receptor responsible for reproduction, because only territorial males reproduce. We have cloned a second GnRH receptor in A. burtoni, GnRH-R1SHS (SHS is a peptide motif in extracellular loop 3), which is up-regulated in pituitaries of territorial males. We have shown that GnRH-R1SHS is expressed in many tissues and specifically colocalizes with LH in the pituitary. In A. burtoni brain, mRNA levels of both GnRH-R1SHS and a previously identified receptor, GnRH-R2PEY, are highly correlated with mRNA levels of all three GnRH ligands. Despite its likely role in reproduction, we found that GnRH-R1SHS has the highest affinity for GnRH2 in vitro and low responsivity to GnRH1. Our phylogenetic analysis shows that GnRH-R1SHS is less closely related to mammalian reproductive GnRH receptors than GnRH-R2PEY. We correlated vertebrate GnRH receptor amino acid sequences with receptor function and tissue distribution in many species and found that GnRH receptor sequences predict ligand responsiveness but not colocalization with pituitary gonadotropes. Based on sequence analysis, tissue localization, and physiological response we propose that the GnRH-R1SHS receptor controls reproduction in teleosts, including A. burtoni. We propose a GnRH receptor classification based on gene sequence that correlates with ligand selectivity but not with reproductive control. Our results suggest that different duplicated GnRH receptor genes have been selected to regulate reproduction in different vertebrate lineages.


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
IN VERTEBRATES, REPRODUCTION is regulated by photoperiod, physical environment, and social status, which act via a final common pathway: the hypothalamic-pituitary-gonadal (HPG) axis. The hypothalamus delivers GnRH (GnRH; see Materials and Methods for peptide nomenclature) to the pituitary causing release of gonadotropins (LH and FSH) that travel to the gonads to regulate gonadal growth and release of sex steroids. The effects of this signaling cascade depend critically on the response of the pituitary GnRH1 receptor. GnRH1 has also been associated with pituitary control of GH (1, 2) and prolactin secretion (3), although the details of these pathways are not as well understood.

Mammalian GnRH (mGnRH1) was first identified in the mammalian hypothalamus, and additional GnRH forms were discovered in nonmammalian vertebrates including teleost fish, which express three GnRH forms: GnRH1, GnRH2, and GnRH3 (4, 5). All species studied to date, including mammals, have two or three forms of GnRH produced by distinct but phylogenetically related genes (5, 6, 7, 8, 9). The amino acid sequences of GnRH2 and GnRH3 are conserved, but the structure of GnRH1 varies considerably across vertebrate species (5, 6, 8). The presence of two or three GnRH peptides within a single organism raises the possibility that GnRH receptors (GnRH-R) may have coevolved with their ligands.

Multiple GnRH-R types have been reported in mammals, birds, fish, and amphibians (8, 10, 11, 12). GnRH-Rs belong to the G protein-coupled receptor family, characterized by seven hydrophobic transmembrane domains linked by hydrophilic extra- and intracellular loops (8, 13, 14, 15). The first GnRH-R sequences were cloned from mammalian pituitary mRNA (8, 13, 14) and designated type I GnRH-Rs because they regulate reproductive function (16, 17). These receptors exhibit high affinity for mGnRH1 and low affinity for other naturally occurring forms of GnRH. A conserved amino acid motif (Ser-Asp/Glu-Pro, SDP), in extracellular loop 3 (EC3) of mammalian pituitary type I GnRH-Rs, determines their specificity for mGnRH1 (18, 19). Receptors that lack the SDP motif in EC3 are not specific for mGnRH1 (8, 10, 20, 21, 22) (see GnRH-R amino acid alignment, supplemental Fig. A, published as supplemental data on The Endocrine Society’s Journals Online web site at http://endo.endojournals.org). GnRH-Rs with a PEY motif in EC3 (GnRH-R1PEY) were initially classified as a nonmammalian type I GnRH-R, and they were predicted to regulate reproduction in response to GnRH1 (10). GnRH-Rs with a PPS motif were originally proposed to be the receptors for GnRH2 (10), and indeed, mammalian type II GnRH-Rs, which are specific for GnRH2, have a PPS motif (8, 23, 24).

Unfortunately, GnRH-Rs have been named numerically in the order of their discovery, without consideration of functional or evolutionary relationships (11, 25, 26, 27). Several GnRH-R nomenclature systems have been proposed that, although based on phylogenetic relationships, contradict one another (supplemental Table A, published as supplemental data on The Endocrine Society’s Journals Online web site at http://endo.endojournals.org). For example, the Astatotilapia burtoni GnRH-R1SHS described here was initially named GnRH-R2, following Troskie et al. (10) but was renamed as a type I GnRH-R to be consistent with Lethimonier et al. (12). Classification of receptors based on function has been hampered by lack of physiological evidence about regulation of the HPG axis in nonmammalian vertebrates.

We used an African cichlid fish species, A. burtoni, in which reproduction is socially regulated, to identify which GnRH-R regulates reproductive function and to clarify the evolutionary relationship among known receptors. A. burtoni males exhibit one of two distinct social phenotypes: dominant (territorial, T) animals, which are larger, brightly colored, and reproductively capable or subordinate (nonterritorial, NT) males, which are smaller and camouflaged, have regressed gonads, and school with females (28, 29). Shifts in social status produce a corresponding change in the gonads via the HPG axis (29, 30). We have previously shown that [Ser8] GnRH1 is the only form of GnRH found in the pituitary of A. burtoni, and it is up-regulated in T males, implying that it must regulate gonadotropin release (31).

We previously cloned a GnRH-R2PEY from A. burtoni and showed that it has a relatively poor response to [Ser8] GnRH1 in vitro (32) and is expressed in the dorsal-anterior and posterior pituitary (33). We also identified a partial sequence corresponding to EC3 of a second GnRH-R, now designated GnRH-R1SHS (32), which is expressed in the ventral anterior and posterior pituitary and up-regulated in T males compared with NT males (33, 34). We now report the complete cloning and tissue distribution of GnRH-R1SHS and compare the functional activities and regulation of GnRH-R2PEY and GnRH-R1SHS to determine whether expression and structure-function relationships support the identification of GnRH-R1SHS as the receptor regulating reproduction in A. burtoni.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Nomenclature of GnRH ligands and receptors
GnRH1 denotes the peptide that regulates reproduction in vertebrates, and we use mGnRH1 to distinguish the mammalian form. For other GnRH peptide sequences, we show amino acid changes relative to the mGnRH1 sequence (35). GnRH peptides used in this study were mGnRH1 (pGlu-His-Trp-Ser-Tyr-Gly-Leu-Arg-Pro-Gly-NH2) and the three forms of GnRH found in A. burtoni, [Ser8] GnRH1 (also called seabream GnRH), GnRH2 ([His5,Trp7,Tyr8] GnRH) (36), and GnRH3 ([Trp7,Leu8] GnRH, also called salmon GnRH) (35, 37). GnRH-Rs are designated as in GenBank with the addition of a three-letter superscript, which is the conserved motif in EC3 (SHS, PEY, PPS, or SDP) to indicate their GnRH-R type classification. Hence, the receptor identified by our experiments becomes GnRH-R1SHS, because it was designated GnRH-R1 in GenBank and contains the SHS sequence in EC3. Correspondingly, the other A. burtoni GnRH-R is named GnRH-R2PEY. In receptors where the EC3 motif is not identical to the conserved motif, the motif in the receptor is used (e.g. GnRH-R1PDY for one of the catfish receptors, which has PDY instead of PEY in EC3).

Animals
A. (Haplochromis) burtoni derived from a wild-caught population were raised in aquaria under conditions similar to their native equatorial habitat in Lake Tanganyika in Africa (pH 8, 28 C) (28). Both male and female fish (3–4 cm length) were used in these experiments. To collect brain tissue, fish were killed via rapid cervical transection. All work was performed in compliance with Stanford University (APLAC) guidelines.

RNA isolation, amplification, sequencing, and analysis of A. burtoni GnRH-Rs
Partial sequences corresponding to the cichlid GnRH-R1SHS were originally amplified using degenerate primers designed to conserved regions of transmembrane TM6 and TM7 (10, 38). We subsequently designed a pair of degenerate primers to conserved amino acid sequences identified by an alignment with amberjack (CAB65407), striped bass (AAF28464) (32), medaka (BAB7056), and European seabass (DLA419594) GnRH-Rs. Sense primer S2 [5'-gtgggcigcica(t/c)(a/t)(g/c)iga(c/t)ggiaa-3'] was designed to the conserved motif WAAHSDGK, and the antisense primer AS1 [5'-gtticc(c/t)tcia(a/g)(g/a)tc(g/a)tcigg(g/a)aa-3'] was designed to the conserved motif FPDDLEGK. These primers were used to amplify a 750-nucleotide product from A. burtoni brain cDNA that was identified as a partial sequence for GnRH-R1SHS. Sequence information from this clone was used to design gene-specific primers for 5' RACE and 3' RACE to clone the full-length GnRH-R1SHS cDNA.

5' RACE and 3' RACE cDNA were prepared according to the manufacturer’s instructions (SMART RACE cDNA Amplification Kit; BD Biosciences, Franklin Lakes, NJ), from total RNA extracted from A. burtoni brain. Antisense primers CH3gsp (5'-tgcaggcaaagtctccggcaagc-3') and CH3ngsp (5'-aatcagcaccctcacgtgggattttcg-3') were used in combination with the 5' universal (UPM) and nested universal (NUP) primers to amplify the 5' end of GnRH-R1 cDNA from the 5' RACE A. burtoni brain cDNA. The 3' end of the cichlid GnRH-R1SHS receptor was amplified using the sense primers 3' gsp1 (5'-gcccgagagcccggatgagaactctg-3') and 3' ngsp1 (5'-gtgattgttctgtctttcatcatctg-3') in combination with the universal (UPM) and nested universal (NUP) from the 3' RACE A. burtoni brain cDNA. Commercial kits (BD Advantage 2 PCR Enzyme System; BD Biosciences) were used according to the manufacturer’s instructions for both 5' and 3' PCR amplification reactions.

Sequence information from the 5' RACE and 3' RACE cDNA clones was used to design a sense and antisense primer that, respectively, included the start and stop codons of the GnRH-R1SHS mRNA transcript, and which also included either an additional BamHI or XhoI site (italic): FLs 5'-gcgcggatccaccatggtggatgggcac-3' and FLas 5'-ggccggctcgagtcataagatgctctcag-3'. The full-length GnRH-R1SHS coding region was amplified from A. burtoni brain cDNA, using this pair of primers and subcloned into a BamHI/XhoI-digested pcDNA3.1(+) expression vector (Invitrogen). Expression plasmids were purified using the PureYield Plasmid Midiprep System according to the manufacturer’s alternative protocol (Promega) and sequenced using T7 and BGH reverse primer (Invitrogen, Carlsbad, CA). The nucleotide sequence of the full-length coding region is in GenBank [AY705931 (GnRH-R1, GnRH-R1SHS)].

Phylogenetic analysis of GnRH-R
We used phylogenetic analysis to situate A. burtoni GnRH-Rs with respect to previously cloned GnRH-Rs. All publicly available (February 14, 2007) vertebrate full-length GnRH-Rs were obtained from GenBank. Multiple GnRH-R sequences were aligned using the predicted sequences of GnRH-R polypeptides (ClustalW program in MEGA3.1) (39). This multiple sequence alignment was used to generate a neighbor-joining tree with bootstrap values using the Jones-Taylor-Thornton matrix. Lamprey GnRH-R was used to root the phylogenetic tree. Full species names and GenBank accession numbers for the receptor cDNAs are listed in the supplemental material (published as supplemental data on The Endocrine Society’s Journals Online web site at http://endo.endojournals.org).

Functional assay of GnRH-R1SHS
For functional tests, the A. burtoni GnRH-R1SHS cDNA expression plasmid was transfected into COS-1 cells using diethylaminoethyl dextran, as previously described (40). The ligand structure-function relationships of the putative GnRH-R1SHS were directly compared with A. burtoni GnRH-R2PEY (32) by measuring ligand-stimulated inositol phosphate (IP) production, as described previously (40), except that IP was extracted with formic acid (10 mM). Ligand concentrations ranged from 10–10 to 10–5 M. Data points were determined in duplicate, and EC50 values are the means of three separate experiments.

For competition binding assays, [His5,D-Tyr6] GnRH1 was radioiodinated as previously described (41). Transfected cells were washed twice with cold PBS and incubated (5 h at 4 C) with 125I-labeled [His5,D-Tyr6] GnRH1 (105 cpm) in the absence (B0) or presence of various concentrations (10–9 to 10–5 M) of unlabeled peptides in 0.5 ml HEPES-buffered DMEM with BSA (0.1%). Unbound 125I-labeled [His5,D-Tyr6]-GnRH was removed by washing twice with PBS, and cell-associated radioactivity was solubilized in 1 ml 1 M NaOH and counted. Nonspecific binding was determined in the presence of excess GnRH2 (10–5 M). Data points were determined in duplicate, and the IC50 values are the means of three separate experiments.

Localization of GnRH-R1SHS mRNA transcripts
GnRH-R1SHS expression was localized using PCR on tissue samples from adult A. burtoni (brain, pituitary gland, retina, gill, heart, intestine, kidney, liver, muscle, ovary, retina, spleen, and testes) and in situ hybridization using antisense receptor-specific probes in the brain and pituitary.

In situ hybridization
In situ hybridization methods are described in detail elsewhere (33), so only essential details are given here. Reproductively active males (n = 3) were killed by rapid cervical transection, and brains were removed and placed in sterile PBS. Tissue was frozen in mounting medium (OCT; Tissue-Tek, Torrance, CA) on dry ice, sectioned coronally, and stored at –80 C until use. Probes for GnRH-R1SHS and GnRH-R2PEY were synthesized, based on GenBank sequences (GnRH-R1SHS, GenBank AY705931; GnRH-R2PEY, AY028476), and plasmids were linearized to generate both sense and antisense templates for each receptor. GnRH-R1SHS and GnRH-R2PEY probes were from 151-1035 bp and from 55–1129 bp, respectively, and were generated by transcribing with SP6 or T7 polymerase in the presence of [35S]UTP (Maxiscript kit; Ambion, Austin, TX). Slides with brain tissue were treated as described previously (33). Photomicrographs were captured digitally (SPOT camera system; Diagnostic Instruments, Sterling Heights, MI) and optimized for clarity (PhotoShop; Adobe, San Jose, CA).

PCR on tissue samples
Tissue from T males (~6 cm; ovaries from sexually mature female ~6 cm standard length) was homogenized and RNA was extracted (RNeasy Micro Kit; QIAGEN Inc., Valencia, CA). cDNA was synthesized from total RNA from each tissue sample (3' RACE, SMART kit; Clontech Laboratories Inc., Palo Alto, CA). Sense (5'-attgttctgtctttcatcatctgctgga-3') and antisense (5'-tgctctcagcacgactctcgt-3') primers specific to A. burtoni GnRH-R1SHS generated a 420-bp PCR product, whereas sense (5'-ggctgctcagttccgagtt-3') and antisense (5'-gtgaggaccctctctggtggacatt-3') primers specific to the A. burtoni GnRH-R2PEY generated a 961-bp PCR product from the cDNA template. The GnRH-R2PEY primer pair spans two introns predicted from the conserved GnRH-R genomic structure in fish (12). Positive controls for primers were performed using brain cDNA because in situ hybridization assays have shown both GnRH-R1SHS and GnRH-R2PEY expression in the brain (33). Negative controls were performed using the same procedure as for the experimental group without adding cDNA from any tissue. PCR was performed with a 68–60 C touchdown protocol as follows: 3 min denaturation at 95 C, followed by 16 cycles of 30 sec denaturation at 95 C, 30 sec annealing (68–60 C), and 15 min extension at 72 C. These reactions yielded a single product, as revealed by gel electrophoresis.

Immunocytochemical staining
To localize the somatotropes and gonadotropes, reproductively active fish (T males; n = 2) were killed by rapid cervical transection, and their brains removed and placed in 4% paraformaldehyde in PBS (pH 7.4) overnight and then transferred into a 30% sucrose solution overnight. Fixed tissue was frozen in OCT on dry ice and stored at –80 C. Three series of sections were cut on a cryostat (Microm; Zeiss, Oberkochen, Germany) in coronal or sagittal planes at 14 µm.

Sections were rehydrated (PBS) and incubated in blocking solution (0.3% Triton X-100, 0.2% BSA, 10% normal goat serum diluted in PBS). Slides were incubated separately in three primary antisera [rabbit anti-GH antibody (lot 8502), which labels somatotropes, and rabbit anti-LH antibody (lot 8506) and rabbit anti-FSH antibody (lot 8707), which label gonadotropes in A. burtoni, all generously provided by Dr. A. Takahashi] diluted 1:1000 in blocking solution at 4 C overnight (42). The other set of slides was incubated in blocking solution instead of primary antiserum for the control group. All of the sections were washed in PBS, and the signal was amplified (ABC kit; Vector Laboratories, Burlingame, CA). Sections were incubated in 3',3'-diaminobenzidine (Sigma-Aldrich, St. Louis, MO) to visualize the bound antibody and counterstained with cresyl violet to visualize cell nuclei. The slides were dehydrated in an ethanol/xylene series and coverslipped with Permount medium for microscope viewing.

Regulation of GnRH-R mRNA abundance
Tissue preparation.
To discover whether expression of GnRH-R mRNA was regulated, either as a function of diurnal time or with respect to GnRH ligand expression patterns, animals (n = 24) were housed in separate aquaria with three animals of mixed sex including one T male in each of eight tanks to provide minimal disturbance of other animals in the room when the animals were killed. Tissue was collected from animals in a single tank at 3-h intervals over a 24-h period beginning at 1000 h (fish were maintained on a 12-h light, 12-h dark cycle with dark-to-light transitions at 0900 h). At night, the aquarium room was accessed through a light-tight door, and all tissue-processing steps were performed in complete darkness using minimal infrared illumination. Fish were killed by rapid cervical transection, and their brains and pituitaries were rapidly removed. Brain and pituitary samples were stored at –80 C until use.

Tissue taken at each time point was separately homogenized in chilled TRIzol (Invitrogen) with a Tissue Tearor (Biospec Products, Bartlesville, OK) and kept frozen in TRIzol reagent at –80 C until RNA extraction.

RNA extraction and PCR sample preparation.
Total RNA was extracted from samples following a standard protocol (TRIzol; Invitrogen). RNA was DNase treated to remove genomic DNA contamination (TURBO DNA-free; Ambion), and 1.0 µg total RNA from each tissue sample was reverse transcribed to cDNA (SuperScript II RNase H reverse transcriptase; Invitrogen).

Quantitative real-time PCR of A. burtoni GnRH ligands and receptors.
Primers for quantitative real-time PCR for [Ser8] GnRH1 were 5'-cagacacactgggcaatatg-3' and 5'-ggccacactcgcaaga-3', generating a 128-bp product; for GnRH2, primers were 5'-tggactcctttggcacatcagaga-3' and 5'-ctctggctaaggcatccagaagaa-3', generating a 126-bp product; for GnRH3, 5'-atggatggctaccaggtggaaaga-3' and 5'-tggatttgggcatttgcctcatcg-3', generating a 11- bp product; for GnRH-R1SHS, 5'-tcagtacagcggcgaaag-3' and 5'-gcatctacgggcatcacgat-3', generating a 187-bp product; and for GnRH-R2PEY, 5'-ggctgctcagttccgagtt-3' and 5'-cgcatcaccaccataccact-3', generating a 220-bp product. Three housekeeping genes, actin, 18S rRNA, and tubulin, cloned in A. burtoni were used to control for sample differences in total cDNA content. PCR were performed (iCycler; Bio-Rad, Hercules, CA), and reaction progress in 30-µl reaction volumes was monitored by fluorescence detection at 490 nm during each annealing step. Reactions contained 2x IQ SYBR Green SuperMix (Bio-Rad), 10 µM of each primer, and 1 ng cDNA (RNA equivalent). Reaction conditions were 1 min at 95 C and then 35 cycles of 30 sec at 95 C, 30 sec at 60 C, and 30 sec at 72 C, followed by a melting curve analysis over the temperature range from 95 C to 4 C. All reactions were run in duplicate.

Quantitative real-time PCR data analysis.
Fluorescence readings for each sample were baseline subtracted, and suitable fluorescence thresholds were automatically determined by the MyiQ software. To determine the number of cycles needed to reach threshold, the original fluorescence reading data were analyzed using a special PCR algorithm (43). To calculate the relative mRNA amount, target gene levels were normalized using the geometric mean of housekeeping gene mRNA values.

Statistical analyses
To test the hypothesis that GnRH ligands and/or receptors were expressed with a diurnal rhythm, we analyzed the mRNA samples taken around the clock for rhythmic variation using generalized additive models (44, 45). This technique identifies nonlinear effects that might produce rhythmicity by analyzing the residuals of a scatterplot smoother fit to the data with a back-fitting algorithm. For the generalized additive models and the linear regressions, we took P < 0.05 as the initial significance level and applied a Bonferroni step-down correction to account for family-wise error due to multiple testing. These analyses were computed using R (R:Statistical computing environment; R Foundation, Vienna, Austria).

To identify potential linear relationships between ligands/receptors in a larger data set, we used robust regression (Siegel repeated means) (46) that allowed for the observed nonconstant variance (heteroskedasticity) in the data without requiring data transformation. As a control, we similarly compared the relationships between the ligands/receptors and the Clock 1a gene (GenBank DQ923857). To look for independent sources underlying the data distribution, we used independent component analysis (ICA), which is a recently developed technique (see survey in Ref. 44) related to latent variable and factor analyses. Unlike more traditional approaches, ICA yields unique (except for the sign) solutions to the problem of finding underlying sources of variation by maximizing the information contained in each independent source of interaction (44). The resulting sources may thus be interpreted as latent variables containing information about independent sources of variation in the data. We examined the relationships of each of the ligands/receptors with these latent variables by regressing each of the ligands/receptors on each of the independent sources, again using robust regression. We used a consensus source matrix from 100 solutions arrived at using the maximum-likelihood ICA package (computed using R). It is important to note that because the ICA signs are arbitrary, we can discover whether ligands/receptors relate to sources in the same or opposite directions or not at all (i.e. the signs can be reversed, but not the relative relationships between variables).


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Cloning of full-length GnRH-R1SHS from A. burtoni
Analysis of the full-length putative GnRH-R cDNA sequence from A. burtoni tissue (supplemental Fig. A) shows it to be a second GnRH-R, with conserved Asp residues in transmembrane segments 2 and 7 and the carboxyl-terminal tail that are characteristic of GnRH-Rs that are not mammalian SDP-type GnRH-Rs (47, 48). There are conserved motifs SCAFVT in transmembrane segment 3 and VSHS in EC3 characteristic of teleost receptors (supplemental Fig. A) as identified by Lethimonier et al. (12). This sequence is most similar to other SHS-type GnRH-Rs especially from teleost fish (e.g. GnRH-R1SHS O. niloticus, GenBank AY705931). Consistent with our proposed nomenclature, we named this receptor GnRH-R1SHS (12, 33, 34). The sequence also has residues conserved among all GnRH-Rs (supplemental Fig. A), Asp107, Trp110, Asn111, Lys130, Asn217, and Tyr298, which are equivalent to residues Asp98, Trp101, Asn102, Lys121, Asn212, and Tyr290 of the human GnRH-R1SDP, known to be important in ligand binding (8).

Phylogenetic analysis
Phylogenetic analysis of vertebrate GnRH-Rs identified two families of GnRH-Rs (a and b), each containing two subfamilies, which are distinguished by conserved motifs in EC3 (Fig. 1Go and supplemental Fig. A). The a1 family includes GnRH-Rs from all nonmammalian vertebrate taxa, including fish, reptiles, birds, and amphibians that are characterized by a conserved TPEYVH motif in EC3. All mammalian GnRH-Rs, known to regulate reproduction, form a distinct subfamily a2, characterized by a consensus VSDPVN motif in EC3 and the absence of a C-terminal tail (supplemental Fig. A). With less bootstrap support, family b GnRH-Rs can be split into two subfamilies that are also characterized by different but related motifs in EC3. Teleost and amphibian receptors in subfamily b2 have a VSHSLT motif, whereas GnRH-Rs from amphibians and some mammals belong to subfamily b1, characterized by the VPPSLS motif. It should be noted that although the GnRH-R2PPS from chicken appears to be an exception, its phylogenetic inclusion in subfamily b2 has a very low bootstrap value. Phylogenetic analysis with alternative algorithms (e.g. minimal evolution) yielded very similar topology (data not shown). These and other data (8, 10, 12) showing that GnRH-R families can be distinguished by conserved amino acid motifs in EC3 allow an unambiguous classification system for vertebrate GnRH-Rs. The phylogenetic tree shows that the two types of GnRH-Rs found in teleosts, GnRH-RSHS and GnRH-RPEY, diverged early in vertebrate evolution and are neither recent duplications nor limited to teleosts.


Figure 1
View larger version (14K):
[in this window]
[in a new window]

 
FIG. 1. Phylogenetic relationship among vertebrate GnRH-Rs. A neighbor-joining phylogenetic tree based on the multiple-sequence alignment of vertebrate GnRH receptor amino acid sequences (supplemental Fig. A). The tree was constructed using the Jones-Taylor-Thornton matrix, and bootstrap values are given at each node. The GenBank entry receptor name is given, followed by the EC3 motif. The GnRH-R sequence identified in the lamprey M. marinus was used as an out-group to root the tree. The two main GnRH-R families referred to in the text are labeled a and b.

 
Function of the A. burtoni GnRH-R1SHS
To discover whether the A. burtoni GnRH-R1SHS cDNA encoded a functional receptor, we analyzed its function using a recombinant system. Because fish cell lines were not available, we transfected COS-1 cells with GnRH-R1SHS and showed specific binding of a radiolabeled GnRH analog, which was displaced by mGnRH1 and each of the three forms of GnRH expressed in A. burtoni (Fig. 2Go and Table 1Go). We found that GnRH-R1SHS had relatively low affinity for the three forms of GnRH that are expressed in A. burtoni. IC50 values ranged from 53.5 nM for GnRH2, which had highest affinity for GnRH-R1SHS, to 20,300 nM for [Ser8] GnRH1 (Table 1Go and Fig. 2Go). Consistent with our previous report, GnRH-R2PEY exhibited similar low affinities for all peptides, and peptides had the same rank order of IC50 values in cells expressing either GnRH-R1SHS or GnRH-R2PEY (GnRH2 < GnRH3 < mGnRH1 < [Ser8] GnRH1; Table 1Go and Fig. 2Go).


Figure 2
View larger version (25K):
[in this window]
[in a new window]

 
FIG. 2. Competition binding and IP production of GnRH-R1SHS and GnRH-R2PEY receptors. COS-1 cells transiently transfected with cichlid GnRH-R1SHS (A and C) or cichlid GnRH-R2PEY (B and D) were incubated with 125I-labeled [His5,D-Tyr6] GnRH (A and B) in the presence of increasing concentrations of mGnRH1 ({diamondsuit}), [Ser8] GnRH1 ({blacktriangleup}), GnRH2 (bullet), or GnRH3 ({blacksquare}) for competition binding assays or (C and D) labeled with [3H]myoinositol and incubated with increasing concentrations of peptides in the presence of LiCl for peptide-stimulated IP production.

 

View this table:
[in this window]
[in a new window]

 
TABLE 1. Summary of ligand binding and intracellular signaling of A. burtoni GnRH-R1SHS and GnRH-R2PEY

 
The four GnRH peptides stimulated IP production in cells transfected with A. burtoni GnRH-R1SHS cDNA, confirming that the cDNA encodes a functional receptor that activates the Gq/11 family of G proteins in COS-1 cells and couples to the IP second-messenger cascade. GnRH2 and GnRH3 exhibited the highest potency (EC50 values 8.11 ± 2.63 and 11.0 ± 2.5 nM, respectively), whereas [Ser8] GnRH1 and mGnRH1 had low potency (Table 1Go and Fig. 2Go). In contrast, only GnRH2 exhibited high potency at the GnRH-R2PEY (EC50 0.208 ± 0.034 nM) (Table 1Go and Fig. 2Go).

In vitro analyses show that the A. burtoni GnRH-Rs have similar affinities for each of the GnRH forms in A. burtoni. However, comparison of IC50 and EC50 values indicates that the two receptors differ in their efficiency in coupling binding of different GnRH peptides to intracellular signaling. This is evident from the similar potencies of GnRH2 and GnRH3 (EC50 values 8.11 ± 2.63 and 11.0 ± 2.5 nM) at the GnRH-R1SHS, even though GnRH3 has considerably lower affinity (IC50, 1143 ± 578 nM) at this receptor than GnRH2 (IC50, 53.5 ± 24.6 nM). To quantify these observations, we have calculated an efficiency quotient (IC50/EC50, Table 1Go). Clearly, the GnRH-R1SHS signaling response is particularly sensitive to GnRH3 and poorly responsive to GnRH2, whereas the GnRH-R2PEY is highly responsive to GnRH2 (Table 1Go). It is notable that neither receptor is highly responsive to [Ser8] GnRH1.

Tissue distribution and localization of A. burtoni GnRH-R1SHS and GnRH-R2PEY
We mapped the distribution of both GnRH-Rs in A. burtoni using RT-PCR for peripheral tissue and in situ hybridization for expression in the brain. GnRH-R1SHS is significantly more widespread in the body than is GnRH-R2PEY, being expressed in muscle, intestine, liver, and heart in addition to brain, retina, pituitary, gill, kidney, testes, and ovary where both forms are expressed (Fig. 3Go).


Figure 3
View larger version (19K):
[in this window]
[in a new window]

 
FIG. 3. Tissue distribution of GnRH-Rs. Distribution of mRNA from both GnRH-R types in A. burtoni tissue was assessed using PCR. Note that both receptor types are distributed throughout the body, but GnRH-R1SHS is found more widely than GnRH-R2PEY, particularly in the muscle, intestine, liver, and heart. Neither form is found in spleen. The first lane contains markers for 517 and 396 bp for GnRH-R1SHS and 1018 bp for GnRH-R2PEY. Material for the control lane was reverse transcribed without RNA and then used as the PCR template.

 
To localize GnRH-Rs relative to the gonadotropes and somatotropes, we compared GnRH cell type locations visualized using immunocytochemistry (Fig. 4Go) with each receptor type mapped using in situ hybridization (Fig. 5Go). The organization of the teleost pituitary requires that it be sectioned parasagittally to see these relationships clearly. In sagittal sections of the pituitary, the distributions of cells containing GH (Fig. 4AGo) or LH (Fig. 4BGo) and a schematic representation (Fig. 4DGo) are shown. GH- and LH-containing cells are in nonoverlapping regions of the pituitary. Distribution of FSH-containing cells was similar to that of LH-containing cells (not shown). Figure 4DGo shows the location of the plane of the coronal section through the pituitary used in Fig. 5Go. Based on in situ hybridization, GnRH-R1SHS is localized in a crescent around the ventral border of the pituitary, and GnRH-R2PEY is located dorsally, directly adjacent to the attachment of the pituitary to the brain, confirming previous results for A. burtoni (33). Comparing the in situ hybridization results with those of GH and LH distribution in the pituitary of A. burtoni, it is clear that the distribution of the GnRH-R1SHS correlates with gonadotropes and GnRH-R2PEY with somatotropes (Fig. 6Go).


Figure 4
View larger version (84K):
[in this window]
[in a new window]

 
FIG. 4. Localization of gonadotropes and somatotropes. Fixed sagittal sections of A. burtoni pituitary from T males were incubated with anti-GH antibody (A) anti-LH (B), or no antibody (C). Bound antibodies were visualized with 3',3'-diaminobenzidine, and nuclei were counterstained with cresyl violet. A schematic illustration (D) shows the distribution of GH cells (horizontal shading) and LH cells (vertical shading). The line labeled 5A-5B shows the position of the coronal sections of in situ staining shown in Fig. 5Go. Scale bar, 100 µm. A, Anterior; D, dorsal; P, posterior; V, ventral; H, hypothalamus; RPD, rostral pars distalis; PN, pars nervosa; PI, pars intermedia.

 

Figure 5
View larger version (72K):
[in this window]
[in a new window]

 
FIG. 5. Localization of GnRH-R1SHS and GnRH-R2PEY in the pituitary of A. burtoni using in situ hybridization with 35S-labeled antisense mRNA probes. Photomicrographs show coronal sections through an A. burtoni pituitary showing localization of GnRH-R2PEY mRNA (A) and GnRH-R1SHS mRNA (B). A1, GnRH-R2PEY mRNA as radioactively labeled silver grains (black spots) in bright field with cresyl violet counterstain (purple), showing that it is concentrated in the dorsal part of the pituitary in the proximal pars distalis (PPD) adjacent to the pituitary stalk; A2, dark-field photomicrograph of the same field of view where silver grains appear silver. Arrowheads mark the boundary of the staining that lies between the proximal pars distalis (PPD) and the pars nervosa (PN). B1, GnRH-R1SHS mRNA as silver grains (black spots) in bright field with cresyl violet counterstain (purple), showing that it is concentrated in ventral crescent of the PPD; B2, dark-field photomicrograph of the same field of view where silver grains appear silver. Arrowheads mark the boundary of the staining that lies along the ventral border of the pituitary. Scale bar, 200 µm.

 

Figure 6
View larger version (49K):
[in this window]
[in a new window]

 
FIG. 6. Schematic illustration of a sagittal section through the A. burtoni pituitary showing the distribution of GnRH-R1SHS ({circ}) and GnRH-R2PEY ({triangleup}) receptor types relative to the GH cells (horizontal lines) and LH cells (vertical lines). Scale bar, 100 µm. Abbreviations are as in Fig. 4Go. This figure summarizes a series of sagittal and coronal sections analyzed as described for Figs. 4Go and 5Go.

 
Diurnal rhythm and coregulation of GnRH-R mRNA levels
In A. burtoni, reproductive behavior in T males and activity in non-T males show distinct daily rhythms (28). To determine whether there might be a corresponding rhythm at the GnRH-R mRNA level, we collected tissue from three animals at 3-h intervals over a 24-h period and performed quantitative RT-PCR to determine the levels of receptor expression (Fig. 7Go). Although previous reports suggested GnRH-R mRNA levels change with a diurnal rhythm (49), we found no statistically significant nonlinearity as would be consistent with a diurnal rhythm in either GnRH-R1SHS or GnRH-R2PEY expression.


Figure 7
View larger version (34K):
[in this window]
[in a new window]

 
FIG. 7. GnRH-R mRNA abundance as a function of time of day over a 24-h period. Panels show the fluctuation of relative GnRH-R1SHS (A) and GnRH-R2PEY (B) mRNA abundance over 24 h. A. burtoni tissue samples were taken at 3-h intervals. mRNA abundance was normalized by mRNA levels of housekeeping genes (see Materials and Methods). Gray shaded area represents darkness. Bars indicate SEM. Lights were on at 0900 h and off at 2100 h (n = 24).

 
However, we did find significant coregulation of particular forms of GnRH and GnRH-R mRNA in data collected around the clock from the brain and pituitary. There is a highly significant positive linear relationship between each receptor and all three ligands. ICA showed that two independent sources contribute to changes in all the variables in the same direction. This means that when values in the first underlying source’s distribution increase, expression of all the mRNAs increases. A third underlying source, however, shows GnRH-R2PEY changing opposite of all other variables. Interestingly, the third underlying source shows GnRH3 is expressed in the same direction as receptor 2 but in the opposite direction of GnRH-R1SHS.

Comparison of vertebrate GnRH-R expression and responsiveness to different GnRH forms
The physiological role(s) of GnRH-Rs have been assessed by expression in appropriate cells of the pituitary, measurement of expression after physiological perturbations, and measurement of high affinity for GnRH1. To determine whether data in other species support the identification of GnRH-R1SHS as the reproductive receptor in A. burtoni, despite its low affinity, we reviewed prior research. We asked whether receptor amino acid sequence correlates with ligand responsiveness and with evidence of expression in gonadotropes or responses to perturbation of the HPG axis in nonmammalian vertebrates. We found no correlation of expression of either SHS- or PEY-type GnRH-Rs with reproductive function (supplemental Table B, published as supplemental data on The Endocrine Society’s Journals Online web site at http://endo.endojournals.org).

Because there are limited data on ligand-binding affinities of most vertebrate GnRH-Rs, we compared EC50 values. However, absolute EC50 values measured in recombinant signaling assay systems are influenced by variables such as receptor expression levels and coupling efficiencies that vary among systems. Because all nonmammalian GnRH-Rs are highly responsive to GnRH2, we compared signaling responses stimulated by various GnRH forms relative to the response stimulated by GnRH2 to calculate a ligand response index (Table 2Go). We included all vertebrate GnRH-Rs that fit the four defined receptor types, SDP, PEY, SHS, and PPS, for which we could find dose-response data for GnRH2 and at least one other GnRH form (Table 2Go). The comparison shows that ligand response is similar within each receptor type, and there is a strong correlation of receptor sequence with ligand-stimulated response relative to the response to GnRH2 (see Table 2Go). As expected, the mammalian, SDP-type GnRH-Rs are highly responsive to mGnRH1 relative to GnRH2 (mean EC50 ratio 0.15, Table 2Go), whereas PPS-type receptors, the classic type II GnRH2 receptors, show high responsiveness to GnRH2 (mean ratio 239, Table 2Go). Of the GnRH-R types that occur in teleosts, the PEY-type receptors are highly selective for GnRH2 and very poorly responsive to all other GnRH forms (mean EC50 ratios 140-2900, Table 2Go). It is notable that chicken and bullfrog PEY-type receptors are exceptions, exhibiting good responses to GnRH1. In contrast, teleost SHS-type receptors show smaller differentials in peptide-stimulated responses (mean EC50 ratios 3.3–58, Table 2Go), meaning they are much less selective for GnRH2 and suggesting that a receptor of this type may be the GnRH1 receptor in teleosts.


View this table:
[in this window]
[in a new window]

 
TABLE 2. Comparison of GnRH ligand structure-function data from GnRH-Rs

 

    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Cloning and genomic analysis have identified multiple GnRH-Rs in a wide variety of nonmammalian vertebrates, including chickens and several species of fish and amphibians, yet it is not known which receptor type regulates reproduction. The GnRH-R responsible for reproduction must 1) be expressed in the gonadotrope cells of the pituitary, 2) exhibit enhanced expression in animals that are reproductively competent, 3) be responsive to the endogenous GnRH ligand that regulates the HPG axis, and 4) be phylogenetically close to receptors that regulate reproductive function in related species.

To identify which GnRH-R regulates reproduction in A. burtoni, we cloned the GnRH-R1SHS receptor, known to be up-regulated in reproductively active T males, and measured its phylogenetic and functional characteristics. We found that GnRH-R1SHS is less closely related than the GnRH-R2PEY to the SDP-type receptors, which are clearly established as reproductive receptors in mammals. Interestingly, neither cichlid GnRH-R has high-affinity interaction with [Ser8] GnRH1, but GnRH-RSHS expression colocalizes with gonadotropes, whereas GnRH-RPEY does not.

Phylogenetic relationships of jawed vertebrate GnRH-Rs
Phylogenetic analysis identified four GnRH-R subfamilies with characteristic motifs in EC3. We have defined GnRH-R types according to their EC3 motif, indicated by a three-letter superscript, to reduce confusion of multiple contradictory nomenclature systems for GnRH-Rs (supplemental Table A) and poor correlation of assigned receptor names with receptor type. This unambiguous classification scheme placing GnRH-Rs into four types, SDP, PEY, PPS, and SHS using the EC3 motif, is based solely on receptor amino acid sequence. It has major advantages in that it includes all but a few closely related SHS- and PPS-type receptors (Fig. 1Go and supplemental Fig. A) and cannot be confused with GenBank assigned receptor names.

It is likely that the four GnRH-R subfamilies identified (a1, a2, b1, and b2, corresponding to types PEY, SDP, PPS, and SHS) arose through gene duplication, an important process in evolution. Duplications allow new gene functions to evolve (50, 51), and there is good evidence in yeast that duplicated genes can provide functional compensation against null mutations (52). In this case, GnRH-RSHS and GnRH-RPEY show similar affinity for [Ser8] GnRH1, but they are clearly not functionally redundant, because their expression patterns are divergent. GnRH-R subfamilies can be understood in terms of the two rounds of genome duplication that occurred early in the vertebrate lineage (known as the 2R hypothesis) (53). Because this hypothesis would predict the presence of four GnRH-Rs in all vertebrates, the phylogenetic tree predicts that GnRH-R genes have been lost independently, on multiple occasions from different branches of the vertebrate lineages. For example, SDP-type (family a2) GnRH-Rs have been identified only in mammals and thus must have been lost independently in fish, amphibians, and reptiles. Similarly, PEY-type (family a1) GnRH-Rs are not found in mammals. Amphibians have GnRH-Rs in both subfamilies of family b (PPS-type and P/SQS variations of the SHS type), whereas all fish GnRH-Rs in family b are SHS type and all bird and mammalian receptors are PPS type. Although GnRH-RPPS receptors have been identified in amphibians, chicken, and nonhuman primates, this subtype is absent from the mouse and rat genome sequence data bases. The human GnRH-RIIPPS lacks a start codon, is truncated at the second transmembrane segment, and is nonfunctional (54, 55, 56). Similarly, a sheep GnRH-RIIPPS ortholog is also nonfunctional (57). Thus, this receptor has been independently lost in rodents and has become a pseudogene in humans and sheep.

Teleost fish have undergone an additional genome duplication since their divergence from tetrapods (58, 59), so it is not unexpected that there are multiple GnRH-Rs of a single receptor type in the Ostariophysi lineage including goldfish (GnRH-RAPEY and GnRH-RBPEY) (38), catfish (GnRH-R1PDY and GnRH-R2PEY) (25), and zebrafish (GnRH-R1PEY, GnRH-R3PEY, GnRH-R2SHS, and GnRH-R4SHS) (12). Similarly, in the Neoteleostei lineage multiple GnRH-Rs are present in pufferfish (GnRH-R1/III-1SHS, GnRH-R1/III-2SHA, and GnRH-R1/III-3SHS) and medaka (GnRH-R1SHS and GnRH-R3SHS).

In the case of A. burtoni, the different expression patterns of the GnRH-R genes suggest that these receptors are in the process of diverging in function. Only GnRH-R2PEY is colocalized on all three groups of neurons that produce GnRH, apparently specialized for feedback control of GnRH production (33). Moreover, GnRH-R1SHS colocalizes with gonadotropes, whereas GnRH-R2PEY colocalizes with somatotropes, revealing different roles in the pituitary. When detailed data about receptor localization and function are available from other species, it will become clear whether these specializations are unique to teleosts and whether evolution has taken different directions in different species.

Role of GnRH-R1SHS as a functional receptor for [Ser8] GnRH1
We have previously shown, in A. burtoni, that [Ser8] GnRH1 is the only GnRH peptide found in the pituitary and that it is up-regulated in reproductively active T males. The low affinity of both A. burtoni receptors for [Ser8] GnRH1 indicates that high concentrations of [Ser8] GnRH1 (in the micromolar range) would be required in the pituitary to activate either receptor. Superficially, this might suggest that neither receptor is the target of [Ser8] GnRH1 but that some other receptor with a higher affinity for [Ser8]-GnRH1 may exist. However, all of the nonmammalian GnRH-Rs that have been functionally analyzed have the highest affinity for GnRH2, and all fish receptors are poorly responsive to GnRH1 forms (11, 12, 25, 27, 38) (Table 2Go). Available information from zebrafish and pufferfish genomes does not reveal additional GnRH-Rs, making it unlikely that receptors with high similarity to currently known GnRH receptors and high affinity for GnRH1 exist.

Although [Ser8] GnRH1 had low potency in stimulating IP production in our recombinant system, the efficiency of [Ser8] GnRH1-stimulated activation of fish G proteins is not known. High-affinity GnRH-R binding has been reported in fish pituitary tissue (60, 61, 62, 63) using synthetic, high-affinity ligands, before the discovery of the teleost forms of GnRH1 (see Ref. 12). However, consistent with results in recombinant systems, catfish GnRH1 was found to have low potency in stimulating gonadotropin release and had very low affinity in GnRH-R binding assays performed with catfish pituitary membranes (60). This suggests the known receptors may regulate reproduction despite their low affinity for the biologically relevant ligands, raising the possibility that synchronous release of GnRH in the pituitary may lead to transiently elevated levels triggering the receptors.

Consistent with their known functions in mammals, all SDP-type GnRH-Rs are highly selective for mGnRH1, and all PPS-type receptors are highly selective for GnRH2 (Table 2Go). This suggests that PPS-type GnRH-Rs may have a conserved function in amphibians and in nonhuman primates. As in A. burtoni, all SHS-type GnRH-Rs exhibit relatively low ligand selectivity, with similar responses to different GnRH peptides, suggesting that SHS-type receptors may have similar functions in teleost and amphibian species. In contrast, all PEY-type GnRH-Rs from fish species are highly selective for GnRH2, suggesting that GnRH2 may be their major physiological ligand. However, PEY-type receptors from chicken and bullfrog are much less selective for GnRH2 and are quite responsive to GnRH1 ([Gln8] GnRH1 or mGnRH1, respectively) suggesting that the functions of PEY-type GnRH-Rs may differ between teleosts and other vertebrates. Indeed, there is physiological evidence that the GnRH-R1PEY regulates reproduction in chickens, because its expression is up-regulated in the pituitaries of castrated cockerels (64). However, these data are not consistent with a recent study in which expression of GnRH-R2SHS was shown to increase with reproductive status in male and female chickens (65).

Which GnRH-R regulates reproduction in A. burtoni?
Several lines of evidence suggest that the GnRH-R1SHS receptor regulates reproduction in A. burtoni. First, GnRH-R1SHS receptor mRNA transcripts are up-regulated in the pituitaries of T males, whereas GnRH-R2PEY receptor mRNA transcripts are not (34). Second, GnRH-R1SHS mRNA transcripts specifically localize with expression of gonadotropes, whereas GnRH-R2PEY transcripts localize with somatotropes. Our results using in situ hybridization contrast with results using receptor immunocytochemistry in the closely related teleost O. niloticus (42). In O. niloticus, antibodies that would be predicted to recognize the GnRH-R1SHS (AB111356) reacted with somatotrope cells, whereas antibodies that would be predicted to recognize GnRH-R2PEY (AB111357) reacted with gonadotropes (42). Because both species are cichlids, this difference is likely to result from the different techniques used. The antibodies were raised against short peptides corresponding to the EC3 sequences of GnRH-Rs. These peptide sequences can be readily identified in diverse proteins by a BLAST search, so it is possible that the antibodies recognize the peptide sequences of other proteins besides GnRH-Rs (57). Thus, although comparing GnRH-R expression patterns in different species (supplemental Table B) does not yield a consensus on whether GnRH-RPEY or GnRH-RSHS receptors regulate release of gonadotropins in nonmammalian vertebrates, our additional expression data confirm that GnRH-R1SHS is the receptor expressed in the gonadotropes of A. burtoni.

Coregulation of GnRH ligand and receptor expression
The widespread distribution of GnRH-Rs (66, 67) has led to speculation that different forms of GnRH may play roles beyond reproduction. Chen and Fernald (33) recently showed that the receptor types differ in their brain distribution in A. burtoni and only GnRH-R2PEY colocalized with the three GnRH-producing cell types in the brain, suggesting direct feedback control of the GnRH production. The coregulation of receptors with GnRH suggests that the GnRH system has some common regulatory underpinnings. The presence of multiple GnRH forms across all vertebrates suggests they might share some functional role in reproduction with GnRH1.

GnRH-R evolution and the usefulness of the EC3 motif classification system
The poor correlation of GnRH-R sequence with a role in reproduction across species (supplemental Table B) makes it impossible to predict which GnRH-R will regulate reproductive function in any particular species. The absence of SHS-type GnRH-Rs in higher vertebrates precludes their role as the receptor responsible for reproduction in these species. Different families of duplicated receptor genes regulate reproductive function in mammals (family a, SDP-type) compared with fish (family b, SHS-type). There is insufficient evidence in other taxa (birds, reptiles, and amphibians) to conclude which receptor is used or when selection for one receptor rather than another might have occurred. Thus, it makes sense to define GnRH-R types by genotype, which, perforce, correlates with EC3 amino acid motif and with ligand selectivity, rather than using a classification system based on physiological function in the face of insufficient and conflicting information.

In summary, we cloned the GnRH-R1SHS receptor, which is up-regulated in reproductively active A. burtoni, and showed that it is coexpressed with LH in the pituitary. Like all teleost GnRH-Rs, it binds and responds poorly to GnRH1. Although this remains puzzling, it suggests that synchronous release of GnRH1 into the pituitary allows delivery of sufficient ligand to activate the receptor. The absence of the SHS-type GnRH-R in many taxa and the poor correlation of GnRH-R sequence with a role in reproduction across species, suggest that different families of duplicated GnRH-R genes have been selected to regulate reproduction during vertebrate evolution.


    Acknowledgments
 
We thank M. Morley for technical help and C. Seoighe for insightful comments on the manuscript.


    Footnotes
 
This work was supported by the National Research Foundation (South Africa), the Medical Research Council (South Africa), the University of Cape Town, Lucille P. Markey Biomedical Research Fellowship to C.-C.C., New York State Hatch Grant NYC-165407 to K.E.W., and National Institutes of Health NS34950 Jacob Javits Investigator Award to R.D.F.

Present address for C.A.F.: School of Physiology, University of the Witwatersrand, Parktown, South Africa.

Present address for K.E.W.: Centro de Neurociencia, Universidad de Valparaiso, Valparaiso, Chile.

Disclosure Summary: The authors have nothing to disclose.

First Published Online June 26, 2007

Abbreviations: EC3, Extracellular loop 3; GnRH-R, GnRH receptor; HPG, hypothalamic-pituitary-gonadal; ICA, independent component analysis; IP, inositol phosphate; mGnRH1, mammalian GnRH1; NT, nonterritorial; T, territorial.

Received October 17, 2006.

Accepted for publication June 15, 2007.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 

  1. Marchant TA, Chang JP, Nahorniak CS, Peter RE 1989 Evidence that gonadotropin-releasing hormone also functions as a growth hormone-releasing factor in the goldfish. Endocrinology 124:2509–2518[Abstract/Free Full Text]
  2. Klausen C, Chang JP, Habibi HR 2002 Time- and dose-related effects of gonadotropin-releasing hormone on growth hormone and gonadotropin subunit gene expression in the goldfish pituitary. Can J Physiol Pharmacol Rev 80:915–924[CrossRef]
  3. Weber GM, Powell JF, Park M, Fischer WH, Craig AG, Rivier JE, Nanakorn U, Parhar IS, Ngamvongchon S, Grau EG, Sherwood NM 1997 Evidence that gonadotropin-releasing hormone (GnRH) functions as a prolactin-releasing factor in a teleost fish (Oreochromis mossambicus) and primary structures for three native GnRH molecules. J Endocrinol 155:121–132[Abstract/Free Full Text]
  4. Okubo K, Amano M, Yoshiura Y, Suetake H, Aida K 2000 A novel form of gonadotropin-releasing hormone in the medaka, Oryzias latipes. Biochem Biophys Res Commun 276:298–303[CrossRef][Medline]
  5. Morgan K, Millar RP 2004 Evolution of GnRH ligand precursors and GnRH receptors in protochordate and vertebrate species. Gen Comp Endocrinol 139:191–197[CrossRef][Medline]
  6. White RB, Eisen JA, Kasten TL, Fernald RD 1998 Second gene for gonadotropin-releasing hormone in humans. Proc Natl Acad Sci USA 95:305–309[Abstract/Free Full Text]
  7. Kasten TL, White SA, Norton TT, Bond CT, Adelman JP, Fernald RD 1996 Characterization of two new preproGnRH mRNAs in the tree shrew: first direct evidence for mesencephalic GnRH gene expression in a placental mammal. Gen Comp Endocrinol 104:7–19[CrossRef][Medline]
  8. Millar RP, Lu ZL, Pawson AJ, Flanagan CA, Morgan K, Maudsley SR 2004 Gonadotropin-releasing hormone receptors. Endocr Rev 25:235–275[Abstract/Free Full Text]
  9. Fernald RD, White RB 1999 Gonadotropin-releasing hormone genes: phylogeny, structure, and functions. Front Neuroendocrinol 20:224–240[CrossRef][Medline]
  10. Troskie B, Illing N, Rumbak E, Sun YM, Hapgood J, Sealfon S, Conklin D, Millar R 1998 Identification of three putative GnRH receptor subtypes in vertebrates. Gen Comp Endocrinol 112:296–302[CrossRef][Medline]
  11. Wang L, Bogerd J, Choi HS, Seong JY, Soh JM, Chun SY, Blomenrohr M, Troskie BE, Millar RP, Yu WH, McCann SM, Kwon HB 2001 Three distinct types of GnRH receptor characterized in the bullfrog. Proc Natl Acad Sci USA 98:361–366[Abstract/Free Full Text]
  12. Lethimonier C, Madigou T, Munoz-Cueto JA, Lareyre JJ, Kah O 2004 Evolutionary aspects of GnRHs, GnRH neuronal systems and GnRH receptors in teleost fish. Gen Comp Endocrinol 135:1–16[CrossRef][Medline]
  13. Tsutsumi M, Zhou W, Millar RP, Mellon PL, Roberts JL, Flanagan CA, Dong K, Gillo B, Sealfon SC 1992 Cloning and functional expression of a mouse gonadotropin-releasing hormone receptor. Mol Endocrinol 6:1163–1169[Abstract/Free Full Text]
  14. Stojilkovic SS, Reinhart J, Catt KJ 1994 Gonadotropin-releasing hormone receptors: structure and signal transduction pathways. Endocr Rev 15:462–499[Abstract/Free Full Text]
  15. Ruf F, Fink MY, Sealfon SC 2003 Structure of the GnRH receptor-stimulated signaling network: insights from genomics. Front Neuroendocrinol 24:181–199[CrossRef][Medline]
  16. Karges B, Karges W, de Roux N 2003 Clinical and molecular genetics of the human GnRH receptor. Hum Reprod Update 9:523–530[Abstract/Free Full Text]
  17. Ulloa-Aguirre A, Janovick JA, Leanos-Miranda A, Conn PM 2004 Misrouted cell surface GnRH receptors as a disease aetiology for congenital isolated hypogonadotrophic hypogonadism. Hum Reprod Update 10:177–192[Abstract/Free Full Text]
  18. Flanagan CA, Becker II, Davidson JS, Wakefield IK, Zhou W, Sealfon SC, Millar RP 1994 Glutamate 301 of the mouse gonadotropin-releasing hormone receptor confers specificity for arginine 8 of mammalian gonadotropin-releasing hormone. J Biol Chem 269:22636–22641[Abstract/Free Full Text]
  19. Fromme BJ, Katz AA, Roeske RW, Millar RP, Flanagan CA 2001 Role of aspartate7.32(302) of the human gonadotropin-releasing hormone receptor in stabilizing a high-affinity ligand conformation. Mol Pharmacol 60:1280–1287[Abstract/Free Full Text]
  20. Fromme BJ, Katz AA, Millar RP, Flanagan CA 2004 Pro7.33(303) of the human GnRH receptor regulates selective binding of mammalian GnRH. Mol Cell Endocrinol 219:47–59[CrossRef][Medline]
  21. Wang C, Yun O, Maiti K, Oh da Y, Kim KK, Chae CH, Lee CJ, Seong JY, Kwon HB 2004 Position of Pro and Ser near Glu7.32 in the extracellular loop 3 of mammalian and nonmammalian gonadotropin-releasing hormone (GnRH) receptors is a critical determinant for differential ligand selectivity for mammalian GnRH and chicken GnRH-II. Mol Endocrinol 18:105–116[Abstract/Free Full Text]
  22. Li JH, Choe H, Wang AF, Maiti K, Wang C, Salam A, Chun SY, Lee WK, Kim K, Kwon HB, Seong JY 2005 Extracellular loop 3 (EL3) and EL3-proximal transmembrane helix 7 of the mammalian type I and type II gonadotropin-releasing hormone (GnRH) receptors determine differential ligand selectivity to GnRH-I and GnRH-II. Mol Pharmacol 67:1099–1110[Abstract/Free Full Text]
  23. Millar R, Lowe S, Conklin D, Pawson A, Maudsley S, Troskie B, Ott T, Millar M, Lincoln G, Sellar R, Faurholm B, Scobie G, Kuestner R, Terasawa E, Katz A 2001 A novel mammalian receptor for the evolutionarily conserved type II GnRH. Proc Natl Acad Sci USA 98:9636–9641[Abstract/Free Full Text]
  24. Neill JD, Duck LW, Sellers JC, Musgrove LC 2001 A gonadotropin-releasing hormone (GnRH) receptor specific for GnRH II in primates. Biochem Biophys Res Commun 282:1012–1018[CrossRef][Medline]
  25. Bogerd J, Diepenbroek WB, Hund E, van Oosterhout F, Teves AC, Leurs R, Blomenrohr M 2002 Two gonadotropin-releasing hormone receptors in the African catfish: no differences in ligand selectivity, but differences in tissue distribution. Endocrinology 143:4673–4682[Abstract/Free Full Text]
  26. Okubo K, Nagata S, Ko R, Kataoka H, Yoshiura Y, Mitani H, Kondo M, Naruse K, Shima A, Aida K 2001 Identification and characterization of two distinct GnRH receptor subtypes in a teleost, the medaka Oryzias latipes. Endocrinology 142:4729–4739[Abstract/Free Full Text]
  27. Okubo K, Ishii S, Ishida J, Mitani H, Naruse K, Kondo M, Shima A, Tanaka M, Asakawa S, Shimizu N, Aida K 2003 A novel third gonadotropin-releasing hormone receptor in the medaka Oryzias latipes: evolutionary and functional implications. Gene 314:121–131[CrossRef][Medline]
  28. Fernald RD, Hirata NR 1977 Field study of Haplochromis burtoni: quantitative behavioural observations. Anim Behav 25:964–975[CrossRef]
  29. Fraley NB, Fernald RD 1982 Social control of developmental rate in the African cichlid fish, Haplochromis burtoni. Zeitschrift Tierpsychol 60:66–82
  30. Fernald RD 2003 How does behavior change the brain? Multiple methods to answer old questions. Integr Comp Biol 43:771–779[Abstract/Free Full Text]
  31. Powell JF, Fischer WH, Park M, Craig AG, Rivier JE, White SA, Francis RC, Fernald RD, Licht P, Warby C, Sherwood NM 1995 Primary structure of solitary form of gonadotropin-releasing hormone (GnRH) in cichlid pituitary; three forms of GnRH in brain of cichlid and pumpkinseed fish. Regul Pept 57:43–53[CrossRef][Medline]
  32. Robison RR, White RB, Illing N, Troskie BE, Morley M, Millar RP, Fernald RD 2001 Gonadotropin-releasing hormone receptor in the teleost Haplochromis burtoni: structure, location, and function. Endocrinology 142:1737–1743[Abstract/Free Full Text]
  33. Chen CC, Fernald RD 2006 Distributions of two gonadotropin-releasing hormone receptor types in a cichlid fish suggest functional specialization. J Comp Neurol 495:314–323[CrossRef][Medline]
  34. Au TM, Greenwood AK, Fernald RD 2006 Differential social regulation of two pituitary gonadotropin-releasing hormone receptors. Behav Brain Res 170:342–346[CrossRef][Medline]
  35. White SA, Kasten TL, Bond CT, Adelman JP, Fernald RD 1995 Three gonadotropin-releasing hormone genes in one organism suggest novel roles for an ancient peptide. Proc Natl Acad Sci USA 92:8363–8367[Abstract/Free Full Text]
  36. Miyamoto K, Hasegawa Y, Nomura M, Igarashi M, Kangawa K, Matsuo H 1984 Identification of the second gonadotropin-releasing hormone in chicken hypothalamus: evidence that gonadotropin secretion is probably controlled by two distinct gonadotropin-releasing hormones in avian species. Proc Natl Acad Sci USA 81:3874–3878[Abstract/Free Full Text]
  37. Sherwood N, Eiden L, Brownstein M, Spiess J, Rivier J, Vale W 1983 Characterization of a teleost gonadotropin-releasing hormone. Proc Natl Acad Sci USA 80:2794–2798[Abstract/Free Full Text]
  38. Illing N, Troskie BE, Nahorniak CS, Hapgood JP, Peter RE, Millar RP 1999 Two gonadotropin-releasing hormone receptor subtypes with distinct ligand selectivity and differential distribution in brain and pituitary in the goldfish (Carassius auratus). Proc Natl Acad Sci USA 96:2526–2531[Abstract/Free Full Text]
  39. Kumar S, Tamura K, Nei M 2004 MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform 5:150–163[Abstract/Free Full Text]
  40. Millar RP, Davidson J, Flanagan C, Wakefield I 1995 Ligand binding and second-messenger assays for cloned Gq/G11-coupled neuropeptide receptors: the GnRH receptor. Methods Neurosci 25:145–162[CrossRef]
  41. Flanagan CA, Fromme BJ, Davidson JS, Millar RP 1998 A high affinity gonadotropin-releasing hormone (GnRH) tracer, radioiodinated at position 6, facilitates analysis of mutant GnRH receptors. Endocrinology 139:4115–4119[Abstract/Free Full Text]
  42. Parhar IS, Soga T, Sakuma Y, Millar RP 2002 Spatio-temporal expression of gonadotropin-releasing hormone receptor subtypes in gonadotropes, somatotropes and lactotropes in the cichlid fish. J Neuroendocrinol 14:657–665[CrossRef][Medline]
  43. Zhao S, Fernald RD 2005 Novel algorithm for quantifying real-time polymerase chain reactions. J Comput Biol 12:1047–1064[CrossRef][Medline]
  44. Hyvärinen A, Oja E 2000 Independent component analysis: algorithms and applications. Neural Networks 13:411–430[CrossRef][Medline]
  45. Hastie T, Tibshirani R 1990 Generalized additive models. London: Chapman and Hall
  46. Siegel AF 1982 Robust regression using repeated medians. Biometrika 69:242–244[Abstract/Free Full Text]
  47. Tensen C, Okuzawa K, Blomenrohr M, Rebers F, Leurs R, Bogerd J, Schulz R, Goos H 1997 Distinct efficacies for two endogenous ligands on a single cognate gonadoliberin receptor. Eur J Biochem 243:134–140[Medline]
  48. Blomenrohr M, Bogerd J, Leurs R, Schulz RW, Tensen CP, Zandbergen MA, Goos HJ 1997 Differences in structure-function relations between nonmammalian and mammalian gonadotropin-releasing hormone receptors. Biochem Biophys Res Commun 238:517–522[CrossRef][Medline]
  49. Schirman-Hildesheim TD, Ben-Aroya N, Koch Y 2006 Daily GnRH and GnRH-receptor mRNA expression in the ovariectomized and intact rat. Mol Cell Endocrinol 252:120–125[CrossRef][Medline]
  50. Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J 1999 Preservation of duplicate genes by complementary, degenerative mutations. Genetics 151:1531–1545[Abstract/Free Full Text]
  51. Ohno S 1970 Evolution by gene duplication. Berlin: Springer-Verlag
  52. Gu Z, Steinmetz LM, Gu X, Scharfe C, Davis RW, Li WH 2003 Role of duplicate genes in genetic robustness against null mutations. Nature 421:63–66[CrossRef][Medline]
  53. Holland PW, Garcia-Fernandez J, Williams NA, Sidow A 1994 Gene duplications and the origins of vertebrate development. Dev Suppl 125–133
  54. Faurholm B, Millar RP, Katz AA 2001 The genes encoding the type II gonadotropin-releasing hormone receptor and the ribonucleoprotein RBM8A in humans overlap in two genomic loci. Genomics 78:15–18[CrossRef][Medline]
  55. Conklin DC, Rixon MW, Kuestner RE, Maurer MF, Whitmore TE, Millar RP 2000 Cloning and gene expression of a novel human ribonucleoprotein. Biochim Biophys Acta 1492:465–469[Medline]
  56. Pawson AJ, Morgan K, Maudsley SR, Millar RP 2003 Type II gonadotrophin-releasing hormone (GnRH-II) in reproductive biology. Reproduction 126:271–278[Abstract]
  57. Gault PM, Morgan K, Pawson AJ, Millar RP, Lincoln GA 2004 Sheep exhibit novel variations in the organization of the mammalian type II gonadotropin-releasing hormone receptor gene. Endocrinology 145:2362–2374[Abstract/Free Full Text]
  58. Amores A, Force A, Yan YL, Joly L, Amemiya C, Fritz A, Ho RK, Langeland J, Prince V, Wang YL, Westerfield M, Ekker M, Postlethwait JH 1998 Zebrafish hox clusters and vertebrate genome evolution. Science 282:1711–1714[Abstract/Free Full Text]
  59. Brunet FG, Crollius HR, Paris M, Aury JM, Gibert P, Jaillon O, Laudet V, Robinson-Rechavi M 2006 Gene loss and evolutionary rates following whole-genome duplication in teleost fishes. Mol Biol Evol 23:1808–1816[Abstract/Free Full Text]
  60. Schulz RW, Bosma PT, Zandbergen MA, Van der Sanden MC, Van Dijk W, Peute J, Bogerd J, Goos HJ 1993 Two gonadotropin-releasing hormones in the African catfish, Clarias gariepinus: localization, pituitary receptor binding, and gonadotropin release activity. Endocrinology 133:1569–1577[Abstract/Free Full Text]
  61. Murthy CK, Wong AO, Habibi HR, Rivier JE, Peter RE 1994 Receptor binding of gonadotropin-releasing hormone antagonists that inhibit release of gonadotropin-II and growth hormone in goldfish, Carassius auratus. Biol Reprod 51:349–357[Abstract]
  62. Pagelson G, Zohar Y 1992 Characterization of gonadotropin-releasing hormone binding to pituitary receptors in the gilthead seabream (Sparus aurata). Biol Reprod 47:1004–1008[Abstract]
  63. Habibi HR, Marchant TA, Nahorniak CS, Van der Loo H, Peter RE, Rivier JE, Vale WW 1989 Functional relationship between receptor binding and biological activity for analogs of mammalian and salmon gonadotropin-releasing hormones in the pituitary of goldfish (Carassius auratus). Biol Reprod 40:1152–1161[Abstract]
  64. Sun YM, Dunn IC, Baines E, Talbot RT, Illing N, Millar RP, Sharp PJ 2001 Distribution and regulation by oestrogen of fully processed and variant transcripts of gonadotropin releasing hormone I and gonadotropin releasing hormone receptor mRNAs in the male chicken. J Neuroendocrinol 13:37–49[CrossRef][Medline]
  65. Shimizu M, Bedecarrats GY 2006 Identification of a novel pituitary-specific chicken gonadotropin-releasing hormone receptor and its splice variants. Biol Reprod 75:800–808[Abstract/Free Full Text]
  66. Jennes L, Eyigor O, Janovick JA, Conn PM 1997 Brain gonadotropin releasing hormone receptors: localization and regulation. Recent Prog Horm Res 52:475–490; discussion 490–491[Medline]
  67. Kakar SS, Grantham K, Musgrove LC, Devor D, Sellers JC, Neill JD 1994 Rat gonadotropin-releasing hormone (GnRH) receptor: tissue expression and hormonal regulation of its mRNA. Mol Cell Endocrinol 101:151–157[CrossRef][Medline]



This article has been cited by other articles:


Home page
J EndocrinolHome page
N. T Joseph, K. Morgan, R. Sellar, D. McBride, R. P Millar, and I. C Dunn
The chicken type III GnRH receptor homologue is predominantly expressed in the pituitary, and exhibits similar ligand selectivity to the type I receptor
J. Endocrinol., July 1, 2009; 202(1): 179 - 190.
[Abstract] [Full Text] [PDF]


Home page
EndocrinologyHome page
R. D. Fernald
Gonadotropin-Releasing Hormone Receptors: Where Did They Come From?
Endocrinology, June 1, 2009; 150(6): 2507 - 2508.
[Full Text] [PDF]


Home page
EndocrinologyHome page
J. A. Tello and N. M. Sherwood
Amphioxus: Beginning of Vertebrate and End of Invertebrate Type GnRH Receptor Lineage
Endocrinology, June 1, 2009; 150(6): 2847 - 2856.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Flanagan, C. A.
Right arrow Articles by Illing, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Flanagan, C. A.
Right arrow Articles by Illing, N.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Endocrinology Endocrine Reviews J. Clin. End. & Metab.
Molecular Endocrinology Recent Prog. Horm. Res. All Endocrine Journals