Next Article in Journal
Novel Microwave-Assisted Synthesis of the Immunomodulator Organotellurium Compound Ammonium Trichloro(dioxoethylene-O,O')tellurate (AS101)
Previous Article in Journal
miR-137: A New Player in Schizophrenia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Bactrocera dorsalis Serine Proteases and Evidence for Their Indirect Role in Insecticide Tolerance

Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University, Chongqing 400716, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2014, 15(2), 3272-3286; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15023272
Submission received: 20 December 2013 / Revised: 9 February 2014 / Accepted: 12 February 2014 / Published: 21 February 2014
(This article belongs to the Section Biochemistry)

Abstract

:
The oriental fruit fly Bactrocera dorsalis (Hendel) causes devastating losses to agricultural crops world-wide and is considered to be an economically important pest. Little is known about the digestive enzymes such as serine proteases (SPs) in B. dorsalis, which are important both for energy supply and mitigation of fitness cost associated with insecticide tolerance. In this study, we identified five SP genes in the midgut of B. dorsalis, and the alignments of their deduced amino acid sequences revealed the presence of motifs conserved in the SP superfamily. Phylogenetic analyses with known SPs from other insect species suggested that three of them were trypsin-like proteases. Analyses of the expression profiles among the different developmental stages showed that all five genes were most abundant in larvae than in other stages. When larvae were continuously fed on diet containing 0.33 μg/g β-Cypermethrin, expression of all five genes were upregulated in the midgut but the larval development was delayed. Biochemical assays were consistent with the increased protease activity exhibited by SPs in the midgut after treatment with β-Cypermethrin. Taken together, these findings provide evidence for the hypothesis that enhanced SP activity may play an indirect role in relieving the toxicity stress of insecticide in B. dorsalis.

Graphical Abstract

1. Introduction

The oriental fruit fly, Bactrocera dorsalis (Hendel), is one of the most destructive pests in the world and causes huge financial losses worldwide [1]. Damages caused by the oriental fruit fly include both punctures into the fruits and vegetables during oviposition and feeding on the fruit pulp by the developing larvae [2]. Since the biological and ecological traits of B. dorsalis include invasive ability (high mobility, dispersive powers, and fecundity), wild distribution, and economic significance, it is ranked high on the quarantine target lists [3].
The insect digestive proteases catalyze the dietary protein to release the free amino acids and supply essential nutrients for growth and development. All insect digestive proteases can be classified into three classes: aspartic (also called carboxyl or acid), cysteine (also called thiol), and serine proteases [4]. Among them, serine proteases (SPs) are considered to be important proteolytic enzymes. Most insects were presumed to use SPs with trypsin or chymotrypsin-like specificities to digest their dietary protein [5]. It is known that SPs represent the major proteolytic enzymes in the midgut of many lepidopterans and dipterans [6]. In some lepidopteran insects, SPs are known to dominate in the larval gut environment and contribute to about 95% of the total digestive activity [7]. They are also dominant proteolytic enzymes in the midgut of Drosophila [8].
SPs are a diverse group of proteolytic enzymes. Almost one-third of all proteases can be classified as SPs, named for the nucleophilic Ser residue at the active site. SPs clan of endoproteolytic enzymes can be further divided into families and subfamilies. The S1.A subfamily includes trypsin, chymotrypsins, elastases, and some recently identified serine collagenases [9], which are involved in a number of critical physiological processes, such as digestion, hemostasis, apoptosis, signal transduction, reproduction, and the immune response [10].
Due to the high effectiveness, low toxicity to non-target organisms, and easy biodegradability, pyrethroids are used more preferably over organochlorines and organophosphates [11]. Most pyrethroids contain cyclopropane carboxylic acid (or a moiety equivalent group) linked to aromatic alcohols through a central ester (or ether) bond. Modifications to this basic structure are designed to increase insecticidal potency or photostability [12]. They act on the insect nervous system by slowing action potential decay, and thereby initiating repetitive discharges in motor and sensory axons. Electrophysiological studies suggested that these phenomena result from modification of gating kinetic of neuronal, voltage-sensitive sodium channels [13]. In insects, the metabolic enzymes that are implicated in the detoxification of pyrethroids include cytochrome P450 monooxygenases (P450s), esterases, and glutathione S-transferases (GSTs) [14]. Our previous studies found that the overexpression of some GSTs and P450 genes in B. dorsalis following exposure to the widely used pyrethroid, β-Cypermethrin [1517], indicating that both GSTs and P450s play important roles in the stress response of B. dorsalis to β-Cypermethrin. In addition, the target enzyme such as acetylcholinesterase (AChE) upon the exposure to β-Cypermethrin in B. dorsalis was also well documented [18]. However, the potential insecticide modulation of digestive and energy metabolic enzyme has been an object of little attention. Previous studies have found that in insects the digestive and metabolic energy enzymes play important roles in mitigating fitness cost associated with insecticide resistance [1922], and it has been confirmed that some xenobiotics could modulate the activity of not only the detoxification enzymes, but also modify the energy and digestive enzymes in an unspecified manner [23,24]. Detoxification is an energy consuming process, it is plausible that the energy utilized originally for normal development and reproduction may be disturbed and channeled toward detoxification. Therefore, we hypothesized that exposure to insecticides may enhance protein digestion in insects resulting in additional energy, and increased availability of amino acids for protein synthesis. This in turn will supply the materials to produce the detoxification apparatus without costing energy for other physiological processes. In this study, we identified five SP genes in the midgut of B. dorsalis and determined their development–specific expressions and transcriptional profiles in larvae upon exposure to β-Cypermethrin. In addition, the activities of SPs in the midgut of B. dorsalis larvae were also measured. The results might provide some evidence for the indirect role of SPs in insecticide tolerance.

2. Results

2.1. β-Cypermethrin Intake Slows Larval Development

Feeding β-Cypermethrin containing diet decreased the average length of 3rd instar larvae from 0.85 ± 0.028 to 0.67 ± 0.021 cm (Figure 1A). Similarly, β-Cypermethrin intake also decreased the average larval body weight from 13.98 ± 0.91 to 11.12 ± 1.09 mg (Figure 1B). In addition, severe effects on overall larval morphology were observed when the insects were fed with β-Cypermethrin containing diets (Figure 1C).

2.2. Effect of β-Cypermethrin on Protease Activities

Activities of SPs were determined using BApNA and L-TAME as substrates for amidolytic and esterolytic activities, respectively. The results showed that in both processes, SP activities in the treatment group were significantly higher than those in the control group. For BApNA hydrolysis, the specific activity of SPs from the treatment group was 450.00 ± 32.97 nmol/mL/min/mg, which was significantly higher than that in the control group (302.67 ± 79.93 nmol/mL/min/mg) (p < 0.05; Table 1). Similarly, for the L-TAME hydrolysis, the specific activity of SPs from the treatment group was 266.00 ± 13.89 nmol/mL/min/mg, which was significantly higher than that in the control group (201.33 ± 16.92 nmol/mL/min/mg) (p < 0.05; Table 1).

2.3. Amino Acid Similarities

The full-length cDNA sequences of all five SPs consisted of 744–798 nucleotides that encoded proteins with 247–265 amino acid residues and were named as BdSP1, BdSP2, BdSP3, BdSP4, and BdSP5, respectively. These sequences were deposited in GenBank and the relative accession numbers were: GAAP01000017, GAAP01000019, GAAP01000020, GAAP01000021, and GAAP01000022.
Protein sequence alignments of BdSPs with other known insect SPs revealed the presence of the hallmark SP feature, which is the catalytic triad consisting of His57, Asp102, and Ser195 (chymotrypsin numbering) with in the highly conserved sequences of TAAHC, DIAL, and GDSGGP, respectively, in most SPs (Table 2). These motifs were easily identified in all the predicted amino acid sequences.
During the analysis we also found that 14 among 21 SP sequences including BdSP1, BdSP2 and BdSP4, and other known insect trypsins contained Asp190, Gly216, and Gly226, while the remaining sequences (BdSP3 and BdSP5) contained Gly, Val and Asp at positions 190, 216, 226, respectively (Figure 2).

2.4. Phylogenetic Analyses

To analyze the sequence homology and phylogenetic relationships, sixteen other insect SPs were downloaded from GenBank and aligned with the five BdSP gene sequences (Table 2). The phylogenetic tree constructed using the maximum likelihood method had two major clades (Figure 3) with all known trypsins along with BdSP1, BdSP2, and BdSP4 clustering on the same clade. BdSP1 was most closely related to DmTR1 and DeTR1 from Drosophila melanogaster and D. erecta, respectively. BdSP4 was in clade A and shared a common branch with BdSP1. Similarly, BdSP3, BdSP5 and other SPs sequences were clustered together in clade B indicting a high sequence similarity. In clade B, BdSP3 and BdSP5 were clustered with AaSP1 from Aedes aegypti and CcSP from Ceratitis capitata, respectively, with 100% bootstrap. Overall the topologies and distances resolved by the maximum likelihood analyses were also supported by the distance-neighbor-joining analyses (data not shown).

2.5. Developmental Expression Profiles of BdSPs

qRT-PCR results revealed that all five SPs were highly expressed in the larvae, followed by the adults. The presence of BdSPs in both eggs and pupae were also detected, but at lower levels. Among the five genes, the expression level of BdSP2 was the highest in larvae and adult (Table 3).

2.6. Temporal Expression of BdSPs in Response to β-Cypermethrin

The transcriptional changes of all five SP genes transcripts in larvae exposed to β-Cypermethrin determined by qRT-PCR which showed significant upregulation of BdSP3, BdSP5, BdSP1, BdSP2, and BdSP4 in larvae following the treatment with β-Cypermethrin, especially for BdSP3. When compared to the control, all five genes registered 10.9-, 2.6-, 1.2-, 1.4- and 1.6-fold increase, respectively (Figure 4). Statistical analyses showed that among the genes, the expression of two genes (BdSP3 and BdSP5) were significantly upregulated compared to the control (p < 0.05).

3. Discussion

SPs are dominant proteolytic enzymes in the insect midgut [8] and can be distinguished by the presence of the Asp-His-Ser “charge relay” system or “catalytic triad” [9]. They are present in all phylogenetic kingdoms including viruses and are involved in many physiological processes [10]. Particularly, the trypsins and chymotrypsins play important roles in insect food digestion, immune defense, and zymogen activation [26]. In addition, SPs and SP homologs constitute the second largest family of genes in the D. melanogaster genome [27]. To date, the genes encoding trypsins and chymotrypsins have been cloned and characterized in several insect species [2729]. This study is the first report on biochemical and molecular characterization of SPs in B. dorsalis.
It is well recognized that toxic substances may induce many biochemical changes preceding cellular and systemic dysfunction in animals [30]. In this study, we found the overexpression of SPs and the enhanced SPs activity in B. dorsalis upon the exposure of β-Cypermethrin, which may probably be explained in this way: the detoxification process is known to draw energy allocated to normal development and therefore, absence of energy supplement may impair the normal development and reproduction in insects [31,32]. To relieve the effect of insecticide, when treated with β-Cypermethrin insects might enhance protein digestion. The enhanced hydrolytic enzyme activity will supply additional energy and amino acids for protein synthesis. This will supply the materials to produce the detoxification apparatus for reducing the cost of energy originally utilized by normal physiological processes.
The similar complementary mechanisms have also been observed in some other insects confronted under environmental stress. It has been reported that there was an increase in carbohydrate metabolism induced by organophosphorus insecticide in silkworm. A significant decrease in pyruvate levels and increase in lactate levels were observed in the hemolymph and fat body of silkworm, suggesting the utilization of pyruvate through the tricarboxylic acid (TCA) cycle to meet the required energy demands to counteract the insecticide toxicity [24]. In a similar study, the sharp elevation in protease activity and free amino acid levels were also observed in the hemolymph and fat body of silkworms when exposed to lethal and sublethal doses of fenitrothion and ethion [33]. Besides, it was found that the enhanced trypsins and cysteine protease activities play important roles in mitigating physiological costs associated with the maintenance of insecticide resistance in studies on maize weevil [19,20]. Taken together, it is reasonable to assume that the higher SP activities in B. dorsalis might play an indirect role in mitigating the stress caused by β-Cypermethrin.
We obtained the full-length cDNA sequences of five SPs on the base of B. dorsalis midgut transcriptome [34]. All the five SPs contained six conserved cysteine residues that correspond to the three disulfide bonds characteristic of arthropod serine proteases [35]. They also had the catalytic triad “Asp-His-Ser”, indicating that these putative proteins were indeed SPs. To predict the substrate specificity of the SPs, we identified the primary substrate–binding pocket (residues 190, 216 and 226 in chymotrypsin) in the aligned sequences [25]. SPs that contain Asp190, Gly216, and Gly/Ala/Ser226 are predicted to have trypsin-like specificity. Therefore we concluded that BdSP1, BdSP2 and BdSP4 are trypsin-like SPs. (Table 4). Although the result of sequence analysis showed that BdSP3 (69%) may be trypsin-like SPs, and BdSP5 (41%) may belong to chymotrypsin-like SPs, our protein alignment did not reveal the hallmark of trypsin (Asp190, Gly216 or Gly/Ala/Ser226) or chymotrypsin (Ser189, Gly216, and Gly226); rather Gly, Val and Asp, were replaced in these positions, respectively. Therefore, we could not unequivocally conclude that all BdSPs are trypsin-like SPs. The final classification of BdSP3 and BdSP5 is required for further clarification.
The determination of SPs’ activities provided some evidences for our hypothesis, but we also determined the expression profiles of five SPs genes after exposure to β-Cypermethrin by qRT-PCR. Consistent with the enhanced SP activities, we observed that there was a significant increase in expression levels of SP transcripts. Although the increase was not dramatic, it still suggested the trend of upregulation of SPs genes, and probably there are other SP genes playing important roles in the process. It was the desirable result that there was no significant difference in larval morphology between control and treatment, which would be more in favor of our hypothesis. However, in fact, after the treatment of β-Cypermethrin, the development of B. dorsalis larvae was delayed. Maybe it can be explained in this way that the stress responses of insects are complicated processes, and the internal changes we found might be not great enough to directly reflect on larval morphology and completely cover the toxicity stress caused by β-Cypermethrin. The five genes were highly expressed in the midgut of larvae suggesting a high energy requirement at this stage, which is consistent with the developmental characteristics of B. dorsalis larvae that feed actively on host fruits and vegetables. The accumulation of nutrition and energy at the larval stage guarantees normal development and reproduction in B. dorsalis. After larval hatching, the enzymes associated with digestion and absorption was highly expressed to satisfy the strong demand for nutrition and energy of B. dorsalis larvae. This is also the reason why we chose the larval stage of B. dorsalis in our study. As expected, the upregulation of these five SPs genes and the increase of serine-protease activity after exposure to β-Cypermethrin demonstrate that SPs do indirectly mitigate the toxicity stress of β-Cypermethrin in B. dorsalis.

4. Experimental Section

4.1. Insects

The laboratory stock of B. dorsalis, originally collected in 2008 from Hainan province, China, was cultured in the insectary at 27 ± 1 °C, 70% ± 5% relative humidity, 14 h light: 10 h dark photoperiod. Larvae were reared on an artificial diet described previously [34].

4.2. Exposure to β-Cypermethrin

About fifty newly emerged B. dorsalis larvae were reared on an artificial diet containing 0.33 μg/g β-Cypermethrin or the control diet with acetone (solvent used to dissolve β-Cypermethrin) from 1st to 3rd instar, which took about five days. At the end of 3rd instar, body length and mass of the surviving larvae were weighed in both control and treated groups. All experimental groups were maintained under the same conditions (27 ± 1 °C, 70% ± 5% relative humidity, and 14 h:10 h (L:D) photoperiod) in culture dishes. The choice of β-Cypermethrin concentration (0.33 μg/g in artificial diet) was based on our previous study [36]. The β-Cypermethrin (95%) was provided by Institute for Control of Agrochemicals of Sichuan Province, Chengdu, China.

4.3. Preparation of Enzyme Extracts

Enzyme extracts were prepared from the midguts of larvae by dissecting on ice under a stereomicroscope (Olympus SZX12, Tokyo, Japan), Treatment and control group both contained three biological replications, and 30 larvae midguts were collected in each sample. The dissected midguts were frozen in liquid nitrogen, homogenized in 0.1 M Tris-HCl (pH 8.2 containing 10 mM (CaCl2) and 2 mL of the homogenate or crude extract was centrifuged at 12 000g for 5 min at 4 °C and filtered through a 0.25 μm filter membrane. The supernatant was centrifuged again at 17,500g for 25 min at 4 °C, and 2 mL of the final supernatant was collected and used as the enzyme source.

4.4. Protein Concentration and Protease Activity

Protein content in the enzyme extracts was determined based on described methods using bovine serum albumin as the standard [37]. Activities of SPs were determined using N-α-benzoyl-dl-arginine p-nitroanilide (BApNA) and N-α-p-tosyl-l-Arg methyl ester (l-TAME) as substrates for amidolytic and esterolytic activities, respectively. Amidolytic activity was measured with BApNA as substrate [32], at a final concentration of 1.7 mM in glycine sodium-hydroxide buffer (pH 10.5) containing 7.4% dimethylformamidine (DMF, v/v) at 37 °C [38,39]. Esterolytic activity was measured with l-TAME; 1 mM as substrate in 0.1 M Tris-HCl buffer (pH 8.2) at 25 °C [40]. Each assay in the control and treated groups was repeated three times.

4.5. RNA Isolation and First-Strand cDNA Synthesis

Total RNA was extracted from the midgut of control and treated larvae using an RNAeasy Micro Kit (Qiagen, Hilden, Germany). Total RNAs from whole bodies of insects were isolated at each developmental stage (eggs, larvae, pupae, and adults) using TRIzol reagent (Biomed, Beijing, China). About 30 individuals were dissected from each of the control and treated groups, and each developmental stage. In addition, three biological replicates were maintained per sample.
The midguts were dissected under a stereomicroscope (Olympus SZX12, Tokyo, Japan) and placed in a 1.5-mL centrifuge tube containing the RNA store reagent (Tiangen, Beijing, China). The collected midguts and insect bodies from different life stages were homogenized immediately in liquid nitrogen. RNA was quantified by measuring the absorbance at 260 nm using a NanoVue UV-Vis spectrophotometer (GE Healthcare Bio-Science, Uppsala, Sweden). The purity of all RNA samples was assessed based on absorbance ratio at OD260/280 and OD260/230, and the integrity of RNA was checked on a 1% agarose gel. The RNA extraction using the RNAeasy Micro Kit included an on-column genomic DNA elimination step that efficiently removed genomic DNA from the total RNA preparations. The total RNA extracted by TRIzol reagent was treated separately by DNase I (Promega, Madison, WI, USA). First strand cDNA was synthesized using a PrimeScript® RT reagent Kit (Takara, Dalian, China).

4.6. Sequence Analysis and Phylogenetic Tree Construction

Five SP gene sequences were obtained from the midgut transcripts of B. dorsalis larvae [34]. DNA sequences were edited with DNAMAN 5.2.2 (Lynnon Pointe-Claire, Quebec, Canada) and the deduced amino acid sequences of all five SPs genes were trimmed and aligned using Clustal X [41]. Homologous SP protein sequences were searched by BLAST in the NCBI database ( http://blast.ncbi.nlm.nih.gov/Blast.cgi). Phylogenetic tree was constructed with MEGA 5.0 using the Maximum likelihood method. Bootstrap values were calculated with 1000 replications. The MsTR of Manduca sexta was used as outgroup in constructing phylogenetic tree.

4.7. Analysis of SP Gene Expression

We analyzed the expression patterns of all five SP genes in each developmental stage (eggs, larvae, pupae, and adults) of B. dorsalis and also their expression in both the control and treatment group by quantitative real-time-PCR (qRT-PCR). Based on standard protocols, α-tubulin and EF1α were selected as reference genes for the analyses [42].
qRT-PCR was conducted on a Stratagene Mx3000P System in a total reaction volume of 20 μL containing 500 ng cDNA, 10 μL GoTaq® qPCR Master Mix (Promega, Madison, WI, USA), 10 pmol of each primer (Table 4), and double distilled water. All reactions were performed under the following conditions: 95 °C for 30 s, 40 cycles of 95 °C for 15 s, 60 °C for 30 s, and 72 °C for 30 s. To verify the specificity of the amplicon for each primer pair, a dissociation curve was included from 60–95°C at the end of each qRT-PCR run. A cDNA dilution series (1, 1/3, 1/9, 1/27, and 1/81) with the sample cDNA was used to construct the standard curve and calculate the efficiency of amplification. The PCR efficiencies of each SP gene and the references were calculated using Mxpro-Mx3000P version 3.20 (Agilent-Stratagene, Santa Clara, CA, USA). Relative expression levels were calculated by the comparative CT method [43].

4.8. Statistical Analysis

Differences in gene expression levels among the developmental stages were analyzed by one-way analysis of variance (ANOVA), followed by Duncan’s multiple range tests in SPSS 17.0 (SPSS Armonk, New York, NY, USA). All other data were analyzed with independent sample t-test for significance (p < 0.05) with SPSS 17.0.

5. Conclusions

In summary, the increases of SP gene expression and enzyme activity after exposure to β-Cypermethrin may reflect the potential role of SPs in relieving the toxicity stress of β-Cypermethrin. Additional studies of SP genes will help elucidate their possible role in different life stages at the molecular level.

Acknowledgments

This research was supported in part by the Natural Science Foundation of Chongqing (CSTC, 2013jjB0176), the Special Fund for Agro-scientific Research in the Public Interest (201203038), the Program for Changjiang Scholars and Innovative Research Teams in Universities (IRT0976), the earmarked fund for Modern Agro-industry (Citrus) Technology Research System, and the Fundamental Research Funds for the Central Universities (XDJK2013A017) of China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bateman, M.A. The ecology of fruit flies. Annu. Rev. Entomol 1972, 17, 493–518. [Google Scholar]
  2. Fletcher, B. The biology of dacine fruit flies. Annu. Rev. Entomol 1987, 32, 115–144. [Google Scholar]
  3. Clarke, A.R.; Armstrong, K.F.; Carmichael, A.E.; Milne, J.R.; Raghu, S.; Roderick, G.K.; Yeates, D.K. Invasive phytophagous pests arising through a recent tropical evolutionary radiation: The Bactrocera dorsalis complex of fruit flies. Annu. Rev. Entomol 2005, 50, 293–319. [Google Scholar]
  4. Wolfson, J.L.; Murdock, L.L. Diversity in digestive proteinase activity among insects. J. Chem. Ecol 1990, 16, 1089–1102. [Google Scholar]
  5. Murdock, L.; Shade, R.; Pomeroy, M. Effects of E-64, a cysteine proteinase inhibitor, on cowpea weevil growth, development, and fecundity. Environ. Entomol 1988, 17, 467–469. [Google Scholar]
  6. Chapman, R.F. The Insects: Structure and Function, 4th ed.; Cambridge University Press: New York, NY, USA, 1998; pp. 54–59. [Google Scholar]
  7. Srinivasan, A.; Giri, A.P.; Gupta, V.S. Structural and functional diversities in lepidopteran serine proteases. Cell Mol. Biol. Lett 2006, 11, 132–154. [Google Scholar]
  8. Li, H.M.; Buczkowski, G.; Mittapalli, O.; Xie, J.; Wu, J.; Westerman, R.; Schemerhorn, B.; Murdock, L.; Pittendrigh, B. Transcriptomic profiles of Drosophila melanogaster third instar larval midgut and responses to oxidative stress. Inset. Mol. Biol 2008, 17, 325–339. [Google Scholar]
  9. Blow, D.M. The tortuous story of Asp-His-Ser: Structural analysis of alpha-chymotrypsin. Trends. Biochem. Sci 1997, 22, 405. [Google Scholar]
  10. Hedstrom, L. Serine protease mechanism and specificity. Chem. Rev 2002, 102, 4501–4524. [Google Scholar]
  11. Kale, M.; Rathore, N.; John, S.; Bhatnagar, D. Lipid peroxidative damage on pyrethroid exposure and alterations in antioxidant status in rat erythrocytes: A possible involvement of reactive oxygen species. Toxicol. Lett 1999, 105, 197–205. [Google Scholar]
  12. Wolansky, M.; Harrill, J. Neurobehavioral toxicology of pyrethroid insecticides in adult animals: A critical review. Neurotoxicol. Teratol 2008, 30, 55–78. [Google Scholar]
  13. Vais, H.; Williamson, M.S.; Devonshire, A.L.; Usherwood, P.N.R. The molecular interactions of pyrethroid insecticides with insect and mammalian sodium channels. PestManag. Sci 2001, 57, 877–888. [Google Scholar]
  14. Li, X.; Schuler, M.A.; Berenbaum, M.R. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu. Rev. Entomol 2007, 52, 231–253. [Google Scholar]
  15. Huang, Y.; Shen, G.M.; Jiang, H.B.; Jiang, X.Z.; Dou, W.; Wang, J.J. Multiple P450 genes: Identification, tissue-specific expression and their responses to insecticide treatments in the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidea). Pestic. Biochem. Physiol 2013, 106, 1–7. [Google Scholar]
  16. Hu, F.; Dou, W.; Wang, J.J.; Jia, F.X.; Wang, J.J. Multiple glutathione S-transferase genes: identification and expression in oriental fruit fly Bactrocera dorsalis. PestManag. Sci 2014, 70, 295–303. [Google Scholar]
  17. Hu, F.; Dou, W.; Wang, J.J.; Jia, F.X.; Wang, J.J. Purification and biochemical characterization of glutathione S-transferases from four field populations of Bactrocera dorsalis (Hendel) (Diptera: Tephritidae). Arch. Insect Biochem. Physiol 2011, 78, 201–215. [Google Scholar]
  18. Hsu, J.C.; Wu, W.J.; Haymer, D.S.; Liao, H.Y.; Feng, H.T. Alterations of the acetylcholinesterase enzyme in the oriental fruit fly Bactrocera dorsalis are correlated with resistance to the organophosphate insecticide fenitrothion. Insect Biochem. Mol. Biol 2008, 38, 146–154. [Google Scholar]
  19. Silva, L.B.; Reis, A.P.; Pereira, E.J.G.; Oliveira, M.G.A.; Guedes, R.N.C. Altered cysteine proteinase activity in insecticide-resistant strains of the maize weevil: purification and characterization. Comp. Bioche. Physiol. B 2010, 157, 80–87. [Google Scholar]
  20. Silva, L.B.; Reis, A.P.; Pereira, E.J.G.; Oliveira, M.G.A.; Guedes, R.N.C. Partial purification and characterization of trypsin-like proteinases from insecticide-resistant and -susceptible strains of the maize weevil Sitophilus zeamais. Comp. Bioche. Physiol. B 2010, 155, 12–19. [Google Scholar]
  21. Philippou, D.; Field, L.; Moores, G. Metabolic enzyme (s) confer imidacloprid resistance in a clone of Myzus persicae (Sulzer)(Hemiptera: Aphididae) from Greece. Pest Manag. Sci 2010, 66, 390–395. [Google Scholar]
  22. Araújo, R.; Guedes, R.N.C.; Oliveira, M.G.A.; Ferreira, G. Enhanced proteolytic and cellulolytic activity in insecticide-resistant strains of the maize weevil Sitophilus zeamais. J. Stored Prod. Res 2008, 44, 354–359. [Google Scholar]
  23. Natsuhara, K.; Shimada, K.; Tanaka, T.; Miyata, T. Phenobarbital induction of permethrin detoxification and phenobarbital metabolism in susceptible and resistant strains of the beet armyworm Spodoptera exigua (Hübner). Pestic. Biochem. Physiol 2004, 79, 33–41. [Google Scholar]
  24. Nath, B.S. Changes in carbohydrate metabolism in hemolymph and fat body of the silkworm, Bombyx mori L., exposed to organophosphorus insecticides. Pestic. Biochem. Physiol 2000, 68, 127–137. [Google Scholar]
  25. Perona, J.J.; Craik, C.S. Structural basis of substrate specificity in the serine proteases. Protein Sci 1995, 4, 337–360. [Google Scholar]
  26. Ge, Z.Y.; Wan, P.J.; Han, Z.J.; Golding, B. Cloning and characterization of trypsin-and chymotrypsin-like genes in the striped rice stem borer Chilo suppressalis. Genome 2012, 55, 281–288. [Google Scholar]
  27. Ross, J.; Jiang, H.; Kanost, M.R.; Wang, Y. Serine proteases and their homologs in the Drosophila melanogaster genome: An initial analysis of sequence conservation and phylogenetic relationships. Gene 2003, 304, 117–131. [Google Scholar]
  28. Zhan, Q.; Zheng, S.; Feng, Q.; Liu, L. A midgut-specific chymotrypsin cDNA (Slctlp1) from Spodoptera litura: Cloning, characterization, localization and expression analysis. Arch. Insect Biochem. Physiol 2011, 76, 130–143. [Google Scholar]
  29. Zhang, C.; Zhou, D.; Zheng, S.; Liu, L.; Tao, S.; Yang, L.; Hu, S.; Feng, Q. A chymotrypsin-like serine protease cDNA involved in food protein digestion in the common cutworm, Spodoptera litura: Cloning, characterization, developmental and induced expression patterns, and localization. J. Insect Physiol 2010, 56, 788–799. [Google Scholar]
  30. Bhosale, S.; Kallapur, V. Changes in the metabolic fuel reserves of the 5th instar, Bombyx mori, following endosulfan treatment. Entomon 1985, 10, 281. [Google Scholar]
  31. Karban, R.; Agrawal, A.A. Herbivore offense. Annu. Rev. Ecol. Syst 2002, 33, 641–664. [Google Scholar]
  32. Appel, H.M.; Martin, M.M. Significance of metabolic load in the evolution of host specificity of Manduca sexta. Ecology 1992, 73, 216–228. [Google Scholar]
  33. Nath, B.S.; Suresh, A.; Varma, B.M.; Kumar, R. Changes in protein metabolism in hemolymph and fat Body of the silkworm, Bombyx mori (Lepidoptera: Bombycidae) in response to organophosphorus insecticides toxicity. Ecotoxicol. Environ. Saf 1997, 36, 169–173. [Google Scholar]
  34. Shen, G.M.; Dou, W.; Huang, Y.; Jiang, X.Z.; Smagghe, G.; Wang, J.J. In silico cloning and annotation of genes involved in the digestion, detoxification and RNA interference mechanism in the midgut of Bactrocera dorsalis [Hendel (Diptera: Tephritidae)]. Insect Mol. Biol 2013, 22, 354–365. [Google Scholar]
  35. Davis, C.A.; Riddell, D.C.; Higgins, M.J.; Holden, J.J.; White, B.N. A gene family in Drosophila melanogaster coding for trypsin-like enzymes. Nucleic. Acids. Res 1985, 13, 6605–6619. [Google Scholar]
  36. Shen, G.M.; Huang, Y.; Jiang, X.Z.; Dou, W.; Wang, J.J. Effect of β-Cypermethrin exposure on the stability of nine housekeeping genes in Bactrocera dorsalis (Diptera: Tephritidae). Fla. Entomol 2013, 96, 442–450. [Google Scholar]
  37. 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]
  38. Hosseininaveh, V.; Bandani, A.; Azmayeshfard, P.; Hosseinkhani, S.; Kazzazi, M. Digestive proteolytic and amylolytic activities in Trogoderma granarium Everts (Dermestidae: Coleoptera). J. Stored Prod. Res 2007, 43, 515–522. [Google Scholar]
  39. Chen, C.C.; Gao, L.T.; Shi, B.J.; Wen, W.D.; Hu, Z.N.; Wu, W.J. Effects of periplocosides P and E from Periploca sepium on the proteinase activities in the midgut of larvae of Mythimna separata and Agrotis ypsilon (Lepidoptera Noctuidae). Acta. Entomol. Sin 2012, 55, 412–419. [Google Scholar]
  40. Hummel, B.C. A modified spectrophotometric determination of chymotrypsin, trypsin, and thrombin. Can. J. Biochem. Physiol 1959, 37, 1393–1399. [Google Scholar]
  41. Jeanmougin, F.; Thompson, J.D.; Gouy, M.; Higgins, D.G.; Gibson, T.J. Multiple sequence alignment with Clustal X. Trends. Biochem. Sci 1998, 23, 403. [Google Scholar]
  42. Shen, G.M.; Jiang, H.B.; Wang, X.N.; Wang, J.J. Evaluation of endogenous references for gene expression profiling in different tissues of the oriental fruit fly Bactrocera dorsalis (Diptera: Tephritidae). BMC Mol. Biol 2010, 11, 76. [Google Scholar]
  43. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic. Acids Res 2001, 29, e45. [Google Scholar]
Figure 1. Body length (A), weight (B) and morphological characteristics (C) of 3rd instar larvae treated with (treatment) or without (control) β-Cypermethrin. * indicates significant difference between treatment and control (Independent sample t-test, p < 0.05).
Figure 1. Body length (A), weight (B) and morphological characteristics (C) of 3rd instar larvae treated with (treatment) or without (control) β-Cypermethrin. * indicates significant difference between treatment and control (Independent sample t-test, p < 0.05).
Ijms 15 03272f1
Figure 2. Protein alignment of insect SPs. Protein alignment of all BdSPs used in this analysis from positions 180–240. Residues of importance are represented as follows: (▼) Ser195 catalytic triad residues, (★) accessory catalytic residues, (●) the third and final disulfide bridge and (○) Asp194 where position 1 (Ile/Val) is buried after activation of the mature peptide.
Figure 2. Protein alignment of insect SPs. Protein alignment of all BdSPs used in this analysis from positions 180–240. Residues of importance are represented as follows: (▼) Ser195 catalytic triad residues, (★) accessory catalytic residues, (●) the third and final disulfide bridge and (○) Asp194 where position 1 (Ile/Val) is buried after activation of the mature peptide.
Ijms 15 03272f2
Figure 3. The maximum likelihood phylogenetic tree of SPs. The amino acid segments aligned in Figure 3 were used in this analysis. Clade A corresponds with the known trypsin and BdSP1, BdTSP2 and BdSP4, clade B represents other SPs along with BdSP3 and BdSP5.
Figure 3. The maximum likelihood phylogenetic tree of SPs. The amino acid segments aligned in Figure 3 were used in this analysis. Clade A corresponds with the known trypsin and BdSP1, BdTSP2 and BdSP4, clade B represents other SPs along with BdSP3 and BdSP5.
Ijms 15 03272f3
Figure 4. Transcription profiles of the five different B. dorsalis SPs (BdSPs) in larvae exposed to β-Cypermethrin. Relative expression was calculated as 2 [ CT ( B d S P C ) - CT ( B d S P T ) ] 2 [ CT ( E F C ) - CT ( E F T ) ]. BdSPT and BdSPC represent target genes from the treatment and control groups, respectively. EF represents reference gene (EF1α). * indicates significant difference in expression level of target genes between control and treatment (Independent sample t-test, p < 0.05).
Figure 4. Transcription profiles of the five different B. dorsalis SPs (BdSPs) in larvae exposed to β-Cypermethrin. Relative expression was calculated as 2 [ CT ( B d S P C ) - CT ( B d S P T ) ] 2 [ CT ( E F C ) - CT ( E F T ) ]. BdSPT and BdSPC represent target genes from the treatment and control groups, respectively. EF represents reference gene (EF1α). * indicates significant difference in expression level of target genes between control and treatment (Independent sample t-test, p < 0.05).
Ijms 15 03272f4
Table 1. Specific activity of SPs in the midgut of 3rd instar larva of B. dorsalis.
Table 1. Specific activity of SPs in the midgut of 3rd instar larva of B. dorsalis.
GroupsSpecific activity (nmol/mL/min/mg)

BApNAL-TAME
Control302.67 ± 79.93201.33 ± 16.92
Treatment450.00 ± 32.97 *266.00 ± 13.89 *
Each value represents the mean (M ± SE) of three replications.
*represents significant difference (Independent sample t-test, p < 0.05).
Table 2. List of serine proteases used in the alignment and phylogenetic analyses.
Table 2. List of serine proteases used in the alignment and phylogenetic analyses.
SpeciesSerine proteaseGenBank Accession numberConserved regions aLengthEnzyme specificity b

TAAHCDIALGDSGGP
Bactrocera. dorsalisBdSP1GAAP01000017SATHCDIGL265T(DGG)
Bactrocera. dorsalisBdSP2GAAP01000019DIAI259T(DGG)
Bactrocera. dorsalisBdSP3GAAP01000020DLAL247E(GVD)
Bactrocera. dorsalisBdSP4GAAP01000021258T(DGG)
Bactrocera. dorsalisBdSP5GAAP01000022257E(GVD)
Glossina morsitans morsitansGmmTR1ACB98719.1DVAV256T(DGG)
Glossina morsitans morsitansGmmTR2ADD19605.1DYSL281T(DGG)
Aedes aegyptiAaSP1XP_004520098.1DLAL248E(GVD)
Aedes aegyptiAaSP2XP_001659492.1247E(GVD)
Aedes aegyptiAaSP3XP_001659851.1DIAV274E(GVD)
Ceratitis capitataCcSPXP_004520096.1DVAL254E(GVD)
Ceratitis capitataCcTRXP_004517776.1DIAV257T(DGG)
Drosophila melanogasterDmTR1AAA17449.1DIVI253T(DGG)
Drosophila melanogasterDmTR2NP_525112.1DIAI256T(DGG)
Drosophila erectaDeTR1XP_001976081.1DIAV256T(DGG)
Drosophila erectaDeTR2XP_001976082.1DVGI262T(DGG)
Drosophila yakubaDyTRXP_002091228.1DIAI256T(DGG)
Phlebotomus papatasiPpTRAAM96943.1DFSL268T(DGW)
Manduca sextaMsTRP35046.1DIAI256T(DGG)
Hypoderma dianaHdSPACF98290.1DVAI256T(DGG)
Stomoxys calcitransSsSPAAC39131.1DVAV254T(DGG)
aIf not listed, sequences are identical to the conserved TAAHC, DIAL, or GDSGGP;
benzyme specificity predicted based on Perona and Craik (1995) [25] (SP, Serine Protease; TR, trypsin).
Table 3. Expression patterns of five BdSPs in different life stages of B. dorsalis.
Table 3. Expression patterns of five BdSPs in different life stages of B. dorsalis.
GenesRelative expression

EggLarvaPupaAdult
BdSP11.07.48 ± 0.930.05 ± 0.021.61 ± 0.51
BdSP21.0189.62 ± 23.110.47 ± 0.16118 ± 28.69
BdSP31.04.04 ± 0.860.45 ± 0.222.21 ± 0.03
BdSP41.06.71 ± 1.720.67 ± 0.251.10 ± 0.19
BdSP51.07.86 ± 1.480.04 ± 0.023.56 ± 0.10
Table 4. Primers used for qRT-PCR.
Table 4. Primers used for qRT-PCR.
GenesGenBank accession numberPrimer name and sequence (5′-3′)PCR efficiency (%)R2
BdSP1GAAP01000017S: ACACACTCGGGTTTTAGCGT
A: GAGGCGCAATCTTCACGTTG
106.250.981
BdSP2GAAP01000019S: AGAGGACGTATTGTTGGCGG
A: CGGCACGGATTCTCAAAACC
96.170.993
BdSP3GAAP01000020S: ATCCCACGGGTCGTGTAGTA
A: GTTGTGCCTGTTTCGACCAC
108.430.986
BdSP4GAAP01000021S: CAACGTGAAGATTGCGCCTC
A: CCTTTTTGGCACAGCCTTCC
91.950.989
BdSP5GAAP01000022S: TTCCCCCATCAGGTCTCAC
A: CCGAGTGTGCGTTGGATAC
94.610.988

Share and Cite

MDPI and ACS Style

Hou, M.-Z.; Shen, G.-M.; Wei, D.; Li, Y.-L.; Dou, W.; Wang, J.-J. Characterization of Bactrocera dorsalis Serine Proteases and Evidence for Their Indirect Role in Insecticide Tolerance. Int. J. Mol. Sci. 2014, 15, 3272-3286. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15023272

AMA Style

Hou M-Z, Shen G-M, Wei D, Li Y-L, Dou W, Wang J-J. Characterization of Bactrocera dorsalis Serine Proteases and Evidence for Their Indirect Role in Insecticide Tolerance. International Journal of Molecular Sciences. 2014; 15(2):3272-3286. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15023272

Chicago/Turabian Style

Hou, Ming-Zhe, Guang-Mao Shen, Dong Wei, Ya-Li Li, Wei Dou, and Jin-Jun Wang. 2014. "Characterization of Bactrocera dorsalis Serine Proteases and Evidence for Their Indirect Role in Insecticide Tolerance" International Journal of Molecular Sciences 15, no. 2: 3272-3286. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms15023272

Article Metrics

Back to TopTop