- Research article
- Open access
- Published:
Transcriptomic analysis of crustacean neuropeptide signaling during the moult cycle in the green shore crab, Carcinus maenas
BMC Genomics volume 19, Article number: 711 (2018)
Abstract
Background
Ecdysis is an innate behaviour programme by which all arthropods moult their exoskeletons. The complex suite of interacting neuropeptides that orchestrate ecdysis is well studied in insects, but details of the crustacean ecdysis cassette are fragmented and our understanding of this process is comparatively crude, preventing a meaningful evolutionary comparison. To begin to address this issue we identified transcripts coding for neuropeptides and their putative receptors in the central nervous system (CNS) and Y-organs (YO) within the crab, Carcinus maenas, and mapped their expression profiles across accurately defined stages of the moult cycle using RNA-sequencing. We also studied gene expression within the epidermally-derived YO, the only defined role for which is the synthesis of ecdysteroid moulting hormones, to elucidate peptides and G protein-coupled receptors (GPCRs) that might have a function in ecdysis.
Results
Transcriptome mining of the CNS transcriptome yielded neuropeptide transcripts representing 47 neuropeptide families and 66 putative GPCRs. Neuropeptide transcripts that were differentially expressed across the moult cycle included carcikinin, crustacean hyperglycemic hormone-2, and crustacean cardioactive peptide, whilst a single putative neuropeptide receptor, proctolin R1, was differentially expressed. Carcikinin mRNA in particular exhibited dramatic increases in expression pre-moult, suggesting a role in ecdysis regulation. Crustacean hyperglycemic hormone-2 mRNA expression was elevated post- and pre-moult whilst that for crustacean cardioactive peptide, which regulates insect ecdysis and plays a role in stereotyped motor activity during crustacean ecdysis, was elevated in pre-moult.
In the YO, several putative neuropeptide receptor transcripts were differentially expressed across the moult cycle, as was the mRNA for the neuropeptide, neuroparsin-1. Whilst differential gene expression of putative neuropeptide receptors was expected, the discovery and differential expression of neuropeptide transcripts was surprising. Analysis of GPCR transcript expression between YO and epidermis revealed 11 to be upregulated in the YO and thus are now candidates for peptide control of ecdysis.
Conclusions
The data presented represent a comprehensive survey of the deduced C. maenas neuropeptidome and putative GPCRs. Importantly, we have described the differential expression profiles of these transcripts across accurately staged moult cycles in tissues key to the ecdysis programme. This study provides important avenues for the future exploration of functionality of receptor-ligand pairs in crustaceans.
Background
Arthropods are characterised by their rigid exoskeleton that, although fundamental to their evolutionary success, must be shed periodically so that they can develop and grow. The moulting process, termed ecdysis, is achieved via an innate sequence of behaviours and stereotyped motor activity, precisely co-ordinated by the interaction of steroid moulting hormones (ecdysteroids) and a complex suite of interacting neuropeptides. The insect ecdysis programme is arguably the best understood endocrine cascade among all invertebrates and a consensus model is now emerging (reviews: [1, 2]). Briefly, new cuticle development is driven by ecdysteroids, which are synthesised and released from the prothoracic glands following stimulation by the brain-derived prothoracicotropic hormone (PTTH). Subsequently, ecdysteroids initiate ecdysis triggering hormone (ETH) production in so-called ‘Inka’ cells (tracheal endocrine cells) and induce central nervous system (CNS) sensitivity to ETH by promoting ETH receptor (ETHR) expression. Two alternatively spliced ETHRs (ETHR-A & ETHR-B) are expressed in discrete populations of neurons that are thought to be differently involved in pre-ecdysis and ecdysis [3,4,5]. ‘Inka’ cells, however only gain secretory competence when ecdysteroid titre declines below a threshold level after which they release ETH, which in turn induces eclosion hormone (EH) release from the CNS. Together, ETH and EH enter a positive feedback loop culminating in a massive release of both neuropeptides into circulation [6, 7]. In the tobacco hornworm, Manduca sexta, the neuropeptide corazonin (CRZ), a paralog of gonadotropin-releasing hormone, initiates the release of ETH from ‘Inka’ cells [8]. However, this is not the case in Drosophila and in some insects corazonin is entirely absent [9]. Elevated hemolymph titres of ETH and EH provides a clear ecdysis signal and commits the animal to the process by evoking release of the cyclic nonapeptide, crustacean cardioactive peptide (CCAP), which starts the ecdysis motor programme and terminates pre-ecdysial behaviours [10, 11]. Post-ecdysis cuticle maturation, including wing inflation and cuticle tanning, is regulated by the peptide bursicon, co-released with CCAP from a sub-set of CCAP neurons in the abdominal ganglia [12,13,14]. Our knowledge on the species-specific actions of ETH, EH, CCAP, and bursicon is under constant review, but the emerging picture is that whilst the functions of these neuropeptides are mostly conserved throughout insects, they may have subtly different roles in different taxonomic groups and at different developmental stages.
Whilst our understanding of the control of insect ecdysis has advanced in recent years, it is relatively poorly described in crustaceans. Nevertheless, the regulation of insect and crustacean ecdysis has a common origin in the neuroendocrine system of the crustacean-like ancestor from which these clades evolved [15, 16] and so crustacean endocrinology has benefitted from comparative work done on insects. In crustaceans, ecdysteroid synthesis occurs in Y-organs (YO) that are homologous with the prothoracic glands of insects, and via a common biosynthesis pathway [17]. Perhaps the defining difference between insect and crustacean cuticular development and ecdysis is that whilst ecdysteroids ultimately co-ordinate the moult cycle in both insects and crustaceans, ecdysteroid biosynthesis in the YO is under negative regulation by the eyestalk-derived neuropeptides, crustacean hyperglycemic hormone (CHH) and moult-inhibiting hormone (MIH) [18]. In crustaceans, CHH and MIH, which are part of the structurally related (but functionally diverse) CHH-superfamily [19] are synthesised in large perikarya concentrated in the so-called X-organ of the eyestalk medulla terminalis and released (in the case of MIH) in a pulsatile fashion by exocytosis from a coalescence of axons emanating from these cells [20], the sinus gland, into the hemolymph. After circulating in the hemolymph, CHH and MIH bind to specific receptors (putatively, GPCRs) on the YO to effect ecdysteroid inhibition [18] during intermoult stages (stages C-D1). Hence the discovery in 1953 in fiddler crabs (Uca) that removal of eyestalks releases this inhibition and can lead to precocious moulting [21]. The temporal release of CHH and MIH implicated in ecdysteroid inhibition in crustaceans is complicated by species specificity and we currently lack a consensus model of the events leading to moulting. This issue is discussed in detail by Webster (2014) [18]. Whatever the precise order of events, YO inhibition must be released in premolt and is followed by elevated ecdysteroid synthesis in stages D1-D3. Downstream of ecdysteroid synthesis our understanding of crustacean moult control is fragmentary at best and there is currently no consensus on the endocrine cascade that culminates in ecdysial behaviour and cuticle tanning. Thus far we know that, during ecdysis (stage E) CHH is dramatically released from paraneurons in the fore and hind-gut whilst CCAP and bursicon are co-released from the pericardial organs, and concomitant with the onset of ecdysis motor programme that precedes post-ecdysial tanning and cuticle hardening [18, 22] (during stages A-C). Thus, the activities of CCAP and bursicon may be generally conserved between crustaceans and insects [18, 22]. Broadly then, our knowledge of the neuropeptides involved in the regulation of crustacean ecdysis, such as CCAP and bursicon, indicate that some elements are conserved across arthropods whilst others are strikingly different (e.g. ecdysteroid regulation). Our inadequate understanding of key elements (and their roles) of the neuroendocrine control of crustacean ecydsis precludes a meaningful evolutionary comparison with the process in insects.
Until recently, the discovery and characterization of neuropeptides relied on the laborious purification of native peptides from animal tissues coupled with manual peptide sequencing or mass spectrometry (e.g. [23, 24]). Today, next generation sequencing (NGS) strategies have become routine and relatively inexpensive and coupled with emerging informatics pipelines have accelerated peptide mRNA transcript discovery and characterization [25]. This new era of whole transcriptome sequencing has generated a surge in neuropeptide discovery in non-model species, including crustaceans, overcoming their previous genetic intractability (e.g. [26,27,28,29]). A consequence of this post-genomic era is a widening gap between the description and characterization of conceptual peptide genes and their cognate peptide products and in vivo functions. Through NGS and transcriptome mining within crustaceans, orthologs of insect neuropeptides have been documented [27], but functional roles for these peptides remain lacking. Similarly, whilst the receptors for numerous neuropeptides have been functionally demonstrated for insects, only a few G protein-coupled receptors (GPCRs) have been de-orphanised among crustaceans [30, 31] and so identification of their GPCRs lags considerably behind. Indeed, only recently have efforts been focussed on revealing the structural identity of putative GPCRs in a crustacean, the spiny lobster Sagmariasus verreauxi [32]. Here, we sought to shed light on the involvement of neuropeptides in crustacean ecdysis by using RNA-sequencing to identify neuropeptides and putative neuropeptide GPCRs in the CNS and YO of the green shore crab, Carcinus maenas and against the insect system as a framework for our exploration. Thus, we describe a comprehensive suite of neuropeptide and receptor transcripts in C. maenas and define their expression profiles across the moult cycle at precisely defined moult stages. This study is the first exploration of a multi-tissue crustacean neurotranscriptome across the moult cycle with high temporal precision and therefore provides the opportunity to elucidate putative neuropeptide signaling during the crustacean ecdysial programme and draw comparisons with the current insect model.
Methods
Carcinus maenas were collected from the Menai Straits, UK using baited traps and held in a flow through aquaria at Bangor University until use. Aquarium conditions were: water temperature = 16 °C, salinity =35 (full seawater), light:dark cycle = 12:12. Aquarium conditions did not replicate ambient tidal cycles. Crabs were fed ad libitum twice per week and were kept for no more than 4 weeks. Crabs were collected and dissected throughout the spring and summer months when moulting is prevalent.
Total RNA extraction, cDNA library preparation, and sequencing
Total RNA was harvested from crabs dissected across the moult cycle at stages A-B, C3–4, D1, D3, D4 (stages according to Phlippen et al. [33]) for n = 5 animals (biological replicates) per moult-stage. Three sets of tissues were collected for transcriptome sequencing and were subsequently analysed independently due to the considerable volume of data which sequencing yielded (see Additional file 1: Table S5, S6, and S7, for further details on sequenced data sets). The first consisted of central nervous system tissues (CNS: eyestalk ganglia [ESG], cerebral ganglion [CG], and ventral ganglion [VG] comprising the fused sub-oesophageal [SOG], thoracic [TG], abdominal ganglia [AG]); the second, pairs of Y-organs (YO); and the third consisted of inter-moult (C3–4) pairs of YOs and epidermis from the branchiostegite region of the carapace (overlying the pleural suture posterior of the YO). These tissues were dissected from ice-anaesthetised C. maenas under ice-cold physiological saline [34], using dedicated instruments for each tissue. Tissues were snap frozen in liquid nitrogen before storage at − 80 °C until use. Tissues were homogenised in TRIzol reagent in a Qiagen TissueLyser LT (Qiagen, UK) for 1 min (CNS and YO tissues) or 3 min (epidermis) at 50 Hz with a 5 mm stainless steel bead and RNA extracted according to the manufacturer’s protocol except additional ethanol washes were performed prior to resuspension of the RNA pellet in DEPC-treated water. RNA integrity was assessed by gel electrophoresis (Additional file 1: Figure S1) whilst yield and purity were measured spectrophotometrically with a NanoDrop ND2000™ (Thermo Scientific, UK). cDNA library preparation was done using Illumina TruSeq RNA sample preparation reagents according to the manufacturer’s protocol. For CNS tissue cDNA library preparation, equal quantities (1 μg) of total RNA taken from separate neural tissues (ESG, CG, VG) were pooled. Library quality was assessed by gel electrophoresis and quantified using a Qubit 2.0 Fluorometer (Invitrogen™, Thermo Fisher, UK) prior to sequencing. Paired-end 126-bp reads were sequenced across multiple lanes on an Illumina HiSeq 2500 platform. Library preparation, quality checks and sequencing were done in-house at the Translational Genomics Facility, IBERS, Aberystwyth University.
Transcriptome assembly, characterization, and differential gene expression analysis
Comparisons of various de novo assembly software packages demonstrate that no single assembler is universally superior [35, 36]. To thoroughly evaluate the potential of the sequenced raw reads, multiple assemblers were used for de novo transcriptome assembly of CNS reads and the resulting transcriptomes compared: Trinity v2.0.6 [37], Oases v0.2.8 [38] (which used Velvet v1.2.10 [39]), IDBA-tran v1.1.1 [40], and BinPacker v1.1 [41]. Prior to assembly, a 13-bp head crop of raw reads was done [42] followed by trimming and in silico read normalisation using the Trinity options --trimmomatic [43] and --normalize_reads. Both options were run with default settings (trimmomatic: phred33, SLIDINGWINDOW:4:5 LEADING:5 TRAILING:5 MINLEN:25, based on [44]). Trimmed and normalised reads were assembled with default settings in each of the four assemblers (k-mer length: Trinity k = 25; Oases k = 23, 27, 31; IDBA-tran k = 20, 30, 40, 50, 60; BinPacker k = 25). Comparisons of these four transcriptome assemblies revealed the Trinity assembly to be superior in terms of the number and completeness of C. maenas neuropeptide transcripts assembled as well as number and completeness of UniProt/Swiss-Prot BLAST hits (for details see Additional file 1: Table S1, Figures S2 and S3). Consequently, Trinity was used for the assembly of transcriptomes in the current study. Transcriptomes for CNS and YO tissues sampled across the moult cycle, and for inter-moult (C3–4) YO and epidermis tissues, were assembled separately. For each transcriptome, Trinity transcripts were clustered into ‘genes’ using Corset v1.06 [45], which clusters transcripts based on sequence similarity and read counts (generated by Bowtie v1.0.1 [46] with multi-mapping enabled). These Corset-clustered Trinity transcriptomes (herein referred to as the Carma_CNS-transcriptome, Carma_YO-transcriptome, and Carma_YO_vs_Epi-transcriptome) were annotated using Trinotate v2.0.1 and associated software (TransDecoder v2.0.1, BLAST v2.2.31+, HMMER v3.1, SignalP v4.1, TMHMM v2.0, RNAMMER v1.2, TrinotateR). The process of Corset clustering, which used Bowtie generated mapping data, outputs a counts file for use in differential expression analysis [45]. For Carma_CNS- and Carma_YO-transcriptomes, differential cluster expression was performed using the edgeR package [47] (with the glm method to enable multiple comparisons) and by comparing inter-moult (C3–4) stage count data with that from each of the other four moult stages (A-B, D1, D3, and D4). For the Carma_YO_vs_Epi-transcriptome, differential cluster expression was performed by comparing count data between YO and epidermis tissues, again using the edgeR package [47]. Bowtie mapping, Corset clustering and counting, and edgeR differential expression analysis followed the example Corset pipeline at: github.com/Oshlack/Corset/wiki. The three transcriptomes assembled as part of this Transcriptome Shotgun Assembly project have been deposited at DDBJ/EMBL/GenBank under the accessions: GFXF00000000 (Carma_CNS-transcriptome), GFYV00000000 (Carma_YO-transcriptome), and GFYW00000000 (Carma_YO_vs_Epi-transcriptome). The transcriptomes have also been made available as publically accessible BLAST databases at: sequenceserver.ibers.aber.ac.uk. Raw read sequence files used to generate these assemblies are archived in the NCBI SRA archive under BioProject PRJNA400568. Trinotate annotation reports for Carma_CNS-, Carma_YO-, and Carma_YO_vs_Epi-transcriptome are provided in Additional file 2 as well as transcripts for neuropeptides and putative neuropeptide receptors mined from Carma_CNS- and Carma_YO-transcriptomes.
Neuropeptide and receptor discovery
The Carma_CNS- and Carma_YO-transcriptomes were mined for contigs coding for putative neuropeptides using local tBLASTn searches with sequences from C. maenas and related species as search terms (performed in BioEdit software, [48]). Short protein motifs in open reading frames (ORFs) identified by TransDecoder were also searched in BioEdit. Peptide sequences used as search terms were taken from the NCBI [49] database and from recent literature [26, 27, 50, 51]. BLAST e-value and sequence identity thresholds considered significant between query peptide sequences and BLAST hits were assessed case-by-case. Contigs mined as putative neuropeptides were translated using the online tool ExPASy Translate (web.expasy.org/translate, [52]) and signal peptides identified with SignalP v.4.1 Server (set to ‘Sensitive’ to reproduce SignalP v.3.0’s sensitivity, [53]). Convertase cleavage sites were predicted following [26, 51, 54] and basic amino acid residues removed by carboxypeptidase were identified. Neuropeptide prepro-hormone processing and deduced mature peptides and precursor-related peptides (PRPs) have been made available in Additional file 3. Cysteine-cysteine disulfide bridges were predicted with the online tool: disulfind.dsi.unifi.it. Contigs coding for putative neuropeptide GPCRs were also mined using local tBLASTn searches with sequences taken from the NCBI database for Drosophila and Tribolium casteneum. Whilst there are studies which have sought to identify crustacean neuropeptide GPCRs based on sequence similarity to known neuropeptide GPCRs (e.g. [29, 32, 51]), until recently no crustacean neuropeptide GPCR had been de-orphanised [30, 31]. For this reason predicted receptors were identified as a result of searching using search terms for Tribolium casteneum and Drosophila receptors, for which functional de-orphanisation has been achieved. GPCR transmembrane helix domains were predicted using TMHMM server v.2.0 (cbs.dtu.dk/services/TMHMM, [55, 56]). Amino acid alignments were done in BioEdit via the ClustalW method [57]. Phylogenetic trees were constructed using the Neighbor-Joining method (1000 bootstrap replicates) [58] in MEGA6 [59]. Drosophila melanogaster metabotropic Glutamate Receptor (CG11144 mGLuR) was used as an outgroup to root phylogenetic trees for class A (rhodopsin-like), class B (secretin-like), and leucine-rich repeat-containing GPCRs (LGR) trees. Drosophila melanogaster trissin (NP_650471.1) was used as an outgroup to root the tree for the ETH, carcikinin, and carcikinin-like peptides phylogenetic tree.
Quantitative RT-PCR
Quantitative RT-PCR assays were developed and performed according to previously published methods [60] using Taqman™ MGB probes and in close adherence to MIQE guidelines [61]. In brief, standard curves were made by generating complementary RNA in vitro with T7 RNA polymerase (MegaShort Script™, Ambion®, UK) using PCR derived DNA templates amplified with T7 phage promotor sequence flanked gene specific primers. All oligonucleotide sequences for cRNA preparations and Taqman assays are detailed in Additional file 1: Table S2. Complementary RNA was purified on 10% 6 M urea PAGE gels and eluted in Elution Buffer (Ambion®, Thermo fisher, UK) overnight at RT before ethanol precipitation. Resulting cRNA was quantified spectrophotometrically, converted to copy number using Avogadro’s constant and diluted in DEPC-treated water and stored at − 80 °C until use. Standard curves were run in the range 109 to 103 copies per reaction. Sample and standard RNA was reverse transcribed using Tetro™ reverse transcriptase (Bioline, UK) according to the manufacturer’s instructions and diluted in water 1:10 for PCR. Assays were run in triplicate on an Applied Biosystems® Quant Studio 12 Flex platform (Thermo-Fisher Scientific, UK) in 10 μL volumes using Sensimix Fast II probe™ qPCR mix with ROX internal reference dye (Bioline, UK). Data were expressed as copies of target mRNA normalized to the geometric mean of the reference genes elongation factor 1-alpha (EF1a) and ubiquitin-conjugating enzyme E2 L3 (UBE2L3) shown to be constitutively expressed across all moult stages with RNAseq analysis (Additional file 1: Figure S4).
End-point PCR
End-point PCRs were performed for cDNAs encoding eclosion hormones 1 and 2, crustacean hyperglycemic hormone 2, and the reference gene elongation factor 1-alpha, to establish the tissue distribution of these transcripts. Tissues were dissected, and total RNA extracted, and RNA reverse transcribed for 5 inter-moult (C4) crabs as outlined above. cDNA of n = 5 individuals was pooled for end-point PCR, which was performed using Bioline myTAQ red mix on a Bio-Rad T100 thermal cycler. PCR conditions were: denature at 95 °C for 3 mins, then 35 cycles of 30 s at 95 °C, 45 s at 60 °C, and 45 s at 72 °C, followed by extension at 72 °C for 5 mins. PCR products were resolved on 1% agarose gels.
Results and discussion
Transcriptome assembly and annotation summary data are detailed in Additional file 1.
Neuropeptide discovery and expression
Carma_CNS-transcriptome
Of the 52 neuropeptide families searched for, representatives of 47 were discovered: 104 transcripts (from 87 clusters) were identified, predicted to encode for 59 full-length peptides or proteins (comprising a ‘start’ methionine, a signal peptide, and a stop codon), 24 C-terminal partial proteins (6 of which contained partial signal peptides), 13 N-terminal partial proteins and 8 internal (lacking both termini) fragments (see Additional files 2 and 3). Analysis of post-translational processing of these conceptual proteins predicts 160 unique mature peptides and 185 precursor-related peptides (PRPs) (see Additional file 3). This is likely an under-estimate of the total number of mature peptides encoded in this transcriptome as prepro-hormones of some neuropeptides contain multiple copies of identical mature peptides (e.g. a HIGSLYamide C-terminal partial protein contained 6 identical mature peptides; see Additional file 3). For Daphnia pulex, 43 genes encoding 73 neuropeptides were predicted in silico; mass spectrometry confirmed 40 of these by mass matching and 30 by fragmentation sequencing [62]. Within the Drosophila genome, some 42 genes encode neuropeptide precursors predicted to produce approximately 75 mature peptides although the translation of many of these transcripts has not been confirmed in vivo [63]. Comparatively, the number of neuropeptide precursors identified, and mature peptides predicted here for C. maenas appears large. However, recent efforts mining decapod transcriptomes have yielded similarly large numbers of transcripts coding for neuropeptide precursors and mature peptides. For example, Christie (2016) identified transcripts from 29 neuropeptide families predicted to produce approximately 112 mature peptides (in addition to 151 PRPs) for C. maenas [26]. Furthermore, Veenstra (2015) identified 58 transcripts encoding neuropeptides from some 41 families for the crayfish Procambarus clarkii [51] and in a substantial transcriptome mining effort, Veenstra (2016) identified 66 ‘genes’ from 47 neuropeptide families for C. maenas as well as comparable values for Scylla paramamosain (59 ‘genes’, 45 neuropeptide families), Eriocheir sinensis (68, 49), Litopenaeus vannamei (70, 49), Macrobrachium rosenbergii (67, 48), Homarus americanus (59, 45), and Procambarus clarkii (63, 49) [27]. Mining efforts for chelicerate transcriptomes and genomes have yielded yet greater numbers of neuropeptide precursors, a result of paralog genes from genome duplications: e.g. 87 and 79 neuropeptide precursors for the spider Stegodyphus mimosarum and the scorpion Mesobuthus martensii, respectively [64].
For C. maenas, the production of 122 mature peptides within the CNS and neuro-endocrine organs has been confirmed by mass spectrometry [65]. Numerous mature peptides predicted in silico here and elsewhere for C. maenas, or confirmed biochemically, are identical or very similar isoforms [26, 65]. Whether such isoforms are physiologically redundant or serve distinct functions is undetermined [63]. For example, the 25 mature allatostatin-As detected in the C. maenas CNS by mass spectrometry comprise the characteristic YXFGLamide C-termini but, differ in N-termini sequence [65]. Here, a single putative allatostatin-A GPCR was discovered (see below, Table 2), suggesting that the 25 or so isoforms may indeed be physiologically redundant, although this has yet to be functionally determined. In this context, it is interesting to note that there is evidence for neuropeptide redundancy in Drosophila [66]. Of course, we should not ignore an alternative scenario that further allatostatin receptors are present but were not predicted in the current analysis. As for all predicted GPCR and their putative ligands, functional de-orphanisation is the only way to clarify the actual situation.
Representatives from allatotropin, androgenic gland specific insulin-like peptide (IAG), DENamide, DXXRLamide, and FXGGXamide families were absent from the Carma_CNS-transcriptome. The discovery of DENamides is so far limited to the cladoceran water flea, D. pulex [62] whilst DXXRLamides and FXGGXamides have been identified only in copepods [67]. Like DENamide, allatotropin is present in D. pulex and the copepods Tigriopus californicus and Lepeophtheirus salmonis and putative allatotropins are present in the transcriptomes of the amphipods, Talitrus saltator (NCBI: PRJNA297565) and Hyalella azteca (NCBI: XM_018171664.1; see Additional file 1: Figure S5). Given that allatotropin is present in molluscs, annelids, and chelicerates, it is presumably an ancient neuropeptide which may have been lost in decapod crustaceans [27]. Androgenic gland specific insulin-like hormone has been identified in a number of decapods, but the same study failed to find it in C. maenas [27]. In studies reported to date [68] IAG was found only in the androgenic gland and it is not therefore, surprising that we did not sequence its transcript in C. maenas and its existence in this species should not be ruled out. We were able to identify a transcript for CHH-MIH-like neuropeptide (the function of which is unknown) previously described in other decapods, but until now absent from C. maenas [27]. The diversity of neuropeptides discovered here from the Carma_CNS-transcriptome is consistent with other studies on C. maenas and comparable with other decapods [26, 27, 65].
Of the 87 neuropeptide contig clusters identified, 9 were excluded from differential gene expression analysis because of very low read numbers (< 1 read per million in > 3 samples) [45]. Our analysis revealed that only 5 neuropeptides were differentially expressed across the moult cycle: carcikinin/ecdysis triggering hormone (ETH), crustacean hyperglycemic hormone-2 (CHH-2), crustacean female sex hormone (CFSH), crustacean cardioactive peptide (CCAP), and neuroparsin-3 (P < 0.05, FDR < 0.05; Table 1).
Carcikinin/ecdysis triggering hormone (ETH)
The transcripts for a neuropeptide, here assigned the name ‘carcikinin’, exhibited the highest magnitude of differential expression of all clusters across the moult cycle in the Carma_CNS-transcriptome (Fig. 1). Three transcripts encoding carcikinin prepro-hormones were identified: two were full prepro-hormone transcripts and one was an N-terminal partial sequence with an N-terminal extension upstream of the signal peptide (Additional file 1: Figure S6; see also Additional files 2 and 3). Such N-terminal extensions were also observed for other prepro-hormones identified here (and that have been reported in other in silico mining studies) though whether these are assembly artefacts or have some biological function is yet undetermined [26]. Of the two deduced full-length precursor peptides, one was 136-aa and the other was truncated at 61-aa (Additional file 1: Figure S6). Both precursors have identical 19-aa signal peptides and 16-aa mature peptides (DAGHFFAETPKHLPRIamide) but differ in their precursor-related peptides (PRPs) (Additional file 1: Figure S6; see also Additional file 3).
In silico expression data revealed that carcikinin transcript abundance was low in post-, inter-, and early pre-moult (moult stages A-B, C3–4, and D1, respectively), but increased dramatically during pre-moult stages, D3 and D4 (Fig. 1). Given the importance of ETH to insect ecdysis, coupled with our intriguing in silico data which showed not only that carcikinin was differentially expressed across the moult cycle but that transcript abundance increased strikingly in late pre-moult, we adopted a more targeted quantitative PCR approach to validate in silico data and measure expression of selected genes of interest within defined regions of the CNS. We observed a considerable increase in carcikinin expression during D3 and D4, consistent with in silico data: two-way ANOVA indicated that there was a significant interaction between tissue type and moult stage on carcikinin expression (F(8, 63) = 2.218, P = 0.0376) with dramatically increased expression in the VG at stages D3 and D4 (Tukey’s post-hoc analysis, P < 0.05; Fig. 1B). This expression pattern across the moult cycle strongly suggests a role for VG derived carcikinin in the regulation of ecdysis within C. maenas and perhaps decapod crustaceans generally (Fig. 1).
Carcikinin has previously been identified as a crustacean orthologue of insect ecdysis triggering hormone (ETH), e.g. [27], and a carcikinin-like neuropeptide identified in the genome of the myriapod, Strigamia maritima, was assigned the name ETH in that study [50]. In a phylogenetic analysis of insect and D. pulex ETHs, crustacean carcikinins, and carcikinin-like peptides mined from the transcriptomes of other arthropods, crustacean carcikinins were clustered in a high confidence branch. A myriapod/chelicerate branch and a separate chelicerate only branch also clustered with high confidence. Insect ETHs were clustered into two high confidence branches, one lepidopteran and the second consisting of dipterans and hymenopterans, which excluded T. casteneum (Fig. 2). The branchiopod crustacean, D. pulex, was excluded from carcikinin/carcikinin-like peptides and ETH branches. Whether differences in the peptide sequences of insect ETHs and carcikinin/carcikinin-like peptides reflect phylogenetic differences between insects and other arthropods, or indicate that ETH and carcikinins/carcikinin-like peptides are distinct peptides is unclear. D. pulex neuropeptides have been found to be more closely related to insect, rather than decapod, homologs and this has been interpreted to confirm that branchiopods are an ancestral group to insects [62]. As such, carcikinin/carcikinin-like peptides may represent a more primitive ETH, common across arthropods. Certainly, amino acid residues within the mature peptide are conserved across arthropods (Fig. 2).
Insect ETH is synthesised and released from peripheral endocrine ‘Inka’ cells located in the trachea and PETH immunochemical localization has revealed peripheral cells in pharate nymphs of the ticks Ixodes ricinus and Rhipicephalus appendiculatus [69]. Given the peripheral localization and expression of these peptides in insects it is possible that a true crustacean ETH homolog has been missed in the current and previous studies that have harvested material from neural tissues; though it is interesting to note that the only ‘ETH’ identified in the myriapod genome is carcikinin-like (Fig. 2). A fully annotated decapod genome is clearly requisite and would be timely.
Eclosion hormones (EHs)
Two eclosion hormone (EH)-like ‘genes’ were identified in our Carma_CNS-transcriptome (EH-1 and EH-2; see Additional files 2 and 3), consistent with other in silico mining studies of this species and other decapods [26, 27]. The presence of two crustacean EHs is in contrast to insects that invariably have only one. In insects, EH is expressed in ventro-median (Vm) cells of the cerebral ganglion, which project axons through the ventral ganglia and their connectives to proctodeal nerves on the hindgut [7]. As alluded to earlier, in the insect ecdysis cascade EH forms a positive feedback loop with ETH, resulting in the near total release of EH and ETH and signaling the start of ecdysis. ETH in turn initiates pre-ecdysis behaviours whilst EH evokes CCAP and bursicon release, thus committing the animal to ecdysis and terminating pre-ecdysial events [11]. However, recent evidence from Drosophila suggests that this view is rather simplistic and that EH may be critical to the initiation of pre-ecdysis behaviours, at least in flies [70].
In the present study, C. maenas EH-2 was represented by a single contig encoding a full length protein but was excluded from differential expression analysis because so few reads were counted for the transcript, indicating extremely low expression levels in the CNS. Endpoint PCR of inter-moult animals revealed that EH-2 mRNA is expressed across multiple tissues at low levels (but see below); EH-1 is also expressed across multiple tissues, but at much higher levels (Additional file 1: Figure S7). EH-1 was not differentially expressed across the moult cycle (F = 1.71, P = 0.18, FDR = 0.70; Fig. 3), but given the importance of EH in insect ecdysis, we thought it worthy of further investigation. In silico gene expression was validated by qPCR which also revealed the primary location of EH-1 expression within the CNS to be the eyestalk ganglia (ESG; Fig. 3B). Two-way ANOVA to test for the effects of moult stage and tissue type on EH-1 expression showed no interaction effect (F(8,63) = 2.04, P = 0.059) but main effects analysis indicated that the mRNA expression in the eyestalk was significantly greater than that in the cerebral and ventral ganglia (F(2,71) = 139.26, P < 0.001; Bonferroni post hoc P < 0.05) (Fig. 3B). The presence of two EHs, coupled with the expression of EH-1 mRNA predominantly in the eyestalks of C. maenas suggests likely differences in the neural architecture of EH producing cells and perhaps EH functioning between insects and crustaceans. Given the fundamental role insect EH plays in insect ecdysis, the localization and functional characterization of crustacean EHs (especially EH-1) and their neural architecture are a priority.
Crustacean cardioactive peptide (CCAP)
Three CCAP transcript variants were discovered and assigned to a single cluster, the translated sequences of which were identical in their first 137-aas (and thus produced identical mature nonapeptides: PFCNAFTGCamide) but differed in their C-termini and therefore the CCAP PRPs that they encode (Additional file 1: Figure S8; see also Additional file 3). CCAP PRPs and other neuropeptide PRPs may serve important physiological functions, but to date little evidence for their roles exists. CCAP mRNA was found to be differentially expressed across the moult cycle (P < 0.05 and FDR < 0.05; Table 1, Fig. 4) and had a strikingly similar expression profile to those of bursicon-α, bursicon-β, and allatostain-CC(part a), which were constitutively expressed across the moult cycle (P-values < 0.05, but FDR > 0.05; Table 1, Fig. 4). Analysis revealed strong positive correlations between CCAP and bursicon-α, −β, and Ast-CC(part a) gene expression (Fig. 4E, F, G). The gene expression profiles observed here for CCAP and bursicon agree with those observed (via qPCR) in other studies of C. maenas [71, 72]. Levels of CCAP and buriscon peptides in the pericardial organs of C. maenas are closely coupled across the moult cycle whilst levels in the fused ventral ganglion are coupled during inter- and pre-moult, but become uncoupled during ecdysis and post-moult [22].
In insects, CCAP release is stimulated by EH and initiates the ecdysis motor programme whilst halting pre-ecdysis behaviours. Bursicon is co-released with CCAP to initiate the ecdysis motor programme and post-ecdysis cuticle tanning and wing inflation [10, 11]. CCAP and bursicon have been reported to perform similar functions in C. maenas [18, 22]. Within C. maenas, CCAP and bursicon are co-localised in paired, segmentally repeating perikarya of the ventral ganglion [22, 72]. Given that CCAP and bursicon expression also appears coupled with that of Ast-CC, it would be interesting to explore whether Ast-CC is co-localised in the CCAP-bursicon neurons of the ventral ganglion. In Drosophila, CAMB neurons (so called because they produce CCAP, Allatostatin-CC, Myoinhibitory peptide, and Bursicon) comprise the ‘master switch’ of ecdysis behaviour following exposure to ETH [73]. Co-localisation of Ast-CC as well as CCAP and bursicon in the paired, segmentally iterated perikarya of the ventral ganglion would indicate that these neurones might be the crustacean equivalent of the fruit fly CAMB neurons and thus, clearly warrant further study.
Crustacean hyperglycemic hormone (CHH) neuropeptide family
The CHH-superfamily comprises structurally related neuropeptides including type-I peptides: crustacean hyperglycemic hormone (CHH) and insect ion transport peptides (ITP); and type-II peptides: moult inhibiting hormone (MIH), vitellogenesis-inhibiting hormone (VIH), and mandibular organ-inhibiting hormone (MOIH) [19, 73,74,75]. CHH neuropeptides are involved in numerous physiological events with reported biological activities including regulation of carbohydrate metabolism, osmo and iono-regulation, moulting, gonad maturation, and methyl farnesoate synthesis by the mandibular organ [19]. Here, two CHHs, an MIH, and (for the first time) a CHH-MIH-like neuropeptide transcript were identified (see Additional files 2 and 3). For each of the CHHs (CHH-1 and CHH-2), two transcript variants were assembled, differing from one-another in their C-terminal regions in the translated sequence (Additional file 1: Figure S9). In a study mining decapod transcriptomes, Veenstra (2016) identified two C. maenas CHHs, naming them CHH-1 and CHH-2. CHH-1 was already known, having been discovered by traditional methods [76], the second (CHH-2) was newly discovered [27]. Here, we have adhered to the naming of CHHs CHH-1 and CHH-2 according to that of Veenstra (2016) [27]. CHH-1 has been the focus of numerous studies investigating the function of CHH within crustaceans (see below). The conceptual CHH-1s identified here correspond to those sequenced from C. maenas pericardial organs (CHH-1a) and the X-organ-sinus gland (CHH-1b) [77] (Additional file 1: Figure S9). The signal peptides and precursor related peptides (CPRP) produced by the two CHH-1 variants correspond to published sequences [27, 77]. To our knowledge, no functional characterisation has been done on CHH-2. Both CHH-2 s have 27-aa signal peptides and 20-aa CHH-precursor-related proteins which have C-terminal amidation signals; however, CHH-2a and b differ in their mature peptides (Additional file 1: Figure S9; see also Additional file 3). Consistent with CHH superfamily peptides, the 83-aa CHH-2a mature peptide has 6 cysteine residues but the 73-aa CHH-2b mature peptide has only 4. The CHH-MIH-like neuropeptide is 128-aa with a 26-aa signal peptide and a 102-aa mature peptide containing 6 cysteine residues (see Additional files 2 and 3).
Three MIH transcript variants were identified from two isoforms, both of which have been reported previously [26]. MIH-a and MIH-b differed from one-another in 9-aa of their signal peptides and 4-aa of their mature peptides (Additional file 1: Figure S10). The second isoform was represented by a C-terminal partial protein (MIH-c; Additional file 1: Figure S10). The cysteines, which form disulfide bridges, are conserved between MIH-a and MIH-b, and the C-terminal fragment of MIH-c (Additional file 1: Figure S10).
Transcript cluster expression analysis revealed CHH-2 to be differentially expressed across the moult cycle with elevated levels of transcript abundance in moult stages A-B and D4 stages (Fig. 5, Table 1). This expression pattern was strikingly different from other CHH super-family neuropeptides which were constitutively expressed across the moult cycle (Fig. 5). Within C. maenas, studies on the function of CHH have concerned CHH-1 only, the two variants of which arise by alternative splicing [19, 20, 77,78,79]. The defining role of CHH is the mobilization of glycogen resulting in hyperglycemia and increased titres of CHH are observed when crustaceans are exposed to physiologically stressful conditions [19]. An elevation in hemolymph CHH titre is also associated with ecdysis: paraneurons in the fore- and hind-gut discharge CHH, which is known to regulate water and ion uptake during ecdysis [79]. CHH also has an inhibitory activity on ecdysteroid synthesis of the YO, though in C. maenas CHH is ~ 10–20 times less potent than MIH [19, 23]. Numerous CHHs have been identified across many decapods species; for example, six were isolated from sinus gland extracts of the kuruma prawn, Marsupenaeus japonicas, five of which showed hyperglycemic activity [80]. Similarly, via a recent transcriptome mining effort, two CHHs were reported for C. maenas (also identified here: CHH-1 and CHH-2) and the brachyuran crabs Scylla paramamosain and Eriocheir sinensis, whilst 3 were found in Macrobrachium rosenbergii and Homarus americanus, and 7 in Litopenaeus vannamei [27]. The physiological functions for numerous CHHs have been studied in a variety of decapods and each appear to perform similar roles, although the potency may vary with isoform (for detailed review, see [19]).
In contrast to CHH-2, gene expression patterns of CHH-1a, CHH-1b, CHH-MIH-like, MIH-a,b, MIH-c were similar to each other (Fig. 5). This could reflect different roles for CHH-2 relative to CHH-1 and coupled with differential expression across the moult cycle, with elevated expression during A-B and D4, may suggest a role for CHH-2 in ecdysis regulation. Endpoint PCR revealed that within the CNS, CHH-2 is expressed in the CG and VG, but not the ESG (Additional file 1: Figure S7). Given that the primary site of CHH-1 expression is the X-organ of the ESG, this crude endpoint tissue distribution has identified likely differences in the neural architecture of CHH-1 and CHH-2 synthesising perikarya. Given the differing expression dynamics identified here between CHH-2 and other CHH-superfamily peptides (CHH-1/CHH-MIH-like/MIH), coupled with likely differences in peptide synthesis location, studies of physiological functions of newly identified CHHs (CHH-like neuropeptides) may find diverse roles that differ from those established for CHH.
Crustacean female sex hormone (CFSH) and Neuroparsin-3
The recently discovered CFSH is implicated in the development of female reproductive characteristics [81]. We identified transcripts coding for CFSH and two CFSH-like neuropeptides in the Carma_CNS-transcriptome (Additional file 1: Figure S11; see also Additional files 2 and 3), the first time the deduced sequences of the CFSH-like neuropeptides have been reported for this species [27]. The deduced C. maenas CFSH prepro-hormone is a 225-aa peptide with a 24-aa signal peptide. It is cleaved at a KR dibasic cleavage site to produce a 32-aa CFSH-PRP and a 166-aa mature CFSH, which has 8 cysteine residues predicted to form 4 disulfide bridges (C46-C150, C80-C112, C105-C119, C107-C148) consistent with that reported for Callinectes sapidus CFSH [81]. CFSH-like 1 prepro-hormone is a 297-aa protein with a 22-aa signal peptide, a 102-aa CFSH-like 1 PRP and a 170-aa mature CFSH-like 1 peptide which has 10 cysteine redisues predicted to form 5 disulfide bridges (C4-C53, C12-C84, C107-C109, C114-C123, C152-C154). This deduced prepro-hormone also has a 12-aa N-terminal extension immediately upstream the start of its signal peptide (see Additional files 2 and 3). The CFSH-like 2 prepro-hormone codes for a 219-aa protein with a 26-aa signal peptide and a 193-aa mature peptide with 8 cysteine residues predicted to form 4 disulfide bridges (C73-C130, C105-C128, C135-C176, C143-C174). CFSH-like 2 prepro-hormone does not contain a PRP. CFSH was differentially expressed across the moult cycle (F = 6.40, P = 0.001, FDR = 0.034, Table 1); however, expression levels were significantly different between males and females, being expressed at only low levels in males -consistent with previous findings [81]. Two-way ANOVA showed that the effects of moult stage (F(4,19) = 2.958, P = 0.047) and sex (F(1,19) = 10.055, P = 0.005) on gene expression were significant, the interaction between these factors could not be investigated due to the unbalanced nature of the data. Hence, we must treat our observation that CFSH expression varies across the moult cycle with due caution.
Neuroparsins were first identified as anti-gonadotropic factors which delay vitellogenesis in insects [82,83,84]. RNA interference experiments confirmed an inhibitory effect of neuroparsin on vitellogenesis and oocyte maturation in female desert locusts, Schistocerca gregaria [85]. In contrast, neuroparsin-like ovary ecdysteroidogenic hormones (OEHs) are gonadotropic and ecdysteroidogenic in the mosquito Aedes aegypti and ovary maturing parsins (OMPs) of Locusta migratoria and Schistocerca gregaria similarly stimulate oocyte growth, probably by ovarian ecdysone synthesis stimulation [86,87,88]. In the shrimp Metapenaeus ensis, RNAi pointed to a role for neuroparsin in the maturation of oocytes by stimulating vitellogenin production in the hepatopancreas, similar to the roles of OEH and OMP [89]. Neuroparsins have also been used as markers of phase change in gregarious locusts, being differentially expressed between phases. In C. maenas neuroparsin-3 was differentially expressed across the moult cycle, elevated expression was evident in D1 (F = 7.03, P = 0.001, FDR = 0.024, Table 1), though this finding is subject to the same caveat as above regarding the unbalanced numbers of males and females per moult stage. Two-way ANOVA was unable to assess interaction because of the unbalanced nature of the data, but indicated that gene expression was affected by moult stage (F(4,19) = 3.931, P = 0.017), but not sex (F(1,19) = 0.379, P = 0.545).
Carma_YO-transcriptome
The YO is an endocrine tissue, derived from the epidermis and as such expression of neuropeptide transcripts was not expected. To our surprise however, transcripts representative of peptides from across 24 families were identified: 37 transcripts (from 34 clusters) were found to encode 19 full-length peptides or proteins (comprising a ‘start’ methionine, a signal peptide, and a stop codon), 10 C-terminal partial proteins, 5 N-terminal partial proteins and 2 internal (lacking both termini) fragments (see Additional file 2). Six of these neuropeptides were expressed at levels > 10 CPM; including neuroparsin-1, CHH-1, inotocin/vasopressin, EH-2, neuroparsin-3, and neuroparsin-4. Notably, neuroparsin-1 mRNA was differentially expressed within the YO across the moult cycle (Fig. 6). Thus, the YO appears to be a source of synthesis (and perhaps release site) for a number of neuroparsins which, as mentioned previously, are known to play roles in reproductive physiology. Neuroparsin-1 expression was low during inter-moult and early pre-moult (C3–4, D1) but increased through late pre-moult, peaking in D4 and remaining elevated during post-moult (A-B, Fig. 6), therefore marking this YO derived transcript as having a putative role in the ecdysis cassette.
Within the Carma_CNS-transcriptome, EH-2 was rarely expressed and end-point PCR using RNA from inter-moult (C4) crabs confirmed low expression across tissues (Additional file 1: Figure S7). However, RNAseq within the YO revealed high expression of EH-2 in pre-moult (D3, D4), though this was not found to be differentially expressed across the moult cycle, presumably due to high variability between biological replicates (Fig. 6). These findings suggest that the two EHs within crustaceans are likely expressed in discrete tissues, perhaps with discrete physiological functions.
Putative neuropeptide receptor discovery and expression
Carma_CNS-transcriptome
A total of 66 putative neuropeptide GPCRs were discovered within the C. maenas CNS: 100 transcripts (from 72 clusters) were identified, encoding 77 full-length proteins, 13 C-terminal partial proteins, nine N-terminal partial proteins, and one likely internal fragment (see Additional file 2). Of the 77 full-length protein coding sequences discovered, 63 were predicted to contain 7 transmembrane helixes (predicted by tmhmm). In addition to these putative GPCRs, one insulin receptor tyrosine kinase (de-orphanised as a receptor for IAG [90] and, since it is found in both sexes, likely also a receptor for other insulin-like ligands [91]) three insulin receptor-like tyrosine kinases, and two venus kinase receptors (recently deorphanised in the mosquito Aedes aegypti as neuroparsin receptors [92]) were identified (see Additional file 2). Within the Drosophila genome, 49 neuropeptide GPCRs are present, 35 of which have been functionally characterised [93]. The highly conserved nature of arthropod (and particularly insect and crustacean) neuropeptidomes enabled us to predict putative C. maenas neuropeptide receptors based on sequence similarity to those de-orphanised in insects [27, 64]. On this basis, the 66 putative C. maenas neuropeptide GPCRs have been assigned putative ligands from 33 neuropeptide families; six C. maenas GPCRs are homologous to orphan insect GPCRs and thus, at present, cannot be assigned a ligand: neuropeptide B3 R, CG33639-like, MOODY, MOODY-like, neuropeptide A47 R, and neuropeptide A47 R-like. Of the 66 putative neuropeptide receptors, 50 were class A (rhodopsin-like), seven were class B (secretin-like), and nine were leucine-rich repeat-containing GPCRs (LGR) (Figs. 7, 8, 9) (see also Additional file 2). Of the neuropeptides present with characterized GPCRs, only agatoxin-like peptide and periviscerokinin lacked putative GPCRs. Owing to the lack of characterized receptors for CCRFamide, CFSH, HIGSLYamide, Hyrg, and Orcokinin, no putative receptors for these were identified in the present study. Red pigment concentrating hormone is related to adipokinetic hormone and is the putative ligand of an AKH receptor; one of which was discovered in the Carma_CNS-transcriptome as well as an AKH/ACP-like receptor, which has similarities to both AKH and ACP receptors. The AKH receptor identified here (AKH/RPCH R Fig. 7) was subsequently de-orphanised and confirmed as the receptor of the native RPCH neuropeptide [89].
G protein-coupled receptors for ion transport peptide (ITP), which are part of the CHH superfamily (type 1 CHHs, see [19]), and ITP-like peptides have been de-orphanised in the silkworm Bombyx mori [94] In the present study two ITPR-like GPCRs were discovered (Fig. 7, Additional file 2). Since these ligand-receptor pairs are likely to be highly conserved, and given that ITPR-like GPCRs have been identified in the decapods [32, 51] and may be broadly represented across crustaceans, this is a potentially very important avenue that must be pursued further in an effort to establish a functional CHH receptor in C. maenas.
The number of putative neuropeptide GPCRs identified for C. maenas exceeds that determined from insect genomes for Drosophila (49), Tribolium castaneum (48), and Bombyx mori (39) and from insect transcriptomes for Nilaparvata lugens (57) and Chilio suppressalis (51) [93, 95,96,97,98], but is fewer than that determined from chelicerate genomes for Stegodyphus mimosarum (120) and Mesobuthus martensii (93) [64]. The relatively large number of neuropeptide GPCRs in chelicerates, reflected also in the number of peptide ligands, results from gene paralogs that have been maintained in most extant chelicerates after ancient genome duplication events [64]. Phylogenetic analysis of putative neuropeptide GPCRs supports the notion that there are paralogous receptor genes for a number of neuropeptide receptors within C. maenas (Figs. 7, 8, 9, Table 2). Four putative CCHa receptors were discovered which cluster with the Drosophila CCHa-1 and -2 receptors (CG30106 and CG14593, respectively: Fig. 7). Two neuropeptide F receptors (NPFR) and two NFPR-like receptors cluster with the Drosophila NPFR, CG1147. Finally, three putative ecdysis triggering hormone receptors (ETHR) were discovered and that cluster with Drosophila ETHR (CG5911) and T. castaneum ETHR (Fig. 7). These are examples of multiple neuropeptide gene paralogs identified in C. maenas, but there are numerous others where two paralogous receptor sequences have been found in the Carma_CNS-transcriptome but which occur as a single receptor gene in Drosophila; for instance, the natalisin receptor, SIFa receptor, sulfakinin receptor, trissin receptor, and proctolin receptor (Fig. 7). A summary of the deduced number of neuropeptides and putative neuropeptide receptor genes is detailed in Table 2. The high confidence clustering of putative neuropeptide receptors identified here for C. maenas with those de-orphanised in Drosophila (and T. castaneum etc.) is strongly suggestive of the functional identity of the putative ligands (as we have done). However, the sequence similarities between GPCRs means that such an approach has potential to be erroneous and so functional de-orphanisation strategies are essential in the identification of GPCR transcripts with their ligands [93].
Only one of the putative neuropeptide receptors discovered was differentially expressed at the transcript level across the moult cycle: proctolin R1 (Fig. 10). Two potential proctolin receptors were discovered (as is the case for other decapods crustaceans, see Table 2 [29, 51]); gene expression analysis showed that R2 was constitutively expressed across the moult cycle whilst R1 was more highly expressed during pre- and post-moult, in stages D4 and A-B, relative to early inter-moult and early pre-moult, D1 (Fig. 10). A single neuropeptide proctolin is present in the Carma_CNS-transcriptome whilst there are two putative receptors (Table 2). The presence of two paralogs, coupled with the differential expression of one across the moult cycle, suggests that the receptors may have specialized roles and that proctolin R1 may be important during the moult cycle. Proctolin is a powerful stimulant of muscle contractions. For instance in Drosophila, proctolin induced sustained muscle contractions in third instar larvae with the CNS removed whilst diminished proctolin receptor expression leads to reduced crawling velocity at high temperatures [99]. In crustaceans, proctolin has been isolated from the pericardial organs and can affect muscle contractions of heart and skeletal muscles and modulate pyloric rhythms [100,101,102]. These reports align with the necessity for stereotyped muscle contractions evident during ecdysis as the animal extricates the old cuticle, and in light of the putative proctolin receptor expression dynamics it is tempting to speculate that this peptide and receptor pair plays a role in this regard. It will be intriguing to investigate tissue specific expression of the two putative proctolin receptors.
Annotation of the Carma_CNS-transcriptome resulted in 221 clusters with the GO:0004930 ‘G protein-coupled receptor activity’. Of these 221 clusters, eight were found to be differentially expressed across the moult cycle (one of which was proctolin R1: annotated as ‘FMRFamide receptor-like’, Cluster-569.17428; Additional file 1: Table S6 and Figure S12). All except for Cluster-569.40579 (annotated as Fizzled-2-like) were elevated during post- and pre-moult relative to inter-moult (Additional file 1: Figure S12). With the exception of proctolin R1 (Cluster-569-17,428), no other GPCRs showed high sequence similarity to known neuropeptide receptors. Given that the moult cycle and the process of ecdysis pervade all levels of crustacean physiology, it is perhaps unsurprising that non-neuropeptide GPCRs are differentially expressed across the moult cycle.
Carma_YO-transcriptome
Thirty putative neuropeptide GPCRs were discovered within the Carma_YO-transcriptome: 43 transcripts (from 35 clusters) were identified encoding 24 full-length proteins, seven C-terminal partial proteins, five N-terminal partial proteins, and seven internal fragments (see Additional file 2). Of these, 25 were predicted to contain seven transmembrane helices (predicted by tmhmm). In addition to GPCRs, two venus kinase receptors were identified. Of the 30 putative neuropeptide GPCRs, 20 were rhodopsin-like, four were secretin-like, and six were LGR type and were assigned putative ligands from 24 neuropeptide families based on sequence similarity to de-orphanised neuropeptide GPCRs (see Additional file 2). Transcripts for the putative neuropeptide receptors allatostatin-A R, CRZR, CCAPR, short neuropeptide F receptor (sNPFR), and bursicon R2, encoded full-length proteins in the YO, whilst those identified from the Carma_CNS-transcriptome were incomplete. In addition, whilst a single calcitonin-like-DH-31 R was identified from the Carma_CNS-transcriptome, an additional three transcript variants were identified in the Carma_YO-transcriptome (see Additional file 2). Of the putative neuropeptide receptors identified in the YO, seven were differentially expressed across the moult cycle (Table 3, Fig. 11). sNPFR was the most differentially expressed cluster over the moult cycle within the Carma_YO-transcriptome, with low expression levels throughout most of the moult cycle, but was upregulated during early pre-moult, D1 (Table 3, Fig. 11A).
Annotation of the Carma_YO-transcriptome annotated 160 clusters with ‘GO:0004930 G protein-coupled receptor activity’. Of these, 18 were found to be differentially expressed across the moult cycle (Additional file 1: Table S9 and Figure S13), one of which was sNPFR and another was allatostatin-B/myoinhibitory peptide R1. Of the remaining 16 GPCRs, which did not have high sequence similarity to known neuropeptide receptors, 12 contained 7 transmembrane helixes (predicted by tmhmm). Differential expression of GPCRs within the YO, the location of ecdysteroid synthesis, is suggestive of a role for these GPCRs in ecydsis control.
Carma_YO_vs_Epi-transcriptome
Given that the YO is pivotal in moult control via its central role in ecdysteroid synthesis, we anticipated that elucidation of GPCRs and their expression dynamics within this tissue might yield insight to receptor-ligand pairs critical to moult regulation. Our strategy was to explore differential expression between the YO and epidermis, from which it is derived, to separate transcripts exclusive to the YO. Differential expression of clusters between YO and epidermis tissues during inter-moult (C3–4) identified 162 clusters annotated as GPCRs (‘GO: GO:0004930 G protein-coupled receptor activity’), 11 of which were upregulated in the YO relative to the epidermis (Table 4, Fig. 12) and of these, eight were predicted to contain 7 transmembrane helices (predicted by TMHMM). Three of the GPCRs upregulated in the YO relative to the epidermis were essentially absent from the epidermis, and two were also differentially expressed in the YO across the moult cycle; rhodopsin G0-coupled and probable G protein-coupled receptor Methuselah (Mth)-like 1 (Table 4, Fig. 12). The GPCR cluster that showed most marked differential expression between YO and epidermis tissues was annotated as a gonadotropin-releasing hormone receptor and, according to our phylogenetic analysis of rhodopsin-like GPCRs (Fig. 7), was a candidate for the corazonin receptor (CRZR) (Table 4, Fig. 13). Corazonin initiates the release of ETH from ‘Inka’ cells in M. sexta, so the identification of a putative CRZR on the YO was exciting. Analysis of the Carma_YO-transcriptome revealed that this receptor was not differentially expressed over the moult cycle (Table 4, Fig. 13B). Differential expression of CRZR between YO and epidermis was validated by qPCR: two-way ANOVA revealed no interaction effect (tissue x moult stage), but main effects revealed a significant difference in CRZR mRNA expression between tissues (being greater in the YO; F(1,46) = 63.3, P < 0.001) but not across moult stage, F(4, 46) = 1.01, P = 0.412 (Fig. 13C, D).
We also performed qPCR on CNS tissues to verify expression patterns of both the peptide CRZ and the putative CRZR (Fig. 14). In silico data for CRZ revealed no differential expression across the moult cycle. Two way ANOVA of qPCR data however, showed that whilst there was no interaction between moult cycle and different CNS tissues, main effects were significant for tissue (F(2,71) = 482.77, P < 0.001) and moult stage (F(4,71) = 4.27, P = 0.004). Pairwise comparisons with Bonferroni correction post hoc analysis indicated that moult stage D1 differed from all other stages (P < 0.05), mRNA abundances being higher during this stage. For tissues, CRZ mRNA abundance was significantly greater within the ESG relative to CG and VG (P < 0.05). For the putative CRZR within the CNS, in silico data showed no differential gene expression across the moult cycle. Two-way ANOVA of CRZR qPCR data revealed no interaction between moult stage and tissue, but main effects were significant for tissue (F(2,71) = 4.75, P = 0.0116); Bonferroni post-hoc analysis showed that VG CRZR mRNA expression was greater than ESG and CG (P < 0.05; Fig. 14D). Expression values for CRZR determined via qPCR show that expression levels within all regions of the CNS were low relative to expression within the YO (Figs. 13D, 14D). This putative CRZR has subsequently been de-orphanised as the receptor to the native CRZ neuropeptide [31].
CRZ qPCR identified the primary location of peptide synthesis within the CNS to be the eyestalk ganglion (ESG, Fig. 14), whilst CRZR was expressed at far higher values in the YO relative to all parts of the CNS. Intriguingly, given the crucial role of the YO to ecdysteroidogenesis and therefore ecdysis, high values for CRZR expression in the YO clearly indicate a hitherto unknown and unexpected role for this ligand-receptor pair in ecdysis regulation.
Conclusions
The Carma_CNS-, Carma_YO-, and Carma_YO_vs_Epi- transcriptomes presented here provide a high quality resource for transcript discovery in relation to crustacean ecdysis regulation. The Carma_CNS-transcriptome represents the most complete crustacean transcriptome published to date in terms of neuropeptide and neuropeptide receptor sequence diversity and completeness. Moreover, the temporal element of this study has enabled the first comprehensive exploration of putative crustacean neuropeptide signaling across the moult cycle and across multiple tissues. Within the CNS of C. maenas, the neuropeptides carcikinin/ETH, CHH-2, and CCAP were differentially expressed across the moult cycle and are, therefore, anticipated to play important roles in the regulation of ecdysis. ETH and CCAP are fundamental to the temporal regulation of ecdysis in insects and release of CCAP during crustacean ecdysis is concomitant with the onset of ecdysis motor programme. The probable roles for both carcikinin/ETH and CCAP in crustacean ecdysis hint at conserved commonalities between insect and crustacean ecdysis cassettes. CHH-2, previously not studied, is now a neuropeptide of interest with respect to ecdysis physiology. The differing gene expression dynamics of CHH-2 relative to CHH-1, MIHs, and CHH-MIH-like peptides found here suggests that these CHHs may fulfil different physiological functions. However, since the translation and processing of functional peptides may be temporally uncoupled from mRNA expression, conclusions peptides availability and bioactivity of these cannot be inferred from mRNA abundance without due caution. A putative proctolin R1 was differentially expressed across the moult cycle, which marks a novel association of proctolin with ecdysis, but no other putative neuropeptide GPCR was found to be differentially expressed across the moult cycle in the CNS of C. maenas.
To our surprise, neuropeptide transcripts were discovered within the YO and amongst these neuroparsin-1 was found to be differentially expressed across the moult cycle, highlighting the YO as a source of neuropeptide synthesis. A number of putative neuropeptide receptors were identified as being differentially expressed across the moult cycle within the YO and, given that the only known function of the YO is synthesis of ecdysteroids, are putative candidates in ecdysis regulation.
Our analysis revealed relatively few GPCRs specific to the YO (compared to the epidermis). The most differentially expressed of these were essentially absent from epidermis and included a receptor annotated as gonadotropin-releasing hormone receptor which, considering phylogenetic grouping was revealed as a candidate corazonin receptor. Tantalisingly, this putative CRZR was strongly expressed in the YO relative to the epidermis, and was lowly expressed in the CNS, suggesting an unexpected and potentially very important role for CRZ-CRZR pair in ecdysis regulation.
The comprehensive list of putative GPCRs discovered here provides a basis from which receptor de-orphanisation may advance, which is absolutely necessary to enable neuropeptide and receptor pair function to be defined.
Abbreviations
- aa:
-
Amino acid
- ACP:
-
Adipokinetic hormone/corazonin-related peptide
- AG:
-
Abdominal ganglia
- AKH:
-
Adipokinetic hormone
- CCAP:
-
Crustacean cardioactive peptide
- cDNA:
-
Complementary deoxyribonucleic acid
- CFSH:
-
Crustacean female sex hormone
- CG:
-
Cerebral ganglia
- CHH:
-
Crustacean hyperglycemic hormone
- CNS:
-
Central nervous system
- cRNA:
-
Complementary ribonucleic acid
- CRZ:
-
Corazonin
- CRZR:
-
Corazonin receptor
- DEPC:
-
Diethyl pyrocarbonate
- DH:
-
Diuretic hormone
- EF1a:
-
Elongation factor 1-alpha
- EH:
-
Eclosion hormone
- ESG:
-
Eyestalk ganglia
- ETH:
-
Ecdysis triggering hormone
- ETHR:
-
Ecdysis triggering hormone receptor
- GPCR:
-
G protein-coupled receptor
- IAG:
-
Androgenic insulin-like peptide
- LGR:
-
Leucine-rich repeat-containing
- MIH:
-
Moult-inhibiting hormone
- MOIH:
-
Mandibular organ inhibitory hormone
- mRNA:
-
Messenger ribonucleic acid
- Mth:
-
Methuselah
- NGS:
-
Next generation sequencing
- NPFR:
-
Neuropeptide F receptor
- OEH:
-
Ovary ecdysteroidogenic hormone
- OMP:
-
Ovary maturing parsin
- ORF:
-
Open reading frame
- PAGE:
-
Polyacrylamide gel electrophoresis
- PRP:
-
Precursor related peptide
- PTTH:
-
Prothoracicotropic hormone
- RNA:
-
Ribonucleic acid
- RPCH:
-
Red pigment concentrating hormone
- RT:
-
Room temperature
- RT-PCR:
-
Reverse transcription polymerase chain reaction
- sNPFR:
-
Short neuropeptide F receptor
- SOG:
-
Sub-oesophageal ganglia
- TG:
-
Thoracic ganglia
- UBE2L3:
-
Ubiquitin-conjugating enzyme E2 L3
- VG:
-
Ventral ganglia
- VIH:
-
Vitellogenesis-inhibiting hormone
- Vm:
-
Ventro-medial
- YO:
-
Y organ
References
White BH, Ewer J. Neural and hormonal control of postecdysial behaviors in insects. Annu Rev Entomol. 2014;59:363–81.
Žitňan D, Adams M. Neuroendocrine regulation of Ecdysis. In: Gilbert LI, editor. Insect Endocrinol. London: Elsevier Academic Press; 2012. p. 253–301.
Kim Y-J, Žitňan D, Cho K-H, Schooley DA, Mizoguchi A, Adams ME. Central peptidergic ensembles associated with organization of an innate behavior. Proc Natl Acad Sci U S A. 2006;103:14211–6.
Kim Y-J, Žitňan D, Galizia CG, Cho K-H, Adams ME. A command chemical triggers an innate behavior by sequential activation of multiple peptidergic ensembles. Curr Biol. 2006;16:1395–407.
Diao F, Mena W, Shi J, Park D, Diao F, Taghert P, et al. The splice isoforms of the Drosophila ecdysis triggering hormone receptor have developmentally distinct roles. Genetics. 2016;202:175–89.
Kingan TG, Gray W, Žitňan D, Adams ME. Regulation of ecdysis-triggering hormone release by eclosion hormone. J Exp Biol. 1997;200:3245–56.
Žitňan D, Hollar L, Spalovská I, Takáč P, Žitňanová I, Gill SS, et al. Molecular cloning and function of ecdysis-triggering hormones in the silkworm Bombyx mori. J Exp Biol. 2002;205:3459–73.
Kim Y-J, Spalovská-Valachová I, Cho K-H, Žitňanová I, Park Y, Adams ME, et al. Corazonin receptor signaling in ecdysis initiation. Proc Natl Acad Sci U S A. 2004;101:6704–9.
Li B, Predel R, Neupert S, Hauser F, Tanaka Y, Cazzamali G, et al. Genomics, transcriptomics, and peptidomics of neuropeptides and protein hormones in the red flour beetle Tribolium castaneum. Genome Res. 2008;18:113–22.
Gammie SC, Truman JW. Neuropeptide hierarchies and the activation of sequential motor behaviors in the hawkmoth, Manduca sexta. J Neurosci. 1997;17:4389–97.
Gammie SC, Truman JW. Eclosion hormone provides a link between ecdysis-triggering hormone and crustacean cardioactive peptide in the neuroendocrine cascade that controls ecdysis behavior. J Exp Biol. 1999;202:343–52.
Taghert PH, Truman JW. Identification of the bursicon-containing neurones in abdominal ganglia of the tobacco hornworm, Manduca sexta. J Exp Biol. 1982;98:385–401.
Kostron B, Kaltenhauser U, Seibel B, Bräunig P, Honegger HW. Localization of bursicon in CCAP-immunoreactive cells in the thoracic ganglia of the cricket Gryllus bimaculatus. J Exp Biol. 1996;199:367–77.
Dewey EM, McNabb SL, Ewer J, Kuo GR, Takanishi CL, Truman JW, et al. Identification of the gene encoding bursicon, an insect neuropeptide responsible for cuticle sclerotization and wing spreading. Curr Biol. 2004;14:1208–13.
Regier JC, Shultz JW, Zwick A, Hussey A, Ball B, Wetzer R, et al. Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences. Nature. 2010;463:1079–83.
Averof M, Akam M. Insect-crustacean relationships: insights from comparative developmental and molecular studies. Philos Trans R Soc London B. 1995;347:293–303.
Lafont R, Dauphin-Villemant C. Ecdysteroid chemistry and biochemistry. In: Gilbert LI, editor. Insect Endocrinol. London: Elsevier Academic press; 2012. p. 106–56.
Webster SG. Endocrinology of molting. In: Chang ES, Thiel M, editors. Nat Hist Crustac vol 4 Physiol Regul. Oxford: Oxford University Press; 2014. p. 1–67.
Webster SG, Keller R, Dircksen H. The CHH-superfamily of multifunctional peptide hormones controlling crustacean metabolism, osmoregulation, moulting, and reproduction. Gen Comp Endocrinol. 2012;175:217–33.
Chung JS, Webster SG. Dynamics of in vivo release of molt-inhibiting hormone and crustacean hyperglycemic hormone in the shore crab, Carcinus maenas. Endocrinology. 2005;146:5545–51.
Passano LM. Neurosecretory control of molting in crabs by the X-organ sinus gland complex. Physiol Comp Oecol. 1953;3:155–89.
Webster SG, Wilcockson DC, Mrinalini, Sharp JH. Bursicon and neuropeptide cascades during the ecdysis program of the shore crab, Carcinus maenas. Gen Comp Endocrinol. 2013;182:54–64.
Webster SG, Keller R. Purification, characterisation and amino acid composition of the putative moult-inhibiting hormone (MIH) of Carcinus maenas (Crustacea, Decapoda). J Comp Physiol B. 1986;156:617–24.
Fernlund P, Josefsson L. Crustacean color-change hormone: amino acid sequence and chemical synthesis. Science. 1972;177:173–5.
Christie AE, Stemmler EA, Dickinson PS. Crustacean neuropeptides. Cell Mol Life Sci. 2010;67:4135–69.
Christie AE. Expansion of the neuropeptidome of the globally invasive marine crab Carcinus maenas. Gen Comp Endocrinol. 2016;235:150–69.
Veenstra JA. Similarities between decapod and insect neuropeptidomes. PeerJ. 2016;4:e2043.
Christie AE. Expansion of the Litopenaeus vannamei and Penaeus monodon peptidomes using transcriptome shotgun assembly sequence data. Gen Comp Endocrinol. 2014;206:235–54.
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.
Marco HG, Verlinden H, Vanden Broeck J, Gäde G. Characterisation and pharmacological analysis of a crustacean G protein-coupled receptor: the red pigment-concentrating hormone receptor of Daphnia pulex. Sci Rep. 2017;7:6851.
Alexander JL, Oliphant A, Wilcockson DC, Audsley N, Down RE, Lafont R, et al. Functional characterization and signaling systems of corazonin and red pigment concentrating hormone in the green shore crab, Carcinus maenas. Front Neurosci. 2018;11:752. https://doi.org/10.3389/fnins.2017.00752.
Buckley SJ, Fitzgibbon GP, 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.
Phlippen MK, Webster SG, Chung JS, Dircksen H. Ecdysis of decapod crustaceans is associated with a dramatic release of crustacean cardioactive peptide into the hemolymph. J Exp Biol. 2000;203:521–36.
Saver MA, Wilkens JL, Syed NI. In situ and in vitro identification and characterization of cardiac ganglion neurons in the crab, Carcinus maenas. J Neurophysiol. 1999;81:2964–76.
Yang Y, Smith SA. Optimizing de novo assembly of short-read RNA-seq data for phylogenomics. BMC Genomics. 2013;14:328.
Wang S, Gribskov M. Comprehensive evaluation of de novo transcriptome assembly programs and their effects on differential gene expression analysis. Bioinformatics. 2017;33:327–33.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29:644–52.
Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012;28:1086–92.
Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.
Peng Y, Leung HCM, Yiu S-M, Lv M-J, Zhu X-G, Chin FYL. IDBA-Tran: a more robust de novo de Bruijn graph assembler for transcriptomes with uneven expression levels. Bioinformatics. 2013;29:i326–34.
Liu J, Li G, Chang Z, Yu T, Liu B, McMullen R, et al. BinPacker: packing-based de novo transcriptome assembly from RNA-seq data. PLoS Comput Biol. 2016;e1004772:12.
Hansen KD, Brenner SE, Dudoit S. Biases in Illumina transcriptome sequencing caused by random hexamer priming. Nucleic Acids Res. 2010;38:1–7.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
MacManes MD. On the optimal trimming of high-throughput mRNA sequence data. Front Genet. 2014;5:1–7.
Davidson NM, Oshlack A. Corset: enabling differential gene expression analysis for de novo assembled transcriptomes. Genome Biol. 2014;15:410.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acid Symp Ser. 1999;41:95–8.
National Center for Biotechnology Information, U.S. National Library of Medicine 8600 Rockville Pike, Bethesda MD, 20894 USA.
Chipman AD, Ferrier DEK, Brena C, Qu J, Hughes DST, Schröder R, et al. The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima. PLoS Biol. 2014;e1002005:12.
Veenstra JA. The power of next-generation sequencing as illustrated by the neuropeptidome of the crayfish Procambarus clarkii. Gen Comp Endocrinol. 2015;224:84–95.
Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012;40:W597–603.
Nielsen H. Predicting secretory proteins with SignalP. In: Kihara D, editor. Protein function prediction, methods in molecular biology, vol. 1611. New York: Humana Press; 2017. p. 59–73.
Veenstra JA. Mono- and dibasic proteolytic cleavage sites in insect neuroendocrine peptide precursors. Arch Insect Biochem Physiol. 2000;43:49–63.
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.
Sonnhammer ELL, Von Heijne G, Krogh A. A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol. 1998;6:175–82.
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.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Sharp JH, Wilcockson DC, Webster SG. Identification and expression of mRNAs encoding bursicon in the plesiomorphic central nervous system of Homarus gammarus. Gen Comp Endocrinol. 2010;169:65–74.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611–22.
Dircksen H, Neupert S, Predel R, Verleyen P, Huybrechts J, Strauss J, et al. Genomics, transcriptomics, and peptidomics of Daphnia pulex neuropeptides and protein hormones. J Proteome Res. 2011;10:4478–504.
Nässel DR, Winther ÅME. Drosophila neuropeptides in regulation of physiology and behavior. Prog Neurobiol. 2010;92:42–104.
Veenstra JA. Neuropeptide evolution: Chelicerate neurohormone and neuropeptide genes may reflect one or more whole genome duplications. Gen Comp Endocrinol. 2016;229:41–55.
Ma M, Bors EK, Dickinson ES, Kwiatkowski MA, Sousa GL, Henry RP, et al. Characterization of the Carcinus maenas neuropeptidome by mass spectrometry and functional genomics. Gen Comp Endocrinol. 2009;161:320–34.
Larsen MJ, Burton KJ, Zantello MR, Smith VG, Lowery DL, Kubiak TM. Type a allatostatins from Drosophila melanogaster and Diplotera puncata activate two Drosophila allatostatin receptors, DAR-1 and DAR-2, expressed in CHO cells. Biochem Biophys Res Commun. 2001;286:895–901.
Christie AE. Prediction of the peptidomes of Tigriopus californicus and Lepeophtheirus salmonis (Copepoda, Crustacea). Gen Comp Endocrinol. 2014;201:87–106.
Ventura T, Rosen O, Sagi A. From the discovery of the crustacean androgenic gland to the insulin-like hormone in six decades. Gen Comp Endocrinol. 2011;173:381–8.
Roller L, Žitňanová I, Dai L, Šimo L, Park Y, Satake H, et al. Ecdysis triggering hormone signaling in arthropods. Peptides. 2010;31:429–41.
Krüger E, Mena W, Lahr EC, Johnson EC, Ewer J. Genetic analysis of eclosion hormone action during Drosophila larval ecdysis. Development. 2015;142:4279–87.
Chung JS, Wilcockson DC, Zmora N, Zohar Y, Dircksen H, Webster SG. Identification and developmental expression of mRNAs encoding crustacean cardioactive peptide (CCAP) in decapod crustaceans. J Exp Biol. 2006;209:3862–72.
Wilcockson DC, Webster SG. Identification and developmental expression of mRNAs encoding putative insect cuticle hardening hormone, bursicon in the green shore crab Carcinus maenas. Gen Comp Endocrinol. 2008;156:113–25.
Kim DH, Han MR, Lee G, Lee SS, Kim YJ, Adams ME. Rescheduling behavioral subunits of a fixed action pattern by genetic manipulation of peptidergic signaling. PLoS Genet. 2015;11:1–25.
Lacombe C, Grève P, Martin G. Overview on the sub-grouping of the crustacean hyperglycemic hormone family. Neuropeptides. 1999;33:71–80.
Chung JS, Zmora N, Katayama H, Tsutsui N. Crustacean hyperglycemic hormone (CHH) neuropeptides family: functions, titer, and binding to target tissues. Gen Comp Endocrinol. 2010;166:447–54.
Kegel G, Reichwein B, Weese S, Gaus G, Peter-Katalinić J, Keller R. Amino acid sequence of the crustacean hyperglycemic hormone (CHH) from the shore crab, Carcinus maenas. FEBS J. 1989;255(1):10–4.
Dircksen H, Böcking D, Heyn U, Mandel C, Chung JS, Baggerman G, et al. Crustacean hyperglycaemic hormone (CHH)-like peptides and CHH-precursor-related peptides from pericardial organ neurosecretory cells in the shore crab, Carcinus maenas, are putatively spliced and modified products of multiple genes. Biochem J. 2001;356:159–70.
Chung JS, Webster SG. Does the N-terminal pyroglutamate residue have any physiological significance for crab hyperglycemic neuropeptides? Eur J Biochem. 1996;240:358–64.
Chung JS, Dircksen H, Webster SG. A remarkable, precisely timed release of hyperglycemic hormone from endocrine cells in the gut is associated with ecdysis in the crab Carcinus maenas. Proc Natl Acad Sci U S A. 1999;96:13103–7.
Yang W-J, Aida K, Nagasawa H. Amino acid sequences and activities of multiple hyperglycemic hormones from the kuruma prawn, Penaeus japonicus. Peptides. 1997;18:479–85.
Zmora N, Chung JS. A novel hormone is required for the development of reproductive phenotypes in adult female crabs. Endocrinology. 2014;155:230–9.
Girardie J, Boureme D, Couillaud F, Tamarelle M, Girardie A. Anti-juvenile effect of neuroparsin A, a neuroprotein isolated from locust corpora cardiaca. Insect Biochem. 1987;17:977–83.
Girardie J, Girardie A, Huet JC, Pernollet JC. Amino acid sequence of locust neuroparsins. FEBS Lett. 1989;245:4–8.
Badisco L, Claeys I, Van Loy T, Van Hiel M, Franssens V, Simonet G, et al. Neuroparsins, a family of conserved arthropod neuropeptides. Gen Comp Endocrinol. 2007;153:64–71.
Badisco L, Marchal E, Van Wielendaele P, Verlinden H, Vleugels R, Van den Broeck J. RNA interference of insulin-related peptide and neuroparsins affects vitellogenesis in the desert locust Schistocerca gregaria. Peptides. 2011;32:573–80.
Girardie J, Girardie A. Lom OMP, a putative ecdysiotropic factor for the ovary in Locusta migratoria. J Insect Physiol. 1996;42:215–21.
Girardie J, Richard O, Girardie A. Detection of vitellogenin in the hemolymph of larval female locusts (Locusta migratoria) treated with the neurohormone, Lom OMP. J Insect Physiol. 1996;42:107–13.
Brown MR, Graf R, Swiderek KM, Fendley D, Stracker TH, Champagne DE, et al. Identification of a steroidogenic neurohormone in female mosquitoes. J Biol Chem. 1998;273:3967–71.
Yang SP, He J-G, Sun CB, Chan SF. Characterization of the shrimp neuroparsin (MeNPLP): RNAi silencing resulted in inhibition of vitellogenesis. FEBS Open Bio. 2014;4:976–86.
Sharabi O, Manor S, Weil E, Aflalo D, Lezer Y, Levy T, Aizen J, Ventura T, Khalaila I, Sagi A. Identification and characterization of an insulin-like receptor involved in crustacean reproduction. Endocrinology. 2016;157(2):928–41.
Chandler JC, Gandhi NS, Mancera RL, Smith G, Elizur A, Ventura T. Understanding insulin endocrinology in decapod crustacea: molecular modelling characterization of an insulin-binding protein and insulin-like peptides in the eastern spiny lobster, Sagmariasus verreauxi. Int J Mol Si. 2017;18:1832. https://doi.org/10.3390/ijms18091832.
Vogel KJ, Brown MR, Strand MR. Ovary ecdysteroidogenic hormone requires a receptor tyrosine kinase to activate egg formation in the mosquito Aedes aegypti. Proc Natl Acad Sci U S A. 2015;112:5057–62.
Caers J, Verlinden H, Zels S, Vandersmissen HP, Vuerinckx K, Schoofs L. More than two decades of research on insect neuropeptide GPCRs: an overview. Front Endocrinol. 2012;3:151.
Nagai C, Mabashi-Asazuma H, Nagasaw H, Nagat S. Identification and characterization of receptors for ion transport peptide (ITP) and ITP-like (ITPL) in the silkworm Bombyx mori. J Biol Chem. 2014;289(46):32166–77.
Hauser F, Cazzamali G, Williamson M, Park Y, Li B, Tanaka Y, et al. A genome-wide inventory of neurohormone GPCRs in the red flour beetle Tribolium castaneum. Front Neuroendocrinol. 2008;29:142–65.
Yamanaka N, Yamamoto S, Žitňan D, Watanabe K, Kawada T, Satake H, et al. Neuropeptide receptor transcriptome reveals unidentified neuroendocrine pathways. PLoS One. 2008;3:e3048.
Tanaka Y, Suetsugu Y, Yamamoto K, Noda H, Shinoda T. Transcriptome analysis of neuropeptides and G-protein coupled receptors (GPCRs) for neuropeptides in the brown planthopper Nilaparvata lugens. Peptides. 2014;53:125–33.
Xu G, Gu G-X, Teng Z-W, Wu S-F, Huang J, Song Q-S, et al. Identification and expression profiles of neuropeptides and their G protein-coupled receptors in the rice stem borer Chilo suppressalis. Sci Rep. 2016;6:28976.
Ormerod KG, LePine OK, Bhutta MS, Jung J, Tattersall GJ, Mercier AJ. Characterizing the physiological and behavioral roles of proctolin in Drosophila melanogaster. J Neurophysiol. 2016;115:568–80.
Hooper SL, Marder E. Modulation of a central pattern generator by two neuropeptides, proctolin and FMRFamide. Brain Res. 1984;305:186–91.
McGaw IJ, Wilkens JL, McMahon BR, Airriess CN. Crustacean cardioexcitatory peptides may inhibit the heart in vivo. J Exp Biol. 1995;198:2547–50.
Hooper S, Marder E. Modulation of the lobster pyloric rhythm by the peptide proctolin. J Neurosci. 1987;7:2097–112.
Acknowledgements
We thank Maximilian Friedersdorff for his help in making the transcriptomes assembled as part of this project available as publically accessible BLAST databases at: sequenceserver.ibers.aber.ac.uk.
Funding
This project was funded by associated BBRSC research grants BB/L021242/1 and BB/L021552/1 awarded to DCW and SGW, respectively. BBSRC played no part in study design, collection nor analysis and interpretation of data nor writing of the manuscript.
Availability of data and materials
The three transcriptomes assembled as part of this Transcriptome Shotgun Assembly project have been deposited at DDBJ/EMBL/GenBank under the accessions: GFXF00000000 (Carma_CNS-transcriptome), GFYV00000000 (Carma_YO-transcriptome), and GFYW00000000 (Carma_YO_vs_Epi-transcriptome). Raw read sequence files used to generate these assemblies are archived in the NCBI SRA archive under BioProject PRJNA400568.
Author information
Authors and Affiliations
Contributions
Sample collections were done primarily by SGW and JLA with help from AO and DCW. RNA extractions were done by AO with library preparation and sequencing done in-house at the Translational Genomics Facility, IBERS, Aberystwyth University. AO conducted the bioinformatics analysis with assistance from MS. DCW developed and conducted qRT-PCR assays with help from AO. AO and DCW wrote the initial manuscript. All authors contributed to the design of the study and to the final manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This work involved animals (crustaceans). Crustaceans are not covered by the UK Home Office Animal Scientific Procedures Act 1986 or Directive 2010/63/EU, but we adhered to our local ethical guidelines, minimized the number of animals used as far as possible, and ice-anaesthetised animals prior to dissections. Carcinus maenas were collected from Menai Strait, UK by hand and using crab pots. Collection of Carcinus maenas from this location is not restricted or subject to local or national licensing or permissions.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
Figures S1-S13 and Tables S1-S9. (DOCX 3567 kb)
Additional file 2:
Trinotate annotation reports and mined sequences for neuropeptides and G protein-coupled receptors. (XLSX 82011 kb)
Additional file 3:
Neuropeptide prepro-hormone processing and deduced mature neuropeptide and precursor-related peptides. (DOCX 88 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Oliphant, A., Alexander, J.L., Swain, M.T. et al. Transcriptomic analysis of crustacean neuropeptide signaling during the moult cycle in the green shore crab, Carcinus maenas. BMC Genomics 19, 711 (2018). https://doi.org/10.1186/s12864-018-5057-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-018-5057-3