Next Article in Journal
Venom in Furs: Facial Masks as Aposematic Signals in a Venomous Mammal
Next Article in Special Issue
Occurrence and Ear Damage of Helicoverpa zea on Transgenic Bacillus thuringiensis Maize in the Field in Texas, U.S. and Its Susceptibility to Vip3A Protein
Previous Article in Journal
Subclinical Lipopolysaccharide from Salmonella Enteritidis Induces Dysregulation of Bioactive Substances from Selected Brain Sections and Glands of Neuroendocrine Axes
Previous Article in Special Issue
The C-Terminal Domain of the Bacillus thuringiensis Cry4Ba Mosquito-Specific Toxin Serves as a Potential Membrane Anchor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Proteomics of Peritrophic Matrix Provides an Insight into its Role in Cry1Ac Resistance of Cotton Bollworm Helicoverpa armigera

1
Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
2
The State Key Laboratory for Biology of Plant Disease and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, West Yuanmingyuan Road, Beijing, 100193, China
*
Author to whom correspondence should be addressed.
Toxins 2019, 11(2), 92; https://doi.org/10.3390/toxins11020092
Submission received: 2 January 2019 / Revised: 26 January 2019 / Accepted: 29 January 2019 / Published: 2 February 2019
(This article belongs to the Special Issue Insecticidal Toxins from Bacillus thuringiensis)

Abstract

:
Crystalline (Cry) proteins from Bacillus thuringiensis (Bt) are widely used in sprays and transgenic crops to control insect pests, but the evolution of insect resistance threatens their long-term use. Different resistance mechanisms have been identified, but some have not been completely elucidated. Here, the transcriptome of the midgut and proteome of the peritrophic matrix (PM) were comparatively analyzed to identify potential mechanism of resistance to Cry1Ac in laboratory-selected strain XJ10 of Helicoverpa armigera. This strain had a 146-fold resistance to Cry1Ac protoxin and 45-fold resistance to Cry1Ac activated toxin compared with XJ strain. The mRNA and protein levels for several trypsin genes were downregulated in XJ10 compared to the susceptible strain XJ. Furthermore, 215 proteins of the PM were identified, and nearly all had corresponding mRNAs in the midgut. These results provide new insights that the PM may participate in Bt resistance.
Key Contribution: Trypsin identified in peritrophic matrix may associated with Bt resistance in Helicoverpa armigera.

Graphical Abstract

1. Introduction

The toxins produced by Bacillus thuringiensis (Bt) have been used widely in transgenic plants for pest control with little or no harm to people and most non-target organisms [1,2,3,4]. The cotton bollworm, Helicoverpa armigera, is a principal cotton pest and inflicts major losses worldwide [5,6,7]. Transgenic cotton lines that produce Bt toxin have been useful in suppressing this polyphagous pest [6], but they often develop resistance to the toxins after long-term use [8,9,10].
Understanding the mode of action and the mechanisms of resistance to Bt proteins can help to enhance and sustain their efficacy against pests. The mode of action of Bt toxicity is complex. Models of the Bt mode of action agree that protoxins (the full-length forms of Cry1Ac proteins) are digested by midgut proteases into activated toxins, which then bind to insect midgut receptors, forming lytic pores in the membrane and leading to cell breakdown [11,12]. During this activation, five hundred amino acids from the carboxyl terminus and 40 amino acids from the amino terminus were removed. The relative molecular mass of protoxins converted from 130 kDa to 55 ~ 65 kDa of activated toxins [12,13]. It has been reported that reduced conversion of protoxin can cause greater resistance to protoxins than activated toxins [14,15,16,17,18].
Blocking any of these steps will lead to resistance. To date, numerous Bt-resistant insect populations have been selected in the laboratory [14,19,20,21,22,23], and various resistance mechanisms have been identified, such as altered activation of midgut digestive proteases, toxin sequestration by glycolipid moieties or esterase, elevated immune response, and reduced binding of Cry toxins [2,11,24,25].
The peritrophic matrix (PM), the acellular, porous tubular lining of the arthropod gut, is secreted by most insects and important for digestion, serving as a semi-permeable barrier between the epithelial cells of gut and the food bolus and protecting the midgut epithelium from infection by pathogens, damage by toxins, and mechanical damage by rough food particles [26,27,28,29]. The PM is composed of chitin fibrils with associated proteoglycans and glycoproteins and is proposed to assist the digestion process and immobilization of digestive enzymes, allowing reuse of hydrolytic enzymes and efficient acquisition of nutrients. Furthermore, the PM may also be a valid target for insect control [30,31,32].
Our previous work with 10 laboratory-selected strains of H. armigera suggested that reduced protease activity is associated with resistance to Cry1Ac [33]. Here we evaluated a resistant strain XJ10 (derived from XJ5) and found that the resistance ratio of Cry1Ac protoxin is much higher than that of the activated Cry1Ac toxin. We hypothesized that this ratio may be associated with the gut digestive proteases that take part in the conversion of the protoxin. In a comparative analysis of transcriptomes (RNA-seq) of midguts and proteomes (iTRAQ) of peritrophic matrix between XJ strain and XJ10 strain, we aimed to identify proteins that may be associated with Bt resistance and in the PM of H. armigera.

2. Results

2.1. Resistance to Cry1Ac Protoxin and Activated Toxin

Bioassay results indicated that the laboratory-selected XJ10 strain of H. armigera was resistant to Cry1Ac protoxin and active toxin compared to its parental strain (XJ, Table 1). The resistant ratios were calculated as the concentration (ng Cry1Ac per g diet) of Cry1Ac killing 50% (LC50) larvae for XJ10 divided by the LC50 for XJ larvae which were 146.8 and 45.0 for protoxin and activated toxin, respectively (Table 1).

2.2. Transcriptomic Analysis of Midguts from XJ and XJ10 Larvae

The mRNA transcript levels for genes from the midguts of strain XJ were compared to those from XJ10 strain for differential transcription. The mapping data analysis revealed 10,495 (75.85%) and 10,617 (76.73%) coding genes in the reference genome, respectively. We used a fold-change ≥ 1.5 and P-value < 0.05 as the threshold to judge that differences in expression were significant. We found that nearly 7.9% of all detected genes (457 up- and 378 down-regulated) were differentially expressed between XJ10 and XJ strain (Fig. 1A). To better understand the functional categories that differed between XJ10 and XJ strain, we used Blast2GO to assign GO categories to the 835 DEGs. The distribution of the GO terms are shown in Figure 1B. Cellular process, metabolic process and single-organism process were the major categories annotated under biological process. Cell, cell part and organelle were the major categories annotated under cellular component. As for molecular function, the major categories were binding and catalytic activity. Furthermore, KEGG analysis showed that 15 pathways were substantially enriched (P < 0.05), including ribosome, metabolic pathways (Table S1).
Trypsin and chymotrypsin are very important digestive enzymes and play important roles in Bt protoxin activation. In this study, as shown in Table 2, 11 trypsin and chymotrypsin were found differentially expressed between XJ10 and XJ. Mechanism of action of Bt toxins is complex, blocking any step may lead to resistance. Besides trypsin, other differentially regulated genes, including ABCs, polycalin, APN and ALP, possibly linked to Bt resistance are also listed in Table 2. We chose Bt resistance-related genes for proteins such as trypsin, esterase, ABCs and Bt receptors for qRT-PCR analysis, which indicated that most of the Bt resistance-related genes had expression patterns similar to those shown by the RNA-seq data (Figure 2).

2.3. Proteomic Analysis of PM from XJ and XJ10 Larvae

In the parallel iTRAQ analysis to compare the proteome of PM from XJ and XJ10 (three biological replicates in each group), 215 proteins were identified and quantified, (see details on the proteins in Supplementary Table S2). Based on the specific functions of the proteins, we divided them into categories (Table 3) that included chitin-related, digestion-related, lipocalins-related, immune-related and so on. Polycalin and APNs, which have been reported as Cry1Ac receptors, were also found in this study. Similar to other findings, most of the identified proteins could not be characterized or had unknown functions [34,35].
In the GO analysis, the 215 identified proteins were enriched in 16 biological process terms. The top 10 processes (Figure 3A) indicate that most proteins are involved in biological processes related to metabolic processes. Catalytic activity and binding were the most abundant molecular function categories (Figure 3B). Additionally, among 12 other cellular component terms, cell part component had the largest group of proteins (Figure 3C). The results of the GO enrichment revealed that the proteins of PM were predominantly binding proteins and have catalytic activity, located in membrane and involved in metabolic processes.
In the analysis of global changes of PM proteins between strain XJ10 and XJ, 12 proteins were classified as differentially expressed proteins (DEPs) (Table 4). Among these DEPs, most downregulated proteins were active digestive hydrolases including five trypsins, two chymotrypsins, carboxypeptidase A and one uncharacterized protein. In XJ10 strain, α-amylase, which is involved in food digestion, was upregulated. Chitin deacetylase and unconventional myosin were also upregulated in XJ10 strain.

2.4. Correlation between Transcriptome and Proteome

Because PM proteins are secreted from the midgut, we analyzed the correlation between the identified proteins of PM and mRNAs of the midgut. The distributions of the ratio of the corresponding protein to the mRNA are shown in Figure 4. Of the 215 identified proteins, 95.8% (205/215) had corresponding mRNAs, indicating that nearly all PM proteins have a corresponding mRNA in the midgut and confirmed that PM proteins are secreted from the midgut. As for the 12 DEPs, their corresponding mRNA expression levels are given in Table 4. For the downregulated DEPs, including four trypsins (XM_021337877.1, XM_021337869.1, XM_021340592.1, and XM_021340602.1) and one uncharacterized protein (XM_021340037.1), the corresponding transcript levels were also down-regulated.

3. Discussion

The laboratory-selected Cry1Ac-resistant H. armigera strain XJ10 had a resistance ratio of 45 for the active Cry1Ac toxin and 146 for Cry1Ac protoxin. Because the resistance ratio for the Cry1Ac protoxin was much higher than for the activated toxin, these results suggest that reduced activation of Cry1Ac may contribute to resistance in XJ10 strain. However, the resistance ratio for the activated toxin indicated the presence of another resistance mechanism in XJ10. In the present study of the transcriptome of the midgut and proteome of the PM between XJ and XJ10 strains using RNA-seq, hundreds of genes were differently expressed, and 12 of the 215 proteins that were identified by iTRAQ varied in abundance between XJ and XJ10 strain, and nearly all the 215 identified proteins in the PM had corresponding mRNAs in the midgut.

3.1. Identification of Potential Resistance Mechanisms in XJ10

The mechanisms of Bt resistance include diverse process, because it could occur by blocking of any of the steps noted above. Among the differentially regulated genes for proteins potentially related to Bt resistance (Table 2 and Table 4), gut proteases play an important role in protoxin activation and in the Bt-resistance mechanism in several lepidopteran species. For example, downregulation of protease gene confers Cry1Ac resistance to H. armigera stain Akola-R from India and strain LF5 from China [7,14]. In contrast, increased activity of gut proteases in Spodoptera littoralis enhances the resistance to Cry1C, possibly due to overdegradation of the toxin [36]. In the present study, many trypsin genes were differentially expressed in transcriptomic profiles and also in proteome profiles of the PM. The PM is proposed to assist the digestion process by partitioning the gut lumen into ectoperitrophic space (between PM and lumen) and endoperitrophic space (between epithelium and PM) and immobilization of digestive enzymes [36]. These immobilized digestive enzymes may play important roles in the activation of Cry1Ac protoxin. Protease-mediated resistance was first demonstrated in a laboratory-selected strain of Plodia interpunctella, for which larval survival was genetically linked with a lack of major trypsin-like gut proteases after Cry1Ac treatment [16]. In other insects, such as Heliothis virescens and Ostrinia nubilalis, reduced protease activity is associated with resistance to Cry1A toxins [37,38].
Sequestration of the toxin by esterases has also been reported to be associated with resistance. Esterase could be responsible for sequestering large quantities of Cry1Ac in the 275-fold resistant H. armigera strain [39]. Here we found that the expression of esterase was upregulated in strain XJ10. Cry toxins are known to have different receptors such as cadherin (Cad), aminopeptidase-N (APN), alkaline phosphatase (ALP) and ABC transporters, which bind to Cry1Ac toxin through a combination of different modes and are important in the receptor-mediated toxicity of Cry toxins. In the case of XJ10, ALP2 and APN1 were upregulated in mRNA level compared with XJ strain. Many studies showed that the mutation of APN1 or downregulation of APN1 or ALP2 is correlated with Bt resistance, but no studies had reported that the upregulated APN1 or ALP2 in mRNA level was related to Bt resistance [40,41,42,43]. APN and ALP are attached to the cell membrane via a glycosyl phosphotidylinositol (GPI) anchor, and the N- and/or O-linked glycans are thought to mediate the binding of Cry toxins. Whether the upregulated expression of ALP2 and APN1 is associated with resistance or a post-translational modification needs further study.

3.2. Identification of PM Proteins

In our comprehensive analysis of cotton bollworm larval PM proteins using iTRAQ, we identified 215 proteins, more than previously reported [29], among several major classes of proteins (Table 3), and thus potentially different functions for the PM of the midgut. In this study, the two mucins and the chitin deacetylase (CDA) protein represented two classes of chitin-binding proteins. CDA is a hydrolytic enzyme that degrades the glycosidic bonds of chitin and may control the rigidity and porosity of chitin-containing PM [44,45]. In XJ10, the protein level of CDA was higher than in XJ. Thus, chitin rigidity in XJ10 may be higher and obstruct Cry1Ac toxin entry into the PM.
Forty-six of the identified proteins were digestive hydrolases of various types, including trypsin, carboxypeptidase, α-amylase, serine protease and lipase. Since the midgut is for digestion and absorption, these digestive enzymes might be immobilized on the PM and present in two forms—either bound to the PM or soluble in the gut lumen. The PM may accelerate digestion in cotton bollworm larvae via PM-bound digestive enzymes. In strain XJ10, the level of α-amylase 2-like was higher than in XJ. α-Amylase is widely distributed among animals, plants and microbes and catalyzes the hydrolysis of starch [46], and may have no association with resistance. Of the 12 proteins we found to be differentially expressed, seven were trypsin, which were all downregulated in XJ10.
Decreased protoxin activation of Cry1Ac or Cry1Ab has been associated with resistance in many insects. Our finding that the resistance ratio of XJ10 for Cry1Ac protoxin was much higher than the resistance ratio of the activated Cry1Ac toxin is similar to previous work in a laboratory-selected strain LF5 of H. armigera, for which the resistance ratio was 2.8 times higher for the Cry1Ac protoxin than the activated toxin. The resistance of LF5 is genetically linked with a trypsin gene HaTryR and reducing the expression of this gene using RNAi increases the survival of susceptible larvae treated with Cry1Ac protoxin [14], implying that trypsin protease is important for activating the Cry1Ac protoxin. We thus suggest that these downregulated trypsins found in the PM may be responsible for the reduced Cry1Ac activation in XJ10.
Serine proteases are important proteolytic enzymes for digestion and insect innate immune systems [47,48]. Although we identified seven serine proteases and two serine protease inhibitors in this study, their protein expression levels did not change. We also identified some heat shock proteins and numerous proteins in the PM with unknown functions. These proteins should not be neglected, because they may play important roles in the molecular architecture of PM.
In our exploration of the gene expression changes in the midgut and protein changes in the PM between XJ10 and XJ strain, both the mRNA and protein levels for several trypsins were downregulated, revealing a potential association with Cry1Ac protoxin activation and the high resistance to Cry1Ac protoxin. However, reduced activation cannot explain the 45-fold resistance ratio for the activated Cry1Ac toxin, indicating that another mechanism(s) also contributes to resistance in XJ10. Furthermore, the identification of proteins of PM enabled us to better understand the nature of the PM and its involvement in digestion and deactivation of toxins such as Bt toxins.

4. Materials and Methods

4.1. Insect Strains

The XJ strain was collected in a cotton field in Xiajin (Shandong Province, China) in 2004 and reared in the laboratory without exposure to Bt toxins or insecticides [33]. The resistant XJ10 strain was derived from strain XJ via selection on Cry1Ac-contaminated artificial diet through a series of progressively more resistant strains, XJ1, XJ5 and XJ10, with the number of each resistant strain corresponding to the concentration of Cry1Ac. Thus, XJ10 was selected on 10 μg Cry1Ac protoxin per ml of diet. During the selection, neonates were reared on the artificial diet with the respective dose of toxin for 7 days, then well-developed larvae were transferred to non-Bt artificial diet until moth emergence. Insects were reared in an insect chamber with a controlled environment (27 ± 2 °C, 75 ± 10% RH, 14L: 10D).

4.2. Insect Bioassays

Bioassays on artificial diets with various concentrations of Cry1Ac protoxin or active toxin in 24-well were performed according to the methods of Liu et al. [14]: a 4-day-old larva was placed in each well, with 24 larvae for each treatment. Mortality was recorded after 7 d. The concentration that killed 50% of larvae (LC50) was determined using a probit analysis.

4.3. RNA Extraction and cDNA Library Construction

Fifth instar larvae form XJ and XJ10 strains were anesthetized on ice and dissected longitudinally to obtain midguts (n=24, 3 replicates). Total RNA were isolated using Trizol (Invitrogen, USA) according to the manufacturer’s protocol. 1% agarose gels were used to check RNA degradation and contamination and Nanophotometer1 spectrophotometer (IMPLEN, CA, USA) was used to check the purity of RNA. Library construction and sequencing using Illumina HiSeqTM 2500 were done by Gene denovo Biotechnology Co., Ltd., Guangzhou, China. Briefly, mRNA was purified from total RNA using poly T oligo-attached magnetic beads.

4.4. Data Processing and Analysis

Clean reads were obtained by removing reads containing adapters, poly N, and low-quality reads form the raw reads. The clean data were then mapped to the reference genome (Accession NO., PRJNA388211) by TopHat 2 (v 2.0.3.12) [49] using the following modification from default parameters: the distance between mate-pair reads, 50bp; the error of distance between mate-pair reads, ± 80bp; up to two mismatches allowed. The expression levels of genes were estimated using RSEM [50] and normalized using the FPKM (fragments per kb of transcript per million mapped reads). Differential expression analysis was performed using the RStudio with package edgeR. The resulting P values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate. The genes were considered differentially expressed if they had an adjusted P-value <0.05 and a change of at least 1.5-fold. Gene Ontology (GO) analysis was conducted for the functional classification of differentially expressed genes (DEGs), and pathway analysis was carried out using Kyoto Encyclopedia of Genes and Genomes (KEGG).

4.5. qRT-PCR Analysis of Gene Expression Levels

Total RNA was extracted from the midgut as described above with three replications of samples, then 1 μg of RNA was reverse transcribed for first-strand cDNA synthesis using TransScript One-step gDNA Removal and cDNA Synthesis SuperMix (Transgen Biotech, China). Primers (Table S3) were designed with Beacon Designer 7.9 software (PREMIER Biosoft International, CA, USA). The reaction mixture contained 10 μL TransStart Tip Green qRCR SuperMix (Transgen Biotech, China), 1.0 μL cDNA, 8.2 μL ddH2O, and 0.4 μL each of the forward and reverse gene-specific primers and run in a Stepone Plus Real-time system with thermal cycling set at 95 °C for 3 min; 40 cycles of 95 °C for 10 s and 55 °C for 60 s. The qRT-PCR for each sample was conducted with three technical replicates and three biological replicates. Actin was used as a reference gene to normalize the content of cDNA. Relative expression levels for target genes, in relation to the most reliable reference gene, were calculated via the 2−∆∆CT method. Levels were analyzed for significant differences using Student’s t-test (SPSS).

4.6. iTRAQ-Based Proteome Determination and Data Analysis

The peritrophic matrix (PM) was excised from 5th instars of strains XJ and XJ10 that had been actively feeding, as described previously [29]. Total protein extracted using Tissue Protein Extraction Kit (Cwbiotech Co. Ltd., Beijing, China) following the manufacturer’s method. For digestion, 100 μg of total protein, quantified by the Bradford method, was incubated overnight with Trypsin Gold (Promega; protein: trypsin = 30:1) at 30 °C for 16 h. After trypsin digestion, the resultant peptides were dried by vacuum centrifugation. iTRAQ labeling was performed using an iTRAQ reagent 8PLEX Multiplex kit (Applied Biosystems, US) according to the manufacturer’s protocol. Samples from XJ and XJ10 were labeled with iTRAQ tags 115 and 119, respectively. Then, the pooled mixtures of iTRAQ-labeled peptides were fractionated by strong cationic exchange (SCX) chromatography using a Shimadzu LC-20AB HPLC Pump system. Collected fractions were pooled into 10 final fractions for nano LC-MS/MS analysis after desalting by Strata XC18 column (Phenomenex) and vacuum dried. The MS/MS analysis was performed using a LTQ-Qrbitrap velos mass spectrometer (ThermoFisher Scientific, Rockford, IL, USA). MS data were acquired using a data-dependent top 10 method, dynamically choosing the most abundant precursor ions from the survey scan (300–1800 m/z) for HCD fragmentation. Determination of the target value based on predictive automatic gain control. The raw data files were converted into MGF format using Proteome Discover 1.2 (Thermo Fisher, Waltham, MA, USA). Proteins were identified and quantified using the MASCOT engine (Matrix Science, Boston, MA, USA) and the reference genome (accession PRJNA388211). Proteins with a 1.2-fold change and P-value ≤ 0.5 between two samples were considered to be significant differentially expressed proteins.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6651/11/2/92/s1, Table S1: KEGG pathways of P-value < 0.05 in XJ10 compared with XJ, Table S2: Proteins were identified and quantified in Peritrophic matrix, Table S3: Primers used for qRT-PCR.

Author Contributions

K.W. and Y.X. designed the present study. M.J. performed all the experiments. K.W., M.J., Y.X., S.C. and C.L. analyzed data and wrote the manuscript. All authors have read and approved the manuscript for publication.

Funding

This work was supported by Key Project for Breeding Genetic Modified Organisms (2016ZX08012004-003); National Natural Science Foundation of China (No. 31601646); Dapeng New District Industry Development Funds (KY20160103); and ShenZhen Science and Technology Project (JCYJ20170303154440838).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mendelsohn, M.; Kough, J.; Vaituzis, Z.; Matthews, K. Are bt crops safe? Nat. Biotechnol. 2003, 21, 1003–1009. [Google Scholar] [CrossRef] [PubMed]
  2. Liliana, P.L.; Mario, S.; Alejandra, B. Bacillus thuringiensis insecticidal three-domain cry toxins: Mode of action, insect resistance and consequences for crop protection. Fems Microbiol. Rev. 2013, 37, 3–22. [Google Scholar]
  3. Van Frankenhuyzen, K. Insecticidal activity of Bacillus thuringiensis crystal proteins. J. Invertebr. Pathol. 2009, 101, 1–16. [Google Scholar] [CrossRef] [PubMed]
  4. Hilbeck, A.; Otto, M. Specificity and combinatorial effects of Bacillus thuringiensis cry toxins in the context of GMO environmental risk assessment. Front. Environ. Sci. 2015, 3, 71. [Google Scholar] [CrossRef]
  5. Wu, K.M.; Guo, Y.Y. The evolution of cotton pest management practices in china. Annu. Rev. Entomol. 2005, 50, 31–52. [Google Scholar] [CrossRef] [PubMed]
  6. Kong-Ming, W.; Yan-Hui, L.; Hong-Qiang, F.; Yu-Ying, J.; Jian-Zhou, Z. Suppression of cotton bollworm in multiple crops in china in areas with Bt toxin-containing cotton. China Basic Sci. 2008, 321, 1676–1678. [Google Scholar]
  7. Raman, R.; Naresh, A.; Swaminathan, S.; Rao, N.G.V.; Nimbalkar, S.A.; Bhatnagar, R.K. Resistance of Helicoverpa armigera to Cry1Ac toxin from Bacillus thuringiensis is due to improper processing of the protoxin. Biochem. J. 2009, 419, 309–316. [Google Scholar]
  8. Tabashnik, B.E. Evolution of resistance to Bacillus thuringiensis. Annu. Rev. Entomol. 1994, 39, 47–79. [Google Scholar] [CrossRef]
  9. Tabashnik, B.E.; Huang, F.; Ghimire, M.N.; Leonard, B.R.; Siegfried, B.D.; Rangasamy, M.; Yang, Y.; Wu, Y.; Gahan, L.J.; Heckel, D.G.; et al. Efficacy of genetically modified Bt toxins against insects with different genetic mechanisms of resistance. Nat. Biotechnol. 2011, 29, 1128–1131. [Google Scholar] [CrossRef] [Green Version]
  10. Tabashnik, B.E.; Brévault, T.; Carrière, Y. Insect resistance to Bt crops: Lessons from the first billion acres. Nat. Biotechnol. 2013, 31, 510–521. [Google Scholar] [CrossRef]
  11. Bravo, A.; Likitvivatanavong, S.; Gill, S.S.; Soberon, M. Bacillus thuringiensis: A story of a successful bioinsecticide. Insect Biochem. Mol. 2011, 41, 423–431. [Google Scholar] [CrossRef] [PubMed]
  12. Gill, S.S.; Cowles, E.A.; Pietrantonio, P.V. The mode of action of Bacillus thuringiensis endotoxins. Annu. Rev. Entomol. 1992, 37, 615–636. [Google Scholar] [CrossRef] [PubMed]
  13. Adang, M.J.; Crickmore, N.; Jurat-Fuentes, J.L. Chapter two–diversity of Bacillus thuringiensis crystal toxins and mechanism of action. Adv. Insect Physiol. 2014, 47, 39–87. [Google Scholar]
  14. Liu, C.; Xiao, Y.; Li, X.; Oppert, B.; Tabashnik, B.E.; Wu, K. Cis-mediated down-regulation of a trypsin gene associated with Bt resistance in cotton bollworm. Sci. Rep. 2014, 4, 7219. [Google Scholar] [CrossRef] [PubMed]
  15. Tabashnik, B.E.; Zhang, M.; Fabrick, J.A.; Wu, Y.; Gao, M.; Huang, F.; Wei, J.; Zhang, J.; Yelich, A.; Unnithan, G.C. Dual mode of action of Bt proteins: Protoxin efficacy against resistant insects. Sci. Rep. 2015, 5, 15107. [Google Scholar] [CrossRef]
  16. Oppert, B.; Kramer, K.J.; Beeman, R.W.; Johnson, D.; Mcgaughey, W.H. Proteinase-mediated insect resistance to Bacillus thuringiensis toxins. J. Biol. Chem. 1997, 272, 23473–23476. [Google Scholar] [CrossRef]
  17. Wu, Y. Detection and mechanisms of resistance evolved in insects to cry toxins from Bacillus thuringiensis. Adv. Insect Physiol. 2014, 47, 297–342. [Google Scholar]
  18. Wei, J.; Liang, G.; Wang, B.; Zhong, F.; Chen, L.; Khaing, M.M.; Zhang, J.; Guo, Y.; Wu, K.; Tabashnik, B.E. Activation of Bt protoxin Cry1Ac in resistant and susceptible cotton bollworm. PLoS ONE 2016, 11, e0156560. [Google Scholar] [CrossRef]
  19. Luo, S.; Wang, G.; Liang, G.; Wu, K.M.; Bai, L.; Ren, X.; Guo, Y. Binding of three Cry1a toxins in resistant and susceptible strains of cotton bollworm (Helicoverpa armigera). Pestic. Biochem. Phys. 2006, 85, 104–109. [Google Scholar] [CrossRef]
  20. Kaur, P.; Dilawari, V.K. Inheritance of resistance to Bacillus thuringiensis Cry1Ac toxin in Helicoverpa armigera (hübner) (lepidoptera: Noctuidae) from India. Pest Manag. Sci. 2011, 67, 1294–1302. [Google Scholar] [CrossRef]
  21. Xinjun, X.; Liangying, Y.; Yidong, W. Disruption of a cadherin gene associated with resistance to Cry1Ac {delta}-endotoxin of bacillus thuringiensis in Helicoverpa armigera. Appl. Environ. Microbiol. 2005, 71, 948–954. [Google Scholar]
  22. Josephakhurst, R.; Janebird, L.; Beard, C. Resistance to the Cry1Ac δ-endotoxin of Bacillus thuringiensis in the cotton bollworm, Helicoverpa armigera (lepidoptera: Noctuidae). J. Econ. Entomol. 2003, 96, 1290–1299. [Google Scholar]
  23. Xiao, Y.; Zhang, T.; Liu, C.; Heckel, D.G.; Li, X.; Tabashnik, B.E.; Wu, K. Mis-splicing of the abcc2 gene linked with bt toxin resistance in Helicoverpa armigera. Sci. Rep. 2014, 4, 6184. [Google Scholar] [CrossRef] [PubMed]
  24. Schnepf, E.; Crickmore, N.; Van Rie, J.; Lereclus, D.; Baum, J.; Feitelson, J.; Zeigler, D.R.; Dean, D.H. Bacillus thuringiensis and its pesticidal crystal proteins. Microbiol. Mol. Biol. Rev. 1998, 62, 775–806. [Google Scholar] [PubMed]
  25. Heckel, D.G.; Gahan, L.J.; Baxter, S.W.; Zhao, J.Z.; Shelton, A.M.; Gould, F.; Tabashnik, B.E. The diversity of bt resistance genes in species of lepidoptera. J. Invertebr. Pathol. 2007, 95, 192–197. [Google Scholar] [CrossRef] [PubMed]
  26. Terra, W.R. The origin and functions of the insect peritrophic membrane and peritrophic gel. Arch. Insect Biochem. 2010, 47, 47–61. [Google Scholar] [CrossRef] [PubMed]
  27. Hegedus, D.; Erlandson, M.; Gillott, C.; Toprak, U. New insights into peritrophic matrix synthesis, architecture, and function. Annu. Rev. Entomol. 2009, 54, 285. [Google Scholar] [CrossRef]
  28. Toprak, U.; Hegedus, D.D.; Baldwin, D.; Coutu, C.; Erlandson, M. Spatial and temporal synthesis of Mamestra configurata peritrophic matrix through a larval stadium. Insect Biochem. Mol. Biol. 2014, 54, 89–97. [Google Scholar] [CrossRef]
  29. Campbell, P.M.; Cao, A.T.; Hines, E.R.; East, P.D.; Gordon, K.H. Proteomic analysis of the peritrophic matrix from the gut of the caterpillar, Helicoverpa armigera. Insect Biochem. Mol. Biol. 2008, 38, 950–958. [Google Scholar] [CrossRef]
  30. Wang, P.; Granados, R.R. Molecular structure of the peritrophic membrane (pm): Identification of potential PM target sites for insect control. Arch. Insect Biochem. 2010, 47, 110–118. [Google Scholar] [CrossRef]
  31. Wang, P.; Granados, R.R. Calcofluor disrupts the midgut defense system in insects. Insect Biochem. Mol. Biol. 2000, 30, 135–143. [Google Scholar] [CrossRef]
  32. Tibor, P.; Allen, C.; Paul, W.; Dawn, S.L. Insect feeding mobilizes a unique plant defense protease that disrupts the peritrophic matrix of caterpillars. Proc. Natl. Acad. Sci. USA 2002, 99, 13319–13323. [Google Scholar] [Green Version]
  33. Cao, G.; Zhang, L.; Liang, G.; Li, X.; Wu, K. Involvement of nonbinding site proteinases in the development of resistance of Helicoverpa armigera (Lepidoptera: Noctuidae) to Cry1Ac. J. Econ. Entomol. 2013, 106, 2514–2521. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, Y.; Xiu, J.F.; Cheng, J.Z.; Luo, M.; Zhao, P.; Shang, X.L.; Wang, T.; Jian Wei, W.U. Proteomic analysis of the peritrophic matrix from the midgut of third instar larvae, Musca domestica. Biomed. Environ. Sci. 2016, 29, 56–65. [Google Scholar] [PubMed]
  35. Wang, L.; Li, F.; Wang, B.; Xiang, J. Structure and partial protein profiles of the peritrophic membrane (PM) from the gut of the shrimp Litopenaeus vannamei. Fish Shellfish Immunol. 2012, 33, 1285–1291. [Google Scholar] [CrossRef] [PubMed]
  36. Keller, M.; Sneh, B.; Strizhov, N.; Prudovsky, E.; Regev, A.; Koncz, C.; Schell, J.; Zilberstein, A. Digestion of delta-endotoxin by gut proteases may explain reduced sensitivity of advanced instar larvae of Spodoptera littoralis to CryIc. Insect Biochem. Mol. Biol. 1996, 26, 365–373. [Google Scholar] [CrossRef]
  37. Karumbaiah, L.; Oppert, B.; Jurat-Fuentes, J.L.; Adang, M.J. Analysis of midgut proteinases from Bacillus thuringiensis -susceptible and -resistant Heliothis virescens (lepidoptera: Noctuidae). Comp. Biochem. Phys. B 2007, 146, 139–146. [Google Scholar] [CrossRef]
  38. Li, H.; Oppert, B.; Higgins, R.A.; Huang, F.; Zhu, K.Y.; Buschman, L.L. Comparative analysis of proteinase activities of Bacillus thuringiensis-resistant and -susceptible Ostrinia nubilalis (Lepidoptera:Crambidae). Insect Biochem. Mol. Biol. 2004, 34, 753–762. [Google Scholar] [CrossRef]
  39. Gunning, R.V.; Dang, H.T.; Kemp, F.C.; Nicholson, I.C.; Moores, G.D. New resistance mechanism in Helicoverpa armigera threatens transgenic crops expressing Bacillus thuringiensis Cry1Ac toxin. Appl. Environ. Microbiol. 2005, 71, 2558–2563. [Google Scholar] [CrossRef]
  40. Juan Luis, J.F.; Gould, F.L.; Adang, M.J. Altered glycosylation of 63- and 68-kilodalton microvillar proteins in Heliothis virescens correlates with reduced Cry1 toxin binding, decreased pore formation, and increased resistance to Bacillus thuringiensis Cry1 toxins. Appl. Environ. Microbiol. 2002, 68, 5711–5717. [Google Scholar]
  41. Zhang, S.; Cheng, H.; Gao, Y.; Wang, G.; Liang, G.; Wu, K. Mutation of an aminopeptidase N gene is associated with Helicoverpa armigera resistance to Bacillus thuringiensis Cry1Ac toxin. Insect Biochem. Mol. Biol. 2009, 39, 421–429. [Google Scholar] [CrossRef]
  42. Herrero, S.; Gechev, T.; Bakker, P.L.; Moar, W.J.; Maagd, R.A.D. Bacillus thuringiensis Cry1Ca-resistant Spodoptera exigua lacks expression of one of four aminopeptidase N genes. BMC Genom. 2005, 6, 96. [Google Scholar] [CrossRef] [PubMed]
  43. Tiewsiri, K.; Wang, P. Differential alteration of two aminopeptidases N associated with resistance to Bacillus thuringiensis toxin Cry1Ac in cabbage looper. Proc. Natl. Acad. Sci. USA 2011, 108, 14037–14042. [Google Scholar] [CrossRef] [PubMed]
  44. Luschnig, S.; Bätz, T.; Armbruster, K.; Krasnow, M.A. Serpentine and vermiform encode matrix proteins with chitin binding and deacetylation domains that limit tracheal tube length in Drosophila. Curr. Biol. 2006, 16, 186–194. [Google Scholar] [CrossRef] [PubMed]
  45. Shen, Z.; Jacobs-Lorena, M. Characterization of a novel gut-specific chitinase gene from the human malaria vector Anopheles gambiae. J. Biol. Chem. 1997, 272, 28895–28900. [Google Scholar] [CrossRef] [PubMed]
  46. Ewelina, C.; Wojciech, B.A.; Monika, B.; Wlodzimierz, G. Cloning, expression, and purification of insect (Sitophilus oryzae) alpha-amylase, able to digest granular starch, in Yarrowia lipolytica host. Appl. Microbiol. Biotechnol. 2015, 99, 2727–2739. [Google Scholar]
  47. Gorman, M.J.; Paskewitz, S.M. Serine proteases as mediators of mosquito immune responses. Insect Biochem. Mol. Biol. 2001, 31, 257–262. [Google Scholar] [CrossRef]
  48. Wolfson, J.L.; Murdock, L.L. Diversity in digestive proteinase activity among insects. J. Chem. Ecol. 1990, 16, 1089–1102. [Google Scholar] [CrossRef]
  49. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. Tophat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14. [Google Scholar] [CrossRef]
  50. Li, B. Rsem: Accurate transcript quantification from rna-seq data with or without a reference genome. BMC Bioinformatics 2011, 12, 323. [Google Scholar] [CrossRef]
Figure 1. Transcriptomic difference between strain XJ10 and XJ. (A) Fold-changes on the x-axis represent the ratio of transcript abundance between strain XJ10 and XJ. Differentially expressed transcripts are highlighted in green (down-regulated) and red (up-regulated), respectively on the Volcano plot. (B) Gene Ontology (GO) classification of genes differentially expressed between XJ10 and XJ strain and grouped into hierarchically structured GO terms biological process, cellular component, and molecular function.
Figure 1. Transcriptomic difference between strain XJ10 and XJ. (A) Fold-changes on the x-axis represent the ratio of transcript abundance between strain XJ10 and XJ. Differentially expressed transcripts are highlighted in green (down-regulated) and red (up-regulated), respectively on the Volcano plot. (B) Gene Ontology (GO) classification of genes differentially expressed between XJ10 and XJ strain and grouped into hierarchically structured GO terms biological process, cellular component, and molecular function.
Toxins 11 00092 g001
Figure 2. qRT-PCR validation of eight Bt-resistance-related genes differentially expressed between XJ10 and XJ strain. The mRNA levels were compared using Student’s t-test (* P < 0.05, ** P < 0.01).
Figure 2. qRT-PCR validation of eight Bt-resistance-related genes differentially expressed between XJ10 and XJ strain. The mRNA levels were compared using Student’s t-test (* P < 0.05, ** P < 0.01).
Toxins 11 00092 g002
Figure 3. GO analysis of the functional categories of proteins identified from PM. (A) Distribution of enriched biological processes (BP). (B) Distribution of cell component (CC) enrichment. (C) Distribution of molecular function (MF) enrichment.
Figure 3. GO analysis of the functional categories of proteins identified from PM. (A) Distribution of enriched biological processes (BP). (B) Distribution of cell component (CC) enrichment. (C) Distribution of molecular function (MF) enrichment.
Toxins 11 00092 g003
Figure 4. Changes in mRNA and cognate protein abundance between midgut and PM. Relative change in abundance (XJ10/XJ) is shown on a log2 scale. Each letter denotes ratio of abundance of mRNA to protein: e, no significant change in both mRNA and protein; c and g, the expression of mRNA and protein with the same trend; and a, b, d, f, h and i, the expression of mRNA and protein with opposite trends.
Figure 4. Changes in mRNA and cognate protein abundance between midgut and PM. Relative change in abundance (XJ10/XJ) is shown on a log2 scale. Each letter denotes ratio of abundance of mRNA to protein: e, no significant change in both mRNA and protein; c and g, the expression of mRNA and protein with the same trend; and a, b, d, f, h and i, the expression of mRNA and protein with opposite trends.
Toxins 11 00092 g004
Table 1. Effects of Cry1Ac protoxin and activated toxin on mortality of H. armigera larvae of resistant strain XJ10 and susceptible strain XJ.
Table 1. Effects of Cry1Ac protoxin and activated toxin on mortality of H. armigera larvae of resistant strain XJ10 and susceptible strain XJ.
StrainForm of Cry1AcLC50 a (95% Confidence Limits)Resistance Ratio b
XJprotoxin8.41 (6.31–11.25)1
XJ10protoxin1233.91 (923.58–1665.62)146.79
XJactivated toxin7.16 (5.25–9.76)1
XJ10activated toxin322.48 (248.42–417.66)45.05
a Concentration killing 50% with 95% fiducial limits in ng Cry1Ac per g diet. b Resistance ratio, LC50 for XJ10 divided by LC50 for XJ.
Table 2. Partial list of Bt-resistance-associated genes that are differentially expressed (DEGs) between XJ10 and XJ strain.
Table 2. Partial list of Bt-resistance-associated genes that are differentially expressed (DEGs) between XJ10 and XJ strain.
Gene IDFC aPadj ValueDescriptionMechanisms
XM_021337885.12.531.44 × 10−6trypsinAltered activation of Cry toxins
XM_021345079.11.543.60 × 10−17trypsin-like protease
XM_021337887.1−2.697.99 × 10−10trypsin
XM_021340602.1−3.565.44 × 10−9trypsin T4
XM_021340597.1−4.293.58 × 10−26trypsin
XM_021342117.1−4.596.80 × 10−11trypsin-like proteinase T2α
XM_021340588.1−8.890.000297197trypsin
XM_021345839.12.416.13 × 10−24chymotrypsin
XM_021345877.11.961.30 × 10−11chymotrypsin-like protease
XM_021334546.1−5.273.40 × 10−24chymotrypsin
XM_021345783.1−5.533.39 × 10−67chymotrypsin-like protease
XM_021337546.12.031.97 × 10−5esterase E4-likeSequestering the toxin
XM_021334360.1−1.705.59 × 10−10ABC transporter G family 23ABCs
XM_021345614.12.028.88 × 10−19ABCC1 protein
XM_021334935.1−1.811.17 × 10−11PolycalinBinding proteins
XM_021337081.11.661.17 × 10−10aminopeptidase N1
XM_021339316.12.246.29 × 10−11alkaline phosphatase 2
a Fold change of DEGs, positive value indicates up-regulation while negative value denotes down-regulation.
Table 3. Proteins identified from the peritrophic matrix of H. armigera.
Table 3. Proteins identified from the peritrophic matrix of H. armigera.
NameNumber of Distinct PeptidesMS/MS NumberSequence Coverage (%)Predicted MW (kDa)Accession Number
Chitin associated
mucin 17-like160.3639.67XM_021332791.1
chitin deacetylase 5a3127.445.297XM_021341184.1
insect intestinal mucin 26327.3121.82XM_021326099.1
Active hydrolases
chymotrypsin-like protease C95821.138.71XM_021335717.1
trypsin, alkaline C-like113.927.167XM_021337856.1
trypsin144.327.749XM_021337869.1
trypsin-like protease2211.926.816XM_021337877.1
trypsin-like protease237.526.86XM_021337879.1
trypsin116.526.217XM_021338513.1
trypsin-7-like102.432.302XM_021340310.1
trypsin T2a3814.327.81XM_021340589.1
trypsin142.727.549XM_021340592.1
trypsin-like protease51323.927.503XM_021340596.1
trypsin-like protease714244.926.916XM_021340599.1
trypsin T4261126.772XM_021340602.1
trypsin 2123.332.037XM_021344969.1
chymotrypsin3314.630.518XM_021345778.1
chymotrypsinogen6221.430.834XM_021345781.1
chymotrypsin-like protease C8146.230.106XM_021345791.1
chymotrypsin-like protease2613.732.376XM_021345818.1
chymotrypsin-like protease128.629.92XM_021345819.1
carboxypeptidase
carboxypeptidase B-like234.448.983XM_021328067.1
carboxypeptidase precursor111.947.903XM_021330765.1
carboxypeptidase B precursor235.148.318XM_021330831.1
carboxypeptidase A122.848.526XM_021330834.1
carboxypeptidase111.642.284XM_021330844.1
carboxypeptidase3510.347.935XM_021330848.1
aminopeptidase N5145114.37XM_021337080.1
aminopeptidase N484.2112.81XM_021337081.1
alpha-amylase11256.009XM_021332568.1
Inactive hydrolases
serine protease inhibitor 5122.845.123XM_021327264.1
serine protease 24226.643.33XM_021330937.1
serine protease inhibitor111.688.323XM_021334028.1
diverged serine protease41313.327.289XM_021338518.1
transmembrane protease serine 9-like2182.579.868XM_021338520.1
serine protease63745.726.952XM_021340600.1
serine protease inhibitor 3 isoform X1111.5106.56XM_021343251.1
diverged serine protease51322.125.554XM_021344422.1
serine protease 52134.727.001XM_021344467.1
serine protease 3, partial3510.730.739XM_021345780.1
lipase3914.230.739XM_021326676.1
neutral lipase113.336.909XM_021331172.1
neutral lipase123.636.476XM_021331174.1
neutral lipase257.236.476XM_021331175.1
pancreatic lipase 2225.738.096XM_021338421.1
lipase112.436.224XM_021338439.1
pancreatic lipase 2218.736.224XM_021343929.1
neutral lipase123.635.29XM_021344621.1
inactive lipase125.730.895XM_021344707.1
Immune-related
tubulin alpha-1 chain-like12249.819XM_021338879.1
Lipocalins
fatty acid-binding protein 34433.314.744XM_021333801.1
fatty acid-binding protein 22216.415.066XM_021341051.1
fatty acid-binding protein 1119.714.986XM_021341061.1
fatty acid-binding protein 2234101.47XM_021341064.1
polycalin120.8101.47XM_021334936.1
Hexamerins
arylphorin8191182.226XM_021340131.1
arylphorin544.982.28XM_021340132.1
heat shock protein
heat shock protein 90221.882.63XM_021341131.1
heat shock protein375.573.029XM_021332476.1
others
Table 4. Proteins differentially expressed between XJ10 and XJ strain with a fold-change >1.5 and P-value < 0.05.
Table 4. Proteins differentially expressed between XJ10 and XJ strain with a fold-change >1.5 and P-value < 0.05.
Gene IDFC aPadj ValuemRNA FCDescription
XM_021332568.121.082260.0334240.24alpha-amylase 2-like
XM_021341184.12.1951510.031540.13chitin deacetylase 5α
XM_021328621.11.5767130.0048−0.45unconventional myosin-XV isoform X2
XM_021337877.1−1.860680.006606−0.91trypsin-like protease
XM_021337869.1−2.204310.018165−7.22trypsin
XM_021340592.1−2.732460.025215−1.36trypsin
XM_021337879.1−4.275660.035587-trypsin-like protease
XM_021330834.1−22.64290.0002191.40carboxypeptidase A
XM_021340037.1−23.74580.011483−1.53uncharacterized protein
XM_021345819.1−24.66710.036741−0.62chymotrypsin-like protease
XM_021345818.1−25.29130.032925−0.02chymotrypsin-like protease
XM_021340602.1−28.74180.017556−3.56trypsin T4
a Fold change of DEGs, positive value indicates up-regulation while negative value denotes down-regulation.

Share and Cite

MDPI and ACS Style

Jin, M.; Liao, C.; Chakrabarty, S.; Wu, K.; Xiao, Y. Comparative Proteomics of Peritrophic Matrix Provides an Insight into its Role in Cry1Ac Resistance of Cotton Bollworm Helicoverpa armigera. Toxins 2019, 11, 92. https://doi.org/10.3390/toxins11020092

AMA Style

Jin M, Liao C, Chakrabarty S, Wu K, Xiao Y. Comparative Proteomics of Peritrophic Matrix Provides an Insight into its Role in Cry1Ac Resistance of Cotton Bollworm Helicoverpa armigera. Toxins. 2019; 11(2):92. https://doi.org/10.3390/toxins11020092

Chicago/Turabian Style

Jin, Minghui, Chongyu Liao, Swapan Chakrabarty, Kongming Wu, and Yutao Xiao. 2019. "Comparative Proteomics of Peritrophic Matrix Provides an Insight into its Role in Cry1Ac Resistance of Cotton Bollworm Helicoverpa armigera" Toxins 11, no. 2: 92. https://doi.org/10.3390/toxins11020092

APA Style

Jin, M., Liao, C., Chakrabarty, S., Wu, K., & Xiao, Y. (2019). Comparative Proteomics of Peritrophic Matrix Provides an Insight into its Role in Cry1Ac Resistance of Cotton Bollworm Helicoverpa armigera. Toxins, 11(2), 92. https://doi.org/10.3390/toxins11020092

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