Next Article in Journal
Prospecting Local Treatments Used in Conjunction with Antivenom Administration Following Envenomation Caused by Animals: A Systematic Review
Next Article in Special Issue
Unveiling the Protein Components of the Secretory-Venom Gland and Venom of the Scorpion Centruroides possanii (Buthidae) through Omic Technologies
Previous Article in Journal
A Novel Multicellular Placental Barrier Model to Investigate the Effect of Maternal Aflatoxin B1 Exposure on Fetal-Side Neural Stem Cells
Previous Article in Special Issue
Application of an Extracellular Matrix-Mimicking Fluorescent Polymer for the Detection of Proteolytic Venom Toxins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Characterization and Functional Analysis of the Dipeptidyl Peptidase IV from Venom of the Ectoparasitoid Scleroderma guani

1
Key Laboratory of Forest Disaster Warning and Control of Yunnan Province, College of Biodiversity Conservation, Southwest Forestry University, Kunming 650224, China
2
Key Laboratory for Forest Resources Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, College of Biodiversity Conservation, Southwest Forestry University, Kunming 650224, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Toxins 2023, 15(5), 311; https://doi.org/10.3390/toxins15050311
Submission received: 13 March 2023 / Revised: 13 April 2023 / Accepted: 25 April 2023 / Published: 27 April 2023
(This article belongs to the Special Issue Animal Venoms: Proteomics, Biochemical Activities and Application)

Abstract

:
Dipeptidyl peptidase IV (DPPIV) is a proline-specific serine peptidase that remains poorly investigated in terms of venom composition. Here, we describe the molecular characteristics and possible functions of DPPIV as a major venom component of the ant-like bethylid ectoparasitoid, Scleroderma guani, named SgVnDPPIV. The SgVnDPPIV gene was cloned, which encodes a protein with the conserved catalytic triads and substrate binding sites of mammalian DPPIV. This venom gene is highly expressed in the venom apparatus. Recombinant SgVnDPPIV, produced in Sf9 cells using the baculovirus expression system, has high enzymatic activity, which can be efficiently inhibited by vildagliptin and sitagliptin. Functional analysis revealed that SgVnDPPIV affects genes related to detoxification, lipid synthesis and metabolism, response to stimuli, and ion exchange in pupae of Tenebrio molitor, an envenomated host of S. guani. The present work contributes towards understanding the role of venom DPPIV involved in the interaction between parasitoid wasp and its host.
Key Contribution: Dipeptidyl peptidase IV (DDPIV) is common in the venom of parasitoid wasps, while its functions remain to be unraveled. In this work, this venom gene was identified in the ectoparasitoid, Scleroderma guani, followed by functional characterization. The results provide insights into the biological role of DPPIV in the venom of parasitoid wasps.

1. Introduction

Dipeptidyl peptidase IV (DPPIV, EC 3.4.14.5) belongs to the subfamily S9B of serine peptidases from peptidase clan SC in the MEROPS database (https://www.ebi.ac.uk/merops/ (accessed on 6 January 2023) [1]. As an aminopeptidase, it preferentially cleaves dipeptides from peptides with a penultimate proline or alanine in the P1 position, such as chemokines, neuropeptides, and peptide hormones [2]. This enzyme is found in a wide range of different living organisms, including mammals, insects, plants, and microorganisms, which has been mostly and best studied in humans [3]. Human DPPIV is widely distributed in cells and tissues to act in various physiological and pathological processes, ranging from glucose homeostasis, immunoregulation, and inflammation to tumorigenesis, through an anchored transmembrane molecule or a soluble circulating enzyme [4,5]. The prominent role of this enzyme is to rapidly inactivate the incretin peptide hormones, glucagon-like peptide-1 and glucose-dependent insulinotropic polypeptide, that are involved in regulating blood glucose levels [6]. Thus, it provides as an ideal pharmacological target to develop drugs for the treatment of type 2 diabetes by potentiating glucose-dependent insulin secretion [7,8].
Compared with mammalian DPPIV, little is known about insect DPPIV, with information derived from limited studies conducted on a few species, including blue blowfly Calliphora vicina [9], cockroach Leucophaea maderae [10], fruit fly Drosophila melanogaster [11], yellow mealworm Tenebrio molitor [3], honey bee Apis mellifera [12], and several wasps [12,13,14,15,16,17]. Insect DPPIV was first partially purified from C. vicina, which has been shown to be able to degrade the ecdysiostatic peptide, i.e., trypsin-modulating oostatic factor, of the grey flesh fly Neobellieria bullata [9]. Partly purified DPPIV activity was characterized from the brain and midgut of L. maderae, and is thought to participate in the inactivation of tachykinin-related peptides [10]. A gene (omega) encoding DPPIV was cloned from D. melanogaster [11]. However, the functions of DPPIV from these insects are not explored. It was found that a soluble DPPIV from the anterior midgut of T. molitor larvae can efficiently hydrolyze gliadins, the main dietary protein of the insect, showing that it has a digestive function [3]. Interestingly, DPPIV is the major component in the venom of hymenopterans. This enzyme from the hornets Vespa basalis and V. magnifica, has been considered as the putative enzyme used to process toxins, including mastoparan B and vespid chemotactic peptide, with antimicrobial activities [13,14,17,18]. It is a major venom allergen of A. mellifera, Polistes dominula, and V. velutina nigrithorax [12,15,16].
Parasitoids likely comprise the largest group of venomous organisms, and are one of the most diverse groups of Hymenoptera [19]. Their venoms comprise a complex cocktail of proteinaceous and non-proteinaceous components with diverse biological functions that facilitate the successful development of its offspring either inside or outside the host [20,21]. Enzymatic proteins, such as serine proteinase, esterase, and metalloprotease, constitute the most abundant venom components of parasitoids [22,23,24]. Of these, DPPIV has been commonly found [25,26]. Unfortunately, the molecular characteristics and functions of DPPIV as a constituent of parasitoid venoms are still unknown.
The ant-like bethylid wasp, Scleroderma guani, is a valuable biological control agent that can parasitize the larvae or pupae of many woodboring beetles [27], and has been successfully used to control the Japanese pine sawyer, a notorious forest pest in Asia known as an efficient vector of the destructive pine wilt disease [28]. A DPPIV (SgVnDPPIV) was identified as a major venom component of this ectoparasitoid [29]. We here describe the molecular characterization and functional analysis of SgVnDPPIV.

2. Results

2.1. Molecular Characteristics of SgVnDPPIV

The cloned SgVnDPPIV gene sequence is 2722 bp in length, containing an open reading frame (ORF) of 2301 bp (Figure S1). The predicted protein contains 767 amino acids, with a calculated molecular weight of 87.30 kDa and isoelectric point of 6.34. There is no signal peptide predicted in the deduced amino acid sequence. SgVnDPPIV shares 47%, 38%, 34%, and 31% identity with DPPIVs of V. basalis, T. molitor, Microplitis demolitor, and Homo sapiens, respectively (Figure 1). The catalytic triads and substrate binding sites of mammalian DPPIV are conserved in SgVnDPPIV. Genes encoding DPPIVs in the genomes of different hymenopteran species vary from 4 to 11 (Tables S1 and S2). In addition to the SgVnDPPIV gene, there are three DPPIV genes identified in the genome of S. guani. Phylogenetic analysis revealed that hymenopteran DPPIVs can be divided into five groups (Figure 2). SgVnDPPIV and the other three DPPIVs of S. guani are grouped into different clades.

2.2. Gene Expression Pattern of SgVnDPPIV

According to the qPCR experimental results (Figure 3A), the SgVnDPPIV gene is most highly expressed in the venom apparatus. The transcriptional level of SgVnDPPIV in the venom apparatus is significantly higher than in abdomens deprived of the venom apparatus and other tissues. Its expressions are dynamic across life stage. Scleroderma guani expresses this venom gene at a very low or even undetectable level in eggs, larvae, and pupae in white and yellow cocoons (Figure 3B). However, its expression then significantly increases in pupae in black cocoons, and reaches to the highest expression level in newly emerged adults. The expression of the SgVnDPPIV gene remains at a high level in other detected adult stages, with an overall decreasing tendency.

2.3. Enzymatic Activity of SgVnDPPIV

Western blot analysis showed that recombinant SgVnDPPIV was successfully produced in baculovirus-infected Spodoptera frugiperda Sf9 cells (Figure S2). The production was found in the cell medium and the soluble fraction of cell lysates. High-purity SgVnDPPIV was purified from concentrated cell culture supernatant (Figure 4A). Analyses of the enzymatic activity of SgVnDPPIV indicated that it has higher activity than that of crude venom from S. guani at their equivalent amounts (Figure 4B). Its activity towards the substrate of Ala-Pro-p-nitroanilide is comparable to the gut extract of T. molitor larvae. Vildagliptin and sitagliptin can efficiently inhibit the activity of SgVnDPPIV (Figure 4C). There are no significant differences between their inhibitory ability towards SgVnDPPIV. However, diprotin A and B have no evident inhibitory effect on SgVnDPPIV activity.

2.4. Functional Analysis of SgVnDPPIV

Transcriptomic sequencing led to the generation of 531555784 raw reads of SgVnDPPIV-injected pupae of T. molitor and their controls (Table S3). After adapter and low-quality nucleotide removal, the raw data were reduced to 79.31 G clean data. Clean reads from each sample were mapped to the genes encoding in the genome of T. molitor to calculate their transcript abundance, expressed in TPM (transcripts per million). Differential gene expression analysis revealed that a total of 455 genes are significantly differentially expressed in comparison with the control (Table S4), which only represents 0.99% of the total 45,780 genes encoding in the genome of T. molitor. Overall, injection with SgVnDPPIV, after 6 h and 24 h, influences the expression of 29 (20 upregulated and 9 downregulated) and 433 (290 upregulated and 143 downregulated) genes of T. molitor pupae, respectively (Figure 5A,B). In Figure 5C, the Venn diagram exhibits that the expression of only seven genes change at both time points upon injection with SgVnDPPIV.
By utilizing GO enrichment analysis at level three, the significantly differentially expressed genes 6 h after injection with SgVnDPPIV were enriched in the categories of oxidoreductase activity (7), organelle subcompartments (7), small molecule metabolic processes (11), and response to chemical stimuli (13), respectively (Figure S3A). Those genes significantly differentially expressed 24 h after injection with SgVnDPPIV are enriched in the GO terms of M band (6), oxidoreductase activity (31), organelle (142), and cellular response to stimuli (109) (Figure S3B). Detailed examination of the functions of the enriched genes revealed that they are mostly related to detoxification, lipid synthesis and metabolism, response to stimuli of different compounds, such as xenobiotics, alkaloids, and alcohols, and ion exchange (Figure 6 and Table S5).

3. Discussion

The present study reports the cloning and sequencing of venom DPPIV cDNA from the ectoparasitoid, S. guani (SgVnDPPIV). Similar to mammalian DPPIV [3], SgVnDPPIV lacks a signal peptide and has a transmembrane domain, indicating that it is a transmembrane enzyme. In contrast, DPPIVs of T. molitor and some other insects are secreted enzymes [3]. Although most venom components are recognized as secreted proteins with a signal peptide, a few of them have been revealed to lack a predicted signal peptide [30]. Unlike most venom proteins, despite the fact that the characteristic signal peptide is present in DPPIV from the venom of the pit viper, Gloydius blomhoffi brevicaudus, it is not cleaved off during biosynthesis [31]. According to the phylogenetic result, SgVnDPPIV is grouped into a separate clade, and the other three non-venomous DPPIVs from S. guani are assigned into other separate groups, suggesting that their functions are different. A high gene expression level of SgVnDPPIV was observed in the venom apparatus of S. guani, and strong DPPIV activity was detected in its venom. According to the qPCR and enzymatic assay results, SgVnDPPIV is an abundant venom component of S. guani, which has been evidenced using the proteomic approach [29]. Abundant expression of the SgVnDPPIV gene occurs in pupae in black cocoons, and its expression levels are highly maintained at early adult stages, indicating that the transcription of this venom gene is correlated with the development of the venom apparatus, which is responsible for the production of venom components, and might be due to the biological rhythm of female wasps [32]. The three catalytic triads (Ser, Asp, and His) for serine peptidases are present in the deduced amino acid sequence of SgVnDPPIV [33]. Core substrate binding sites that are important for the activity of mammalian DPPIVs are highly conserved in the amino acid sequence of this venom protein [34,35]. These findings demonstrate that SgVnDPPIV is an active venomous enzyme. In agreement with this, the active recombinant SgVnDPPIV was produced in Sf9 cells for enzymatic and functional characterization.
The activity of SgVnDPPIV is sensitive to the inhibitors of vildagliptin and sitagliptin, but diprotin A and B have no obvious inhibitory effect on it. This is similar to DPPIVs of humans and T. molitor, which are not efficiently inhibited by diprotin A or B [3]. Vildagliptin and sitagliptin are the most effective inhibitors of human DPPIV [35]. Vildagliptin strongly inhibits T. molitor DPPIV, but its activity is not apparently influenced by sitagliptin [3]. These differences might be a result of the differences in the structures of the S1 and S2 subsites of DPPIVs from various species, which bind to different types of inhibitors with different modes of action.
DPPIV is an important enzyme used to degrade proline-containing peptides [2,9]. Recently, it has been found to participate in the hydrolysis of dietary gliadin, which is rich in proline residues, performing a digestive function in the larval gut of T. molitor [3]. As a venom component of the honey bee, hornet, and paper wasp, DPPIV was characterized as an allergen, and may be involved in processing peptidic toxins [12,13,14,15,16,17]. Venom of these wasps commonly causes allergic reactions among humans, and is used to defend their colonies from vertebrate predators; this is different from the venom of parasitoids, which is responsible for paralyzing hosts and manipulating their parasitized hosts for the survival and development of their progeny [21,36,37,38]. Additionally, peptides are often abundant in the venom of Aculeate bees and wasps, while they are rarely identified in the venom of parasitoids [39,40]. According to these differences in the venoms of Aculeate groups and parasitoids, DPPIV in the venom of parasitoids may act in biological processes other than allergic reactions or the maturation of proline-rich peptidic toxins, and might be shaped in response to both ecological and behavioral influences during the evolution of parasitoids.
In order to reveal the putative function of SgVnDPPIV, we investigated the influence of this venom protein on the expression of genes encoding in the genome of T. molitor. Only a very small portion of genes of T. molitor pupae changed their transcriptions after this beetle was envenomated with SgVnDPPIV. This is similar to venom Y, a major component of Nasonia vitripennis venom, in that the presence of this venom can only lead to changes of 1.8% of the total genes of its host, flesh fly Sarcophaga bullata [41]. The majority of them are detoxification, immunity, neural, and fat-body-specific genes [41]. Among the genes affected by SgVnDPPIV, many of them are involved in oxidation reduction, lipid synthesis and metabolism, response to stimuli, and ion exchange, indicating the corresponding function of SgVnDPPIV. These results pave the way towards a better understanding of the role of DPPIV in parasitoid venom; however, the mechanisms by which this is achieved remain to be unraveled.

4. Conclusions

Herein, venom DPPIV of S. guani (SgVnDPPIV) was molecularly characterized and functionally analyzed. The results indicate that this venom protein has conserved catalytic triads and substrate binding sites of mammalian DPPIV, exhibiting high enzymatic activity that can be inhibited by vildagliptin and sitagliptin. It can function by affecting the expression of hosts genes related to detoxification, lipid synthesis and metabolism, stimulus response, and ion exchange. These findings will facilitate the understanding of the biological role of DPPIV in the venom of parasitoids.

5. Materials and Methods

5.1. Insects

The parasitoid, S. guani, and its host, T. molitor, used in this study were derived from colonies maintained in our laboratory, as described by Zhu et al. [42]. In brief, newly emerged adult wasps were fed daily with fresh 20% honey solution in water at room temperature and 70% relative humidity under a 14 h light/10 h dark diurnal cycle. Pupae of T. molitor were supplied as hosts for parasitization by S. guani. These beetles were reared using wheat bran containing dry instant yeast as an artificial diet, and kept in the dark at about 55% relative humidity.

5.2. Gene Cloning and Sequence Analysis

The venom apparatus was dissected from about 150 female adults of S. guani. Total RNA was isolated from the venom apparatus using Trizol reagent (TaKaRa, Dalian, China). Its purity and concentration were measured using a NanoDrop ND 1000 spectrometer (PeqLab, Erlangen, Germany). Its quality was verified using 1.0% agarose gel electrophoresis. Aliquots of total RNA (5 µg) were used to synthesize the cDNA using a RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) with the 3′-RACE CDS primer (5′-AAGCAGTGGTATCAACGCAGAGTAC(T)30-3′). The gene sequence of SgVnDPPIV was derived from the venom apparatus transcriptome of S. guani [29]. Based on this, the gene-specific forward primer (5′-GGAACAGTATCTGGACAGACTAAGGAAT-3′) was designed with Primer 5.0 (PREMIER Biosoft International, Palo Alto, CA, USA). This primer and the 3′-RACE CDS primer were used to clone the SgVnDPPIV gene. PCR (polymerase chain reaction) conditions were as follows: 3 min of pre-denaturation at 94 °C, 40 cycles of 30 s at 94 °C, 1 min at 60 °C, 2 min 30 s at 72 °C, and the final extension at 72 °C of 10 min. PCR product was detected using 1.0% agarose gel electrophoresis, and then sent to Biosune Biotechnology Limited Company (Shanghai, China) for sequencing. Signal peptide was predicted with SignalP 6.0 (https://services.healthtech.dtu.dk/service.php?SignalP-6.0 accessed on 20 October 2022) [43]. Genes encoding DPPIVs of various hymenopteran species were identified using the BlastP in TBtools v1.0987663 [44] from their genomes (Tables S1 and S2). Multiple sequence alignment of the amino acid sequences was performed with ClustalX 2.1 [45]. A phylogenetic tree was constructed using IQ-TREE 2 based on alignment with 1000 ultra-fast bootstrap replicates [46].

5.3. Quantitative Real-Time PCR

Selected tissues, including the venom apparatus, abdomens deprived of the venom apparatus, thoraxes, and heads from female adults and parasitoids from different developmental stages, including eggs, early larvae, late larvae, mature larvae, spinning mature larvae, pupae in white cocoon, pupae in yellow cocoon, pupae in black cocoon, and adults after 1–25 days emergence, were collected. Total RNA was isolated, and its purity and quality were verified as described above. Total RNA from these samples was used to synthesize cDNA with the PrimeScript RT Reagent Kit with gDNA Eraser (TaKaRa, Dalian, China). Gene-specific primers (5′-CGGATACAAGTGCAACGTCAG-3′ and 5′-TTGGTG GTCAAGTAGGAGTGG-3′) were designed using Primer 5.0 (PREMIER Biosoft International, Palo Alto, CA, USA) according to the corresponding gene sequence. The 18S ribosomal RNA gene was used as an endogenous control. Quantitative real-time PCR (qPCR) was performed using SYBR Premix Ex Taq II (TaKaRa, Dalian, China) to quantitatively evaluate the level of gene transcription in triplicate samples. Amplification reactions were 95 °C for 2 min followed by 45 cycles of 95 °C for 15 s and 58 or 60 °C for 30 s, finishing with elongation at 72 °C for 30 s. The 2−ΔΔCT method was used to analyze the qPCR data [47].

5.4. Recombinant Production of Proteins

According to the sequence of the SgVnDPPIV gene described above, its cDNA containing Nco I and Hind III ends was synthesized via a PCR-based accurate synthesis protocol [48]. After validation by Sanger sequencing, it was cut with Nco I and Hind III and then cloned into the pFast-bac1 vector provided in the Bac-to-Bac® Vector Kit (Thermo Fisher Scientific, Waltham, MA, USA), and the final vector was verified by sequencing. This was used to generate bacmid following the manufacture’s instruction. Generated bacmid was used to transfect the Sf9 cells cultivated in Sf-900 II SFM (Thermo Fisher Scientific, Waltham, MA, USA) media at a density of about 9 × 105 cells/mL using Cellfectin II Reagent (Thermo Fisher Scientific, Waltham, MA, USA). Infected cells were maintained in Sf-900 II SFM at 27 °C and 100 rpm to amplify the baculovirus. At 72 h post infection, cells and medium were harvested to determine the expression of the recombinant protein using 10% SDS-PAGE. Recombinant SgVnDPPIV bacmid-derived baculoviruses were amplified in an additional round of infection for high levels of protein expression. The supernatant of cell medium was collected to purify the recombinant proteins by affinity separation with Ni-IDA -Sepharose CL-6B (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. Bound protein fractions were eluted with 300 mM imidazole. Fractions containing the SgVnDPPIV were pooled and concentrated for subsequent 10% SDS-PAGE analysis and stained with Coomassie Brilliant Blue (TaKaRa, Dalian, China). Western blotting was performed using anti-His Tag mouse monoclonal antibodies (TaKaRa, Dalian, China). The concentration of purified proteins was measured by the Quick Start Bradford Protein Assay (Bio-Rad, Hercules, CA, USA). They were stored at −80 °C prior to use.

5.5. Enzyme Assay

The enzymatic activity of SgVnDPPIV was assayed with substrate Ala-Pro-p-nitroanilide (Sigma-Aldrich (Shanghai) Trading Co., Ltd., Shanghai, China) following the method described by Tereshchenkova et al. [3]. The substrate was initially dissolved in dimethylformamide, and its final reaction concentration was 0.25 mM. SgVnDPPIV (1 µg) was combined with 2.5 µL substrate. The freshly prepared PBS (phosphate-buffered saline) (0.1 M, PH 8.0) was added to the final volume of 200 µL. The mixture was incubated at 40 °C for 30 min. Its absorbance was measured periodically using an Infinite F50 Plus and Infinite F50 Robotic microplate reader (TECAN, Männedorf, Switzerland) at 405 nm for 60 min. Venom of S. guani and gut extract from the T. molitor larvae were used as the positive controls. The effects of the inhibitors, including vildagliptin, sitagliptin, diprotin A, and diprotin B (Sigma-Aldrich (Shanghai) Trading Co., Ltd., Shanghai, China), on the enzymatic activity of SgVnDPPIV were tested. The inhibitor concentration was 0.1 mM in the incubate assay with the SgVnDPPIV before its reaction with the substrate, as described above. All assays were performed in three independent biological replicates. The enzymatic activity was expressed in units (U).

5.6. Transcriptomic Sequencing and Analysis

SgVnDPPIV (0.25 µg in 2 µL PBS, PH 7.4) was microinjected into the pupae of T. molitor after larvae molt into pupae 24 h. Pupae injected with equal volume of PBS were set as controls. Subsequently, 6 h and 24 h after treatment, pupae were collected and grinded in liquid nitrogen with a mortar and pestle. Samples from at least five pupae were pooled together and set as one biological replicate. Three biological replicates were used for each treatment at each time point. Total RNA was extracted using Trizol reagent (TaKaRa, Dalian, China). Its quantity and quality were determined with the Agilent 2100 Bioanalyzer (Agilent Technologies, technologies, Santa Clara, CA, USA) and using 1% agarose gel electrogenesis. mRNA was isolated from the total RNA with oligo d(T) attached magnetic beads, and was used to construct the library using the TruSeq RNA Access Library Prep Kit (Illumina, San Diego, CA, USA). Libraries were sequenced on an Illumina HiSeq 2000 platform (Illumina, San Diego, CA, USA) and 150 bp paired-end reads were generated. Raw data were deposited into the National Center for Biotechnology Information (NCBI) Sequence Read Archive under BioProject accession PRJNA954966 and PRJNA955130. After quality control to remove residual adapters and low-quality nucleotides, clean reads were mapped to the gene repertoire encoded in the genome of T. molitor using the Kallisto Super Wrapper function in TBtools v1.0987663 [44]. This gene repertoire was generated by combining genes annotated in different genome versions of T. molitor sequenced by our laboratory [49] and two other groups [50,51], followed by reducing the redundant sequences using CD-HIT [52]. The number of mapped reads to each gene was used to measure the gene expression in TPM using DESeq2 [53]. The final set of differentially expressed genes was defined as those with a significant adjusted p-value below 0.05 and an absolute fold change above 1.5. Gene ontology (GO) enrichment for the differentially expressed genes with blast hits of an e-value ≤ 1 × 10−5 against the NCBI Nr protein database was performed with the GO Enrichment function in TBtools v1.0987663 [44,54]. Significantly enriched GO categories are defined with an adjusted p-value < 0.05.

5.7. Statistical Analysis

All data were analyzed by one-way analysis of variance (ANOVA) with the least significant difference (LSD) method using Data Processing System (DPS) software v20.00 [55]. Statistical significance was set at the p < 0.05 level.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/toxins15050311/s1. Figure S1: The nucleotide and deduced amino acid sequences of the SgVnDPPIV gene. Figure S2: Western blot analysis of the expression of recombinant SgVnDPPIV in Sf9 cells using the baculovirus expression system (1, cell supernatant; 2, cell lysis supernatant). Figure S3: Gene ontology (GO) enrichment analysis at three levels of differentially expressed genes in Tenebrio molitor pupae envenomated by SgVnDPPIV (A, enriched GO terms and differentially expressed gene number 6 h post SgVnDPPIV injection; B, enriched GO terms and differentially expressed gene number 24 h post SgVnDPPIV injection). Table S1: Number of DPPIV genes encoding in the genomes of various hymenopteran species. Table S2: Amino acid sequences of DPPIVs identified in the genomes of various hymenopteran species. Table S3: Statistics of the transcriptomic sequencing. Table S4: Detailed list of differentially expressed genes in Tenebrio molitor pupae envenomated by SgVnDPPIV. Table S5: Detailed results of the gene ontology (GO) enrichment analysis of the differentially expressed genes in Tenebrio molitor pupae envenomated by SgVnDPPIV.

Author Contributions

Investigation, methodology, software, data curation, formal analysis, validation, visualization, writing—original draft, C.W.; investigation, methodology, software, data curation, formal analysis, visualization, C.Y.; software, data curation, formal analysis, validation, visualization, Y.W.; software, formal analysis, visualization, J.W.; conceptualization, resources, supervision, project administration, funding acquisition, writing—original draft, review, and editing, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Science and Technology Innovation Talent Program in Forestry and Grassland for Young Top-Notch Talents (2019132615), the National Natural Science Foundation of China (32060126), the CAS “Light of West China” Program (Western Youth Scholars “A”), the Xing Dian Talents Support Program of Yunnan Province (YNWR-QNBJ-2018-393), and the Open Fund of State Key Laboratory for Conservation and Utilization of Bio-Resources in Yunnan to J.Z.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the information is at this document.

Conflicts of Interest

The authors declare no conflict of interests.

References

  1. Rawlings, N.D.; Waller, M.; Barrett, A.J.; Bateman, A. MEROPS: The database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2014, 42, D503–D509. [Google Scholar] [CrossRef] [PubMed]
  2. Borloo, M.; De Meester, I. Dipeptidyl peptidase IV: Development, design, synthesis and biological evaluation of inhibitors. Verh. K. Acad. Geneeskd. Belg. 1994, 56, 57–88. [Google Scholar] [PubMed]
  3. Tereshchenkova, V.F.; Goptar, I.A.; Kulemzina, I.A.; Zhuzhikov, D.P.; Serebryakova, M.V.; Belozersky, M.A.; Dunaevsky, Y.E.; Oppert, B.; Filippova, I.Y.; Elpidina, E.N. Dipeptidyl peptidase 4—An important digestive peptidase in Tenebrio molitor larvae. Insect Biochem. Mol. Biol. 2016, 76, 38–48. [Google Scholar] [CrossRef] [PubMed]
  4. Mulvihill, E.E.; Drucker, D.J. Pharmacology, physiology, and mechanisms of action of dipeptidyl peptidase-4 inhibitors. Endocr. Rev. 2014, 35, 992–1019. [Google Scholar] [CrossRef]
  5. Yang, Q.; Fu, B.; Luo, D.; Wang, H.B.; Cao, H.Y.; Chen, X.; Tian, L.; Yu, X.J. The multiple biological functions of dipeptidyl peptidase-4 in bone metabolism. Front. Endocrinol. 2022, 13, 856954. [Google Scholar] [CrossRef] [PubMed]
  6. Juillerat-Jeanneret, L. Dipeptidyl peptidase IV and its inhibitors: Therapeutics for type 2 diabetes and what else? J. Med. Chem. 2014, 57, 2197–2212. [Google Scholar] [CrossRef]
  7. Dávalos Terán, I.; Imai, K.; Lacroix, I.M.E.; Fogliano, V.; Udenigwe, C.C. Bioinformatics of edible yellow mealworm (Tenebrio molitor) proteome reveal the cuticular proteins as promising precursors of dipeptidyl peptidase-IV inhibitors. J. Food Biochem. 2020, 44, e13121. [Google Scholar] [CrossRef]
  8. Przezak, A.; Bielka, W.; Pawlik, A. Incretins in the therapy of diabetic kidney disease. Int. J. Mol. Sci. 2021, 22, 12312. [Google Scholar] [CrossRef]
  9. Martensen, I.; Koolman, J.; Mentlein, R. Proline-specific dipeptidyl peptidase from the blue blowfly Calliphora vicina hydrolyzes in vitro the ecdysiostatic peptide trypsin-modulating oostatic factor (Neb-TMOF). Arch. Insect Biochem. Physiol. 1998, 37, 146–157. [Google Scholar] [CrossRef]
  10. Nässel, D.R.; Mentlein, R.; Bollner, T.; Karlsson, A. Proline-specific dipeptidyl peptidase activity in the cockroach brain and intestine: Partial characterization, distribution, and inactivation of tachykinin-related peptides. J. Comp. Neurol. 2000, 418, 81–92. [Google Scholar] [CrossRef]
  11. Chihara, C.J.; Song, C.Y.; LaMonte, G.; Fetalvero, K.; Hinchman, K.; Phan, H.; Pineda, M.; Robinson, K.; Schneider, G.P. Identification and partial characterization of the enzyme of omega: One of five putative DPP IV genes in Drosophila melanogaster. J. Insect Sci. 2005, 5, 26. [Google Scholar] [PubMed]
  12. Blank, S.; Seismann, H.; Bockisch, B.; Braren, I.; Cifuentes, L.; McIntyre, M.; Rühl, D.; Ring, J.; Bredehorst, R.; Ollert, M.W.; et al. Identification, recombinant expression, and characterization of the 100 kDa high molecular weight Hymenoptera venom allergens Api m 5 and Ves v 3. J. Immunol. 2010, 184, 5403–5413. [Google Scholar] [PubMed]
  13. Jang, S.; Chung, T.Y.; Shin, J.; Lin, K.L.; Tzen, J.T.; Li, F.Y. Docking study of the precursor peptide of mastoparan onto its putative processing enzyme, dipeptidyl peptidase IV: A revisit to molecular ticketing. J. Comput. Aided Mol. Des. 2010, 24, 213–224. [Google Scholar] [CrossRef] [PubMed]
  14. Lee, V.S.Y.; Tu, W.C.; Jinn, T.R.; Peng, C.C.; Lin, L.J.; Tzen, J.T.C. Molecular cloning of the precursor polypeptide of mastoparan B and its putative processing enzyme, dipeptidyl peptidase IV, from the black-bellied hornet, Vespa basalis. Insect Mol. Biol. 2007, 16, 231–237. [Google Scholar] [CrossRef]
  15. Pretre, G.; Asturias, J.A.; Lizaso, M.T.; Tabar, A.I. Dipeptidyl peptidase IV of the Vespa velutina nigrithorax venom is recognized as a relevant allergen. Ann. Allergy Asthma Immunol. 2022, 129, 101–105. [Google Scholar]
  16. Schiener, M.; Hilger, C.; Eberlein, B.; Pascal, M.; Kuehn, A.; Revets, D.; Planchon, S.; Pietsch, G.; Serrano, P.; Moreno-Aguilar, C.; et al. The high molecular weight dipeptidyl peptidase IV Pol d 3 is a major allergen of Polistes dominula venom. Sci. Rep. 2018, 8, 1318. [Google Scholar]
  17. Yu, H.N.; Yang, H.L.; Ma, D.Y.; Lv, Y.; Liu, T.G.; Zhang, K.Y.; Lai, R.; Liu, J.Z. Vespid chemotactic peptide precursor from the wasp, Vespa magnifica (Smith). Toxicon 2007, 50, 377–382. [Google Scholar] [CrossRef]
  18. Chen, W.H.; Yang, X.B.; Yang, X.L.; Zhai, L.; Lu, Z.K.; Liu, J.Z.; Yu, H.N. Antimicrobial peptides from the venoms of Vespa bicolor Fabricius. Peptides 2008, 29, 1887–1892. [Google Scholar] [CrossRef]
  19. Mrinalini; Werren, J.H. Parasitoid wasps and their venoms. In Evolution of Venomous Animals and Their Toxins; Starr, C., Ed.; Toxinology; Springer: Dordrecht, The Netherlands, 2017; pp. 1–26. [Google Scholar]
  20. Beckage, N.E.; Gelman, D.B. Wasp parasitoid disruption of host development: Implications for new biologically based strategies for insect control. Annu. Rev. Entomol. 2004, 49, 299–330. [Google Scholar]
  21. Moreau, S.J.M.; Asgari, S. Venom proteins from parasitoid wasps and their biological functions. Toxins 2015, 7, 2385–2412. [Google Scholar] [CrossRef]
  22. Gatti, J.L.; Belghazi, M.; Legeai, F.; Ravallec, M.; Frayssinet, M.; Robin, S.; Aboubakar-Souna, D.; Srinivasan, R.; Tamò, M.; Poirié, M.; et al. Proteo-trancriptomic analyses reveal a large expansion of metalloprotease-like proteins in atypical venom vesicles of the wasp Meteorus pulchricornis (Braconidae). Toxins 2021, 13, 502. [Google Scholar] [PubMed]
  23. Moreau, S.J.M.; Guillot, S. Advances and prospects on biosynthesis, structures and functions of venom proteins from parasitic wasps. Insect Biochem. Mol. Biol. 2005, 35, 1209–1223. [Google Scholar] [CrossRef] [PubMed]
  24. Yang, Y.; Ye, X.H.; Dang, C.; Cao, Y.S.; Hong, R.; Sun, Y.H.; Xiao, S.; Mei, Y.; Xu, L.; Fang, Q.; et al. Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation. BMC Biol. 2021, 19, 145. [Google Scholar] [CrossRef] [PubMed]
  25. Dorémus, T.; Urbach, S.; Jouan, V.; Cousserans, F.; Ravallec, M.; Demettre, E.; Wajnberg, E.; Poulain, J.; Azéma-Dossat, C.; Darboux, I.; et al. Venom gland extract is not required for successful parasitism in the polydnavirus-associated endoparasitoid Hyposoter didymator (Hym. Ichneumonidae) despite the presence of numerous novel and conserved venom proteins. Insect Biochem. Mol. Biol. 2013, 43, 292–307. [Google Scholar] [PubMed]
  26. Liu, N.Y.; Xu, Z.W.; Yan, W.; Ren, X.M.; Zhang, Z.Q.; Zhu, J.Y. Venomics reveals novel ion transport peptide-likes (ITPLs) from the parasitoid wasp Tetrastichus brontispae. Toxicon 2018, 141, 88–93. [Google Scholar] [CrossRef]
  27. Chen, J.; Cheng, H.Z. Advances in applied research on Scleroderma spp. Chin. J. Biol. Control 2000, 16, 166–170. [Google Scholar]
  28. Li, X.R.; Lu, D.G.; Liu, X.X.; Zhang, Q.W.; Zhou, X.G. Ultrastructural characterization of olfactory sensilla and immunolocalization of odorant binding and chemosensory proteins from an ectoparasitoid Scleroderma guani (Hymenoptera: Bethylidae). Int. J. Biol. Sci. 2011, 7, 848–868. [Google Scholar] [CrossRef] [PubMed]
  29. Zhu, J.Y. Deciphering the main venom components of the ectoparasitic ant-like bethylid wasp, Scleroderma guani. Toxicon 2016, 113, 32–40. [Google Scholar]
  30. Martinson, E.O.; Mrinalini, W.J.H.; Kelkar, Y.D.; Chang, C.H.; Werren, J.H. The evolution of venom by co-option of single-copy genes. Curr. Biol. 2017, 27, 2007–2013.e8. [Google Scholar]
  31. Ogawa, Y.; Mamura, Y.; Murayama, N.; Yanoshita, R. Characterization and cDNA cloning of dipeptidyl peptidase IV from the venom of Gloydius blomhoffi brevicaudus. Comp. Biochem. Physiol. 2006, 145, 35–42. [Google Scholar] [CrossRef]
  32. Qian, C.; Fang, Q.; Wang, L.; Ye, G.Y. Molecular cloning and functional studies of two kazal-type serine protease inhibitors specifically expressed by Nasonia vitripennis venom apparatus. Toxins 2015, 7, 2888–2905. [Google Scholar] [CrossRef] [PubMed]
  33. Misumi, Y.; Hayashi, Y.; Arakawa, F.; Ikehara, Y. Molecular cloning and sequence analysis of human dipeptidyl peptidase IV, a serine proteinase on the cell surface. Biochim. Biophys. Acta 1992, 1131, 333–336. [Google Scholar] [CrossRef] [PubMed]
  34. Meadows, S.A.; Edosada, C.Y.; Mayeda, M.; Tran, T.; Quan, C.; Raab, H.; Wiesmann, C.; Wolf, B.B. Ala657 and conserved active site residues promote fibroblast activation protein endopeptidase activity via distinct mechanisms of transition state stabilization. Biochemistry 2007, 46, 4598–4605. [Google Scholar] [CrossRef]
  35. Nabeno, M.; Akahoshi, F.; Kishida, H.; Miyaguchi, I.; Tanaka, Y.; Ishii, S.; Kadowaki, T. A comparative study of the binding modes of recently launched dipeptidyl peptidase IV inhibitors in the active site. Biochem. Biophys. Res. Commun. 2013, 434, 191–196. [Google Scholar] [CrossRef] [PubMed]
  36. Lee, S.H.; Baek, J.H.; Yoon, K.A. Differential properties of venom peptides and proteins in solitary vs. social hunting wasps. Toxins 2016, 8, 32. [Google Scholar] [CrossRef]
  37. Wertheim, B. Adaptations and counter-adaptations in Drosophila host-parasitoid interactions: Advances in the molecular mechanisms. Curr. Opin. Insect Sci. 2022, 51, 100896. [Google Scholar] [CrossRef]
  38. Yoon, K.A.; Kim, K.; Kim, W.J.; Bang, W.Y.; Ahn, N.H.; Bae, C.H.; Yeo, J.H.; Lee, S.H. Characterization of venom components and their phylogenetic properties in some aculeate bumblebees and wasps. Toxins 2020, 12, 47. [Google Scholar] [CrossRef]
  39. Becchimanzi, A.; Avolio, M.; Bostan, H.; Colantuono, C.; Cozzolino, F.; Mancini, D.; Chiusano, M.L.; Pucci, P.; Caccia, S.; Pennacchio, F. Venomics of the ectoparasitoid wasp Bracon nigricans. BMC Genom. 2020, 21, 34. [Google Scholar] [CrossRef]
  40. Dos Santos-Pinto, J.R.A.; Perez-Riverol, A.; Lasa, A.M.; Palma, M.S. Diversity of peptidic and proteinaceous toxins from social Hymenoptera venoms. Toxicon 2018, 148, 172–196. [Google Scholar] [CrossRef]
  41. Martinson, E.O.; Siebert, A.L.; He, M.; Kelkar, Y.D.; Doucette, L.A.; Werren, J.H. Evaluating the evolution and function of the dynamic venom Y protein in ectoparasitoid wasps. Insect Mol. Biol. 2019, 28, 499–508. [Google Scholar] [CrossRef]
  42. Zhu, J.Y.; Yang, P.; Zhang, Z.; Wu, G.X.; Yang, B. Transcriptomic immune response of Tenebrio molitor pupae to parasitization by Scleroderma guani. PLoS ONE 2013, 8, e54411. [Google Scholar] [CrossRef] [PubMed]
  43. Teufel, F.; Almagro Armenteros, J.J.; Johansen, A.R.; Gíslason, M.H.; Pihl, S.I.; Tsirigos, K.D.; Winther, O.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat. Biotechnol. 2022, 40, 1023–1025. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.H.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  45. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef]
  46. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef]
  47. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  48. Xiong, A.S.; Yao, Q.H.; Peng, R.H.; Duan, H.; Li, X.; Fan, H.Q.; Cheng, Z.M.; Li, Y. PCR-based accurate synthesis of long DNA sequences. Nat. Protoc. 2006, 1, 791–797. [Google Scholar] [CrossRef]
  49. Zhu, J.Y.; Song, Q.S.; Stanley, D. Genomic and transcriptomic analysis of Tenebrio molitor: Preface. Arch. Insect Biochem. Physiol. 2022, 111, e21969. [Google Scholar] [CrossRef]
  50. Eleftheriou, E.; Aury, J.M.; Vacherie, B.; Istace, B.; Orcid, C.B.; Noel, B.; Moret, Y. Chromosome-scale assembly of the yellow mealworm genome. Open Res. Eur. 2022, 1, 94. [Google Scholar] [CrossRef]
  51. Kaur, S.; Stinson, S.A.; DiCenzo, G.C. Whole genome assemblies of Zophobas morio and Tenebrio molitor. G3 Genes 2023, jkad079. [Google Scholar] [CrossRef]
  52. Li, W.Z.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [PubMed]
  53. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  54. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef] [PubMed]
  55. Tang, Q.Y.; Zhang, C.X. Data Processing System (DPS) software with experimental design, statistical analysis and data mining developed for use in entomological research. Insect Sci. 2013, 20, 254–260. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Multiple alignment of SgVnDPPIV and DPPIVs from other insects and humans. Sgua, Scleroderma guani; Tmol, Tenebrio molitor (CAH1368919); Vbas, Vespa basalis (ABG57265); Mdem, Microplitis demolitor (XP_008553451); Hsap, Homo sapiens (P27487). The transmembrane domains of S. guani, V. basalis, M. demolitor, and H. sapiens DPPIVs are marked in red boxes. The signal peptides of T. molitor and V. basalis DPPIVs are indicated in green boxes. The catalytic triads and substrate binding sites are indicated with triangles and solid circles, respectively. The identical and highly conserved residues are shaded in black and grey, respectively.
Figure 1. Multiple alignment of SgVnDPPIV and DPPIVs from other insects and humans. Sgua, Scleroderma guani; Tmol, Tenebrio molitor (CAH1368919); Vbas, Vespa basalis (ABG57265); Mdem, Microplitis demolitor (XP_008553451); Hsap, Homo sapiens (P27487). The transmembrane domains of S. guani, V. basalis, M. demolitor, and H. sapiens DPPIVs are marked in red boxes. The signal peptides of T. molitor and V. basalis DPPIVs are indicated in green boxes. The catalytic triads and substrate binding sites are indicated with triangles and solid circles, respectively. The identical and highly conserved residues are shaded in black and grey, respectively.
Toxins 15 00311 g001
Figure 2. Phylogenetic tree of DPPIVs from various hymenopteran species. Amino acid sequences of DPPIVs used to construct the phylogenetic tree are listed in Table S2.
Figure 2. Phylogenetic tree of DPPIVs from various hymenopteran species. Amino acid sequences of DPPIVs used to construct the phylogenetic tree are listed in Table S2.
Toxins 15 00311 g002
Figure 3. Gene expression profiles of SgVnDPPIV detected using quantitative real-time PCR. (A) Expression profiles of SgVnDPPIV gene in various tissues of female adults. He, head; Th, thorax; Ab, abdomen derived of venom apparatus; VA, venom apparatus. (B) Expression profiles of SgVnDPPIV gene at different developmental stages. E, egg; EL, early larvae; LL, late larvae; ML, mature larvae; SL, spinning mature larvae; WP, pupae in white cocoon; YP, pupae in yellow cocoon; BP, pupae in back cocoon; A1–25, adults after 1–25 days emergence. Different superscript letters indicate statistically significant differences among different samples.
Figure 3. Gene expression profiles of SgVnDPPIV detected using quantitative real-time PCR. (A) Expression profiles of SgVnDPPIV gene in various tissues of female adults. He, head; Th, thorax; Ab, abdomen derived of venom apparatus; VA, venom apparatus. (B) Expression profiles of SgVnDPPIV gene at different developmental stages. E, egg; EL, early larvae; LL, late larvae; ML, mature larvae; SL, spinning mature larvae; WP, pupae in white cocoon; YP, pupae in yellow cocoon; BP, pupae in back cocoon; A1–25, adults after 1–25 days emergence. Different superscript letters indicate statistically significant differences among different samples.
Toxins 15 00311 g003
Figure 4. Enzymatic assay of recombinant SgVnDPPIV. (A) Expression of SgVnDPPIV in Sf9 cells using the baculovirus expression system. The recombinant SgVnDPPIV was purified and detected with Coomassie Brilliant Blue. (B) Enzymatic activity of gut extract of Tenebrio molitor larvae, venom of Scleroderma guani, and recombinant SgVnDPPIV. (C) Effect of inhibitors on the activity of recombinant SgVnDPPIV. Different superscript letters indicate statistically significant differences among different samples or treatments.
Figure 4. Enzymatic assay of recombinant SgVnDPPIV. (A) Expression of SgVnDPPIV in Sf9 cells using the baculovirus expression system. The recombinant SgVnDPPIV was purified and detected with Coomassie Brilliant Blue. (B) Enzymatic activity of gut extract of Tenebrio molitor larvae, venom of Scleroderma guani, and recombinant SgVnDPPIV. (C) Effect of inhibitors on the activity of recombinant SgVnDPPIV. Different superscript letters indicate statistically significant differences among different samples or treatments.
Toxins 15 00311 g004
Figure 5. Differentially expressed genes in Tenebrio molitor pupae envenomated by SgVnDPPIV. (A) Volcano plot of differentially expressed genes identified 6 h post SgVnDPPIV injection. (B) Volcano plot of differentially expressed genes identified 24 h post SgVnDPPIV injection. (C) Venn diagram of differentially expressed genes identified 6 h and 24 h post SgVnDPPIV injection.
Figure 5. Differentially expressed genes in Tenebrio molitor pupae envenomated by SgVnDPPIV. (A) Volcano plot of differentially expressed genes identified 6 h post SgVnDPPIV injection. (B) Volcano plot of differentially expressed genes identified 24 h post SgVnDPPIV injection. (C) Venn diagram of differentially expressed genes identified 6 h and 24 h post SgVnDPPIV injection.
Toxins 15 00311 g005
Figure 6. Gene ontology (GO) enrichment analysis at a broad level of differentially expressed genes in Tenebrio molitor pupae envenomated by SgVnDPPIV. (A) Enriched GO terms and differentially expressed gene number 6 h post SgVnDPPIV injection. (B) Enriched GO terms and differentially expressed gene number 24 h post SgVnDPPIV injection.
Figure 6. Gene ontology (GO) enrichment analysis at a broad level of differentially expressed genes in Tenebrio molitor pupae envenomated by SgVnDPPIV. (A) Enriched GO terms and differentially expressed gene number 6 h post SgVnDPPIV injection. (B) Enriched GO terms and differentially expressed gene number 24 h post SgVnDPPIV injection.
Toxins 15 00311 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, C.; Yang, C.; Wang, Y.; Wang, J.; Zhu, J. Molecular Characterization and Functional Analysis of the Dipeptidyl Peptidase IV from Venom of the Ectoparasitoid Scleroderma guani. Toxins 2023, 15, 311. https://doi.org/10.3390/toxins15050311

AMA Style

Wu C, Yang C, Wang Y, Wang J, Zhu J. Molecular Characterization and Functional Analysis of the Dipeptidyl Peptidase IV from Venom of the Ectoparasitoid Scleroderma guani. Toxins. 2023; 15(5):311. https://doi.org/10.3390/toxins15050311

Chicago/Turabian Style

Wu, Chaoyan, Cheng Yang, Yuqin Wang, Jun Wang, and Jiaying Zhu. 2023. "Molecular Characterization and Functional Analysis of the Dipeptidyl Peptidase IV from Venom of the Ectoparasitoid Scleroderma guani" Toxins 15, no. 5: 311. https://doi.org/10.3390/toxins15050311

APA Style

Wu, C., Yang, C., Wang, Y., Wang, J., & Zhu, J. (2023). Molecular Characterization and Functional Analysis of the Dipeptidyl Peptidase IV from Venom of the Ectoparasitoid Scleroderma guani. Toxins, 15(5), 311. https://doi.org/10.3390/toxins15050311

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop