Transcriptome analysis reveals salivary gland-specific neuropeptide signaling genes in the predatory stink bug, Picromerus lewisi

Predatory stink bugs derive from phytophagous stink bugs and evolved enhanced predation skills. Neuropeptides are a diverse class of ancient signaling molecules that regulate physiological processes and behavior in animals, including stink bugs. Neuropeptide evolution might be important for the development of predation because neuropeptides can be converted to venoms that impact prey. However, information on neuropeptide signaling genes in predatory stink bugs is lacking. In the present study, neuropeptide signaling genes of Picromerus lewisi, an important predatory stink bug and an effective biological agent, were comprehensively identified by transcriptome analysis, with a total of 59 neuropeptide precursor genes and 58 potential neuropeptide receptor genes found. In addition, several neuropeptides and their receptors enriched in salivary glands of P. lewisi were identified. The present study and subsequent functional research contribute to an in-depth understanding of the biology and behavior of the predatory bugs and can provide basic information for the development of better pest management strategies, possibly including neuropeptide receptors as insecticide targets and salivary gland derived venom toxins as novel killing moleculars.


Introduction
Neuropeptides are a large diverse class of signaling molecules with key roles in insect physiology and behavior (Schoofs et al., 2017).Neuropeptides mediate their effects mainly through G protein-coupled receptors (GPCRs), receptor guanylyl cyclases (RGCs), and receptor tyrosine kinases (RTKs) (Caers et al., 2012).Given their essential regulatory functions and high specificity, neuropeptide signaling systems have been considered promising targets for "green" pest control (Verlinden et al., 2014;Audsley and Down, 2015).In order to develop a "green" insecticide based off this system, scientists should seek out unique neuropeptide signaling pathways that are absent in beneficial predatory insects.To this end, recent improvements in genome, transcriptome and proteome analysis have led to the discovery of neuropeptides and their receptors in a number of insects, including several beneficial species, which will provide valuable information for the development of novel insecticides with high selectivity.
The predatory stink bug Picromerus lewisi Scott (Hemiptera: Pentatomidae) is widely distributed in China and other Asian regions (Lin et al., 2000).It has been selected as an excellent biological control agent for a wide range of agricultural and forest insect pests, such as Lepidoptera larvae (Mu et al., 2022).Currently, P. lewisi has been successfully mass produced by natural enemy factories in China and has been shown to have good control of many important pests, such as Spodoptera frugiperda (Wang et al., 2019b;Tang et al., 2019).Predatory stink bugs use their salivary venom to paralyze prey and initiate extra-oral digestion.The salivary glands with venom components (also called venom glands) of predatory stink bugs are likely to be derived from the salivary glands of related non-venomous stink bug species (Cohen, 1990;Cohen, 1995;Walker et al., 2016).Insect salivary systems are influenced by neuropeptides and neurohormones secreted from endocrine organs and various neuronal cells (Spit et al., 2012).For example, in the cockroach, Periplanea americana, the SMYamide neuropeptide gene is specifically expressed in the neurons innervating the salivary glands and functions as a hormone to generate action potentials during feeding (Veenstra, 2020).In addition, insect salivary glands also have endocrine functions, producing neuropeptides and neurohormones to regulate other tissues and organs (Li et al., 2022b;Titos et al., 2023).A salivary gland-derived peptide, Sgsf, has been identified in Drosophila as an endocrine factor secreted into the hemolymph that systemically regulates larval growth (Li et al., 2022b).What's more, there are a few interesting cases where invertebrate neuropeptides have been recruited into salivary glands as venom toxins to affect prey, such as TKs from the cephalopods Eledone moschata and Octopus vulgaris (Champagne and Ribeiro, 1994;Kanda et al., 2003).
Given the state of the field described above, our fundamental hypotheses are that: 1) Insect neuropeptides might be repurposed into novel insecticides; 2) finding a neuropeptide pathway that is different or absent in beneficial predators, could lead to insecticides that do not hurt beneficial predators; 3) natural enemy evolution might repurpose neuropeptides as venom toxins, causing them to be expressed in salivary glands.To begin testing these hypotheses, we identified neuropeptides and their receptors in P. lewisi using transcriptome data analysis.In addition, the expression patterns of these genes in heads, guts and salivary glands were analyzed.The fundamental genetic information of the neuropeptide signaling system will be helpful to better understand the biology and ecology of the predatory bugs and to develop better pest management strategies utilizing bothe natural enemies and novel synthetic chemistries.

De novo assembly of the P. lewisi transcriptome
In our previous study, full-length transcriptome and RNA-seq transcriptome analysis of P. lewisi were jointly performed (Li et al., 2022a).However, no more than 20 neuropeptide precursor transcripts were found in the full-length transcriptome, which might be because a number of short precursor transcripts were polished during sequence clustering of raw reads obtained from Iso-seq sequencing.Herein, for the identification of neuropeptide signaling genes of P. lewisi, we performed de novo re-assembly of the RNA-seq transcriptome data that were generated from the previous study (Li et al., 2022a).Briefly, 100 fifth-instar nymphs starved for about 6 h were sampled for total RNA extraction from salivary glands (SG), guts (G), antennae (A), legs (L), and heads without antennae and salivary glands (H), with three independent biological replicates for each tissue/part sample.RNA samples were extracted using TRIzol reagent (Invitrogen, USA) and about 3 μg RNA per sample was used as the template for cDNA synthesis and Illumina sequencing using Ilumina NovaSeq 6000 (Illumina, USA).A total of fifteen libraries with approximately 20 million sequence clean reads each were generated and deposited in the National Center for Biotechnology Information (NCBI) Short Read Archive (SRA) under accession numbers of SRR20681617~SRR20681631 (Wang et al., 2019a).All libraries were pooled and then assembled into unigenes using the Trinity software (v2.6.6), with min_kmer_cov set to 2 by default and all other parameters set default.A Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis was performed to assess the completeness of the assembly (Simão et al., 2015).All unigenes were de novo annotated against Nr/Nt (NCBI non-redundant protein/nucleotide sequences), Pfam (Protein family), KOG/ COG (Clusters of Orthologous Groups of proteins), Swiss-Prot (A manually annotated and reviewed protein sequence database), KEGG (Kyoto Encyclopedia of Genes and Genomes), and GO (Gene Ontology) databases.

Differential expression analysis
Differential expression analysis of transcripts from SG, G and H with three biological replicates each was performed using the DESeq2 R package (1.20.0).The Fragments per kilobase of transcript per million fragments mapped (FPKM) value was calculated to estimate the expression level of each transcript from each library.The Benjamini-Hochberg approach was used to adjust p values (P_adj) to control the false discovery rate.Genes with P_adj values <0.05 and | log2 (fold change) | > 1 were designated as differentially expressed genes (DEGs).Heat maps showing gene expression profiles were constructed based on Log10 (FPKM +1) values.
To validate the transcriptome data, qRT-PCR analysis of eighteen genes had been performed (Li et al., 2022a).Briefly, eight cytochrome P450 monooxygenase (CYP) genes, six carboxyl/choline esterase (CCE) genes, and four glutathione S-transferase (GST) genes were selected, and the EF1A gene was used as the candidate reference gene.The RNA samples for qRT-PCR were the same as those for RNA-seq sequencing.Significant consistency was found between the expression profiles obtained by qRT-PCR and RNA-Seq (Li et al., 2022a).

Peptide structural prediction
A well-established workflow was used to predict the potential active peptides of P. lewisi.Briefly, the presence of signal peptides was predicted using the online program SignalP 5.0 (https://services.healthtech.dtu.dk/services/SignalP-5.0/).Prohormone cleavage sites and post-translational modifications were identified based on the information presented in Veenstra (Veenstra, 2000) and/or by homology to known arthropod peptides.

Sequence alignment and phylogenetic analysis
Multiple alignments of amino acid sequences were performed using the online program MAFFT version 7 (https://mafft.cbrc.jp/alignment/server/, "G-INS-1" progressive method setting).For the alignments of neuropeptides, the putative active peptides or the amino acid sequences removing the putative signal peptides were adopted.For the alignments of neuropeptide receptors, full-length sequences were used when available, or partial sequences were used.Amino acid identity was subsequently determined and the alignment result was visualized using GeneDoc version 2 (Nicholas et al., 1997).Sequence logos of the aligned peptide sequences were generated using the online program WebLogo (http://weblogo.berkeley.edu/logo.cgi)(Crooks et al., 2004).
Maximum likelihood (ML) phylogenetic analysis was conducted using IQ-TREE (Minh et al., 2020).In a first run, the ModelFinder function was employed to determine the best-fit model using Bayesian Information Criterion (Darriba et al., 2019).In a second run, ML phylogenetic tree was constructed using a ultrafast (UF) bootstrap test with 3000 replicates and default settings to reduce overestimation of bootstrap support (-bnni) (Hoang et al., 2017).Phylogenetic trees were visualized using the Interactive Tree Of Life (iTOL) web server (https:// itol.embl.de/)(Letunić and Bork, 2021).

Cloning neuropeptide transcripts
Two species-specific neuropeptide transcripts, Crustacean cardioactive peptide (CCAP) and Orcokinin B (OKB), were amplified by PCR and sequenced to verify their reliability.Primers were designed to amplify the coding sequences (Supplementary Table S1).Five fifth-instar nymphs of P. lewisi used for cloning were supplied by the Fenggang County Natural Enemy Breeding Center of the Guizhou Tobacco Company Zunyi Branch, Zunyi, Guizhou Province, China.Total RNA was isolated using TRIzol reagent (Invitrogen, United States).RNA quality and quantity were determined using the RNA Nano 6000 Assay Kit in the Bioanalyzer 2,100 system (Agilent Technologies, United States).Total RNA was treated with DNase I (Invitrogen, United States) to remove any residual genomic DNA.cDNA were synthesized from 1 μg of the total RNA using a SMARTer PCR cDNA Synthesis Kit (Takara Bio United States, Inc, United States).Neuropeptide precursor transcripts were amplified in a 50-uL reaction system using 2 × Phanta Flash Master Mix (Vazyme, China).The PCR procedure was set as follows: pre-denaturation at 98 °C for 30 s, 35 cycles of 98 °C for 10 s, 55 °C for 5 s and 72 °C for 5 s, and 3 Results

Transcript assembly
The raw sequences of 15 biosamples from five tissues/parts with three biological repeats each were re-assembled using the de novo assembly procedure, resulting in a total of 62,183 unigenes with a mean length of 1,203 bp, and an N50 length of 2,040 bp (Supplementary Table S2).BUSCO analysis showed a high degree of completeness (87.0%) in the assembly (Supplementary Figure S1).41.9% of the assembled unigenes were annotated in at least one database and 64.1% of the annotated unigenes were best matched to their H. halys homologs (Supplementary Table S2).

Identification of neuropeptide precursors in P. lewisi
A total of 59 neuropeptide precursors were identified in the P. lewisi transcriptome, including several novel neuropeptides, CNP, PTH, PaOGS36577, and RFLa (Table 1; Supplementary Data S1).Glycoprotein hormone beta 5 (GPB5), allatostatin C (AST-C), and trissin were not found in the P. lewisi transcriptome.Although AST-C was not identified, potential homologs of AST-CC and AST-CCC were found in P. lewisi and other hemipterans (Table 2).GPB5 was identified from H. halys and R. prolixus but still not found in N. viridula and L. hesperus (Table 2).Adipokinetic hormone (AKH) and sulfakinin (SK) were not found in our custom transcriptome but were identified in NCBI SRA data (SRR10134979).50 neuropeptide transcripts have full-length sequences and the remaining nine nonfull-length sequences include AKH, ALP1, ALP2, argininevasopressin-like peptide (AVLP), natalisin (NTL), OKA, OKB,   PTH and SK (Supplementary Data S1).Most of the predicted P. lewisi active neuropeptides are identical or highly similar to their homologous peptides from H. halys or other bugs (Table 1), whereas a few peptides have low similarity to their homologs, such as OKB, CCAP, prothoracicotropic hormone (PTTH) and neuroparsin A5 (NPA5) (with no more than 70% identity) (Figure 1).To increase the confidence in the putative P. lewisi-specific neuropeptide precursors, the coding sequences of the OKB and CCAP transcripts were targeted for RT-PCR amplification.Amplimers of the expected sizes were obtained for these two transcripts (730 bp and 195 bp for OKB and CCAP, respectively; Supplementary Figure S1).All of the cloned products were found to have 100% nucleotide identity with the in silico sequences from the transcriptomic data.The P. lewisi OKB precursor can encode eight potential active peptides that share the DXI/LGGG consensus sequence (Figure 1A).The H. halys OKB precursor can also encode eight distinct active peptides.Only one of the P. lewisi OKB active peptides shares the identical sequence (NLDTIGGGHLV) with an H. halys OKB peptide, whereas other P. lewisi OKB peptides had 45%~90% identity with their H. halys orthologues (Figures 1A~1B).
CCAP is known as the most conserved arthropod neuropeptide, with an identical amino acid sequence, PFCNAFTGCamide, found in all insects examined.Unexpectedly, the mature CCAP peptide predicted from P. lewisi (SFCNSFTGC) has two variant residues and lacks a C-terminal glycine residue (Figure 1C).The same CCAP peptide sequence was also found in P. lewisi RNA-seq SRA reads released in NCBI (e.g., SRR10134979.5654305.2).
The P. lewisi PTTH transcript has features typical of other insect PTTHs, encoding a putative active peptide containing seven cystine residues that form inter-and intra-chain disulfide bonds.PTTH has also been found in another heteropteran species, L. hesperus, but it is absent in a number of heteropteran species such as H. halys, N. viridula, and R. prolixus (Table 2).The P. lewisi PTTH peptide shares 39% amino acid identity with that of L. hesperus (Figure 1D).
A total of seven P. lewisi NPA transcripts were found (Supplementary Data S1).The P. lewisi NPA precursors can encode potential active peptides with a length of approximately 80 amino acid residues and 12~14 cystine residues forming 6~7 intrachain disulfide bridges (Figure 1E).Sequence alignment of the NPA peptides showed that Picle_NPA5 shares a relatively low identity (68%) with its H. halys orthologue, whereas other P. lewisi NPA peptides had 71%~96% identity with their H. halys orthologs (Figure 1E).

Identification of salivary gland-specific neuropeptides in P. lewisi
A heatmap based on FPKM values of 57 neuropeptide genes (excluding AKH and SK, which were identified from the NCBI SRA database) in heads, salivary glands, and guts of P. lewisi, is shown in Table 1.A total of 22 neuropeptide genes were found to be expressed in salivary glands using a criterion of more than 1 FPKM in at least one repeat (Supplementary Table S3).Among them, six neuropeptides, OKB, CCAP, OKA, diuretic hormone 31 (DH31), ALP2, and ecdysis triggering hormone (ETH), showed FPKM values higher than ten in at least one repeat (Supplementary Table S3).Compared to heads, only two neuropeptide genes (OKB and CCAP) were significantly more highly expressed in salivary glands, with log2 (fold) >1 and P_adj <0.05 (Table 1).OKB had the highest FPKM value (3228.1) in the salivary glands among all the identified neuropeptide genes and its relative expression level in the salivary glands compared to heads was 190.6-fold.CCAP was the second most highly expressed neuropeptide gene in the salivary glands (FPKM value = 109.3)and it was also more highly expressed in the salivary glands than heads, with an expression fold of 134.9.OKB, CNMa, AST-CC and CCHamide 2 (CCHa2) were significantly upregulated in guts than heads of P. lewisi (Table 1).

Identification of salivary gland-specific neuropeptide receptor genes in P. lewisi
A heatmap based on FPKM values of 57 neuropeptide receptor genes (excluding AKHR, which was identified from the trinity sequences rather than the unigene sequences of our custom transcriptome) in heads, guts and salivary glands is shown in  Frontiers in Physiology frontiersin.orgTable 3.A total of 15 out of 57 neuropeptide receptors were expressed in P. lewisi salivary glands with the FPKM values higher than one in at least one repeat, among which two receptor genes (EHR and NPLPR) presented the FPKM value higher than 10 in one repeat (Supplementary Table S4).Compared to heads, receptors for EH, NPLP, CAPA, GPA2/ GPB5, Ele, DH44, and ILP, and one orphan RTK receptor were more highly expressed in the salivary glands than heads, with log 2 (fold) >1 and P_adj <0.05 (Table 3).Receptors for CAPA, NPLP, and GPA2/GPB5, and one orphan RTK receptor were significantly more highly expressed in guts than heads (Table 3).

Discussion
In the present study, neuropeptide signaling genes were systematically identified in the predatory stink bug P. lewisi, with a total of 59 neuropeptide precursor genes and 58 potential neuropeptide receptor genes found.The number of neuropeptide precursors and their receptors identified in P. lewisi is similar to that found in H. halys and higher than those in other heteropterans such as N. viridula, L. hesperus, and R. prolixus (Christie et al., 2016;Lavore et al., 2018;Hull et al., 2021).In P. lewisi, nearly all of the neuropeptide signaling genes identified in other heteropterans were found, with the exception of GPB5 and RYaR.GPB5 and RYaR were not detected in the P. lewisi transcriptome, likely due to sample bias inherent to transcriptomes as opposed to genomes.
Tissue-specific expression profiles of neuropeptides and their receptors in P. lewisi provide basic information for in-depth studies of the biology and ecology of this important biological control species.The majority of neuropeptide precursor genes and their receptor genes in P. lewisi were found to be predominantly expressed in heads, indicating that most of neuropeptide signaling systems may act mainly as neuroregulators in the predatory bugs and/or derive from tissues in the head.A few neuropeptide signaling genes were also expressed in salivary glands and guts.These may play critical roles in the regulation of feeding and digestion.Neuropeptides in insect guts are involved in the regulation of feeding processes, including food choice, appetite, digestion, metabolism and excretion (Abou El Asrar et al., 2020).DH31, CCHa-1, CCHa-2, NPF, and CNMa are well-known insect neuropeptides associated with the gut-brain axis, regulating appetite, sleep and even courtship (Sano et al., 2015;Kim et al.,   Frontiers in Physiology frontiersin.org2021; Lin et al., 2022;Malita et al., 2022;Titos et al., 2023).In P. lewisi, three neuropeptides (CNMa, ASTCC and CCHa2) and five receptor genes were significantly upregulated in the guts when compared to heads, suggesting that they may have functions in feeding regulation.
In P. lewisi, three neuropeptide genes (CCAP and OKB) and eight receptors were enriched in salivary glands.In all these cases, the physiological functions of only a few neuropeptides and their receptors in insect salivary glands have been determined.More research is needed on this topic to describe specific functions.
Several insect neuropeptides have been implicated in the neural control of salivary production or secretion, such as FMRFa from the blow fly Calliphora vomitoria (Duve et al., 1992) and the kissing bug R. prolixus (Ons, 2017).In R. prolixus, AST-B, Pro, CCAP, and CCAPR have also been detected in processes innervating salivary glands, indicating their involvement in the hormonal control of salivary production or secretion (Ons, 2017).In ticks, several neuropeptides have been identified to be expressed in innervations of salivary glands, such as AST-B, DH31, Ele, ELeR, NPLP1, OKA, PDF, SIFa, and SIFaR (Sterkel et al., 2011;Simo et al., 2012;Ladislav et al., 2015;Kim et al., 2018;Vancová et al., 2019;FIGURE 2 Phylogenetic tree of neuropeptide GPCR receptors from P. lewisi (gene IDs beginning "Picle"), A. lucorum (Apolu) and N. lugens (Nl).The ultrafast (UF) bootstrap value with more than 99% was marked in the tree.The P. lewisi gene names are marked with a blur color.Guerrib et al., 2023).CCAP, OKA, DH31, and NPLPR, which were highly expressed in P. lewisi, indicated their possible functions in neural control of the salivary system, as consistent with previous studies, however, expression of OKB and EHR in venom gland or salivary glands has not been reported in other insects so far.
Limited studies provide evidence for the functions of neuropeptides derived from invertebrate salivary glands as endocrine factors that regulate other tissues and organs.In addition, some of these salivary peptides might be venom toxins which alter the physiology of another species.A attractive study demonstrated that a peptide Sgsf expressed in Drosophila salivary glands can be secreted into the hemolymph and regulate Dilp2 secretion in the brain (Li et al., 2022b).Nässel et al. (2019) reviewed that TKs produced by salivary glands of mosquitos and cephalopods have been identified as exogenous vasodilators or paralyzing agents that can be delivered to prey.Determining whether neuropeptides expressed in the salivary glands of P. lewisi could be endogenous endocrine factors or exogenous venom peptides is important because venom peptides might be utilized as novel insecticidal peptides and/or signaling neuropeptides that their receptors might become useful targets for insecticides.
Our study presented species-and tissue-specific expression patterns of the neuropeptide signaling system in the predatory bugs, which will be used to generate testable functional genetic hypotheses in future studies.The most noteworthy neuropeptide is P. lewisi OKB, which was the most highly expressed neuropeptide in P. lewisi salivary glands.The expression pattern of OKB in the nervous system and intestine has also been reported in D. melanogaster, Bombyx mori, T. castaneum, R. prolixus, and Blattella germanica (Chen et al., 2015;Jiang et al., 2015;Ons et al., 2015;Wulff et al., 2017;Wang et al., 2019a).However, expression of OKB in salivary glands has not been reported in other insect species to date.OKs were first discovered with the myotropic activity in the crayfish Orconectes limosus (Stangier et al., 1992).One study suggested a role of OKB awakening behavior in T. castaneum (Ladislav et al., 2015).There is plenty of evidence for the recruitment of neuropeptides into animal venoms (Nässel et al., 2019;Sachkova et al., 2020;Goudarzi et al., 2023).Therefore, we would like to determine whether OKBs could be recruited as venom peptides in the salivary glands of P. lewisi and injected into prey to regulate myotropic activity or behavior.From an evolutionary perspective, the high genetic diversity of OKB peptide sequences possibly supports the repurposing of neuropeptides into venom peptides.A well-known example is TKs, another highly diverse and pleiotropic neuropeptide class, which have been convergently recruited into the venom or salivary glands of venomous invertebrates to affect prey (Nässel et al., 2019).This interesting evolutionary path of neuropeptide recruitment for novel toxins has also been recently revealed in the sea anemone Nematostella vectensis (ShK-like peptides) and the caterpillar Acharia stimulea (RF-amide peptides) (Sachkova et al., 2020;Goudarzi et al., 2023).Nevertheless, we cannot exclude the possibility that OKB peptides are endogenous endocrine factors.The OK receptor has not yet been identified in any species, therefore efforts to identify the OK receptor in insects will be crucial to understanding the role of OKBs in the salivary glands of P. lewisi.
The most unexpected result in the present study was finding a very atypical CCAP in P. lewisi.To our knowledge, CCAP is the identical neuropeptide (PFCNAFTGC-NH 2 ) in all examined insect species.CCAP from P. lewisi (SFCNSFTGC) identified in the present study has two variant amino acid residues and no amidation at the C-terminus.Thus our findings indicate a highly novel primary structure for this peptide which we will seek to confirm through further studies utilizing liquid chromatography tandem mass spectrometry analysis.CCAP is mainly expressed in the central nervous system and/or in the intestine in most of the determined insects, such as D. melanogaster and R. prolixus (Lee et al., 2011;Shi et al., 2019).In R. prolixus, CCAP and its receptor have also been detected to be present in salivary glands (Lee et al., 2011;Lee and Lange, 2011;Lee et al., 2013).However, salivary gland-specific expression of CCAP in P. lewisi has not been reported in other insects to date.CCAP plays a crucial role in numerous biological and physiological processes in insects, mainly including the regulation of heart contraction, ecdysis, and feeding (Sakai et al., 2006;Estevez-Lao et al., 2013;Lee et al., 2013;Marco et al., 2018;Shi et al., 2019;Shen et al., 2021;Verbakel et al., 2021;Shi et al., 2022a), but its potential physiological function in the salivary glands of R. prolixus is still unknown.In P. lewisi, three putative CCAP receptors were identified, which were highly identical to their homologs (>90% identities) and highly expressed in heads, indicating that salivary gland-derived CCAP of P. lewisi could be considered as an endogenous endocrine factor.It is worth mentioning that the presence of the atypical CCAP in P. lewisi venom as a toxin cannot be ruled out given an intriguing example of a CCAPrelated peptide discovered in the venom of Conus villepinii and having the activity of decreasing the heart frequency in Drosophila larvae (Möller et al., 2010).The insect neuropeptidergic system has been considered as an ideal target for the development of greener pest control strategies.Comparative genomics and transcriptomics provide useful information for appropriate design strategies to develop target-specific insecticidal molecules that could successfully control pests while protecting beneficial species.To date, insecticidal activity and biosafety have been demonstrated for a few insect neuropeptides and their analogues, such as kinins, proctolin, CAPA, and TKs (Shi et al., 2022b).For example, insect kinins and their analogues were determined to exhibit high efficacy against aphids, however, they showed safety to an aphid predator, the common green lacewing Chrysoperla carnea, based on the transcriptome analysis information that insect kinins were not found in C. carnea (Shi et al., 2022b).Based on the comparative analysis of gene sets of known neuropeptides and their receptors between P. lewisi and other heteropterans, P. lewisi harbors almost all kinds of neuropeptide signaling system identified in other heteropteran species.Because of the conservation of this signaling system, it seems difficult to develop green pest control strategies based on neuropeptide systems putatively lost in benefical predatory bugs.Although sequence alignment of mature peptides between P. lewisi and its close herbivorous bug species showed several neuropeptides like OKB, CCAP, and PTTH, which are highly diverse in bugs, these neuropeptide signaling classes could be potential targets for the development of highly selective insecticidal agents needs to be further determined.

Conclusion
In the present study, a total of 59 neuropeptide precursors and 58 potential neuropeptide receptor genes were identified through transcriptomic analysis.The present study also revealed a set of neuropeptides and their receptors that were enriched in the salivary glands of P. lewisi, providing basic information for in-depth study on repurposing neuropeptides and their receptors into insecticides and targets.Li et al. 10.3389/fphys.2023.1270751 halys] 100% a Neuropeptides from P. lewisi are classified into those encoding a single and short active peptide (S1, the length of an encoded peptide is less than 25 amino acid residues), multiple and short active peptides (Sx, x as the number of encoded peptides, the number in parentheses indicates the distinct peptide numbers), and long active peptides (L1~L2, the length of an encoded peptide is more than 25 amino acid residues).b Identity values of predicted active peptides of P. lewisi and their homologous peptides were calculated based on pairwise alignemnts or multiple alignments.c Expression levels of two neuropeptide genes (AKH, and SK) are missing because they were identified from the public SRA, database.d The potential short active peptides of CNP, and NPLP1 of P. lewisi are not clear.*indicates that the expression level of a gene was significantly higher in G or SG than that in H, with log2 (fold) >1 and P_adj<0.05.The heatmap scale among heads (H), guts (G) and salivary glands (SG) was based on Log10 (FPKM+1) values:

FIGURE 1
FIGURE 1 Multiple sequence alignment of four neuropeptides from P. lewisi and other bugs.(A,B) OKB active peptides from P. lewisi and H. halys.(C) CCAP precursors cutting signal peptides from P. lewisi, H. halys and N. viridula.(D) PTTH precursors cutting signal peptides from P. lewisi and L. hesperus.(E) NPA precursors cutting signal peptides from P. lewisi and H. halys.Full conservative residues are shaded in black background.Predicted convertase cleavage sites are shaded in red.Glycine residues are shaded in pink.Conserved cysteine residues are shaded in yellow.Identity values between the P. lewisi precursor and its ortholog are shown at the end.

a
Identity values of neuropeptide receptor genes of P. lewisi and their homologous genes were calculated based on pairwise alignments or multiple alignments.b Expression level of AKHR, is missing because it was identified from the public SRA, database.c Not full length.*indicates that the expression level of one gene in G or VG was significantly higher than that in H, with log2 (fold) >1 and P_adj<0.05.The heatmap scale among heads (H), guts (G) and salivary glands (SG) was based on Log10 (FPKM+1) values:
(Continued on following page)Frontiers in Physiology frontiersin.orgtermination at 72 °C for 1 min.The PCR products were cut from the gel, purified and sequenced by Beijing Tsingke Biotech Co., Ltd.

TABLE 2
Neuropeptide precursors identified in P. lewisi and four other heteropteran species.

TABLE 2 (
Continued) Neuropeptide precursors identified in P. lewisi and four other heteropteran species.