Candidate chemoreceptor subfamilies differentially expressed in the chemosensory organs of the mollusc Aplysia
© Cummins et al; licensee BioMed Central Ltd. 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 mollusc Aplysia 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 the Aplysia rhinophore. Analysis of the A. 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 novel Aplysia 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 in Aplysia involves species- or lineage-specific families of chemosensory receptors.
All animals must recognize and respond to chemosensory information in their environment. Although the marine mollusc Aplysia 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 in Drosophila 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 in Aplysia, conspecific and congener attraction is mediated by a remarkable cocktail of water-borne protein pheromones [4, 22]. In Aplysia, 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 in Aplysia 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 the Aplysia 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 in Aplysia 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 for Aplysia 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 the A. 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 as Aplysia Gαq, Gαi, and Gαo, confirmed their expression in sensory tissues, with localization to the outer sensory epithelium.
Identification of genes encoding rhodopsin G-protein coupled receptor-like proteins in the Aplysia californica genome
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 of A. 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.
Phylogenetic construction and analysis of identified rhodopsin G-protein coupled receptor-like genes
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 published Aplysia EST neuronal transcriptome . They show only distant similarity with known molluscan multi-transmembrane receptors, including well-characterized Aplysia 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 candidate Aplysia californica chemosensory receptors (AcCRs) subfamilies a to c.
Candidate Aplysia californicachemosensory receptors subfamily a
Candidate Aplysia californicachemosensory receptors subfamily b (AcCRb)
Candidate Aplysia californicachemosensory receptors subfamily c (AcCRc)
Tissue specificity of expression
Scanning electron microscopy and molecular identification of candidate chemoreceptors
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 gene AcCR29b). 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 gene AcCR2c). 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.
Immunofluorescent localization of a candidate chemoreceptor
Distribution of Gαq, Gαi and Gαoimmunoreactivity in rhinophore
Summary of candidate Aplysia 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 identified Aplysia 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 novel Aplysia 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 identified Aplysia 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 to Aplysia 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 of Aplysia 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 various Aplysia 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 of A. californica (that is, AcCR1a, AcCR5a, AcCR15a, AcCR32b, AcCR5c, AcCR9c, AcCR11c, AcCR15c, AcCR16c and AcCR20c, 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 for Aplysia chemoreceptors in oocyte recognition, possibly because both are activated by identical or structurally similar hormones and pheromones. In addition to the exocrine albumen gland, the Aplysia pheromone seductin is known to be expressed in ovotestis tissue .
It will be of interest to perform a comparative gene analysis to determine whether A. californica candidate chemoreceptors are also present in other Aplysia species. This finding would suggest that these genes are highly similar throughout Aplysia species and strongly imply a selective pressure for conservation. We have already established that several Aplysia 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 the Aplysia 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 the A. 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 the Drosophila melanogaster genome . Upon genome completion, a more comprehensive analysis of GPCR gene families in Aplysia 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 in C. 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 some Drosophila ORs  and gustatory receptors , as well as Anopheles 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.
Database mining for multi-transmembrane rhodopsin G-protein coupled receptor-like genes
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 the Aplysia 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. Putative Aplysia 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 the Aplysia 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 .
Phylogenetic and gene analysis
Selection of outgroup sequences was performed by GenBank tBLASTn searches. The approach provided homologous counterparts from C. elegans for the Aplysia 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 of C. elegans and Strongylocentrotus 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).
Animal and sample preparation
Adult A. 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.
Reverse transcription-polymerase chain reaction
Total RNA was extracted from rhinophore, oral tentacle, skin, pooled central nervous system, albumen gland, large hermaphroditic duct, and ovotestis tissues of A. 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 Red Taq 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.
Scanning electron microscopy of rhinophore and oral tentacle
Adult Aplysia 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.
Laser capture microdissection and molecular identification of candidate A. californicachemosensory receptors
The location of the rhinophore groove, glomeruli underlying the sensory epithelium, and rhinophore ganglia in Aplysia 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 incorporated EcoRI 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 the EcoRI 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 Red Taq 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.
Molecular identification of candidate Aplysia chemosensory receptors within subfamily c
Total RNA was extracted from oral tentacle tissue of A. 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 Institute Aplysia Genome Initiative for making the partial genome sequence of A. 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).
- Antonov I, Antonova I, Kandel ER, Hawkins RD: Activity-dependent presynaptic facilitation and hebbian LTP are both required and interact during classical conditioning in Aplysia. Neuron. 2003, 37: 135-147.View ArticlePubMedGoogle Scholar
- Balakrishnan R: Learning from a sea snail: Eric kandel. Resonance. 2001, 6: 86-90.View ArticleGoogle Scholar
- Nambu JR, Scheller RH: Egg-laying hormone genes of Aplysia: evolution of the ELH gene family. J Neurosci. 1986, 6: 2026-2036.PubMedGoogle Scholar
- Painter SD, Cummins SF, Nichols AE, Akalal DB, Schein CH, Braun W, Smith JS, Susswein AJ, Levy M, de Boer PA, ter Maat A, Miller MW, Scanlan C, Milberg RM, Sweedler JV, Nagle GT: Structural and functional analysis of Aplysia attractins, a family of water-borne protein pheromones with interspecific attractiveness. Proc Natl Acad Sci USA. 2004, 101: 6929-6933.PubMed CentralView ArticlePubMedGoogle Scholar
- Buck L, Axel R: A novel multigene family may encode odorant receptors: a molecular basis for odor recognition. Cell. 1991, 65: 175-187.View ArticlePubMedGoogle Scholar
- Liberles SD, Buck LB: A second class of chemosensory receptors in the olfactory epithelium. Nature. 2006, 442: 645-650.View ArticlePubMedGoogle Scholar
- Dulac C, Axel R: A novel family of genes encoding putative pheromone receptors in mammals. Cell. 1995, 83: 195-206.View ArticlePubMedGoogle Scholar
- Ryba NJ, Tirindelli R: A new multigene family of putative pheromone receptors. Neuron. 1997, 19: 371-379.View ArticlePubMedGoogle Scholar
- Li X, Staszewski L, Xu H, Durick K, Zoller M, Adler E: Human receptors for sweet and umami taste. Proc Natl Acad Sci USA. 2002, 99: 4692-4696.PubMed CentralView ArticlePubMedGoogle Scholar
- Matsunami H, Montmayeur JP, Buck LB: A family of candidate taste receptors in human and mouse. Nature. 2000, 404: 601-604.View ArticlePubMedGoogle Scholar
- Touhara K: Molecular biology of peptide pheromone production and reception in mice. Adv Genet. 2007, 59: 147-171.View ArticlePubMedGoogle Scholar
- Zufall F, Leinders-Zufall T: Mammalian pheromone sensing. Curr Opin Neurobiol. 2007, 17: 483-489.View ArticlePubMedGoogle Scholar
- Cao Y, Oh BC, Stryer L: Cloning and localization of two multigene receptor families in goldfish olfactory epithelium. Proc Natl Acad Sci USA. 1998, 95: 11987-11992.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang X, Firestein S: The olfactory receptor genesuperfamily of the mouse. Nat Neurosci. 2002, 5: 124-133.PubMedGoogle Scholar
- Clyne PJ, Warr CG, Freeman MR, Lessing D, Kim J, Carlson JR: A novel family of divergent seven-transmembrane proteins: candidate odorant receptors in Drosophila. Neuron. 1999, 22: 327-338.View ArticlePubMedGoogle Scholar
- Vosshall LB, Amrein H, Morozov PS, Rzhetsky A, Axel R: A spatial map of olfactory receptor expression in the Drosophila antenna. Cell. 1999, 96: 725-736.View ArticlePubMedGoogle Scholar
- Robertson HM, Wanner KW: The chemoreceptor superfamily in the honey bee, Apis mellifera: expansion of the odorant, but not gustatory, receptor family. Genome Res. 2006, 16: 1395-1403.PubMed CentralView ArticlePubMedGoogle Scholar
- Benton R, Sachse S, Michnick SW, Vosshall LB: Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLoS Biol. 2006, 4: e20-PubMed CentralView ArticlePubMedGoogle Scholar
- Sato K, Pellegrino M, Nakagawa T, Nakagawa T, Vosshall LB, Touhara K: Insect olfactory receptors are heteromeric ligand-gated ion channels. Nature. 2008, 452: 1002-1006.View ArticlePubMedGoogle Scholar
- Raible F, Tessmar-Raible K, Arboleda E, Kaller T, Bork P, Arendt D, Arnone MI: Opsins and clusters of sensory G-protein-coupled receptors in the sea urchin genome. Dev Biol. 2006, 300: 461-475.View ArticlePubMedGoogle Scholar
- Mobley AS, Mahendra G, Lucero MT: Evidence for multiple signaling pathways in single squid olfactory receptor neurons. J Comp Neurol. 2007, 501: 231-242.View ArticlePubMedGoogle Scholar
- Cummins SF, Schein CH, Xu Y, Braun W, Nagle GT: Molluscan attractins, a family of water-borne protein pheromones with interspecific attractiveness. Peptides. 2005, 26: 121-129.View ArticlePubMedGoogle Scholar
- Cummins SF, Nichols AE, Amare A, Hummon AB, Sweedler JV, Nagle GT: Characterization of Aplysia enticin and temptin, two novel water-borne protein pheromones that act in concert with attractin to stimulate mate attraction. J Biol Chem. 2004, 279: 25614-25622.View ArticlePubMedGoogle Scholar
- Cummins SF, Nichols AE, Warso CJ, Nagle GT: Aplysia seductin is a water-borne protein pheromone that acts in concert with attractin to stimulate mate attraction. Peptides. 2005, 26: 351-359.View ArticlePubMedGoogle Scholar
- Levy M, Blumberg M, Susswein AJ: The rhinophores sense pheromones regulating multiple behaviors in Aplysia fasciata. Neurosci Lett. 1997, 225: 113-116.View ArticlePubMedGoogle Scholar
- Gobbeler K, Klussmann-Kolb A: A comparative ultrastructural investigation of the cephalic sensory organs in Opisthobranchia (Mollusca, Gastropoda). Tissue Cell. 2007, 39: 399-414.View ArticlePubMedGoogle Scholar
- Emery DG, Audesirk TE: Sensory cells in Aplysia. J Neurobiol. 1978, 9: 173-179.View ArticlePubMedGoogle Scholar
- Wertz A, Rossler W, Obermayer M, Bickmeyer U: Functional neuroanatomy of the rhinophore of Aplysia punctata. Front Zool. 2006, 3: 6-PubMed CentralView ArticlePubMedGoogle Scholar
- Cummins SF, De Vries MR, Hill KS, Boehning D, Nagle GT: Gene identification and evidence for expression of G protein alpha subunits, phospholipase C, and an inositol 1,4,5-trisphosphate receptor in Aplysia californica rhinophore. Genomics. 2007, 90: 110-120.View ArticlePubMedGoogle Scholar
- Moroz LL: Localization of putative nitrergic neurons in peripheral chemosensory areas and the central nervous system of Aplysia californica. J Comp Neurol. 2006, 495: 10-20.View ArticlePubMedGoogle Scholar
- Clifford KT, Gross L, Johnson K, Martin KJ, Shaheen N, Harrington MA: Slime-trail tracking in the predatory snail, Euglandina rosea. Behav Neurosci. 2003, 117: 1086-1095.View ArticlePubMedGoogle Scholar
- Moroz LL, Edwards JR, Puthanveettil SV, Kohn AB, Ha T, Heyland A, Knudsen B, Sahni A, Yu F, Liu L, Jezzini S, Lovell P, Lannucculli W, Chen M, Nguyen T, Sheng H, Shaw R, Kalachikov S, Panchin YV, Farmerie W, Russo JJ, Ju J, Kandel ER: Neuronal transcriptome of Aplysia: neuronal compartments and circuitry. Cell. 2006, 127: 1453-1467.PubMed CentralView ArticlePubMedGoogle Scholar
- Mulder NJ, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Buillard V, Cerutti L, Copley R, Courcelle E, Das U, Daugherty L, Dibley M, Finn R, Fleischmann W, Gough J, Haft D, Hulo N, Hunter S, Kahn D, Kanapin A, Kejariwal A, Labarga A, Langendijk-Genevaux PS, Lonsdale D, Lopez R, Letunic I, Madera M, Maslen J, et al: New developments in the InterPro database. Nucleic Acids Res. 2007, 35: D224-228.PubMed CentralView ArticlePubMedGoogle Scholar
- Freitag J, Ludwig G, Andreini I, Rossler P, Breer H: Olfactory receptors in aquatic and terrestrial vertebrates. J Comp Physiol [A]. 1998, 183: 635-650.View ArticleGoogle Scholar
- Katada S, Tanaka M, Touhara K: Structural determinants for membrane trafficking and G protein selectivity of a mouse olfactory receptor. J Neurochem. 2004, 90: 1453-1463.View ArticlePubMedGoogle Scholar
- Scott K, Brady R, Cravchik A, Morozov P, Rzhetsky A, Zuker C, Axel R: A chemosensory gene family encoding candidate gustatory and olfactory receptors in Drosophila. Cell. 2001, 104: 661-673.View ArticlePubMedGoogle Scholar
- Rosenkilde MM, Kledal TN, Schwartz TW: High constitutive activity of a virus-encoded seven transmembrane receptor in the absence of the conserved DRY motif (Asp-Arg-Tyr) in transmembrane helix 3. Mol Pharmacol. 2005, 68: 11-19.PubMedGoogle Scholar
- Acharya S, Karnik SS: Modulation of GDP release from transducin by the conserved Glu134-Arg135 sequence in rhodopsin. J Biol Chem. 1996, 271: 25406-25411.View ArticlePubMedGoogle Scholar
- Ballesteros J, Kitanovic S, Guarnieri F, Davies P, Fromme BJ, Konvicka K, Chi L, Millar RP, Davidson JS, Weinstein H, Sealfon SC: Functional microdomains in G-protein-coupled receptors. The conserved arginine-cage motif in the gonadotropin-releasing hormone receptor. J Biol Chem. 1998, 273: 10445-10453.View ArticlePubMedGoogle Scholar
- Pathe-Neuschäfer-Rube A, Neuschäfer-Rube F, Püschel GP: Role of the ERC motif in the proximal part of the second intracellular loop and the C-terminal domain of the human prostaglandin F2alpha receptor (hFP-R) in G-protein coupling control. Biochem J. 2005, 388: 317-324.PubMed CentralView ArticlePubMedGoogle Scholar
- Wodicka LM, Morse DE: cDNA sequences reveal mRNAs for two Gα signal transducing proteins from larval cilia. Biol Bull. 1991, 180: 318-327.View ArticleGoogle Scholar
- Shinohara H, Asano T, Kato K: Differential localization of G-proteins Gi and Go in the accessory olfactory bulb of the rat. J Neurosci. 1992, 12: 1275-1279.PubMedGoogle Scholar
- Jia C, Goldman G, Halpern M: Development of vomeronasal receptor neuron subclasses and establishment of topographic projections to the accessory olfactory bulb. Brain Res Dev Brain Res. 1997, 102: 209-216.View ArticlePubMedGoogle Scholar
- Berghard A, Buck L: Sensory transduction in vomeronasal neurons: evidence for G alpha o, G alpha i2, and adenylyl cyclase II as major components of a pheromone signaling cascade. J Neurosci. 1996, 16: 909-918.PubMedGoogle Scholar
- Vassar R, Chao SK, Sitcheran R, Nunez JM, Vosshall LB, Axel R: Topographic organization of sensory projections to the olfactory bulb. Cell. 1994, 79: 981-991.View ArticlePubMedGoogle Scholar
- Ferrand N, Pessah M, Frayon S, Marais J, Garel JM: Olfactory receptors, Golf alpha and adenylyl cyclase mRNA expressions in the rat heart during ontogenic development. J Mol Cell Cardiol. 1999, 31: 1137-1142.View ArticlePubMedGoogle Scholar
- Parmentier M, Libert F, Schurmans S, Schiffmann S, Lefort A, Eggerickx D, Ledent C, Mollereau C, Gerard C, Perret J, Grootegoed A, Vassart G: Expression of members of the putative olfactory receptor gene family in mammalian germ cells. Nature. 1992, 355: 453-455.View ArticlePubMedGoogle Scholar
- Fukuda N, Yomogida K, Okabe M, Touhara K: Functional characterization of a mouse testicular olfactory receptor and its role in chemosensing and in regulation of sperm motility. J Cell Sci. 2004, 117: 5835-5845.View ArticlePubMedGoogle Scholar
- Branscomb A, Seger J, White RL: Evolution of odorant receptors expressed in mammalian testes. Genetics. 2000, 156: 785-797.PubMed CentralPubMedGoogle Scholar
- Krieger J, Klink O, Mohl C, Raming K, Breer H: A candidate olfactory receptor subtype highly conserved across different insect orders. J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2003, 189: 519-526.View ArticlePubMedGoogle Scholar
- Hill CA, Fox AN, Pitts RJ, Kent LB, Tan PL, Chrystal MA, Cravchik A, Collins FH, Robertson HM, Zwiebel LJ: G protein-coupled receptors in Anopheles gambiae. Science. 2002, 298: 176-178.View ArticlePubMedGoogle Scholar
- Go Y, Niimura Y: Similar numbers but different repertoires of olfactory receptor genes in humans and chimpanzees. Mol Biol Evol. 2008, 25: 1897-1907.View ArticlePubMedGoogle Scholar
- Bohbot J, Pitts RJ, Kwon HW, Rutzler M, Robertson HM, Zwiebel LJ: Molecular characterization of the Aedes aegypti odorant receptor gene family. Insect Mol Biol. 2007, 16: 525-537.PubMed CentralPubMedGoogle Scholar
- Troemel ER, Chou JH, Dwyer ND, Colbert HA, Bargmann CI: Divergent seven transmembrane receptors are candidate chemosensory receptors in C. elegans. Cell. 1995, 83: 207-218.View ArticlePubMedGoogle Scholar
- Xie SY, Feinstein P, Mombaerts P: Characterization of a cluster comprising approximately 100 odorant receptor genes in mouse. Mamm Genome. 2000, 11: 1070-1078.View ArticlePubMedGoogle Scholar
- Nozawa N, Masatochi N: Evolutionary dynamics of olfactory receptor genes in Drosophila species. Proc Natl Acad Sci USA. 2007, 104: 7122-7127.PubMed CentralView ArticlePubMedGoogle Scholar
- Stewart MK, Clark NL, Merrihew G, Galloway EM, Thomas JH: High genetic diversity in the chemoreceptor superfamily of Caenorhabditis elegans. Genetics. 2005, 169: 1985-1996.PubMed CentralView ArticlePubMedGoogle Scholar
- McBride C: Rapid evolution of smell and taste receptor genes during host specialization in Drosophila sechellia. Proc Natl Acad Sci USA. 2007, 104: 4996-5001.PubMed CentralView ArticlePubMedGoogle Scholar
- Robertson HM, Warr CG, Carlson JR: Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster. Proc Natl Acad Sci USA. 2003, 100 (Suppl 2): 14537-14542.PubMed CentralView ArticlePubMedGoogle Scholar
- Liman ER, Innan H: Relaxed selective pressure on an essential component of pheromone transduction in primate evolution. Proc Natl Acad Sci USA. 2003, 100: 3328-3332.PubMed CentralView ArticlePubMedGoogle Scholar
- NCBI trace database. [http://www.ncbi.nlm.nih.gov/sutils/genom_table.cgi?organism=euk]
- Aplysia genome sequencing, Broad Institute at MIT. [http://www.genome.gov/Pages/Research/Sequencing/SeqProposals/AplysiaSeq.pdf]
- HMMTOP v2. [http://www.enzim.hu/hmmtop/index.html]
- TMHMM program. [http://www.enzim.hu/hmmtop/html/submit.html]
- Pfam. [http://pfam.janelia.org/search]
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302: 205-217.View ArticlePubMedGoogle Scholar
- Stamatakis A: RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574.View ArticlePubMedGoogle Scholar
- Tavare S: Some probabilistic and statistical problems on the analysis of DNA sequences. Lect Math Life Sci. 1986, 17: 57-86.Google Scholar
- Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18: 691-699.View ArticlePubMedGoogle Scholar
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.