- Research article
- Open Access
Deep sequencing of transcriptomes from the nervous systems of two decapod crustaceans to characterize genes important for neural circuit function and modulation
BMC Genomicsvolume 17, Article number: 868 (2016)
Crustaceans have been studied extensively as model systems for nervous system function from single neuron properties to behavior. However, lack of molecular sequence information and tools have slowed the adoption of these physiological systems as molecular model systems. In this study, we sequenced and performed de novo assembly for the nervous system transcriptomes of two decapod crustaceans: the Jonah crab (Cancer borealis) and the American lobster (Homarus americanus).
Forty-two thousand, seven hundred sixty-six and sixty thousand, two hundred seventy-three contigs were assembled from C. borealis and H. americanus respectively, representing 9,489 and 11,061 unique coding sequences. From these transcripts, genes associated with neural function were identified and manually curated to produce a characterization of multiple gene families important for nervous system function. This included genes for 34 distinct ion channel types, 17 biogenic amine and 5 GABA receptors, 28 major transmitter receptor subtypes including glutamate and acetylcholine receptors, and 6 gap junction proteins – the Innexins.
With this resource, crustacean model systems are better poised for incorporation of modern genomic and molecular biology technologies to further enhance the interrogation of fundamentals of nervous system function.
Despite their status as important economic species, their important place in understanding the evolution and phylogeny of arthropods, and as models for neurobiology research, crustaceans have been largely overlooked in the rush to apply modern molecular biology and high throughput sequencing approaches to work in “non-genetic” systems. Even among arthropods they are fairly poorly represented, with insects dominating the ranks of those with available genome and transcriptome assemblies. Two indicators of this are seen in the number of relative sequence read archive (SRA) and GEO profile publicly available in NCBI: at the time of this article, there were 2,323 crustacean and 46,866 insect SRAs, and 4,608 crustacean and 1,275,029 insect GEO profiles. To date, only two crustacean genomes have been made publicly available [1, 2], with the first, the water flea Daphnia pulex, coming only in 2011  – a full 11 years after the first arthropod genome was sequenced . As a microcrustacean, Daphnia is still a far cry from the large decapod crustaceans that are common models in neuroscience research, such as crabs, lobsters, crayfish and shrimp. Thus even with some very recent additions to decapod crustacean transcriptome data [4, 5], there is a strong need to add to our sequence knowledge of these species.
Many fundamental findings in neuroscience were made with crustacean preparations. To mention only a small subset of these, command fibers , electrical coupling  and presynaptic inhibition  were first described using crustacean preparations. Work on crayfish and lobsters established GABA as an inhibitory transmitter [9, 10], and allowed early studies of the relevance of the fast outward current, IA, for action potential generation and propagation [11, 12]. The first intracellular fluorescent dye-fills were pioneered with crustaceans [13, 14], and crustacean systems were used early on to understand the organization of circuits in behavior [15–18].
Several crustacean circuits, including the stomatogastric nervous system and the cardiac ganglion, continue to provide important new insights into circuit dynamics and modulation [19–22], but this work has been partially hampered by the lack of extensive molecular sequence knowledge in crustaceans. In this study, we generated de novo transcriptome assemblies from central nervous system tissue for two commonly used species in neuroscience research: the Jonah crab (Cancer borealis) and the American lobster (Homarus americanus). We focus on an initial identification, curation and comparison of genes that will have the most profound impact on our understanding of circuit function in these species, namely channels and receptors, with the hopes of fostering new avenues of research for these preparations that continue to be valuable assets in our understanding of nervous system dynamics. Additionally, such sequence information provides a valuable resource for comparative molecular neuroscience approaches across phyla.
Tissue collection and RNA preparation
Adult lobsters, H. americanus, and crabs, C. borealis, were obtained from The Fresh Lobster Company (Gloucester, Massachusetts, USA) and maintained in chilled (12 °C) artificial seawater tanks until experiments were performed. Lobsters and crabs were anesthetized on ice for at least 30 minutes prior to dissection. The brain, abdominal nerve cord, cardiac ganglion and complete stomatogastric nervous system (STNS) (including the commissural, esophageal, and stomatogastric ganglia) was dissected out of two lobsters and pinned out in a Sylgard (Dow Corning)–coated dish containing chilled (12–13 °C) physiological saline. From two crabs we dissected out brain, complete STNS, and cardiac ganglia. Connective tissue and muscle were removed to the fullest extent possible, and the tissues were rinsed multiple times in chilled physiological saline (Lobster saline composition in mM/l: 479.12 NaCl, 12.74 KCl, 13.67 CaCl2, 20.00 MgSO4, 3.91 Na2SO4, 11.45 Trizma base, and 4.82 maleic acid [pH = 7.45]; Crab saline composition in mM/l: 440.0 NaCl, 11.0 KCl, 13.0 CaCl2, 20.00 MgCl2, 11.2 Trizma base, and 5.1 maleic acid [pH = 7.45]) in ultrapure, RNase-free water. After dissection, tissues for each species were homogenized in Trizol (Invitrogen). The resulting combined pool of RNA therefore consisted of mixed nervous system tissue. Insoluble tissues were pelleted by centrifugation, and the supernatant stored at −80 ° C until RNA extraction. Total RNA was isolated as per the manufacturer’s protocol (Invitrogen), and treated with DNase (Zymo Research) prior to library construction.
Library construction, sequencing, and de novo assembly
Library construction and RNA-sequencing were performed as a fee-for-service by GENEWIZ, Inc. (South Plainfield, New Jersey, USA). Briefly, quantification of RNA samples was performed using a Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, California, USA) and RNA quality checked with an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, California, USA). Illumina TruSeq RNA library prep, clustering, and sequencing reagents were used throughout the process as specified by the manufacturer (Illumina, San Diego, California, USA). mRNAs were purified using oligo-attached poly-T magnetic beads. The mRNAs were fragmented and first and second strand cDNAs were synthesized and end-repaired. cDNA templates were enriched by PCR following adaptor ligation after adenylation at the 3′ends. cDNA libraries were validated using an Agilent 2100 Bioanalyzer with a High Sensitivity Chip. cDNA library yield was quantified with a Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, California, USA) and by qPCR. After clustering on a flow cell using the cBOT, the samples were loaded on an Illumina HiSeq 2000 instrument for sequencing with 2x100 paired-end reads.
Raw reads were converted into fastq files and de-multiplexed using Illumina CASSAVA 1.8.2. Fastq files were imported into CLC Genomics Workbench Server 5.0.1. Sequence reads were trimmed to remove bases with low quality ends. De novo assembly was conducted with the trimmed reads utilizing the CLC Genomics Server. The total length of the assembled transcripts was 66,058,464 bp for crab and 99,847,148 bp for lobster (see Table 1). To ensure that the CLC Genomics assembly was of high quality, we performed a second round of de novo assembly using the SeqMan NGen assembler from the DNAstar software suite (SeqMan NGen®. Version 13.0. DNASTAR. Madison, WI.). Following assembly, quality of assembled contigs was investigated by comparison with species-specific sequences contained within GenBank that were previously generated largely by Sanger sequencing approaches. We used BLAST+ command line application (Version 2.2.30+) to perform blastn comparisons of the curated GenBank sequence versus transcriptome contigs, and calculate percent nucleotide identity for the top hit of each sequence for both species.
BUSCO transcriptome quality assessment
To analyze the completeness of our transcriptomes, a reference-based alignment was performed using Benchmarking Universal Single-Copy Orthologs (BUSCO) software (Version 1.22). The arthropod BUSCO reference contains 2675 orthologous genes found within >90 % of the 38 arthropod species’ genomes used to construct the reference . The four transcriptomes we assembled were aligned against the arthropod reference, resulting in percentages of the reference genes found as complete, fragmented, or missing from our transcriptomes. “Complete” genes are those which align to a reference gene with a mean length within two standard deviations (i.e. 95 %) of the reference value. Genes that only partially align are deemed ”fragmented”, and those present in the reference with no match found in the transcriptome are classified as “missing”.
For gene ontology (GO) term analysis, the Blast2GO software package (Version 3.1.3)  was used for functional annotation of the assembled transcriptomes. A blastx search with an E-value threshold of 10−5 was carried out against the NCBI non-redundant (nr) protein sequence database. Assignment of gene names to each contig was based on the highest scoring BLAST hit. Scoring of the annotated sequences utilized a threshold score of ≥ 55. The top 10 significant hits for each query extracted from the blastx search were used for further gene annotation. Query sequences were categorized into three broad ontological classifications: molecular function, cellular component, and biological process. GO annotation filters included: E-value-Hit-Filter of 1.0e-6, Annotation CutOff of 55, and GO Weight of 5.
Whole-transcriptome alignment comparison
The software VennBLAST  was used to compare the whole C. borealis and H. americanus transcriptomes against the Daphnia pulex (GCA_000187875.1) protein sequences from Ensembl Metazoa. Protein sequence database for D. pulex was chosen as a common subject to query against the C. borealis and H. americanus transcriptomes. Initially, a local blastx of C. borealis or H. americanus contigs against D. pulex protein sequences was performed with the BLAST+ command line application (Version 2.2.30+). This output was run through the VennBLAST Merge tool with the InterGroup Option: Use Subject to quantify the relative overlap of C. borealis and H. americanus with D. pulex. A second layer of filtering was performed using the VennBLAST Filter tool with an Identity percent of 70 and an E-value threshold of 1.0e-5, and this output was subsequently merged in the same manner mentioned previously.
Ion channel and receptor sequence identification and alignment analysis
We identified putative orthologs of channels and receptors from the transcriptomes of crab and lobster as follows. We created local blast databases from the assembled contigs of each transcriptomes. Because channels and receptors are fairly well conserved across diverse taxa, and because the mouse research community has agreed upon a well-curated systematic naming system for channel and receptor genes, we used mouse reference mRNA sequence for each gene of interest as the query in a tblastn search of each transcriptome database. We used initially stringent e-value cutoffs (1e-100 to 1e-50) for our searches to find very high sequence similarity matches. Top contig matches from these blast searches were then compared with the results of the remainder of the blast queries for a given gene family. Often multiple mouse input sequences resulted in the same top hit from the crustacean transcriptomes, indicating that fewer members of the gene family were present in our invertebrate sequence than the mammalian gene families. Once a complete gene family search was obtained, all putative orthologs were then blasted against the Non-Redundant Protein (NR) Sequence Database hosted at NCBI via blastx. This allowed us to look for conserved sequence across all taxa and confirm a given gene identification. Once gene families were obtained from the C. borealis transcriptome, the process was repeated with H. americanus as the subject database. We additionally used the crab sequences as queries to find the direct ortholog for a given gene in H. americanus. These sequences were confirmed in the same way via blastx against the NR database, and moved forward into sequence alignment as described below. As is to be expected from manually performed sequence-by-sequence discovery and curation such as this, at times other searches and sequence comparisons were performed on a case-by-case basis with comparator species such as Daphnia pulex, Drosophila melanogaster, or other insect species to gain insight or clarification as to the best possible identification for a given transcript.
The web-based software tool Biology Workbench (Version 3.2)  was used for sequence analysis of putative ion channels and receptors from the assembled transcriptomes. Coding regions were determined based on the longest open reading frame (ORF) from the SIXFRAME tool in Biology Workbench. ClustalW was utilized (default parameters) to perform the multiple sequence alignment (MSA) for ion channel and receptor family subtypes based on amino acid sequences from predicted coding regions. The rooted phylogenetic trees were constructed from the output of the MSA from ClustalW. The data matrix for all phylogenetic trees was deposited into TreeBASE (Study Accession URL: http://purl.org/phylo/treebase/phylows/study/TB2:S19948).
We used blastp to generate percent identities and similarities for predicted amino acid sequences of orthologs between species. We used only sequences that were full-length or those that were near full length. Sequences were assumed to be full length coding sequences if they met three criteria: the sequence began with a start codon, was approximately the same length as similar sequences in the non-redundant database based on a blastx search, and the sequence ended with a stop codon. Sequences were considered to be close to full-length if they were at least 80 % the length of similar sequences on the NR database, regardless of the presence of start and stop codons. In addition, the two sequences from crab and lobster had to be at least 80 % of the length of one another. This generated 42 pairwise comparisons for orthologous protein sequence between crabs and lobsters.
We note that all nomenclature for transcripts described in this study will conform to one nomenclature convention: transcripts from Cancer borealis, will be noted with the species prefix Cb-, and the species prefix Ha- will be used for Homarus americanus gene products. All curated gene sequences described below were submitted to GenBank and assigned individual accession numbers as noted in Tables 2, 3, 4 and 5.
Mixed nervous system transcriptome sequencing and de novo assembly
Constructing RNA-seq libraries from nervous tissues of adult crustaceans, a total of 414,978,768 and 452,237,240 raw reads were obtained from the paired - end sequencing of C. borealis and H. americanus, respectively. The average read length for both species was approximately 97 bp, as expected for 100 bp paired-end Illumina Sequencing. Following quality checks removing adaptors, contaminating sequences, and low-quality sequences, 391,060,790 (94.2 %) clean reads were found for C. borealis and 426,712,238 (94.4 %) for H. americanus. These high-quality cleaned reads were subsequently assembled de novo into contigs using two different assemblers: CLC Genomics and Seqman NGen. For C. borealis, CLC assembly resulted in 42,766 contigs with an average length of 1544 bp and an N50 length of 2178 bp, while SeqMan assembly resulted in 67,380 contigs with an average length of 1076 bp and N50 of 1239 (Table 1; Fig. 1a). For H. americanus, CLC assembly resulted in 60,273 contigs with an average length of 1657 bp and N50 length of 2357 bp, while SeqMan NGen resulted in 45,043 contigs with an average length of 1799 bp and N50 of 2258 (Table 1; Fig. 1a).
To compare the quality of our transcriptomes from multiple assemblies, a reference-based alignment was performed using BUSCO . The arthropod BUSCO reference contains 2675 orthologous genes expected within most arthropod species that were compared against the gene content of our transcriptomes. The alignment of our transcriptomes against the arthropod reference resulted in similar percentages of the reference genes found as complete (C), fragmented (F), or missing (M) across our transcriptomes (Fig. 1b). The C. borealis metrics were C: 59.0 %, F: 13.5 %, and M: 27.5 % for the CLC Genomics assembly and C: 58.4 %, F: 20.1 %, and M: 21.5 % for the SeqMan NGen assembly. The H. americanus metrics were C: 56.1 %, F: 16.4 %, and M: 27.5 % for the CLC Genomics assembly and C: 66.5 %, F: 13.4 %, and M: 20.1 % for the SeqMan NGen assembly. These results were compared against the arthropod transcriptome reference scores provided in the BUSCO supplementary materials, a recently published Homarus americanus nervous system transcriptome assembled using Trinity , and a recent transcriptome of the freshwater crayfish Astacus astacus . Our results are comparable with the Astacus transcriptome in completeness and an apparent extension of the published Homarus americanus transcriptome . One possible explanation for the missing arthropod genes from our transcriptomes can be explained by the fact that our sequences were derived solely from nervous system tissue, while the references were built from arthropod genomic sequences.
Using the NCBI BLAST+ suite to perform a blastn of 28 Cancer borealis sequences already contained within GenBank against our assembled contigs, we found an average sequence identity from the CLC assembly of 99.2 %, with the lowest identity score 96 %. SeqMan NGen assembly for C. borealis transcripts had an average sequence identity of 99.03 %, and the lowest percent identity was 95 %. For H. americanus, 75 GenBank sequences were aligned against our transcriptome, resulting in an average sequence identity to CLC assembled sequences of 99.3 % with the lowest being 94.5 %. SeqMan NGen assembly for H. americanus transcripts had an average sequence identity of 98.97 %, and the lowest percent identity was 87 %.
Based on the relative similarity in many of the metrics for these two assembly methods, the somewhat better performance of CLC contigs when compared with Sanger sequencing generated orthologs, and the fact that portions of the H. americanus transcriptome based on the CLC assembly have previously been published , we chose to perform the remaining representative analysis of these sequence data based on the CLC assembled contigs. The H. americanus Transcriptome Shotgun Assembly (TSA) project has been deposited at GenBank under Accession No. GEBG00000000 (BioProject No. PRJNA300643; BioSample No. SAMN04230440). The C. borealis Transcriptome Shotgun Assembly (TSA) project has been deposited at GenBank under the Accession No. GEFB00000000 (BioProject No. PRJNA310325; BioSample No. SAMN04450329). The versions described in this paper represent the first versions, GEBG01000000 and GEFB01000000 respectively.
Annotation and gene ontology mapping
Entrez Gene IDs were obtained for both transcriptomes using blastx against the NCBI non-redundant (NR) protein database. These annotations consisted of 9489 unique proteins among C.borealis transcripts, and 11,061 among H, americanus transcripts. Mapping these gene IDs to Gene Ontology (GO) categories yielded 9351 (22 %) of the C. borealis and 6191 (10 %) of the H. americanus contigs successfully identified (Fig. 2a). Similar percentages have been observed in other de novo transcriptome analyses [27, 28]. From the functional annotation, transcripts were classified into three broad categories: cellular compartment (CC), molecular function (MF), and biological process (BP) . Within these broad categorizations, the highest abundance GO terms of H. americanus and C. borealis were compared against each other, which included the top 9 CCs, 18 MFs, and 16 BPs for both species (Fig. 2b). The arrangement of GO terms was based on the highest abundance H. americanus terms, in descending order. The only notable exception to this order was the BP GO term “RNA-dependent DNA Replication” ontology due to its high abundance in C. borealis but relatively low abundance in H. americanus. These same GO terms were compared between the two species using the relative percentage of each GO term for its broad GO classification (CC, MF, BP) (Fig. 3). The most striking differences between the GO ontologies of each species include a much higher incidence of “protein binding” terms for C. borealis MF than that of H. americanus, a much higher incidence of “metabolic process” in H. americanus BP, and a prominent difference between the “RNA-dependent DNA Replication” term for BP. The source of these differences could be attributed to factors including, but not limited to, the variation in tissue types (abdominal nerve cord was used in H. americanus, but not C. borealis), depth of sequencing, or natural variation in transcript abundance.
Species comparisons: distribution and VennBLAST analysis
Using the Blast2GO software suite, the number of species that the C. borealis and H. americanus neural transcriptomes align with was determined from a blastx against the NCBI non-redundant database. The species distribution for both C. borealis and H. americanus gave similar top species, such as Tribolium castaneum, Daphnia pulex, and Strongylocentrotus purpuratus within the top 5 species hits (Fig. 4). The absence of termite (Zootermopsis nevadensis) from the C. borealis species distribution of blast hits is due to the fact that the Z. nevadensis protein sequences had yet to be uploaded to the NCBI non-redundant database at the time of the blast analysis of the C. borealis transcriptome, while the H. americanus analysis was performed after the Z. nevadensis reference became available.
Venn diagrams were generated (Fig. 5a) using the software VennBLAST  to compare the whole C. borealis and H. americanus transcriptomes against the Daphnia pulex (GCA_000187875.1) protein sequences from Ensembl Metazoa. The protein sequence database for D. pulex was chosen as a common subject to query the C. borealis and H. americanus transcriptomes against due to the mutual high top-hit species distribution (Fig. 4), as well as the well-annotated genome of the crustacean D. pulex . Initially, a local blastx of C. borealis or H. americanus contigs against D. pulex protein sequences resulted in 17,343 and 14,818 hits, respectively. Upon overlapping these hits with the VennBLAST Merge tool, 11,258 hits from C. borealis and H. americanus were found to have the same top hit for D. pulex. A second analysis with increased stringency was performed using the VennBLAST Filter tool to retain only high-quality matches, leaving C. borealis with 7,460 and H. americanus with 7,268 hits to D. pulex. Subsequent merger of these filtered hits resulted in 6,226 common top-hits for D. pulex, resulting in an increased percentage (from 54 % overlap to 73 %) of common top-hits.
For the remainder of our transcriptome analysis, we identified and characterized sequences for 6 different Innexin proteins (gap junctions), 34 distinct ion channel types, 17 biogenic amine receptors, 5 GABA receptors, and 28 major transmitter receptor subtypes including glutamate and acetylcholine receptors. These are described in detail below. These receptor groups consisted of 27 different ligand-gated channel subunits (ionotropic receptors) and 23 metabotropic receptor types. To better quantify the similarity across lobster and crab, we performed analyses of percent amino acid identity and similarity between orthologs of a subset of genes (see Methods). Overall, the sequence similarity is very high between these species, as one might expect for members of the same Order (Fig. 5b); across 42 genes surveyed, there was a mean ± SD of 85.27 % ± 8.46 % amino acid identity between genes. However, we also noticed a significant trend across different types of gene products: there was significantly lower amino acid identity and similarity for metabotropic receptors than for the other classes of genes (Fig. 5b). In particular, the G-protein coupled (GPCRs) receptors that we analyzed were some of the most divergent between crab and lobster. For example, Octβ-R1, 3, and 4 shared 72, 72, and 74 % amino acid identity respectively. Conversely, the most highly conserved genes were in the Shaker family of voltage-gated K+ channels. Shaker, Shal, and Shaw1were 98, 98, and 96 % identical between crabs and lobsters. From these results we would predict more conservation in channel function and physiology across species than that of the GPCRs.
For our initial analysis of these crustacean transcriptomes, we decided to focus on some of the most critical proteins involved in nervous system function. We therefore first conducted an analysis of ion channel subtypes. Putative ion channels were identified based on tblastx searches utilizing the transcriptomes as a reference database and querying with known channel protein largely consisting of sequences from mouse (Mus musculus) and Drosophila melanogaster. A 100 % overlapping set of ion channels were found to be present in both C. borealis and H. americanus (Table 2; Fig. 6). We specifically hand-curated and annotated these channel sequences, and the full list is available in Table 2, including putative current types carried by each channel. We used the multiple sequence alignment (MSA) output from CLUSTALW  to develop a fairly comprehensive ion channel tree based on amino acid sequence similarity, allowing us to cluster channels by type to effectively interrogate the nervous system channel content of these crustaceans.
Our analysis of the crab and lobster transcriptomes led us to identify and characterize 34 distinct channel subtypes representing several gene families (Table 2; Fig. 6). Each identified channel transcript was found in both crab and lobster transcriptomes. Our analysis confirmed the presence of 3 major voltage-dependent calcium channel subtypes  corresponding one each to the L-type (CaV1), P/Q-type (CaV2), and T-type (CaV3) families of calcium channels. In addition, we identified a single member of the NaV-type voltage-gated Na+ channel representing the para type channel identified in other species. Finally, in both species we identified a non-selective sodium leak channel (NALCN) thought to underlie TTX-resistant Na+ conductance important in baseline neuronal excitability . One other major family of non-selective cation channels we identified was the cyclic-nucleotide-gated channels of the HCN/CNG type. Both species contained one member of the hyperpolarization-activated cyclic nucleotide-gated (HCN) channel family, the channels that give rise to IH type currents in crustacean neurons . In addition, we identified 3 α -subunits and 1 β-subunit of the cyclic nucleotide gated ion channel types (CNG), which are activated by the binding of cAMP and cGMP to carry a non-selective cation current .
The pore-forming α-subunits of K+ channels can be sub-divided into voltage-dependent subunits (Kv), inward rectifiers (Kir), two-pore subunits (K2P), and those activated by intracellular calcium (KCa) or sodium (KNa) ions. We identified a diverse array of voltage-dependent potassium channel subtypes in the transcriptomes of both the crab and lobster (Table 2; Fig. 5). The best already characterized of these channels in crustaceans are the Shaker family of channels, having been identified in crab [35, 36] and spiny lobsters (Panulirus interruptus, [37, 38]), with the latter having an extensive characterization via expression in oocyte systems [39, 40]. Previously, 4 members of this family were already known from Cancer borealis: shaker, shal, shab, and shaw. We found orthologs to each of these in H. americanus as well. We further discovered that in both species there actually were 2 distinct shaw-related channel transcripts (Fig. 6). The newly re-named “Shaw1” transcript from this analysis is a perfect match with the previously identified shaw transcript from Cancer borealis (Accession #EF089569), while the newly identified shaw-like transcript is presented as Shaw2. We also identified two members of the KCNQ family of K+ channels. KCNQ genes encode a family of six transmembrane domain K+ channel alpha-subunits that have a wide range of physiological roles, including likely underlying the slow voltage-gated M-type currents . Rounding out the voltage-dependent K+ channel subtypes are 3 members of the ether-a-go-go/KCNH family. In addition to voltage-dependent K+ channels, we also identified one Kir channel (IRK), two members of the K2P family (KCNK), one sodium-activated K+ channel (KCNT), and two calcium-activated K+ channel types. These KCa channels had previously been identified in C. borealis  and correspond to one BK- and one SK- type channel.
Transient receptor potential (TRP) channels have been implicated as a primary channel for generation of sensations including temperature, taste, pain, pressure, and vision. In our analysis, we found various TRP subfamilies within both C. borealis and H. americanus (Table 2; Fig. 6). These subfamilies included TRPV (vanilloid), mediating odor and pain sensations; TRPA (ankyrin), associated with mechanical stress receptors; TRPM (melastatin), associated with magnesium reabsorption ; and TRP pyrexia, a thermal sensing receptor . In crustaceans, TRP channels have been primarily studied for their role in olfactory reception  and stretch reception . Not all orthologs of TRP channels were identified in both species. We did not identify in this data set orthologs of TRP-V6 and pyrexia from C. borealis and an ortholog of TRP-M1 was not identified in H. americanus. It is most likely that these “missing” orthologs are due to limitations in the sequence depth, although we cannot rule out the possibility that these two species have distinct complements of TRP channel genes. No sequences were found in either species that represent the TRPC, TRPP, TRPL, or TRPN subfamilies. These results are consistent with found in a previously published transcriptome of H. americanus that identified 2 TRPA, one pyrexia, and two TRPM type channels . We extend these results to include identification of the TRPV family of channels in both H. americanus and C. borealis.
Biogenic amine receptors
Biogenic amine neuromodulators were some of the first modulatory compounds to be thoroughly studied in the crustacean nervous system [46–49], and specifically in the stomatogastric nervous system [47, 50–52] where some of the most comprehensive understanding of the multiple targets and modulatory impacts of these compounds on neural circuits has been described [53–59]. Therefore, we decided it would be valuable to provide a thorough characterization of these receptor subtypes as well to complement the extensive and elegant physiology work that has been going on for decades.
Dopamine has been perhaps the most extensively characterized biogenic amine from a functional and biochemical perspective in crustaceans. Previous work [60, 61] identified 3 subtypes of dopamine receptors in the nervous system of the spiny lobster, Panulirus interruptus: D1αPan (Type 1A DAr), D1βPan (Type 1B DAr), and D2αPan (Type 2 DAr). We found clear orthologs to all three of these receptor subtypes in both C. borealis and H. americanus (Table 3; Fig. 7), and our transcriptome search protocol did not come up with any other putative DAr subunit transcripts. Therefore, it is likely that these three receptor subtypes represent the complete complement of dopamine receptors in these decapod crustaceans. In deference to the extensive characterization of these receptors in the closely related spiny lobster, we conform the naming of these channels to match with the Panulirus nomenclature: for example, Cb-D1αR, Cb-D1βR, and Cb-D2αR (Fig. 7).
Serotonin receptors are less described in crustaceans than the dopamine receptors. Previous reports describe two distinct subtypes of serotonin receptors in the P. interruptus  as well as the crayfish Procambarus clarkii : one type-1 and one type-2 serotonin receptor. Our analysis of the transcriptomes of C. borealis and H. americanus found clear orthologs to both of these receptor subtypes, and based on homology with mouse and Drosophila sequences we identified two novel putative serotonin receptor subtypes as well. Figure 6 uses the existing P. interruptus and P. clarkii sequences to root the new sequences in a tree representing these crustacean serotonin receptors. Our analysis suggests that the previous type-1 5HTr subtypes identified are most similar to mammalian 5HT1A (HTR1A) receptors, while the crustacean type-2 5HTr subunit is most similar to 5HT2B (HTR2B) receptors. We also identified a putative type-1B receptor (HTR1B) as well as a putative type-7 receptor (HTR7). These identities are assigned based on mouse query sequences used in our tblastx protocols that generated the strongest hits (i.e. lowest e-values) when queried against the crustacean transcriptomes. We follow the mammalian classification and nomenclature guidelines for these 5HT receptors in assigning gene names (Table 3; Fig. 7), as these are well defined and organized relative to the invertebrate nomenclature: e.g. Cb-HTR1A, Cb-HTR1B, Cb-HTR2B, and Cb-HTR7.
Octopamine receptors are virtually undescribed in crustaceans, with the sole decapod receptor described as a tyramine/octopamine receptor from the freshwater prawn, Macrobrachium rosenbergii . Thorough work with crustacean octopamine receptors is found in the barnacle, Balanus improvisus, where one alpha- and four beta-like receptor subtypes have been very nicely characterized . Our analysis identified the same distribution of receptor types in C. borealis and H. americanus as was described in the barnacle – one alpha- and four beta-like subunits (Table 3; Fig. 7). However, there were no particularly conserved motifs that resulted in a clustering of decapod and barnacle receptor subtypes to converge on a common nomenclature for these receptors; the four beta-like receptors in barnacle most closely related one another rather than subtypes across species. As a result, we have simply named these β-like octopamine receptors subtypes with ascending numbers and in the style of the descriptions given to those identified in B. improvisus (Bi): e.g. Cb-Oct-αR, Cb-Octβ-R1, Cb-Octβ-R2, Cb-Octβ-R3, and Cb-Octβ-R4. However, we do not mean to imply direct orthology between BiOctβ-R1 and Cb-Octβ-R1, for example. Finally, we also identified a single putative tyramine receptor in both of our species (Cb-Tyr-R and Ha-Tyr-R), and these sequences are most similar to the tyramine/octopamine receptor reported from Macrobrachium rosenbergii, which is most effectively activated by tyramine [64, 66].
Histamine is also a major neurotransmitter in invertebrates, but sequence information for histamine receptors in crustaceans has yet to be described. In vertebrates, these receptors are part of a protein family of G-protein-coupled metabotropic receptors. In invertebrates, histamine acts exclusively through ionotropic histamine receptors, as no metabotropic histamine receptors appear to be present in invertebrates . We discovered putative orthologs for 3 distinct ionotropic histamine receptors in both C. borealis as well as H. americanus (HisR1-3) and a fourth sequence was found in H. americanus that did not have an obvious match from C. borealis (Fig. 7). In sequence comparisons with Drosophila melanogaster, the closest orthologs are described histamine-gated chloride channels. This is consistent with the known physiological characterization of crustacean histamine gating of chloride channels [68, 69]. Because naming conventions are lacking for these receptors in crustaceans, we have simply assigned these transcripts ascending numbers (1–4) to distinguish between different putative histamine receptor subtypes (Fig. 7; Table 3).
Metabotropic glutamate receptors
Metabotropic glutamate receptors (mGluRs) are seven-transmembrane domain proteins coupled to G-proteins capable of controlling many cellular processes through signaling cascades. mGluRs can be classified into three primary classes : Class I consists of mGluR1 and mGluR5 and are associated with phospholipase-C and utilize intracellular calcium signaling cascades; Class II consists of mGluR2 and mGluR3; and Class III consists of mGluR4, mGluR6, mGluR7, and mGluR8 and are negatively coupled with adenylyl cyclase activity. These receptors exist as either homo- or heterodimers on the cell surface, and it is the associated G-protein alpha-subunit that determines which class the mGluR falls under (e.g. Class I associates with Gq and G11). In crustaceans, mGluRs have a relatively short history of study, mostly owing to the fact that the metabotropic form of glutamate receptors had not been characterized until the late 1980s . In crustacean preparations, it has been found that mGluRs play a role in rhythm generation in the stomatogastric ganglion [72, 73]. In our analysis of these transcriptomes, we found six mGluR sequences for each species, which covered all three primary classes of mGluRs (Table 4; Fig. 8). We did not find mGluR6 and mGluR8 orthologs in either species. It should also be noted that Ha-mGluR2 and Ha-mGluR4 aligned more closely to one another than to their C. borealis counterparts. This discrepancy could be due to the relatively short partial sequence found for Cb-mGluR2; that is, only the first 200 amino acids were found for Cb-mGluR2, while the Ha-mGluR2 sequence found is 1027 amino acids long.
Ionotropic glutamate receptors
The most common excitatory neurotransmitter found in crustaceans at the neuromuscular synapse is glutamate . The three primary ionotropic glutamate receptors are NMDA, AMPA, and kainate receptors, named respectively after the agonists N-methyl-D-aspartate, α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid, and kainic acid that activate them. One of the initial characterizations of NMDA receptors from crustaceans was performed in the crayfish optic lobe . Since then, NMDA receptors have been studied in crustaceans for their role in memory , axon-to-glial signaling , and central pattern generation . We were able to identify 4 separate NMDA receptors from both C. borealis and H. americanus, which fell into two primary categories: 1-like and 2-like (Table 4; Fig. 8). We further identified these receptors into 1A, 1B, 2A, and 2B based on their pairing, but this is not meant to imply orthology with other naming schemes in other species. We also found one NMDA receptor in both C. borealis and H. americanus that we did not find a highly similar pair for, which we have included as Cb-NMDA-2-like and Ha-NMDA-2-like.
AMPA receptors have been virtually undescribed in crustacean preparations, which coincides with our results. In our analysis, we did not find any receptors that most closely resembled AMPA receptors. Known AMPA receptors blasted against our transcriptomes aligned best against the putative kainite-type receptors, which is unsurprising considering both AMPA and kainate are considered non-NMDA receptors. Kainate receptors have been implicated as modulators of synaptic transmission and excitability . In crustacean systems, kainic acid historically has been shown to elicit depolarizations at the crab neuromuscular junction , as well as the crayfish neuromuscular junction . Five pairs of kainite-like receptors were found between the C. borealis and H. americanus transcriptomes (Table 4; Fig. 8), falling into two discrete categories: kainate type-1-like and 2-like. Sequences were further subdivided into A, B, and C on the simple basis of pairing between the two species and not based on any specific orthology to other species.
Beyond the excitatory glutamate-gated cation channels (NMDA- and kainate-like), an inhibitory glutamate-gated chloride channel (GluCl) was also found for both C. borealis and H. americanus (Table 4). Soon after their discovery as extrajunctional receptors in locust muscles , GluCls were described as postsynaptic receptors in the crustacean stomatogastric ganglion [50, 83]. In our analysis, we found a single GluCl transcript for each species, named Cb-GluCl and Ha-GluCl. The finding of a single channel is consistent with some other invertebrates, with a single GluCl also found in most insects . Because of their distinct characteristics, the GluCl channels we identified were not placed on any of the trees shown in the Figures.
The neurotransmitter γ-amino butyric acid (GABA) has been studied in crustacean species for decades for its role in synaptic transmission and neural inhibition [85–88]. Interestingly, in invertebrates several different neural responses to GABA have been found that have distinct profiles from that of vertebrate GABA receptors [83, 89–91]. GABA receptors are classified into two major groups: GABAA type receptors, comprising receptor complexes that are part of a ligand-gated ion channels, or GABAB type receptors, G-protein-coupled receptors that act via metabotropic signaling systems. GABAA receptors are pentameric transmembrane receptors responsible for fast, usually inhibitory synaptic currents, and heteromultimers of the individual subunit types can form distinct channel properties in invertebrates . We identified orthologs of three GABAA type receptor subunits from both H. americanus and C. borealis (Fig. 9; Table 3), including orthologs of Drosophila LCCH3-, RDL-, and GRD-like receptor subunits – and we have preserved naming conventions for these subtypes. Two GABA receptor subunits previously have been cloned from H. americanus , and very recently a GABAA type receptor was identified in the crayfish, Procambarus clarkii . Sequence comparison reveals these previously described sequences to be orthologous to the RDL-like receptor from our data. The metabotropic GABAB type receptors are GPCRs responsible largely for slower inhibitory synaptic effects, and functional GABAB receptors are heterodimers formed by GABAB1 and GABAB2 subunits . We identified orthologs of both GABAB1 and GABAB2 subunits in both C. borealis and H. americanus (Table 3; Fig. 9).
Acetylcholine receptors are classified into two family subtypes: nicotinic receptors (nAChRs), which are ligand-gated ion channels that are activated by nicotine; and muscarinic receptors (mAChRs), which are metabotropic GPCRs that respond to the agonist muscarine. Muscarinic acetylcholine receptors are further classified into subtypes based on the specific G-protein associated with the receptor . In our analysis, we discovered two discrete subtypes of mAChRs (Table 4; Fig. 9) in both C. borealis and H. americanus, which is consistent with other arthropods . The A- and B-type mAChR are defined based on their differential sensitivity to muscarine (A-type is 1000x more sensitive than B-type), as well as the antagonist binding properties (atropine, scopolamine, and QNB block A-type but not B-type) that each receptor exhibits .
Nicotinic acetylcholine receptors are common throughout the invertebrate central nervous system, mediating largely fast excitatory neurotransmission [98–101]. For nAChRs in nervous systems of C. borealis and H. americanus, each species was found to have 1 β-subtype and 8 α-subtypes (Fig. 9). A thorough characterization of invertebrate nAChRs has been performed for the snail, Lymnaea stagnalis . Our analysis revealed that few direct orthologous sequences occurred from snail (Mollusca) to crustaceans that would allow us to adopt the nomenclature put forward in Lymnaea. Therefore, the crustacean receptor subtypes were named based on their most similar mammalian counterpart, with the exception of the α16 subunit found, named so due to its similarity to acr-16 found in C. elegans . This subunit is most comparable to α7 in humans, but we found other sequences nAChRs more similar to human α7 than that of the putative α16.
Gap junction proteins (Innexins)
The proteins responsible for gap junctions in invertebrates are the Innexins . A family of Innexins have previously been described for both C. borealis and H. americanus . We include them here as characterized through the transcriptome analysis for completeness. Three full-length sequences were named Innexins 1–3 (Fig. 10; Table 5) based on significant sequence similarity to coding sequences of Innexins from multiple other organisms, including other decapod crustaceans . Innexins 4–6 (Fig. 10; Table 5) were subsequently identified from our transcriptome analysis of C. borealis nervous system. We identified clear orthologs from lobster and crab for Innexins 1–4 and 6; but two Innexins showed enough dissimilarity to be classified separately, and these are classified as Innexin 5 in C. borealis and Innexin 7 in H. americanus (Fig. 10). All of these identified Innexin sequences have the signature motif YYQWV in the second TM domain as well as a series of other conserved amino acid residues considered hallmarks of Innexins .
The era of modern genomics and high-throughput sequencing has revolutionized the study of neuroscience, and provided an opportunity for classic physiology systems in the study of neural circuit properties to experience a renewed level of impact. In particular, invertebrate model systems that historically have been invaluable to our understanding of basic circuit properties, dynamics, systems neuroscience, and neuromodulation now present themselves as novel contributors to molecular neuroscience. In particular, classic preparations such as the Tritonia swim system , Aplysia feeding circuits as well as the classic gill and siphon withdrawal reflex [108, 109], and crustacean stomatogastric systems  revolutionized our understanding of neural circuitry. Each of these systems is the renewed focus of genomic and transcriptomic approaches [4, 5, 110, 111] – including this study – that promise to merge the unparalleled experimental accessibility on the neurophysiology end of the spectrum with new molecular tools to understand and manipulate these circuits. Decapod crustacean systems have also been foundational in the understanding of modulation of behavioral states. The earliest work implicating serotonin broadly in aggression can be traced to seminal work in lobster behavioral studies [48, 112]. GABA was first identified as an inhibitory transmitter in these same decapod lobster species [9, 113]. Finally, the crayfish escape behavioral response has been a paradigm for the true integration of neuroethological work across single neurons, neural circuitry, behavior, modulation, and social status . Therefore, the accessibility of a molecular perspective and toolset will allow researchers to revisit these seminal works with greater potential to understand integrated nervous system function.
It can be challenging to stay up to date with all of the sequence data being published. To the best of our knowledge, there are a relatively small number of published transcriptome projects with decapod crustaceans as model systems. Sequence discovery projects in decapod crustaceans began with expressed sequence tag (EST) analyses over a decade ago [115, 116], and these have been used to examine olfactory receptor expression in the lobster system . Since then, a mixed tissue transcriptome sequencing projects have been performed from the spider crab, Hyas araneus [118, 119] and spiny lobster Sagmariasus verreauxi . Some of the most proactive work in crustacean transcriptomics lies in the area of neuropeptidomics for the discovery of neuroactive peptides and their receptors [121–124], including a study that makes use of the H. americanus sequence data contained reported here . However, none of these studies provides analysis of ion channel or receptor sequence as contained within our study.
The most directly comparable transcriptome study was carried out in H. americanus, and used RNA-seq to examine differences in gene expression across tissues, including nervous system, as well as responses to hormonal treatment and temperature perturbations . While this study uses GO and other annotation approaches to identify channel and receptor groups that are differentially expressed, identification of these gene products is at the level of gene family. For example, McGrath et al.  identify glutamate receptors as a broad group that are differentially expressed in tissues, but make no further distinction as to receptor subtype –reporting that they have 126 transcripts that map to a global glutamate receptor subtype group. This high number is likely due to identification of multiple short contigs of the same receptor sequence, as opposed to distinct receptor sequences, and may be reflected by the lower N50 value of this transcriptome (1,289). In contrast, the N50 of our H. americanus transcriptome is 2,357, leading to longer contigs and often the identification of full-length coding sequence that greatly aided our ability to characterize and identify distinct channel and receptor subtypes. Nevertheless, to the extent these studies can be compared, we saw results largely consistent with the overall identifications provided . For example, both studies identified 5 octopamine receptor subunits, including at least 3 beta-type receptors. Further, all reported differentially expressed voltage-dependent ion channel and transmitter receptor subtypes  were present in our data set. While there is no direct overlap in the goals of these studies, the characterization provided in our study will provide a reference by which subsequent expression analyses can be combined with thorough annotation to provide more insight into gene-specific changes related to neural function.
In this study we sequenced the nervous system transcriptomes for two highly utilized species in invertebrate neuroscience research: the Jonah crab (Cancer borealis) and the American lobster (Homarus americanus). Our sequencing, assembly, and annotation efforts have yielded an extensive set of sequence information from which we can begin to mine gene products critical to fundamental nervous system output: channels and receptors. This study represents the first attempt to characterize to this extent these critical building blocks of circuit function from these model systems. In doing so we have identified for the first time in these species previously undescribed channel and receptor families, as well as added to the incomplete characterization of amine receptors known to modulate both circuit function and behavior in these animals. This sequence information opens up these target proteins for use in gene manipulation techniques such as overexpression  or RNA-interference mediate knockdown  to deeply interrogate circuit function. Finally, the stomatogastric system has been used extensively in computational studies that have revolutionized our understanding of circuit fundamentals, dynamics, and the role of variability in neuronal parameters in circuit function [127–130]. These models have relied on biological data for identification of likely membrane conductances present in the networks. Molecular screening and quantitative assays of channel expression can effectively be used in concert with computational modeling [131, 132] to generate better and more biologically realistic models with which to uncover fundamental aspects of neural circuit dynamics.
Benchmarking universal single-copy orthologs
Expressed sequence tag
Glutamate-gated chloride channel
G-protein coupled receptor
Muscarinic acetylcholine receptor
Metabotropic glutamate receptor
Multiple sequence alignment
Nicotinic acetylcholine receptor
Open reading frame
Sequence read archive
Stomatogastric nervous system
Transient receptor potential
Transcriptome shotgun assembly
Colbourne JK, Pfrender ME, Gilbert D, Thomas WK, Tucker A, Oakley TH, et al. The ecoresponsive genome of Daphnia pulex. Science. 2011;331:555–61.
Zeng V, Villanueva KE, Ewen-Campen BS, Alwes F, Browne WE, Extavour CG. De novo assembly and characterization of a maternal and developmental transcriptome for the emerging model crustacean Parhyale hawaiensis. BMC Genomics. 2011;12:581.
Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, et al. The genome sequence of Drosophila melanogaster. Science. 2000;287:2185–95.
Christie AE, Chi M, Lameyer TJ, Pascual MG, Shea DN, Stanhope ME, et al. Neuropeptidergic signaling in the american lobster homarus americanus: New insights from high-throughput nucleotide sequencing. PLoS One. 2015;10:e0145964.
McGrath LL, Vollmer SV, Kaluziak ST, Ayers J. De novo transcriptome assembly for the lobster Homarus americanus and characterization of differential gene expression across nervous system tissues. BMC Genomics. 2016;17:63.
Wiersma CA, Ikeda K. Interneurons commanding swimmeret movements in the crayfish, Procambarus clarkii (girard). Comp Biochem Physiol. 1964;12:509–25.
Furshpan EJ, Potter DD. Transmission at the giant motor synapses of the crayfish. J Physiol. 1959;145:289–325.
Dudel J, Kuffler S. The quantal nature of transmission and spontaneous miniature potentials at the crayfish neuromuscular junction. J Physiol. 1961;155:514–29.
Otsuka M, Iversen LL, Hall ZW, Kravitz EA. Release of gamma-aminobutyric acid from inhibitory nerves of lobster. Proc Natl Acad Sci U S A. 1966;56:1110–5.
Otsuka M, Kravitz EA, Potter DD. Physiological and chemical architecture of a lobster ganglion with particular reference to gamma-aminobutyrate and glutamate. J Neurophysiol. 1967;30:725–52.
Connor JA, Walter D, McKown R. Neural repetitive firing: modifications of the Hodgkin-Huxley axon suggested by experimental results from crustacean axons. Biophys J. 1977;18:81–102.
Connor JA. Neural repetitive firing: a comparative study of membrane properties of crustacean walking leg axons. J Neurophysiol. 1975;38:922–32.
Remler M, Selverston A, Kennedy D. Lateral giant fibers of cray fish: location of somata by dye injection. Science. 1968;162:281–3.
Stretton AO, Kravitz EA. Neuronal geometry: determination with a technique of intracellular dye injection. Science. 1968;162:132–4.
Hartline D. Impulse identification and axon mapping of the nine neurons in the cardiac ganglion of the lobster Homarus americanus. J Exp Biol. 1967;47:327–41.
Maynard DM. Simpler networks. Ann N Y Acad Sci. 1972;193:59–72.
Mulloney B, Selverston AI. Organization of the stomatogastric ganglion of the spiny lobster - III. Coordination of the two subsets of the gastric system. J Comp Physiol. 1974;91:53–78.
Heitler WJ. Coupled motoneurones are part of the crayfish swimmeret central oscillator. Nature. 1978;275:231–4.
Marder E, Bucher D. Understanding circuit dynamics using the stomatogastric nervous system of lobsters and crabs. Annu Rev Physiol. 2007;69:291–316.
Marder E. Neuromodulation of neuronal circuits: Back to the future. Neuron. 2012;76:1–11.
Diehl F, White RS, Stein W, Nusbaum MP. Motor circuit-specific burst patterns drive different muscle and behavior patterns. J Neurosci. 2013;33:12013–29.
Williams AH, Kwiatkowski MA, Mortimer AL, Marder E, Lou ZM, Dickinson PS. Animal-to-animal variability in the phasing of the crustacean cardiac motor pattern: an experimental and computational analysis. J Neurophysiol. 2013;109:2451–65.
Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31:3210–2.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.
Zahavi T, Stelzer G, Strauss L, Salmon AY, Salmon-Divon M. VennBLAST—whole transcriptome comparison and visualization tool. Genomics. 2015;105:131–6.
Subramaniam S. The Biology Workbench--a seamless database and analysis environment for the biologist. Proteins. 1998;32:1–2.
Elmer KR, Fan S, Gunter HM, Jones JC, Boekhoff S, Kuraku S, et al. Rapid evolution and selection inferred from the transcriptomes of sympatric crater lake cichlid fishes. Mol Ecol. 2010;19:197–211.
Gu XC, Zhang YN, Kang K, Dong SL, Zhang LW. Antennal transcriptome analysis of odorant reception genes in the red turpentine beetle (RTB), Dendroctonus valens. PLoS One. 2015;10:e0125159.
Theissinger K, Falckenhayn C, Blande D, Toljamo A, Gutekunst J, Makkonen J, et al. De Novo assembly and annotation of the freshwater crayfish Astacus astacus transcriptome. Mar Genomics. 2016;28:7–10.
Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Ransdell JL, Temporal S, West NL, Leyrer ML, Schulz DJ. Characterization of inward currents and channels underlying burst activity in motoneurons of crab cardiac ganglion. J Neurophysiol. 2013;110:42–54.
Lu B, Su Y, Das S, Liu J, Xia J, Ren D. The neuronal channel NALCN contributes resting sodium permeability and is required for normal respiratory rhythm. Cell. 2007;129:371–83.
Ouyang Q, Goeritz M, Harris-Warrick RM. Panulirus interruptus Ih-channel gene PIIH: modification of channel properties by alternative splicing and role in rhythmic activity. J Neurophysiol. 2007;97:3880–92.
Podda MV, Grassi C. New perspectives in cyclic nucleotide-mediated functions in the CNS: the emerging role of cyclic nucleotide-gated (CNG) channels. Pflügers Arch. Eur J Physiol. 2014;466:1241–57.
Schulz DJ, Goaillard J-M, Marder E. Variable channel expression in identified single and electrically coupled neurons in different animals. Nat Neurosci. 2006;9:356–62.
Schulz DJ, Goaillard J-M, Marder EE. Quantitative expression profiling of identified neurons reveals cell-specific constraints on highly variable levels of gene expression. Proc Natl Acad Sci U S A. 2007;104:13187–91.
Baro DJ, Coniglio LM, Cole CL, Rodriguez HE, Lubell JK, Kim MT, et al. Lobster shal: comparison with Drosophila shal and native potassium currents in identified neurons. J Neurosci. 1996;16:1689–701.
Baro DJ, Cole CL, Harris-Warrick RM. The lobster shaw gene: Cloning, sequence analysis and comparison to fly shaw. Gene. 1996;170:267–70.
Kim M, Baro DJ, Lanning CC, Doshi M, Moskowitz HS, Farnham J, et al. Expression of Panulirus shaker potassium channel splice variants. Receptors Channels. 1998;5:291–304.
Baro DJ, Quiones L, Lanning CC, Harris-Warrick RM, Ruiz M. Alternate splicing of the shal gene and the origin of IA diversity among neurons in a dynamic motor network. Neuroscience. 2001;106:419–32.
Robbins J. KCNQ potassium channels: physiology, pathophysiology, and pharmacology. Pharmacol Ther. 2001;90:1–19.
Wissenbach U, Niemeyer BA, Flockerzi V. TRP channels as potential drug targets. Biol Cell. 2004;96:47–54.
Lee Y, Lee Y, Lee J, Bang S, Hyun S, Kang J, et al. Pyrexia is a new thermal transient receptor potential channel endowing tolerance to high temperatures in Drosophila melanogaster. Nat Genet. 2005;37:305–10.
Bobkov YV, Ache BW. Pharmacological properties and functional role of a TRP-related ion channel in lobster olfactory receptor neurons. J Neurophysiol. 2005;93:1372–80.
Rydqvist B, Purali N. Transducer properties of the rapidly adapting stretch receptor neurone in the crayfish (Pacifastacus leniusculus). J Physiol. 1993;469:193–211.
Dudel J. Facilitatory effects of 5-hydroxy-tryptamine on the crayfish neuromuscular junction. Naunyn Schmiedebergs Arch Exp Pathol Pharmakol. 1965;249:515–28.
Barker DL, Kushner PD, Hooper NK. Synthesis of dopamine and octopamine in the crustacean stomatogastric nervous system. Brain Res. 1979;161:99–113.
Livingstone MS, Harris-Warrick RM, Kravitz EA. Serotonin and octopamine produce opposite postures in lobsters. Science. 1980;208:76–9.
Glanzman DL, Krasne FB. Serotonin and octopamine have opposite modulatory effects on the crayfish’s lateral giant escape reaction. J Neurosci. 1983;3:2263–9.
Marder E, Eisen JS. Transmitter identification of pyloric neurons: electrically coupled neurons use different transmitters. J Neurophysiol. 1984;51:1345–61.
Flamm RE, Harris-Warrick RM. Aminergic modulation in lobster stomatogastric ganglion. I. Effects on motor pattern and activity of neurons within the pyloric circuit. J Neurophysiol. 1986;55:847–65.
Flamm RE, Harris-Warrick RM. Aminergic modulation in lobster stomatogastric ganglion. II. Target neurons of dopamine, octopamine, and serotonin within the pyloric circuit. J Neurophysiol. 1986;55:866–81.
Beltz B, Eisen JS, Flamm R, Harris-Warrick RM, Hooper SL, Marder E. Serotonergic innervation and modulation of the stomatogastric ganglion of three decapod crustaceans (Panulirus interruptus, Homarus americanus and Cancer irroratus). J Exp Biol. 1984;109:35–54.
Cleland TA, Selverston AI. Dopaminergic modulation of inhibitory glutamate receptors in the lobster stomatogastric ganglion. J Neurophysiol. 1997;78:3450–2.
Ayali A, Johnson BR, Harris-Warrick RM. Dopamine modulates graded and spike-evoked synaptic inhibition independently at single synapses in pyloric network of lobster. J Neurophysiol. 1998;79:2063–9.
Johnson BR, Kloppenburg P, Harris-Warrick RM. Dopamine modulation of calcium currents in pyloric neurons of the lobster stomatogastric ganglion. J Neurophysiol. 2003;90:631–43.
Johnson BR, Schneider LR, Nadim F, Harris-Warrick RM. Dopamine modulation of phasing of activity in a rhythmic motor network: contribution of synaptic and intrinsic modulatory actions. J Neurophysiol. 2005;94:3101–11.
Harris-Warrick RM, Johnson BR. Checks and balances in neuromodulation. Front Behav Neurosci. 2010;4:1–9.
Kiehn O, Harris-Warrick RM. 5-HT modulation of hyperpolarization-activated inward current and calcium-dependent outward current in a crustacean motor neuron. J Neurophysiol. 1992;68:496–508.
Clark MC, Baro DJ. Molecular cloning and characterization of crustacean type-one dopamine receptors: D1alphaPan and D1betaPan. Comp Biochem Physiol B Biochem Mol Biol. 2006;143:294–301.
Clark MC, Khan R, Baro DJ. Crustacean dopamine receptors: Localization and G protein coupling in the stomatogastric ganglion. J Neurochem. 2008;104:1006–19.
Clark MC, Dever TE, Dever JJ, Xu P, Rehder V, Sosa MA, et al. Arthropod 5-HT2 receptors: a neurohormonal receptor in decapod crustaceans that displays agonist independent activity resulting from an evolutionary alteration to the DRY motif. J. Neurosci. 2004;24:3421–35.
Sosa MA, Spitzer N, Edwards DH, Baro DJ. A crustacean serotonin receptor: Cloning and distribution in the thoracic ganglia of crayfish and freshwater prawn. J Comp Neurol. 2004;473:526–37.
Reyes-Colón D, Vázquez-Acevedo N, Rivera NM, Jezzini SH, Rosenthal J, Ruiz-Rodríguez EA, et al. Cloning and distribution of a putative octopamine/tyramine receptor in the central nervous system of the freshwater prawn Macrobrachium rosenbergii. Brain Res. 2010;1348:42–54.
Lind U, Alm Rosenblad M, Hasselberg Frank L, Falkbring S, Brive L, Laurila JM, et al. Octopamine receptors from the barnacle Balanus improvisus are activated by the alpha2-adrenoceptor agonist medetomidine. Mol Pharmacol. 2010;78:237–48.
Jezzini SH, Reyes-Colón D, Sosa MA. Characterization of a prawn OA/TA receptor in Xenopus oocytes suggests functional selectivity between octopamine and tyramine. PLoS One. 2014;9:e111314.
Roeder T. Metabotropic histamine receptors—nothing for invertebrates? Eur J Pharmacol. 2003;466:85–90.
McClintock TS, Ache BW. Histamine directly gates a chloride channel in lobster olfactory receptor neurons. Proc Natl Acad Sci U S A. 1989;86:8137–41.
Claiborne BJ, Selverston AI. Histamine as a neurotransmitter in the stomatogastric nervous system of the spiny lobster. J Neurosci. 1984;4:708–21.
Ferraguti F, Shigemoto R. Metabotropic glutamate receptors. Cell Tissue Res. 2006;326:483–504.
Hirono C, Ito I, Yamagishi S, Sugiyama H. Characterization of glutamate receptors induced in Xenopus oocytes after injection of rat brain mRNA. Neurosci Res. 1988;6:106–14.
Krenz WD, Nguyen D, Pérez-Acevedo NL, Selverston AI. Group I, II, and III mGluR compounds affect rhythm generation in the gastric circuit of the crustacean stomatogastric ganglion. J Neurophysiol. 2000;83:1188–201.
Pérez-Acevedo NL, Krenz WD. Metabotropic glutamate receptor agonists modify the pyloric output of the crustacean stomatogastric ganglion. Brain Res. 2005;1062:1–8.
Takeuchi A, Takeuchi N. The effect on crayfish muscle of iontrophoretically applied glutamate. J Physiol. 1964;170:296–317.
Pfeiffer-Linn C, Glantz RM. An arthropod NMDA receptor. Synapse. 1991;9:35–42.
Troncoso J, Maldonado H. Two related forms of memory in the crab Chasmagnathus are differentially affected by NMDA receptor antagonists. Pharmacol Biochem Behav. 2002;72:251–65.
Gafurov BS, Urazaev AK, Grossfeld RM, Lieberman EM. Mechanism of NMDA receptor contribution to axon-to-glia signaling in the crayfish medial giant nerve fiber. Glia. 2002;38:80–6.
Cleland TA, Selverston AI. Glutamate-gated inhibitory currents of central pattern generator neurons in the lobster stomatogastric ganglion. J Neurosci. 1995;15:6631–9.
Contractor A, Mulle C, Swanson GT. Kainate receptors coming of age: milestones of two decades of research. Trends Neurosci. 2011;34:154–63.
King AE, Wheal HV. The excitatory actions of kainic acid and some derivatives at the crab neuromuscular junction. Eur J Pharmacol. 1984;102:129–34.
Shinozaki H, Ishida M. Electrophysiological studies of kainate, quisqualate, and ibotenate action on the crayfish neuromuscular junction. Adv Biochem Psychopharmacol. 1981;27:327–36.
Cull-Candy SG. Two types of extrajunctional L-glutamate receptors in locust muscle fibres. J Physiol. 1976;255:449–64.
Marder E, Paupardin-Tritsch D. The pharmacological properties of some crustacean neuronal acetylcholine, gamma-aminobutyric acid, and L-glutamate responses. J Physiol. 1978;280:213–36.
Semenov EP, Pak WL. Diversification of Drosophila chloride channel gene by multiple posttranscriptional mRNA modifications. J Neurochem. 1999;72:66–72.
Kravitz E, Potter DD, Van Gelder N. Gamma-aminobutyric acid and other blocking substances extracted from crab muscle. Nature. 1962;194:382–3.
Iversen LL, Kravitz EA, Otsuka M. Release of gamma-aminobutyric acid (GABA) from lobster inhibitory neurones. J Physiol. 1967;188:21P–2.
Dudel J, Gryder R, Kaji A, Kuffler SW, Potter DD. Gamma-aminobutyric acid and other blocking compounds in crustacea. I. Central nervous system. J Neurophysiol. 1963;26:721–8.
Florey E, Chapman DD. The non-identity of the transmitter substance of crustacean inhibitory neurons and gamma-aminobutyric acid. Comp Biochem Physiol. 1961;3:92–8.
Fuchs PA, Getting PA. Ionic basis of presynaptic inhibitory potentials at crayfish claw opener. J Neurophysiol. 1980;43:1547–57.
Swensen AM, Golowasch J, Christie AE, Coleman MJ, Nusbaum MP, Marder E. GABA and responses to GABA in the stomatogastric ganglion of the crab Cancer borealis. J Exp Biol. 2000;203:2075–92.
Beg AA, Jorgensen EM. EXP-1 is an excitatory GABA-gated cation channel. Nat Neurosci. 2003;6:1145–52.
Gisselmann G, Plonka J, Pusch H, Hatt H. Drosophila melanogaster GRD and LCCH3 subunits form heteromultimeric GABA-gated cation channels. Br J Pharmacol. 2004;142:409–13.
Xu F, Hollins B, Landers TM, McClintock TS. Molecular cloning of a lobster Gbeta subunit and Gbeta expression in olfactory receptor neuron dendrites and brain neuropil. J Neurobiol. 1998;36:525–36.
Jiménez-Vázquez EN, Díaz-Velásquez CE, Uribe RM, Arias JM, García U. Molecular cloning and expression of a GABA receptor subunit from the crayfish Procambarus clarkii. J Neurosci Res. 2016;94:190–203.
Padgett CL, Slesinger PA. GABAB receptor coupling to G-proteins and ion channels. Adv Pharmacol. 2010;58:123–47.
Ren GR, Folke J, Hauser F, Li S, Grimmelikhuijzen CJP. The A- and B-type muscarinic acetylcholine receptors from Drosophila melanogaster couple to different second messenger pathways. Biochem Biophys Res Commun. 2015;462:358–64.
Collin C, Hauser F, Gonzalez de Valdivia E, de Valdivia EG, Li S, Reisenberger J, et al. Two types of muscarinic acetylcholine receptors in Drosophila and other arthropods. Cell. Mol. Life Sci. 2013;70:3231–42.
Florey E. Acetylcholine in Invertebrate Nervous Systems. Can J Biochem Physiol. 1963;41:2619–26.
Hildebrand JG, Townsel JG, Kravitz EA. Distribution of acetylcholine, choline, choline acetyltransferase and acetylcholinesterase in regions and single identified axons of the lobster nervous system. J Neurochem. 1974;23:951–63.
Marder E, Paupardin-Tritsch D. The pharmacological profile of the acetylcholine response of a crustacean muscle. J Exp Biol. 1980;88:147–59.
Marder E. Cholinergic motor neurones in the stomatogastric system of the lobster. J Physiol. 1976;257:63–86.
van Nierop P, Bertrand S, Munno DW, Gouwenberg Y, van Minnen J, Spafford JD, et al. Identification and functional expression of a family of nicotinic acetylcholine receptor subunits in the central nervous system of the mollusc Lymnaea stagnalis. J Biol Chem. 2006;281:1680–91.
Touroutine D, Fox RM, Von Stetina SE, Burdina A, Miller DM, Richmond JE. acr-16 encodes an essential subunit of the levamisole-resistant nicotinic receptor at the Caenorhabditis elegans neuromuscular junction. J Biol Chem. 2005;280:27013–21.
Phelan P. Innexins: Members of an evolutionarily conserved family of gap-junction proteins. Biochim Biophys Acta. 2005;1711:225–45.
Shruti S, Schulz DJ, Lett KM, Marder E. Electrical coupling and innexin expression in the stomatogastric ganglion of the crab Cancer borealis. J Neurophysiol. 2014;112:2946–58.
Ducret E, Alexopoulos H, Le Feuvre Y, Davies JA, Meyrand P, Bacon JP, et al. Innexins in the lobster stomatogastric nervous system: cloning, phylogenetic analysis, developmental changes and expression within adult identified dye and electrically coupled neurons. Eur J Neurosci. 2006;24:3119–33.
Getting PA. Mechanisms of pattern generation underlying swimming in Tritonia. II Network reconstruction. J Neurophysiol. 1983;49:1017–35.
Castellucci V, Kandel E. Presynaptic facilitation as a mechanism for behavioral sensitization in Aplysia. Science. 1976;194:1176–8.
Baxter DA, Byrne JH. Feeding behavior of Aplysia: A model system for comparing cellular mechanisms of classical and operant conditioning. Learn Mem. 2006;13:669–80.
Moroz LL, Edwards JR, Puthanveettil SV, Kohn AB, Ha T, Heyland A, et al. Neuronal transcriptome of Aplysia: Neuronal compartments and circuitry. Cell. 2006;127:1453–67.
Senatore A, Edirisinghe N, Katz PS. Deep mRNA sequencing of the Tritonia diomedea brain transcriptome provides access to gene homologues for neuronal excitability, synaptic transmission and peptidergic signalling. PLoS One. 2015;10:e0118321.
Kravitz EA. Serotonin and aggression: insights gained from a lobster model system and speculations on the role of amine neurons in a complex behavior. J Comp Physiol A. 2000;186:221–38.
Grundfest H, Reuben JP, Rickles WH. The electrophysiology and pharmacology of lobster neuromuscular synapses. J Gen Physiol. 1959;42:1301–23.
Krasne FB, Edwards DH. Modulation of the crayfish escape reflex--physiology and neuroethology. Integr Comp Biol. 2002;42:705–15.
Stepanyan R, Day K, Urban J, Hardin DL, Shetty RS, Derby CD, et al. Gene expression and specificity in the mature zone of the lobster olfactory organ. Physiol Genomics. 2006;25:224–33.
Towle DW, Smith CM. Gene discovery in Carcinus maenas and Homarus americanus via expressed sequence tags. Integrative Comp Biol. 2006;46:912–18.
McClintock TS, Ache BW, Derby CD. Lobster olfactory genomics. Integr Comp Biol. 2006;46:940–7.
Harms L, Frickenhaus S, Schiffer M, Mark FC, Storch D, Held C, et al. Gene expression profiling in gills of the great spider crab Hyas araneus in response to ocean acidification and warming. BMC Genomics. 2014;15:789.
Harms L, Frickenhaus S, Schiffer M, Mark FC, Storch D, Pörtner HO, et al. Characterization and analysis of a transcriptome from the boreal spider crab Hyas araneus. Comp Biochem Physiol Part D Genomics Proteomics. 2013;8:344–51.
Ventura T, Fitzgibbon QP, Battaglene SC, Elizur A, Medina M, Anger K, et al. Redefining metamorphosis in spiny lobsters: molecular analysis of the phyllosoma to puerulus transition in Sagmariasus verreauxi. Sci Rep. 2015;5:13537.
Ventura T, Cummins SF, Fitzgibbon Q, Battaglene S, Elizur A. Analysis of the central nervous system transcriptome of the Eastern rock lobster Sagmariasus verreauxi reveals its putative neuropeptidome. PLoS One. 2014;9:e97323.
Christie AE, Chi M. Prediction of the neuropeptidomes of members of the Astacidea (Crustacea, Decapoda) using publicly accessible transcriptome shotgun assembly (TSA) sequence data. Gen Comp Endocrinol. 2015;224:38–60.
Christie AE. Prediction of Scylla olivacea (Crustacea; Brachyura) peptide hormones using publicly accessible transcriptome shotgun assembly (TSA) sequences. Gen Comp Endocrinol. 2016;230–231:1–16.
Buckley SJ, Fitzgibbon QP, Smith GG, Ventura T. In silico prediction of the G-protein coupled receptors expressed during the metamorphic molt of Sagmariasus verreauxi (Crustacea: Decapoda) by mining transcriptomic data: RNA-seq to repertoire. Gen Comp Endocrinol. 2016;228:111–27.
MacLean JN, Zhang Y, Johnson BR, Harris-Warrick RM. Activity-independent homeostasis in rhythmically active neurons. Neuron. 2003;37:109–20.
Krenz W-D, Parker AR, Rodgers E, Baro DJ. Monoaminergic tone supports conductance correlations and stabilizes activity features in pattern generating neurons of the lobster. Panulirus interruptus. Front Neural Circuits. 2015;9:63.
Marder E. Variability, compensation, and modulation in neurons and circuits. Proc Natl Acad Sci U S A. 2011;108(Suppl):15542–8.
Hamood AW, Marder E. Animal-to-animal variability in neuromodulation and circuit function. Cold Spring Harb Symp Quant Biol. 2014;79:21–8.
Marder E, O’Leary T, Shruti S. Neuromodulation of circuits with variable parameters: Single neurons and small circuits reveal principles of state-dependent and robust neuromodulation. Annu Rev Neurosci. 2014;37:329–46.
Temporal S, Lett KM, Schulz DJ. Activity-dependent feedback regulates correlated ion channel mRNA levels in single identified motor neurons. Curr Biol. 2014;24:1899–904.
Ball JM, Franklin CC, Tobin A-E, Schulz DJ, Nair SS. Coregulation of ion channel conductances preserves output in a computational model of a crustacean cardiac motor neuron. J Neurosci. 2010;30:8637–49.
Tobin AE, Cruz-Bermúdez ND, Marder E, Schulz DJ. Correlations in ion channel mRNA in rhythmically active neurons. PLoS One. 2009;4:2–9.
This work was supported by NIH grant MH46742 (EM and DJS), NIH grant NIGMS 5T32GM008396 (support for AJN), and a grant from the University of Missouri Research Board (DJS).
Availability of data and materials
All sequence data has been deposited with GenBank, and Accession numbers for the TSA, SRA, and individual sequence submissions are provided in the manuscript and accompanying tables. Matrices used to generate trees have been deposited in TreeBase.
Conceived the study: DJS, EM. Performed wet lab experiments: AJN, KML, CMD, BJL, VBG. Design of RNA-seq: DJS. Bioinformatic analysis: AJN, DJS. Data analysis: AJN, DJS. Interpretation of results: AJN, DJS, EM. Wrote the manuscript: AJN, DJS, EM. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate