Next Article in Journal
Impacts of Molecular Structure on Nucleic Acid–Protein Interactions
Previous Article in Journal
Monochromatic Green Light Stimulation during Incubation Alters Hepatic Glucose Metabolism That Improves Embryonic Development in Yangzhou Goose Eggs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prophenoloxidase of Odontotermes formosanus (Shiraki) (Blattodea: Termitidae) Is a Key Gene in Melanization and Has a Defensive Role during Bacterial Infection

1
Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
2
College of Forestry, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(1), 406; https://doi.org/10.3390/ijms24010406
Submission received: 27 September 2022 / Revised: 22 December 2022 / Accepted: 22 December 2022 / Published: 26 December 2022
(This article belongs to the Section Biochemistry)

Abstract

:
Melanization mediated by the prophenoloxidase (PPO)-activating system is an important innate immunity to fight pathogens in insects. In this study, the in vitro time-dependent increase in the intensity of melanization and phenoloxidase (PO) activity from the hemolymph of Odontotermes formosanus (Shiraki) challenged by pathogenic bacteria was detected. PPO is one of the key genes in melanization pathway, whereas the molecular characteristics and functions of O. formosanus PPO are unclear. The OfPPO gene was cloned and characterized. The open reading frame of OfPPO is 2085 bp in length and encodes a 79.497 kDa protein with 694 amino acids. A BLASTx search and phylogenetic analyses revealed that OfPPO shares a high degree of homology to the Blattodea PPOs. Moreover, real-time fluorescent quantitative PCR analysis showed that OfPPO is ubiquitously expressed in all castes and tissues examined, with the highest expression in workers and variable expression patterns in tissues of different termite castes. Furthermore, the expression of OfPPO was significantly induced in O. formosanus infected by pathogenic bacteria. Intriguingly, in combination with silencing of OfPPO expression, pathogenic bacteria challenge caused greatly increased mortality of O. formosanus. These results suggest that OfPPO plays a role in defense against bacteria and highlight the novel termite control strategy combining pathogenic bacteria application with termite PPO silencing.

1. Introduction

Termites are important agricultural pests that can damage to nearly anything made of cellulose, leading to an annual loss of 40 billion dollars worldwide [1,2]. They are eusocial insects with definite castes, such as king, queen, soldiers, and workers, within colonies. The Odontotermes formosanus (Shiraki) (Blattodea: Termitidae) is a notorious termite species that is distributed in China, India, Japan, Myanmar, Vietnam, and Thailand [3,4]. O. formosanus not only damages crops, plantations, and forests but is also a major concern for disasters if they build numerous hidden chambers in dams, leading to structural instability [5,6]. O. formosanus populations can be controlled by applications of pesticides, but the excessive use of pesticides can seriously affect the human health and environment safety [7]. Therefore, there is an urgent need to develop environmentally friendly biological control. Serratia marcescens and Bacillus thuringiensis have promising applications in the biological control of termites [8,9,10,11,12]. S. marcescens, a pathogenic bacterium recently isolated from dead O. formosanus, has a pathogenic effect on O. formosanus [13]. Laboratory studies have demonstrated that B. thuringiensis, the mostly and widely used bacterial insecticide, has a good control effect on termite [11,12]. Thus, these two insect pathogens have great potential to be good microbial control pathogens against O. formosanus.
Insects combat infection by pathogens through mounting powerful immune responses (cellular and humoral immunity), which reduce the virulence of pathogens, such as apoptosis, phagocytosis, autophagy, melanization, the JAK/STAT pathway, the Toll pathway and the immune deficiency (Imd) pathway, etc. [14,15,16,17]. Currently, immune genes are regarded as potential target genes for RNA interference (RNAi)-based control strategies. The RNAi technique has potential practical value for the development of new tools for the management of insect pests [18,19]. For example, the Imd pathway is involved in the antifungal response of Anthonomus grandis and interference with the relish gene by RNAi can increase the susceptibility of A. grandis to Metarhizium anisopliae [20]. Inhibiting apoptosis related genes of Colorado potato beetle is an effective RNAi-mediated control measure [21]. The silencing of the locust apoptosis-related gene can reduce insect immunity and increase insect susceptibility to fungal pathogens [22]. The knockdown of several genes in the Toll pathway also can increase pest lethality, which provides new targets for pest control [23].
In insects, melanization reaction is one of the most important immune processes. Mounting evidence suggests that the antimicrobial effect of melanization is an indispensable component of the insect defense mechanism to protect the insects against microbes [24,25]. Phenoloxidase (PO) catalyzes the conversion of phenols to quinones via an oxidation reaction and promotes the formation of melanin which directly encapsulates and clears away foreign pathogens [26,27,28]. Active PO is formed by prophenoloxidase (PPO) in a specific protease system [29,30]. Previous research suggested that the PPO-mediated melanization pathway was indispensable for insects to defend against pathogens, in which PPO played a critical role in regulating immune responses [31]. According to the information mentioned above, interfering with melanization pathway by RNAi technology can be an effective method to control termites. The potentially critical target gene is PPO in the RNAi-mediated control strategy. However, its characteristics and function have not been defined yet.
In this study, we used pathogenic bacteria and O. formosanus as an infection model to investigate the effects of pathogenic bacterial infection on the hemolymph melanization of O. formosanus. We cloned OfPPO gene and analyzed the function of OfPPO in O. formosanus. OfPPO was determined to serve as an important component of hemolymph melanization and to defend insects against pathogenic bacteria. This study not only characterizes the function of PPO gene of O. formosanus, but also provides insights for the integrated control of O. formosanus with RNAi and pathogenic bacteria mediated control strategies.

2. Results

2.1. Melanization Accelerated and PO Activity Induced by SM1 and Bt in O. formosanus

Melanization of hemolymph from O. formosanus infected with SM1 or Bt was investigated. As shown in Figure 1A, compared to the uninfected groups, accelerated melanization was observed in O. formosanus treated by SM1 and Bt for 0, 1, 2, 3, and 4 h. Moreover, the elevated PO activity occurred in the bacteria-infected groups at most of the time points post treatment (Figure 1B,C).

2.2. Cloning and Characterization of OfPPO in O. formosanus

We successfully cloned the complete cDNA of the OfPPO gene (accession number: OP828923) and found that OfPPO contains an ORF of 2085 bp encoding a 694 amino acid protein with a molecular weight of 79.497 kDa and a predicted isoelectric point (pI) of 6.01. Alignment of amino acid sequences of OfPPO with the other nine insect PPOs showed that six motifs were extremely well-conserved in OfPPO (Figure 2). Among the six motifs, two were annotated as potential proteolytic cleavage sites, another two were copper binding sites (CuA and CuB), one was a thiol-ester-like motif, and the other was a conserved C-terminal motif. Comparing the amino acid sequence using BLASTx, the OfPPO protein exhibited high degree of identity ranged between 59.68% and 93% with the PPO of other insect species (Table 1). The OfPPO protein sequence is highly similar to insect PPO sequences of the Blattella. OfPPO shares 92.8% pairwise identity with Coptotermes formosanus PO2 (CfPO2) at the amino acid level and 87.9% and 79.34% pairwise identity with Zootermopsis nevadensis PO2 (ZnPO2) and Blattella germanica PO1 (BgPO1), respectively (Table 1). Based on the deduced OfPPO amino acid sequence, 81 invertebrate PPO (PO) sequences were selected, and a high-accuracy phylogenetic tree was generated using MEGA X. The phylogenetic tree was divided into six different groups. In particular, OfPPO was most clustered with Blattodea PPO (Figure 3). This result showed that OfPPO was highly similar to the PPO of Blattodea.

2.3. Caste- and Tissue-Specific Expression Profiles of OfPPO in O. formosanus

The mRNA levels of OfPPO from different castes of O. formosanus were analyzed using real-time fluorescent quantitative PCR (qPCR). The expression of OfPPO in the workers was significantly higher than that in soldiers, larval instars, and dealates (Figure 4A).
The expression of OfPPO was examined in different body parts as well. For female dealates, OfPPO was expressed at relatively high levels in the hemolymph and head (Figure 4B). For male dealates, OfPPO was expressed at the highest level in the head followed by the hemolymph (Figure 4C). For workers, the expression pattern of OfPPO was similar to that of male dealates (Figure 4D). For larval instars, OfPPO was expressed at high levels in the hemolymph and gut (Figure 4E). For soldiers, OfPPO was expressed at the highest level in the gut, followed by the leg (Figure 4F).

2.4. The OfPPO Expression Induced by SM1 and Bt in O. formosanus

Next, we investigated whether OfPPO is involved in defense against bacterial infection. We measured the mRNA levels of OfPPO in O. formosanus challenged by SM1 and Bt, respectively. As shown in Figure 5A, the OfPPO expression was significantly increased after SM1 infection at all treated time points. Moreover, the transcript levels of OfPPO were significantly upregulated when treated by Bt for 6, 12, and 24 h, respectively (Figure 5B). In total, the results showed that SM1 and Bt challenges could induce the expression of OfPPO in O. formosanus.

2.5. The Function of OfPPO in Resisting SM1 and Bt Infection in O. formosanus

To evaluate the effectiveness of two dsRNAs from two important OfPPO domains, two in vitro synthesized dsRNAs (dsOfPPO1 and dsOfPPO2) were designed (Figure 6A). We found that dsOfPPO1 and dsOfPPO2 suppressed the transcript levels of OfPPO by 92.7% and 85.1% at 6 h, respectively (Figure 6B,C).
The effect of OfPPO on bacterial infection was examined in vivo. After feeding O. formosanus for 6 h with sterilized water (negative control), GFP dsRNA, and dsOfPPO1 or dsOfPPO2, O. formosanus was exposed to SM1. The mortality of O. formosanus treated with dsOfPPO1 and SM1 was significantly increased compared with that of O. formosanus treated with dsGFP and SM1 at 12 and 24 h (Figure 7B,C). Interestingly, the mortality of O. formosanus treated with dsOfPPO1 and SM1 was almost 100% at 48 h (Figure 7D). In addition, the mortality of O. formosanus treated with dsOfPPO1 and Bt was significantly increased at 24, 48, and 72 h compared with that of O. formosanus treated with dsGFP and Bt (Figure 8A–C). Furthermore, the mortality of O. formosanus treated with dsOfPPO1 and Bt was almost 100% at 96 h (Figure 8D). Similar results were obtained when O. formosanus was treated with dsOfPPO2 and SM1 or dsOfPPO2 and Bt (Figures S1 and S2).

3. Discussion

The development of entomopathogens as classical, safe, eco-friendly, and sustainable biological control agents has achieved many successes in recent years [32]. Microbial control techniques involving entomopathogens have been developed as a nonchemical alternative for pest control [33]. Meanwhile, entomopathogens have received attention as sustainable environmentally friendly control agents for the control of the insects [8,9,12]. To combat pathogen infection, insects have evolved potent immune defense mechanisms, including cellular immunity and humoral immunity. In insect immunity, melanization is an important part of humoral immunity. It is well known that the melanization reaction is an important anti-pathogen strategy that attacks pathogenic bacteria in insects, thus reducing the damage caused by microorganisms to insects [29,34,35,36]. Growing evidence suggests that different insects have developed PPO-mediated melanization pathways to defend against infection [35,37,38]. In this study, we observed that melanization in O. formosanus was accelerated by SM1 and Bt infection. It demonstrated that termite melanization is activated in response to biocontrol bacterial infection, which is one of the major obstacles to the use of pathogenic bacteria for termite control. Thus, suppressing the melanization of termites can improve the effectiveness of biocontrol pathogen-mediated pest control.
PPO, as a zymogen, is an important component in the melanization pathway, which produces PO through a specific protease cascade. PO plays an important role in invertebrates, including wound healing [39], hemolymph clotting [40], and encapsulation or melanization of foreign pathogens [17,27]. In Aedes aegypti, the expression of PPOs could be induced after pathogen infection, and the melanization pathway was activated [41,42]. In Spodoptera exigua larval hemolymph, the expression of PPO was significantly induced, and melanin formation occurred after Heliothis virescens ascovirus infection [43]. The pea aphid defended against bacterial and fungal infection through upregulating the expression of PPO and improving the PO activity [30]. PO activity was prominent in the gut of Reticulitermes flavipes, which was useful for maintaining the digestive balance between the insect and the symbiont bacteria [44]. Due to the higher pathogenic pressure in the soil environment, the PO activity of Nasutitermes acajutlae was higher compared to arboreal termites [45]. In this study, we found that the expression of OfPPO and the activity of PO in O. formosanus infected by pathogenic bacteria gradually increased, suggesting that the PPO-mediated melanization has a conserved role in the defense against pathogens in insects.
Very little is known about the characteristics of PPO in O. formosanus. In this research, we first obtained a 2085 bp full-length DNA sequence of OfPPO, which contained specific conserved motifs (proteolytic cleavage sites, conserved copper binding sites, thiol ester-like motifs and C-terminal conserved motifs), which are regarded as expressed sequence tags for PPOs [46,47]. Therefore, we initially speculate that OfPPO has the same function as other insect PPOs, e.g., antibacterial function. In addition, the results of phylogenetic analysis of eighty-one PPO proteins indicated that OfPPO shares a common ancestor with C. formosanus, Z. nevadensis, and B. germanica, as these species cluster on adjacent branches. The results of gene structure and phylogenetic analysis indicated that OfPPO was relatively conserved among different invertebrates. In Apis mellifera, the expression levels of PPO in older pupae and adults were higher than those in larvae and younger pupae [47]. In Hyphantria cunea, the expression of PPO was detected in most stages of development [48]. In Anopheles gambiae, six PPO genes were expressed in different life stages [49]. Here, OfPPO was differentially expressed with respect to castes. Interestingly, workers showed the highest OfPPO expression in all castes. We hypothesized that workers often foraged outside of their colony and suffered diverse stresses, and thus had a higher immunity than soldiers and dealates. With respect to tissue expression, melanization mainly occurred in insect hemolymph, and the expression of OfPPO was relatively high in the hemolymph, which was consistent with a previous report [50,51].
RNAi is a promising new sustainable strategy for ecologically friendly pest control [52]. This technology achieves the inhibition of specific endogenous gene expression by double-stranded RNA (dsRNA) sequences. The effectiveness of RNAi has been proven in many species [52,53]. For example, immune function of Spodoptera littoralis larvae was destroyed by dsRNA molecules, which led to an increase in its susceptibility to Bt [54]. The regulatory genes in the Toll pathway were identified as new lethal targets for pest control [23]. Inhibitors of apoptosis were found to be potent target genes for RNAi-based control of pests [21,22]. In our study, we chose the OfPPO gene, the essential immune gene in melanization, as an RNAi target gene in O. formosanus. When O. formosanus treated by dsOfPPOs was inoculated with SM1 or Bt, the lethality of the pathogenic bacteria to O. formosanus was significantly increased. These results demonstrated that the PPO gene is a good candidate RNAi target for controlling O. formosanus. Some RNAi delivery systems have been developed to control pests, including the use of bacterial, viral, and fungal agents and genetically modified plants or the use of nanoparticles, as well as direct spraying or feeding of dsRNA [22,55]. OfPPO may control O. formosanus via these dsRNA delivery systems.
In summary, PPOs serve important roles in the resistance of insects to pathogenic bacterial infection. Our results defined the OfPPO expression distribution in castes and tissues. RNAi-mediated silencing of OfPPO demonstrated that knockdown of OfPPO expression facilitated pathogenic bacterial infection in O. formosanus. Therefore, PPO genes involved in the melanization reaction can be used as candidate targets for RNAi control.

4. Materials and Methods

4.1. Insects and Bacteria

Odontotermes formosanus colonies were obtained from Nanjing Forestry University (Nanjing, China) and kept in a dark environment at 25 ± 1 °C with 75% relative humidity (RH) that was maintained in a feeding device consisted of colony and foraging areas. The larval instars and soldiers from the colonies were utilized for caste- and tissue-specific gene expression analysis. The naïve female and male dealates were collected for caste- and tissue-specific gene expression analysis. The Serratia marcescens (SM1) and the B. thuringiensis (Bt) strains were stored at −80 °C in our laboratory. For experiments, SM1 was cultured in fermentation medium for 30 h in a shaker incubator with a temperature set at 28 ℃ and a rotation set at 180 rpm. Bt was cultured with lysogen broth under the same conditions for 6 days.

4.2. Melanization Assay and PO Activity

For the melanization assay, twenty O. formosanus workers were placed into culture dishes (9 cm) with moist sterilized filter paper. One microliter of 8.75 × 1010 CFU/mL SM1 or 2.14 × 1010 CFU/mL Bt was placed on each worker’s pronotum (treatment group) (n = 20 workers per replicate) using an injection needle (PB600-1 Repeating Dispenser, Hamilton Co., Reno, NV, USA). For controls, one microliter of culture medium was placed on each worker’s pronotum (n = 20 workers per replicate). The experiment was replicated three times independently. The SM1 and Bt concentrations were predetermined prior to the experiments using dilution coating method. The hemolymph of O. formosanus workers was collected on ice at 24 h post infection (p.i.). The hemolymph and 100 µL PBS were mixed in Eppendorf tubes, then incubated at 28 °C for 4 h and then observed the qualitative level of melanization.
For PO activity assay, the termites of the 3 h treatment groups were treated by one microliter of 8.75 × 1010 CFU/mL SM1 or 2.14 × 1010 CFU/mL Bt for 3 h (n = 20 workers per replicate), while the termites of the 3 h control groups were treated by one microliter of culture medium for 3 h (n = 20 workers per replicate). In addition, the 6 h, 12 h or 24 h treatment groups and corresponding controls were collected in the same way. The experiment was replicated three times independently. Each sample was then homogenized manually on ice with 500 µL phosphate buffer solution (PBS, 0.1 M, pH = 7.5) using a plastic grinding pestle (TIANGEN, Beijing, China). Each homogenate was centrifuged by a 5471R centrifuge (Eppendorf, Hamburg, Germany) at 5900 g for 5 min at 4 °C, and then the supernatant was transferred into a new Eppendorf tube and centrifuged at 13,300× g for 30 min at 4 °C. The final supernatant was used to test PO activity. According to a previous study [56], bovine serum albumin (BSA) was set as a standard for determining the protein content of the supernatant, and the absorbance was detected at 595 nm in a Model 680 microplate reader (Bio-Rad, Hercules, CA, USA). PO activity was measured according to a previous study [57]. This mixture contained 30 µL of the supernatant, 135 µL PBS, and 135 µL of 0.0375 moL catechol. The absorbance value was detected at 420 nm for 2 min with a Model 680 microplate reader. A mixture that contained only catechol (CK) and PBS was used as a control and measured. The experiment was independently repeated three times. The results are indicated as the absorbance value of enzyme activity per mg protein (ΔOD420 nm·min−1·mg protein−1).

4.3. Cloning of the Full-Length Prophenoloxidase Gene of O. formosanus

Total RNA from non-treated O. formosanus workers was extracted using TsingZol Total RNA Extraction Reagent (TSP401, Tsingke Biotechnology Beijing, Co., Ltd., Beijing, China) according to the manufacturer’s protocol. The integrity and quantity of RNA were checked using 1% agarose gel electrophoresis and a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). O. formosanus cDNA was synthesized using an Evo M-MLV Plus 1st Strand cDNA Synthesis Kit (AG11615, Accurate Biotechnology, Hunan, Co., Ltd., Changsha, China) according to the manufacturer’s protocol, and cDNA was used as template to clone OfPPO gene. The PCR primers specific to OfPPO were designed using Premier 5.0 based on the full-length sequence of O. formosanus PPO obtained from the full-length transcriptome data (Table 2). PCR was executed using PrimeSTAR Max DNA Polymerase (Takara, Dalian, Liaoning, China). The PCR program was set as follows: 98 °C for 3 min; 35 cycles of 98 °C for 10 s, 52 °C for 15 s and 72 °C for 60 s; and an extension cycle of 72 °C for 5 min. The OfPPO DNA was subcloned into the TA clone vector and sequenced by Sangon Biotech (Shanghai) Co., Ltd., Shanghai, China.

4.4. Characteristics of OfPPO and Phylogenetic Analysis

The amino acid sequences were deduced from the NCBI Open Reading Frame (ORF) finder (https://www.ncbi.nlm.nih.gov/orffinder/, accessed on 20 June 2022). Then, their molecular weight and isoelectric point were calculated by the Expert Protein Analysis System (EXPASY) proteomics server (http://www.expasy.org/, accessed on 20 June 2022). The amino acid sequences of OfPPO were aligned with several known characterized PPOs from C. formosanus, Z. nevadensis, B. germanica, B. mori, P. xylostella, S. frugiperda, B. dorsalis, T. castaneum and T. madens to analyze conserved OfPPO motifs using DNAMAN multiple sequence alignment analysis (Version 6, Lynnon Biosoft Co., San Ramon, CA, USA). To construct the phylogenetic tree and explore the relationship between OfPPO and other PPOs, we chose prophenoloxidase as a keyword to search the nonredundant database from the NCBI website (https://www.ncbi.nlm.nih.gov/, accessed on 20 June 2022). Multiple amino acid sequence alignment analysis was performed using MEGA X (version 10.1) and Clustal X software (version 2.1). The phylogenetic tree was constructed using the neighbor-joining (NJ) method in MEGA X with 1000 bootstrap replicates.

4.5. OfPPO Expression in O. formosanus Using qPCR

Total RNA was extracted from each sample using TsingZol Total RNA Extraction Reagent, and mRNA was first reverse transcribed to cDNA using PrimeScript™ RT Master Mix (Takara, Dalian, Liaoning, China) according to the manufacturer’s protocol and subsequently verified by qPCR using 2×TSINGKE Master qPCR Mix (TSE201). In general, qPCR was carried out in a 20 µL volume consisting of 1 µL of cDNA (100 ng/uL) and 10 µL of 2×TSINGKE Master qPCR Mix, 0.4 µL of ROX Reference Dye II, 0.8 µL each of 10 µM forward and reverse primer and 7 µL of ddH2O using ViATM 7 Real-time PCR system (Applied Biosystems, Foster City, CA, USA). The qPCR program was 95 °C for 60 s, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s, after which the melting curve was analyzed. Each primer for qPCR was designed using Premier 5.0 software according to the gene sequence. GAPDH and RPS18 were used as internal reference genes. All primers are listed in Table 2. The amplification efficiency of primers was tested using Linreg PCR software (Version: 2016.1), according to the method described by Ramakers [58]. All primer sets utilized in this study had mean amplification efficiency values ranged between 1.7 and 2.0. Caste- and tissue-specific expression experiments were calculated using the 2−∆Ct method, and other experiments were calculated using the 2−∆∆Ct method [59,60]. Three replicates were performed for each treatment independently.

4.6. In Vitro Synthesis of dsOfPPOs (dsOfPPO1 and dsOfPPO2)

Two dsOfPPO sequences were designed based on the copper-binding site B motif and C-terminal conserved motif to evaluate the RNAi efficiency. Two dsOfPPOs were amplified from the OfPPO subcloned vector using dsOfPPO1 and dsOfPPO2 primers, and the restriction enzyme cutting sites were added to each primer. The PCR products of dsOfPPOs were digested by restriction enzyme, purified using a DNA Gel Extraction Kit (TSP601-50, Tsingke Biotechnology Co., Ltd., Beijing, China), and cloned into vector L4440 to produce dsOfPPOs.
The recombinant plasmids of dsOfPPO1 and dsOfPPO2 were transformed into E. coli HT115 competent cells [61,62]. The transformed bacteria were cultured in LB medium containing 100 μg/mL ampicillin and 12.5 μg/mL tetracycline at 37 ℃ for 12 h, with continuous shaking (200 rpm). A 50 μL aliquot of cultured bacteria was added to 50 mL of fresh LB medium containing 100 μg/mL ampicillin and 12.5 μg/mL tetracycline and allowed to grow until OD600 = 0.6–0.7. In order to induce the synthesis of dsRNA, isopropyl-beta-D-thiogalactopyranoside (IPTG, final concentration of 0.8 mM) was added to the bacteria culture and the mixture was cultured continuously at 37 °C for 12 h. Fifty microliter of the resulting bacterial culture was utilized to isolate dsRNA using TsingZol Total RNA Extraction kit. The integrity and quantity of RNA were analyzed using 1% agarose gel electrophoresis and a NanoDrop spectrophotometer. The dsRNA isolates were used for subsequent RNAi experiments. Double stranded RNA of GFP was used as a control dsRNA treatment and an expression plasmid containing the GFP gene was utilized to produce dsGFP using the same procedures that generated dsOfPPOs.

4.7. RNAi Efficiency Evaluation after Treatment with dsRNA

For triplicate dsRNA treatment experiments, twenty O. formosanus workers in a culture dish covered with wet filter paper were fed with 0.4 mL solution containing 1 μg/μL dsOfPPO and 10% (w/v) Nile blue (to recognize whether O. formosanus workers fed the fluid). The positive control group was fed with 0.4 mL dsGFP and Nile blue, and the negative control group was fed with 0.4 mL sterilized water containing Nile blue. Each treatment was repeated three times. At 3, 6, 12, and 24 h post treatment, the workers whose intestines turned blue by visual observation were collected and subjected for qPCR evaluation for RNAi interference efficiency on OfPPO gene expression.

4.8. The Role of OfPPOs in Resisting SM1 and Bt Infection

A number of O. formosanus workers for this experiment were placed into culture dishes with 0.4 mL of sterilized water and Nile blue (negative control), dsGFP and Nile blue (positive control), dsOfPPOs (dsOfPPO1 or dsOfPPO2) and Nile blue (treatment group) for 6 h. Then, workers from the negative control, positive control and treatment groups were divided into two groups (uninfected and bacterial infection groups), respectively. One microliter of culture medium was placed on each worker’s pronotum (control group), and 1 μL of SM1 (8.75 × 1010 CFU/mL) or 1 μL of Bt (2.14 × 1010 CFU/mL) was placed on each worker’s pronotum (treatment group). Each treatment contained 20 workers and was run in triplicates. The experiments were monitored hourly until one of the following conditions was met: (1) all members in an infected group died, and the mortality rate in the uninfected group was less than 20%. (2) The mortality rate in the uninfected group was at least 20% regardless of the mortality rate in the infected group.

4.9. Statistical Analysis

InStat software version 3.05 (GraphPad, San Diego, CA, USA) was used to analyze the data. The significant differences between the two samples in the PO activity and OfPPO expression assays were analyzed using a t test. One-way analysis of variance and Tukey’s multiple comparisons were used to analyze caste- and tissue-specific gene expression and RNA interference assays. A p value < 0.05 indicated a significant difference. For the survival assay, GraphPad Prism version 8.0.2 (GraphPad Software, San Diego, CA, USA) was used to analyze the data according to the log-rank (Mantel–Cox) test.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24010406/s1.

Author Contributions

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

Funding

This research was supported by The Natural Science Foundation of the Jiangsu Higher Education Institutions of China [grant Number 20KJA220003], Postgraduate Research & Practice Innovation Program of Jiangsu Province [grant numbers KYCX21_0924] and the Priority Academic Program Development Fund of Jiangsu Higher Education Institutions.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and supplementary material.

Acknowledgments

The authors thank Yinjuan Zhao (Nanjing Forestry University) for providing Bacillus thuringiensis strain and Jing Wei (Soochow University) for editing and revising the English language. The authors also thank the two anonymous reviewers for their helpful comments and suggestions for improving the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rust, M.K.; Su, N.Y. Managing social insects of urban importance. Annu. Rev. Entomol. 2012, 57, 355–375. [Google Scholar] [CrossRef] [PubMed]
  2. Ahmad, F.; Fouad, H.; Liang, S.Y.; Hu, Y.; Mo, J.C. Termites and Chinese agricultural system: Applications and advances in integrated termite management and chemical control. Insect. Sci. 2021, 28, 2–20. [Google Scholar] [CrossRef] [PubMed]
  3. Husseneder, C.; Garner, S.P.; Huang, Q.; Booth, W.; Vargo, E.L. Characterization of microsatellites for population genetic analyses of the fungus-growing termite Odontotermes formosanus (Isoptera: Termitidae). Environ. Entomol. 2013, 42, 1092–1099. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Appel, A.G.; Hu, X.P.; Zhou, J.; Qin, Z.; Zhu, H.; Chang, X.; Wang, Z.; Liu, X.; Liu, M. Observations of the biology and ecology of the black-winged termite, Odontotermes formosanus Shiraki (Termitidae: Isoptera), in camphor, Cinnamomum camphora (L.) (Lauraceae). Psyche 2012, 2012, 123102. [Google Scholar]
  5. Huang, Q.Y.; Lei, C.L.; Xue, D. Field evaluation of a fipronil bait against subterranean termite Odontotermes formosanus (Isoptera: Termitidae). J. Econ. Entomol. 2006, 99, 455–461. [Google Scholar] [CrossRef]
  6. Huang, Q.; Sun, P.; Zhou, X.; Lei, C. Characterization of head transcriptome and analysis of gene expression involved in caste differentiation and aggression in Odontotermes formosanus (Shiraki). PLoS ONE 2012, 7, e50383. [Google Scholar] [CrossRef] [Green Version]
  7. Tudi, M.; Daniel, R.H.; Wang, L.; Lyu, J.; Sadler, R.; Connell, D.; Chu, C.; Phung, D.T. Agriculture Development, Pesticide Application and Its Impact on the Environment. Int. J. Environ. Res. Public Health 2021, 18, 1112. [Google Scholar] [CrossRef]
  8. Aggarwal, C.; Paul, S.; Nain, V.; Tripathi, V.; Paul, B.; Aslam, K.M. Comparative response of Spodoptera litura challenged per os with Serratia marcescens strains differing in virulence. J. Invertebr. Pathol. 2021, 183, 107562. [Google Scholar] [CrossRef]
  9. Wang, Z.; Feng, K.; Tang, F.; Xu, M. Activation of the host immune response in Hyphantria cunea (Drury) (Lepidoptera: Noctuidae) induced by Serratia marcescens Bizio. Insects 2021, 12, 983. [Google Scholar] [CrossRef]
  10. Niu, H.; Wang, N.; Liu, B.; Xiao, L.; Wang, L.; Guo, H. Synergistic and additive interactions of Serratia marcescens S-JS1 to the chemical insecticides for controlling Nilaparvata lugens (Hemiptera: Delphacidae). J. Econ. Entomol. 2018, 111, 823–828. [Google Scholar] [CrossRef]
  11. Wang, C.; Henderson, G.; Gautam, B.K. Lufenuron suppresses the resistance of Formosan subterranean termites (Isoptera: Rhinotermitidae) to entomopathogenic bacteria. J. Econ. Entomol. 2013, 106, 1812–1818. [Google Scholar] [CrossRef] [PubMed]
  12. Wright, M.S.; Cornelius, M.L. Mortality and repellent effects of microbial pathogens on Coptotermes formosanus (Isoptera: Rhinotermitidae). BMC Microbiol. 2012, 12, 291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Fu, R.; Luo, J.; Feng, K.; Lu, X.; Tang, F. Termite-killing components in Serratia marcescens (SM1). J. For. Res. 2021, 32, 1739–1744. [Google Scholar] [CrossRef]
  14. Hoffmann, J.A. Innate immunity of insects. Curr. Opin. Immunol. 1995, 7, 4–10. [Google Scholar] [CrossRef]
  15. Hoffmann, J.A. The immune response of Drosophila. Nature 2003, 426, 33–38. [Google Scholar] [CrossRef]
  16. Ferrandon, D.; Imler, J.L.; Hetru, C.; Hoffmann, J.A. The Drosophila systemic immune response: Sensing and signalling during bacterial and fungal infections. Nat. Rev. Immunol. 2007, 7, 862–874. [Google Scholar] [CrossRef]
  17. Hillyer, J.F. Insect immunology and hematopoiesis. Dev. Comp. Immunol. 2016, 58, 102–118. [Google Scholar] [CrossRef] [Green Version]
  18. Caccia, S.; Astarita, F.; Barra, E.; Di Lelio, I.; Varricchio, P.; Pennacchio, F. Enhancement of Bacillus thuringiensis toxicity by feeding Spodoptera littoralis larvae with bacteria expressing immune suppressive dsRNA. J. Pest Sci. 2020, 93, 303–314. [Google Scholar] [CrossRef] [Green Version]
  19. Fan, Y.; Abbas, M.; Liu, X.; Wang, Y.; Song, H.; Li, T.; Ma, E.; Zhu, K.Y.; Zhang, J. Increased RNAi efficiency by dsEGFP-induced up-regulation of two core RNAi pathway genes (OfDicer2 and OfAgo2) in the Asian Corn Borer (Ostrinia furnacalis). Insects 2022, 13, 274. [Google Scholar] [CrossRef]
  20. Moreira-Pinto, C.E.; Coelho, R.R.; Leite, A.; Silveira, D.A.; de Souza, D.A.; Lopes, R.B.; Macedo, L.; Silva, M.; Ribeiro, T.P.; Morgante, C.V.; et al. Increasing Anthonomus grandis susceptibility to Metarhizium anisopliae through RNAi-induced AgraRelish knockdown: A perspective to combine biocontrol and biotechnology. Pest. Manag. Sci. 2021, 77, 4054–4063. [Google Scholar] [CrossRef]
  21. Maximo, W.; Howell, J.L.; Mogilicherla, K.; Basij, M.; Chereddy, S.; Palli, S.R. Inhibitor of apoptosis is an effective target gene for RNAi-mediated control of Colorado potato beetle, Leptinotarsa decemlineata. Arch. Insect. Biochem. Physiol. 2020, 104, e21685. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, W.; Keyhani, N.O.; Zhang, H.; Cai, K.; Xia, Y. Inhibitor of apoptosis-1 gene as a potential target for pest control and its involvement in immune regulation during fungal infection. Pest. Manag. Sci. 2020, 76, 1831–1840. [Google Scholar] [CrossRef] [PubMed]
  23. Bingsohn, L.; Knorr, E.; Billion, A.; Narva, K.E.; Vilcinskas, A. Knockdown of genes in the Toll pathway reveals new lethal RNA interference targets for insect pest control. Insect. Mol. Biol. 2017, 26, 92–102. [Google Scholar] [CrossRef]
  24. Zhao, P.; Li, J.; Wang, Y.; Jiang, H. Broad-spectrum antimicrobial activity of the reactive compounds generated in vitro by Manduca sexta phenoloxidase. Insect Biochem. Mol. Biol. 2007, 37, 952–959. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kanost, M.R.; Jiang, H.; Yu, X.Q. Innate immune responses of a lepidopteran insect, Manduca sexta. Immunol. Rev. 2004, 198, 97–105. [Google Scholar] [CrossRef] [PubMed]
  26. Lavine, M.D.; Strand, M.R. Surface characteristics of foreign targets that elicit an encapsulation response by the moth Pseudoplusia includens. J. Insect. Physiol. 2001, 47, 965–974. [Google Scholar] [CrossRef]
  27. Cerenius, L.; Lee, B.L.; Soderhall, K. The proPO-system: Pros and cons for its role in invertebrate immunity. Trends Immunol. 2008, 29, 263–271. [Google Scholar] [CrossRef]
  28. Prabu, S.; Jing, D.; Shabbir, M.Z.; Yuan, W.; Wang, Z.; He, K. Contribution of phenoloxidase activation mechanism to Bt insecticidal protein resistance in Asian corn borer. Int. J. Biol. Macromol. 2020, 153, 88–99. [Google Scholar] [CrossRef]
  29. Wang, Y.; Jiang, H.; Cheng, Y.; An, C.; Chu, Y.; Raikhel, A.S.; Zou, Z. Activation of Aedes aegypti prophenoloxidase-3 and its role in the immune response against entomopathogenic fungi. Insect Mol. Biol. 2017, 26, 552–563. [Google Scholar] [CrossRef]
  30. Xu, L.; Ma, L.; Wang, W.; Li, L.; Lu, Z. Phenoloxidases are required for the pea aphid’s defence against bacterial and fungal infection. Insect Mol. Biol. 2019, 28, 176–186. [Google Scholar] [CrossRef]
  31. Zhang, J.; Huang, W.; Yuan, C.; Lu, Y.; Yang, B.; Wang, C.Y.; Zhang, P.; Dobens, L.; Zou, Z.; Wang, C.; et al. Prophenoloxidase-mediated ex vivo immunity to delay fungal infection after insect ecdysis. Front. Immunol. 2017, 8, 1445. [Google Scholar] [CrossRef] [PubMed]
  32. Lacey, L.A.; Grzywacz, D.; Shapiro-Ilan, D.I.; Frutos, R.; Brownbridge, M.; Goettel, M.S. Insect pathogens as biological control agents: Back to the future. J. Invertebr. Pathol. 2015, 132, 1–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chouvenc, T.; Su, N.; Grace, J.K. Fifty years of attempted biological control of termites—Analysis of a failure. Biol. Control 2011, 59, 69–82. [Google Scholar] [CrossRef]
  34. Rodriguez-Andres, J.; Rani, S.; Varjak, M.; Chase-Topping, M.E.; Beck, M.H.; Ferguson, M.C.; Schnettler, E.; Fragkoudis, R.; Barry, G.; Merits, A.; et al. Phenoloxidase activity acts as a mosquito innate immune response against infection with Semliki Forest virus. PLoS Pathog. 2012, 8, e1002977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Yokoi, K.; Hayakawa, Y.; Kato, D.; Minakuchi, C.; Tanaka, T.; Ochiai, M.; Kamiya, K.; Miura, K. Prophenoloxidase genes and antimicrobial host defense of the model beetle, Tribolium castaneum. J. Invertebr. Pathol. 2015, 132, 190–200. [Google Scholar] [CrossRef] [PubMed]
  36. Yassine, H.; Kamareddine, L.; Osta, M.A. The mosquito melanization response is implicated in defense against the entomopathogenic fungus Beauveria bassiana. PLoS Pathog. 2012, 8, e1003029. [Google Scholar] [CrossRef]
  37. Liu, L.; Li, G.; Sun, P.; Lei, C.; Huang, Q. Experimental verification and molecular basis of active immunization against fungal pathogens in termites. Sci. Rep. 2015, 5, 15106. [Google Scholar] [CrossRef] [Green Version]
  38. Kuraishi, T.; Hori, A.; Kurata, S. Host-microbe interactions in the gut of Drosophila melanogaster. Front. Physiol. 2013, 4, 375. [Google Scholar] [CrossRef] [Green Version]
  39. Bidla, G.; Hauling, T.; Dushay, M.S.; Theopold, U. Activation of insect phenoloxidase after injury: Endogenous versus foreign elicitors. J. Innate. Immun. 2009, 1, 301–308. [Google Scholar] [CrossRef]
  40. Karlsson, C.; Korayem, A.M.; Scherfer, C.; Loseva, O.; Dushay, M.S.; Theopold, U. Proteomic analysis of the Drosophila larval hemolymph clot. J. Biol. Chem. 2004, 279, 52033–52041. [Google Scholar] [CrossRef] [Green Version]
  41. Zou, Z.; Shin, S.W.; Alvarez, K.S.; Kokoza, V.; Raikhel, A.S. Distinct melanization pathways in the mosquito Aedes aegypti. Immunity 2010, 32, 41–53. [Google Scholar] [CrossRef] [PubMed]
  42. Zou, Z.; Shin, S.W.; Alvarez, K.S.; Bian, G.; Kokoza, V.; Raikhel, A.S. Mosquito RUNX4 in the immune regulation of PPO gene expression and its effect on avian malaria parasite infection. Proc. Natl. Acad. Sci. USA 2008, 105, 18454–18459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Li, Z.Q.; Song, X.H.; Wang, M.; Wang, S.; Huang, G.H. Melanization induced by Heliothis virescens ascovirus 3h promotes viral replication. Insect Sci. 2021, 28, 472–484. [Google Scholar] [CrossRef] [PubMed]
  44. Tartar, A.; Wheeler, M.M.; Zhou, X.; Coy, M.R.; Boucias, D.G.; Scharf, M.E. Parallel metatranscriptome analyses of host and symbiont gene expression in the gut of the termite Reticulitermes flavipes. Biotechnol. Biofuels 2009, 2, 25. [Google Scholar] [CrossRef] [Green Version]
  45. Rosengaus, R.B.; Reichheld, J.L. Phenoloxidase activity in the infraorder Isoptera: Unraveling life-history correlates of immune investment. Naturwissenschaften 2016, 103, 14. [Google Scholar] [CrossRef]
  46. Yang, N.; Ding, T.; Chu, D. Silencing of the prophenoloxidase gene BtPPO1 increased the ability of acquisition and retention of Tomato chlorosis virus by Bemisia tabaci. Int. J. Mol. Sci. 2022, 23, 6541. [Google Scholar] [CrossRef]
  47. Lourenco, A.P.; Zufelato, M.S.; Bitondi, M.M.; Simoes, Z.L. Molecular characterization of a cDNA encoding prophenoloxidase and its expression in Apis mellifera. Insect Biochem. Mol. Biol. 2005, 35, 541–552. [Google Scholar] [CrossRef]
  48. Park, D.S.; Shin, S.W.; Kim, M.G.; Park, S.S.; Lee, W.J.; Brey, P.T.; Park, H.Y. Isolation and characterization of the cDNA encoding the prophenoloxidase of fall webworm, Hyphantria cunea. Insect Biochem. Mol. Biol. 1997, 27, 983–992. [Google Scholar] [CrossRef]
  49. Muller, H.M.; Dimopoulos, G.; Blass, C.; Kafatos, F.C. A hemocyte-like cell line established from the malaria vector Anopheles gambiae expresses six prophenoloxidase genes. J. Biol. Chem. 1999, 274, 11727–11735. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Zheng, H.; Li, L.; Xu, Q.; Zou, Q.; Tang, B.; Wang, S. Gene cloning and expression patterns of two prophenoloxidases from Catantops pinguis (Orthoptera: Catantopidae). Bull. Entomol. Res. 2013, 103, 393–405. [Google Scholar] [CrossRef] [PubMed]
  51. Asano, T.; Ashida, M. Cuticular pro-phenoloxidase of the silkworm, Bombyx mori. Purification and demonstration of its transport from hemolymph. J. Biol. Chem. 2001, 276, 11100–11112. [Google Scholar] [CrossRef] [PubMed]
  52. Khajuria, C.; Velez, A.M.; Rangasamy, M.; Wang, H.; Fishilevich, E.; Frey, M.L.; Carneiro, N.P.; Gandra, P.; Narva, K.E.; Siegfried, B.D. Parental RNA interference of genes involved in embryonic development of the western corn rootworm, Diabrotica virgifera virgifera LeConte. Insect Biochem. Mol. Biol. 2015, 63, 54–62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Terenius, O.; Papanicolaou, A.; Garbutt, J.S.; Eleftherianos, I.; Huvenne, H.; Kanginakudru, S.; Albrechtsen, M.; An, C.; Aymeric, J.L.; Barthel, A.; et al. RNA interference in Lepidoptera: An overview of successful and unsuccessful studies and implications for experimental design. J. Insect Physiol. 2011, 57, 231–245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Caccia, S.; Di Lelio, I.; La Storia, A.; Marinelli, A.; Varricchio, P.; Franzetti, E.; Banyuls, N.; Tettamanti, G.; Casartelli, M.; Giordana, B.; et al. Midgut microbiota and host immunocompetence underlie Bacillus thuringiensis killing mechanism. Proc. Natl. Acad. Sci. USA 2016, 113, 9486–9491. [Google Scholar] [CrossRef] [Green Version]
  55. Joga, M.R.; Zotti, M.J.; Smagghe, G.; Christiaens, O. RNAi Efficiency, systemic properties, and novel delivery methods for pest insect control: What we know so far. Front. Physiol. 2016, 7, 553. [Google Scholar] [CrossRef] [Green Version]
  56. Bradford, M.M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248–254. [Google Scholar] [CrossRef]
  57. BENJAMIN, N.D.; Montgomery, M.W. Polyphenoloxidase of Royal Ann cherries: Purification and characterization. J. Food Sci. 1973, 38, 799–806. [Google Scholar] [CrossRef]
  58. Ramakers, C.; Ruijter, J.M.; Deprez, R.H.; Moorman, A.F. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci. Lett. 2003, 339, 62–66. [Google Scholar] [CrossRef]
  59. Muller, P.Y.; Janovjak, H.; Miserez, A.R.; Dobbie, Z. Processing of gene expression data generated by quantitative real-time RT-PCR. Biotechniques 2002, 32, 1372–1374, 1376, 1378–1379. [Google Scholar]
  60. Pfaffl, M.W.; Horgan, G.W.; Dempfle, L. Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30, e36. [Google Scholar] [CrossRef]
  61. Newmark, P.A.; Reddien, P.W.; Cebria, F.; Sanchez, A.A. Ingestion of bacterially expressed double-stranded RNA inhibits gene expression in planarians. Proc. Natl. Acad. Sci. USA 2003, 100 (Suppl. 1), 11861–11865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Timmons, L.; Court, D.L.; Fire, A. Ingestion of bacterially expressed dsRNAs can produce specific and potent genetic interference in Caenorhabditis elegans. Gene 2001, 263, 103–112. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The melanization and enzymatic activity of PO in O. formosanus treated by SM1 and Bt. (A): The level of melanization in the hemolymph of O. formosanus infected with bacteria treated for 0, 1, 2, 3 and 4 h. The treatment groups were treated with SM1 or Bt. The control group was treated with culture medium (CK). The melanization reaction was recorded by photography at different time points from 0 to 4 h. (B): PO activity in O. formosanus infected by SM1. The control group was treated with culture medium (CK). (C): PO activity in O. formosanus infected by Bt. The control group was treated with culture medium (CK). The data are presented as the mean ± SD of three replicates. The asterisks indicate significant differences between the treatments groups and the control groups (* p < 0.05, ** p < 0.01).
Figure 1. The melanization and enzymatic activity of PO in O. formosanus treated by SM1 and Bt. (A): The level of melanization in the hemolymph of O. formosanus infected with bacteria treated for 0, 1, 2, 3 and 4 h. The treatment groups were treated with SM1 or Bt. The control group was treated with culture medium (CK). The melanization reaction was recorded by photography at different time points from 0 to 4 h. (B): PO activity in O. formosanus infected by SM1. The control group was treated with culture medium (CK). (C): PO activity in O. formosanus infected by Bt. The control group was treated with culture medium (CK). The data are presented as the mean ± SD of three replicates. The asterisks indicate significant differences between the treatments groups and the control groups (* p < 0.05, ** p < 0.01).
Ijms 24 00406 g001
Figure 2. Multiple amino acid sequence alignment of insect PPOs. Protein sequences of from Odontotermes formosanus PPO (OfPPO), Coptotermes formosanus PO2 (CfPO2, AHB39936.1), Zootermopsis nevadensis PO2 (ZnPO2, XP_021921811.1), Blattella germanica PO1 (BgPO1, PSN38098.1), Bombyx mori PPO1 (BmPPO1, AAG09304), Plutella xyostella PPO2 (PxPPO2, ACS36209.1), Spodoptera frugiperda PPO1 (SfPPO1, ABB92834), Bactrocera dorsalis PPO1 (BdPPO1, AFP81887), Tribolium castaneum PPO2 (TcPPO2, NP_001034522), and Tenebrio madens PO1 (TmPO1, XP_044260894.1) were aligned using DNAMAN. The two potential proteolytic cleavage sites were shown by red and green lines, and the predicted cleavage bonds were marked with black arrows. CuA and CuB were underlined with blue and black lines respectively, and the conserved six histidine residues are marked with black dots. The thiol ester sites and C-terminal conserved motif were underlined with pink and gray lines, respectively.
Figure 2. Multiple amino acid sequence alignment of insect PPOs. Protein sequences of from Odontotermes formosanus PPO (OfPPO), Coptotermes formosanus PO2 (CfPO2, AHB39936.1), Zootermopsis nevadensis PO2 (ZnPO2, XP_021921811.1), Blattella germanica PO1 (BgPO1, PSN38098.1), Bombyx mori PPO1 (BmPPO1, AAG09304), Plutella xyostella PPO2 (PxPPO2, ACS36209.1), Spodoptera frugiperda PPO1 (SfPPO1, ABB92834), Bactrocera dorsalis PPO1 (BdPPO1, AFP81887), Tribolium castaneum PPO2 (TcPPO2, NP_001034522), and Tenebrio madens PO1 (TmPO1, XP_044260894.1) were aligned using DNAMAN. The two potential proteolytic cleavage sites were shown by red and green lines, and the predicted cleavage bonds were marked with black arrows. CuA and CuB were underlined with blue and black lines respectively, and the conserved six histidine residues are marked with black dots. The thiol ester sites and C-terminal conserved motif were underlined with pink and gray lines, respectively.
Ijms 24 00406 g002
Figure 3. Phylogenetic analysis of the protein sequences of 82 PPOs, including 38 insects and 16 crustaceans. The PPO gene name was shown as the Latin name of the species, and the NCBI gene accession number was added. The O. formosanus PPO was marked in red. The branches specific for Lepidoptera, Blattodea, Coleoptera, Hymenoptera, Crustacea, and Diptera PPOs were shaded in different colors. Branch colors represent bootstrap values.
Figure 3. Phylogenetic analysis of the protein sequences of 82 PPOs, including 38 insects and 16 crustaceans. The PPO gene name was shown as the Latin name of the species, and the NCBI gene accession number was added. The O. formosanus PPO was marked in red. The branches specific for Lepidoptera, Blattodea, Coleoptera, Hymenoptera, Crustacea, and Diptera PPOs were shaded in different colors. Branch colors represent bootstrap values.
Ijms 24 00406 g003
Figure 4. Expression profiles of OfPPO in different castes and tissues of O. formosanus. (A) OfPPO mRNA expression levels in five castes. (B) OfPPO mRNA expression levels in five tissues of female dealates. (C) OfPPO mRNA expression levels in five tissues of male dealates. (D) OfPPO mRNA expression levels in five tissues of worker. (E) OfPPO mRNA expression levels in five tissues of larval instars. (F) OfPPO mRNA expression levels in five tissues of soldiers. The data are presented as the mean ± SD of three replicate samples. The same letters on the error bars indicate no significant difference in OfPPO expression (p > 0.05).
Figure 4. Expression profiles of OfPPO in different castes and tissues of O. formosanus. (A) OfPPO mRNA expression levels in five castes. (B) OfPPO mRNA expression levels in five tissues of female dealates. (C) OfPPO mRNA expression levels in five tissues of male dealates. (D) OfPPO mRNA expression levels in five tissues of worker. (E) OfPPO mRNA expression levels in five tissues of larval instars. (F) OfPPO mRNA expression levels in five tissues of soldiers. The data are presented as the mean ± SD of three replicate samples. The same letters on the error bars indicate no significant difference in OfPPO expression (p > 0.05).
Ijms 24 00406 g004
Figure 5. Elevated OfPPO expression in O. formosanus workers infected by SM1 (A) and Bt (B). The data are presented as the mean ± SD of three replicate samples. The asterisks indicate significant differences between the treatment groups and the control groups (** p < 0.01).
Figure 5. Elevated OfPPO expression in O. formosanus workers infected by SM1 (A) and Bt (B). The data are presented as the mean ± SD of three replicate samples. The asterisks indicate significant differences between the treatment groups and the control groups (** p < 0.01).
Ijms 24 00406 g005
Figure 6. Analysis of the silencing efficiency of OfPPO in O. formosanus. (A) Conceptual diagrams illustrating the interference regions for designing dsRNA fragments (including dsOfPPO1 and dsOfPPO2). Black arrows indicate the regions for detecting the expression level of OfPPO by qPCR. (B) Relative expression levels of OfPPO in O. formosanus fed on dsOfPPO1. (C) Relative expression levels of OfPPO in O. formosanus fed on dsOfPPO2. The data are presented as the mean ± SD of three replicates. The same letters indicate no significant difference in the OfPPO expression (p > 0.05).
Figure 6. Analysis of the silencing efficiency of OfPPO in O. formosanus. (A) Conceptual diagrams illustrating the interference regions for designing dsRNA fragments (including dsOfPPO1 and dsOfPPO2). Black arrows indicate the regions for detecting the expression level of OfPPO by qPCR. (B) Relative expression levels of OfPPO in O. formosanus fed on dsOfPPO1. (C) Relative expression levels of OfPPO in O. formosanus fed on dsOfPPO2. The data are presented as the mean ± SD of three replicates. The same letters indicate no significant difference in the OfPPO expression (p > 0.05).
Ijms 24 00406 g006
Figure 7. Bioassay test of SM1-challenged O. formosanus treated with dsOfPPO1. (A) The mortality rate of O. formosanus in each group at 6 h. (B) The mortality rate of O. formosanus in each group at 12 h. (C) The mortality rate of O. formosanus in each group at 24 h. (D) The mortality rate of O. formosanus in each group at 48 h. (E) The survival rate of SM1-challenged O. formosanus treated with dsOfPPO1. The data are presented as the mean ± SD of three replicates. In these figures, the same letters indicate no significant differences in the mortality of O. formosanus (p > 0.05).
Figure 7. Bioassay test of SM1-challenged O. formosanus treated with dsOfPPO1. (A) The mortality rate of O. formosanus in each group at 6 h. (B) The mortality rate of O. formosanus in each group at 12 h. (C) The mortality rate of O. formosanus in each group at 24 h. (D) The mortality rate of O. formosanus in each group at 48 h. (E) The survival rate of SM1-challenged O. formosanus treated with dsOfPPO1. The data are presented as the mean ± SD of three replicates. In these figures, the same letters indicate no significant differences in the mortality of O. formosanus (p > 0.05).
Ijms 24 00406 g007
Figure 8. Bioassay test of Bt-challenged O. formosanus treated with dsOfPPO1. (A) The mortality rate of O. formosanus in each group at 24 h. (B) The mortality rate of O. formosanus in each group at 48 h. (C) The mortality rate of O. formosanus in each group at 72 h. (D) The mortality rate of O. formosanus in each group at 96 h. (E) The survival rate of Bt-challenged O. formosanus treated with dsOfPPO1. The data are presented as the mean ± SD of three replicates. In these figures, the same letters indicate no significant differences in the mortality of O. formosanus (p > 0.05).
Figure 8. Bioassay test of Bt-challenged O. formosanus treated with dsOfPPO1. (A) The mortality rate of O. formosanus in each group at 24 h. (B) The mortality rate of O. formosanus in each group at 48 h. (C) The mortality rate of O. formosanus in each group at 72 h. (D) The mortality rate of O. formosanus in each group at 96 h. (E) The survival rate of Bt-challenged O. formosanus treated with dsOfPPO1. The data are presented as the mean ± SD of three replicates. In these figures, the same letters indicate no significant differences in the mortality of O. formosanus (p > 0.05).
Ijms 24 00406 g008
Table 1. Percentage of similarity among six PPO amino acid sequences calculated using MEGA X.
Table 1. Percentage of similarity among six PPO amino acid sequences calculated using MEGA X.
NameOfPPOCfPO2ZnPO2BgPO1CsPO1TcPPO2
OfPPO100.092.8087.9079.3460.8463.40
CfPO2 100.089.0580.5660.1262.96
ZnPO2 100.080.5659.6861.93
BgPO1 100.059.8263.70
CsPO1 10076.83
TcPPO2 100
Of: Odontotermes formosanus; Cf: Coptotermes formosanus; Zn: Zootermopsis nevadensis; Bg: Blattella germanica; Cs: Coccinella septempunctata; Tc: Tribolium castaneum.
Table 2. Primers used for this study.
Table 2. Primers used for this study.
Primer NamesForward Primer Sequences (5′-3′)Reverse Primer Sequences (5′-3′)Annealing TemperatureExperiments
OfPPOATGTCAAGACCAGCACAATCTCATGAGAGCCGGTTTG52 °CFull-length
amplification
dsOfPPO1CGAGCTCGGTACAACTTCGAACGGCTGTCCCCCGGGGCTTGCCTCCATCAGGTTAC57 °CRNAi
dsOfPPO2CGAGCTCTCCAACAACGATGATGATAAGATATCCCCCGGGGTTTGGCTCATCTCTACCCTG57 °C
dsGFPATCGGAGCTCTAGTTGAACGGATCCATCTTCACCCAAGCTTAGAACTTTTCACTGGA54 °C
OfPPOGTTGATCTGTCCCGAGGTATTGGCTGTTCTCAGCCTGAATCTTA60 °CqPCR
GAPDHTCGTATTGGCCGTCTTGTGCAGCGACCATGGGTGGAATCAT60 °C
RPS18ATGGCAAACCCCCGTCAGTACATACCACGATGCGCACGAA60 °C
Note: The underline within primers indicated the restriction sites; GAPDH and RPS18 were the internal reference genes in this study.
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

Wang, Z.; Luo, J.; Feng, K.; Zhou, Y.; Tang, F. Prophenoloxidase of Odontotermes formosanus (Shiraki) (Blattodea: Termitidae) Is a Key Gene in Melanization and Has a Defensive Role during Bacterial Infection. Int. J. Mol. Sci. 2023, 24, 406. https://doi.org/10.3390/ijms24010406

AMA Style

Wang Z, Luo J, Feng K, Zhou Y, Tang F. Prophenoloxidase of Odontotermes formosanus (Shiraki) (Blattodea: Termitidae) Is a Key Gene in Melanization and Has a Defensive Role during Bacterial Infection. International Journal of Molecular Sciences. 2023; 24(1):406. https://doi.org/10.3390/ijms24010406

Chicago/Turabian Style

Wang, Zhiqiang, Jian Luo, Kai Feng, Yujingyun Zhou, and Fang Tang. 2023. "Prophenoloxidase of Odontotermes formosanus (Shiraki) (Blattodea: Termitidae) Is a Key Gene in Melanization and Has a Defensive Role during Bacterial Infection" International Journal of Molecular Sciences 24, no. 1: 406. https://doi.org/10.3390/ijms24010406

APA Style

Wang, Z., Luo, J., Feng, K., Zhou, Y., & Tang, F. (2023). Prophenoloxidase of Odontotermes formosanus (Shiraki) (Blattodea: Termitidae) Is a Key Gene in Melanization and Has a Defensive Role during Bacterial Infection. International Journal of Molecular Sciences, 24(1), 406. https://doi.org/10.3390/ijms24010406

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