Genome-wide identification and expression profiling of serine proteases and homologs in the diamondback moth, Plutella xylostella (L.)

Identification of seminal fluid proteins and reproductive function of trypsin-1 in male Plutella xylostella

Insect seminal fluid proteins (SFPs) are primary factors affecting physiology and behavior in both sexes, making them valuable targets for pest control. However, SFPs have not been fully characterized in the Plutella xylostella, a global pest that attacks cruciferous crops. Here, 75 putative SFPs were identified in P. xylostella, compared to 10 orthologs in Drosophila melanogaster, 10 in Nilaparvata lugens, 5 in Apis mellifera, and 43 in Heliconius melpomene. Analyses of Ka/Ks suggested that SFPs had high evolution rates. Proteases (22/75, 29.3 %) accounted for the highest proportion of P. xylostella SFPs, including 16 trypsins. The phylogenetic analysis showed that most trypsins from P. xylostella and H. melpomene belonged to the same cluster. SFP04 (trypsin-1) was orthologous to the SFP ADJ58550.1 in H. melpomene. PxTry1 was specifically expressed in adult males and their accessory glands but was also detected in females after mating. A CRISPR/Cas9-induced PxTry1 homozygous mutant strain with a 22-base pair nucleotides insertion was generated. PxTry1 deletion resulted in swollen testes, smaller spermatophores, and abnormal sperm, thus reducing the P. xylostella egg production and hatching. These results clarify the role of insect SFPs in evolution and reproduction and identify a promising target for pest control based on genetic regulation.

Trypsin-encoding gene function of efficient star polycation nanomaterial-mediated dsRNA feeding delivery system of Grapholita molesta

BACKGROUND

Grapholita molesta is an important and harmful fruit pest worldwide, with widespread feeding hosts. Trypsin, an indispensable hydrolytic digestive protease in the insect gut, is crucial in digestion, growth and development. We analyzed the characteristics of the trypsin-encoding genes, screened for the optimal dose of RNAi mediated by nanocarriers, and investigated various indices of larval growth and development of G. molesta.

RESULTS

Gut content (GC) and RNase A degraded double-stranded RNA (dsRNA), with a faster degradation rate at higher concentrations. Star polycation (SPc) nanomaterials protected dsGFP from degradation by anion-cation binding and did not migrate through agarose gel. The key conserved motifs of the trypsin-encoding genes were similar, exhibiting high homology with those in other lepidopteran insects. An interference efficiency of ≈70% was achieved with SPc nanomaterial-mediated RNA interference with 0.05 μg dsRNA. The efficiency of continuous interference was stable. Trypsin activity, body weight of 8-day-old larvae, pupal weight and emergence rate were significantly reduced, and the larval stage was significantly prolonged.

CONCLUSION

The investigated trypsin gene is a key target gene in the growth and development of G. molesta. We investigated the efficiency and convenience of feeding SPc nanomaterials in a functional study of insects. Our results provide valuable data for the development of efficient trypsin-targeting pesticides. © 2024 Society of Chemical Industry.

A sperm-activating trypsin-like protease from the male reproductive tract of Spodoptera litura: Proteomic identification, sequence characterization, gene expression profile, RNAi and the effects of ionizing radiation

Like other lepidopteran insects, males of the tobacco cutworm moth, Spodoptera litura produce two kinds of spermatozoa, eupyrene (nucleate) and apyrene (anucleate) sperm. Formed in the testis, both kinds of sperm are released into the male reproductive tract in an immature form and are stored in the duplex region of the tract. Neither type of sperm is motile at this stage. When stored apyrene sperm from the duplex are treated in vitro with an extract of the prostatic region of the male tract, or with mammalian trypsin, they become motile; activation is greater and achieved more rapidly with increasing concentration of extract or enzyme. The activating effect of prostatic extract is blocked by soybean trypsin inhibitor (SBTI), also in a dose-dependent way. These results suggest that the normal sperm-activating process is due to an endogenous trypsin-like protease produced in the prostatic region. Proteomic analysis of S. litura prostatic extracts revealed a Trypsin-Like Serine Protease, TLSP, molecular weight 27 kDa, whose 199-residue amino acid sequence is identical to that of a predicted protein from the S. litura genome and is highly similar to predicted proteins encoded by genes in the genomes of several other noctuid moth species. Surprisingly, TLSP is only distantly related to Serine Protease 2 (initiatorin) of the silkmoth, Bombyx mori, the only identified lepidopteran protein so far shown to activate sperm. TLSP has features typical of secreted proteins, probably being synthesized as an inactive precursor zymogen, which is later activated by proteolytic cleavage. cDNA was synthesized from total RNA extracted from the prostatic region and was used to examine TLSP expression using qPCR. tlsp mRNA was expressed in both the prostatic region and the accessory glands of the male tract. Injection of TLSP-specific dsRNA into adult males caused a significant reduction after 24 h in tlsp mRNA levels in both locations. The number of eggs laid by females mated to adult males that were given TLSP dsRNA in 10 % honey solution, and the fertility (% hatched) of the eggs were reduced. Injecting pupae with TLSP dsRNA caused the later activation of apyrene sperm motility by adult male prostatic extracts to be significantly reduced compared to controls. Exposure of S. litura pupae to ionizing radiation significantly reduced expression of tlsp mRNA in the prostatic part and accessory gland of irradiated males in both the irradiated generation and also in their (unirradiated) F1 progeny. The implications of these findings for the use of the inherited sterility technique for the control of S. litura and other pest Lepidoptera are discussed.

The Set of Serine Peptidases of the Tenebrio molitor Beetle: Transcriptomic Analysis on Different Developmental Stages

Serine peptidases (SPs) of the chymotrypsin S1A subfamily are an extensive group of enzymes found in all animal organisms, including insects. Here, we provide analysis of SPs in the yellow mealworm Tenebrio molitor transcriptomes and genomes datasets and profile their expression patterns at various stages of ontogeny. A total of 269 SPs were identified, including 137 with conserved catalytic triad residues, while 125 others lacking conservation were proposed as non-active serine peptidase homologs (SPHs). Seven deduced sequences exhibit a complex domain organization with two or three peptidase units (domains), predicted both as active or non-active. The largest group of 84 SPs and 102 SPHs had no regulatory domains in the propeptide, and the majority of them were expressed only in the feeding life stages, larvae and adults, presumably playing an important role in digestion. The remaining 53 SPs and 23 SPHs had different regulatory domains, showed constitutive or upregulated expression at eggs or/and pupae stages, participating in regulation of various physiological processes. The majority of polypeptidases were mainly expressed at the pupal and adult stages. The data obtained expand our knowledge on SPs/SPHs and provide the basis for further studies of the functions of proteins from the S1A subfamily in T. molitor.

Distinct enzyme activities of serine protease p37k in silkworm midgut and molting fluid

Serine proteases possess various biological functions. The serine protease p37k exhibits gelatinolytic activity in the silkworm midgut and degrades cuticular proteins in the molting fluid. In this study, we analyzed the activity changes of recombinant p37k (re-p37k) and p37k in the midgut and molting fluid of Bombyx mori. Firstly, in vitro-expressed re-p37k was activated when a 22 kDa band was observed by western blot. Re-p37k exhibits strong gelatinolytic activity, with the highest activity observed at pH 7.0-9.0 and 45 °C. Compared to p37k in the midgut, re-p37k loses thermal stability but can be restored by midgut extract or ions. E64, AEBSF, and an inhibitor cocktail inhibited the hydrolytic activity of re-p37k on epidermal proteins but did not inhibit the gelatinolytic activity. Subsequently, zymography showed that the positions of gelatinolytic band produced by p37k in the midgut and molting fluid were different, 35 kDa and 40 kDa, respectively. Finally, when heated midgut extract was added to re-p37k or molting fluid, the gelatinolytic band shifted from 40 kDa to 35 kDa, and the proteolytic activity of p37k in the molting fluid was inhibited. Collectively, our results demonstrate that p37k exhibits different activities in various tissues, suggesting its distinct tissue-specific functions during insect metamorphosis.

MiR-2b-3p Downregulated PxTrypsin-9 Expression in the Larval Midgut to Decrease Cry1Ac Susceptibility of the Diamondback Moth, Plutella xylostella (L.)

Crystal (Cry) toxins, produced by Bacillus thuringiensis, are widely used as effective biological pesticides in agricultural production. However, insects always quickly evolve adaptations against Cry toxins within a few generations. In this study, we focused on the Cry1Ac protoxin activated by protease. Our results identified PxTrypsin-9 as a trypsin gene that plays a key role in Cry1Ac virulence in Plutella xylostella larvae. In addition, P. xylostella miR-2b-3p, a member of the micoRNA-2 (miR-2) family, was significantly upregulated by Cry1Ac protoxin and targeted to PxTrypsin-9 downregulated its expression. The mRNA level of PxTrypsin-9, regulated by miR-2b-3p, revealed an increased tolerance of P. xylostella larvae to Cry1Ac at the post-transcriptional level. Considering that miR-2b and trypsin genes are widely distributed in various pest species, our study provides the basis for further investigation of the roles of miRNAs in the regulation of the resistance to Cry1Ac and other insecticides.

Parasitoid Serpins Evolve Novel Functions to Manipulate Host Homeostasis

Parasitoids introduce various virulence factors when parasitism occurs, and some taxa generate teratocytes to manipulate the host immune system and metabolic homeostasis for the survival and development of their progeny. Host-parasitoid interactions are extremely diverse and complex, yet the evolutionary dynamics are still poorly understood. A category of serpin genes, named CvT-serpins, was discovered to be specifically expressed and secreted by the teratocytes of Cotesia vestalis, an endoparasitoid of the diamondback moth Plutella xylostella. Genomic and phylogenetic analysis indicated that the C. vestalis serpin genes are duplicated and most of them are clustered into 1 monophyletic clade. Intense positive selection was detected at the residues around the P1-P1' cleavage sites of the Cv-serpin reactive center loop domain. Functional analyses revealed that, in addition to the conserved function of melanization inhibition (CvT-serpins 1, 16, 18, and 21), CvT-serpins exhibited novel functions, i.e. bacteriostasis (CvT-serpins 3 and 5) and nutrient metabolism regulation (CvT-serpins 8 and 10). When the host-parasitoid system is challenged with foreign bacteria, CvT-serpins act as an immune regulator to reprogram the host immune system through sustained inhibition of host melanization while simultaneously functioning as immune effectors to compensate for this suppression. In addition, we provided evidence that CvT-serpin8 and 10 participate in the regulation of host trehalose and lipid levels by affecting genes involved in these metabolic pathways. These findings illustrate an exquisite tactic by which parasitoids win out in the parasite-host evolutionary arms race by manipulating host immune and nutrition homeostasis via adaptive gene evolution and neofunctionalization.

Differences of functionalized graphene materials on inducing chronic aquatic toxicity through the regulation of DNA damage, metabolism and oxidative stress in Daphnia magna

Graphene can be modified with functional groups when released into the environment. However, very little is known about molecular mechanisms of chronic aquatic toxicity induced by graphene nanomaterials with different surface functional groups. By using RNA sequencing, we investigated the toxic mechanisms of unfunctionalized graphene (u-G), carboxylated graphene (G-COOH), aminated graphene (G-NH₂), hydroxylated graphene (G-OH) and thiolated graphene (G-SH) to Daphnia magna during 21-day exposure. We revealed that alteration of ferritin transcription levels in the "mineral absorption" signaling pathway is a molecular initiating event leading to potential of oxidative stress in Daphnia magna by u-G, while toxic effects of four functionalized graphenes are related to several metabolic pathways including the "protein digestion and absorption" pathway and "carbohydrate digestion and absorption" pathway. The transcription and translation related pathways were inhibited by G-NH₂ and G-OH, which further affected the functions of proteins and normal life activities. Noticeably, detoxifications of graphene and its surface functional derivatives were promoted by increasing the gene expressions related to chitin and glucose metabolism as well as cuticle structure components. These findings demonstrate important mechanistic insights that can potentially be employed for safety assessment of graphene nanomaterials.

Clip-SP1 cleavage activates downstream prophenoloxidase activating protease (PAP) in Plutella xylostella

Melanization is a component of the humoral immune defense of insects and is induced by serine protease-mediated phenoloxidase (PO) catalysis. Prophenoloxidase (PPO) in the midgut of Plutella xylostella is activated by the CLIP domain serine protease (clip-SP) in response to Bacillus thuringiensis (Bt) infection, but the detailed signaling cascade following this activation is unknown. Here, we report that activation of clip-SP enhances PO activity in the P. xylostella midgut by cleaving three downstream PPO-activating proteases (PAPs). First, the expression level of clip-SP1 was increased in the midgut after Bt8010 infection of P. xylostella. Then, purified recombinant clip-SP1 was able to activate three PAPs - PAPa, PAPb and PAP3 - which in turn enhanced their PO activity in the hemolymph. Furthermore, clip-SP1 showed a dominant effect on PO activity compared to the individual PAPs. Our results indicate that Bt infection induces the expression of clip-SP1, which is upstream of a signaling cascade, to efficiently activate PO catalysis and mediate melanization in the midgut of P. xylostella. And it provides a basis for studying the complex PPO regulatory system in the midgut during Bt infection.

Identification and gene expression analysis of serine proteases and their homologs in the Asian corn borer Ostrinia furnacalis

Serine proteases (SPs) and their homologs (SPHs) are among the best-characterized gene families. They are involved in several physiological processes, including digestion, embryonic development and immunity. In the current study, a total of 177 SPs-related genes were characterized in the genome of Ostrinia furnacalis. The activation site of SPs/SPHs and enzyme specificity of SPs were identified, and the findings showed that most of the SPs analyzed possessed trypsin substrate specificity. Several SPs/SPHs with similar simple gene structures had tandem repeat-like distributions on the scaffold, indicated that gene expansion has occurred in this large family. Furthermore, we constructed 30 RNA sequencing libraries including four with developmental stage and four middle larval stage tissues to study the transcript levels of these genes. Differentially upregulated and downregulated genes were obtained via data analysis. More than one-quarter of the genes were specifically identified as highly expressed in the midgut in compared to the other three tissues evaluated. In the current study, the domain structure, gene location and phylogenetic relationship of genes in O. furnacalis were explored. Orthologous comparisons of SPs/SPHs between model insects and O. furnacalis indicated their possible functions. This information provides a basis for understanding the functional roles of this large family.

An ancient enzyme finds a new home: Prevalence and neofunctionalization of trypsin in marine phytoplankton

Trypsin is an ancient protease best known as a digestive enzyme in animals, and traditionally believed to be absent in plants and protists. However, our recent studies have revealed its wide presence and important roles in marine phytoplankton. Here, to gain a better understanding on the importance of trypsin in phytoplankton, we further surveyed the distribution, diversity, evolution and potential ecological roles of trypsin in global ocean phytoplankton. Our analysis indicated that trypsin is widely distributed both taxonomically and geographically in marine phytoplankton. Furthermore, by systematic comparative analyses we found that algal trypsin could be classified into two subfamilies (trypsin I and trypsin II) and exhibited highly duplicated and diversified during evolution. We also observed markedly different domain sequences and organizations between and within the subfamilies, suggesting potential neofunctionalization. Diatoms contain both subfamilies of trypsin, with higher numbers of genes and more environment-responsive expression of trypsin than other lineages. The duplication and subsequent neofunctionalization of the trypsin family may be important in diatoms for adapting to dynamical environmental conditions, contributing to diatoms' dominance in the coastal oceans. This work advances our knowledge on the distribution and neofunctionalization of this ancient enzyme and creates a new window of research on phytoplankton biology.

Molting-related proteases in the brown planthopper, Nilaparvata lugens

Digestion and absorption of old cuticles during insect molting are necessary for new cuticle formation, during which complicated enzyme catalysis is essential. To date, a few carboxypeptidases, aminopeptidases and serine proteases (mostly trypsins) connected with cuticle digestion, zymogen activation and histological differentiation during the ecdysis of lepidopteran, dipteran and hymenopteran insects have been identified. However, little is known about these proteins in hemimetabolous insects. In this study, we identified 33 candidate trypsin and trypsin-like homologs, 14 metallocarboxypeptidase and 32 aminopeptidase genes in the brown planthopper Nilaparvata lugens, a hemipteran rice pest. Among the proteins encoded by these genes, 9 trypsin-like proteases, 3 metallocarboxypeptidases and 1 aminopeptidase were selected as potential procuticle hydrolases by bioinformatics analysis and in vivo validation. RNA interference targeting these genes demonstrated that 3 trypsin-like proteases (NlTrypsin-8, NlTrypsin-29 and NlTrypsin-32) genes and 1 metallocarboxypeptidase (NlCpB) gene were found to be essential for ecdysis in N. lugens; specifically, gene silencing led to incomplete cuticle degradation and arrested ecdysis, causing lethal morphological phenotype acquisition. Spatiotemporal expression profiling by quantitative PCR and western blotting revealed their specific expression in the integument and their periodic expression during each stadium, with a peak before ecdysis and eclosion. Transmission electron microscopy demonstrated corresponding ultrastructural defects after RNAi targeting, with NlCpB-silenced specimens having the most undigested old procuticles. Immunohistochemical staining revealed that NlTrypsin-8, NlTrypsin-29 and NlCpB were predominantly located in the exuvial space. This research further adds to our understanding of proteases and its potential role in insect ecdysis.

Characterization and expression profiling of serine protease inhibitors in the yellow mealworm Tenebrio molitor

Serine protease inhibitors (SPIs) act in diverse biological processes in insects such as immunity, development, and digestion by preventing the unwanted proteolysis. So far, the repertoire of genes encoding SPIs has been identified from few insect species. In this study, 62 SPI genes were identified from the genome of the yellow mealworm, Tenebrio molitor. According to their modes of action, they were classified into three families, serpin (26), canonical SPI (31), and α-macroglobulins (A2M) (5). These SPIs feature eight domains including serpin, Kazal, TIL, Kunitz, WAP, Antistasin, pacifastin, and A2M. In total, 39 SPIs contain a single SPI domain, while the others encode at least two inhibitor units. Based on the amino acids in the cleaved reactive sites, the abilities of these SPIs to inhibit trypsin, chymotrypsin, or elastase-like enzymes are predicted. The expression profiling based on the RNA-seq data showed that these genes displayed stage-specific expression patterns during development, suggesting to us their significance in development. Some of the SPI genes were exclusively expressed in particular tissues such as hemocyte, fat body, gut, ovary, and testis, which may be involved in biological processes specific to the indicated tissues. These findings provide necessary information for further investigation of insect SPIs.

Identification and expression profiling of serine protease-related genes in Tenebrio molitor

In insects, serine proteases and serine protease homologs (SPs/SPHs) are involved in a variety of physiological processes including digestion, development, and immunity. Here, we identified 112 SP and 88 SPH genes in the genome of the yellow mealworm, Tenebrio molitor. Based on the features of domain structure, they were divided into "S" group containing single Tryp-SPc or Tryp-SPHc domain, "C" group containing 1-4 CLIP domain (CLIPA-D) and "M" group containing the CBD, CUB, EGF, Fz, Gd, LDLa, PAN, SEA, SR, Sushi, and TSP domains, and have 115, 48, and 37 gene members, respectively. According to the active sites in the catalytic triad, the putative trypsin, chymotrypsin, or elastase-like enzyme specificity of the identified SPs/SPHs were predicted. Phylogenetic and genomic location analyses revealed that gene duplication exists in the large amount of SPs/SPHs. Gene expression profiling using RNA-seq data along with real time reverse transcription-polymerase chain reaction analysis showed that most SP/SPH genes display life stage specific expression patterns, indicating their important roles in development. Many SP/SPH genes are specifically or highly expressed in the gut, salivary gland, fat body, hemocyte, ovary, and testis, suggesting that they participate in digestion, immunity, and reproduction. The findings lay the foundation for further functional characterization of SPs/SPHs in T. molitor.

RNA Sequencing Reveals the Potential Adaptation Mechanism to Different Hosts of Grapholita molesta

Grapholita molesta is an important fruit tree worldwide pest which feeds on hosts extensively and does serious harm. In this paper, the growth and development parameters and protease activities of G. molesta fed on different hosts were compared. Using Illumina RNA sequencing technology, 18 midgut samples from five different hosts (apple, pear, plum, peach and peach shoots) and artificial diet were sequenced and compared with the reference genome, resulting in 15269 genes and 2785 predicted new genes. From 15 comparative combinations, DEGs were found from 286 to 4187 in each group, with up-regulated genes from 107 to 2395 and down-regulated genes from 83 to 2665. KEGG pathway analysis showed that DEGs were associated with amino acid metabolism, starch and sucrose metabolism, carbohydrate metabolism, and hydrolase activity. A total of 31 co-expression gene modules of different hosts were identified by WGCNA. qRT-PCR showed that the expression pattern of the trypsin gene was consistent with RNA sequencing. In this study, growth and development parameters, protease activity, DEGs, enrichment analysis and qRT-PCR were combined to reveal the adaptation process to different hosts of G. molesta in many aspects. The results of this study provide a basis for further exploration of the molecular mechanism of host adaptation of G. molesta.

Differential Transcriptional Responses in Two Old World Bemisia tabaci Cryptic Species Post Acquisition of Old and New World Begomoviruses

Begomoviruses are transmitted by several cryptic species of the sweetpotato whitefly, Bemisia tabaci (Gennadius), in a persistent and circulative manner. Upon virus acquisition and circulative translocation within the whitefly, a multitude of molecular interactions occur. This study investigated the differentially expressed transcript profiles associated with the acquisition of the Old World monopartite begomovirus, tomato yellow leaf curl virus (TYLCV), and two New World bipartite begomoviruses, sida golden mosaic virus (SiGMV) and cucurbit leaf crumple virus (CuLCrV), in two invasive B. tabaci cryptic species, Middle East-Asia Minor 1 (MEAM1) and Mediterranean (MED). A total of 881 and 559 genes were differentially expressed in viruliferous MEAM1 and MED whiteflies, respectively, compared with their non-viruliferous counterparts, of which 146 genes were common between the two cryptic species. For both cryptic species, the number of differentially expressed genes (DEGs) associated with TYLCV and SiGMV acquisition were higher compared with DEGs associated with CuLCrV acquisition. Pathway analysis indicated that the acquisition of begomoviruses induced differential changes in pathways associated with metabolism and organismal systems. Contrasting expression patterns of major genes associated with virus infection and immune systems were observed. These genes were generally overexpressed and underexpressed in B. tabaci MEAM1 and MED adults, respectively. Further, no specific expression pattern was observed among genes associated with fitness (egg production, spermatogenesis, and aging) in viruliferous whiteflies. The weighted gene correlation network analysis of viruliferous B. tabaci MEAM1 and MED adults identified different hub genes potentially implicated in the vector competence and circulative tropism of viruses. Taken together, the results indicate that both vector cryptic species and the acquired virus species could differentially affect gene expression.

Transcriptome analysis of Holotrichia oblita reveals differentially expressed unigenes related to reproduction and development under different photoperiods

Holotrichia oblita (Faldermann) (Coleoptera: Scarabaeidae) is an insect whose feeding and mating behaviors occur at night. A scotophase is necessary for H. oblita reproduction. We used RNA sequencing (RNA-seq) to compare the expression patterns of H. oblita at five photoperiods (0:24, 8:16, 12:12, 16:8, and 24:0 h) (L:D). Compared to the control (24:0) (L:D), 161-684 differentially expressed unigenes (DEUs) were found in female samples, while 698-2322 DEUs were found in male samples. For all DEUs, a total of 92-1143 DEUs were allocated to 116-662 categories of gene ontology (GO), and 81-1116 DEUs were assigned into 77-286 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The iPath diagram showed that the DEUs generated by comparing female and male samples with photoperiods of 0:24 and 24:0, respectively, involved multiple metabolic pathways, such as carbohydrate metabolism, lipid metabolism, nucleotide metabolism, purine metabolism and glutathione metabolism. Most of these DEUs were upregulated. Finally, 13 DEUs related to reproduction and development were selected to confirm the consistency of relative expression between RNA-Seq and quantitative real-time polymerase chain reaction (qRT-PCR). Most of these comparison results agreed well, except for some qRT-PCR results that were not detected in male samples due to their low expression. These results provide useful information for understanding the dark-induced reproduction of H. oblita.

Mutation of Serine protease 1 Induces Male Sterility in Bombyx mori

Serine proteases are important in reproduction, embryonic development, cell differentiation, apoptosis, and immunity. The genes encoding some serine proteases are essential for male fertility in both humans and rodents and are functionally conserved among metazoan. For example, the Serine protease 1 (Ser1) gene determines male reproductive success in the model lepidopteran insect Bombyx mori. In this study, we explored the function of BmSer1 through transgenic CRISPR/Cas9 technology-mediated mutations in silkworm. We found that the mutation of BmSer1 gene resulted in male sterility but had no effect on female fertility. Male mutants produce normal eupyrene sperm bundles, but the sperm bundles do not dissociate into single sperm. Male sterility caused by the BmSer1 gene mutation was inherited stably through female individuals. Therefore, the serine protease encoded by BmSer1 is essential for male reproductive success in lepidopterans and is a potential target gene for biological reproductive regulation.

Homeodomain proteins POU-M2, antennapedia and abdominal-B are involved in regulation of the segment-specific expression of the clip-domain serine protease gene CLIP13 in the silkworm, Bombyx mori

Clip-domain serine proteases (CLIPs) play important roles in insect innate immunity and development. Our previous studies indicated that CLIP13, an epidermis-specific gene, was involved in cuticle remodeling during molting and metamorphosis in the silkworm, Bombyx mori. However, the transcriptional regulatory mechanism and regulatory pathways of CLIP13 remained unclear. In the present study, we investigated CLIP13 expression and the regulation pathway controlled by 20-hydroxyecdysone (20E) in the silkworm. At the transcriptional level, expression of CLIP13 exhibited pronounced spatial and temporal specificity in different regions of the epidermis; homeodomain transcription factors POU-M2, antennapedia (Antp), and abdominal-B (Abd-B) showed similar expression change trends as CLIP13 in the head capsule, thorax, and abdomen, respectively. Furthermore, results of cell transfection assays, electrophoretic mobility shift assays, and chromatin immunoprecipitation demonstrated that POU-M2, Antp, and Abd-B were involved in the transcriptional regulation of CLIP13 by directly binding to their cis-response elements in CLIP13 promoter. RNA interference-mediated silencing of POU-M2, Antp, and Abd-B led to a decrease of CLIP13 expression in the head capsule, the epidermis of the 1st to 3rd thoracic segments and the 7th to 10th abdominal segments, respectively. Consistent with CLIP13, 20E treatment significantly upregulated expression of POU-M2, Antp, and Abd-B in the silkworm epidermis. Taken together, these data suggest that 20E positively regulates transcription of CLIP13 via homeodomain proteins POU-M2, Antp, and Abd-B in different regions of the silkworm epidermis during metamorphosis, thus affecting the molting process. Our findings provide new insight into the functions of homeodomain transcription factors in insect molting.

Genome-wide identification, characterization, and expression profiling of ATP-binding cassette (ABC) transporter genes potentially associated with abamectin detoxification in Cydia pomonella

The codling moth Cydia pomonella L. (Lepidoptera: Tortricidae) is one of the most notorious pests of pome fruits and walnuts worldwide, which has developed resistance to almost all classes of insecticides, including abamectin (ABM). ATP-binding cassette (ABC) transporters are thought to play a vital roles in insecticide detoxification by reducing the toxic concentrations of insecticides in an organism tissues. Despite the tremendous progress in understanding the detoxification mechanisms at the molecular level, the physiological functions of ABC transporters in insects have been poorly investigated. In this study, we found that the ABC inhibitor verapamil synergized significantly the toxicity of ABM, suggesting a potential role of ABC in detoxification. A total of 54 ABC genes were identified in the third-instar larvae of C. pomonella after treatment with sublethal doses (LD10 and LD30) of ABM. The expression profile of these genes in ABM-treated larvae at different time points (24, 48, 72 hr) using transcriptomic analysis (RNA-seq) was also investigated. The results showed that the expression of about 30 ABC genes was significantly co-upregulated after treatment. Several specific genes were up-regulated at 48 hr after treatment of larvae with LD10 ABM. Among these up-regulated genes, we found that the relative expression level of the CPOM19553 was 29.7-fold and 16.0-fold higher when larvae were exposed to ABM at the LD10 and LD30 doses compared to control, respectively. Unlike other ABC genes, only CPOM08323 exhibited significant expression levels in the head and cuticle of the third-instar larvae of C. pomonella exposed to the two sublethal doses of ABM, with no expression was observed in the detoxification tissues such as midgut and Malpighian tubule. This study suggests that these up-regulated genes may be involved in ABM resistance in C. pomonella. Our findings will provide an additional information required for further analysis of ABC transporter genes associated with xenobiotic metabolism in C. pomonella.

Evaluation of Optimal Reference Genes for qRT-PCR Analysis in Hyphantria cunea (Drury)

The relative quantification of gene expression is mainly achieved through reverse transcription-quantitative PCR (qRT-PCR); however, its reliability and precision rely on proper data normalization using one or more optimal reference genes. Hyphantria cunea (Drury) has been an invasive pest of forest trees, ornamental plants, and fruit trees in China for many years. Currently, the molecular physiological role of reference genes in H. cunea is unclear, which hinders functional gene study. Therefore, eight common reference genes, RPS26, RPL13, UBI, AK, RPS15, EIF4A, β-actin, α-tub, were selected to evaluate levels of gene expression stability when subjected to varied experimental conditions, including developmental stage and gender, different tissues, larvae reared on different hosts and different larval density. The geNorm, BestKeeper, ΔCt method, and NormFinder statistical algorithms were used to normalize gene transcription data. Furthermore, the stability/suitability of these candidates was ranked overall by RefFinder. This study provides a comprehensive evaluation of reference genes in H. cunea and could help select reference genes for other Lepidoptera species.

The Draft Genome of Yellow Stem Borer, an Agriculturally Important Pest, Provides Molecular Insights into Its Biology, Development and Specificity Towards Rice for Infestation

Simple Summary

Yellow stem borer (YSB), is the most destructive and widely occurring pest that attacks rice throughout the growing season. Rice (Oryza sativa L.) is a major staple cereal worldwide, providing essential caloric requirements for more than half of the world’s population. Annual losses to rice borers are approximately 5–10%, but losses in individual fields may reach up to 50–60%. The use of traditional pest management strategies in controlling YSB is somewhat challenging due to its unique internal feeding habit. Genome sequence information of economically important crop pests is important for designing or developing pest-resistant rice varieties. In an approach to achieve this, we present our first-ever study on the draft genome sequence of YSB. The information provided from our current study might be useful in developing genome-based approaches for the management of pest species.

Abstract

Yellow stem borer (YSB), Scirpophaga incertulas (Walker) (Lepidoptera: Crambidae), a major monophagous insect pest of rice, causes significant yield losses. The rice–YSB interaction is very dynamic, making it difficult for management. The development of resistant lines has been unsuccessful as there are no effective resistant sources in the germplasm. Genome information is necessary for a better understanding of interaction with rice in terms of its recognition, response, and infestation mechanism. The draft genome of YSB is predicted to have 46,057 genes with an estimated size of 308 Mb, being correlated with the flow cytometry analysis. The existence of complex metabolic mechanisms and genes related to specific behavior was identified, being conditioned by a higher level of regulation. We deciphered the possible visual, olfactory, and gustatory mechanisms responsible for its evolution as a monophagous pest. Comparative genomic analysis revealed that YSB is unique in the way it has evolved. The obvious presence of high-immunity-related genes, well-developed RNAi machinery, and diverse effectors provides a means for developing genomic tools for its management. The identified 21,696 SSR markers can be utilized for diversity analysis of populations across the rice-growing regions. We present the first draft genome of YSB. The information emanated paves a way for biologists to design novel pest management strategies as well as for the industry to design new classes of safer and specific insecticide molecules.

Role of digestive protease enzymes and related genes in host plant adaptation of a polyphagous pest, Spodoptera frugiperda

The evolutionary success of phytophagous insects depends on their ability to efficiently exploit plants as a source of energy for survival. Herbivorous insects largely depend on the efficiency, flexibility, and diversity of their digestive physiology and sophistication of their detoxification system to use chemically diverse host plants as food sources. The fall armyworm, Spodoptera frugiperda (J.E. Smith), is a polyphagous pest of many commercially important crops. To elucidate the ability of this insect pest to adapt to host plant mechanisms, we evaluated the impact of primary (corn) and alternate (rice) host plants after 11 generations on gut digestive enzymatic activity and expression profiles of related genes. Results indicated that the total protease and class-specific trypsin- and chymotrypsin-like protease activity of S. frugiperda significantly differed among host plant treatments. The class-specific protease profiles greatly differed in S. frugiperda midguts upon larval exposure to different treatments with inhibitors compared with treatments without inhibitors. Similarly, the single and cumulative effects of the enzyme-specific inhibitors TLCK, TPCK, and E-64 significantly increased larval mortality and reduced larval growth/mass across different plant treatments. Furthermore, the quantitative reverse transcription polymerase chain reaction results revealed increased transcription of two trypsin (SfTry-3, SfTry-7) and one chymotrypsin gene (Sfchym-9), which indicated that they have roles in host plant adaptation. Knockdown of these genes resulted in significantly reduced mRNA expression levels of the trypsin genes. This was related to the increased mortality observed in treatments compared with the dsRED control. This result indicates possible roles of S. frugiperda gut digestive enzymes and related genes in host plant adaptation.

Mechanisms of Insecticidal Action of Metarhizium anisopliae on Adult Japanese Pine Sawyer Beetles (Monochamus alternatus)

Pine wilt disease, caused by Bursaphelenchus xylophilus (pine wood nematode), leads to severe environmental and economic damage. Here, we report the results of experiments on the biological control of pine wilt disease through termination of the insect vector of the nematode and the mechanism of the insecticidal action of Metarhizium anisopliae JEF-279 against Monochamus alternatus (Japanese pine sawyer). A combined treatment with a fungal conidia suspension and a fungal protease-containing culture filtrate caused 75.8% mortality of the insect vector. Additionally, the presence of destruxins was confirmed in the dead Japanese pine sawyer adults, and half of the 10 protein spots in proteomic analysis were identified as an actin related to muscle contraction. Based on proteomic and microscopic analyses, the infection cycle of the Japanese pine sawyer by M. anisopliae JEF-279 was inferred to proceed in the following sequence: (1) host adhesion and germination, (2) epicuticle degradation, (3) growth as blastospore, (4) killing by various fungal toxins (insecticidal metabolites), (5) immune response as defense mechanism, and (6) hyphal extrusion and conidiation. Consequently, the combined fungal conidia suspension and protease-containing culture filtrate treatment may be applied as an insecticidal agent, and flaccid paralysis is likely a major mechanism underlying the insecticidal action of M. anisopliae JEF-279 on host insects.

Immune responses to Bacillus thuringiensis in the midgut of the diamondback moth, Plutella xylostella

The diamondback moth, Plutella xylostella, is the first insect to develop resistance to Bacillus thuringiensis (Bt) in the field. To date, little is known about the molecular mechanism of the interaction between Bt and midgut immunity in P. xylostella. Here, we report immune responses in the P. xylostella midgut to Bt strain Bt8010 using a combined approach of transcriptomics and quantitative proteomics. Many genes in the Toll, IMD, JNK and JAK-STAT pathways and antimicrobial peptide genes were activated at 18 h post-infection. In the prophenoloxidase (PPO) cascade, four serpin genes were activated, and the PPO1 gene was suppressed by Bt8010. Inhibition of the two PPO proteins was observed at 18 h post-infection. Feeding Bt8010-infected larvae recombinant PPOs enhanced their survival. These results revealed that the Toll, IMD, JNK and JAK-STAT pathways were triggered and participated in the immune defence of the midgut against Bt8010, while the PPO cascade was inhibited and played an important role in this process.

Identification, characterization, and expression analysis of clip-domain serine protease genes in the silkworm, Bombyx mori

Clip-domain serine proteases (CLIPs), characterized by regulatory module clip domains, constitute an important serine protease family identified in insects and other arthropods. They participate in host immune response and embryonic development in a cascade-activated manner. Here, we present a genome-wide identification and expression analysis of CLIP genes in the silkworm, Bombyx mori. A total of 26 CLIP genes were identified in the silkworm genome. Bioinformatics analysis indicated that these CLIPs clustered into four subfamilies (CLIPA-D), and exhibit a close evolutionary relationship with CLIPs of Manduca sexta. Tissue expression profiling revealed that silkworm CLIP genes are mainly expressed in the integument, head, fat body, and hemocytes. Temporal expression profiles showed that 15 CLIP genes were predominantly expressed during the fifth-instar larval stage, early and later period of the pupal stage, and adult stage, whereas 10 CLIP genes were mainly expressed in the wandering stage and middle to later period of the pupal stage in the integument. Pathogens and 20-hydroxyecdysone (20E) induction analysis indicated that 14 CLIP genes were positively regulated by 20E, 9 were negatively regulated by 20E but positively regulated by pathogens, and 5 were positively regulated by both factors in the integument. Together, these results suggested that silkworm CLIP genes may play multiple functions in integument development, including melanization of new cuticle, molting and immune defense. Our data provide a comprehensive understanding of CLIP genes in the silkworm integument and lays a foundation for further functional studies of CLIP genes in the silkworm.

Clip-Domain Serine Protease Gene (LsCLIP3) Is Essential for Larval-Pupal Molting and Immunity in Lasioderma serricorne

Clip-domain serine proteases (CLIPs) play crucial roles in insect development and innate immunity. In this study, we identified a CLIP gene (designated LsCLIP3) from the cigarette beetle Lasioderma serricorne. LsCLIP3 contains a 1,773-bp open reading frame (ORF) encoding a 390-amino-acid protein and shows a conserved clip domain and a trypsin-like serine protease domain. Phylogenetic analysis indicated that LsCLIP3 was orthologous to the CLIP-B subfamily. LsCLIP3 was prominently expressed in larva, pupa, and early adult stages. In larval tissues, it was highly expressed in the integument and fat body. The expression of LsCLIP3 was induced by 20-hydroxyecdysone. A similar induction was also found by peptidoglycans from Escherichia coli and Staphylococcus aureus. RNA interference (RNAi)-mediated silencing of LsCLIP3 disrupted larval–pupal molting and specifically reduced the expression of genes in 20-hydroxyecdysone synthesis and signaling pathway. The chitin amounts of LsCLIP3 RNAi larvae were greatly decreased, and expressions of six chitin metabolic-related genes were significantly reduced. Knockdown of LsCLIP3 increased larval sensitivity to Gram-negative and Gram-positive bacteria. There was significantly decreased expression of four antimicrobial peptide (AMP) genes. The results suggest that LsCLIP3 is an important component of the larva to pupa molt and for the immunity of L. serricorne.

Genome-wide profiling of the alternative splicing provides insights into development in Plutella xylostella

Background

The diamondback moth (DBM), Plutella xylostella (L.), is a major pest of cruciferous crops worldwide. While the species has become a model for genomics, post-transcriptional mechanisms associated with development and sex determination have not been comprehensively studied and the lack of complete structure of mRNA transcripts limits further research.

Results

Here, we combined the methods of single-molecule long-read sequencing technology (IsoSeq) and RNA-seq to re-annotate the published DBM genome and present the genome-wide identification of alternative splicing (AS) associated with development and sex determination of DBM. In total, we identified ~ 13,900 genes (~ 77%) annotated in the DBM genome (version-2), resulting in the correction of 1586 wrongly annotated genes and identification of 78,000 previously unannotated transcripts. We also identified 1804 genes showing alternative splicing (AS) in each of the developmental stages and sexes, suggesting that AS events are ubiquitous in DBM. Comparative analyses showed that these AS events were rarely shared among developmental stages, indicating that they may play key specific roles in regulation of insect development. Further, we found 156 genes showing different AS events and expression patterns between males and females, linking them to potential functions in sex determination.

Conclusion

Overall, the P. xylostella transcriptome provides the significant information about regulatory alternative splicing events, which are shown to be involved in development and sex determination. Our work presents a solid foundation to better understand the mechanism of post-transcriptional regulation, and offers wider insights into insect development and sex determination.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5838-3) contains supplementary material, which is available to authorized users.

Building a platform for predicting functions of serine protease-related proteins in Drosophila melanogaster and other insects

Serine proteases (SPs) and serine protease homologs (SPHs) play essential roles in insect physiological processes including digestion, defense and development. Studies of insect genomes, transcriptomes and proteomes have generated a vast amount of information on these proteins, dwarfing the biological data acquired from a few model species. The large number and high diversity of homologous sequences makes it a challenge to use the limited functional information for making predictions across a broad taxonomic group of insects. In this work, we have extensively updated the framework of knowledge on the SP-related proteins in Drosophila melanogaster by identifying 52 new SPs/SPHs, classifying the 257 proteins into four groups (CLIP, gut, single- and multi-domain SPs/SPHs), and detecting inherent connections among phylogenetic relationships, genomic locations and expression profiles for 99 of the genes. Information on the existence of specific proteins in eggs, larvae, pupae and adults is presented to facilitate future research. More importantly, we have developed an approach to reveal close homologous or orthologous relationships among SPs/SPHs from D. melanogaster, Anopheles gambiae, Apis mellifera, Manduca sexta, and Tribolium castaneum thus inspiring functional studies in these and other holometabolous insects. This approach is useful for tackling similar problems on large and diverse protein families in other groups of organisms.

Insights into the Temporal Gene Expression Pattern in Lymantria dispar Larvae During the Baculovirus Induced Hyperactive Stage

Baculoviruses are effective biological control agents for many insect pests. They not only efficiently challenge the host immune system but also make them hyperactive for better virus dispersal. Some investigations have focused on the viral mechanisms for induction of such altered response from the host. However, there are no current studies monitoring changes in gene expression during this altered phenotype in infected larvae. The L. dispar multiple nucleopolyhedrovirus (LdMNPV) induces hyperactivity in third instar L. dispar larvae at 3-days post infection (dpi), to continued till 6 dpi. The transcriptome profiles of the infected and uninfected larvae at these time points were analyzed to provide new clues on the response of the larvae towards infection during hyperactivity. Gene ontology enrichment analysis revealed, most of the differentially expressed genes (DEGs) were involved in proteolysis, extracellular region, and serine-type endopeptidase activity. Similarly, Kyoto Encyclopedia of Genes and Genome enrichment analysis showed maximum enrichment of 487 genes of the signal transduction category and neuroactive ligand-receptor interaction sub-category with 85 annotated genes. In addition, enrichment map visualization of gene set enrichment analysis showed the coordinated response of neuroactive ligand-receptor interaction genes with other functional gene sets, as an important signal transduction mechanism during the hyperactive stage. Interestingly all the DEGs in neuroactive ligand-receptor interactions were serine proteases, their differential expression during the hyperactive stage correlated with their conceivable involvement in disease progression and the resulting altered phenotype during this period. The outcome provides a basic understanding of L. dispar larval responses to LdMNPV infection during the hyperactive stage and helps to determine the important host factors involved in this process.

A serine protease homologue Bombyx mori scarface induces a short and fat body shape in silkworm

Body shape is one of the most prominent and basic characteristics of any organism. In insects, abundant variations in body shape can be observed both within and amongst species. However, the molecular mechanism underlying body shape fine-tuning is very complex and has been largely unknown until now. In the silkworm Bombyx mori, the tubby (tub) mutant has an abnormal short fat body shape and the abdomen of tub larvae expands to form a fusiform body shape. Morphological investigation revealed that the body length was shorter and the body width was wider than that of the Dazao strain. Thus, this mutant is a good model for studying the molecular mechanisms of body shape fine-tuning. Using positional cloning, we identified a gene encoding the serine protease homologue, B. mori scarface (Bmscarface), which is associated with the tub phenotype. Sequence analysis revealed a specific 312-bp deletion from an exon of Bmscarface in the tub strain. In addition, recombination was not observed between the tub and Bmscarface loci. Moreover, RNA interference of Bmscarface resulted in the tub-like phenotype. These results indicate that Bmscarface is responsible for the tub mutant phenotype. This is the first study to report that mutation of a serine protease homologue can induce an abnormal body shape in insects.

The genomic and transcriptomic analyses of serine proteases and their homologs in an endoparasitoid, Pteromalus puparum

In insects, serine proteases (SPs) and serine protease homologs (SPHs) constitute a large family of proteins involved in multiple physiological processes such as digestion, development, and immunity. Here we identified 145 SPs and 38 SPHs in the genome of an endoparasitoid, Pteromalus puparum. Gene duplication and tandem repeats were observed in this large SPs/SPHs family. We then analyzed the expression profiles of SP/SPH genes in response to different microbial infections (Gram-positive bacterium Micrococcus luteus, Gram-negative bacterium Escherichia coli, and entomopathogenic fungus Beauveria bassiana), as well as in different developmental stages and tissues. Some SPs/SPHs also displayed distinct expression patterns in venom gland, suggesting their specific physiological functions as venom proteins. Our finding lays groundwork for further research of SPs and SPHs expressed in the venom glands.

Hemolymph proteins of Anopheles gambiae larvae infected by Escherichia coli

Anopheles gambiae is a major vector of human malaria and its immune system in part determines the fate of ingested parasites. Proteins, hemocytes and fat body in hemolymph are critical components of this system, mediating both humoral and cellular defenses. Here we assessed differences in the hemolymph proteomes of water- and E. coli-pricked mosquito larvae by a gel-LC-MS approach. Among the 1756 proteins identified, 603 contained a signal peptide but accounted for two-third of the total protein amount on the quantitative basis. The sequence homology search indicated that 233 of the 1756 may be related to defense. In general, we did not detect substantial differences between the control and induced plasma samples in terms of protein numbers or levels. Protein distributions in the gel slices suggested post-translational modifications (e.g. proteolysis) and formation of serpin-protease complexes and high Mr immune complexes. Based on the twenty-five most abundant proteins, we further suggest that major functions of the larval hemolymph are storage, transport, and immunity. In summary, this study provided first data on constitution, levels, and possible functions of hemolymph proteins in the mosquito larvae, reflecting complex changes occurring in the fight against E. coli infection.

Optimization of recombinant bacteria expressing dsRNA to enhance insecticidal activity against a lepidopteran insect, Spodoptera exigua

Double-stranded RNA (dsRNA) has been applied to control insect pests due to its induction of RNA interference (RNAi) of a specific target gene expression. However, developing dsRNA-based insecticidal agent has been a great challenge especially against lepidopteran insect pests due to variations in RNAi efficiency. The objective of this study was to screen genes of chymotrypsins (SeCHYs) essential for the survival of the beet armyworm, Spodoptera exigua, to construct insecticidal dsRNA. In addition, an optimal oral delivery method was developed using recombinant bacteria. At least 7 SeCHY genes were predicted from S. exigua transcriptomes. Subsequent analyses indicated that SeCHY2 was widely expressed in different developmental stages and larval tissues by RT-PCR and its expression knockdown by RNAi caused high mortality along with immunosuppression. However, a large amount of dsRNA was required to efficiently kill late instars of S. exigua because of high RNase activity in their midgut lumen. To minimize dsRNA degradation, bacterial expression and formulation of dsRNA were performed in HT115 Escherichia coli using L4440 expression vector. dsRNA (300 bp) specific to SeCHY2 overexpressed in E. coli was toxic to S. exigua larvae after oral administration. To enhance dsRNA release from E. coli, bacterial cells were sonicated before oral administration. RNAi efficiency of sonicated bacteria was significantly increased, causing higher larval mortality at oral administration. Moreover, targeting young larvae possessing weak RNase activity in the midgut lumen significantly enhanced RNAi efficiency and subsequent insecticidal activity against S. exigua.

Genome-wide Identification and Expression Analysis of Amino Acid Transporters in the Whitefly, Bemisia tabaci (Gennadius)

The whitefly (Bemisia tabaci) is a cosmopolitan and devastating pest of agricultural crops and ornamentals. B. tabaci causes extensive damage by feeding on phloem and by transmitting plant viruses. Like many other organisms, insects depend on amino acid transporters (AATs) to transport amino acids into and out of its cells. We present a genome-wide and transcriptome-wide investigation of the following two families of AATs in B. tabaci biotype B: amino acid/auxin permease (AAAP) and amino acid/polyamine/organocation (APC). A total of 14 putative APCs and 25 putative AAAPs were identified, and a 10-paralog B. tabaci-specific expansion of AAAPs was found by maximum likelihood phylogeny. Detailed gene structure information revealed that 9 members of the B. tabaci-specific AAAP family expansion closely situated on a same scaffold. Expression profiling of the B. tabaci B APC and AAAP genes as affected by stage and plant host showed diverse expression patterns. The analysis of evolutionary rates indicated that purifying selection can explain the B. tabaci-specific AAAP expansion. RNA interference (RNAi)-mediated suppression of two AAAP genes (BtAAAP15 and BtAAAP21) significantly increased the mortality of B. tabaci B adults. The results provide a foundation for future functional analysis of APC and AAAP genes in B. tabaci.

Characterization and expression profiling of serine protease inhibitors in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

BACKGROUND

Serine protease inhibitors (SPIs) have been found in all living organisms and play significant roles in digestion, development and innate immunity. In this study, we present a genome-wide identification and expression profiling of SPI genes in the diamondback moth, Plutella xylostella (L.), a major pest of cruciferous crops with global distribution and broad resistance to different types of insecticides.

RESULTS

A total of 61 potential SPI genes were identified in the P. xylostella genome, and these SPIs were classified into serpins, canonical inhibitors, and alpha-2-macroglobulins based on their modes of action. Sequence alignments showed that amino acid residues in the hinge region of known inhibitory serpins from other insect species were conserved in most P. xylostella serpins, suggesting that these P. xylostella serpins may be functionally active. Phylogenetic analysis confirmed that P. xylostella inhibitory serpins were clustered with known inhibitory serpins from six other insect species. More interestingly, nine serpins were highly similar to the orthologues in Manduca sexta which have been demonstrated to participate in regulating the prophenoloxidase activation cascade, an important innate immune response in insects. Of the 61 P.xylostella SPI genes, 33 were canonical SPIs containing seven types of inhibitor domains, including Kunitz, Kazal, TIL, amfpi, Antistasin, WAP and Pacifastin. Moreover, some SPIs contained additional non-inhibitor domains, including spondin_N, reeler, and other modules, which may be involved in protein-protein interactions. Gene expression profiling showed gene-differential, stage- and sex-specific expression patterns of SPIs, suggesting that SPIs may be involved in multiple physiological processes in P. xylostella.

CONCLUSIONS

This is the most comprehensive investigation so far on SPI genes in P. xylostella. The characterized features and expression patterns of P. xylostella SPIs indicate that the SPI family genes may be involved in innate immunity of this species. Our findings provide valuable information for uncovering further biological roles of SPI genes in P. xylostella.