Identification of the neuropeptide precursor genes potentially involved in the larval settlement in the Echiuran worm Urechis unicinctus
BMC Genomics volume 21, Article number: 892 (2020)
In marine invertebrate life cycles, which often consist of planktonic larval and benthonic adult stages, settlement of the free-swimming larva to the sea floor in response to environmental cues is a key life cycle transition. Settlement is regulated by a specialized sensory–neurosecretory system, the larval apical organ. The neuroendocrine mechanisms through which the apical organ transduces environmental cues into behavioral responses during settlement are not fully understood yet.
In this study, a total of 54 neuropeptide precursors (pNPs) were identified in the Urechis unicinctus larva and adult transcriptome databases using local BLAST and NpSearch prediction, of which 10 pNPs belonging to the ancient eumetazoa, 24 pNPs belonging to the ancient bilaterian, 3 pNPs belonging to the ancient protostome, 9 pNPs exclusive in lophotrochozoa, 3 pNPs exclusive in annelid, and 5 pNPs only found in U. unicinctus. Furthermore, four pNPs (MIP, FRWamide, FxFamide and FILamide) which may be associated with the settlement and metamorphosis of U. unicinctus larvae were analysed by qRT-PCR. Whole-mount in situ hybridization results showed that all the four pNPs were expressed in the region of the apical organ of the larva, and the positive signals were also detected in the ciliary band and abdomen chaetae. We speculated that these pNPs may regulate the movement of larval cilia and chaeta by sensing external attachment signals.
This study represents the first comprehensive identification of neuropeptides in Echiura, and would contribute to a complete understanding on the roles of various neuropeptides in larval settlement of most marine benthonic invertebrates.
Most marine benthic invertebrates have planktonic larvae during their life cycle. After going through a pelagic period, these planktonic larvae settle to the bottom and metamorphose into benthonic individuals (crawling, attaching, fixing and burrowing) [1, 2]. The larval settlement is the key event for their development and survival, which commonly includes the cessation of swimming and the appearance of substrate exploratory behavior [3,4,5,6]. This is a complex process determined by the interaction of biotic and abiotic factors at different temporal and spatial scale [7, 8]. The apical organ, a cluster of sensory neurons at the anterior of the larva in diverse groups as phoronids, polychaetes and chitons, has been implicated to be the site of perception cues for settlement and metamorphosis [9, 10]. Researchers found that neuropeptides expressed in chemosensory–neurosecretory cells of the apical organ can innervate ciliary bands, and suggested that they may play a role in the regulation of larval locomotion [11,12,13,14], which contribute to the larval settlement behavior .
Neuropeptides are considered to be the oldest neuronal signaling molecules in metazoans , and participate in the control of neural circuits and physiology [16,17,18]. They are generated from inactive precursor proteins by proteolytic cleavage and further modification such as C-terminal alpha-amidation and N-terminal pyroglutamination [19, 20], and then released into the hemolymph as hormones or the synapses as nerotransmitters to regulate the physiological activities of target cells . Studies on the marine invertebrate larval neuropeptide systems have mainly been focused on lophotrochozoan including Mollusca and Annelida. For example, in mollusca, they mainly focus on larval development, larval feeding behavior, larval muscle innervation and muscular contractions [22,23,24,25]. In annelid Platynereis, neuropeptides have been indicated to involve in the ciliary beating, some neuropeptides (RYa, FVMa, DLa, FMRFa, FVa, LYa, YFa SPY and L11) for the larval upward swimming and others (FLa and WLD) for the downward swimming . Furthermore, MIPs (myoinhibitory peptide) have been experimentally verified to play a role in regulating the larval settlement of marine annelid . So far, other neuropeptides related to the larval settlement remain to be explored.
Urechis unicinctus is a representative species in Echirua inhabiting the U-shaped burrows in the coastal mud flats, and is also a commercial echiuran worm in China, Japan and Korea. The worm has a typical free-swimming trochophore larva beginning with the early trochophore stage (ET, 2 days post-fertilization; dpf) and the planktonic larva settles to the bottom during the segmentation larva stage (SL, 35 dpf; also called competent larva, CL), and then burrows the sediment and metamorphoses into the benthic worm (worm-shaped larva, WL, 42 dpf). Previous studies indicated that the SL stage larvae will delay metamorphosis and their mortality rate will increase if they do not find the adaptive substrate [3, 26, 27]. In this study, to provide a basic profile of the neuropeptide precursors for investigating the role of neuropeptides in U. unicinctus larval settlement, we screened the neuropeptide precursors potentially involved in the larval settlement from the U. unicinctus larval and adult transcriptomes. Furthermore, expression characteristics of the candidate genes were validated by qRT-PCR and whole-mount in situ hybridization. To map the candidate genes to the nerve cells at the special sites, nervous system in U. unicinctus larvae was analyzed using the fluorescence immunohistochemistry. The aim of this study was to identify neuropeptide precursors potentially involved in the larval settlement in the U. unicinctus and to provide new insights in larval settlement of marine benthic invertebrates.
Results and discussion
Overview of the neuropeptide precursors in U. unicinctus
We performed BLAST search and NpSearch prediction to screen the neuropeptide precursors in the transcriptomes of U. unicinctus. A total of 54 neuropeptide precursors (pNPs) were identified, 7 from BLAST search, 5 from NpSearch prediction, and 42 from both methodologies (Fig. 1a and Supplementary Table S2). Among them, 49 pNPs had been reported in other species, and the remaining 5 pNPs were first identified in U. unicinctus and we named them FxFamide, FILamide, FW, FRWamide and ASYY according to their conserved amino acid residues. In the U. unicinctus transcriptomes, most neuropeptide precursor sequences contained the full-length open reading frame (ORF) with a signal peptide (SP), except pedal peptide 1 and FVRIamide. The sequence characteristics of U. unicinctus neuropeptide precursors for the SP presence, the conserved peptide motifs and other hallmarks of bioactive peptides, e.g. amidation C-terminal Gly, Cys-containing stretches, mono- or dibasic cleavage sites were summarized in Fig. 1a and Supplementary Fig. S1.
Due to the inherent difficulties of analyzing highly diverse and repetitive pNPs, the relationships among different families are often elusive. Therefore, Jékely  and Conzelmann , using similarity-based clustering and sensitive similarity searches, obtained a global view of metazoan pNP diversity and evolution based on a curated dataset of 6225 pNPs from 10 phyla. This approach was also useful for analyzing the phylogenetic distribution of U. unicinctus pNPs and we classified the pNP families using the same methodology. The results showed that ten pNPs in U. unicinctus were categorized as the ancient eumetazoan families (Fig. 1a), which are the repertoire neuropeptides with the short amidated peptides, such as R [F/Y] amide, Wamide, insulin-related peptide and the glycoprotein hormones [28, 31, 32]. Then, twenty-four pNPs in U. unicinctus were categorized as the ancient bilaterian families (Fig. 1a), which belong to 17 neuropeptide families . Three members of the ancient protostome neuropeptide precursor families were present in U. unicinctus, including prohormone-3, myomodulin, and whitnin-2 (Fig. 1a). Moreover, we identified nine pNPs in U. unicinctus (Fig. 1a), which were proposed to be the lophotrochozoan-specific families . Three pNPs in U. unicinctus had recognizable orthologs only in annelids, including DLamide, SLRFamide and QERAS (Fig. 1a). In addition, five pNPs did not have recognizable orthologs outside Urechis, and were temporarily classified as neuropeptides unique to U. unicinctus, including the FxFamide, FILamide, FW, FRWamide and ASYY (Fig. 1a).
Traditionally Echiura was ranked as a phylum, but recent studies, especially on molecular phylogenetic analysis  and morphological observation [34, 35], have generated an increasing body of evidence that they actually are derived annelids and provide strong support for a sister group relationship between Echiura and Capitellidae. This is consistent with our study in which we find three Annelid-specific pNPs were presented in U. unicinctus.
Screen of the neuropeptide precursors potentially involving in the larval settlement
We performed a hierarchical clustering of the neuropeptide precursors based on their stage-specific expression (FPKM values) from the U. unicinctus transcriptomes (Fig. 1b and Supplementary Fig. S2). The results showed that most of the neuropeptide precursors were expressed at multiple stages, and the expression levels were significantly different. We found when the larvae developed from LT to ES, a process which the larvae initially transited from upper to middle layer in the water column, expression levels of NPF-3 and DH44 decreased significantly (p < 0.05), while that of bursicon-A2 and NPF-1 was significantly increased (p < 0.05) (Fig. 1b and Supplementary Fig. S2). During the development progress from ES to SL, a period that the larvae move from the middle to the bottom of water layer, eleven pNP genes (MIP, bursicon-A2, NPF-2, RGWamide, 7B2, pedal peptide 1, myomodulin, FVRIamide, FxFamide, FILamide and FRWamide) were significantly up-regulated (p < 0.05) (Fig. 1b and Supplementary Fig. S2). However, eight genes (except 7B2, pedal peptide 1 and FVRIamide) among the eleven pNPs above were again down-regulated (p < 0.05) when the larvae developed from SL to WL, which is the period that the larvae begin to explore the suitable substrate and finally became benthic larvae (Fig. 1b and Supplementary Fig. S2). As MIP have been confirmed to regulate larvae settlement behavior , we speculated preliminarily the eight pNPs with similar expression pattern were considered to be most likely pNPs involved in the regulation of larval settlement and metamorphosis in U. unicinctus.
Sequence characteristics of the selected pNPs that may be involved in the regulation of larval settlement in U. unicinctus
Four interesting pNPs, including previously reported MIP  and three Uu-specific pNPs (FxFamide, FILamide and FRWamide), were selected for further analysis.
MIPs (Myoinhibitory peptides) are pleiotropic neuropeptides first described in insects as inhibitors of muscle contractions [18, 36, 37]. In some insect species, MIPs modulate juvenile hormone synthesis and reduce food intake, and they are also referred to as allatostatin-B or WWamide [38,39,40,41]. In Platynereis the MIPs have been confirmed to regulate larvae settlement behavior  and feeding behavior . They are characterized by a conserved domain containing two Trp residues which are usually separated by five to eight amino acid residues in insects, molluscs and annelids [43,44,45]. In U. unicinctus transcriptomes, we identified a neuropeptide precursor which is an orthologue of arthropod MIP (Fig. 2). The Uu-MIP precursor contains 11 mature peptides, the number of mature peptides in U. unicinctus is the same as that of the annelid Platynereis dumerilii, while differs from the mollusc Patinopecten yessoensis and the arthropod Megabalanus volcano which have 10 mature peptides (Fig. 2a). Sequence alignment of the bioactive peptides revealed that the sequence similarity among the different mature MIPs in U. unicinctus was higher than that in P. dumerilii, P. yessoensis and M. volcano (Fig. 2b). Moreover, the MRVWamide motif in C-terminal of the mature MIPs is present in U. unicinctus and P. dumerilii, but not in P. yessoensis and M. volcano (Fig. 2b). The above results show that the characteristics of the MIP precursor sequence of U. unicinctus are closer to that of P. dumerilii, which is consistent with the classic species evolution.
In this study, five potential neuropeptide precursors were for the first time identified in U. unicinctus (Fig. 1a), and three of them (FRWamide, FILamide and FxFamide) were predicted to play a role in regulating U. unicinctus larvae settlement based on the significant differences in mRNA level from the segmented larvae to worm-shaped larvae (Fig. 1b). FRWamide precursor is comprised of 202 amino acids which contains a 25-residue signal peptide and 7 copies of neuropeptides with FRWamide motif in the C-terminal (Fig. 3a and d). FILamide precursor is comprised of 336 amino acids which contains a 27-residue signal peptide and 8 copies of neuropeptides with FILamide motif in the C-terminal (Fig. 3b and e). FxFamide precursor is comprised of 509 amino acids which contains a 33-residue signal peptide and 15 copies of neuropeptides with FxFamide motif in the C-terminal (Fig. 3c and f). These newly discovered neuropeptide precursors enrich the intension of neuropeptide composition.
Spatio-temporal expression of the selected pNPs during the larval settlement
To verify the expression of the four pNP transcripts (MIP, FILamide, FxFamide and FRWamide), U. unicinctus larvae including late-trochophore (LT), pre-competent larva (PL), competent larva (CL), post-competent larva (POL) and worm-shaped larva (WL) were employed for qRT-PCR analysis (Fig. 4a). The results showed that the mRNA levels of the four pNP genes increased through larval development, with the highest expression in CL, and then significant decrease in POL and WL (Fig. 4b). These results are consistent with the transcriptome data (Fig. 1b and Supplementary Fig. S2). During the developmental progression from LT to CL, the U. unicinctus larvae move from the upper to the middle layer in water, and gradually acquire the ability to explore a suitable substrate in CL, finally become benthic larvae in WL. Thus, we suggested that these four genes may be involved in the biological activities of the larvae exploring the substrate for settlement in U. unicinctus.
To map the expression of these pNPs (MIP, FxFamide, FILamide and FRWamide) to nerve cells at the special sites, nervous system in U. unicinctus larvae was analyzed using fluorescence immunohistochemistry with an anti-5HT antibody (Fig. 5). The results showed that, in trochophore up to an age of approximately 15 days, only a few structures of the nervous system are labeled with antibodies against 5-HT. In the episphere of the larvae, the circumoesophageal connectives (CC) and two nerve rings innervating the prototroch and metatroch are visible (Fig. 5a). In the hyposphere of the larva, two longitudinal nerves (LN) merge after a short distance forming a median nerve named ventral nerve cord (VNC). Two pairs of perikarya are discernible in the anterior region, directly behind the slit-shaped mouth opening (Fig. 5a, c) and on the telotroch nerve ring (Fig. 5a, b, c). In dorsal view of the larva, 3–4 LNs can be seen in the episphere which connect to the prototroch and metatroch nerve rings (Fig. 5b). As development proceeds up to the competent larva, in which the anterior chaetae have already been formed, the paired longitudinal nerve tracts of the VNC are fused in the ventral midline (Fig. 5c). In addition, the metatroch nerve ring is disappeared and two labeled perikaryas are visible on the dorsal side of the larva just under the prototroch nerve ring (Fig. 5d, e). The apical organ of the larva is shown in Fig. 5f. Fluorescence immunohistochemistry were also used to study the development of the nervous system in various other Echiuran species, such as Bonellia viridis [46, 47] and Urechis caupo . Our results are consistent with those previous studies which have proven to be informative in the study of neurogenesis in neuronal structures of Echiurans. Besides, we revealed several previously unreported details — eight nerve fibers and six large labeled perikaryas are visible in the apical organ in Echiuran worm (Fig. 5f), which is similar to that reported in P. dumerilii [9, 21, 30, 48] and especially in Capitella teleta .
Next, location of four pNP mRNAs including MIP, FxFamide, FILamide and FRWamide were detected by Whole-mount mRNA in situ hybridization (WISH) (Fig. 6 and Figure S3). The results showed that a positive MIP signal was first observed in the central region of the episphere in the early-trochophore larva (Fig. 6a) which is similar to that of the apical organ in C. teleta and P. dumerilii [9, 48]. As the development proceeds, four positive cells are exclusively located in the apical organ of the late-trochophore larva (Fig. 6b). Until the competent larvae, the obvious positive signals were located in four regions, including the apical organ (the 4–6 cells), above the abdomen chaetae (the two cells), the prototroch in the dorsal side of the larvae (the two cells) and the both side of the telotroch (the two cells) (Fig. 6c). The expression patterns of FRWamide, FxFamide and FILamide were similar to that of MIP (Fig. 6), except no visible FRWamide positive signal was observed in U. unicinctus early-trochophore (Fig. 6d), in competent larvae FxFamide was detected only in apical organ (the six stained cells) and above the abdomen chaetae (the two stained cells) (Fig. 6i), while the positive signal of FILamide in the competent larva was only in the four positive cells of apical organ (Fig. 6l). Since the WISH experiment was observed after being sealed with resin, it was difficult to observe the apical view of the larva.
In many marine invertebrates the transition from free-swimming larvae to bottom-dwelling juveniles is regulated by neuroendocrine signals (including neuropeptides and hormones) [49,50,51]. In diverse ciliated marine larvae, the apical organ, has been implicated in the detection of cues for the initiation of larval settlement [9, 10, 52,53,54,55]. Previous studies suggest that several neuropeptides expressed in distinct sensory neurons (apical organ) innervate locomotor cilia, which contribute to the larvae swimming depth [9, 21]. The alternation of active upward swimming and passive sinking, together with swimming speed and sinking rate, is thought to determine vertical distribution in the water . In Platynereis, neuropeptides including RYa, FVMa, DLa, FMRFa, FVa, LYa, YFa, FLa, MIP, GWa et al. can alter ciliary beat frequency and the rate of calcium-evoked ciliary arrests [9, 21], which eventually may be involved in regulating the larvae settlement. In our study, the MIP and FRWamide were detected in the dorsal side of the prototroch nerve ring and the both side of the telotroch nerve ring (Fig. 6c, f), indicating that they may play a role in U. unicinctus ciliary beating and eventually cause the larvae sinking to the bottom. MIP is the only neuropeptide that has been shown to be involved in larval settlement in Platynereis , which is expressed in chemosensory–neurosecretory cells in the annelid larval apical organ. The researchers found that synthetic MIPs can induce the settlement of P. dumerilii larvae, and they demonstrate by morpholino-mediated knockdown that MIP signals via a G protein-coupled receptor to trigger settlement . These results reveal a role for a conserved MIP receptor–ligand pair in regulating marine annelid settlement. In this study, we revealed that MIP, FxFamide, FILamide and FRWamide all localized in the region of the apical organ (Fig. 6), like MIPs expression pattern in P. dumerilii, indicating that these neuropeptides may also be involved in triggering larval settlement. In addition, the expressions of the MIP, FRWamide and FxFamide were also detected at base of the abdomen chaetae (Fig. 6c, f and i). The chaetae are important in locomotion, stabilization during peristalsis, and sensing the environment in annelids , and have been implicated in assisting movement and stabilizing body segments within the tube for worms living in burrows or tubes [58, 59]. However, in sediment-dwellers they contact the inside walls of tubes or burrows. The cantilever nature of capillary chaetae and their astounding breadth of flexural stiffness suggest that they could be very effective at transmitting specific mechanical information about their surroundings to the body of the worm [34, 57]. Therefore, we propose that MIP, FRWamide and FxFamide located in base of the abdomen chaetae may help U. unicinctus larvae crawl on the sediment surface or explore the bottom and eventually contribute to larval settlement. Meanwhile, there are also some limitations to this study which need to be acknowledged. Firstly, this study only uses transcriptomic data to identify neuropeptides, therefore these predicted pNPs have not been confirmed by mass spectrometry to show that they are definitely released as signaling peptides in the worm. Secondly, we only used WISH technology to explore the location of the pNP genes at the mRNA level. Therefore, some issues including mass spectrometry, immunohistochemistry, western blot and other functional studies remained to be investigated in the future.
In this study, we identified 54 pNP genes in U. unicinctus larvae and adult transcriptome databases based on BLAST and NpSearch prediction, and suggested that the neuropeptide system of U. unicinctus is very close to that of annelids according to their phylogenetic distribution. Based on the expression data of pNP genes in the transcriptome of U. unicinctus larvae, four pNPs that may be associated with larval settlement were selected to further investigation. qRT-PCR results showed that the four pNPs were indeed highly expressed in competent larvae. WISH results indicated all the four pNPs were expressed in the region of the apical organ of the larva, and the positive signals were also detected in the ciliary bands and abdomen chaetae. These results imply that the four pNPs may be involved in sensing settlement signals and regulating larval ciliary locomotor and the movement of chaeta, and eventually may play an important role in the settlement of U. unicinctus larvae. Our findings provide some basic data for investigate the complex regulatory mechanisms in larval settlement of marine benthic invertebrates.
Animals and sampling
Adult U. unicinctus were obtained from Jiutian aquatic products market in Zhifu District of Yantai, China. Sperm and ova were obtained by dissecting the nephridia (gonaduct) of male and female, respectively. Artificial insemination was conducted through mixing the sperm and ova with a ratio of 10: 1 in filtered sea water (FSW). The fertilized eggs were reared in FSW (17 °C, pH 7.7, and salinity 30 PSU), and the hatched larvae were fed with single-cell algae (Isochrysis galbana, Chlorella vulgaris and Chaetoceros muelleri). The larvae at different stages were sampled and fixed in 4% paraformaldehyde for 15 h at 4 °C, and then dehydrated with serial methanol (25, 50, 75 and 100%) and stored in 100% methanol at − 30 °C for immunofluorescent histochemistry and whole-mount in situ hybridization analysis. The larvae from five developmental stages, late-trochophore (LT, 25 dpf), pre-competent larva (PL, 32 dpf) correspond to early-segmented larva (ES), competent larva (CL, 35 dpf) correspond to segmented larva (SL) which is the fully developed larvae prior to settlement and metamorphosis , post-competent larva (POL, 38 dpf) correspond to late-segmented larva (LS), and worm-shaped larva (WL, 42 dpf) were collected (Fig. 4a), frozen with liquid nitrogen immediately and then stored at − 80 °C, respectively for total RNA extraction. Three biological replicates from each developmental stage were prepared.
Identification, classification and sequence alignment of neuropeptide precursors
Data from U. unicinctus larval transcriptomes [29, 60] and adult transcriptome  were used in this study. To search for transcripts encoding putative neuropeptides or peptide hormone precursor proteins in U. unicinctus, the homologous sequences previously identified in annelids (Platynereis dumerilii , Capitella capitate  and Helobdella robusta ), molluscs (Patinopecten yessoensis , Pinctata fucata , Lottia gigantea , Crassostrea gigas  and Deroceras reticulatum ) and echinoderms (Asterias rubens , Ophionotus Victoria , Strongylocentrotus purpuratus  and Apostichopus japonicus ) were downloaded from NCBI and used as queries in tBLASTn searches of the assembled U. unicinctus transcriptome database using BioEdit v7.0.9 with an E-value cutoff of 1e-5. Open reading frames (ORFs) in these mRNA sequences of potential neuropeptides were identified using DNASTAR v7.1. The resultant protein sequences were further evaluated based on (i) the presence of a putative N-terminal signal peptide identified by SignalP 4.1 and Signal-3 L 2.0 [69, 70]; (ii) the presence of putative monobasic or dibasic cleavage sites flanking the putative bioactive peptides according to the existing neuropeptide cleavage motifs ; (iii) the presence of a C-terminal glycine residue which is a putative amidation site, and (iv) the presence of cysteine residues which can form disulfide linkages.
Furthermore, we used a neuropeptide-prediction tool NpSearch (https://github.com/wurmlab/npsearch) to identify the putative neuropeptide precursors with low sequence similarity to known precursors based on various characteristics (signal peptide, cleavage sites, C-terminal glycine and cysteine residues). Functional annotation of the identified neuropeptide precursors was finally conducted by searching against NCBI non-redundant protein sequence (nr) database using BLASTx algorithm with the E-value of 1e-5.
The classification of U. unicinctus neuropeptides is mainly based on the researches of Jékely  and Conzelmann , and we have also updated the classification status of LFRYamide [72, 73], Cerebrin [28, 63, 74, 75] and RGWamide [30, 44, 45, 62,63,64, 76, 77] according to recent studies.
The neuropeptide precursor homologous sequences from other species were collected from GenBank. Multiple alignments were conducted using ClustalW , and the results were annotated with GeneDoc (https://genedoc.software.informer.com). The frequency of each amino acid in the alignment result was presented using the online tool WebLogo . The hierarchical clustering of the pNP genes according to their FPKM values in the U. unicinctus larval transcriptome  was performed by an online tool (https://www.omicshare.com/tools/).
RNA isolation, cDNA synthesis and quantitative real-time PCR (qRT-PCR)
Total RNA from each stored larval sample was isolated using MicroElute®Total RNA Kit (Omega, Norcross, USA) according to the manufacture’s instruction. The cDNA was synthesized for each sample using Prime-Script™ RT reagent Kit (TaKaRa, Dalian, China) with the gene specific primers (Supplementary Table S1) designed using Primer Premier 5.0 according to their predicted CDS sequences. The amplifications were performed with SYBR Premix Ex Taq kit (TaKaRa, Dalian, China) in LightCycler 480 Real-Time PCR system. The PCR mixture consisted of 10 μl SYBR Premix Ex Taq II, 1 μl template cDNA, 1 μl forward primer (10 μM), 1 μl reverse primer (10 μM) and 7 μl ddH2O. The qRT-PCR condition was: denature at 95 °C for 30 s, followed by 39 cycles of 5 s at 95 °C, and 60 °C for 30 s. Each sample was run in 3 technical replicates. The relative expression levels were normalized to the reference gene ATPase , and expression ratios were calculated using the 2–ΔΔCt method. The experimental data were presented as mean ± standard deviation from three samples with three parallel repetitions, and all RT-PCR assays were validated in compliance with “the MIQE guidelines” . Significant differences between means were tested using one-way analysis of variance (ANOVA) followed by Tukey’s HSD test with SPSS software 18.0 (SPSS Inc., Chicago, USA). The significance level was set at p < 0.05.
Whole-mount in situ hybridization (WISH)
Specific fragments from cDNA of each neuropeptide genes were amplified using the primers with T7 or Sp6 promoter sequence at their 5′ ends (Supplementary Table S1). DIG-labeled RNA probes were prepared using the DIG RNA Labeling Kit SP6/T7 (Roche, Basel, Switzerland) with the PCR products as templates. WISH was carried out with the following protocols: the larvae were rehydrated with PBT (PBS + 0.1% Tween-20), and treated with 200 ng/μl proteinase K for 30 min to optimize hybridization; pre-hybridization was carried out for 6 h at 60 °C, then the probe was added and incubated at last 16 h at 60 °C; after the excess probe was removed by several rinses in hybridization buffer (50% formamide, 5 × SSC, 0.1% Tween, 9.2 mM citric acid for adjustment to pH 6.0, 50 μg/mL heparin, 500 μg/mL yeast RNA), the non-specific binding sites in the larval cells were blocked using the blocking buffer; the samples were incubated with the Anti-DIG-AP antibody (Roche, Basel, Switzerland) for 16 h at 4 °C; finally, the samples were stained in an NBT/BCIP staining solution (Roche, Basel, Switzerland) and kept in the dark for 1 h, and then dehydrated. The results were observed and photographed using a Nikon E80i microscope (Nikon, Tokyo, Japan). Drawings and final panels were designed using Adobe Photoshop (Adobe, San Jose, CA, USA).
The larva samples were rehydrated in a gradient methanol series (100, 75, 50 and 25%), and treated with 3% bovine serum albumin (BSA) (Shanghai biotechnology Technology, Shanghai, China) diluted by PBT (pH 7.4). Subsequently, the samples were transferred into primary antibody (Anti-5-hydroxytryptamine, an antibody produced in rabbit, Sigma, Jaffrey, USA) diluted 1: 200 in BSA and incubated overnight at 4 °C on a nutator. Afterwards, the samples were rinsed in PBT for 2 h and incubated subsequently with secondary fluorochrome conjugated antibodies (donkey anti-rabbit Alexa Fluor 488, Invitrogen, CA, USA) diluted 1: 300 in PBT for 2 h. At last, the larvae were washed six times in PBT and incubated in PBT with 2.5% DAPI (Solarbio, Beijing, China) in the dark for 15 min to label cell nuclei. Negative controls were obtained by pre-immune serum in order to check for antibody specificity. All the samples were analyzed with the confocal laser-scanning microscope Nikon A1RSi (Nikon, Tokyo, Japan). Drawings and final panels were designed using Adobe Photoshop (Adobe, San Jose, CA, USA).
Availability of data and materials
The datasets analysed during the current study are available in the NCBI SRA repository (SRX397931, SRX398497, SRX4526076, SRX4526077, SRX4526078, SRX2999430, SRX4526079, SRX4526072, SRX4526073, SRX4526074, SRX4526075, SRX4526080, SRX2999431, SRX4526070, SRX4526071, SRX4526081).
Quantitative real-time PCR
Open reading frame
Early- segmented larva
Ventral nerve cord
Whole-mount mRNA in situ hybridization
Filtered sea water
Practical salinity units
One-way analysis of variance
Phosphate belanced solution
Hadfield MG. Why and how marine-invertebrate larvae metamorphose so fast. Semin Cell Dev Biol. 2000;11:437–43.
Hodin J. Expanding networks: signaling components in and a hypothesis for the evolution of metamorphosis. Integr Comp Biol. 2006;46:719–42.
Feng D, Ke C, Lu C, Li S. The influence of temperature and light on larval pre-settlement metamorphosis: a study of the effects of environmental factors on pre-settlement metamorphosis of the solitary ascidian Styela canopus. Mar Freshw Behav Physiol. 2010;43:11–24.
Hadfield MG, Koehl MAR. Rapid behavioral responses of an invertebrate larva to dissolved settlement cue. Biol Bull. 2004;207:28–43.
Qian PY. Larval settlement of polychaetes. Hydrobiologia. 1999;402:239–53.
Walters LJ, Miron G, Bourget E. Endoscopic observations of invertebrate larval substratum exploration and settlement. Mar Ecol Prog Ser. 1999;182:95–108.
Kornijów R, Pawlikowski K, Drgas A, Rolbiecki L, Rychter A. Mortality of post-settlement clams Rangia cuneata (Mactridae, Bivalvia) at an early stage of invasion in the Vistula Lagoon (South Baltic) due to biotic and abiotic factors. Hydrobiologia. 2018;811:207–19.
Rodriguez SR, Ojeda FP, Inestrosa NC. Settlement of benthic marine invertebrates. Mar Ecol Prog Ser. 1993;97:193–207.
Conzelmann M, Williams EA, Tunaru S, Randel N, Shahidi R, Asadulina A, et al. Conserved MIP receptor–ligand pair regulates Platynereis larval settlement. Proc Natl Acad Sci U S A. 2013;110:8224–9.
Hadfield MG, Meleshkevitch EA, Boudko DY. The apical sensory organ of a gastropod veliger is a receptor for settlement cues. Biol Bull. 2000;198:67–76.
Beer AJ, Moss M, Thorndyke M. Development of serotonin-like and SALMFamide-like Immunoreactivity in the nervous system of the sea urchin Psammechinus miliaris. Biol Bull. 2001;200:268–80.
Dickinson AJG, Nason J, Croll RP. Histochemical localization of FMRFamide, serotonin and catecholamines in embryonic Crepidula fornicata (Gastropoda, Prosobranchia). Zoomorphology. 1999;119:49–62.
Gruhl A. Serotonergic and FMRFamidergic nervous systems in gymnolaemate bryozoan larvae. Zoomorphology. 2009;128:135–56.
Voronezhskaya EE, Tsitrin EB, Nezlin LP. Neuronal development in larval polychaete Phyllodoce maculata (Phyllodocidae). J Comp Neurol. 2003;455:299–309.
Watanabe H, Fujisawa T, Holstein TW. Cnidarians and the evolutionary origin of the nervous system. Develop Growth Differ. 2009;51:167–83.
Bargmann CI. Beyond the connectome: how neuromodulators shape neural circuits. Bioessays. 2012;34:458–65.
Gruber CW. Physiology of invertebrate oxytocin and vasopressin neuropeptides. Exp Physiol. 2014;99:55–61.
Nässel DR, Winther ÅME. Drosophila neuropeptides in regulation of physiology and behavior. Prog Neurobiol. 2010;92:42–104.
Chun JY, Korner J, Kreiner T, Scheller RH, Axel R. The function and differential sorting of a family of aplysia prohormone processing enzymes. Neuron. 1994;12:831–44.
Eipper BA, Stoffers DA, Mains RE. The biosynthesis of neuropeptides: peptide alpha-amidation. Annu Rev Neurosci. 1992;15:57–85.
Conzelmann M, Offenburger SL, Asadulina A, Keller T, Münch TA, Jékely G. Neuropeptides regulate swimming depth of Platynereis larvae. Proc Natl Acad Sci U S A. 2011;108:1174–83.
Braubach OR, Dickinson AJ, Evans CC, Croll RP. Neural control of the velum in larvae of the gastropod, Ilyanassa obsoleta. J Exp Biol. 2006;209:4676–89.
Dickinson AJG, Croll RP. Development of the larval nervous system of the gastropod Ilyanassa obsoleta. J Comp Neurol. 2003;466:197–218.
Dyachuk V, Odintsova N. Development of the larval muscle system in the mussel Mytilus trossulus (Mollusca, Bivalvia). Develop Growth Differ. 2009;51:69–79.
Kempf SC, Masinovsky B, Willows AD. A simple neuronal system characterized by a monoclonal antibody to SCP neuropeptides in embryos and larvae of Tritonia diomedea (Gastropoda, Nudibranchia). J Neurobiol. 1987;18:217–36.
Hadfield MG, Carpizo-Ituarte EJ, Del Carmen K, Nedved BT. Metamorphic competence, a major adaptive convergence in marine invertebrate larvae. Am Zool. 2001;41:1123–31.
Huan P, Wang H, Liu B. A label-free proteomic analysis on competent larvae and juveniles of the pacific oyster Crassostrea gigas. PLoS One. 2015;10:e0135008.
Jékely G. Global view of the evolution and diversity of metazoan neuropeptide signaling. Proc Natl Acad Sci U S A. 2013;110:8702–7.
Hou X, Wei M, Li Q, Zhang T, Zhou D, Kong D, et al. Transcriptome analysis of larval segment formation and secondary loss in the Echiuran worm Urechis unicinctus. Int J Mol Sci. 2019;20:1806.
Conzelmann M, Williams EA, Krug K, Franz-Wachtel M, Macek B, Jékely G. The neuropeptide complement of the marine annelid Platynereis dumerilii. BMC Genomics. 2013;14:906.
Watthanasurorot A, Söderhäll K, Jiravanichpaisal P, Söderhäll I. An ancient cytokine, astakine, mediates circadian regulation of invertebrate hematopoiesis. Cell Mol Life Sci. 2011;68:315–23.
Aravind L, Koonin EV. A colipase fold in the carboxy-terminal domain of the Wnt antagonists–the Dickkopfs. Curr Biol. 1998;8:477–8.
Struck TH, Schult N, Kusen T, Hickman E, Bleidorn C, Mchugh D, et al. Annelid phylogeny and the status of Sipuncula and Echiura. BMC Evol Biol. 2007;7:57.
Hessling R. Metameric organisation of the nervous system in developmental stages of Urechis caupo (Echiura) and its phylogenetic implications. Zoomorphology. 2002;121:221–34.
Tilic E, Lehrke J, Bartolomaeus T. Homology and evolution of the chaetae in Echiura (Annelida). PLoS One. 2015;10:e0120002.
Schoofs L, Holman GM, Hayes TK, Nachman RJ, De Loof A. Isolation, identification and synthesis of locustamyoinhibiting peptide (LOM-MIP), a novel biologically active neuropeptide from Locusta migratoria. Regul Pept. 1991;36:111–9.
Schoofs L, Veelaert D, Vanden Broeck J, De Loof A. Immunocytochemical distribution of Locustamyoinhibiting peptide (Lom-MIP) in the nervous system of Locusta migratoria. Regul Pept. 1996;63:171–9.
Hua YJ, Tanaka Y, Nakamura K, Sakakibara M, Nagata S, Kataoka H. Identification of a prothoracicostatic peptide in the larval brain of the silkworm, Bombyx mori. J Biol Chem. 1999;274:31169–73.
Lorenz MW, Kellner R, Hoffmann KH. A family of neuropeptides that inhibit juvenile hormone biosynthesis in the cricket, Gryllus bimaculatus. J Biol Chem. 1995;270:21103–8.
Minakata H, Ikeda T, Muneoka Y, Kobayashi M, Nomoto K. WWamide-1, −2 and −3: novel neuromodulatory peptides isolated from ganglia of the African giant snail, Achatina fulica. FEBS Lett. 1993;323:104–8.
Williamson M, Lenz C, Winther ME, Nässel DR, Grimmelikhuijzen CJ. Molecular cloning, genomic organization, and expression of a B-type (cricket-type) allatostatin preprohormone from Drosophila melanogaster. Biochem Bioph Res Co. 2001;281:544–50.
Williams EA, Conzelmann M, Jékely G. Myoinhibitory peptide regulates feeding in the marine annelid Platynereis. Front Zool. 2015;12:1.
Huybrechts J, Bonhomme J, Minoli S, Prunier-Leterme N, Dombrovsky A, Abdel-Latief M, et al. Neuropeptide and neurohormone precursors in the pea aphid, Acyrthosiphon pisum. Insect Mol Biol. 2010;19:87–95.
Veenstra JA. Neuropeptide evolution: neurohormones and neuropeptides predicted from the genomes of Capitella teleta and Helobdella robusta. Gen Comp Endocrinol. 2011;171:160–75.
Zhang M, Wang Y, Li Y, Li W, Li R, Xie X, et al. Identification and characterization of neuropeptides by Transcriptome and proteome analyses in a bivalve Mollusc Patinopecten yessoensis. Front Genet. 2018;9:197.
Hessling R. Novel aspects of the nervous system of Bonellia viridis (Echiura) revealed by the combination of immunohistochemistry, confocal laser-scanning microscopy and three-dimensional reconstruction. Hydrobiologia. 2003;496:225–39.
Hessling R, Westheide W. Are Echiura derived from a segmented ancestor? Immunohistochemical analysis of the nervous system in developmental stages of Bonellia viridis. J Morphol. 2002;252:100–13.
Conzelmann M, Jékely G. Antibodies against conserved amidated neuropeptide epitopes enrich the comparative neurobiology toolbox. Evodevo. 2012;3:23.
Heyland A, Moroz LL. Signaling mechanisms underlying metamorphic transitions in animals. Integr Comp Biol. 2006;46:743–59.
Laudet V. The origins and evolution of vertebrate metamorphosis. Curr Biol. 2011;21:726–37.
Truman JW, Riddiford LM. The origins of insect metamorphosis. Nature. 1999;401:447–52.
Kempf SC, Page LR, Pires A. Development of serotonin-like immunoreactivity in the embryos and larvae of nudibranch mollusks with emphasis on the structure and possible function of the apical sensory organ. J Comp Neurol. 1997;386:507–28.
Rentzsch F, Fritzenwankerv JH, Scholz CB, Technau U. FGF signalling controls formation of the apical sensory organ in the cnidarian Nematostella vectensis. Development. 2008;135:1761–9.
Voronezhskaya EE, Khabarova MY. Function of the apical sensory organ in the development of invertebrates. Dokl Biol Sci. 2003;390:231–4.
Voronezhskaya EE, Khabarova MY, Nezlin LP. Apical sensory neurones mediate developmental retardation induced by conspecific environmental stimuli in freshwater pulmonate snails. Development. 2004;131:3671–80.
Chia FS, Buckland-Nicks J, Young CM. Locomotion of marine invertebrate larvae: a review. Can J Zool. 1984;62:1205–22.
Merz RA, Woodin SA. Polychaete chaetae: function, fossils, and phylogeny. Integr Comp Biol. 2006;46:481–96.
Knight-Jones P, Fordy MR. Setal structure, functions and interrelationships in Spirorbidae (Polychaeta, Sedentaria). Zool Scr. 1979;8:119–38.
Sendall KA, Fontaine AR, O’Foighil D. Tube morphology and activity patterns related to feeding and tube building in the polychaete Mesochaetopterus taylori Potts. Can J Zool. 1995;73:509–17.
Park C, Han YH, Lee SG, Ry KB, Oh J, Kern EMA, et al. The developmental transcriptome atlas of the spoon worm Urechis unicinctus (Echiurida: Annelida). GigaScience. 2018;7:1–7.
Liu X, Zhang L, Zhang Z, Ma X, Liu J. Transcriptional response to sulfide in the Echiuran worm Urechis unicinctus by digital gene expression analysis. BMC Genomics. 2015;16:829.
Stewart MJ, Favrel P, Rotgans BA, Wang T, Zhao M, Sohail M, et al. Neuropeptides encoded by the genomes of the Akoya pearl oyster Pinctata fucata and Pacific oyster Crassostrea gigas: a bioinformatic and peptidomic survey. BMC Genomics. 2014;15:840.
Veenstra JA. Neurohormones and neuropeptides encoded by the genome of Lottia gigantea, with reference to other mollusks and insects. Gen Comp Endocrinol. 2010;167:86–103.
Ahn SJ, Martin R, Rao S, Choi MY. Neuropeptides predicted from the transcriptome analysis of the gray garden slug Deroceras reticulatum. Peptides. 2017;93:51–65.
Semmens DC, Mirabeau O, Moghul I, Pancholi MR, Wurm Y, Elphick MR. Transcriptomic identification of starfish neuropeptide precursors yields new insights into neuropeptide evolution. Open Biol. 2016;6:150224.
Zandawala M, Moghul I, Yañez Guerra LA, Delroisse J, Abylkassimova N, Hugall AF, et al. Discovery of novel representatives of bilaterian neuropeptide families and reconstruction of neuropeptide precursor evolution in ophiuroid echinoderms. Open Biol. 2017;7:170129.
Rowe ML, Elphick MR. The neuropeptide transcriptome of a model echinoderm, the sea urchin Strongylocentrotus purpuratus. Gen Comp Endocrinol. 2012;179:331–44.
Rowe ML, Achhala S, Elphick MR. Neuropeptides and polypeptide hormones in echinoderms: new insights from analysis of the transcriptome of the sea cucumber Apostichopus japonicus. Gen Comp Endocrinol. 2014;197:43–55.
Petersen TN, Brunak S, Von Heijne G, Nielsen HH. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Zhang YZ, Shen HB. Signal-3L 2.0: a hierarchical mixture model for enhancing protein signal peptide prediction by incorporating residue-domain cross-level features. J Chem Inf Model. 2017;57:988–99.
Liu F, Baggerman G, D'Hertog W, Verleyen P, Schoofs L, Wets G. In silico identification of new secretory peptide genes in Drosophila melanogaster. Mol Cell Proteomics. 2006;5:510–22.
Fadda M, Hasakiogullari I, Temmerman L, Beets I, Zels S, Schoofs L. Regulation of feeding and metabolism by neuropeptide F and short neuropeptide F in invertebrates. Front Endocrinol. 2019;10:64.
Yañez-Guerra LA, Zhong X, Moghul I, Butts T, Zampronio CG, Jones AM, et al. Echinoderms provide missing link in the evolution of PrRP/sNPF-type neuropeptide signalling. Elife. 2020;9:e57640.
Nässel DR, Zandawala M. Recent advances in neuropeptide signaling in Drosophila, from genes to physiology and behavior. Prog Neurobiol. 2019;179:101607.
Mirabeau O, Joly JS. Molecular evolution of peptidergic signaling systems in bilaterians. Proc Natl Acad Sci U S A. 2013;110:2028–37.
Bauknecht P, Jékely G. Large-scale combinatorial deorphanization of Platynereis neuropeptide GPCRs. Cell Rep. 2015;12:684–93.
Thiel D, Franz-Wachtel M, Aguilera F, Hejnol A. Xenacoelomorph neuropeptidomes reveal a major expansion of neuropeptide systems during early bilaterian evolution. Mol Biol Evol. 2018;35:2528–43.
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.
Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14:1188–90.
Wei M, Lu L, Wang Q, Kong D, Zhang T, Qin Z, et al. Evaluation of suitable reference genes for normalization of RT-qPCR in Echiura (Urechis unicinctus) during developmental process. Russ J Mar Biol. 2019;45:464–9.
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.
We would like to thank the anonymous reviewers for their kind and helpful comments on the original manuscript.
This work was support by Fundamental Research Funds for the Central Universities (202064006), the National Natural Science Foundation of China (31572601), China Postdoctoral Science Foundation (2020 M680095) and Shandong province science outstanding Youth Fund. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The collection and handing of the U. unicinctus and experimental procedures were performed in accordance with the Guidelines for Experimental Animals of the Ministry of Science and Technology (Beijing, China) and approved by the Institutional Animal Welfare Committee of the College of Marine Life Sciences, Ocean University of China (Approval number: 2016016).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Specific primers used in this research.
. Detailed information of the identified neuropeptide precursors in U. unicinctus.
Structures of U. unicinctus pNPs and identified repetitive peptide motifs.
Expression trends of the neuropeptide precursor genes in U. unicinctus larval transcriptome.
The negative controls of MIP, FRWa, FxFa and FILa in U. unicinctus larvae detected by whole-mount in situ hybridization.
About this article
Cite this article
Hou, X., Qin, Z., Wei, M. et al. Identification of the neuropeptide precursor genes potentially involved in the larval settlement in the Echiuran worm Urechis unicinctus. BMC Genomics 21, 892 (2020). https://doi.org/10.1186/s12864-020-07312-4
- Urechis unicinctus
- Neuropeptide precursor
- Larval settlement