|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Biology, University of Victoria, Victoria BC, Canada V8W 3N51 Peptide Biology Laboratory, The Salk Institute, 10010 N. Torrey Pines Road, La Jolla, California 92037, USA
(Correspondence should be addressed to N M Sherwood; Email: nsherwoo{at}uvic.ca)
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
|
Of the ten hormones in humans, PACAP is the most conserved member. The two biologically active forms, PACAP27 and PACAP38, are identical as the latter is a C-terminal extension of PACAP27 (Miyata et al. 1989, 1990). Zebrafish also have PACAP27 and PACAP38, but the gene encoding the peptides has been duplicated producing two forms of each peptide (Fradinger & Sherwood 2000, Wang et al. 2003). PACAP has a wide range of functions as a neuromodulator, neurotropin, smooth muscle relaxant, and a releasor of pituitary hormones (Sherwood et al. 2000, Chow et al. 2003). In fish, PACAP additionally can release growth hormone (Parker et al. 1997, Wong et al. 1998). The receptors that activate these physiological actions are PAC1R, which preferentially binds PACAP with high affinity compared with VIP, and VPAC1R that shows similar potency with PACAP or VIP (Laburthe et al. 2003). Both PAC1R and VPAC1R are identified in humans and zebrafish, but humans also have a VPAC2R that is related in structure and function to VPAC1R. Fish have receptors that resemble VPAC2Rs in structure, but they have not been tested for binding to any peptides.
Although PRP and PACAP share the same precursor, the function of PRP in mammals is unknown. In human, PRP is 52% identical in amino acids to GHRH27, whereas in zebrafish one of the two PRP peptides, known as GHRH-like peptide1 (GHRH-LP1) has a higher identity (59%) with the zebrafish true GHRH27. Also, release of GH has been observed in some fish experiments but not in others after application of GHRH-LP in vivo or in vitro (see Tam et al. 2007 for review). A PRP receptor has not been identified for any mammal but a full-length sequence (PRPR or alternatively GHRH-LPR) in goldfish (Chan et al. 1998) and partial sequence (transmembrane domains 5–7) in zebrafish (Fradinger et al. 2005) have been reported.
Another two neuropeptides, PHM and VIP, are listed in Table 1; they are cleaved from a single precursor resulting in PHM27 and VIP28 in humans with 44% amino acid identity. VIP is a potent peptide that mediates smooth muscle relaxation, stimulates endocrine and exocrine secretion, regulates circadian rhythms and neuromodulates, whereas PHM is thought to be a less potent version of VIP that uses the same receptors, VPAC1R and VPAC2R. The name of PHM in human becomes PHI in other vertebrates, as the C-terminal peptide is isoleucine (I) rather than methionine (M). In zebrafish, a PHI/VIP cDNA has not been cloned. The receptors for VIP are shared with PACAP as noted earlier. However, a receptor for PHM/PHI has not been identified for any mammal and only one PHIR (goldfish) has been characterized to date (Tse et al. 2002).
Mammalian GHRHs act on the pituitary to release GH. A true GHRH cDNA was identified only in 2007 for chicken (Wang et al. 2007), goldfish, and zebrafish (Lee et al. 2007), and shown to be a homolog with the mammalian form. Relevant here, the zebrafish and goldfish GHRH receptors have been shown to be homologs of the mammalian receptor using synteny and response to fish true GHRH in vivo and in vitro (Lee et al. 2007).
In the present experiments, our aim was to use zebrafish to identify a PHIR and determine which signaling pathway is coupled to this receptor. To test the biological activity of the PHIR, we cloned the PHI/VIP cDNA precursor and synthesized the two peptides. The zebrafish peptides of PHI and VIP along with zebrafish PACAP1 and 2 were used to test the signaling path of the zebrafish PHIR. To further expand our understanding of the evolution of the secretin superfamily, we also cloned and expressed GHRH-LP receptor and characterized its signaling properties. Human and zebrafish peptides were tested also on human VPAC2R and GHRHR. Our hypothesis is that the evolution of the secretin superfamilies of hormones and receptors in vertebrates has resulted in the loss or alteration of certain receptors in mammals with subsequent changes in ligand binding. Phylogenetic analysis and synteny maps were used to evaluate the evolution of the receptors.
| Materials and methods |
|---|
|
|
|---|
Wild-type zebrafish (Danio rerio) were purchased from a local pet store and maintained in our laboratory with a 14h light:10h darkness cycle at 28.5 °C. Procedures were approved by the Animal Care Committee at the University of Victoria. The tissues of adult zebrafish were dissected and frozen in liquid nitrogen.
Peptides
The sequences of seven zebrafish peptides used to test receptor activity are shown in Table 2. All peptides were synthesized at the Salk Institute (La Jolla, CA, USA) using a CS-biopeptide synthesizer (model CS536, CS Bio Co. Inc., San Carlos, CA, USA) on a methyl benzhydrylamine resin using Boc-strategy. The peptides were cleaved with hydrofluoric acid, concomitantly deprotected, and then purified as described in Adams et al. (2003). The purity of the peptides was characterized by capillary zone electrophoresis performed on a Beckman P/ACE System 2050 connected to a ChromJet integrator (Spectra Physics, San Jose, CA, USA). Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry of the peptides was measured on an ABI Perseptive DE-STR instrument (PE Applied Biosystems, Foster City, CA, USA). In addition, three human peptides were tested: PHM27 and VIP28 (Bachem, Torrance, CA, USA), and GHRH44 (Salk Institute).
|
Total RNA was extracted from 11 tissues of adult zebrafish: brain, eye, heart, testis, ovary, spleen, kidney, gill, skin, gut, and swimbladder according to the manufacturer's protocol (RNeasy Mini Kit, Qiagen Inc.) The concentration of total RNA was measured and 900 ng total RNA from each sample was reverse transcribed with Oligo (dT) and Superscript II (Invitrogen) according to the manufacturer's instructions.
RACE reactions for zebrafish phir, ghrh-lpr, and vip (PHI/VIP)
Total RNA (500–600 ng) from adult zebrafish brains was used for the synthesis of cDNA with the First Choice RLM-RACE kit (Ambion) following the manufacturer's instructions. The 5' and 3' RACE primers for phir were designed from a published partial sequence for vipr2 (VPAC2R) (Wang et al. 2003) and ghrh-lpr (Fradinger et al. 2005). A partial predicted zebrafish VIP peptide sequence was found by performing a BLAST search with VIP of rat/mouse, guinea pig, chick, frog, trout/bowfin, cod, dogfish in the Vega Multi Blast View program (http://vega.sanger.ac.uk/Multi/blastview?species=D._rerio); this partial predicted zebrafish VIP peptide was identical to trout VIP. Primers for 3' RACE of vip (PHI/VIP) were designed to this predicted nucleotide sequence. Later, 5' RACE primers were designed from our 3'RACE PCR sequence. RACE-PCRs were carried out for 35 cycles at 94 °C for 30 s, 55 °C annealing temperature for 30 s, and 72 °C for 1 min with a 7-min extension at 72 °C on the last cycle. All primers are listed in Table 3. The PCR products were separated by electrophoresis on a 1.5% agarose gel and stained with ethidium bromide. Then, PCR products were selected and isolated with QIAquick Gel Extraction Kit (Qiagen), and ligated into PGEM-T vector (Promega Corp.)
|
The NCBI database was used to search for PAC1R, VPAC1R, VPAC2R, GHRHR and GHRH-LPR in human, chicken, frog, goldfish, fugu, and zebrafish genomes. The amino acid sequence of the open reading frame of zebrafish PHI receptor with its seven transmembrane domains was used separately to search for homologous receptors in the same genomes as above and in those of medaka (Oryzias latipes), stickleback (Gasterosteus aculeatus), Tetraodon nigroviridis, and salmonids (Salmonidae family) not only in NCBI website (www.ncbi.nlm.nih.gov) but also in Ensembl (www.ensembl.org) and the Salmon Genome Project database (www.salmongenome.no). The deduced amino acid sequences were aligned using the ClustalW alignment program (Thompson et al. 1994) with BLOSUM62. A neighbor-joining (NJ) phylogenetic tree was generated based on the ClustalW alignment of the amino acids from transmembrane domains 1-7 by MEGA3.0 (Kumar et al. 2004) to produce the tree using JTT (Jones–Taylor–Thornton) substitution. In addition, a maximum-likelihood (ML) tree was generated with 100 bootstraps and the JTT substitution using PHYML (Guindon & Gascuel 2003) and its online web server (http://atgc.lirmm.fr/phyml/menu.html; Guindon et al. 2005).
The GenBank accession numbers for the receptors isolated here are zebrafish PHIR (EU150381) and zebrafish GHRH-LPR (EU150382). For other sequences, the numbers are: human (Homo sapiens) GHRHR (NP000814); chicken GHRHR (ABB84385); goldfish GHRHR (ABJ55978); zebrafish GHRHR (ABJ55981); chicken PRP receptor (PRPR; CAC82589); goldfish GHRH-LPR (AAC15698); fugu GHRH-LPR (CAC82589); human VPAC2R (P41587); chicken VPAC2R (NP_001014970); fugu VPAC2RA (AJ408877); fugu VPAC2RB (AJ296143); trout (AY706217), human VPAC1R (P32241); marsh frog VPAC1R (Q9YHC6); goldfish VPAC1R (AAB05459); zebrafish VPAC1R (NP_001013371.1 and for the second receptor ENSDARP00000046126); fugu VPAC1RA (CAC82588); fugu VPAC1RB (CAC82587); human PAC1R (P41586); goldfish PAC1R (O73769); zebrafish PAC1R-short (AAW65134 [GenBank] .1 and for the second receptor XP682980); fugu PAC1RA (CAD35690); and fugu PAC1RB (CAD33842). The Ensembl accession numbers are: medaka VPAC2RA (ENSORLP00000023740), medaka VPAC2RB (ENSORLP0000007394), Tetraodon VPAC1R (GSTENP00016553001), stickleback VPAC2RA (ENSGACP0000002397), and stickleback VPAC2RB (ENSGACP00000023187).
Ensembl also lists two accession numbers for zebrafish VPAC2R, one of which (ENSDARG00000012353) is the full-length (441 amino acids) receptor that we call the PHIR and the other (ENSDARG00000067568) is a truncated (318 amino acids) receptor; the shorter receptor lacks an extracellular N-terminal domain and its 318 amino acids are 100% identical to the longer receptor. The truncated receptor is not considered further. The duplicate receptor for zebrafish VPAC1R was not included in the ML or NJ trees as the receptor had low sequence similarity in TMD 3–6, probably due to inaccurate exon predictions. However, we have included in the ML tree a second zebrafish GHRH-LPR (ENSDARP00000070262), which lacks 100 amino acids in the N-terminal domain, because the remaining amino acids include all 7 TM regions and are distinct from the first form of the receptor; the truncated receptor has not been cloned to date but the genomic sequence shows a 26 amino acid insertion compared with other GHRH-LPRs (Fig. 2). The zebrafish glucagon receptor (GLUCR) (XM_685886) was used as the outgroup for tree construction. Goldfish PHIR was reported earlier but does not have an accession number (Tse et al. 2002). Two receptors with accession numbers are not considered further as second forms in our study, as they are almost identical (99–99.5%) to ones listed in Fig. 1 or 2; one of these receptors is a second form of fugu GHRH-LPR (ENSTRUP00000032622) and the other is a second form of trout VPAC2R (CU069615), which is only a fragment (residues 255–440) of the full-length receptor.
|
|
Tissue-specific distribution of each mRNA was performed by RT-PCR on 11 tissues of zebrafish. Primers for the adcyap1r1-short isoform, vipr1, and phir were designed according to the NCBI Gene Bank database (accession No.: AAW65134.1; NP_001013371.1; NP_571854). ghrh-lpr and vip primers were designed according to our 5' and 3' RACE sequences. A control reaction was prepared using 1 µl of cDNA amplified with zebrafish actin primers. All PCRs were carried out on 1 µl cDNA for 35 cycles at different annealing temperatures: 58 °C (bactin1), 60 °C (adcyap1r1, phir), 62 °C (vip, vipr1, ghrh-lpr, ghrhr). All PCR products spanned at least two exons as a control check for potential genomic contamination. PCR products of the expected size were sequenced and confirmed.
Transfection and expression of zebrafish PHIR, PAC1R-short, VPAC1R, GHRH-LPR and human VPAC2R and GHRHR
Constructs of cDNA full-length coding sequences for receptors PAC1R-short and VPAC1R were a gift from Javier Tello as previously reported (Fradinger et al. 2005). The cDNAs for PHIR and GHRH-LPR containing full-length open reading frames with flanking restriction enzymes sites and a Kozak sequence were subcloned into pcDNA3.1(–) (Invitrogen). Two rounds of PCRs were conducted with platinum Taq polymerase high fidelity (Invitrogen) as outlined by the manufacturer. Primers are listed in Table 3. Human vasoactive intestinal peptide receptor 2 (VPAC2R) and GHRHR cDNA clones were obtained from University of Missouri-Rolla (UMR) cDNA Resource Center (Clone ID: VPAC2R00000, GHRHR00000; GB Acc. No.: NM_003382, AY557192). All six receptors were separately transfected into COS7 cells (American Type Culture Collection, Manassas, VA, USA) as described (Tello et al. 2005) with a few modifications: 1) cells were plated at a density of 60 000 cells/well, 2) medium was replaced after 14 h of transfection with fresh VP-SFM (Invitrogen) for later cAMP assay or replaced with labeling medium (Medium 199, Invitrogen) plus 0.3% bovine albumin (Sigma–Aldrich) for the inositol phosphate (IP) assay, and 3) labeling medium with isobutylmethylxanthine (IBMX) and 10 mM LiCl was used to incubate cells during ligand stimulation. The cAMP assay was done with an Amersham cAMP enzyme immunoassay Biotrak system (GE Healthcare, Little Chalfont Buckinghamshire, UK). In the total IP accumulation assay, quantitation was performed with cell extracts by the multi-well filtration method (Chengalvala et al. 1999).
Data analysis
All IP and cAMP samples were measured in triplicate within each assay and each experiment was independently repeated at least thrice. Data analysis was performed using PRISM3 software (GraphPad Software Inc., San Diego, CA, USA) with nonlinear regression (curve fit) and sigmoidal dose–response curves. The EC50 values (dose of peptide stimulating half-maximal IP or cAMP response) of each peptide were calculated from the means of at least three independent experiments. The data were analyzed by one-way ANOVA followed by Tukey's multiple comparison test. P<0.05 was considered as statistically significant.
| Results |
|---|
|
|
|---|
Full-length cDNAs obtained with the use of total RNA from zebrafish brains and RACE reactions were sequenced in both the forward and reverse directions. The zebrafish PHIR cDNA was 2088 bp in length (GenBank accession No. EU150381) with an open reading frame of 1323 bp (from 169–1491) encoding a 441 amino acid protein (Fig. 1A and B). A Kyte–Doolittle hydrophobicity plot (data not shown) of the deduced protein revealed that it has seven transmembrane domains. Comparison of the zebrafish PHIR deduced amino acid sequence of the full receptor revealed that its highest match is with goldfish PHIR. Also, the PHIR shares high sequence identity with trout VPAC2R (65%), other teleost VPAC2Rs (47–61%), chicken VPAC2R (59%), and human VPAC2R (56%; Table 4). For other members of the secretin receptor family, zebrafish PHIR shares sequence identity with zebrafish PAC1R-short (42%), VPAC1R (42%), GHRH-LPR (37%), and GHRHR (35%). The percent identity plus the functional study suggests that the present cDNA encodes VPAC2R a PHIR, not a VPAC2R. Hereafter, we refer to the receptor as PHIR.
|
The zebrafish GHRH-LPR cDNA is 1747 bp in length (GenBank, accession No. EU150382) with an open reading frame of 1296 bp (from 181 to 1476) encoding a 432 amino acid protein (Fig. 2). A Kyte–Doolittle hydrophobicity plot (data not shown) of the deduced protein was used to identify seven putative transmembrane domains. The amino acid sequence of the zebrafish GHRH-LPR shares the highest level of identity with goldfish GHRH-LPR (88%), followed by fugu (64%), and chicken (43%) as seen in Table 5. On the other hand, zebrafish GHRH-LPR only shares 39–40% sequence identity with GHRHR of zebrafish, goldfish, chicken, and human. The intron and exon organization was not revealed by a BLAST search of the full-length cDNA of GHRH-LPR in the zebrafish genome database (Vega Multi Blast View) as the genome data are not complete. Sequence alignment of GHRH-LPRs from various species (Fig. 2) revealed that zebrafish GHRH-LPR is very similar and structurally related to receptors of other fish. Thirteen conserved cysteines and three putative glycosylation sites, and signature motifs like RLAK, FQGBBVXXBYCFXNXEXQ, and SQR were observed in zebrafish, fugu, and goldfish receptors.
|
3584 bp on chromosome 13 and contains six exons (Fig. 3A and B). The zebrafish mature PHI peptide shares high sequence identity with goldfish PHI A (85%) and B (89%), chicken PHI (67%), and human PHM (78%) peptides (Tse et al. 2002; Fig. 3C). The zebrafish VIP is identical with the mature VIP peptide of trout. Both PHI and VIP in the zebrafish precursor cDNA are flanked at the C-terminus by a GKR, indicating that the peptides are amidated and have potential proteolytic cleavage sites (Fig. 3B).
|
Phylogenetic analysis involved zebrafish PHIR and GHRH-LPR with other members of the secretin receptor superfamily. Analysis of all ClustalW multiple aligned receptors (TMD 1–7) was conducted by the ML method and the NJ method with the GLUCR as the outgroup. The ML tree (Fig. 4) showed a similar result as the NJ tree (Supplementary Figure 1, see Supplementary data in the online of version of the Journal of Endocrinology at http://jme.endocrinology-journals.org/content/vol41/issue/5). The analysis indicates that the five major groups of secretin superfamily receptors included the GHRHR, GHRH-LPR, PAC1R, VPAC1R, and the PHIR/VPAC2R; each group had a bootstrap value of 99–100. The identified zebrafish GHRH-LPR was in a subgroup with like receptors and the potential for grouping with GHRHRs (bootstrap value of 66 in NJ tree but only 46 in ML tree). The PHIR/VPAC2R group had three subgroups: 1) zebrafish and goldfish PHIRs were in a subgroup with five other teleost VPAC2RAs (bootstrap=100), 2) three teleost duplicate VPAC2Rs (receptor B, bootstrap=75) for fugu, medaka, and stickleback were together, and 3) chicken and human VPAC2Rs formed the third subgroup (bootstrap=96).
|
|
To evaluate whether the identified PHIR and PHI peptide are able to transduce a physiological signal, COS7L cells were separately transfected to express zebrafish PHIR, VPAC2R, PAC1R-short, and human VPAC2R. The ability of zebrafish peptides to activate adenylyl cyclase and phospholipase C pathways was assessed by the cAMP and [3H] IP assay (Fig. 6). In COS7L cells expressing the PAC1 receptor, three peptides but not PHI were able to activate the cAMP and IP accumulation in a dose-dependent manner with different half-maximal response (EC50) values (Fig. 6A and B). There was no significant difference (P<0.05) in EC50 values of PAC1R stimulation by the two different PACAP peptides in either the cAMP or IP pathway. But the EC50 values of the PAC1R response to VIP were significantly higher (P<0.05) than those of the PACAP peptides in either the cAMP or IP pathway. Hence, zebrafish PACAP was more potent than zebrafish VIP on the zebrafish PAC1R.
|
|
The physiological characteristics of the zebrafish receptors were compared with human VPAC2R using our zebrafish VIP and human VIP, two peptides with 82% identical amino acids. Zebrafish VIP was able to activate the human VPAC2R in both the cAMP (Fig. 6G) and IP (Fig. 6H) signaling pathways in a dose-dependent manner. It was striking in the cAMP pathway that human VIP and zebrafish VIP produced similar EC50 values with the human VPAC2R (7.6x10–9 M in response to human VIP and 5.4x10–9 M in response to zebrafish VIP; Table 6). Unexpectedly, zebrafish VIP also stimulated the human VPAC2R in the IP signaling pathway (Fig. 6H), as this VPAC2R primarily activates the cAMP path. VIP did not stimulate the zebrafish VPAC1R in the IP path as expected. Hence, zebrafish VIP is a potent stimulus for the human VPAC2R.
Also, the human VPAC2R was stimulated by zebrafish PHI. The EC50 was 7.4x10–9 M, which was not significantly different from the EC50 for VIP suggesting a potent response (Fig. 6G). Part of the effect may be explained by the high sequence identity (78%) between zebrafish PHI and human PHM peptides, although zebrafish PHI was even more potent than human PHM (EC50=4.7x10–8) with the human receptor. In short, zebrafish PHI and VIP are potent stimulators of the human VPAC2R.
Functional analysis of zebrafish GHRH-LPR and human GHRHR
To confirm the functional identity of zebrafish GHRH-LPR, we used three peptides: zebrafish GHRH-LP1, GHRH-LP2, and true GHRH. Each peptide was tested on the zebrafish GHRH-LPR expressed in COS7L cells. For comparison of the physiological profiles, human GHRHR was also expressed in COS7L cells. Both cAMP and IP production were measured after stimulation with different doses of peptides. Zebrafish GHRH-LPR responded to zebrafish GHRH-LP1 (EC50=3.3x10–8 M) and to zebrafish GHRH (1.2x10–7 M) with no significant difference (P>0.05); the responses in the cAMP signaling pathway were dose dependent (Fig. 7) and similar to those in the chicken or goldfish system (Lee et al. 2007, Wang et al. 2007). However, there was no response of zebrafish GHRH-LP receptor to zebrafish GHRH-LP2 stimulation; likewise the goldfish receptor only responded to one of its two GHRH-like peptides (Lee et al. 2007). We found that human GHRHR had no response to either zebrafish GHRH-LP1 or zebrafish GHRH-LP2 but did respond to human GHRH as a positive control.
|
To investigate the expressions of zebrafish PHIR, GHRH-LPR, and PHI/VIP mRNAs in comparison with zebrafish PAC1R, VPAC1R, and GHRHR mRNAs, RT-PCR was used for 11 tissues (Fig. 8). Zebrafish β-actin cDNA (1014 bp) was used as a control to verify the quality of the first-strand cDNA synthesized. The most widely distributed receptor is PHIR mRNA, which is expressed strongly in all tissues. The VPAC1R-short mRNA (arrow in Fig. 8) was also widely expressed in adult zebrafish with stronger expression in the eye, brain, testis, gill, skin, and ovary; lower levels were detected in the heart, kidney, swimbladder, skin, and gut. The longer bands (arrow head in Fig. 8) are from PAC1R-hop mRNA, which has a hop cassette inserted in intracellular loop 3. PAC1R-short mRNA is more widely distributed than other forms of PAC1R mRNA in zebrafish (Fradinger et al. 2005). Also, PAC1R-short is the predominant isoform in goldfish (Kwok et al. 2006). VPAC1R mRNA was expressed in all 11 tissues with a lower level of expression in the gill, kidney, skin, and ovary. PHI/VIP mRNA was more strongly expressed in the eye, heart, brain, testis, gill, gut, and skin than in the swimbladder or ovary. GHRH-LPR and GHRHR mRNAs are distributed to only a few tissues. GHRH-LPR mRNA was only strongly expressed in eye, brain, testis, and kidney with very weak expression in swimbladder and heart. The GHRHR mRNA is mainly distributed in the brain, testis, and gut with very faint bands in the eye, kidney, and swimbladder. β-actin mRNA was expressed in a relatively similar manner in the 11 tissues.
|
| Discussion |
|---|
|
|
|---|
Initially, we thought our newly isolated receptor of 2088 bp was a VPAC2R because it encoded a small section that was identical to a published fragment of VPAC2R (274 bp; Wang et al. 2003). Also, alignment of our deduced amino acid sequence with human and chicken VPAC2R showed a sequence identity of 56 and 59% respectively. Signature motifs characteristic of human and chicken VPAC2R were present in the zebrafish sequence: Thirteen conserved cysteines including seven in the N-terminal extracellular domain, PDV, and RLAK binding motifs and other highly conserved amino acid residues for VIP binding (Fig. 1). Our phylogenetic analysis and chromosome linkage study also grouped this sequence with other VPAC2Rs.
In contrast to the mammalian VPAC2Rs, our receptor was not significantly stimulated by VIP even with a 10–5 M concentration analyzed in both cAMP and IP signaling pathways. Also, in the N-terminal extracellular domain that is fundamental for VIP binding (Laburthe et al. 2003), the zebrafish sequence is missing a potential glycosylation site, which is expressed in the human, chicken, and frog VPAC2R (Hoo et al. 2001). N-glycosylation sites are critical for ligand binding and correct delivery to plasma membrane of the human VPAC1R (Couvineau et al. 1995, 1996). However, the missing glycosylation site is not strong evidence for establishing the identity of our receptor, as the human VPAC1R functions correctly as long as two of the glycosylation sites are intact (Couvineau et al. 1996).
Finally, we discovered that our receptor has 80% amino acid sequence identity with goldfish PHIR; the goldfish sequence had been published (Tse et al. 2002) but not entered into the NCBI database. The goldfish PHIR is the only one characterized in vertebrates so far as we know. The sequence identity of our receptor with goldfish PHI receptor is higher than that for all reported VPAC2Rs (up to 65% in trout). Phylogenetic analysis also grouped the two fish PHIRs together. Moreover, in the N-terminal glycosylation site, goldfish has isoleucine in the same position as zebrafish (Fig. 1B; position 107, open circle) instead of asparagine needed for glycosylation. Based on functional studies (specific response to PHI and no response to VIP or PACAP), sequence identity, and phylogenetic analysis, we conclude that this receptor is zebrafish PHIR. Our strongest evidence comes from functional analysis with PHI in which the receptor was activated in a range expected for the natural ligand (EC50=1.77x10–9 M). Further analysis of fish (e.g. fugu, medaka, and stickleback) VPAC2RA and VPAC2RB for PHIR characteristics is essential. First, to use fugu as an example, its VPAC2Rs only have two glycosylation sites each (Fig. 1B). Second, fugu VPAC2Rs share relatively high sequence identity with zebrafish PHIR: 55 and 62% of VPAC2RA and VPAC2RB respectively. Third, in NJ and ML phylogenetic trees, fugu VPAC2RA and VPAC2RB cluster more closely with zebrafish PHIR than with chicken and human VPAC2R. However, the sequence similarity may not represent functional equivalence as addressed by Cardoso et al. (2005). Whether fugu VPAC2RA and VPAC2RB are PHI receptors will be defined by functional assays.
One hypothesis is that the tetrapod VPAC2R and the teleost PHIR shared a common origin. Support for this hypothesis includes the close structural identity between the two types of receptors (Fig. 1), the chromosomal location and linked genes (Fig. 5), and the activation of the human VPAC2R in which zebrafish PHI (EC50=7.4x10–9 M) is more potent than human PHM (EC50=4.7x10–8 M; Fig. 6E). Common ancestry is also suggested in that zebrafish PHI is a potent stimulator of both the zebrafish PHIR (1.8x10–9 M) and human VPAC2R (7.4x10–9 M; Fig. 6E and H). However, a phylogenetic map for peptide affinity changes cannot be prepared until a number of other so-called VPAC2Rs are tested with a species-specific PHI. This is especially true for the teleost fish where two forms of VPAC2Rs have been annotated in the genome; the duplicates are presumably due to a large scale or whole genome duplication. But the duplicate VPAC2Rs in fish has not yet been cloned and their start sites are ambiguous (Fig. 1). If the duplicates are functional, it is possible that one receptor is a PHI receptor and the other has a different function, possibly as a VPAC2R. Mapping of the receptor changes using affinity values to PHI, VIP, and PACAP should help to elucidate the proposed receptor status.
Zebrafish PAC1R-short has unconventional response to VIP
It was generally accepted that PAC1R specifically binds PACAP with a very low affinity to VIP and that only VPAC1R and VPAC2R have similar affinity to both PACAP and VIP (Vaudry et al. 2000). This is true for PAC1R-short, which is the basic form, and for the PAC1Rs that have added cassettes in the third intracellular loop. However, we show here that the PAC1R-short isoform in zebrafish has a significant dose-dependent response to PACAPs and VIP in both cAMP and IP assays. Although VIP is not as potent as the PACAP (EC50=10–9 M) in activating the PAC1R, the EC50 value of zebrafish VIP is a physiological concentration (EC50=4.7x10–8 M) and there was no significant difference (P<0.05) between PAC1R and VPAC1R in response to VIP in the cAMP path. Support for the VIP action on the PAC1R comes from a recent paper on sea bream in which human VIP activated the sea bream PAC1R (EC50=10–7 M), although less than the PACAPs (EC50=10–9 M; Cardoso et al. 2007b). Relevant also is a report that there are at least 14 variants of the PAC1R in human neuroblastoma cells due to alternative splicing and some of these variants also respond to VIP stimulation in a dose-dependent manner in the cAMP and IP signaling pathways (Lutz et al. 2006). Likewise, human fetal brain has PAC1R variants that respond to VIP (Lutz et al. 2006). It is possible that certain variant PAC1Rs in both fish and human retained the function of responding to VIP. Moreover, human VPAC2R also responds well to zebrafish VIP with no significant difference (P>0.05) compared with zebrafish PAC1R in both the cAMP and IP signaling pathways. The zebrafish VIP peptide shares 82% amino acid identity with human VIP.
Zebrafish GHRH-LPR responds to both zebrafish GHRH-LPR and GHRH
In the present study, we identified a new receptor in zebrafish, which we argue is GHRH-LPR because it has 88% sequence identity with goldfish GHRH-LPR but only 40% identity with zebrafish and human GHRHR (Lee et al. 2007). Sequence alignment with goldfish and fugu GHRH-LPRs reveals closely related structures: 13 conserved cysteines, three N-glycosylation sites, and RLAK for the coupling to Gs and FQGBBVXXBYCFXNXEXQ (B is a hydrophobic residue and X is any residue), characteristic of secretin family receptors. Also phylogenetic analysis grouped zebrafish and goldfish GHRH-LPRs along with other GHRH-LPRs with a bootstrap value of 100. The second form of zebrafish GHRH-LPR also clustered with the same group, as it has seven transmembrane domains but it lacks 100 amino acids at the N-terminus.
Functional assays revealed that the GHRH-LP receptors in both zebrafish and goldfish are most sensitive to one of their two GHRH-LPs but both also respond to their GHRH with about threefold less potency in zebrafish and 30-fold less in goldfish (Lee et al. 2007). The activation of the zebrafish GHRH-LPR was concentration dependent in the cAMP pathway and the potency of zebrafish GHRH-LP1 and zebrafish GHRH was not significantly different (P>0.05). Similarly, Lee et al. (2007) found that the goldfish GHRH-LPR (PRPR) was stimulated by one of two goldfish GHRH-LPs, and by fish GHRH. The EC50 values for fish GHRH tested in goldfish (1.8x10–7 M) and zebrafish (1.15x10–7 M) are nearly the same. Nonetheless, a more potent response occurs with GHRH-LP1; the EC50 values are 5.8x10–9 M for goldfish and 3.41x10–8 M for zebrafish. Moreover, neither zebrafish nor goldfish GHRH-LPR responded to its second form of GHRH-LP. We did not detect any response of human GHRHR to either form of zebrafish GHRH-LP. We conclude that zebrafish GHRH-LPR (PRP) has at least two specific ligands: zebrafish GHRH-LP1 and zebrafish GHRH in the cAMP path. However, zebrafish GHRH receptor only responds to fish GHRH, not to goldfish GHRH-LP stimulation (Lee et al. 2007).
The homologous receptor in humans (PRPR) has been lost from the human genome (Cardoso et al. 2007a, Lee et al. 2007). This implies that the function of GHRH-LP is distinct for fish and other non-mammals compared with mammals.
Cleavage of zebrafish prepro-PHI/VIP is predicted to produce two amidated peptides or extended forms without amidation
PHI and VIP are both transcribed from the same gene in vertebrates. PHI has an amino acid cleavage site of GKR, which implies that PHI is amidated at position 27. Another zebrafish peptide, PHM, is a C-terminal extension of PHI to a methionine at position 42 but is not amidated, as it is followed by only a KR cleavage site. VIP also has a GKR cleavage site after position 28. The zebrafish VIP C-terminus can be extended, which is similar to human VIP. Both PHI and VIP are amidated at the C-terminus but goldfish and chicken PHIs are not amidated because of the lack of a glycine residue (Talbot et al. 1995, Tse et al. 2002). Also, a VIP extension is absent in goldfish because VIP28 is directly followed by a stop codon. The primary structure of zebrafish PHI has been conserved (78% amino acid identity) compared with the human PHM27. In general, PHI/PHM activates all known mammalian VIP receptors (VPACR) with moderate to weak activation (Moriarty et al. 1984, Palle et al. 1989, Lelièvre et al. 1998). However, Lelièvre et al. (1998) suggested that PHI had its own specific receptor that preferentially binds PHI over VIP. In their study, PHI had a 100-fold more potent anti-proliferative action than VIP using radio-labeled PHI as a tracer for binding studies on neuroblastoma cells. The study also revealed that there were high-affinity receptors that selectively bound PHI over VIP, although neuroblastoma cells are known to express many alternatively spliced receptors with different binding profiles (Lutz et al. 2006) The existence of zebrafish and goldfish PHIR and the high conservation of PHI in fish and human suggest that both VIP and PHI play an important role in fish and possibly other vertebrates.
Location of receptors indicates diverse target tissues for peptides
To compare the location of zebrafish PHIR and GHRH-LPR with other related peptide receptors, we performed an RT-PCR study. PHIR is one of the most widely distributed receptors and its mRNA highly expressed in all 11 tissues. The co-expression of PHIR with PHI in the eye, heart, brain, testis, gill, gut, swimbladder, skin, and ovary indicates that PHI may act in an autocrine/paracrine manner. Goldfish PHIR expression was quantified by real-time PCR (Tse et al. 2002). High expression levels were found in the brain, heart, testis, and gut, which are similar to zebrafish.
Zebrafish PAC1R-short mRNA is expressed widely except in the spleen. This result is similar to the observation of Fradinger et al. (2005), except we found additional weak expression in the heart and swimbladder. Goldfish PAC1R-short mRNA was not detected in the spleen either but has strong expression in the heart measured by real-time PCR (Kwok et al. 2006). The PAC1R-short isoform is one of the most widely distributed receptors tested here, along with PHIR in goldfish and zebrafish. In frog and mammal, the distribution was largely investigated by RT-PCR or in situ hybridization in the brain and pituitary (Hashimoto et al. 1996, Shioda et al. 1997, Hu et al. 2000) where PAC1Rs are abundantly expressed.
By contrast, the distribution of zebrafish mRNA for GHRH-LPR and GHRHR are limited compared with other receptors studied here. This probably relates to tissue-specific functions. Zebrafish GHRH mRNA was only strongly expressed in brain, testis, and gut with weak signals in eye, kidney, swimbladder, skin, and ovary. GHRH-LPR mRNA was only detected in eye, brain, testis, and kidney with weak signals in the heart. But GHRH-LPR was expressed strongly in the eye compared with GHRHR. The distribution of zebrafish GHRH-LPR is comparable with a previous study in our laboratory (Fradinger et al. 2005).
| Conclusions |
|---|
|
|
|---|
| Declaration of interest |
|---|
|
|
|---|
| Funding |
|---|
|
|
|---|
| Acknowledgements |
|---|
| References |
|---|
|
|
|---|
Cardoso JC, Power DM, Elgar G & Clark MS 2004 Duplicated receptors for VIP and PACAP (VPAC1R and PAC1R) in a teleost fish, Fugu rubripes. Journal of Molecular Endocrinology 33 411–428.
Cardoso JC, Clark MS, Viera FA, Bridge PD, Gilles A & Power DM 2005 The secretin G-protein-coupled receptor family: teleost receptors. Journal of Molecular Endocrinology 34 753–765.
Cardoso JCR, Vieira FA, Gomes AS & Power DM 2007a PACAP, VIP and their receptors in the metazoa: insights about the origin and evolution of the ligand-receptor pair. Peptides 28 1902–1919.[CrossRef][Web of Science][Medline]
Cardoso JCR, de Vet ECJM, Louro B, Elgar G, Clark MS & Power DM 2007b Persistence of duplicated PAC1 receptors in the teleost, Sparus auratus. BMC Evolutionary Biology 7 221[CrossRef][Medline]
Chan KW, Yu KL, Rivier J & Chow BK 1998 Identification and characterization of a receptor from goldfish specific for a teleost growth hormone-releasing hormone-like peptide. Neuroendocrinology 68 44–56.[CrossRef][Web of Science][Medline]
Chengalvala M, Kostek B & Frail DE 1999 A multi-well filtration assay for quantitation of inositol phosphates in biological samples. Journal of Biochemical and Biophysical Methods 38 163–170.[CrossRef][Web of Science][Medline]
Chow BK, Yuen TT & Chan KW 1997 Molecular evolution of vertebrate VIP receptors and functional characterization of a VIP receptor from goldfish Carassius auratus. General and Comparative Endocrinology 105 176–185.[CrossRef][Web of Science][Medline]
Chow BK, Pang RT & Ng SSMolecular evolution of PACAP precursor and PACAP receptorsH Vaudry & A Arimura Pituitary Adenylate Cyclase-Activating Polypeptide ,2003Norwell, MA:Kluwer Academic Publishers:23–47.
Couvineau A, Gaudin P, Maoret JJ, Rouyer-Fessard C, Nicole P & Laburthe M 1995 Highly conserved aspartate 68, tryptophane 73 and glycine 109 in the N-terminal extracellular domain of the human VIP receptor are essential for its ability to bind VIP. Biochemical and Biophysical Research Communications 206 246–252.[CrossRef][Web of Science][Medline]
Couvineau A, Fabre C, Gaudin P, Maoret JJ & Laburthe M 1996 Mutagenesis of N-glycosylation sites in the human vasoactive intestinal peptide 1 receptor. Evidence that asparagine 58 or 69 is crucial for correct delivery of the receptor to plasma membrane. Biochemistry 35 1745–1752.[CrossRef][Web of Science][Medline]
Foord SM, Bonner TI, Neubig RR, Rosser EM, Pin JP, Davenport AP, Spedding M & Harmar AJ 2005 International Union of Pharmacology. XLVI. G protein-coupled receptor list. Pharmacological Reviews 57 279–288.
Fradinger EA & Sherwood NM 2000 Characterization of the gene encoding both growth hormone-releasing hormone (GRF) and pituitary adenylate cyclase-activating polypeptide (PACAP) in the zebrafish. Molecular and Cellular Endocrinology 165 211–219.[CrossRef][Web of Science][Medline]
Fradinger EA, Tello JA, Rivier JE & Sherwood NM 2005 Characterization of four receptor cDNAs: PAC1, VPAC1, a novel PAC1 and a partial GHRH in zebrafish. Molecular and Cellular Endocrinology 231 49–63.[CrossRef][Web of Science][Medline]
Guindon S & Gascuel O 2003 A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology 52 696–704.
Guindon S, Lethiec F, Duroux P & Gascuel O 2005 PHYML Online – a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Research 33 W557–W559.
Hashimoto H, Nogi H, Mori K, Ohishi H, Shigemoto R, Yamamoto K, Matsuda T, Mizuno N, Nagata S & Baba A 1996 Distribution of the mRNA for a pituitary adenylate cyclase-activating polypeptide receptor in the rat brain: an in situ hybridization study. Journal of Comparative Neurology 371 567–577.[CrossRef][Web of Science][Medline]
Hoo RL, Alexandre D, Chan SM, Anouar Y, Pang RT, Vaudry H & Chow BK 2001 Structural and functional identification of the pituitary adenylate cyclase-activating polypeptide receptor VPAC2 from the frog Rana tigrina rugulosa. Journal of Molecular Endocrinology 27 229–238.[Abstract]
Hu Z, Lelievre V, Tam J, Cheng JW, Fuenzalida G, Zhou X & Waschek JA 2000 Molecular cloning of growth hormone-releasing hormone/pituitary adenylyl cyclase-activating polypeptide in the frog Xenopus laevis: brain distribution and regulation after castration. Endocrinology 141 3366–3376.
Kumar S, Tamura K & Nei M 2004 MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Briefings in Bioinformatics 5 150–163.
Kwok YY, Chu JY, Vaudry H, Yon L, Anouar Y & Chow BK 2006 Cloning and characterization of a PAC1 receptor hop-1 splice variant in goldfish (Carassius auratus). General and Comparative Endocrinology 145 188–196.[CrossRef][Web of Science][Medline]
Laburthe M, Couvineau A & Nicole PMolecular pharmacology and structure-function analysis of PACAP/VIP receptorsH Vaudry & A Arimura Pituitary Adenylate Cyclase-Activating Polypeptide ,2003Norwell, MA:Kluwer Academic Publishers:69–93.
Lee LT, Siu FK, Tam JK, Lau IT, Wong AO, Lin MC, Vaudry H & Chow BK 2007 Discovery of growth hormone – releasing hormones and receptors in nonmammalian vertebrates. PNAS 104 2133–2138.
Lelièvre V, Pineau N, Du J, Wen CH, Nguyen T, Janet T, Muller JM & Wascheck JA 1998 Differential effects of peptide histidine isoleucine (PHI) and related peptides on stimulation and suppression of neuroblastoma cell proliferation. Journal of Biological Chemistry 273 19685–19690.
Lok S, Kuijper JL, Jelinek LJ, Kramer JM, Whitmore TE, Sprecher CA, Mathewes S, Grant FJ, Biggs SH, Rosenberg GB et al. 1994 The human glucagon receptor encoding gene: structure, cDNA sequence and chromosomal localization. Gene 40 203–209.[CrossRef]
Lutz EM, Ronaldson E, Shaw P, Johnson MS, Holland PJ & Mitchell R 2006 Characterization of novel splice variants of the PAC1 receptor in human neuroblastoma cells: consequences for signaling by VIP and PACAP. Molecular and Cellular Neurosciences 31 193–209.[CrossRef][Web of Science][Medline]
McGonnell IM & Fowkes RC 2006 Fishing for gene function – endocrine modelling in the zebrafish. Journal of Endocrinology 189 425–439.
Miyata A, Arimura A, Dahl RR, Minamino N, Uehara A, Jiang L, Culler MD & Coy DH 1989 Isolation of a novel 38 residue-hypothalamic polypeptide which stimulates adenylate cyclase in pituitary cells. Biochemical and Biophysical Research Communications 164 567–574.[CrossRef][Web of Science][Medline]
Miyata A, Jiang L, Dahl RD, Kitada C, Kubo K, Fujino M, Minamino N & Arimura A 1990 Isolation of a neuropeptide corresponding to the N-terminal 27 residues of the pituitary adenylate cyclase activating polypeptide with 38 residues (PACAP 38). Biochemical and Biophysical Research Communications 170 643–648.[CrossRef][Web of Science][Medline]
Moriarty KJ, Hegarty JE, Tatemoto K, Mutt V, Christofides ND, Bloom SR & Wood JR 1984 Effect of peptide histidine isoleucine on water and electrolyte transport in the human jejunum. Gut 25 624–628.
Ohkubo S, Kimura C, Ogi K, Okazaki K, Hosoya M, Onda H, Miyata A, Arimura A & Fujino M 1992 Primary structure and characterization of the precursor to human pituitary adenylate cyclase activating polypeptide. DNA and Cell Biology 11 21–30.[Web of Science][Medline]
Okazaki K, Kimura C, Kosaka T, Watanabe T, Ohkubo S, Ogi K, Kitada C, Onda H & Fujino M 1992 Expression of human pituitary adenylate cyclase activating polypeptide (PACAP) cDNA in CHO cells and characterization of the products. FEBS Letters 298 49–56.[CrossRef][Web of Science][Medline]
Palle C, Ottesen B, Jorggensen J & Fahrenkrug J 1989 Peptide histidine methionine and vasoactive intestinal peptide: occurrence and relaxant effect in the human female reproductive tract. Biology of Reproduction 41 1103–1111.[Abstract]
Parker DB, Power ME, Swanson P, Rivier J & Sherwood NM 1997 Exon skipping in the gene encoding pituitary adenylate cyclase-activating polypeptide in salmon alters the expression of two hormones that stimulate growth hormone release. Endocrinology 138 414–423.
Sherwood NM, Krueckl SL & McRory JE 2000 The origin and function of the pituitary adenylate cyclase-activating polypeptide (PACAP)/glucagon superfamily. Endocrine Reviews 21 619–670.
Shioda S, Yada T, Nakajo S, Nakaya K, Nakai Y & Arimura A 1997 Pituitary adenylate cyclase-activating polypeptide (PACAP): a novel regulator of vasopressin-containing neurons. Brain Research 765 81–90.[CrossRef][Web of Science][Medline]
Solano RM, Langer I, Perret J, Vertongen P, Juarranz MG, Robberecht P & Waelbroeck M 2001 Two basic residues of the h-VPAC1 receptor second transmembrane helix are essential for ligand binding and signal transduction. Journal of Biological Chemistry 276 1084–1088.
Talbot RT, Dunn IC, Wilson PW, Sang HM & Sharp PJ 1995 Evidence for alternative splicing of the chicken vasoactive intestinal polypeptide gene transcript. Journal of Molecular Endocrinology 15 81–91.
Tam JKV, Lee LTO & Chow BKC 2007 PACAP-related peptide (PRP)-Molecular evolution and potential functions. Peptides 28 1920–1929.[CrossRef][Web of Science][Medline]
Tello JA, Rivier JE & Sherwood NM 2005 Tunicate gonadotropin-releasing hormone (GnRH) peptides selectively activate Ciona intestinalis GnRH receptors and the green monkey type II GnRH receptor. Endocrinology 146 4061–4073.
Thompson JD, Higgins DG & Gibson TJ 1994 CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Research 22 4673–4680.
Tse DL, Pang RT, Wong AO, Chan SM, Vaudry H & Chow BK 2002 Identification of a potential receptor for both peptide histidine isoleucine and peptide histidine valine. Endocrinology 143 1327–1336.
Ulloa-Aguirre A & Conn PMG protein-coupled receptors and G proteinsPM Conn & R Means Principles of Molecular Regulation ,2000Totowa, NJ:Humana Press:3–25.
Vaudry D, Gonzalez BJ, Basille M, Yon L, Fournier A & Vaudry H 2000 Pituitary adenylate cyclase-activating polypeptide and its receptors: from structure to functions. Pharmacological Reviews 52 269–324.
Wang Y, Wong AO & Ge W 2003 Cloning, regulation of messenger ribonucleic acid expression, and function of a new isoform of pituitary adenylate cyclase-activating polypeptide in the zebrafish ovary. Endocrinology 144 4799–4810.
Wang Y, Li J, Wang CY, Kwok AH & Leung FC 2007 Identification of the endogenous ligands for chicken growth hormone-releasing hormone (GHRH) receptor: evidence for a separate gene encoding GHRH in submammalian vertebrates. Endocrinology 148 2405–2416.
Wong AOL, Leung MY, Shea WLC, Tse LY, Chang JP & Chow BKC 1998 Hypophysiotropic action of pituitary adenylate cyclase-activating polypeptide (PACAP) in the goldfish: immunohistochemical demonstration of PACAP in the pituitary, PACAP stimulation of growth hormone release from pituitary cells, and molecular cloning of pituitary type I PACAP receptor. Endocrinology 139 3465–3479.
Received in final form 15 August 2008
Accepted 29 August 2008
Made available online as an Accepted Preprint 29 August 2008
This article has been cited by other articles:
![]() |
D. Vaudry, A. Falluel-Morel, S. Bourgault, M. Basille, D. Burel, O. Wurtz, A. Fournier, B. K. C. Chow, H. Hashimoto, L. Galas, et al. Pituitary Adenylate Cyclase-Activating Polypeptide and Its Receptors: 20 Years after the Discovery Pharmacol. Rev., September 1, 2009; 61(3): 283 - 357. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |