Function and distribution of 5-HT2 receptors in the honeybee (Apis mellifera).

BACKGROUND: Serotonin plays a pivotal role in regulating and modulating physiological and behavioral processes in both vertebrates and invertebrates. In the honeybee (Apis mellifera), serotonin has been implicated in division of labor, visual processing, and learning processes. Here, we present the cloning, heterologous expression, and detailed functional and pharmacological characterization of two honeybee 5-HT2 receptors.
METHODS: Honeybee 5-HT2 receptor cDNAs were amplified from brain cDNA. Recombinant cell lines were established constitutively expressing receptor variants. Pharmacological properties of the receptors were investigated by Ca(2+) imaging experiments. Quantitative PCR was applied to explore the expression patterns of receptor mRNAs.
RESULTS: The honeybee 5-HT2 receptor class consists of two subtypes, Am5-HT2α and Am5-HT2β. Each receptor gene also gives rise to alternatively spliced mRNAs that possibly code for truncated receptors. Only activation of the full-length receptors with serotonin caused an increase in the intracellular Ca(2+) concentration. The effect was mimicked by the agonists 5-methoxytryptamine and 8-OH-DPAT at low micromolar concentrations. Receptor activities were blocked by established 5-HT receptor antagonists such as clozapine, methiothepin, or mianserin. High transcript numbers were detected in exocrine glands suggesting that 5-HT2 receptors participate in secretory processes in the honeybee.
CONCLUSIONS: This study marks the first molecular and pharmacological characterization of two 5-HT2 receptor subtypes in the same insect species. The results presented should facilitate further attempts to unravel central and peripheral effects of serotonin mediated by these receptors.

Introduction

The biogenic amine serotonin (5-hydroxytryptamine, 5-HT) has important physiological functions including the regulation of energy balance and food intake, gastrointestinal and endocrine function, and cardiovascular and pulmonary physiology (for a recent review, see: [1]). In the central nervous system (CNS), serotonin contributes to the regulation of appetite, mood, sleep, and cognitive functions, including learning and memory. In humans, impaired serotonin signaling has been implicated in diseases such as migraine, schizophrenia, and depression (for reviews, see: [2], [3]). At present, 14 distinct genes encoding 5-HT receptors have been identified in total in mammals. These proteins have been assigned to seven receptor classes (5-HT1 to 5-HT7; [4], [5]). Except for 5-HT3, which is a ligand-gated ion channel, all mammalian 5-HT receptors belong to the family of G protein-coupled receptors (GPCRs). The 5-HT1 and 5-HT5 receptors couple preferentially to Gi/o proteins and inhibit cAMP synthesis. The 5-HT2 receptors activate Gq/11 proteins, which mediate the hydrolysis of inositol phosphates and cause a subsequent increase in cytosolic Ca2+ ([Ca2+]i). The 5-HT4, 5-HT6, and 5-HT7 receptors couple to Gs proteins and promote cAMP formation. Orthologous G protein-coupled 5-HT receptors with conserved signaling pathways have been discovered in protostomes including nematodes [6], crustaceans [7], [8], and insects [9]-[11]. To date, four 5-HT receptor subtypes have been characterized in the fruit fly (Drosophila melanogaster). Phylogenetically, they cluster within the mammalian 5-HT1A [12], 5-HT2 [13], and 5-HT7 [14] subfamilies.

In addition to the genetic model organism D. melanogaster, the honeybee (Apis mellifera) has greatly contributed to improving the general understanding of serotonergic control and modulation of various types of behavior [11], [15]–[19]. The distribution of serotonergic neurons has been mapped precisely in the CNS of this social insect [20], [21], and the serotonin content in the brain has been determined at various developmental stages [22], [23]. Serotonin levels are not constant during a bee’s lifespan, with older foraging bees having higher serotonin contents than younger bees working in the hive [23]. In single-cohort colonies, which consist of same-aged bees, significant differences in serotonin levels have been found in the antennal lobes of forager compared with nurse bees [24]. Together, these results suggest that differences in serotonin levels are related to specific task(s) that bees perform, rather than to their age.

Several aminergic receptors have been characterized from honeybees in recent years [25]–[31]. However, our present knowledge about the molecular and pharmacological properties of 5-HT receptors remains limited. So far, only a 5-HT1A [32] and a 5-HT7 [33] receptor have been uncovered. The aim of this study was the molecular identification and detailed pharmacological characterization of honeybee 5-HT2 receptors.

So far, not much is known about the functions of 5-HT2 receptors in insects. In D. melanogaster, the Dm5-HT2α receptor (CG1056) regulates cellular movements during germband extension and cuticular formation during early embryogenesis [34], [35]. Furthermore, 5-HT2 receptors have been implicated in the modulation of the daily activity pattern, anticipatory behavior, and aggression (for a review, see: [19]). Recently, the Cv5-HT2α receptor of the blowfly, Calliphora vicina, has been shown to be expressed not only in the brain but also in the salivary gland [36]. Here, its activation by serotonin leads, via an increase in [Ca2+]i, to an elevation of the Cl− permeability of both the basolateral and apical membrane and thus facilitates Cl− movement from the haemolymph into the lumen of the gland [36].

Inspection of the completely sequenced honeybee genome [37] has revealed the existence of two candidate genes encoding 5-HT2 receptors: Am5-HT2α and Am5-HT2β [10], [11]. After heterologous expression, both receptors cause increases in ([Ca2+]i) upon stimulation with nanomolar concentrations of serotonin. These responses are efficiently blocked by 5-HT receptor antagonists, but with subtype-specific patterns of efficacy and potency. Because of their preferential expression in glandular tissues, both 5-HT2 receptor subtypes are likely candidates for the control or modulation of important secretory processes in the honeybee.

Materials and Methods

Cloning of Am5-ht2 cDNAs

Single Apis mellifera drone brains were used to prepare poly(A)+ RNA with the Micro-FastTrack™ 2.0 Kit (Invitrogen, Karlsruhe, Germany). Drones possess a haploid genome and, therefore, single nucleotide substitutions in cDNA clones cannot be due to allelic polymorphisms. Synthesis of cDNA employed the AccuScript™ High Fidelity 1st Strand cDNA Synthesis Kit (Stratagene, Amsterdam, Netherlands). Specific primers (Table S1) allowed the entire coding region of the receptors to be amplified. The polymerase chain reaction (PCR) was carried out for 2.5 min at 94°C (1 cycle), followed by 35 cycles of 40 s at 94°C, 40 s at 54°C (Am5-ht2α) or 64°C (Am5-ht2β), 150 s at 72°C, and a final extension of 10 min at 72°C. PCR products were cloned into pGEM-T vector (Promega, Mannheim, Germany) and subsequently sequenced (AGOWA, Berlin, Germany). The nucleotide sequences of Am5-ht2α and Am5-ht2β have been submitted to the European Bioinformatics Institute (EBI) database (accession nos. FR727107 and FR727108, respectively).

Multiple sequence alignment and phylogenetic analysis

Amino-acid sequences used for phylogenetic analysis were identified by protein-protein Basic Local Alignment Search Tool (BLAST) searches of the National Center for Biotechnology Information (NCBI) database with the deduced amino acid sequence of Am5-ht2α (Am5-HT2α) as “bait”. Values for identity (ID) and similarity (S) were calculated by using the BLOSUM62 substitution matrix in BioEdit 7.0.5. MEGA 4 [38] was used to calculate the genetic distances between the core sequences and to construct maximum parsimony trees with 2000-fold bootstrap re-sampling. The D. melanogaster rhodopsin 1 (ninaE) and FMRFamide receptor sequences were used as out-groups.

Quantitative real-time PCR

Samples from various tissues of individual honeybee workers and drones were collected, immediately frozen in liquid nitrogen, and stored at -80°C until use. Total RNA was extracted by using the RNeasy Mini Kit (Qiagen, Hilden, Germany) and served as the template for cDNA synthesis. From each sample, two independent cDNA syntheses were performed by using SuperScriptIII (Invitrogen) according to the manufacturer’s instructions. Quantitative real-time PCR (qPCR) was carried out on a Rotor Gene Q (Qiagen) by using TaqMan technology with various fluorescent dyes to allow duplex measurements of receptor and reference gene expression. The fluorescent dyes used as 5’modifications were 6-FAM-phosphoramidite, Cy5, and Cy5.5. BlackBerry quencher was attached to the 3’-end of TaqMan probes. The sequences of the primers and TaqMan probes are presented in Fig. S1. The PCR was performed with an initial step at 60°C for 1 min and a denaturation step at 95°C for 5 min, followed by 45–55 cycles at 95°C for 20 s and at 60°C for 60 s. Tissue samples of at least five individual bees were examined in triplicate. Mean copy numbers were calculated by using Rotor Gene Q software (Qiagen). Samples from each tissue were tested in at least two independent qPCR runs. Receptor transcript levels were normalized to elongation factor 1α (Amef-1α) transcript levels ( =  100%) by using the standard curve method. The standards covered copy numbers from 103 – 107. To detect statistically significant differences in gene expression, one-way ANOVA followed by Bonferroni’s multiple comparison test was applied on log(x + 1) transformed data as implemented in SPSS19 (Rel. 19.0.0, IBM SPSS Statistics).

Construction of expression vectors

Expression-ready constructs of Am5-ht2α and Am5-ht2αΔIII were generated in pcDNA3.1 vector (Invitrogen). PCR was performed with specific primers (Table S1). PCR products were digested with NheI and EcoRI and subcloned. The resulting constructs were named pcAm5-ht2α-HA and pcAm5-ht2αΔIII-HA.

In order to generate an expression vector for Am5-ht2β, a PCR with specific primers (Fig. S1) was performed. The PCR product was digested with HindIII and AgeI and subcloned into pcDNA6/myc-His A vector (Invitrogen) yielding pcAm5-ht2β-His. All insert fragments were checked by DNA sequencing.

Functional expression of honeybee 5-HT2 receptors

Approximately 8 µg of pcAm5-ht2α-HA, pcAm5-ht2αΔIII-HA, and pcAm5-ht2β-His were transfected into exponentially growing HEK 293 cells (∼4×105 cells per 5-cm Petri dish) by a modified calcium phosphate method [39]. Stably transfected cells were selected in the presence of the antibiotic G418 (0.8 mg/ml; pcDNA3.1 vector) or blasticidin (0.1 mg/ml; pcDNA6 vector). Isolated foci were propagated and analyzed for the expression of the receptor proteins either by immunocytochemistry and Western blotting (see Fig. S2) or by functional Ca2+ imaging upon receptor activation.

Functional analysis of honeybee 5-HT2 receptors

The ability of the honeybee receptors to activate Gq proteins was assessed by monitoring changes in [Ca2+]i with the Ca2+-sensitive fluorescent dye Fluo-4 (Invitrogen). Non-transfected HEK 293 cells and cells expressing Am5-HT2α, Am5-HT2αΔIII, or Am5-HT2β were grown in 96-well plates to a density of ∼3×104 per well. Cells were loaded at room temperature with Fluo-4 as described earlier [40]. Each 96-well plate was transferred into a fluorescence reader (FLUOstar Galaxy/Optima; BMG Labtech, Offenburg, Germany) to monitor Fluo-4 fluorescence. The excitation wavelength was 485 nm, and fluorescence emission was detected at 520 nm. Various concentrations of biogenic amines or receptor ligands were added, once Fluo-4 fluorescence had reached a stable value in each well. The changes in Fluo-4 fluorescence were recorded automatically. Dose-response curves for putative agonists/antagonists were established in at least two independent experiments.

Results

Cloning and sequence analysis of 5-HT2 receptors from A. mellifera

Two candidate genes encoding putative serotonin 5-HT2 receptors were identified in the completely sequenced honeybee genome (Am16 and Am17; [10]). Here, we used this sequence information and applied a PCR-based strategy to amplify the full-length cDNAs. The Am5-ht2α cDNA (Am16 in [10]) contains an open reading frame (ORF) of 1,962 bp and encodes a protein of 653 amino-acid residues (Am5-HT2α) with a calculated molecular weight of 71.5 kDa. The hydropathy profile and topology predictor Phobius [41] suggest seven trans-membrane (TM) domains, which are characteristic of GPCRs. The TM segments are flanked by an extracellular N-terminus of 140 and an intracellular C-terminus of 96 residues. An alignment of the cDNA with the genomic database (release Amel4.5) revealed that the Am5-ht2α gene contains seven exons and six introns (Fig. 1), and that it is located on chromosome LG9. Am5-ht2α is the ortholog of the D. melanogaster 5-ht2α gene (CG1056) with which it has three introns in common (Fig. S3A+C). A remarkable difference exists between the cloned cDNA and the annotated Am5-ht2α gene sequence (accession no.: XM394798). Exon V of the cDNA consists of 270 bp, whereas the annotated exon V contains only 220 bp (Fig. 1, Vannot). This discrepancy originates from the usage of a different 5’ splice site during the in silico assembly of the gene, such that the annotated exon V starts at position 1,202 of the Am5-ht2α ORF. The 3’ splice sites of both the cloned and the annotated sequences are completely conserved (Fig. 1). Notably, in all PCR experiments performed, we never obtained a fragment as predicted by the annotation.

Figure 1

Genomic organization of Am5-ht2α and Am5-ht2β genes.

The Apis mellifera genomic database (Amel_4.5) was screened with the Am5-ht2α and Am5-ht2β cDNA sequences. The figure shows the splicing pattern of the Am5-ht2α (upper part) and Am5-ht2β (lower part) genes. Nucleotide sequences of exons are given in upper-case letters and nucleotide sequences of introns are given in lower-case letters. Full-length and truncated receptor variants are displayed schematically. Boxes represent exons, and bent lines represent introns. The position of the seven trans-membrane domains is highlighted by black bars. Exons containing a frameshift are indicated in gray. The seven exons of the Am5-ht2α gene are completely covered in the genomic contig NW_001253562. The six exons of the Am5-ht2β gene are dispersed in three genomic contigs (A: NW_001262451, B: NW_001252966, und C: NW_001262048).

However, we isolated a splice variant of Am5-ht2α mRNA (Am5-ht2αΔIII). In Am5-ht2αΔIII, exon III (225 bp) is absent, which results in a shortened ORF of 1,737 bp (Fig. 1). This mRNA variant can potentially be translated into a protein of 578 amino acid residues. Compared with the full-length receptor, this protein will lack TM3, TM4, and the entire CPL2 (see Fig. 1).

In addition to Am5-ht2α, three partial mRNAs that might encode proteins have been annotated in silico. However, with 256 bp (XM_001119970), 522 bp (XM_001122856), and 1,053 bp (XM_624894), these sequences are too short to code for complete GPCRs. Nevertheless, the deduced amino-acid sequences showed striking sequence similarity to a 5-HT receptor of the spiny lobster, Panulirus interruptus (Pan5-HT2β; [7]; accession no.: AY550910). We used the crustacean sequence to design primers annealing in close proximity to the start and the stop codons of the hypothetical honeybee ortholog. A single cDNA, which contained an ORF of 2,202 bp, was amplified on honeybee brain cDNA (Am5-ht2β). The ORF is distributed on three genomic contigs (Fig. 1). Contig NW_001252966 has been allocated to chromosome LG10, whereas the other two, i.e., NW_001262451 and NW_001262048, have not been assigned to any chromosome, as yet. The coding sequence of the gene is interrupted by five introns (Fig. 1) with two of them being in common with the orthologous gene in D. melanogaster (CG42796, Fig. S3B+C). The ORF of the Am5-ht2β gene encodes a protein (Am5-HT2β) of 733 amino-acid residues with a calculated molecular weight of 80.7 kDa. The Am5-HT2β receptor contains an extremely long CPL3 of 399 residues and a relatively short C-terminus of 25 residues.

Similar to Am5-ht2α, we identified a splice variant of the Am5-ht2β mRNA. In this variant (Am5-ht2βΔII), exon II is completely missing. This results in a frameshift of the ORF (ORF1) with an early stop codon at position 1,075 - 1,077 (Fig. 1). ORF1 is predicted to code for a protein of 358 amino-acid residues and a calculated molecular weight of 39.5 kDa. The sequence from Met1 to Gly132 including TM1, TM2, and the extracellular loop (ECL) 1 is identical to the full-length Am5-HT2β protein. Because of the frameshift, the consecutive sequence (Gly133 to Pro357) differs completely from the full-length receptor protein. A second ORF (ORF2, 1,605 bp) starts within exon III and lasts to the stop codon of the full-length mRNA. This reading frame is predicted to encode a protein of 491 residues with a calculated molecular weight of 54.1 kDa. The amino-acid sequence is completely identical to that of full-length Am-5-HT2β from Met243 (ECL2 prior to TM5) to Arg733 but obviously misses a large part of the N-terminal region encompassing TM domains 1-4 of an intact GPCR (Fig. 1).

Sequence motifs, which are essential for the three-dimensional structure, ligand binding, and signal transduction of the receptor, are well conserved in both full-length Am5-HT2α and Am5-HT2β proteins (Fig. 2). An exception is the Asp residue within the D-R-Y motif required for G protein coupling [42]; this is substituted by Gly in Am5-HT2α (G237-R238-Y239, Fig. 2). Five and four consensus motifs for potential N-glycosylation (N-x-[S/T]) are located in the extracellular N-terminus of Am5-HT2α and Am5-HT2β, respectively (Fig. 2). Cysteine residues in the C-terminus (Cys570, Cys572, and Cys594 in Am5-HT2α, Cys722 in Am5-HT2β, Fig. 2) are possible sites for post-translational palmitoylation. Whereas no consensus site for phosphorylation by protein kinase A (PKA, R-[R/K]-x-[S/T]) or protein kinase C (PKC, [S/T]-x-[R/K]) is found in Am5-HT2α, five PKA consensus sites and four PKC consensus sites are present within CPL3 of Am5-HT2β (Fig. 2).

Figure 2

Amino-acid sequence alignment of Am5-HT2α, Am5-HT2β, and orthologous receptors from Drosophila melanogaster (Dm5-HT2; accession no. CAA57429) and Panulirus interruptus (Pan5-HT2β; AAS57919).

Identical residues between all four receptors are shown as white letters against black, whereas conservatively substituted residues are shaded. Putative trans-membrane domains (TM1–7) are indicated by gray bars. Putative consensus sites for post-translational modifications in Am5-HT2α and Am5-HT2β are indicated by open and filled symbols, respectively: N-glycosylation sites are shown as triangles, PKA phosphorylation sites as diamonds, PKC phosphorylation sites as circles, palmitoylation sites as squared boxes. Underlined letters represent the region within the CPL3 of Am5-HT2α used to raise specific antibodies (see Fig. S2). The amino acid position is indicated on the right.

A comparison of Am5-HT2α and Am5-HT2β amino-acid sequences with NCBI databases identified several orthologous protostomian and deuterostomian 5-HT2 receptors. For Am5-HT2α, the highest amino acid identity/similarity (ID/S) existed to the 5-HT2 receptor of D. melanogaster (Dm5-HT2α; [13]; ID 28.4%, S 34.8%). For Am5-HT2β, homology was more pronounced to 5-HT2 receptors from the crustaceans Panulirus interruptus (Pan5-HT2β; [7]; ID 31.2%, S 42.8%) and Procambarus clarkii (Pro5-HT2; [8]; ID 33.2%, S 40.4%). In phylogenetic tree analyses, both Am5-HT2α and Am5-HT2β always grouped with protostomian and deuterostomian 5-HT2 receptors (Fig. 3).

Figure 3

Phylogenetic relationship between various 5-HT receptors from vertebrates, mollusks, and arthropods.

Alignments were performed with BioEdit (version 7.0.5) by using the core amino-acid sequences lacking the variable regions of the N- and C-terminus and the third cytoplasmic loop. Maximum parsimony analysis was calculated with MEGA4. The receptor sequences followed by their accession numbers are listed in the order illustrated: Homo sapiens (Hs5-HT1B, acc. no. NP_000854), Danio rerio (Dr5-HT1B, AAI63698), Hs5-HT1D (NP_000855), Hs5-HT1E (NP_000856), Hs5-HT1F (NP_000857), Dr5-HT1F (NP_001076425), Hs5-HT1A (NP_000515), Dr5-HT1A (NP_001139238), Aplysia californica (Ac5-HT1, AAM46088, AAC28786), Panulirus interruptus (Pi5-HT1, AY528822), Apis mellifera (Am5-HT1, CBI75449), Drosophila melanogaster (Dm5-HT1A, CAA77570), Dm5-HT1B (CAA77571), Hs5-HT5A (NP_076917), Dr5-HT5A, (NP_001119882), Am5-HT7 (AM076717), Dm5-HT7 (NP_524599), Aplysia kurodai (Ak5-HT7, ACQ90247), Hs5-HT7A (NP_000863), Hs5-HT7D (NP_062873), Hs5-HT7B (NP_062874), Hs5-HT4 (NP_001035259), Ac5-HT2 (Q16951), Ac5-HT1 (Q16950), Hs5-HT6 (NP_000862), Hs5-HT2B (NP_000858), Dr5-HT2B (ABI18978), Hs5-HT2A (NP_000612), Dr5-HT2A (CAQ15355), Hs5-HT2C (NP_000859), Dr5-HT2C (CAX14715), Am5-HT2β (FR727108), Dm5-HT2β (NM_141548/NM_169229), Pi5-HT2β (AY550910), Am5-HT2α (FR727107), and Dm5-HT2α (CAA57429). The divergent D. melanogaster ninaE-encoded rhodopsin 1 (DmninaE, NM_079683) and D. melanogaster FMRFamide receptor (DmFR, AAF47700) were used as out-groups. The numbers at the nodes of the branches represent the percentage bootstrap support (2,000 replications) for each branch.

Tissue-specific expression patterns of the Am5-HT2 genes

The tissue-specific distribution of Am5-ht2α and Am5-ht2β mRNA was determined by qPCR in tissue samples from pollen foragers and drones (Fig. 4, Fig. S4). Application of splice-variant-specific primers (Fig. S1) allowed us to unravel the expression pattern of the different splice variants for each receptor gene. PCR fragments originating from full-length and splice variants of Am5-ht2α and Am5-ht2β were found in the CNS, the hypopharyngeal glands, the salivary glands, and the Malpighian tubules of pollen foragers (Fig. 4; see Fig. S4 for drone tissues). Expression of Am5-ht2α was higher than Am5-ht2β expression, e.g., up to ten times in the nervous system and ∼200 times in the hypopharyngeal glands. Interestingly, in the brain and in the hypopharyngeal gland, the expression of the Am5-ht2αΔIII transcript was lower than the full-length version, whereas the opposite held true for the Am5-ht2β and Am5-ht2βΔII transcripts (Fig. 4). All transcripts showed similar expression patterns in the CNS of drones compared with female honeybees (Fig. 4A+B; Fig. S4).

Figure 4

Tissue-specific expression patterns of Am5-ht2 receptor genes in pollen foragers determined by quantitative real-time PCR.

Transcript levels were normalized to Amef-1α. (A) Relative levels of receptor-gene mRNAs in the brain (n = 11). (B) Relative levels of receptor-gene mRNAs in the ventral nerve cord (Am5-ht2α and Am5-ht2αΔIII: n = 11; Am5-ht2β and Am5-ht2βΔII: n = 10) (C) Relative levels of receptor-gene mRNAs in the hypopharyngeal gland (Am5-ht2α: n = 6; Am5-ht2αΔIII: n = 5; Am5-ht2β: n = 7; Am5-ht2βΔII: n = 6). (D) Relative levels of receptor-gene mRNAs in the head salivary gland (Am5-ht2α and Am5-ht2αΔIII: n = 5; Am5-ht2β and Am5-ht2βΔII: n = 7). (E) Relative levels of receptor-gene mRNAs in the thoracic salivary gland (each: n = 5). (F) Relative levels of receptor-gene mRNAs in Malpighian tubules (Am5-ht2α: n = 5; Am5-ht2αΔIII: n = 6; Am5 ht2β: n = 8; Am5-ht2βΔII: n = 7). Groups that differed significantly from one another in relative mRNA levels are indicated with different letters above the box plots (p < 0.05, Bonferroni’s multiple comparison test).

Functional analyses of Am5-HT2 receptors in HEK 293 cells

5-HT2 receptors are known to induce an increase in [Ca2+]i. We generated HEK 293 cell lines constitutively expressing either Am5-HT2α or Am5-HT2β to examine whether these honeybee receptors also induced serotonin-dependent Ca2+ signals. In a first series of experiments, we analyzed the ligand specificity of both receptors. Only serotonin (10 µM) increased [Ca2+]i in Am5-HT2α- and Am5-HT2β-expressing cells, whereas no signals were observed after the application of 10 µM tyramine, octopamine, or dopamine (Fig. 5). Non-transfected cells did not respond to any of these biogenic amines (Fig. 5). A third cell line expressing the alternatively spliced variant Am5-HT2αΔIII did not respond to the application of any biogenic amine, including serotonin. Therefore, this construct most likely does not code for a bona fide GPCR.

Figure 5

Biogenic amine (10 µM each) modulation of intracellular Ca2+ concentrations in HEK 293 cells constitutively expressing Am5-HT2α or Am5-HT2ß receptors and in non-transfected HEK 293 cells.

Cells were loaded with the Ca2+-sensitive dye Fluo-4 and changes in [Ca2+]i were measured fluorimetrically. All values of HEK 293 cells stably expressing Am5-HT2β (black bars) and the corresponding non-transfected HEK 293 cells (gray bars) were normalized to the serotonin response of the Am5-HT2β-expressing cell line. The values obtained for the Am5-HT2α cell line (white bar) and non-transfected HEK 293 cells (striped bar) were normalized analogously. Data represent the mean ± SD of octuplicate determinations from a representative experiment. The values for serotonin-stimulated Am5-HT2α- and Am5-HT2β-expressing cells are significantly different from all other values (one-way ANOVA, followed by Dunnett’s multiple comparison test, p < 0.05).

For both full-length receptors, the concentration-response relationship for serotonin was examined with concentrations ranging from 1 nM to 300 µM. In cell lines expressing either Am5-HT2α or Am5-HT2β, the serotonin effect was concentration-dependent and saturable, resulting in sigmoidal concentration-response curves (Fig. 6A). The half-maximal effective concentrations (EC50) were similar at 25.7 nM and 32.5 nM for Am5-HT2α and Am5-HT2β, respectively. Maximal increase in [Ca2+]i was observed at serotonin concentrations of ≥1 µM.

Figure 6

Concentration-dependent effects of agonists and antagonists on [Ca2+]i in Am5-HT2α- and Am5-HT2β-expressing HEK 293 cells.

[Ca2+]i is depicted as relative fluorescence given in percent. Data represent the mean ± SEM of eight replicates from one experiment representative of at least two similar experiments. (A) Am5-HT2α- or Am5-HT2β-expressing cells were incubated with increasing concentrations of serotonin. Values were normalized to the response measured at the highest serotonin concentration ( = 100%) for each cell line. (B) Cells expressing Am5-HT2α or Am5-HT2β were incubated with increasing concentrations of the agonists 5-methoxytryptamine (5-MeOT) and 8-hydroxy-DPAT (8-OH-DPAT). Values were normalized to the response measured at the highest 5-MeOT concentration ( = 100%) for each cell line. (C, D) Cells expressing Am5-HT2α or Am5-HT2β were incubated with 50 nM serotonin and increasing concentrations of the antagonist mianserin (C) or clozapine (D). Values were normalized to the serotonin response in the absence of antagonist ( = 100%).

To establish pharmacological profiles for both 5-HT2 receptors, the effects of various 5-HT receptor agonists and antagonists were tested. Two agonists, namely 5-methoxytryptamine and 8-hydroxy-2-(di-n-propylamino)tetralin (8-OH-DPAT), caused an increase in [Ca2+]i in both receptor-expressing cell lines (Fig. 6B). At the Am5-HT2α receptor, the EC50 values were 70 nM and 55.9 µM for 5-methoxytryptamine and 8-OH-DPAT, respectively. At the Am5-HT2β receptor, the EC50 for 5-methoxytryptamine was 60.4 nM, and for 8-OH-DPAT, the EC50 was 561.5 nM. Neither 5-methoxytryptamine nor 8-OH-DPAT had an effect on non-transfected HEK 293 cells. In contrast to the previously mentioned substances, 5carboxamidotryptamine and 2,5-dimethoxy-4-iodoamphetamine even affected [Ca2+]i in non-transfected cells and thus were not examined with the receptor-expressing cell lines.

The potency of 5-HT-receptor antagonists to inhibit serotonin-induced Ca2+ responses in Am5-HT2α- and Am5-HT2β-expressing cells was examined in another series of experiments. The calculated IC50 and maximal inhibitory effect (Imax) of the substances are summarized in Table 1, and examples of concentration-response relationships for mianserin and clozapine are shown in Figures 6C+D. Most antagonists were able to inhibit the activation of the Am5-HT2α receptor by 30-60%. At the Am5-HT2α receptor, antagonists showed decreasing inhibitory potency with: SB200646  =  mianserin > cyproheptadine > methiothepin > clozapine > methysergide.

Table 1

Antagonist profile of Apis mellifera 5-HT2 receptors.

antagonist	Am5-HT2αIC50 [µM]	Am5-HT2αmaximal inhibition	Am5-HT2βIC50 [µM]	Am5-HT2βmaximal inhibition
clozapine	1.7	56%	6.1	95%
cyproheptadine	0.7	64%	0.2	77%
ketanserin	-	no effect	1.5	89%
methiothepin	0.6	61%	-	no effect
methysergide	2.2	45%	-	no effect
mianserin	0.2	51%	2.4	76%
SB-200646	0.1	32%	-	no effect

IC50 values (potency) and relative efficacy were calculated from concentration-response curves for each drug. Efficacy is given as the maximal inhibition (%) of fluorescence compared with 50 nM serotonin in the absence of antagonist and calculated by curve fitting. Values are means of the results of two independent experiments, each performed in octuplicates.

Interstingly, SB-200646, methiothepin, and methysergide showed no effects at the Am5-HT2β receptor, whereas clozapine, cyproheptadine, ketanserin, and mianserin were strongly inhibitory with a maximal inhibitory effect of ≥77%. In comparison with Am5-HT2α, the rank order of inhibitory potency at Am5-HT2β was different with: cyproheptadine > ketanserin > mianserin > clozapine.

Discussion

The intense arborization of serotonergic neurons in the honeybee CNS suggests that serotonin controls the activity of many neuronal circuits [11]. In order to obtain further insight into serotonergic signaling in the honeybee, the molecular and functional identity of the 5-HT receptor repertoire of the honeybee has to be unraveled. Concomitantly with the release of the honeybee genomic sequence [37], several 5-HT receptor genes were annotated [10]. These in silico analyses, however, lack experimental data related to pharmacological properties, expression patterns, and physiological functions of the identified receptors. So far, only a 5-HT1A receptor [32] and a 5-HT7 receptor [33] have been thoroughly characterized. With the analysis of the two paralogous 5-HT2 receptors, Am5-HT2α and Am5-HT2β, the current work provides an important step to complement our knowledge of the serotonergic receptor system in the honeybee.

Genomic organization and gene prediction

The coding sequences of Am5-ht2α and Am5-ht2β are heavily spliced and consist of seven and six exons, respectively. Two splice sites are located in similar positions, one in ECL1 and a second in ECL2. These two splice sites are not only conserved in orthologous genes of the fruit fly (Fig. S3) but also represent the positions from where the splice variants originate (Fig. 1).

A comparison of cloned and annotated Am5-ht2-receptor gene sequences reveals several discrepancies. Whereas the overall exon-intron structure is similar, the annotated Am5-ht2α sequence contains a shorter exon V (Fig. 1). Obviously, native honeybee spliceosomes not only use exactly conserved splice donor and acceptor sites, which, however, are central elements in gene prediction algorithms [43]. The situation for the Am5-ht2β gene is even more complicated. Alignment of the cloned cDNA with the honeybee genome sequence showed that: (1) the Am5-ht2β gene contains extremely large introns (>100,000 bp), which are hardly recognized by gene prediction software [44]; (2) the Am5-ht2β gene is dispersed over three genomic contigs; (3) the ORF of Am5-ht2β is oriented in either a 5’ to 3’ or 3’ to 5’ orientation compared with deposited contig sequences. This obviously prevents a computer-based assembly of the correct ORF. As exemplified by Am5-ht2β, gene annotations provide important information, but molecular cloning of cDNA remains a prerequisite to determine the functional properties of a protein.

Phylogenetic analysis of 5-HT2 receptors

Phylogenetic analyses of various 5-HT receptors revealed a separation of two main groups (Fig. 3). One group comprises receptors coupled to the cAMP pathway, i.e., 5-HT1 (Gi-coupled) and 5-HT7 receptors (Gs-coupled). The other group comprises 5-HT2 receptors preferentially coupled to Gq proteins. Interestingly, deuterostomian and protostomian 5-HT receptors do not form distinct branches suggesting common ancestors for each receptor class. In arthropods, usually one or two receptor subtypes exist per 5-HT receptor class. In D. melanogaster, the Dm5-HT1A and Dm5-HT1B receptors result from a recent duplication event [12]. So far, one 5-HT1 receptor has been characterized in the honeybee (Am5-HT1A, [32]). In addition, one 5-HT7 receptor has been found in D. melanogaster and A. mellifera [14], [33]. Here, we have added two functional 5-HT2 receptors to the list of honeybee 5-HT receptors. A receptor orthologous to Am5-HT2α has been characterized in D. melanogaster (CG1056; [13]). The D. melanogaster genome contains another gene (CG42796) potentially encoding an ortholog of Am5-HT2β. This hypothesis, however, remains to be established experimentally.

Am5-ht2 genes give rise to several transcripts by alternative splicing

Despite the apparent functional diversity afforded by the existence of multiple 5-HT2 receptor subtypes, the existence of functional splice variants of these receptors adds to the potential repertoire of 5-HT receptors an additional mode of regulation. This phenomenon was observed for several 5-HT2 receptor subtypes (Ascaris suum 5-HT2: [45]; 5-HT2A: [46]; 5-HT2C: [47], [48]). Splice variants differ in distribution, ligand-binding properties and functional signaling (for a review, see: [49]).

For both honeybee 5-HT2 receptor genes, only one splice product contains information for a complete GPCR (Fig. 1). Shorter versions are generated by exon skipping. Interestingly, Am5-ht2αΔIII and Am5-ht2βΔII have lost the exon coding for a central part of the receptor including TM3 and TM4 (Fig. 1). Several 5-HT2 receptor splice variants affecting the same region were also reported in mammalian genes [46], [47]. A deletion covering the region encoding CPL2 and TM4 leads to a frameshift in the downstream sequence and a premature stop codon of the rat 5-HT2C receptor [47]. Similarly truncated forms were also detected in cDNA samples from mouse and human brain [47]. A variant of the human 5-HT2A receptor contains an 118 bp insertion which produces a frameshift and a premature stop codon [46]. Unfortunately, nothing is known about the function(s) of truncated 5-HT2 receptors in mammals.

When heterologously expressed, the Am5-HT2αΔIII receptor is synthesized as a glycosylated and non-glycosylated protein suggesting that it is processed in the Golgi apparatus (Fig. S2A). Attempts to monitor functional receptor activity from this construct were unsuccessful. This might be attributable to the absence, in TM3, of critical residues that interact with the amino group of biogenic amines [50]. Although the protein is not functional on its own, it might assemble with a full-length receptor and thereby modulate the binding and/or signal transduction properties in a heterodimeric complex. Such modes of regulation have been observed for a Caenorhabditis elegans dopamine receptor variant lacking TM6 and TM7 [51] and for mutants of the human dopamine D2 receptor [52].

Deletion of exon II in Am5-ht2βΔII leads to a frameshift in the ORF, which is terminated by an early stop codon. Thus, ORF1 of Am5-ht2βΔII codes for a protein that shares, with full-length Am5-HT2β, the N-terminal part down to ECL1 but that thereafter contains a completely divergent sequence. A second ORF of this splice variant codes for a protein identical to the full-length receptor from TM5 down to the C-terminus (Fig. 1). Whether any of these proteins can be heterologously expressed and eventually modify the properties of the full-length receptor remains to be examined in a forthcoming study.

Functional characterization of Am5-HT2 receptors

Am5-HT2α and Am5-HT2β were functionally expressed in HEK 293 cells. Like their mammalian orthologs, both receptors couple to the Gq signaling pathway. These results support the concept that the signaling pathway is conserved for a given receptor class, whereas pharmacological profiles might differ.

With EC50 values of 25.7 nM and 32.5 nM, respectively, Am5-HT2α and Am5-HT2β displayed similar potencies for serotonin as those also reported for 5-HT2α receptors from D. melanogaster (Dm5-HT2α, 16 nM; [13]) and Calliphora vicina (Cv5-HT2α, 24 nM; [36]). In addition to serotonin, 5-methoxytryptamine also activates honeybee 5-HT2 receptors but with 2-fold higher EC50 values (Fig. 6B). In the blowfly C. vicina, nanomolar concentrations of 5-methoxytryptamine activate only Cv5-HT2α (EC50 = 67 nM) but not Cv5-HT7 [36]. However, at higher concentrations (EC50 = 55.9 µM and 560 nM, respectively), both honeybee receptors are also activated by 8-OH-DPAT, as is Cv5-HT2α (EC50 = 62 µM). In conclusion, both 5-methoxytryptamine and 8-OH-DPAT can be considered as non-selective 5-HT receptor agonists in the honeybee, because 5methoxytryptamine also activates the Am5-HT1A receptor [32], and because 8-OH-DPAT is a poor agonist at the Am5-HT7 receptor [33].

In general, the pharmacological profiles of both Am5-HT2 receptors are different for antagonists, especially for SB-200646, methiothepin, and methysergide, which display high potency at Am5-HT2α but do not block Am5-HT2β. To be considered as an Am5-HT2α-specific antagonist, however, these substances need to be re-examined at the Am5-HT1 and Am5-HT7 receptors. Clozapine displays similar potencies for both 5-HT2 receptors. In contrast, mianserin has a 24fold higher potency at the Am5-HT2α receptor (Table 1) and thus is a promising candidate for an Am5-HT2α–specific antagonist. Interestingly, ketanserin seems to be a specific antagonist for Am5-HT2β. However, these substances also remain to be tested with Am5-HT1 and Am5-HT7.

We have demonstrated that many drugs are active at several 5-HT receptor subtypes of the honeybee ([32], [33]; this study). Thus, when applying an agonist to elicit a physiological or behavioral “serotonin effect”, the possible activation of multiple 5-HT receptors or even other aminergic receptors simultaneously should not be overlooked.

Localization of Am5-HT2 receptors

The expression patterns of only two 5-HT receptors have been reported in the honeybee. The Am5-HT1A receptor is restricted to neuronal tissues [32]. Transcripts of Am5-HT7 have been detected in neural and peripheral tissues [33]. Based on our qPCR data, both Am5-ht2 receptor genes share the widespread expression pattern with Am5-ht7.

Transcript levels for Am5-ht2α were clearly higher than those for Am5-ht2β, especially in the hypopharyngeal gland (Fig. 4C). The hypopharyngeal gland plays an important role in the control of eusocial behavior and displays ‘functional flexibility’ in secreting either royal jelly compounds or digestive enzymes depending on the age and task of the worker bee [53], [54]. Queen caste development strongly depends on royal jelly [55], [56]. A tempting speculation is that glandular activity is modulated by serotonin, and thus, signaling via Am5-HT2α is crucial for generating the female reproductive caste in honeybees.

The Am5-ht2β gene is expressed in the nervous system of worker bees and drones. In D. melanogaster, the orthologous receptor (CG42796) is preferentially expressed in the male CNS [57]. Thus, 5-HT2β receptors might be involved in different functional circuits in these insects.

With the present work, we extend the basic knowledge of the serotonergic system of the honeybee by adding two receptors, Am5-HT2α und Am5-HT2β, to the growing list of well-characterized honeybee GPCRs. The detailed knowledge of the pharmacological properties of a receptor is essential for designing and conducting targeted pharmacological experiments (e.g., [32]). In combination with expression data, knowledge of the pharmacological profiles of these receptors should now facilitate such studies for individual 5-HT receptor subtypes in the honeybee.

Sequences of primers and TaqMan probes (including 5’- and 3’-modifications; see material and methods) used for qPCR assays and the expected length of the resulting amplicons. The positions of primers on the cDNAs are schematically shown by open arrows, TaqMan probes are indicated as asterisks.

(PDF)

Click here for additional data file.

Western blot and immunocytochemical analyses of Am5-HT Anti-Am5-HT2α antibodies were raised against a fusion protein containing part of the third cytoplasmic loop (CPL3; amino acid Arg340 to Glu452; see Fig. 1). The cDNA fragment was amplified by PCR with specific primers (Table S1). The fragment was cloned into pMAL-c2X vector (New England Biolabs, Frankfurt, Germany). The fusion protein was over-expressed in E. coli BL21 CP and purified by amylose affinity-chromatography (New England Biolabs). In collaboration with the Nachwuchsgruppe Antikörper-Technologien (University of Potsdam, Germany), the fusion protein was used to immunize mice and to raise monoclonal antibodies. A second fusion protein containing a His-tag attached to the same receptor fragment was expressed from pET-30a vector (Novagen, Darmstadt, Germany) and used for testing the specificity of the monoclonal antibodies. Membrane proteins (10 µg protein per lane) of human embryonic kidney cells (HEK 293) expressing Am5-HT2α-HA and Am5-HT2αΔIII-HA receptors (see below) were isolated as previously described (Thamm et al., 2010). Proteins were separated by SDS polyacrylamide gel electrophoresis on 10% or 12% gels and transferred to polyvinylidene fluoride membranes (Roth, Karlsruhe, Germany). These membranes were blocked with 5% (w/v) dry milk in Tris-buffered saline containing Tween 20 (TBS-T, 10 mM Tris-HCl, pH 7.5, 150 mM NaCl, 0.01% Tween 20) for 30 min at room temperature, incubated either with specific anti-HA antibodies (Anti-HA High Affinity, Roche, Penzberg, Germany; dilution 1∶5,000) or with receptor-specific antibodies (dilution 1∶100) in TBS-T, washed with TBS-T, and finally incubated with secondary antibodies (1∶5,000, anti-rat-HRP; American Qualex, La Mirada, USA; 1∶200, anti-mouse Alexa568; Invitrogen) for 1 h. Signals were visualized by enhanced chemiluminescence. (A) Western blot analyses of membrane proteins (10 µg per lane) of non-transfected HEK 293 cells (nt) and HEK 293 cells expressing either full-length Am5-HT2α-HA (full) or Am5-HT2αΔIII-HA (ΔIII) proteins. Both anti-HA (dilution 1:5,000, left) and anti-Am5-HT2α (culture supernatant 1∶10, right) antibodies recognize bands of identical size in protein preparations from transfected cells. No bands were detected in protein preparations from non-transfected cells. (B) Similar staining patterns were observed in receptor-expressing cell lines with both antibodies (anti-HA and anti-Am5-HT2α). Scale bar 40 µm). Unfortunately, the anti-Am5-HT2α antibody did not work with native tissue, neither on Western Blots nor on fixed tissue sections. This is probably due to the very low endogenous expression level of Am5-HT2α, especially in nervous tissue.

(PDF)

Click here for additional data file.

Comparison of splice sites in genes (coding regions) of and . Splice donor and acceptor sites were identified using the Splign alignment tool (http://www.ncbi.nlm.nih.gov/sutils/splign/splign.cgi). In the mRNA sequence alignments (A+B), which were created using the program BioEdit version 7.1.3.0, splice sites are indicated by +++ and highlighted in yellow and green for 5-ht2 genes of the honeybee and D. melanogaster, respectively. (A) Am5-ht2α is the orthologue of the Dm5-ht2α gene (CG1056) with which it has three introns in common. (B) Am5-ht2β is the orthologue of the D. melanogaster CG42796 gene with which it has two introns in common. (C) A pictogram of a GPCR with its seven transmembrane segments (green bars) is displayed. Arrows point to the relative positions of splice sites in the primary structures of the receptors. For each exon, the last amino acid residue is indicated and numbered according to its position in the deduced amino acid sequence. One splice site is conserved in all four 5-ht2 genes (red arrow). Two additional splices sites are conserved in the two 5-ht2α genes only (green arrows) whereas the two 5-ht2β genes have one additional splice site in common (blue arrow).

(PDF)

Click here for additional data file.

Tissue-specific expression patterns of receptor genes in drones determined by quantitative real-time PCR. Transcript levels were normalized to Amef-1α. (A) Relative levels of receptor-gene mRNAs in the brain of drones (Am5-ht2α and Am5-ht2αΔIII: n = 10; Am5-ht2β: n = 6; Am5-ht2βΔII: n = 5). (B) Relative levels of receptor-gene mRNAs in the ventral nerve cord of drones (Am5-ht2α, Am5-ht2αΔIII, and Am5-ht2β: n = 7; Am5-ht2βΔII: n = 6). Groups that differed significantly in relative mRNA levels within a given tissue are indicated with different letters above the box plots (p < 0.05, Bonferroni’s multiple comparison test).

(PDF)

Click here for additional data file.

Sequences of primers used for full-length cloning of cDNAs, construction of expression vectors, and production of monoclonal antibodies against Am5-HT

(DOCX)

Click here for additional data file.