Persistence of duplicated PAC1 receptors in the teleost, Sparus auratus
© Cardoso et al; licensee BioMed Central Ltd. 2007
Received: 06 June 2007
Accepted: 12 November 2007
Published: 12 November 2007
Duplicated genes are common in vertebrate genomes. Their persistence is assumed to be either a consequence of gain of novel function (neofunctionalisation) or partitioning of the function of the ancestral molecule (sub-functionalisation). Surprisingly few studies have evaluated the extent of such modifications despite the numerous duplicated receptor and ligand genes identified in vertebrate genomes to date. In order to study the importance of function in the maintenance of duplicated genes, sea bream (Sparus auratus) PAC1 receptors, sequence homologues of the mammalian receptor specific for PACAP (Pituitary Adenylate Cyclase-Activating Polypeptide), were studied. These receptors belong to family 2 GPCRs and most of their members are duplicated in teleosts although the reason why both persist in the genome is unknown.
Duplicate sea bream PACAP receptor genes (sbPAC1A and sbPAC1B), members of family 2 GPCRs, were isolated and share 77% amino acid sequence identity. RT-PCR with specific primers for each gene revealed that they have a differential tissue distribution which overlaps with the distribution of the single mammalian receptor. Furthermore, in common with mammals, the teleost genes undergo alternative splicing and a PAC1Ahop1 isoform has been characterised. Duplicated orthologous receptors have also been identified in other teleost genomes and their distribution profile suggests that function may be species specific. Functional analysis of the paralogue sbPAC1s in Cos7 cells revealed that they are strongly stimulated in the presence of mammalian PACAP27 and PACAP38 and far less with VIP (Vasoactive Intestinal Peptide). The sbPAC1 receptors are equally stimulated (LOGEC50 values for maximal cAMP production) in the presence of PACAP27 (-8.74 ± 0.29 M and -9.15 ± 0.21 M, respectively for sbPAC1A and sbPAC1B, P > 0.05) and PACAP38 (-8.54 ± 0.18 M and -8.92 ± 0.24 M, respectively for sbPAC1A and sbPAC1B, P > 0.05). Human VIP was found to stimulate sbPAC1A (-7.23 ± 0.20 M) more strongly than sbPAC1B (-6.57 ± 0.14 M, P < 0.05) and human secretin (SCT), which has not so far been identified in fish genomes, caused negligible stimulation of both receptors.
The existence of functionally divergent duplicate sbPAC1 receptors is in line with previously proposed theories about the origin and maintenance of duplicated genes. Sea bream PAC1 duplicate receptors resemble the typical mammalian PAC1, and PACAP peptides were found to be more effective than VIP in stimulating cAMP production, although sbPAC1A was more responsive for VIP than sbPAC1B. These results together with the highly divergent pattern of tissue distribution suggest that a process involving neofunctionalisation occurred after receptor duplication within the fish lineage and probably accounts for their persistence in the genome. The characterisation of further duplicated receptors and their ligands should provide insights into the evolution and function of novel protein-protein interactions associated with the vertebrate radiation.
Increased gene number and complexity are generally assumed to have contributed to the success of vertebrates. The evolutionary driving forces behind this are still under debate however gene and/or genome duplications and exon shuffling events are proposed to have been of fundamental importance [1–5]. The increased complexity of metazoan genomes have been attributed to rounds of gene or whole genome duplication [1, 6–9]. Analysis of metazoa genomes reveals that a remarkable percentage of duplicated genes exist [10–13] and whilst some genes decay to non-functionality and are subsequently eliminated from the genome, others are maintained either through the acquisition of novel functions (neofunctionalisation) or by partitioning the function of the ancestral molecule between the duplicated isoforms (subfunctionalisation).
The secretin family of G-protein coupled receptors (GPCRs) (a.k.a. family 2 GPCRs) is a large hormone and neuropeptide receptor gene family present in metazoan genomes. Members of this family have been identified in both protostomes and deuterostomes [14–16] and their conserved sequence and gene organisation has led to the proposal that they evolved from a common ancestral gene as a consequence of total, or partial genome duplication . Vasoactive Intestinal Peptide (VIP) and Pituitary Adenylate Cyclase Activating Polypeptide (PACAP) receptors (VPAC and PAC1, respectively) are closely related members of family 2 GPCRs. They are important pharmaceutical targets as their ligands, the brain-gut peptides VIP and PACAP, control a number of important physiological functions in mammals [17, 18].
In humans three receptors exist, PAC1, VPAC1 and VPAC2 and binding studies reveal that VPACs are able to bind the ligands, VIP and PACAP with similar affinities, while PAC1 preferentially binds PACAP [18, 19]. In vertebrates, activation of PAC1/VPAC receptors occurs via three intracellular transduction mechanisms. These involve either cyclic AMP (cAMP) production, IP turnover via PLC or/and calcium mobilization as a consequence of the activation of a G-protein complex [18, 20–22]. In mammals, five PAC1 isoforms, which result from the insertion of one or two, 28 (hip or hop1 variant) or 27 (hop2 variant) amino acid cassettes in intracellular loop 3 (IL3) have been identified [23, 24]. In addition, two VPAC1 and VPAC2 isoforms have been recently described which lack TM6 and TM7 . In other vertebrates PAC1 splice isoforms have been isolated, however no alternative splice forms of VPAC have yet been identified and in vertebrates these latter receptors do not activate the IP3 signalling pathway [26–28].
Teleosts, which diverged from the tetrapod lineage approximately 450 million years ago (MYA), represent one of the most successful vertebrate groups with over 25,000 species. The existence of a variety of teleost genome sizes and ploidy levels has made them very useful for studies of gene evolution and function . Recently duplicated PAC1 and VPAC receptor genes have been identified in teleosts using in silico approaches  and in Takifugu they have a differential tissue distribution . In the present study duplicate PAC1 cDNAs were isolated from the marine teleost, sea bream (Sparus auratus) and their tissue expression and functional profile characterised using the peptides VIP, PACAP27, PACAP38 and SCT (secretin). The persistence of duplicate sbPAC1 receptors in teleost genomes is discussed in relation to the current proposed theories for gene evolution in vertebrates.
Sea bream duplicate PAC1 receptors
A microsatellite sequence is present in the 5'UTR region of both receptors and sbPAC1A contains an imperfect (AC)32 dinucleotide repeat and sbPAC1B an imperfect (GA)28 repeat upstream of the initiation codon (Figures 1 and 2). Genotyping of the microsatellites using genomic DNA from sea bream caught at different geographic locations and a family panel, revealed both microsatellites are polymorphic and the locus sbPAC1B scores more alleles (8 alleles) than sbPAC1A (4 alleles) (Additional file 1).
Sequence comparisons and phylogenetic analysis
Accession numbers of the putative teleost PAC1, VPAC, PRPR and GHRHR identified in silico.
Tissue expression of sbPAC1 receptors
Functional studies of duplicate sbPAC1
The sbPAC1A and sbPAC1B were successfully expressed in Cos7 and Hek293 cell lines as confirmed by immunofluorescence and Western blot analysis. Cos7 cell extracts subject to Western blot contained a specific immunoreactive fusion protein of approximately 52 kDa in transfected cells, which corresponds to the estimated molecular weight in silico of the fusion proteins, (T7PAC1A is 52.75 kDa and T7PAC1B is 52.32 kDa). Cos7 and Hek293 cell lines expressing the recombinant vector were activated by Forskolin and gave maximal cAMP production. Negative control experiments in which Cos7 and Hek293 cells were transfected with pcDNA3 without insert and incubated with the maximum concentration used of test peptide (10-6 M), revealed that Hek293 cells are responsive to PACAP and VIP peptides. This suggests the existence of endogenous receptors in Hek293 which is confirmed by available proteome data . For this reason only the results obtained with transfected Cos7 cells are presented.
In the present study, duplicate sea bream PAC1 genes (sbPAC1A and sbPAC1B) have been isolated and functionally characterised. Structural motif identification and expression assays reveal they are functional family 2 GPCR members. The sbPAC1B is predominantly found in brain and pituitary while sbPAC1A has a widespread distribution but is absent from brain. The overall tissue distribution of the sbPAC1 receptors is similar to the single mammalian PAC1 receptor and suggests they may have similar functions in the endocrine, nervous, gastrointestinal and reproductive systems [17, 18]. In common with the mammalian PAC1 receptor both sbPAC1 receptors are highly stimulated by PACAP but poorly by VIP and SCT. The tissue distribution of PAC1 receptors in sea bream is different from that in Takifugu in which PAC1B has a widespread tissue distribution and PAC1A distribution is restricted to brain, gill and gonads . Therefore, despite the high sequence conservation between teleost PAC1A or PAC1B genes (91% and 86%, respectively), the paralogue genes seems to have evolved in a species specific manner in relation to tissue distribution and possibly function.
Duplicate PAC1 genes have been identified in other teleosts and also for other family 2 GPCR members [14, 15] as expected in light of the proposed partial or full genome duplication event suggested to have occurred within the teleost lineage [32, 33, 35–37].
Analysis of teleost genomes indicates that duplicate PAC1 genes are linked with the recently reclassified GHRH-like receptor which based upon ligand binding characteristics has been reassigned as a teleost PRP receptor (PRPR) [38, 39]. In Takifugu, PAC1A and PRPR genes are localised on scaffold N000080 and PAC1B/PRPR in scaffold N002399 and in the zebrafish genome on chromosome 10 and 2, respectively [40, 41]. In terrestrial vertebrates, with the exception of mammals, co-localization of both receptors is also observed suggesting that PAC1 and PRPR genes arose by tandem gene duplication prior to the teleost divergence and that PRPR gene was subsequently lost in the mammalian lineage (Figure 6).
Stimulation of cAMP production and not peptide affinity or activation of alternative signalling pathways has been investigated and revealed that sbPAC1 receptors are highly stimulated and have identical potency profiles for the mammalian PACAP27 and PACAP38 peptides compared to VIP. The latter peptide was found to be more potent for sbPAC1A in comparison to sbPAC1B. Although, the identification in teleosts of duplicate genes for the ligands and the existence of two potentially active PACAP peptides raises further issues in relation to PAC1 receptor activation and function. Two copies of PACAP are present in Takifugu (DQ659331 and DQ659332), Tetraodon (Q4RN19 and Q4RH43) and zebrafish (NW_652622 and NW_634478), but functional studies with the duplicate receptors and their ligands are scarce. The ligand binding characteristics of the duplicated zebrafish zfPACAP27 peptides for the zebrafish PAC1 receptor, the sequence homologue of sbPAC1A, was tested and both peptides strongly stimulate in a similar way the cAMP and phospholipase pathways [28, 42] suggesting conservation of function. A further measure of complexity is also introduced by the identification of duplicate VIP genes in teleost genomes [38, 43] and it will be of importance in future to establish the affinity of the duplicate peptides for the duplicate receptors, as well as compare their tissue distribution in teleosts.
Moreover, a microsatellite is identified for the first time in the 5'UTR of PAC1 receptors. Genotyping analysis of the microsatellite reveals it is polymorphic in sea bream and raises intriguing questions about its potential influence on gene expression. Analysis of homologue regions in other teleosts  also reveals the presence of a microsatellite in the 5'UTR. For example, in the stickleback PAC1A and PAC1B respectively, two perfect microsatellite repeats (TG)21 and (CA)34 and an imperfect (CA)21 dinucleotide repeat are present. The PAC1A gene in medaka contains a (TG)8 repeat, although no microsatellite is present in the paralogue gene (Additional file 3). So far no microsatellites have been described or identified in the mammalian, chicken, frog and goldfish homologue receptors [28, 44–46]. The importance of PACAP in growth and development [17, 18] and the presence of a microsatellite in its receptor suggests it may be a useful tool for genomic and functional analysis. A previous study of early growth variation in the Artic charr (Salvelinus alpinus, ) revealed a strong marker-trait association with a single nucleotide polymorphism (G/A) of the 18th base pair of the intronic region between the exons that code for PACAP-related peptide and PACAP in the PRP/PACAP gene precursor.
The human SCT peptide failed to significantly stimulate cAMP production by the sbPAC1 receptor which may indicate either a failure to activate the cAMP signalling pathway or a failure to bind the receptor. The preceding results are at odds with the physiological role attributed to secretin in fish, pancreatic stimulation via cholecystokinin and oxyntomodulin . However, SCT has only been identified by immunohistochemistry in the gastrointestinal tract of the flower fish (Pseudophoxinus antalyae) . The failure to identify in teleosts a gene for the SCT receptor or its ligand suggests they probably evolved subsequent to teleost divergence or were lost during the teleost radiation and the significant sequence differences between SCT and VIP or PACAP probably explains the absence of activity of SCT in the present study .
Gene duplicates and generation of alternative splice isoforms are major contributors to functional diversity of the vertebrate proteome. In teleosts, the existence of gene duplicates is proposed to reduce the incidence of gene isoforms since they are assumed to generate functional redundancy of single-copy gene splice isoform , and this may explain reduced number of PAC receptor splice isoforms detected in the sea bream. In sea bream, a PAC1Ahop receptor isoform was identified with an overlapping tissue distribution with the shorter PAC1A receptor transcript, but it was different from the tissue distribution reported for alternative splice isoforms in zebrafish and goldfish [26, 28]. Recently in zebrafish, two novel IL3 splice isoforms were characterised, a hop2 isoform (insertion of 87 bp) present in ovaries and a novel skip isoform (resulting in a truncated protein) in the gills . Characterisation of a recently isolated PAC1 and hop1 receptor isoform in the goldfish reveals that they have similar activation profiles raising question about their functional role . The mounting evidence for the presence of hop1 receptor isoforms in all vertebrates suggests that a homologue transcript was probably also expressed by the PAC1 receptor gene in ancestral metazoa . In order to understand the functional relevance of duplicate receptors and isoforms it will be important to establish if they are co-localised in vivo in order to establish possible interactions.
Expression studies of the sbPAC1 genes demonstrated they are functional family 2 GPCRs and PACAP and VIP peptides were able to stimulate cAMP production in a dose dependent manner. Such assays also indicate that both sea bream receptors are specific PACAP receptors and have different activation profiles for VIP although, as only cAMP production was measured, it remains to be established if a similar activation profile occurs for the alternative IP3 signalling pathway. One hypothesis for the apparent functional divergence and differential expression of sea bream paralogue receptors may be related to functional specialisation related to the existence in teleosts of duplicate ligands which might also explain the significant amino acid changes (48% sequence identity) observed in the N-terminal region. In mammals, important amino acid motifs involved in receptor binding have been identified and mutation studies reveal the importance of the N-terminus . Amino acid motifs involved in ligand binding are conserved between tetrapod and teleost homologue receptors and include the motifs W-D, G-W-S and the following amino acids W, P and P (Figure 3). However, the amino acid motifs that account for potential ligand selectivity of the teleost duplicated receptors are still unknown, future receptor mutation studies should help to clarify this issue.
An intriguing aspect about the sbPAC1 genes is their divergent tissue distribution and this may be the basis of their specific physiological functions and persistence in the genome. Understanding the biological function of receptors is complex as it is not only the availability and concentration of receptors but also of ligands and accessory factors at a given site which will determine receptor preference/activity and ultimately biological function. Physiological studies of the ligands, PACAP and VIP, in teleosts are not very numerous and certainly do not encompass all the actions assigned in mammals making it difficult currently to assign possible biological roles to the duplicate receptors. However, one major function attributed to PACAP is its role in GH-release [46, 51, 52]. In common with mammals, PACAP38 is the predominant isoform in teleost brain and this peptide is found to have a more potent stimulatory effect on fish GH secretion by pituitary cells when compared to GHRH and GnRH (potent mammalian GH releasing factors). In contrast, VIP has little or no effect on GH release [46, 51, 52] and this has been taken to suggest that teleost GH-release is mediated via PAC1 pituitary receptors and the sea bream PAC1B may play a central role in this process. Relatively few studies have been carried out to characterise the biological activity of fish VIP and as in mammals it is proposed to have an important role in the gastrointestinal system [51, 53, 54]. In cod, Gadus morhua, VIP stimulates gastric and pancreatic secretion [55, 56] and in tilapia it is involved in ion and water absorption by the intestine . In sea bream, tissue distribution of PACAP and VIP transcripts indicate that PACAP is mainly restricted to nervous tissue whilst VIP is abundant in the gastrointestinal system but is also present in a wide range of other tissue (unpublished data). The widespread distribution in non nervous tissue of sbPAC1A indicates that this receptor may have a broader physiological role. Unquestionably much more work is required to elucidate the physiological relevance of duplicate sbPAC1 and alternative approaches such as ligand mutational studies; gene knock-down strategies and physiological experiments will be needed.
Duplicate PAC1 receptors genes are present in the majority of teleost genomes, although the reason for their persistence is not yet clearly established. The present study with duplicate sbPAC1 receptors suggests that their maintenance may be due to a process of neofuncionalisation as a consequence of the accumulation of mutations in the ligand binding domain of the receptors after duplication. Such a proposal is supported by the poor sequence conservation (48%) in the N-terminal ligand binding domain of the receptor. The divergent tissue distribution of the receptors, with one form predominantly found in nervous tissue and the other with a more widespread distribution is highly suggestive of functional divergence. The isolation and characterisation of ligands for the receptors in teleosts will be an important step in establishing receptor function, as will improved characterisation of the tissue distribution of both receptors and ligands and the factors regulating their expression.
Sea bream cDNA library screening
Primer sequences used for PCR amplification reactions.
RT-PCR analysis of receptor tissue distribution
Tissue distribution and analysis of PAC1 splice variants was carried out by RT-PCR. Total RNA was extracted from sea bream pituitary, brain, kidney, gills, gut, heart, gonads, liver and skin using TRI reagent (Sigma-Aldrich, Spain). Mature sea bream of 11–13 months old weighing approximately 350–500 g were sacrificed by decapitation and their tissues collected and immediately frozen at -80°C. cDNA was synthesised using 1 μg of sea bream total RNA and a reverse transcription system kit (Promega, Spain) and hexameric oligonucleotides following the manufacture's instructions. The quality of cDNA obtained was verified by PCR with sea bream EF1-α (a housekeeping gene) using the following thermocycle; 94°C for 2 minutes; 25 cycles (94°C for 1 minute, 58°C for 1 minute and 72°C for 1 minute) and a final extension step at 72°C for 5 minutes. Specific primers for each sbPAC1 gene were designed (Table 2) and PCR was performed as previously described using the following cycle: 94°C for 2 minutes; 34 cycles (94°C for 1 minute, 62°C for 1 minute and 72°C for 1 minute) and a final step at 72°C for 5 minutes. The products obtained were sequenced to confirm their identity.
The variability of a microsatellite identified in the 5' UTR of sbPAC1A and sbPAC1B was assessed using specific primers (Table 2) and sea bream genomic DNA from fish of diverse geographic origins (Morocco, Portugal, France and Adriatic sea) and from a genomic panel composed of parents and 50 first generation progeny. PCR reactions were performed using fluorescently labelled forward primers (0.3 μM; 6-FAM and TET; Metabion International AG) and unlabelled reverse primers (0.3 μM) using the following thermocycle: 95°C for 2 minutes; 35 cycles (95°C for 30 seconds, 57°C for 30 seconds and 72°C for 30 seconds) and a final step at 72°C for 5 minutes. PCR products were separated using high resolution 6% Long Ranger acrylamide gels (Cambrex, USA) on a automated ABI 377 sequencer (Applied Biosystems, USA) and data was analysed with the GenScan software (Applied Biosystems, USA).
Database searches, sequence alignments and phylogenetic analysis
The conserved amino acid sequence of sbPAC1 transmembrane (TM) domains was used to search for homologous receptors in teleost genomes. Briefly, the amino acid sequences of the seven TM domains were extracted, concatenated and used in BLAST sequence similarity searches  against the Tetraodon nigroviridis , medaka (Oryzias latipes), stickleback(Gasterosteus aculeatus), zebrafish (Danio renio)  and atlantic salmon (Salmon salar)  genome databases and NCBI EST database . The amino acid sequence of the TM domains of a total of 52 receptors were concatenated and a multiple sequence alignment was produced using the ClustalX vs1.83 (Blosum matrix and Gap opening penalty 10 and Gap extension 0.2)  and percentage of identity/similarity calculated using Genedoc . The alignment produced (length 170, with 162 informative sites) did not require the insertion of gaps and was used to construct phylogenetic trees using both maximum parsimony and neighbour joining methods  with 1000 bootstrap replicates, complete gap deletion and Poisson correction in MEGA 3.1 phylogenetic programme .
Short-range linkage analysis
The gene environment of the Takifugu, Xenopus, chicken and human PAC1 homologue regions were compared using a sequence similarity approach. The human and chicken gene environments were accessed using the NCBI Mapview interfaces  and Xenopus using the Ensembl database . The gene environment of the Takifugu scaffolds (release17/05) was accessed using NIX annotation  and the neighbouring genes were used to search for orthologues in human, chicken and Xenopus genomes using the tblastn algorithm .
Construction of the recombinant expression vector
Specific primers for each sbPAC1 gene were designed to amplify the mature receptor sequence (Table 2) and were cloned into pcDNA3 vector (Invitrogen, UK) containing a signal peptide of CD33 and a T7-epitope tag (CD33-T7-pcDNA3, ). In order to facilitate cloning, restriction digestion sites for BglII and EcoRI enzymes were incorporated in the forward and reverse primers respectively. Template amplification was carried out using the thermocyle; 95°C for 2 minutes; 35 cycles of (95°C for 1 minute, 68°C for 1 minute and 70°C 2 minutes); 72°C for 10 minutes with a proof-reading Pfu DNA polymerase (Promega, Spain) in a 25 μl reaction volume containing 1 × PCR buffer (Promega, Spain), 1,5 mM MgCl2 (Promega, Spain), 0,2 mM dNTPs (Amersham), 1 mM of each primer and 1,5 U Pfu DNA polymerase (3 U/μl)) and DNase Free water (Sigma-Aldrich, Spain). The amplified PCR products were cloned into pGEMT-easy vector (Promega, Spain), sequenced and subcloned into CD33-T7-pcDNA3 vector in frame with the T7-epitope tag. The recombinant constructs produced (CD33-T7-pcDNA3+sbPAC1A and CD33-T7-pcDNA3+sbPAC1B) were sequenced and used to transfect mammalian Cos7 and Hek293 cells lines. The CD86 (NM_175862) membrane protein of the immunoglobulin superfamily cloned in CD33-T7-pcDNA3 expression vector was used as positive control for cell transfection (de Vet et al., 2001).
Mammalian cell transfections
Mammalian Cos7 and Hek293 cell lines were transfected using the Effectene transfection kit (Qiagen, Germany) and the success of transfection was assessed by western blot and immunofluorescence assays using antisera raised against the T7-epitope tag protein . One day prior to transfection, approximately 60,000 to 70,000 Cos7 and Hek293 cells were seeded and transfections carried out using approximately 0.4 μg of recombinant construct (pcDNA3+sbPAC1A or B) and cells were grown for 2 days at 37°C. The viability of transfected cells was determined by dye exclusion using Trypan blue (0,4% solution, Sigma) and the success of transfection calculated by fluorescent confocal microscopy using DAPI (4',6-Diamidino-2-phenylindole) and FITC (Fluorescein) staining for Cos7 cells and the FACS (Fluorescence Activated Cell Sorting) method for Hek293 cells. The percentage of cell transfection was estimated to be approximately 30% in both cell lines.
Immunofluorescence and western blot assays
An anti-T7 epitope monoclonal antibody (Novagen, UK) was used in immunofluorescence and western blots assays to assess the success of transfection and production of recombinant fusion protein (T7+sbPAC1) as described in de Vet et al., 2001 . Briefly, for the immunofluorescence assay, cells were grown in 6 well plates (Greiner, Germany) on sterile cover slips. Approximately 120,000 cells were used per well and cell transfections were carried out using 0.4 μg of DNA as described above. Immunofluorescence localization studies were performed  and DAPI and FITC staining examined using a microscope linked to a confocal imaging system (Bio-Rad, UK). Western blots were carried out by lysing transfected cells (30 μl) in 1× Laemmli SDS-PAGE loading buffer  and proteins were fractionated on a 10% SDS-PAGE gel with a constant current of 35 mA. The fractionated proteins were transferred to a nitrocellulose membrane (GE healthcare, UK) and blocked (2% milk powder in 1 × PBS/0.1%Tween) for 1 hour at room temperature. Incubations were carried out with anti-T7 tag monoclonal antibody (Novagen, UK) and a secondary antibody coupled to Horse-Radish peroxidase as described in de Vet., 2001 . Immunoreactive proteins were detected using the ECL system (PerkinElmer Life Sciences, UK). The size of the immunoreactive protein was comparable to the estimated sizes of the recombinant receptor protein determined in silico using the Swiss Prot database interface .
Ligand binding studies were performed two days after transfections, in three independent experiments. Transfected cells were incubated in duplicate with different concentrations (10-6 to 10-11M) of the human VIP, PACAP38 and secretin (SCT) peptides and ovine PACAP27 peptide for 30 minutes (Sigma-Aldrich, Spain). Homologous sea bream peptides are unavailable but their predicted amino acid sequences are 82%, 100%and 92% identical to human VIP, PACAP27 and PACAP38 respectively and important amino acids at N-terminus potential involved in receptor binding are totally conserved (unpublished data). Prior to ligand-binding assays, Cos7 and Hek293 cells were washed 3 times with DMEM medium without fetal calf serum and incubated in a CO2 incubator with 1 mM of IBMX (3-Isobutyl-1-Methylxantine) for 30 minutes. Ligand-binding assays were performed by incubating the cells in fresh medium containing the peptides in the presence of 1 mM IBMX for 30 minutes at 37°C in the CO2 incubator. After incubation, cells were lysed and stored at -80°C until required. Mammalian Cos7 and Hek293 cells transfected with wild type pcDNA3 were used as negative controls which were incubated with the highest tested peptide concentration (10-6mM). Prior to assays the responsiveness of transfected Cos7 and Hek293 cells was tested by stimulating cAMP production using Forskolin (10 mM, 1 mM and 0.1 mM; Sigma-Aldrich).
Radioimmunoassay (RIA) and statistical analyses
The quantification of cAMP produced was determined by radioimmunoassay using the TRK432 kit (GE Healthcare, UK) following the manufacturer's instructions. Cos7 and Hek293 cells were lysed by sonication, centrifuged and the supernatant heat denatured for 10 minutes at 100°C. The concentration of cAMP produced (pmol/well) was determined in duplicate for each sample and calculations performed based on a linear regression curve constructed using standard concentrations of labelled (3H) cAMP. The cAMP data was normalized as a percentage of stimulation above basal levels and plotted as a percentage of cAMP production per well (%cAMP/well). The results are presented as the mean ± SEM of three independent experiments in duplicate and analysis was performed using the SigmaPlot9.01 programme. Data was analyzed by comparing the potency profile of each sbPAC1 receptor in the presence of different concentrations of test peptides and by comparing the potency profile of both receptors in the presence of identical peptide concentrations. The presence of significant differences in cAMP production was assessed with the SigmaStat3.11 programme using two way Anova and the Holm-Sidak method for pairwise multiple comparisons (P < 0.05 and P < 0.001 were considered statistically significant).
This work was supported by a FCT/CCMAR pluriannual grant. JCRC was funded by FCT grant BPD/14449/03. The authors would like to acknowledge the technical assistance of Elisabete Caldas for performing part of the teleost sequence database searches and to Elsa Couto for helping with the cAMP radioimmunoassays.
- Ohno S: Evolution by gene duplication. 1970, Berlin, New York, Springer-VerlagView ArticleGoogle Scholar
- Iwabe N, Kuma K, Miyata T: Evolution of gene families and relationship with organismal evolution: rapid divergence of tissue-specific genes in the early evolution of chordates. Mol Biol Evol. 1996, 13 (3): 483-493.View ArticlePubMedGoogle Scholar
- Patthy L: Genome evolution and the evolution of exon-shuffling--a review. Gene. 1999, 238 (1): 103-114. 10.1016/S0378-1119(99)00228-0.View ArticlePubMedGoogle Scholar
- Kondrashov FA, Rogozin IB, Wolf YI, Koonin EV: Selection in the evolution of gene duplications. Genome Biol. 2002, 3 (2): RESEARCH0008-10.1186/gb-2002-3-2-research0008.PubMed CentralView ArticlePubMedGoogle Scholar
- Seoighe C: Turning the clock back on ancient genome duplication. Curr Opin Genet Dev. 2003, 13 (6): 636-643. 10.1016/j.gde.2003.10.005.View ArticlePubMedGoogle Scholar
- Ferris SD, Whitt GS: Loss of duplicate gene expression after polyploidisation. Nature. 1977, 265 (5591): 258-260. 10.1038/265258a0.View ArticlePubMedGoogle Scholar
- Holland PW, Garcia-Fernandez J, Williams NA, Sidow A: Gene duplications and the origins of vertebrate development. Dev Suppl. 1994, 125-133.Google Scholar
- Holland PW: Gene duplication: past, present and future. Semin Cell Dev Biol. 1999, 10 (5): 541-547. 10.1006/scdb.1999.0335.View ArticlePubMedGoogle Scholar
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, Nicaud S, Jaffe D, Fisher S, Lutfalla G, Dossat C, Segurens B, Dasilva C, Salanoubat M, Levy M, Boudet N, Castellano S, Anthouard V, Jubin C, Castelli V, Katinka M, Vacherie B, Biemont C, Skalli Z, Cattolico L, Poulain J, De Berardinis V, Cruaud C, Duprat S, Brottier P, Coutanceau JP, Gouzy J, Parra G, Lardier G, Chapple C, McKernan KJ, McEwan P, Bosak S, Kellis M, Volff JN, Guigo R, Zody MC, Mesirov J, Lindblad-Toh K, Birren B, Nusbaum C, Kahn D, Robinson-Rechavi M, Laudet V, Schachter V, Quetier F, Saurin W, Scarpelli C, Wincker P, Lander ES, Weissenbach J, Roest Crollius H: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431 (7011): 946-957. 10.1038/nature03025.View ArticlePubMedGoogle Scholar
- Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations. Genetics. 1999, 151 (4): 1531-1545.PubMed CentralPubMedGoogle Scholar
- Betran E, Long M: Expansion of genome coding regions by acquisition of new genes. Genetica. 2002, 115 (1): 65-80. 10.1023/A:1016024131097.View ArticlePubMedGoogle Scholar
- Amores A, Suzuki T, Yan YL, Pomeroy J, Singer A, Amemiya C, Postlethwait JH: Developmental roles of pufferfish Hox clusters and genome evolution in ray-fin fish. Genome Res. 2004, 14 (1): 1-10. 10.1101/gr.1717804.PubMed CentralView ArticlePubMedGoogle Scholar
- Hughes AL: Adaptive evolution after gene duplication. Trends Genet. 2002, 18 (9): 433-434. 10.1016/S0168-9525(02)02755-5.View ArticlePubMedGoogle Scholar
- Cardoso JC, Pinto VC, Vieira FA, Clark MS, Power DM: Evolution of secretin family GPCR members in the metazoa. BMC Evol Biol. 2006, 6: 108-10.1186/1471-2148-6-108.PubMed CentralView ArticlePubMedGoogle Scholar
- Cardoso JC, Power DM, Elgar G, Clark MS: Duplicated receptors for VIP and PACAP (VPAC1R and PAC1R) in a teleost fish, Fugu rubripes. J Mol Endocrinol. 2004, 33 (2): 411-428. 10.1677/jme.1.01575.View ArticlePubMedGoogle Scholar
- Cardoso JC, Clark MS, Viera FA, Bridge PD, Gilles A, Power DM: The secretin G-protein-coupled receptor family: teleost receptors. J Mol Endocrinol. 2005, 34 (3): 753-765. 10.1677/jme.1.01730.View ArticlePubMedGoogle Scholar
- Sherwood NM, Krueckl SL, McRory JE: The origin and function of the pituitary adenylate cyclase-activating polypeptide (PACAP)/glucagon superfamily. Endocr Rev. 2000, 21 (6): 619-670. 10.1210/er.21.6.619.PubMedGoogle Scholar
- Vaudry D, Gonzalez BJ, Basille M, Yon L, Fournier A, Vaudry H: Pituitary adenylate cyclase-activating polypeptide and its receptors: from structure to functions. Pharmacol Rev. 2000, 52 (2): 269-324.PubMedGoogle Scholar
- Schafer H, Zheng J, Morys-Wortmann C, Folsch UR, Schmidt WE: Structural motifs of pituitary adenylate cyclase-activating polypeptide (PACAP) defining PAC1-receptor selectivity. Regul Pept. 1999, 79 (2-3): 83-92. 10.1016/S0167-0115(98)00147-5.View ArticlePubMedGoogle Scholar
- Laburthe M, Couvineau A, Marie JC: VPAC receptors for VIP and PACAP. Receptors Channels. 2002, 8 (3-4): 137-153. 10.1080/10606820213680.View ArticlePubMedGoogle Scholar
- Couvineau A, Lacapere JJ, Tan YV, Rouyer-Fessard C, Nicole P, Laburthe M: Identification of cytoplasmic domains of hVPAC1 receptor required for activation of adenylyl cyclase. Crucial role of two charged amino acids strictly conserved in class II G protein-coupled receptors. J Biol Chem. 2003, 278 (27): 24759-24766. 10.1074/jbc.M301916200.View ArticlePubMedGoogle Scholar
- Dickson L, Aramori I, McCulloch J, Sharkey J, Finlayson K: A systematic comparison of intracellular cyclic AMP and calcium signalling highlights complexities in human VPAC/PAC receptor pharmacology. Neuropharmacology. 2006, 51 (6): 1086-1098. 10.1016/j.neuropharm.2006.07.017.View ArticlePubMedGoogle Scholar
- Spengler D, Waeber C, Pantaloni C, Holsboer F, Bockaert J, Seeburg PH, Journot L: Differential signal transduction by five splice variants of the PACAP receptor. Nature. 1993, 365 (6442): 170-175. 10.1038/365170a0.View ArticlePubMedGoogle Scholar
- Journot L, Spengler D, Pantaloni C, Dumuis A, Sebben M, Bockaert J: The PACAP receptor: generation by alternative splicing of functional diversity among G protein-coupled receptors in nerve cells. Semin Cell Biol. 1994, 5 (4): 263-272. 10.1006/scel.1994.1032.View ArticlePubMedGoogle Scholar
- Bokaei PB, Ma XZ, Byczynski B, Keller J, Sakac D, Fahim S, Branch DR: Identification and characterization of five-transmembrane isoforms of human vasoactive intestinal peptide and pituitary adenylate cyclase-activating polypeptide receptors. Genomics. 2006, 88 (6): 791-800. 10.1016/j.ygeno.2006.07.008.View ArticlePubMedGoogle Scholar
- Kwok YY, Chu JY, Vaudry H, Yon L, Anouar Y, Chow BK: Cloning and characterization of a PAC1 receptor hop-1 splice variant in goldfish (Carassius auratus). Gen Comp Endocrinol. 2006, 145 (2): 188-196. 10.1016/j.ygcen.2005.08.011.View ArticlePubMedGoogle Scholar
- Alexandre D, Vaudry H, Grumolato L, Turquier V, Fournier A, Jegou S, Anouar Y: Novel splice variants of type I pituitary adenylate cyclase-activating polypeptide receptor in frog exhibit altered adenylate cyclase stimulation and differential relative abundance. Endocrinology. 2002, 143 (7): 2680-2692. 10.1210/en.143.7.2680.View ArticlePubMedGoogle Scholar
- Fradinger EA, Tello JA, Rivier JE, Sherwood NM: Characterization of four receptor cDNAs: PAC1, VPAC1, a novel PAC1 and a partial GHRH in zebrafish. Mol Cell Endocrinol. 2005, 231 (1-2): 49-63. 10.1016/j.mce.2004.12.002.View ArticlePubMedGoogle Scholar
- Venkatesh B: Evolution and diversity of fish genomes. Curr Opin Genet Dev. 2003, 13 (6): 588-592. 10.1016/j.gde.2003.09.001.View ArticlePubMedGoogle Scholar
- Sarropoulou E, Franch R, Louro B, Power DM, Bargelloni L, Magoulas A, Senger F, Tsalavouta M, Patarnello T, Galibert F, Kotoulas G, Geisler R: A gene-based radiation hybrid map of the gilthead sea bream Sparus aurata refines and exploits conserved synteny with Tetraodon nigroviridis. BMC Genomics. 2007, 8: 44-10.1186/1471-2164-8-44.PubMed CentralView ArticlePubMedGoogle Scholar
- Chow BK, Yuen TT, Chan KW: Molecular evolution of vertebrate VIP receptors and functional characterization of a VIP receptor from goldfish Carassius auratus. Gen Comp Endocrinol. 1997, 105 (2): 176-185. 10.1006/gcen.1996.6818.View ArticlePubMedGoogle Scholar
- Christoffels A, Koh EG, Chia JM, Brenner S, Aparicio S, Venkatesh B: Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes. Mol Biol Evol. 2004, 21 (6): 1146-1151. 10.1093/molbev/msh114.View ArticlePubMedGoogle Scholar
- Taylor JS, Braasch I, Frickey T, Meyer A, Van De Peer Y: Genome duplication, a trait shared by 22,000 species of ray-finned fish. Genome Res. 2003, 13: 382-390. 10.1101/gr.640303.PubMed CentralView ArticlePubMedGoogle Scholar
- HEK293 Cell Database. [http://www.mbi.ufl.edu/~shaw/293.html]
- Taylor JS, Van de Peer Y, Braasch I, Meyer A: Comparative genomics provides evidence for an ancient genome duplication event in fish. Philos Trans R Soc Lond B Biol Sci. 2001, 356: 1661-1679. 10.1098/rstb.2001.0975.PubMed CentralView ArticlePubMedGoogle Scholar
- Wittbrodt J, Meyer A, Schartl M: More genes in fish? . BioEssays. 1998, 20: 511-515. 10.1002/(SICI)1521-1878(199806)20:6<511::AID-BIES10>3.0.CO;2-3.View ArticleGoogle Scholar
- Meyer A, Van de Peer Y: From 2R to 3R: evidence for a fish-specific genome duplication (FSGD). Bioessays. 2005, 27 (9): 937-945. 10.1002/bies.20293.View ArticlePubMedGoogle Scholar
- Lee LT, Siu FK, Tam JK, Lau IT, Wong AO, Lin MC, Vaudry H, Chow BK: Discovery of growth hormone-releasing hormones and receptors in nonmammalian vertebrates. Proc Natl Acad Sci U S A. 2007, 104 (7): 2133-2138. 10.1073/pnas.0611008104.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Li J, Ying Wang CY, Ho Yan Kwok AJ, Leung FC: Identification of the endogenous ligands for chicken growth hormone-releasing hormone (GHRH) receptor: evidence for a separate gene encoding GHRH in submammalian vertebrates. Endocrinology. 2007, 148 (5): 2405-2416. 10.1210/en.2006-1013.View ArticlePubMedGoogle Scholar
- Fugu Genome database. [http://fugu.biology.qmul.ac.uk]
- Ensembl database. [http://www.ensembl.org]
- Wang Y, Wong AO, Ge W: Cloning, regulation of messenger ribonucleic acid expression, and function of a new isoform of pituitary adenylate cyclase-activating polypeptide in the zebrafish ovary. Endocrinology. 2003, 144 (11): 4799-4810. 10.1210/en.2003-0501.View ArticlePubMedGoogle Scholar
- Cardoso J, Vieira FA, Gomes A, Power D: PACAP, VIP and their receptors in the metazoa: insights about the origin and evolution of the ligand-receptor pair. Peptides. 2007Google Scholar
- Aino H, Hashimoto H, Ogawa N, Nishino A, Yamamoto K, Nogi H, Nagata S, Baba A: Structure of the gene encoding the mouse pituitary adenylate cyclase-activating polypeptide receptor. Gene. 1995, 164 (2): 301-304. 10.1016/0378-1119(95)00391-I.View ArticlePubMedGoogle Scholar
- Couvineau A, Maoret JJ, Rouyer-Fessard C, Carrero I, Laburthe M: Cloning and functional characterization of the human VIP1/PACAP receptor promoter. Ann N Y Acad Sci. 1998, 865: 59-63. 10.1111/j.1749-6632.1998.tb11163.x.View ArticlePubMedGoogle Scholar
- Wong AO, Leung MY, Shea WL, Tse LY, Chang JP, Chow BK: 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. 1998, 139 (8): 3465-3479. 10.1210/en.139.8.3465.PubMedGoogle Scholar
- Tao WJ, Boulding EG: Associations between single nucleotide polymorphisms in candidate genes and growth rate in Arctic charr (Salvelinus alpinus L.). Heredity. 2003, 91 (1): 60-69. 10.1038/sj.hdy.6800281.View ArticlePubMedGoogle Scholar
- Nelson LE, Sheridan MA: Gastroenteropancreatic hormones and metabolism in fish. Gen Comp Endocrinol. 2006, 148 (2): 116-124. 10.1016/j.ygcen.2006.01.011.View ArticlePubMedGoogle Scholar
- Cinar K, Senol N, Ozen MR: Immunohistochemical study on distribution of endocrine cells in gastrointestinal tract of flower fish (Pseudophoxinus antalyae). World J Gastroenterol. 2006, 12 (42): 6874-6878.PubMed CentralPubMedGoogle Scholar
- Su Z, Wang J, Yu J, Huang X, Gu X: Evolution of alternative splicing after gene duplication. Genome Res. 2006, 16 (2): 182-189. 10.1101/gr.4197006.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong AO, Li WS, Lee EK, Leung MY, Tse LY, Chow BK, Lin HR, Chang JP: Pituitary adenylate cyclase activating polypeptide as a novel hypophysiotropic factor in fish. Biochem Cell Biol. 2000, 78 (3): 329-343. 10.1139/bcb-78-3-329.View ArticlePubMedGoogle Scholar
- Montero M, Yon L, Kikuyama S, Dufour S, Vaudry H: Molecular evolution of the growth hormone-releasing hormone/pituitary adenylate cyclase-activating polypeptide gene family. Functional implication in the regulation of growth hormone secretion. J Mol Endocrinol. 2000, 25 (2): 157-168. 10.1677/jme.0.0250157.View ArticlePubMedGoogle Scholar
- Holmgren S, Vaillant C, Dimaline R: VIP-, substance P-, gastrin/CCK-, bombesin-, somatostatin- and glucagon-like immunoreactivities in the gut of the rainbow trout, Salmo gairdneri. Cell Tissue Res. 1982, 223 (1): 141-153. 10.1007/BF00221505.View ArticlePubMedGoogle Scholar
- Fouchereau-Peron M, Laburthe M, Besson J, Rosselin G, Le Gal Y: Characterization of the vasoactive intestinal polypeptide (VIP) in the gut of fishes. Comp Biochem Physiol. 1980, 65A: 489-492. 10.1016/0300-9629(80)90064-X.View ArticleGoogle Scholar
- Thwaites DT, Young J, Thorndyke MC, Dimaline R: The isolation and chemical characterization of a novel vasoactive intestinal peptide-related peptide from a teleost fish, the cod, Gadus morhua. Biochim Biophys Acta. 1989, 999 (2): 217-220.View ArticlePubMedGoogle Scholar
- Holstein B, Humphrey CS: Stimulation of gastric acid secretion and suppression of VIP-like immunoreactivity by bombesin in the Atlantic codfish, Gadus morhua. Acta Physiol Scand. 1980, 109 (2): 217-223.View ArticlePubMedGoogle Scholar
- Mainoya JR, Bern HA: Influence of vasoactive intestinal peptide and urotensin II on the absorption of water and NaCl by the anterior intestine of the tilapia, Sarotherodon mossambicus. Zool Sci. 1984, 1 (100-105):Google Scholar
- Santos CR, Brinca L, Ingleton PM, Power DM: Cloning, expression, and tissue localisation of prolactin in adult sea bream (Sparus aurata). Gen Comp Endocrinol. 1999, 114 (1): 57-66. 10.1006/gcen.1998.7228.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Salmon Genome Project. [http://www.salmongenome.no]
- NCBI database. [http://www.ncbi.nih.gov]
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Genedoc software. [http://www.psc.edu/biomed/genedoc]
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5 (2): 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
- de Vet EC, Aguado B, Campbell RD: G6b, a novel immunoglobulin superfamily member encoded in the human major histocompatibility complex, interacts with SHP-1 and SHP-2. J Biol Chem. 2001, 276 (45): 42070-42076. 10.1074/jbc.M103214200.View ArticlePubMedGoogle Scholar
- Aguado B, Campbell RD: Characterization of a human lysophosphatidic acid acyltransferase that is encoded by a gene located in the class III region of the human major histocompatibility complex. J Biol Chem. 1998, 273 (7): 4096-4105. 10.1074/jbc.273.7.4096.View ArticlePubMedGoogle Scholar
- Laemmli UK: Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970, 227 (5259): 680-685. 10.1038/227680a0.View ArticlePubMedGoogle Scholar
- ExPASy Proteomics Server. [http://www.expasy.ch]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.