Candidate chemoreceptor subfamilies differentially expressed in the chemosensory organs of the molluscAplysia
© Cummins et al. 2009
Received: 06 November 2008
Accepted: 04 June 2009
Published: 04 June 2009
Skip to main content
© Cummins et al. 2009
Received: 06 November 2008
Accepted: 04 June 2009
Published: 04 June 2009
Marine molluscs, as is the case with most aquatic animals, rely heavily on olfactory cues for survival. In the molluscAplysia californica, mate-attraction is mediated by a blend of water-borne protein pheromones that are detected by sensory structures called rhinophores. The expression of G protein and phospholipase C signaling molecules in this organ is consistent with chemosensory detection being via a G-protein-coupled signaling mechanism.
Here we show that novel multi-transmembrane proteins with similarity to rhodopsin G-protein coupled receptors are expressed in sensory epithelia microdissected from theAplysia rhinophore. Analysis of theA. californica genome reveals that these are part of larger multigene families that possess features found in metazoan chemosensory receptor families (that is, these families chiefly consist of single exon genes that are clustered in the genome). Phylogenetic analyses show that the novelAplysia G-protein coupled receptor-like proteins represent three distinct monophyletic subfamilies. Representatives of each subfamily are restricted to or differentially expressed in the rhinophore and oral tentacles, suggesting that they encode functional chemoreceptors and that these olfactory organs sense different chemicals. Those expressed in rhinophores may sense water-borne pheromones. Secondary signaling component proteins Gαq, Gαi, and Gαo are also expressed in the rhinophore sensory epithelium.
The novel rhodopsin G-protein coupled receptor-like gene subfamilies identified here do not have closely related identifiable orthologs in other metazoans, suggesting that they arose by a lineage-specific expansion as has been observed in chemosensory receptor families in other bilaterians. These candidate chemosensory receptors are expressed and often restricted to rhinophores and oral tentacles, lending support to the notion that water-borne chemical detection inAplysia involves species- or lineage-specific families of chemosensory receptors.
All animals must recognize and respond to chemosensory information in their environment. Although the marine molluscAplysia has been a valuable model to investigate the molecular basis of behavior [1, 2] and reproduction [3, 4], our knowledge of how they recognize and respond to environmental signals is limited. In particular, it is unknown how they distinguish and bind water-soluble molecules and transfer exogenous information intracellularly. In contrast, the molecular components and mechanisms of chemical detection in a range of vertebrates and other invertebrates have been well studied.
Vertebrate chemoreception is made possible by six distinct classes of multi-transmembrane receptors: (i) olfactory receptors (ORs) , (ii) trace amine-associated receptors , vomeronasal receptors (iii) type 1 and (iv) type 2 [7, 8] and taste receptors (v) type 1 and (vi) type 2 [9, 10]. Besides binding chemical molecules, all share the common traits of seven transmembrane (7-TM) domains, G-protein signaling and precise sensory cell expression. In mammals, non-volatile pheromone perception is thought to act primarily through the vomeronasal organ sensory epithelium  and be mediated intracellular via the interaction of chemical molecules with vomeronasal receptors located on the dendrites of vomeronasal sensory neurons . However, in teleost fishes who do not have a vomeronasal organ, the vomeronasal receptors are found in the main olfactory epithelium . It appears that genes involved in an animal's response to its environment are subject to extensive gene duplication, gene loss and lineage-specific expansion over time, leading to large gene families such as those observed in the OR and vomeronasal receptor repertoire. In fact, OR genes represent the largest mammalian gene family .
Chemoreception through 7-TM domain receptors appears to have evolved multiple times independently, as vertebrate chemoreceptors are not closely related to those known in insects and nematodes. Recognition of external chemicals inDrosophila is accomplished by families of 130 genes encoding 7-TM domain receptors [15, 16], including OR (60) and gustatory receptors (70). Gustatory receptors are greatly reduced in the honeybee . Insect chemoreceptors do not belong to the G-protein coupled receptor (GPCR) family due to a unique inverse membrane topology . Rather, they use an alternative, non-G protein-based signaling pathway where receptors not only detect chemicals but can also act as ion channels . In support of this, heterologous cells expressing silkmoth, fruitfly or mosquito heteromeric OR complexes showed G-protein independent extracellular calcium influx and cation-non-selective ion conductance upon stimulation with odorant . Nevertheless, chemical detection is still mediated by a large and divergent family of 7-TM domain receptors.
A central issue that has not been adequately addressed is how water-borne chemicals are detected at the molecular level by the huge diversity of invertebrates that inhabit marine environments. In marine invertebrates, chemosensory abilities are essential for almost all aspects of their life, from feeding to predator avoidance and reproduction. A recent bioinformatic survey of the sea urchin genome resulted in the identification of a remarkable diversity of chemoreceptors, expressed specifically and differentially in adult sensory structures . Meanwhile, there have been important findings forthcoming from research into the molluscan group. Olfactory studies of squid have shown that both phospholipase C (PLC) and cAMP-mediated pathways may be involved in olfactory sensory neurons activation . In support of this, immunolocalization experiments revealed the presence of G proteins involved in both cAMP (Gαo) and PLC (Gαq) pathways which are clearly co-expressed in certain cell types.
Aplysia possesses many advantages necessary for chemical communication research, such as an extensive knowledge of its anatomy, a detailed understanding of the molecular and cellular basis of behavior, and now considerable genomic and expressed sequence tag (EST) resources. Moreover, we have found that inAplysia, conspecific and congener attraction is mediated by a remarkable cocktail of water-borne protein pheromones [4, 22]. InAplysia, freshly laid egg cordons are considered to be a source of both water-borne and contact pheromones that attract conspecifics and closely related species to the area and induce them to mate and lay eggs. Egg laying results in the release of at least four proteinaceous attraction pheromones, including the 58-residue attractin [4, 23, 24]. T-maze bioassays have demonstrated that binary blends of attractin with either enticin, temptin or seductin are sufficient to attract potential mates .
At the anatomical level,Aplysia chemosensory detection is achieved by the rhinophore , specialized anterior sensory organs on the dorsal surface of the head. Rhinophore are retractile and primarily used for distance chemoreception and rheoreception (response to water current), whereas the oral tentacles, which are found more ventrally, are possibly involved in contact chemoreception and mechanoreception . The neuroanatomical organization of rhinophores includes a rhinophore groove where most of the sensory cells appear to be concentrated. Its sensory epithelium contains sensory neurons that project axons back to rhinophore ganglia and dendrites that end in either a surface-exposed cilium or a small protuberance [26–28]. Consistent with a potential role in chemical transduction, gene transcripts encoding G protein, PLC or inositol 1,4,5-trisphosphate receptor were found to be expressed inAplysia rhinophore sensory epithelium . The involvement of nitric oxide as a potential chemosensory processing component has also been implicated in molluscan chemoreception based on cytological NADH-diaphorase histochemistry of theAplysia rhinophore . Of significance, was the finding that nitric oxide synthase is present in epithelial sensory-like cells that had multiple apical ciliated processes exposed to the environment. This is consistent with findings demonstrating that inhibition of nitric oxide synthase disrupts slime trail following, suggesting a role for nitric oxide in neural processing of stimuli in snails .
The presence of G protein mRNAs inAplysia sensory epithelium suggested that multi-transmembrane GPCR-like proteins could play an important role in chemosensory detection. With the availability of a 2× genome coverage forAplysia californica, we expected that it would provide an excellent and first opportunity to investigate the molecular basis of chemical detection in a mollusc. Here, we performed iterative Basic Local Alignment Search Tool (BLAST) searches to identify genes similar to rhodopsin GPCR genes encoding 7-TM domains from theA. californica genome. We identified genes representing three unique monophyletic families that show rhinophore, oral tentacle and ovotestis expression. Based on their expression, these may encode chemosensory proteins, including pheromone and gustatory receptors. Antisera directed against a conserved region of a candidate chemosensory receptor, as well asAplysia Gαq, Gαi, and Gαo, confirmed their expression in sensory tissues, with localization to the outer sensory epithelium.
We performed iterative tBLASTn for closely related novel genes and discovered a large number of genes encoding rhodopsin GPCR-like proteins. Using this approach, we successfully identified a total of 90 genes encoding proteins belonging to the GPCR superfamily. Of these, 72 were predicted to contain 7-TM domains. It was not possible to annotate the full-length sequence of all genes, especially in the 5'-regions, and 18 genes encoding six transmembrane domains were considered partial-length. Note that these numbers are the minimal estimates, because the genome sequencing ofA. californica had not been fully completed. We expect that further rhodopsin GPCR-like genes will be found beyond the ones we have identified. Also, several of the multi-transmembrane gene models appeared to be pseudogenes with various defects, including the insertion of stop codons and frame-shifting indels leading to premature termination of the coding region. These were not included in the final data set.
Besides structural predictions that place these genes within the 7-TM superfamily, which includes rhodopsin GPCRs, these genes have little amino acid identity (<10%) to any known genes, and do not appear in the publishedAplysia EST neuronal transcriptome . They show only distant similarity with known molluscan multi-transmembrane receptors, including well-characterizedAplysia neurotransmitter GPCRs. GenBank tBLASTn searches also reveal most amino acid identity with regions of orphan GPCRs of the sea urchin,Strongylocentrotus purpuratus (E value 1e-11), various ghrelin receptors (for example,Rattus norvegicus, E value 2e-07) and a candidate GPCR (Caenorhabditis elegans, E value 4e-04) for subfamilies a, b and c, respectively. Based on the gene characteristics described below and observed tissue distribution (see Results, Tissue specificity of expression), we subsequently called them candidateAplysia californica chemosensory receptors (AcCRs) subfamilies a to c.
PCR amplification of a rhinophore epithelium LCM library using degenerate primers that were selective for AcCRb sequences (primer combinations B1 to B8) generated amplification products of 744 bp. Several amplicons were successfully cloned and sequenced, revealing multiple partial-length AcCRb genes. Subsequently, the full length of one gene sequence was identified by 5'- and 3'-RACE, containing 1483 bp and encoding a protein of 354 amino acids (Figure 7b – corresponding most closely, 99.4%, to geneAcCR29b). The sequence data has been submitted to the GenBank database under accession number EU808014. It possesses potentially six PKC phosphorylation sites (S11SK, T15QK, T147PR, T249NR, S322SK, S343ER) and four N-linked glycosylation sites (N5NS, N65IT, N184VT and N310ST) within the predicted N-terminal region and intracellular loop domains. Kyte-Doolittle hydropathy profiles of the deduced amino acid sequence showed that it contained seven hydrophobic transmembrane segments that were each composed of 25 residues.
AcCRc genes could not be PCR-amplified from a rhinophore LCM library. However, transcripts could be obtained from oral tentacle cDNA preparations (Figure 7c – corresponding most closely, 98%, to identified geneAcCR2c). PCR amplification of oral tentacle cDNAs using degenerate primers selective for AcCRc (gene combination C1) generated an amplification product of 824 bp. The amplicon was successfully cloned and sequenced. Subsequently, the full-length gene sequence was identified by 5'- and 3'-RACE, containing 1752 bp and encoding a protein of 398 amino acids. The sequence data has been submitted to the GenBank database under accession number EU808013. It possesses 10 PKC phosphorylation sites (T7ER, T28LR, T150FK, T188TR, S195SK, S275RR, S295NK, S308AK, T332SR, S385YR, a cAMP phosphorylation site (K234KSS) and six N-linked glycosylation sites (N5ET, N16IS, N54IT, N187TT, N280IS, N329TS) within the N-terminal region and intracellular loop domains. Kyte-Doolittle hydropathy profiles of the deduced amino acid sequence showed that it contained seven hydrophobic transmembrane segments that were composed of between 20 and 25 residues.
Summary of candidateAplysia chemoreceptor genes selected by genome mining
Number of intact genes
Size range (aa's)
Predicted TM domains
Sites of expression**
rhino, ot, ovo
Consistent with known chemosensory receptors (for example, insect, rodent, fish), all selected genes that were considered full length encoded 7-TM regions and semi-conserved glycosylation sites, as well as several common cysteine residues and amino acid sequence motifs. Conserved amino acids and post-translational modifications would likely contribute to the correct folding and functioning within the plasma membrane so that they may bind chemical stimuli and couple to appropriate secondary signaling molecules. For example, Katada and colleagues  demonstrated that in rodents, N-terminal glycosylation is critical for proper targeting of ORs to the plasma membrane. While proteins encoded by AcCRa and AcCRb genes share notably high sequence identity, comparative analysis within AcCRc shows that they share as little as 19% amino acid identity. As a consequence, there are very few defining sequence motifs which are retained throughout. However, it is not uncommon for chemoreceptors, and in particular gustatory receptors, to be divergent; similarity between most insect and mammalian gustatory receptor pairs is only 15% and 25% or less at the amino acid sequence level, respectively [10, 36].
Of the major GPCR superfamily groups, the identifiedAplysia genes categorize most closely to the rhodopsin GPCR family (based on Interpro database ). As is the case with many other rhodopsin family GPCRs, these genes largely lack introns. Moreover, all encoded proteins have a short N-terminus and a highly conserved arginine (R) residue located at the cytosolic end of the third transmembrane domain. This residue is typically associated with the DRY (Asp-Arg-Tyr) motif, crucial for controlling agonist-dependent receptor activation. Of the three residues constituting DRY, arginine is the most conserved residue, and appears to be essential for forming intramolecular interactions that constrain receptors in either the inactive or activated conformation . Consistent with this, receptors lacking the arginine side chain fail to activate G-protein signaling [38, 39]. In the novelAplysia proteins, however, this has been replaced by an ERC motif, a feature also observed in the human prostaglandin F2α receptor and most other prostanoid receptors . Studies show that substitution of the glutamic acid to a threonine residue leads to full constitutive activation and implicates the region in agonist-dependent G-protein coupling control . We predict that this may also be essential to receptor activation in the identifiedAplysia receptors.
The canonical model of GPCR activation is via an interaction with intracellular heterotrimeric G-protein signaling components. The genes identified in this study show little amino acid identity to GPCRs found in the Metazoa and we have not shown that they directly interact with G-proteins. Despite this, our study indicates that Gα proteins are present in the rhinophore sensory epithelium, possibly in close association with transmembrane GPCRs. The presence of sensory tissue G-protein immunoreactivity adds further support to studies of other marine invertebrate olfactory systems implicating G-proteins in sensory transduction [21, 41]. Moreover, the existence of multiple G-type proteins in sensory epithelium suggests that multiple signal transduction pathways may be activated following ligand stimulus. In squid, for example, the pattern of immunolabeling implies that a G protein coupled to a PLC pathway (Gαq and Gαo) may be present in similar cells as those coupled to a cAMP pathway (Gαi). As suggested by Mobley and colleagues , overlapping G-protein pathways could facilitate discrimination between odorants detected by the same neuron. This contrasts rodent models where the role of G proteins in olfactory transmembrane signaling at the dendrites has been studied extensively. Researchers have demonstrated spatially restricted patterns of expression of respective G proteins [42–44]. Gαo and Gαi2 are highly expressed by separate subsets of neurons that are located in different regions of the vomeronasal neuroepithelium .
Of particular relevance toAplysia chemosensory studies is the rhinophore epithelium, where water-borne molecules such as pheromones presumably bind and initiate activation of pheromone-receptive neurons. In the rhinophore groove, receptor cells with a suspected chemosensory role have been described previously in molluscs [26–28] and their presence was further supported by our SEM analysis. The rhinophore groove ciliary aggregations are likely necessary in the separation and circulation of fluid throughout the groove, and may also be directly involved in detection of external chemical stimuli. It is from this precise location that we isolated cDNAs encoding identified novel GPCR-like proteins. Their presence raised the question as to whether their expression was specific to sensory tissues. Subsequently, representatives of each subfamily were found to be restricted to or differentially expressed in the rhinophore, oral tentacles and ovotestis, suggesting that they encode functional receptors and that these olfactory organs sense different chemicals. In the rhinophore, spatial expression of cells immunoreactive to candidate chemoreceptor AcCR29b was most prominent in the tip and outer epithelium, peripheral to the groove. Chemosensory detection could likely benefit from this broad distribution, whereby stimulation may activate sensory fibers that extend to higher brain centers. Although this finding clearly indicates a sensory role, a more extensive study ofAplysia sensory organs at higher magnification is required to delineate the precise distribution of this receptor, as well as other receptor subfamily members.
Interestingly, we found gene expression within the ovotestis, and our preliminary analysis of variousAplysia neuronal EST databases indicate that a relatively small fraction of these genes may be expressed in the central nervous system. Deep sequencing of neuronal transcripts has resulted in identification of tags for 13 different genes in the central ganglia ofA. californica (that is,AcCR1a,AcCR5a,AcCR15a,AcCR32b,AcCR5c,AcCR9c,AcCR11c,AcCR15c,AcCR16c andAcCR20c, see Additional file 1) as well as several orphan receptors similar to vertebrate ghrelin and histamine receptors (L Moroz, unpublished). Some of those are associated with centrally located sensory neurons that send neuronal processes to the periphery and therefore may be involved in chemoreception. Other neuronal cell types were previously described as motorneurons. Taken together, these findings imply that external chemical detection may not be their sole function, which is consistent with that described of chemoreceptors in mammalian olfactory bulb , cardiac muscle  and vertebrate germ tissues [47–49]. The functional significance of such expression is currently unknown. We speculate that there may be a role forAplysia chemoreceptors in oocyte recognition, possibly because both are activated by identical or structurally similar hormones and pheromones. In addition to the exocrine albumen gland, theAplysia pheromone seductin is known to be expressed in ovotestis tissue .
It will be of interest to perform a comparative gene analysis to determine whetherA. californica candidate chemoreceptors are also present in otherAplysia species. This finding would suggest that these genes are highly similar throughoutAplysia species and strongly imply a selective pressure for conservation. We have already established that severalAplysia species share a comparable attraction pheromone blend [4, 23], and therefore cognate receptor binding sites are likely to be similar. The identification of ligands for chemosensory receptors is often problematic and therefore it is an advantage that we have these attraction pheromones for use in future functional studies to definitively link identified genes to chemosensory detection. As these are the first such novel GPCRs identified in molluscs, it will also be of interest to see if analogous receptors are found in evolutionarily distant molluscs. However, apart from the highly conserved insect receptor Or83b , it is generally accepted that chemoreceptors seem to be very divergent with little sequence conservation within and across orders [51, 52].
Our analysis of theAplysia genome noted that AcCRa genes as well as various AcCRb genes are clustered, a common feature of fast-evolving genes such as chemoreceptor genes [15, 53–55]. We also found that, although the assembly of the genome used was incomplete, some genes contain mutations that introduce stop codons to encode truncated proteins, one of which appears to be expressed in theA. californica rhinophore. Hominoids, in particular, are known to possess a high pseudogene content (50%) among their ORs, whereas only 20% of OR genes are pseudogenes in the mouse  and less in theDrosophila melanogaster genome . Upon genome completion, a more comprehensive analysis of GPCR gene families inAplysia will be necessary to determine the precise pseudogene number. Indeed, some of the pseudogenes identified may in fact be 'flatliners', that is, genes whose functional versus pseudogene status is unclear . As demonstrated inC. elegans, many of these genes have apparently functional alleles in one or more wild isolates and therefore are not pseudogenes. Evidence for this has also been shown for someDrosophila ORs  and gustatory receptors , as well asAnopheles gambiae gustatory receptors . Although pseudogenes are generally accepted as nonfunctional and therefore not transcribed, occasionally it has been shown that such pseudogenes can be transcribed ; however, there is no evidence of the functional relevance.
Aplysia is an excellent model animal for studying the molecular mechanism of chemical communication in the marine environment. In this study we have isolated a novel group of genes encoding multi-transmembrane rhodopsin GPCR-like proteins that show expression in chemosensory tissues. This expression pattern and observed genomic clustering provide strong evidence that these have arisen via gene duplication and may be used to discriminate the large diversity of water-soluble molecules. The expression of some of these in rhinophore suggests that they are excellent candidates to be involved in pheromone detection. Further knowledge of the receptor gene genome organization, characterization of their developmental and spatial expression profile, secondary signaling and their evolutionary relationship to other molluscan species would be the next significant steps towards defining the logic behind how chemical communication in molluscs, and potentially other marine animals, operates.
A. californica genome contig sequences were procured from the NCBI trace database http://www.ncbi.nlm.nih.gov/sutils/genom_table.cgi?organism=euk. This was a preliminary genome assembly from 2× coverage of the genome by the Broad Institute at MIT http://www.genome.gov/Pages/Research/Sequencing/SeqProposals/AplysiaSeq.pdf. An iterative tBLASTn strategy was adopted to identify multi-transmembrane rhodopsin GPCR-like genes in theAplysia genome. Selection criteria included that receptors would be encoded by a family of related genes; at least some receptor genes would be clustered at the same genetic loci; receptors would have limited amino acid identity to members of known GPCR superfamilies; and a full-length coding region would encode multiple transmembrane domains. A search was initially performed using molluscan GPCR protein sequences already submitted to GenBank, as queries; a non-stringent expectation value cutoff of 1e-4 was employed. During this search we retrieved two sequences on contig AASC01105652 with genes encoding hydrophobic multi-transmembrane domains with no significant amino acid identity to other proteins. PutativeAplysia transmembrane receptors were in turn employed in searches to find more genes in an iterative tBLASTn process. A candidate rhodopsin GPCR-like gene having a complete open reading frame (methionine, 7-TM domains, three extracellular domains, three intracellular domains, and a stop codon) was considered intact and probably functional. To be conservative, genes that were >98% identical in amino acid sequence were considered allelic variants. As theAplysia genome has not yet been fully assembled and consists of only contigs, this method does not identify splice-variants or provide a comprehensive analysis of gene clustering. Pseudogenes were identified by premature stop codons and we have arbitrarily chosen to disregard those genes that encode less than six transmembrane domains. Computer analyses of sequences were performed using BLAST and CLUSTALW for nucleotide alignment. Transmembrane helix domain and topology of predicted receptors was performed using HMMTOP version 2.0 program; http://www.enzim.hu/hmmtop/index.html. Hydrophilicity plots were generated using the TMHMM program at http://www.enzim.hu/hmmtop/html/submit.html. In order to categorize identified genes, we used http://pfam.janelia.org/search and the Interpro database .
Selection of outgroup sequences was performed by GenBank tBLASTn searches. The approach provided homologous counterparts fromC. elegans for theAplysia rhodopsin GPCR-like genes representing subfamilies a and b. We did not find any significant counterparts for subfamily c. The sequences were aligned with t-coffee  under default settings. Alignments have subsequently been improved by eye. Character positions, which could not be aligned unambiguously, were not considered for the phylogenetic analyses in order to avoid conflicting phylogenetic signal. The inclusion of outgroup sequences ofC. elegans andStrongylocentrotus purpuratus resulted in a large proportion of only ambiguously alignable positions and consequently required deletion of a relatively high number of characters. Therefore we performed two separate analyses: first, we aligned and analyzed the ingroup sequences only in order to reconstruct their evolution under a maximal number of informative characters. Second, we aligned and analyzed the dataset with the non-Aplysia sequences included in order to infer the polarity to the phylogeny.
Phylogenetic analyses were conducted on a multi-processor Linux-Cluster under the likelihood criterion using RAxML v. 7.0  for Maximum Likelihood and MrBayes v. 3.1.2  for the Bayesian Inference. MrBayes analyses were performed in two runs of eight MCMCMC chains and under the GTR+G+I Model . Chains ran for 10,000,000 generations or were stopped when the standard deviation of split frequencies between both runs fell below 0.01. RAxML bootstrap analyses on 1,000 replicates have been performed under the PROTMIX algorithm with the WAG amino acid substitution model . In PROTMIX the tree inference is performed under the PROTCAT model followed by the final tree evaluation under the PROTGAMMA model in order to obtain stable likelihood values (see the RAxML manual for further details).
AdultA. californica (100 to 500 g) were obtained from Marine Research and Educational Products (Escondido, CA, USA). Animals were anesthetized in isotonic MgCl2 (337 mM) equivalent to 50% of their weight, relevant tissues dissected out and either (1) embedded in optimal cutting temperature (OCT) compound for LCM or sectioning, (2) snap frozen in liquid nitrogen for RNA and protein isolation, or (3) prepared for SEM.
Total RNA was extracted from rhinophore, oral tentacle, skin, pooled central nervous system, albumen gland, large hermaphroditic duct, and ovotestis tissues ofA. californica using a Tripure Isolation Reagent. Any contaminating genomic DNA was removed by treatment with DNase I followed by lithium chloride/ethanol precipitation. First-strand cDNA synthesis was performed in a 20 μl reverse transcription mixture containing oligo d(T)12–18 and 200 U Superscript III RNase H- reverse transcriptase, following the manufacturer's instructions. PCR was performed using 1 μl of prepared cDNA using subfamily-specific primers (see Table S1 of Additional file 3: primer combinations A3, B3 and C1). Each reaction was performed in a final concentration of 1× PCR Buffer, 1.5 mM MgCl2, 200 μM of each dNTP, 0.5 μM of sense and antisense primer, 1.25 units of RedTaq polymerase and ddH2O. Negative controls contained no template cDNA. PCR using actin-specific primers (sense, 5'-GCTTCACCACCACTGCCGAGAG-3' and antisense, 5'-ACCAGCAGATTCCATACCCAGG-3') were used to ensure the integrity of each tissue cDNA sample. Reactions were heated at 94°C for 2 min and amplified for 36 cycles (94°C, 60 s; 50 to 55°C, 30 s; 72°C, 60 s). Following PCR, 15 μl of reaction mix was fractionated on 2% agarose gels and visualized by ethidium bromide staining. Based on the primer design, the expected amplification sizes were 756 bp (AcCRa), 832 bp (AcCRb), 512 bp (AcCRc) and 224 bp (actin). PCR products were cloned into pGEM-T vector and sequenced.
AdultAplysia rhinophore and oral tentacle were fixed with 2.5% glutaraldehyde in phosphate buffer (pH 7.2 to 7.4) for 3 days at 4°C. Secondary fixation was in 1% OsO4 (osmium tetroxide) in 0.1M sodium cacodylate. This material was dehydrated in a graded series of ethanol (20% to 100%). The samples were dried using hexamethyldisilazane and platinum-coated using Eiko IB-5 Sputter Coater. The specimens were viewed using a Jeol 6300 Field Emission Scanning Electron Microscope.
The location of the rhinophore groove, glomeruli underlying the sensory epithelium, and rhinophore ganglia inAplysia have been described previously [26–28]. To examine whether selected genes are expressed in rhinophore sensory epithelial cells, a combination of LCM, total RNA isolation, library construction and RT-PCR were performed. Rhinophore tissue that had been embedded in OCT compound was sectioned (10 μm) onto slides and dehydrated. LCM was performed using a PixCell II laser capture microscope with an infrared diode laser (Arcturus Engineering Inc., Mountain View, CA, USA) and a laser spot size of 15 μm. Cells were marked, captured on CapSure HS caps (Arcturus), and total RNA was isolated using the Picopure Isolation Kit (Arcturus) including DNase I incubation. RNA quality was determined by measurement of absorbance ratio at 260 nm/280 nm.
First-strand cDNA was synthesized according to the SMART cDNA Library Construction Kit protocol (Clontech, Palo Alto, CA, USA), with the minor modification of incorporatedEcoRI restriction sites during cDNA synthesis (5'-AAGCAGTGGTATCAACGCAGAGTGAATTCACGCGGG-3' and 5'-AAGCAGTGGTATCAACGCAGAGTGAATTCT30VN-3'). Amplification of cDNA was performed by PCR using 5' and 3' PCR Primer mix (5'-CTAATACGACTCACTATAGGGCAAGCAGTGGTATCAACGCAGAGT-3' and 5'-AAGCAGTGGTATCAACGCAGAGTGAATTCT30VN-3'); samples were heated at 95°C for 1 min and amplified for 30 cycles (95°C, 1 min; 50°C, 30 s; 68°C, 4 min). Reaction volumes of 50 μl were treated with 2 μl of proteinase K (20 μg/μl) at 45°C for 20 min, and PCR products were precipitated. Dried samples were resuspended in 80 μl of deionized water.EcoRI digestion was performed and size fractionation was achieved using a CHROMA SPIN-400 column (Stratagene, La Jolla, CA, USA). Products were purified by precipitation and the dried pellet resuspended in 7 μl deionized water.EcoRI cDNA was cloned into theEcoRI sites of Lambda ZAP II vector, and the cDNA library (complexity 1 × 106) amplified once.
The sequences of oligonucleotide primers (Sigma-Genosys, Australia) used for library PCR are located in Table S1 of Additional file 3. For AcCRa and AcCRb genes, PCR was performed using degenerate sense and antisense primer combinations A1 to A5 (AcCRa) and B1 to B8 (AcCRb). Each reaction was performed using RedTaq polymerase (Sigma) following the manufacturer's instructions. Samples were heated at 94°C for 3 min and amplified for 36 cycles (94°C, 60 s; 45°C, 30 s; 72°C, 60 s), followed by a 7-min extension at 72°C. PCR products were cloned into the TA vector pGEM-T (Promega) and sequenced as previously described . To obtain 5' and 3' sequences, PCR was performed using gene-specific primers. For 3'-RACE, gene-specific sense primers (A3' and B3') were used in combination with vector primer T3. For 5'-RACE, gene-specific antisense primers (A5' and B5') were used in combination with vector primer T7. Samples were heated at 94°C for 2 min and amplified for 36 cycles (94°C, 60 s; 50°C, 30 s; 72°C, 60 s), followed by a 7-min extension at 72°C. PCR products were cloned into pGEM-T vector and sequenced.
Total RNA was extracted from oral tentacle tissue ofA. californica using a Tripure Isolation Reagent (Roche), and any contaminating genomic DNA was removed by treatment with DNase I (Invitrogen). First strand cDNA synthesis was performed using 1 μg of total RNA in a 20 μl reverse transcription mixture containing oligo d(T)12–18 and 200 U Superscript™ III RNase H- reverse transcriptase, following the manufacturer's instructions. The sequences of oligonucleotide primers used for PCR are located in Table S1 of Additional file 3. PCR was performed using degenerate sense and antisense primer combinations C1-C3. 3'- and 5'-RACE was performed with the SMART RACE amplification kit (BD Biosciences) and using gene-specific primers sense and antisense primers (C3' and C5'). Samples were heated at 94°C for 2 min and amplified for 36 cycles (94°C, 60 s; 50°C, 30 s; 72°C, 60 s), followed by a 7-min extension at 72°C. PCR products were cloned into a pGEM-T vector and sequenced.
A rabbit polyclonal antibody was generated to the N-terminal region of candidate chemoreceptor 29b, corresponding to N6SQARSSKSTQKGL (GenScript Corporation). This region was chosen due to its lack of significant amino acid identity to other receptors. Details of the immunohistochemical protocol have been described [23, 24]. Briefly, tissue cryostat sections of rhinophore were incubated overnight at 4°C in either affinity-purified 29b antibody (0.6 mg/ml, 1:500), Gαq, Gαi or Gαo antisera (Chemicon, 1:500 dilution), rinsed in phosphate buffered saline (PBS), incubated in fluorescein isothiocyanate (FITC)-conjugated goat anti-rabbit Ig (Sigma-Aldrich, St. Louis, MO, USA) for 1 h at 22°C, rinsed in PBS, and then mounted in FITC mounting media (90% glycerol/100 mM Tris pH 8.0). Preparations were examined using an Olympus FluoView confocal microscope (Leeds Precision Instruments, Inc., Minneapolis, USA), and the images captured on a spot-cooled charged coupled device camera. Sections were counterstained with 4',6-diamidino-2-phenylindole or propidium iodide at 1 μg/ml in water. As a control, the primary antiserum was replaced with no primary antibody. For 29b, a control also included using the primary antibody that had been preabsorbed against its antigenic peptide (20 μg/ml).
Aplysia californica chemosensory receptor
basic local alignment search tool
expressed sequence tag
G-protein coupled receptor
laser capture microdissection
optimal cutting temperature
reverse transcription-polymerase chain reaction
We acknowledge the Broad InstituteAplysia Genome Initiative for making the partial genome sequence ofA. californica available. We thank Dr Darren Boehning for comments on an earlier version of the manuscript. We acknowledge the assistance of the UTMB Protein Chemistry Lab, Steve LePage (MREP) and Erica Lovus (Institute of Molecular Biology). SFC is supported by a University of Queensland Fellowship. This research was supported by grants to BMD from the Australian Research Council, LLM from NIH and NSF, and to GTN from the NSF (grant IBN-0314377).
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.