Identification of Neuropeptides and Their Receptors in the Ectoparasitoid, Habrobracon hebetor

Neuropeptides are a group of signal molecules that regulate many physiological and behavioral processes by binding to corresponding receptors, most of which are G-protein-coupled receptors (GPCRs). Using bioinformatic methods, we screened genomic and transcriptomic data of the ectoparasitoid wasp, Habrobracon hebetor, and annotated 34 neuropeptide candidate precursor genes and 44 neuropeptide receptor candidate genes. The candidate neuropeptide genes were found to encode all known insect neuropeptides except allatotropin, neuropeptide F, pigment dispersing factor, and CCHamides. When compared with the endoparasitic wasp Pteromalus puparum and the ectoparasitic wasp Nasonia vitripennis, trissin and FMRFamide were found only in H. hebetor. A similar result held for the neuropeptide receptor genes, for the receptors were found in H. hebetor except the receptors of CCHamides and neuroparsin. Furthermore, we compared and analyzed the differences in neuropeptides in eight Braconidae wasps and identified natalisin in H. hebetor, Diachasma alloeum, Fopius arisanus and Microplitis demolitor, but not in the other wasps. We also analyzed the transcriptome data and qRT-PCR data from different developmental stages and tissues to reveal the expression patterns of the neuropeptides and their receptors. In this study, we revealed composition of neuropeptides and neuropeptide receptors in H. hebetor, which may contribute to future neurobiological studies.


INTRODUCTION
Neuropeptides and peptide hormones function as key signals in orchestrating the regulation of numerous physiological processes and behaviors (Schmitt et al., 2015;Wang L. et al., 2018). Neuropeptides constitute a large and diverse class of signaling molecules that are produced by various types of neurons, neurosecretory cells, endocrine cells and other cells (Nassel et al., 2019). Being a kind of endogenous active substance, neuropeptides widely occur in multicellular biological nerve tissues, where they are involved in the functioning of the nervous system, mostly acting upon neuropeptide receptors (Schoofs et al., 2017).
In insects, neuropeptides and their receptors also play significant roles in controlling myriad physiological processes, including development, reproduction, feeding, homeostasis, courtship, circadian rhythm, olfaction, movement, water, and ion homeostasis, among others (Schoofs et al., 2017;Dickinson et al., 2019;Pandit et al., 2019). Due to their specificity and high activity at low concentrations, neuropeptides have also been investigated as potential leads for developing new environmentally-friendly pest control agents (Xiong et al., 2019). Importantly, many neuropeptides are capable of pleiotropic actions such as enabling them to function as neuromodulators, co-transmitters or circulating hormones, some of effects may be exerted simultaneously (Nassel et al., 2019). Neuropeptides are produced from larger precursor proteins called pre-propeptides. A given insect species can harbor 50 neuropeptide precursor-encoding genes, with some species featuring greater complementarity of precursors than others (Yeoh et al., 2017).
Neuropeptides mediate their biological actions via interactions with specific receptors on the cell surface (Ma et al., 2017). Most of these receptors are G-protein-coupled receptors (GPCRs), which have a similar structure characterized by seven transmembrane domains that are highly conserved through evolution, and they constitute the largest superfamily of cell membrane-spanning proteins (Jastrzebska, 2017). Functionally, these proteins are able to recognize extracellular transmitter molecules (Bao et al., 2018).
Habrobracon hebetor is an ectoparasitoid distributed worldwide, with an extensive host range and a rapid life cycle (Ghimire and Phillips, 2010). It provides excellent biological control services in terrestrial agroecosystems and has become a specific agent for effectively managing multiple lepidopteran pests (Ghimire and Phillips, 2010). For example, H. hebetor is the natural enemy of Heliocheilus albipunctella, parasitizing the latter's larvae at a rate of 50-78% when released in pearl millet fields (Gahukar and Ba, 2019). With the development of next-generation sequencing technology, here we identified a complete set of neuropeptide genes and interrelated receptor genes via genomic and transcriptomic sequencing of H. hebetor. Building on this, we then compared them with those annotated in representative species of different insect orders, including the wasp parasitoids Nasonia vitripennis (Hauser et al., 2010) and Pteromalus puparum , the bee Apis mellifera (Hummon et al., 2006), the fly Drosophila melanogaster (Broeck, 2001), the moth Bombyx mori (Roller et al., 2008), the beetle Tribolium castaneum (Li et al., 2008), and an hemipteran insect: Nilaparvata lugens (Tanaka et al., 2014). In this study, we identified 34 neuropeptide candidate genes and 44 neuropeptide candidate receptor genes in H. hebetor. The expression profiles of neuropeptides and neuropeptide receptors at different developmental stages and various tissues were determined. These results allow us to compare the neuropeptidergic signaling systems across different insect species, to provide relevant information for further functional studies in H. hebetor. The findings also provide us novel practical insights for developing new, environmentally friendly insecticides with high effectiveness to against insect pests but special while remaining safety to for their specific parasitoid wasps.

Sample Collection and Sequencing
In a recent study, we carrier out the sequencing and assembly of the H. hebetor genomem, its assembled genome is 131.6 Mb in size with a contig N50 of 1.63 Mb (Ye et al., 2020). Here, we prepared embryos, larvae, and other samples from the wasps for RNA-Seq analysis. Embryos spawned by mating female wasps were obtained from the newly parasitized fifth-instar larvae of P. interpunctella within 2 h. Larvae spawned by mating female wasps were isolated from parasitized host larvae at the third instar. White pupae and newly emerged wasp adults of both sexes were collected. Moreover, we also prepared samples of salivary glands, venom glands, the carcass without salivary glands of larvae, and the carcass without the venom gland of female adults. All samples were collected, washed with phosphate buffer saline (10 mM of 1× PBS; pH 7.4), and then ground with TRIzol reagent (Invitrogen, USA) for subsequent experiments. The construction of complementary DNA (cDNA) libraries was performed by Nextomics Biosciences (Nextomics, Wuhan, China), and then sequenced using an Illumina HiSeq TM system (Illumina, NEB, United States). We sequenced two 100-bp paired-end lanes. We used the NGSQCToolKit (v2.3.3) for data filtering (Patel and Jain, 2012), and the FastQC (http://www.bioinformatics.babraham.ac. uk/projects/fastqc/) for quality control of data. The Q30 of all the samples were above 94%, and we then generated ∼144.3 million clean reads, for a total of 21556.9 MB for all sequences. The two un-filtered paired-end lanes of each sequence have been deposited as a series, under the accession number PRJNA642006 at NCBI's GEO database or at the NCBI Short Read Archive under submission number SUB7674593.

Identification of Neuropeptides and Neuropeptide Receptors
We used the already-reported neuropeptides and neuropeptide receptors of P. puparum, N. vitripennis, D. melanogaster, A. mellifera, B. mori, and other insects as queries to search for candidate sequences in H. hebetor genome database. TBLASTN were used to search for and acquire candidate sequences with the criteria of identity ≥30% and a threshold E-value of 10 −3 . We amplified 5000-bp upstream and downstream of the corresponding genome regions into high scoring segment pairs (HSPs). Gene predictions from the HSPs were carried out using the GeneWise tools (Birney and Durbin, 2000), and then confirmed through a BLASTP search with the non-redundant protein sequence (NR) at NCBI. After gene identification, signal peptides of putative neuropeptides were identified using the SignalP v4.0 software tool (http://www.cbs.dtu.dk/services/ SignalP/). We identified candidate genes with signal peptides as precursor neuropeptide genes. The transmembrane domains of putative neuropeptides receptors were verified by HMMER (Finn et al., 2011) in the Pfam (http://pfam.xfam.org/) database. Identification of neuropeptides of the other seven wasps was conducted applying the same procedure.

Sequence Alignments and Phylogenetic Analysis
The sequences of neuropeptide precursors and neuropeptide receptors from H. hebetor and other insects were aligned by ClustalX2 software (http://www.clustal.org/clustal2/) and edited with GeneDoc v2.7 software (https://genedoc.software.informer. com/). According to an earlier study, the location of the disulfide bond is distinguishable in the results of the multiple sequence alignment (Kono et al., 1990). A phylogenetic analysis of neuropeptide receptors of H. hebetor and other insect species was used to study their evolutionary relationships, for which neighbor-joining trees were constructed in MEGA 7.0 with n = 1,000 bootstrap replicates (Kumar et al., 2016). The circle trees were drawn with FigTree v1.4.3 software and all trees were rooted by the D. melanogaster metabotropic glutamate receptor (DmCG11144).

Expression Profiling of the Neuropeptides and Neuropeptide Receptors
Gene expression levels for H. hebetor samples from different stages of both sexes were estimated via RNA-Seq using the Expectation Maximization (RSEM) (Li and Dewey, 2011), with expression levels calculated and conveyed as fragments per kilobase of transcript per million (FPKM) values. The expression profiles of each neuropeptide and neuropeptide receptor were visualized by using TBtools (Chen et al., 2020).

Quantitative Real-Time Polymerase Chain Reaction (qPCR) Analyses
We performed qPCR using SYBR R Premix Ex Taq TM II (Tli RNaseH Plus) (Takara, Japan) on the CFX Connect TM Real Time Detection System (Bio-Rad). Total RNA was extracted with TRIzol and quantified by a NanoDrop (Thermo Scientific) spectrophotometer, after which the cDNA was synthesized using a PrimeScript TM One Step RT-PCR Kit (Takara, Japan). Then, in a 25-µL reaction volume, 1 ng cDNA was used as a template for which primers were designed by Primer3 web v4.1.0 and Primer-Blast (http://www.ncbi.nlm.nih.gov/tools/ primer-blast/). Initially, we used the reverse transcription polymerase chain reaction (RT-PCR) and sequencing of its product to correct the nucleotide sequences of the selected genes for the qPCR verification. Multiple specific primers were designed for each tested gene. To determine the efficacy of the primers, we performed qPCR for the templates with serial dilutions that ranged from 10 to 10,000, respectively, and calculated their efficiency values. A dissociation curve was drawn spanning 60-95 • C at the end of each qPCR reaction, to verify the specificity of the used qPCR primers. Based on the results from these specificity and efficiency verifications, the suitable primers for gene expression profile determination were selected. We chose glyceraldehyde-3-phosphate dehydrogenase (GAPDH) to sever the reference gene (Li et al., 2019). Statistical analyses were implemented in GraphPad Prism 6 (https://www. graphpad.com/). To quantify the mRNA expression levels the 2 − Ct method was utilized (Livak and Schmittgen, 2001). The correlation of R-values was calculated by using data processing system (DPS) software (Tang and Zhang, 2013).
In identifying the genes encoding neuropeptides and their receptors in insects, it became apparent that at least three neuropeptide genes seemed prevalent: those encoding for corazonin (Crz), LK, and allatostatin A (AstA) (Veenstra, 2019a). We also annotated these three genes in H. hebetor in the present study ( Table 1). The number of neuropeptide precursor genes in H. hebetor was lower than that in non-parasitic insect species, such as B. mori, C. suppressalis, D. melanogaster, N. lugens, T. castaneum, and A. mellifera, probably due to its parasitic life cycle (Supplementary Table 1). Furthermore, we identified the lowest number of neuropeptides in M. cingulum ( Table 2). With the exception of F. arisanus, all the braconid wasps lacked an adipokinetic hormone/corazonin-related peptide (ACP) ( Table 2). In stark contrast, ACP, along with AT, inotocin, CCHa, and PDF, were absent in H. hebetor compared with P. puparum and N. vitripennis (Supplementary Table 1). This disparity likely exists because these three wasps are from two different families. Yet many differences were discerned in the comparison within Braconidae among C. vestalis, D. alloeum, D. collaris, F. arisanus, M. cingulum, M. demolitor, T. pretiosum, and H. hebetor. For example, we could not identify inotocin in H. hebetor and M. cingulum, but inotocin was found in the six other wasps we   (Wei et al., 2010), D. alloeum (Lowe and Eddy, 1997), D. collaris (Li et al., 2015), F. arisanus (Lowe and Eddy, 1997), M. cingulum (Yin et al., 2018), M. demolitor (Burke et al., 2014), and T. pretiosum (Lowe and Eddy, 1997). investigated ( Table 2). We speculate such differences may also relate to the species' host specificity and parasitic mode. Apart from H. hebetor, the other seven wasps are endoparasitoids, and both F. arisanus and T. pretiosum parasitize host eggs, whereas the other species parasitize larvae. The hosts of F. arisanus and D. alloeum belong to Diptera, while the other wasps can only parasitize Lepidoptera.

Multiple Sequence Alignments of Some Neuropeptides in H. hebtor
We annotated the eclosion hormone (EH) in H. hebetor ( Table 1). The alignment of EH sequences from various insects showed that H. hebetor EH contains all three disulfide bridges usually present in EHs (Figure 1). Notably, the sequences of some candidate neuropeptide genes were found to be incomplete, in that they did not harbor signal peptide, as in the case of the bursicon alpha subunit (Burα) and CNMamide (CNMa). Nevertheless, those sequences were able to align with the same genes in other insects (Figures 2, 3). The complete genes of the bursicon beta subunit (Burβ) and a bursicon receptor were found in H. hebetor; hence, we deemed Hheb01190 candidate neuropeptide gene. CNMa is a member of a conserved insect neuropeptide family recently identified in the genomes of D. melanogaster and other species, including Rhodnius prolixus (Traverso et al., 2016) so named after its C-terminal consensus motif (Jung et al., 2014).

Receptors for Neuropeptides in H. hebetor G Protein-Coupled Receptors (GPCRs) for Neuropeptides
A total of 44 putative neuropeptide receptors were identified in the transcriptomes of H. hebetor (Supplementary Material File 2). Of these receptors, 32 GPCRs belong to the A-family and another five to the B-family, two are leucine-rich repeat-containing GPCRs (LGRs), three are receptor guanylyl cyclase (RGCs), two are receptor tyrosine kinases (RTKs). If a likely receptor gene is identified in H. hebetor, this gives additional support for the presence of these  neuropeptides in H. hebetor. Those neuropeptides, such as CCHa, NPF, and PDF, which we could not identify, yet found a corresponding putative receptor in H. hebetor, are probably missing from this wasp. However, according to Roller et al. (2008), the genes encoding CCHamides and their receptors can be found in most sequenced insect genomes (Roller et al., 2008). Another reason for why some were missing in H. hebetor could be the large differences between their sequences in Braconidae vis-à-vis those in other species. Similarly, neuroparsin (NP) was annotated in H. hebetor, but its receptors were not found. These missing neuropeptides may require further experimental verification.

A-Family GPCRs
The A-family GPCRs are also known as the rhodopsin family (Jung et al., 2014). Here we identified 32 A-family GPCRs in H. hebetor (Figure 4). For example, one subassemblage includes the receptors for AKH, ACP, and Crz. Since those three and their corresponding GPCRs do not always occur together in all insects (Li et al., 2008). We tried to annotate all of them, but the receptor of ACP found by BLAST in H. hebetor did not cluster with the ACPs of B. mori, N. vitripennis, and P. puparum, and being positioned instead on the outskirts of this phylogenetic branch (Figure 4).
There are two AstC receptors positioned in H. hebetor one of the Cluster (Figure 4), which differs from P. puparum, namely Hheb023220.1 and Hheb044860.1. We also found two CAPA receptors in H. hebetor, whereas neither CAPA or its receptor occurred in N. vitripennis, and the receptors of CAPA and pyrokinin were clustered together in the phylogenetic tree. There are other orthologs evident as well, including RPRC004565, Hheb55420.1, Hheb008990.1, Hheb008990.1, and Hheb07740.1 (Figure 4), all of which are orphan receptors. Some of these A-family GPCRs in H. hebetor, N. vitripennis and P. puparum were clustered together, such as Hheb077040.1, the receptor of FMRF.

B-Family GPCRs
Family B (secretin-like) GPCRs include the calcitonin-like receptor (CTR), CRF-like diuretic hormone receptor (DHR) and pigment dispersing factor receptor (PDFR) . We found 5 GPCRs belonging to B-family GPCRs in H. hebetor, including one diuretic hormone 44 (DH44) receptor, one PDF receptor, one diuretic hormone 31 (DH31) receptor and two other B-family GPCRs (Figure 5). In the phylogenetic tree, Hheb039610.1 is closely related to PpB1 and NvB1, while Hheb104700.1 is closely related to both PpDHR and NvDHR (Figure 5). The two GPCRs were orthologs of B6 B7, while the orthologs of B3 and B4 are not present in H. hebetor (Figure 5). LGRs LGRs are receptors with key functions in organismal development and reproduction, a prime example being the bursicon receptor, which triggers the cuticle hardening and tanning in newly emerged insects (Van Hiel et al., 2012). Three distinct types of LGR (type A-C)are known to exist, distinguishable by their number of leucine-rich repeats (LRRs), their type-specific hinge region, and the presence or absence of a low density lipoprotein receptor domain class A (LDLa) motif (Van Hiel et al., 2012). In H. hebetor, a type C1 LRR was identified, it being an ortholog of the fruit fly's LGR3 and LGR4 ( Figure 6); likewise, a type B LRR was also identified (Figure 6) this being receptor of bursicon. However, the orthologs of Type A and Type C2 were not present in H. hebetor.

RGCs
RGCs are conserved homodimeric membrane proteins with an intracellular protein kinase and guanylyl cyclase domains that catalyze the formation of cGMP (Potter, 2011). Here we identified three neuropeptide receptors which belong to the RGC family. In addition, similar to the circumstances of P. puparum and N. vitripennis, we found an EH receptor and a NPLP receptor in H. hebetor. Notably, Hheb087320.1 is an ortholog of orphan receptor guanylyl cyclase 4 (OGC4) (Figure 7).

RTKs
The receptors of prothoracicotropic hormone (PTTH), ILP, and NP all belong to the family of RTKs. In this study, we identified two receptors in H. hebetor that belonged to the RTKs (Figure 8), both being orthologs of the ILP receptor (InR). A phylogenetic analysis indicated that the NP receptor and InR are closely related .

Expression Profiles of Neuropeptides and Neuropeptide Receptors
Based on H. hebetor's RNA-Seq data, the expression profiles of its neuropeptides and neuropeptide receptors at different developmental stages (embryo, larva, pupa, adult) were characterized (Figures 9, 10). The genes AstA, AstCC, FMRF, Crz, NTL, ILP1, and LK were expressed at low levels in nearly all stages, whereas the ion transport peptide (ITP) was highly expressed in every stage, especially the embryonic stage (Figure 9). Both ITG and tachykinin (TK) were highly expressed through development except in the embryonic stage (Figure 9). Expression of ecdysis triggering hormone (ETH) was greater in the pupal stage than other developmental stages (Figure 9). Stage-specific expression profiles of the neuropeptides indicated that mRNAs encoding ITP undergo a high level of expression in the embryos, in contrast to NP, RY and ETH, whose expression was higher in the pupal stage (Figure 9). Concerning OKA, ITG and DH31, they were expressed more in the adult stage (Figure 9). Finally, these expression profiles also revealed that TK was preferentially expressed in the male pupae of H. hebetor (Figure 9).
Using for inference the RNA-Seq of the wasp tissues salivary glands, venom glands, the carcass of larvae without salivary glands, and the carcass of female adults without the venom glands, we characterized their expression profiles. This showed that TK and ITP were expressed at higher levels in both the venom and salivary glands than in the corresponding residual tissue (Supplementary Figure 1). Most expressed in the venom gland tissues was ILP1, whereas it was barely expressed in salivary glands (Supplementary Figure 1). We also determined that PDFR was highly expressed in all developmental stages and tissues examined (Figure 10, Supplementary Figure 2), but we failed to annotate PDF in H. hebetor. The mRNAs encoding AstAR underwent a high level of expression in male adults ( Figure 10); conversely, it was the mRNAs encoding CAPAR2 that were most expressed in female adults (Figure 10). Many receptors, however, were not at all expressed in the salivary and venom glands (Supplementary Figure 2). From the cluster perspective, females and males differed little, indicating negligible sex specificity, but the adults were significantly separate from the pre-eclosion period, suggesting the expression of these neuropeptides before and after the wasp's emergence is different.
To verify the reliability of our RNA-Seq data, we randomly selected some genes for a qPCR analysis. These results revealed that the expression patterns of those genes generally agreed with that registered by RNA-Seq. For example, the developmental expression profiles ETH and NP obtained by qPCR were consistent with those registered by RNA-Seq (Figure 11). Accordingly, the positive correlation was extremely strong (Rvalues > 0.80), for the four genes: ETH (R = 0.9181, p = 0.0098), NP (R = 0.9420, p = 0.0049), DH31 (R = 0.9099, p = 0.0118), and ILP1 (R = 0.8381, p = 0.0372). The other selected genes had nearly identical expression trends recorded via both analyses (R-values > 0.3) (Figure 11). In addition, the qPCR results for AstCCC, TK and AstCCR each had the same peak in the expression profiles as determined by the RNA-Seq at different developmental stages (Supplementary Figure 3). Overall, these validation results indicate that our RNA-Seq data were robust.

DISCUSSION
As ancient molecules, neuropeptides and peptide hormones mediate cell-to-cell communication in all multicellular animals  (Li et al., 2008;Garczynski et al., 2019). To date, many neuropeptides have been identified and sequenced in several insect species (Yeoh et al., 2017). Insect neuropeptides show heterogeneous similarities to mammals, yet display clear similarity among insects taxa, despite the large phylogenetic distance across these huge and various classes (Yeoh et al., 2017). Indeed, it is well-established that the structure of the genes encoding neuropeptides and their receptors are highly conserved during evolutionary history, and this is not surprising because they are important regulators of various vital physiological processes (Li et al., 2008;Veenstra, 2019a,b). Neuropeptide precursors share few features in common, however. The only common feature is the presence of an aminoterminal signal peptide that directs the ribosomes synthesizing  neuropeptide precursors to the endoplasmic reticulum (ER). Most neuropeptides can interact with GPCRs, thereby generating an intracellular response (Clynen et al., 2010). Neuropeptides and their receptors regulate fundamental events in the insect life cycle, hence, they have been proposed as potential insecticides or targets to replace or complement the use of neurotoxic compounds against pests (Verlinden et al., 2014). In addition, potential mutations to the neuropeptide precursor or receptor genes are usually harmful, so resistance is not easy to evolve (Ons, 2017). Insect neuropeptides are interesting because their receptors may be reasonable targets for a new generation of insecticides (Veenstra, 2019a), and we should also study these related genes in natural enemy insects (predators, parasitoids) to avoid harming them with insecticides.
In this study, we annotated 33 neuropeptide genes in H. hebetor by Nr-annotation and homology searches. Inotocin was first discovered from the central nervous system (CNS) of Locusta migratoria, and later also identified in T. castaneum, N. vitripennis, P. puparum, C. vestalis, D. alloeum, D. collaris, F. arisanus, M. demolitor, T. pretiosum, and N. lugens ( Table 2,  Supplementary Table 1). But it is apparently the lost from A. mellifera (Hummon et al., 2006) and H. hebetor. The genes encoding two well-known neuropeptides, PDF and NPF were not found in the H. hebetor genome ( Table 1). PDF was first isolated from crustaceans, earning its name because of its pigment celldispersing activity (Hauser et al., 2010). Similar peptides were later found in insects (Rao et al., 1987), and these participate in the regulation of circadian rhythm (Renn et al., 1999). In D. melanogaster, PDF is involved in maintaining behavioral rhythms (Lin et al., 2004), while NPF has the function of regulating olfactory responses to food odor (Lee et al., 2020). In D. melanogaster, NPF is a 36 amino acid peptide whose C-terminal sequence is RVRFamide (Hauser et al., 2010), it is structurally related to the vertebrate neuropeptide Y family and affects food intake and feeding behavior (Wu et al., 2003). Interestingly, we found the highly expressed PDF receptor in H. hebetor (Figure 10 and Supplementary Figure 2) but could not identify the PDF neuropeptide. We speculate this absence arose because we failed to annotate it in the genomic and transcriptome data, rather than it not existing per se in H. hebetor. This would fit with the fact that we did not obtain the transcriptome of the neural tissue and so could have easily missed those genes represented by small fragments and low expression when we annotated them. In this context, we also could not find NPF in H. hebetor, however, it is a very conservative neuropeptide, and it also went unfound in D. alloeum, D. collaris, F. arisanus, M. cingulum, and M. demolitor in this study. This result is similar to findings reported by Chang et al. (2018), who also did not find NPF in F. arisanus, D. collaris, and D. alloeum (Chang et al., 2018). Thus, we speculate there are the large differences between the sequences in Braconidae and those in other species. This would explain the failed annotations made by using BLASTP and TBLASTN. Nonetheless, we did succeed in identifying the receptor of NPF, but at a low expression level (Figure 10, Supplementary Figure 2), the receptor of inotocin was never found. Therefore, we speculate that inotocin neuropeptides may have been lost from H. hebetor.
CCHa members and their receptors were also lost from all braconids analyzed in this study, except D. collaris and T. pretiosum ( Table 2). More than 10 years ago, Roller et al. (2008) discovered a novel neuropeptide CCHa in B. mori, and recently, two CCHa genes (CCHa1 and CCHa2) were identified in the tsetse fly Glossina morsitans . This type of neuropeptides contains two highly conserved cysteines and an amidated histidine residue at the C-terminus (Roller et al., 2008). The CCHa pre-prohormone is expressed in several small neurons in the CNS and in the midgut endocrine cells of B. mori larvae, but the biological function of the peptide remains unknown (Roller et al., 2008).
In contrast to N. vitripennis (Hauser et al., 2010) and P. puparum , we identified NTL in H. hebetor (Table 1). NTL is named for its function in promoting reproduction and modulating sexual activity. It was recently discovered and characterized in D. melanogaster, T. castaneum, and B. mori (Jiang et al., 2013). Being closely related to motif of tachykinin-related peptides (TKRPs), NTL is an arthropodspecific neuropeptide (Figure 12), that has a tachykinin-like signaling system (Jiang et al., 2013;Xu et al., 2016). In our study, both allatostatin CC (AstCC) and allatostatin CCC (AstCCC) were located to the same genomic scaffold. Among various insects, AstA inhibits the corpora allata to produce juvenile hormones or block muscle contraction in the gut (Hauser et al., 2010). The neuropeptides from insects' allatostatin family C are myoinhibitory in nature (Pandit et al., 2018). Allatostatins were discovered as rapid and reversible inhibitors of juvenile hormone synthesis in insects, though AstC is a cyclic neuropeptide that is structurally unrelated to AstA (Stay and Tobe, 2007). AstCC is a paralog of AstC, and AstCCC was also regarded as an AtsC in previous studies (Veenstra, 2009;Xu et al., 2016). AstCC and AstCCC are likely generated by gene duplication, and the similarity between the precursors and receptors of these genes suggests they have a common ancestor (Xu et al., 2016).
Most insect neuropeptides are peptide hormones involved in modulating physiology and behavior, such as the AKHs, EH, and ETH. In H. hebetor, we annotated an AKH ( Table 1). The AKHs are insect neuropeptides with 8-10 amino acid residues that are expressed in endocrine cells from the two neurohemal organs, corpora cardiaca, that are often fused and closely situated to the insect brain (Hauser et al., 2010;Li et al., 2016). AKH is a well-studied neuropeptide, thought to be similar to mammalian glucagon, which acts antagonistically to insulin by activating glycogen phosphorylase and mobilizing carbohydrates (Bacci et al., 2004). Some AKHs mobilize lipids and sugar mobilization from insect fat body during energy-consuming activities such as flight and intense locomotion (Li et al., 2008;Hauser et al., 2010). One axon branch from a pair of glucose-sensing neurons projects toward insulin-producing cells, thereby triggering the release of Drosophila insulin-like peptide 2 (dilp2), while the other branch extends to AKH-producing cells to inhibit the secretion of AKH, the analog of glucagon in flies (Oh et al., 2019). EH is known for its involvement in ecdysis behavior and is produced by a pair of brain neurons every time an insect molts in its lifetime (Zitnan et al., 2007;Li et al., 2008). In the CNS, EH induces the release of crustacean cardio-active peptide (CCAP), though it can also act peripherally on Inka cells in the epitracheal glands to induce the release of ETH (Hauser et al., 2010). A membrane-bound guanylyl cyclase was shown to function as an EH receptor in Inka cells (Chang et al., 2009). Respectively, located in different scaffolds in H. hebetor, ETH, EH and CCAP are the main players in the peptidergic circuit that controls ecdysis in insects . EH is expressed in CNS, ETH production occurs in endocrine cells, and the CCAP is expressed in neurons (Kim et al., 2006). In insects, development and metamorphosis are coordinated by a class of ecdysteroid hormones, of which 20hydroxyecdysone (20E) is the main one (Gilbert et al., 2002). Both the production and release of 20E from prothoracic glands are regulated by the brain-derived PTTH (Gilbert et al., 2002). We annotated PTTH in H. hebetor's genome. PTTH is secreted by the brain's neurosecretory cells and released into hemolymph at a specific time and in particular developmental stages (Yamanaka et al., 2010). The amino acid sequences of various insect PTTHs are not very well-conserved, but they are all believed to operate by forming homodimers (Hauser et al., 2010).
GPCRs mediate almost all physiological and behavioral processes, by delivering signals from a series of extracellular ligands, including neuropeptide hormones and biogenic amines, across the cell membrane to signal pathways that elicit specific responses (Audsley and Down, 2015). Some GPCRs can mediate extracellular ligand signals to intracellular signal transduction proteins (Ja et al., 2007), while others recognize various ligands, including odorants, photons, neurotransmitters, lipids, hormones, peptides, and other small molecules, enabling the regulation of intracellular responses when adapting to a changing environment (Zhu and Roth, 2014). In our study, we annotated 44 GPCRs. From the phylogenetic analysis, evidently some GPCRs in H. hebetor, P. puparum, N. vitripennis are clustered together, but some are not. This evolutionary divergence could be due to different taxonomic families that H. hebetor and the other two wasps belong to. If a neuropeptide gene is genuinely missing from a species, one should expect its receptor to be disabled and no longer affected by positive selection (Veenstra, 2019a). Hence, when both neuropeptide and its unique receptor are absent in a genome assembly, it is a reliable indication that the particular neuropeptide signaling system has been lost from the species in question (Veenstra, 2019a). As mentioned above, inotocin may fit this situation.
In our research, we performed RNA-Seq at different developmental stages in different organs including the venom and salivary glands of wasps. This expression data set provides us with basic characteristic information on the neuropeptides and their receptors, which is required to functionally address their biological significance in ectoparasitoid wasps. As revealed by the expression profile analysis, ETH is expressed higher in the pupal than other developmental stages (Figure 9), which is likely related to its function. Both ETH and EH regulate the release of CCAP from central CCAP neurons, thus inhibiting pre-ecdysis (Wang L. et al., 2018). Accordingly, it should be expressed more at the pupal stage. In parallel, we found that the peak expression of ATR and SIF differ (Figure 11). SIFamide (SIFa) influences sexual behavior, feeding, pupal mortality and sleep regulation in holometabolous insects (Ayub et al., 2020), and AT has multiple functions, including the regulation of juvenile hormone synthesis, growth, development, and reproduction (Zhang et al., 2019). We speculate that this variation in the expression level of SIFa and AT receptor could drive from differences between individual insects such as in their body size.
To conclude, through this comprehensive study we gained some timely insight into the neuropeptide precursor and neuropeptide receptor genes harbored by H. hebetor. We anticipate that our data will provide some basic yet practically useful information applicable to future pesticide development and other research.

DATA AVAILABILITY STATEMENT
The transcriptome data has been deposited in the Gene Expression Omnibus (accession: PRJNA642006).

AUTHOR CONTRIBUTIONS
KY, SX, GX, XY, HY, FW, QF, QS, and GY conceived and designed the experimental plan, and KY performed the experiments. KY and SX analyzed and interpreted the sequence and experimental data. GX, XY, HY, FW, QF, QS, and GY contributed the reagents, materials, and analysis tools. KY and SX prepared all the figures and tables. GX, QF, QS, and GY authored or reviewed drafts of the paper. GY revised the final draft. All authors read and approved the final manuscript.

ACKNOWLEDGMENTS
We are grateful to the Program for Chinese Outstanding Talents in Agricultural Scientific Research of the Ministry of Agriculture and Rural Affairs of the People's Republic of China. The authors sincerely thank the Guangzhou Institute of Cereal Science Research for providing the H. hebetor insects and Professor Hongxia Hua at the Huazhong Agricultural University for providing us with P. interpunctella.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys. 2020.575655/full#supplementary-material Supplementary Figure 1 | Expression profiles of neuropeptide precursor genes in H. hebetor across its different tissues. Log2 FPKM (fragments per kilobase of transcript per million) values for neuropeptide precursor genes are presented by in colored bars, where darker red denotes higher expression values, and darker blue denotes lower expression values. VG, venom glands; SG, salivary glands; VGC, carcass of female adults without their venom gland; SGC, carcass of larvae without salivary glands.
Supplementary Figure 2 | Expression profiles of neuropeptide precursor receptor genes in H. Hebetor across its different tissues. Log2 FPKM (fragments per kilobase of transcript per million) values for neuropeptide precursor genes are presented by bar colors where darker red denotes higher expression values, and darker blue denotes lower expression values. VG, venom glands; SG, salivary glands; VGC, carcass of female adults without their venom gland; SGC, carcass of larvae without salivary glands.
Supplementary Figure 3 | The qPCR analysis of relative expression levels of neuropeptide genes and neuropeptide receptor genes in different tissues. VG, venom glands; SG, salivary glands; VGC, carcass of female adults without their venom gland; SGC, carcass of larvae without salivary glands. The pink dotted line represents the FPKM (fragments per kilobase of transcript per million) value, and the relative expression levels from the qPCR are depicted in represented by the light blue histogram.