ORIGINAL RESEARCH article
iTRAQ-Based Comparative Proteomic Analysis of Larval Midgut From the Beet Armyworm, Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae) Challenged With the Entomopathogenic Bacteria Serratia marcescens
- Key Laboratory of Bio-Pesticide Innovation and Application of Guangdong Province, College of Agriculture, South China Agricultural University, Guangzhou, China
Entomopathogenic bacteria Serratia marcescens is widely used as an environmentally friendly biocontrol agent against various pests, including Spodoptera exigua. Understanding the immune defense mechanism of S. exigua through comparative proteomic analysis can identify the key proteins expressed in response to the microbial infection. Here, we employed the as isobaric tags for relative and absolute quantification (iTRAQ) technique to investigate the effects of S. marcescens on the proteomic expression of S. exigua. Based on the molecular functional analysis, the differentially expressed proteins (DEPs) were mainly involved in the binding process and catalytic activities. Further bioinformatics analysis revealed important DEPs that played a crucial role in innate immunity of S. exigua with recognition (C-type lectin), melanization (propanol oxidase 3, serine protease, Serine-type carboxypeptidase activity, clip domain serine protease 4), antimicrobial activity (lysozyme, lysozyme-like, gloverin, cecropin B), detoxification (acetyl-CoA C-acetyltransferase, 3-dehydroecdysone 3-alpha-reductase, glucuronosyltransferase, glutathione S-transferase) and others. The Quantitative real-time PCR (qRT-PCR) results further indicated the significant upregulation of the immune-related genes in Spodoptera exigua following S. marcescens infection. To the best of our knowledge, this is the first iTRAQ based study to characterize S. marcescens mediated proteomic changes in S. exigua and identified important immune-related DEPs. The results of this study will provide an essential resource for understanding the host-pathogen interactions and the development of novel microbial biopesticides against various pests.
Environmentally friendly entomopathogenic microbes used as an alternative to chemicals for controlling different pests. Over the few decades, there has been a sustained research activity in different entomopathogenic microbes and their roles for agricultural pest controls (D’Argenio et al., 2001; Huszar and Imler, 2008; Labropoulou et al., 2008; Vonkavaara et al., 2008; Pascual et al., 2012). However, with increasing evidence of the development of resistance to existing microbial pesticides, scientists are forced to look for alternative biocontrol agents having a novel mode of action. S. marcescens is a gram-negative entomopathogenic bacterium, belongs to the family Enterobacteriaceae and is commonly found in different insects (Pineda-Castellanos et al., 2015; Bidari et al., 2018). It has been successfully applied to control S. litura (Aggarwal et al., 2015), Helicoverpa armigera (Mohan et al., 2011), Phyllophaga Blanchardi (Pineda-Castellanos et al., 2015), Heliothis virescens (Sikorowski et al., 2001), and other pests. The primary metabolite chitinases play a vital role in the insecticidal activities of S. marcescens by hydrolyzing the body wall and peritrophic membrane structure of the insects. Tandem-mass spectrometry (LC-MS/MS) also showed that active insecticidal protein of S. marcescens has high similarity with the Serralysin-like protein (Pineda-Castellanos et al., 2015). It has been reported that serralysin can suppress the cellular immunity of the insect by lowering the adhesive characteristics of the immunosurveillance cells (Ishii et al., 2014). However, the effectiveness and stabilities of microbial insecticides were hindered due to the development of resistance by the insects (Cory and Franklin, 2012). The efficiency of microbial insecticides may vary due to the influences of various biotic and abiotic factors (Bonaterra et al., 2012). Therefore, understanding the host defense mechanisms is essential to the development of effective entomopathogen-based biopesticides to manage pests in different environments.
The beet armyworm, S. exigua (Hübner) (Lepidoptera: Noctuidae), is an important worldwide pest that causes severe damage to the agricultural and ornamental plants (Xiu and Dong, 2007). The larvae of this devastating pest can feed both foliage and fruit and damages many crops, such as potato, corn, onion, cotton, soybean, and tomato (Greenberg et al., 2001; Saeed et al., 2010). Evidence from previous studies indicated that larvae of S. exigua consume more cabbage foliage than diamondback moth Plutella xylostella, but less than cabbage looper Trichoplusia ni (Hübner) (East et al., 1989). So far, the control strategy of this pest is achieved by chemical pesticides, which induce the development of resistance against a wide range of insecticides, including organophosphates, carbamates, and pyrethroids (Moulton et al., 2002). Furthermore, the intensive use of pesticides leads to environmental as well as food safety problems. Therefore, it is necessary to develop alternative pest control strategies, especially biopesticides against the resistance of this pest (Pascual et al., 2012). Among the biocontrol agent, microbial-based biopesticide is gaining importance due to smooth delivery, formulations, and target specific. They are also less toxic and effective in small quantities (Chattopadhyay et al., 2017). The larval stage of S. exigua consumes various plant materials, including pathogenic microorganism that comes into direct contact with the midgut. This activates gut immunity to protect themselves from the infection and to guarantee the normal development of the host (Lemaitre and Miguel-Aliaga, 2013; Wu et al., 2016). Toll pathway played a major role in the defense mechanism against the gram positive bacteria and fungi where as the immune deficiency (Imd) pathway used against the Gram negative bacteria. Viral pathogens are targeted through RNAi, apoptosis and autophagy. Other pathways, such as JNK, JAK/STAT, and p38 pathways are also played an important role in the insect innate immunity (Clem, 2005; Chen et al., 2010; Nakamoto et al., 2012; Wu et al., 2016). Thus analyzing differential expressed midgut proteins might reveal the host defense or detoxification mechanism of S. exigua.
With the advances of high throughput sequencing technology, millions of DNA molecules obtained from mRNA can be sequenced in parallel without the needs of bacterial clones, which significantly improve our understanding of the differentially expressed genes at the mRNA level (Cloonan and Grimmond, 2008; Morozova and Marra, 2008). However, the genomic analysis fails to analyze the functions of some critical proteins that depend on post-translational and protein-protein interactions. Proteomics offers a direct interpretation of protein responses under the influence of various biological perturbations, such as disease or drug treatment (Anderson and Anderson, 1998). Advanced proteomics method such as isobaric tags for relative and absolute quantification (iTRAQ) has been introduced to identify a higher number of proteins with a and more reliable quantitative information. This method used as an isotope labeling approach and successfully applied to study the changes in protein expression under different conditions (Song et al., 2018). Here we investigated whether S. exigua can raise a specific immune response to S. marcescens infections by comparing the gut proteomes of control and treated insects. We predicted that S. exigua might trigger its immune system in the gut to defense against the pathogen. Identifying and analyzing the differentially expressed proteins (DEPs) associated immunity will significantly improve our understanding of the molecular mechanism of the host defense system.
Materials and Methods
Insect Rearing and Bacterial Sample Preparation
A susceptible population of S. exigua was obtained from the College of Agriculture, South China Agricultural University, China, and kept in an insecticide-free environment for 10 generations. Adults were fed 10% honey solution, and larvae were reared on an artificial diet as described by the previous researcher (Ren et al., 2013). All populations were maintained at 25 ± 1°C under a photoperiod of 14:10 h (light:dark) and 68 ± 5% relative humidity (RH). The highly entomopathogenic strain S. marcescens PS-1 was kindly provided by Professor Zhang Maoxin of the College of Agriculture, South China Agricultural University, China. Bacteria from a glycerol stock (stored at −80°C) were inoculated into Luria-Bertani (LB) and grown 10 h at 30°C. The bacterial suspension was inoculated into fresh LB medium, incubated, and centrifuged at 5,000 r/min for 10 min. The pallet was resuspended in phosphate-buffered saline (PBS) in different concentrations and used for the bioassay experiments.
Insect Treatment With S. marcescens
Firstly, to test the toxicity of this S. marcescens and find the best concentration to treat the insect getting the midgut, we did the following experiment. The 5th instar larvae of S. exigua were selected and divided into two groups, each with 5 concentrations, each concentrations containing 24 insects, 3 repetitions. The treatment group includes larvae, which were fed for 120 h with an artificial diet containing different concentration of S. marcescens (1.0 × 104 CFU/mL, 1.0 × 105 CFU/mL, 1.0 × 106 CFU/mL, 1.0 × 107 CFU/mL, 1.0 × 108 CFU/mL, 1.0 × 109 CFU/mL) and the number of dead larvae was recorded at 24 h intervals. Whereas the Control group was fed an artificial diet devoid of S. marcescens. All larvae used in the feeding experiment were starved for 2 h before the exposure to S. marcescens. The percentage mortality was calculated using the standard procedure (Abbott, 1925; Thanigaivel et al., 2019). The lethal dose (LC50) was calculated using Probit analysis using Minitab®17 software. The bioassay experiment was repeated with the bacterial lethal dose obtained after 72 h of post-infection. After treatment of 24, 48, and 72 h, larval midguts were dissected from the treated and control group and proceeds for the isolation of protein and RNA. The concentrations of RNA were analyzed using Nanodrop (Bio-Rad, United States), and its integrity was determined on Agilent 2100 Bioanalyzer (Agilent, United States). The sample was frozen in liquid nitrogen and stored at −80°C for further use.
Midgut samples were placed in lysis buffer 3 containing 8 M Urea, 40 mM Tris- HCl, 2 mM EDTA, one mM dithiothreitol (DTT), and one mM phenylmethanesulfonyl fluoride (PMSF). The mixture was transferred to a tissue laser for 2 min at 50 Hz, followed by centrifuged with 25,000 g at 20 min. The supernatant was transferred into a new tube, reduced with 10 mM DTT at 56°C for 1 h, and alkylated by 55 mM iodoacetamide for 45 min. The protein concentration was measured using Bradford assay with bovine serum albumin (BSA) as a standard (Bradford, 1976).
iTRAQ Labeling and Strong Cation Exchange (SCX) Chromatography
Protein (100 μg) obtained from each sample was digested using Trypsin Gold (Promega, United States) at an enzyme to substrate ratio of 1: 50 at 37°C for 16 h. The peptides were desalted with C18 cartridge (Thermo Fisher Scientific, United States) to remove the high urea, and desalted peptides were dried by vacuum centrifugation. A peptide mixture of each sample was labeled using iTRAQ reagent according to the manufacturer’s instructions (Applied Biosystems). iTRAQ labeled peptides were fractionated by the SCX chromatography using the AKTA Purifier system (GE Healthcare). The dried peptide mixture was reconstituted and acidified with buffer A (10 mM KH2PO4 in 25% of ACN, pH 3.0) followed by loaded onto a PolySULFOETHYL column (5 μm, 200 Å, PolyLC Inc, Maryland, United States). The peptides were further eluted at a flow rate of 1 ml/min with a gradient of 0–8% buffer B (500 mM KCl, 10 mM KH2PO4 in 25% of ACN, pH 3.0) for 22 min. The collected fractions were finally desalted on the C18 Cartridges (EmporeTM SPE Cartridges C18 (standard density), bed I.D.7 mm, volume 3 ml, Sigma), and concentrated by vacuum centrifugation (Fan et al., 2018).
Each fraction was injected for nano LC-MS/MS analysis. The peptide mixture was loaded onto a reverse-phase trap column connected to the C18 reversed-phase analytical column in buffer A (0.1% Formic acid) and divided with a linear gradient of buffer B (84% acetonitrile and 0.1% Formic acid). LC-MS/MS analysis was carried out on a Q Exactive mass spectrometer (Thermo Fisher Scientific) that was coupled to Easy nLC (Proxeon Biosystems, Thermo Fisher Scientific) for 60/120/240 min. The mass spectrometer was operated in positive ion mode. MS data were obtained using a data-dependent top10 method dynamically by selecting the most abundant precursor ions present in the survey scan (300–1,800 m/z) for HCD fragmentation. The automatic gain control (AGC) target was set to 3e6, maximum injection time to 10 ms; dynamic exclusion duration was 40.0 s. Survey scans were acquired at a resolution of 70,000 at m/z 200, and resolution for the HCD spectra was set to 17,500 at m/z 200, and the isolation width was 2 m/z (Wang et al., 2016; Fan et al., 2018).
iTRAQ MS Raw Data Processing
The original mass spectrometry file (.raw) generated by Q Exactive was converted to a. mgf file using Proteome Discoverer 1.4 (Thermo Fisher Scientific) software and submitted to the MASCOT2.2 server for database retrieval through the built-in tools of the software. Then, through the Proteome Discoverer 1.4, the library file (. dat file) formed on the MASCOT server is transmitted back to the software, and the database is screened according to the FDR < 0.01 standard to obtain highly reliable qualitative results. The protein library used in this experiment was derived from the previously generated transcriptome in our laboratory.
Screening of DEPs
To compare the differences between different samples, FDR is often used to determine the domain value of p-value, and the differential expression multiple of the gene between different samples is calculated according to the expression level of the protein (RPKM value). The smaller the FDR value, the larger the difference multiplier, and the more significant the difference in expression. We finally screened DEPs with a protein whose FDR could not exceed 0.01 and a fold difference of no less than 1.2 times and sorted them for further analysis and verification.
In the present study, the functional annotations of the DEPs sequences were performed by utilizing the Blast2go_v2.5 program against the NCBInr and UniProt/SwissProt database. Functional analysis of DEPs was carried out using gene ontology (GO) analysis, and the DEPs were grouped according to the biological process, molecular function, and cellular component (Conesa et al., 2005). The Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis was used to classify and group the DEPs (Kanehisa and Goto, 2000).
Quantitative Real-Time PCR (qRT-PCR) Analysis
Fifteen well-known immune-related genes (JNK, NOX1, Lectin, IMD, Relish, TAB, IKK, Defensin, Gioverin, Lebocin, Lysozyme, PGRP-LC, Moricin, Hemolin, and Cecropin), as well as six, randomly selected DEPs were investigated by qRT-PCR at the transcription level. Total RNA was extracted by E.Z.N.A.® Total RNA Kit II method as manufacturer’s description (Omega Bio-tek, Inc.), and reverse-transcribed cDNA from equal amounts (1.0 μg) of total RNA using the PrimeScriptTM RT Master Mix (Perfect Real Time) (TaKaRa Biomedical Technology (Beijing) Co., Ltd.). List of Primers and their sequences used in RT-qPCR is shown in Table 1. qRT-PCR was carried out in a BioRad CXF96 Real-Time PCR Detection System (Bio-Rad Laboratories, Inc.) using TB Green® Premix Ex TaqTM II (Tli RNaseH Plus) (TaKaRa) according to the instructions of the manufacturer. The reaction program was set as initial denaturation at 95°C for 30 s, 40 cycles of 95°C for 5 s and 60°C for 30 s. The expression level of genes was calculated by the 2–ΔΔCt method, and the value stood for an n-fold difference relative to the calibrator (L10). Each experiment was performed in triplicate. All data were given in terms of relative mRNA expression as mean ± SE.
Figures were prepared using the tool GraphPad Prism 7.0. All experiments were performed in triplicates, and the statistical analysis was conducted using the SPSS software (SPSS 26.0, IBM). Data from mortality experiments were subjected to analysis of variance (ANOVA of arcsine, logarithmic and square root transformed percentages), with data expressed as a mean of five replicates. Significant differences between treatment groups were analyzed by Tukey’s multiple range test (significance at p < 0.05) using Minitab®17 program. One-way analysis of variation (ANOVA) followed by Least-significant difference (LSD) multiple comparisons was used to compare the relative expression of immune genes in different time intervals. Differences were considered statistically significant at a p < 0.05.
The experimental larvae were inoculated with a different bacterial concentration to characterize the interaction between S. exigua and S. marcescens. Feeding bioassays demonstrated that the oral ingestion of S. marcescens is lethal to the larvae of S. exigua. The mortality rate was significantly higher post 120 h at the maximum dosage of 1 × 109 [F(4, 20 = 15.11, P ≤ 0.001] and it is statistically significant with other treatment dosage and control. Despite there is no significant difference between 1 × 109 and 1 × 108 [F(4, 20 = 21.22, P ≤ 0.001]. The mortality rate is higher when the treatment dosage and the treatment hours gets increased. Similarly, the percentage of mortality was statistically different between the treatment dosages post 96 h and it was significantly higher in 1 × 109 [F(4, 20 = 10.57, P ≤ 0.001] as compared to other dosages and control. Though, there is no statistical difference between 1 × 105 and 1 × 104 [F(4, 20) = 17.44, P ≤ 0.001]. However, all the treatment dosage were statistically different with control. Despite there is no statistical difference between the treatment dosage of 1 × 109 and 1 × 108 [F(4, 20) = 20.88, P ≤ 0.001] post 72 h of treatment. Similarly, there is no statistical difference between 105 and 1 × 104 [F(4, 20) = 19.11, P ≤ 0.001]. Apart from that, the mortality rate was significantly different between the control and all treatment dosages [F(4, 20) = 12.22, P ≤ 0.001]. Similar trends were observed in the 48 h of treatment as the mortality rate significantly reduced as compared to 120 and 96 h of treatment. As there is no significant difference in the treatment dosage of 1 × 109 and 1 × 108 [F(4, 20) = 14.18, P ≤ 0.001] also between 10–5 and 1 × 10–4 [F(4, 20) = 13.36, P ≤ 0.001]. Although there is no statistical difference between all the treatment dosage post 24 h of treatment [F(4, 20 = 10.18, P ≤ 0.001]. The mortality rate was maximum recorded at 1 × 109 [F(4, 20) = 10.18, P ≤ 0.001], however, there is no statistical difference between the other treatment dosage. Though, all the treatment dosages were significant with control [F(4, 20) = 15.11, P ≤ 0.001] (Figure 1A). Probit analysis showed that the lethal concentration (LC50) at 72 h was 1.0 × 108 CFU/mL (Figure 1B).
Figure 1. Bioassay of S. exigua. (A) Larval mortality of fifth instar larvae of S. exigua post-treatment with different concentrations (CFU/ml) of S. marcescens. Each point represents mean ± SE of three replications from different sets of experiments of identical environmental conditions. Same letters above bars indicate no significant difference (P ≤ 0.05) according to a Tukey test. Y axis showed the mortality rate (%), and X-axis represent the induction time (h). (B) Probit analysis showing the lethal concentrations (LC50) of S. marcescens against fifth instars larvae of S. exigua. Y axis showed the mortality rate (%) at 72 h, and X-axis represent the concentration (log) of S. marcescens.
All MS/MS spectra were processed using Mascot software. iTRAQ analysis of midgut proteome showed that a total of 291,113 spectra, including 49,292 unique spectra, 24,135 identified peptides, 23,224 unique peptides resulted in 4,559 protein (Supplementary Figure S1 and Supplementary Table S1).
Identification of Differentially Expressed Proteins (DEPs)
In the present study, 211 proteins (cut off of 1.2-fold change and p < 0.05) were found to be differentially expressed in S. exigua after exposure with S. marcescens. Comparative proteomic analysis showed that 86 proteins were upregulated, and 125 were downregulated after exposure to S. marcescens with a fold change of 0.24–2.08 (Figure 2). The top 10 up- and downregulated proteins by fold change in response to bacterial infection are listed in Table 2. Top up-regulated proteins belonged to cubilin, nicotinic acetylcholine receptor subunit beta 2, G-protein coupled receptor Mth2-like, Fibroblast growth factor receptor 1-A-like, and lactase-phlorizin hydrolase nesprin-1. The top 10 down-regulated proteins included NACHT, LRR, and PYD domain-containing protein 3-like, double-stranded RNA-specific disease 1-like, ADP, ATP carrier protein-like, SEC14-like protein, etc. (Table 2).
Figure 2. Differentially expressed proteins (DEPs) in S. exigua responding to S. marcescens. (A) Statistics of up- and downregulated DEPs in S. exigua. (B) Distribution of fold changes and molecular masses for DEPs. Red circles indicate upregulated proteins, and green circles indicate downregulated proteins.
Functional Classification of DEPs
Based upon Gene ontology (GO) classification of the DEPs were classified into three categories; biological process, molecular function, and cellular component. Molecular functional classification showed that the proteins identified by iTRAQ analysis were classified into different categories; the majority of them contain catalytic and binding activities (Figure 3). KEGG1 ontology assignments were used to classify functional annotations of the identified DEPs to further understand their biological functions and identified various pathways associated with Phagosome (3.8%), Metabolism of xenobiotics by cytochrome P450 (3.2%), Chemical carcinogenesis (3.2%); Lysosome (3.2%), Carbon metabolism (2.5%), Drug metabolism – cytochrome P450 (2.5%), Pentose and glucuronate interconversions (2.5%), Glutathione metabolism (2.5%), Ribosome (2.5%), and so on (Figure 4).
Figure 3. Gene ontology (GO) analysis of all proteins in midguts identified by iTRAQ analysis. Shown above is the classification of these proteins in different categories based on biological process, molecular function, cellular component.
Figure 4. KEGG functional classification of the identified proteins in the midgut of S. exigua identified by iTRAQ analysis.
DEPs Involved in Carbohydrate and Energy Metabolism
S. marcescens infection alters the expression of many proteins involved in carbohydrate metabolism in S. exigua, which was mainly included the citrate cycle (TCA cycle) metabolism of Fructose and mannose, Galactose, Ascorbate and aldarate, Pyruvate, Glyoxylate and dicarboxylate, Propanoate, Butanoate and, Inositol phosphate. Among these DEPs, six proteins were up-regulated subsequent to S. marcescens treatment, with the highest ratio 1.50 ± 0.15 for phosphatidylinositol 4-phosphatase followed by glucuronosyltransferase 1.43 ± 0.08. Other identified upregulated proteins were 2-oxoglutarate dehydrogenase E1 component (1.35 ± 0.09); fumarate hydratase, class II (1.24 ± 0.06); aldehyde reductase (1.21 ± 0.05); acetyl-CoA C-acetyltransferase (1.25 ± 0.05). The only downregulated protein involved in the carbohydrate metabolism was synaptojanin (0.59 ± 0.08) (Table 3). DEPs involved in the energy metabolism were also identified in S. exigua infected by S. marcescens and were related to Oxidative phosphorylation, Carbon fixation pathways in prokaryotes, and Nitrogen metabolism. Five identified upregulated proteins associated with energy metabolism were V-type H+-transporting ATPase 21 kDa proteolipid subunit; NADH dehydrogenase (ubiquinone) 1 beta subcomplex subunit 1; acetyl-CoA C-acetyltransferase; fumarate hydratase, class II; glutamate dehydrogenase with the average ratio ranged from 1.24 ± 0.06 to 1.34 ± 0.07. However, none of the downregulated proteins associated with the energy metabolism were identified (Table 4).
DEPs Involved in Immunity
Based on the proteomic analysis, 27 DEPs were identified associated with the immune defense system of S. exigua (Table 5). These proteins were involved in recognition, melanization, antimicrobial activity, detoxification, and other functions. In the present study, the level of three identified recognition proteins belongs to C-type lectin and galectin-8-like were upregulated, with the ratio ranged from 1.07 ± 0.02 to 1.13 ± 0.07. Identified DEPs involved in the melanization process were belongs to PPO3, a serine protease, Serine-type carboxypeptidase activity, clip domain serine protease 4, and their expression was ranged with the ratio from 0.68 ± 0.04 to 1.32 ± 0.07. Moreover, four AMPs and seven detoxifying proteins were identified with increased levels from DEPs, with the ratio ranged from 1.06 ± 0.30 to 1.25 ± 0.09, and 1.20 ± 0.02 to 1.43 ± 0.07, respectively.
Transcriptional (qRT-PCR) Analysis of the Immune-Related Genes
Quantitative real-time PCR was used to detect the expression levels of 15 immune-related genes and six following post-infection of S. marcescens. Results illustrate that most of the studied immune genes were significantly upregulated in S. exigua following the post-infection at different time intervals. The expression level of JNK, NOX1, PGRP-LC, Unigene29765, and Unigene19139 were significantly higher (P < 0.05) in the treatment group than those in the control group, and Lebocin, Lysozyme, Gioverin, Defensin, Moricin, Cecropin, and Hemolin was significantly (P < 0.01) higher than that in the control group following 72 h of post-infection. However, the expression levels of IKK, Relish in the treated group was highly significantly (P < 0.01) lesser than that in the control group following 72 h of post-infection. Overall, similar trends were observed between the expression of the gene at both transcriptomics and proteomics levels among the treated larvae of S. exigua (Figure 5).
Figure 5. Quantitative real-time PCR (qRT-PCR) analysis of 21 genes (A–U for JNK, NOX1, Lectin, IMD, Relish, TAB, IKK, Defensin, Gioverin, Lebocin, Lysozyme, PGRP-LC, Moricin, Hemolin, Cecropin, CL1044, Unigene 14488, Unigene 19139, Unigene 4854, Unigene50 and Unigene 29765, respectively) of the midgut of S. exigua, which was fed with an artificial diet containing S. marcescens (1.0 × 108 CFU/mL) for 24, 48, and 72 hours. L10 was used as an internal control. Error bars represent mean ± SD from three independent experiments. *p < 0.05, **p < 0.01. Same letters above bars indicate no significant difference (P ≤ 0.05).
Understanding the host defense system is an essential aspect of improving existing pest management strategies. Despite recent advances in insect science, the molecular mechanisms underlying the host defense system of most of the destructive pests not has been fully understood. iTRAQ based proteomic technique can identify numerous proteins and quantify them more reliably compared to the traditional two-dimensional electrophoresis (Karp et al., 2010). This method has been applied successfully to screen the DEPs to understand the immune defense mechanism of the host (Malécot et al., 2011; Gao et al., 2017; Grassl et al., 2017). However, limited studies are available on the proteomic analysis of insect pests. To the best of our knowledge, this is the first iTRAQ study on the molecular mechanisms of the host defense system in S. exigua against bacterial infections. The differentially expressed immune proteins of S. exigua induced by the entomopathogenic bacteria S. marcescens is discussed.
Changes in proteomic of the host following infection illustrate important clues to host defense systems, the pathogenesis mechanism, and suggests target candidate genes or pathways for developing control strategies (Munoz-Gomez et al., 2014; Yuan et al., 2015; Gao et al., 2017). In the present study, 211 DEPs of all detected 4,559 proteins, were up- or downregulated in S. exigua induced by entomopathogenic bacteria. These DEPs were associated with various biological functions and metabolic pathways associated with the with Phagosome, Metabolism of xenobiotics by cytochrome P450, Chemical carcinogenesis; Lysosome, Carbon metabolism, Drug metabolism – cytochrome P450, Pentose and glucuronate interconversions, Glutathione metabolism, Ribosome, etc., indicating that bacterial infection strongly influences the host physiology. During the microbial infection in insects, partial energy originally used for host protein synthesis will be assigned for the reproduction of the microorganism. In the present study, several proteins involved in carbohydrate and energy metabolism appeared to be differentially expressed during S. marcescens infection, including 2-oxoglutarate dehydrogenase E1 component; fumarate hydratase, class II; aldehyde reductase; glucuronosyltransferase; acetyl-CoA C-acetyltransferase; synaptojanin and phosphatidylinositol 4-phosphatase; V-type H+-transporting ATPase 21 kDa proteolipid subunit; NADH dehydrogenase (ubiquinone) 1 beta subcomplex subunit one and glutamate dehydrogenase. The upregulation of the metabolic pathways of the host is necessary to meet the basic materials and energy requirements for the growth and reproduction of the host during bacterial infection (Yu et al., 2015). Taken together, these results suggest that energetic metabolic pathways, such as carbohydrate and energy metabolism, are vigorously triggered in the host cell during the infection process to survive against the pathogen.
Once infected by the pathogen, the host activates its defense system produces a rapid and effective response to locate and neutralize the multiplication of the pathogen. Several DEGs associated with the host immune defense mechanism were also identified from different insect species (Malécot et al., 2011; Gao et al., 2017; Grassl et al., 2017). In this study, differentially expressed immune proteins involved in recognition, modulation, signaling, and other immune defense mechanisms were identified after infection with S. marcescens. During infection in insects, multiple Pattern Recognition Receptors (PRR) recognize conservative determinants such as peptidoglycan, lipopolysaccharide, b 1-3 glucan, and lipoteichoic acid) on the pathogen surface to trigger defense reactions (Jiang, 2008). In the present study, C-type lectin (CTL) was identified in the DEP dataset. CTLs are a large family of calcium-dependent carbohydrate-binding proteins and are mainly associated with cell adhesion and immunity by recognition of various glycoconjugates (Shen et al., 2018). Upregulation of this protein (CL5489.Contig1) following 24 h of post-infection may help in the innate immunity of S. exigua to efficiently recognize and eliminate pathogens (Adachi et al., 2004).
It has been reported that antimicrobial peptides (Hampson, 1892), such as cecropins, glove rings, attacins, or REsponse to PAThogen (REPAT), were differentially expressed in the host upon pathogen induction (Ji and Kim, 2004; Freitak et al., 2007; Herrero et al., 2007; Van Munster et al., 2007; Wu et al., 2010; Hwang and Kim, 2011; Bang et al., 2012). Similarly, the present study also observed a significant upregulation of AMP-Cecropins, which is known to exhibit broad-spectrum antimicrobial activity against both Gram-positive and Gram-negative bacteria (Wu et al., 2018). The DEP data mining also led to the identification of upregulated gloverin (Unigene18150) protein in S. exigua, which was reported to be associated with the targeting of the Gram-negative bacteria and fungi (Yi et al., 2014). This indicates that the upregulation of the expression of this peptide is a way of protecting against bacterial virulence factors. The AMP lysozyme is involved in the disruption of bacterial cells, which is considered as an essential function in the defense against the pathogen (Leippe, 1999; Schulenburg and Boehnisch, 2008). Two AMP proteins, such as lysozyme (CL1493.Contig1_All) and lysozyme-like protein (CL4703.Contig1_All), showed increased expression in S. exigua. Similar to the iTRAQ result, significant upregulation was also observed among AMPs such as Defensin, Gloverin, Lebocin, Lysozyme, PGRP-LC, Moricin, Hemolin, and Cecropin at the transcription level, which further indicates that the antimicrobial peptide plays a vital role in the defense system of S. exigua. A similar observation was observed in T. ni and H. armigera after the ingestion of B. thuringiensis (Tamez-Guerra et al., 2008; Wang et al., 2010).
Melanization is another crucial process of the insect defense system and is triggered by trace transduction factors when encountering foreign invaders (Eleftherianos and Revenis, 2011; Nakhleh et al., 2017). The present study identified the upregulation of propanol oxidase (PPO), an essential immune protein associated with melanization, and intermediates formed during this process might allow the host to kill the pathogen. Bacteria produce a large number of toxins inside the host to combat the host defense system. Glutathione S-transferase has been reported to be associated with the antiviral and antibacterial activity. They help in the maturation and formation of the phagosome and subsequent clearance of the pathogens through innate immune responses (Mahlapuu et al., 2016). The upregulation of glucuronosyltransferase in the present study might be associated with the detoxification of the S. marcescens toxin. A similar observation was reported in Bombyx mori infected by Bacillus bombyseptieus (Molina-Heredia et al., 2006; Huang et al., 2009). The present study also revealed other immune-related DEPs involved in the immune response of S. exigua. For example, Integrin was significantly upregulated, while Hexamerin was downregulated. Integrins transmembrane receptors that play a vital function in cell-to-cell and cell-to-extracellular matrix interactions (Gramkow et al., 2010; Zhang et al., 2017). Integrin from A. gambiae appears to function as a phagocytic receptor for Gram-negative bacteria (Moita et al., 2006). C-type lectin can interact with β-integrin, which leads to the increased hemocyte encapsulation in Helicoverpa armigera (Wang et al., 2017). However, the hexameric protein is mainly involved in the accumulation of amino acids during insect metamorphosis but may also have an essential function in immunity (Martins and Bitondi, 2012).
The present study applied the iTRAQ technique for the first time to understand the effect of S. marcescens on the defense system of S. exigua at different time points. A total of 4,559 proteins were identified with 86 up-regulated and 125 down-regulated proteins and classified into various GO and KEGG functional groups. Our research indicated that the infection of S. marcescens could affect the different metabolic pathways of S. exigua associated with the carbohydrate metabolism, energy metabolism as well as immune defense system involving recognition, melanization, antimicrobial activity, detoxification and others. Taken together, this study identified several important DEPs in S. exigua infected by S. marcescens, which can provide the foundation for further functional studies aimed at illustrating the molecular mechanism of the host defense system in S. exigua.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
FJ, SD, and XX conceived and designed the experiments. BL, MS, YL, XX, and SD performed the experiments. SD, XX, and YL analyzed the data. FJ and XX contributed to reagents, materials, and analysis tools. SD and XX wrote the manuscript. SD and FJ revised the manuscript.
This work was supported by a grant from the National Natural Science Foundation of China (31972345), National Natural Science Foundation of Guangdong, China (2019A1515011221, 2018A030313402), and Guangzhou Science and Technology project (201803020011).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.00442/full#supplementary-material
FIGURE S1 | Basic information statistics of the proteome by iTRAQ. Total spectra: the secondary mass spectrums; Spectra: the secondary mass spectrums after quality control; Peptide: the identified peptide after quality control; Unique peptide: the identified peptides which belong only to a group of proteins; Protein: identified by Mascot 2.3.02 software.
TABLE S1 | Protein identification list.
Aggarwal, C., Paul, S., Tripathi, V., Paul, B., and Khan, M. A. (2015). Chitinolytic activity in Serratia marcescens (strain SEN) and potency against different larval instars of Spodoptera litura with effect of sublethal doses on insect development. BioControl 60, 631–640.
Bang, K., Park, S., Yoo, J. Y., and Cho, S. (2012). Characterization and expression of attacin, an antibacterial protein-encoding gene, from the beet armyworm, Spodoptera exigua (Hübner)(Insecta: Lepidoptera: Noctuidae). Mol. Biol. Rep. 39, 5151–5159. doi: 10.1007/s11033-011-1311-3
Bidari, F., Shams−Bakhsh, M., and Mehrabadi, M. (2018). Isolation and characterization of a Serratia marcescens with insecticidal activity from Polyphylla olivieri (Col.: Scarabaeidae). J. Appl. Entomol. 142, 162–172.
Bonaterra, A., Badosa, E., Cabrefiga, J., Francés, J., and Montesinos, E. (2012). Prospects and limitations of microbial pesticides for control of bacterial and fungal pomefruit tree diseases. Trees 26, 215–226. doi: 10.1007/s00468-011-0626-y
Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1006/abio.1976.9999
Chattopadhyay, P., Banerjee, G., and Mukherjee, S. (2017). Recent trends of modern bacterial insecticides for pest control practice in integrated crop management system. 3 Biotech 7:60. doi: 10.1007/s13205-017-0717-6
Chen, J., Xie, C., Tian, L., Hong, L., Wu, X., and Han, J. (2010). Participation of the p38 pathway in Drosophila host defense against pathogenic bacteria and fungi. Proc. Natl. Acad. Sci.U.S.A. 107, 20774–20779. doi: 10.1073/pnas.1009223107
Clem, R. (2005). “The role of apoptosis in defense against baculovirus infection in insects,” in Role of Apoptosis in Infection, ed. D. E. Griffin (Berlin: Springer), 113–129. doi: 10.1007/3-540-27320-4_5
Conesa, A., Götz, S., García-Gómez, J. M., Terol, J., Talón, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
East, D., Edelson, J., and Cartwright, B. (1989). Relative cabbage consumption by the cabbage looper (Lepidoptera: Noctuidae), beet armyworm (Lepidoptera: Noctuidae), and diamondback moth (Lepidoptera: Plutellidae). J. Econ. Entomol. 82, 1367–1369.
Fan, W., Ge, G., Liu, Y., Wang, W., Liu, L., and Jia, Y. (2018). Proteomics integrated with metabolomics: analysis of the internal causes of nutrient changes in alfalfa at different growth stages. BMC Plant Biol. 18:78. doi: 10.1186/s12870-018-1291-8
Freitak, D., Wheat, C. W., Heckel, D. G., and Vogel, H. (2007). Immune system responses and fitness costs associated with consumption of bacteria in larvae of Trichoplusia ni. BMC Biol. 5:56. doi: 10.1186/1741-7007-5-56
Gao, K., Deng, X.-Y., Shang, M.-K., Qin, G.-X., Hou, C.-X., and Guo, X.-J. (2017). iTRAQ-based quantitative proteomic analysis of midgut in silkworm infected with Bombyx mori cytoplasmic polyhedrosis virus. J. Proteom. 152, 300–311. doi: 10.1016/j.jprot.2016.11.019
Gramkow, A. W., Perecmanis, S., Sousa, R. L. B., Noronha, E. F., Felix, C. R., Nagata, T., et al. (2010). Insecticidal activity of two proteases against Spodoptera frugiperda larvae infected with recombinant baculoviruses. Virol. J. 7:143. doi: 10.1186/1743-422X-7-143
Grassl, J., Peng, Y., Baer-Imhoof, B., Welch, M., Millar, A. H., and Baer, B. (2017). Infections with the sexually transmitted pathogen Nosema apis trigger an immune response in the seminal fluid of honey bees (Apis mellifera). J. Proteom. Res. 16, 319–334. doi: 10.1021/acs.jproteome.6b00051
Greenberg, S., Sappington, T., Legaspi, B., Liu, T.-X., and Setamou, M. (2001). Feeding and life history of Spodoptera exigua (Lepidoptera: Noctuidae) on different host plants. Ann. Entomol. Soc. Am. 94, 566–575.
Herrero, S., Ansems, M., Van Oers, M. M., Vlak, J. M., Bakker, P. L., and de Maagd, R. A. (2007). REPAT, a new family of proteins induced by bacterial toxins and baculovirus infection in Spodoptera exigua. Insect Biochem. Mol. Biol. 37, 1109–1118. doi: 10.1016/j.ibmb.2007.06.007
Huang, L., Cheng, T., Xu, P., Cheng, D., Fang, T., and Xia, Q. (2009). A genome-wide survey for host response of silkworm, Bombyx mori during pathogen Bacillus bombyseptieus infection. PLoS One 4:e8098. doi: 10.1371/journal.pone.0008098
Hwang, J., and Kim, Y. (2011). RNA interference of an antimicrobial peptide, gloverin, of the beet armyworm, Spodoptera exigua, enhances susceptibility to Bacillus thuringiensis. J. Invertebr. Pathol. 108, 194–200. doi: 10.1016/j.jip.2011.09.003
Ishii, K., Adachi, T., Hamamoto, H., and Sekimizu, K. (2014). Serratia marcescens suppresses host cellular immunity via the production of an adhesion-inhibitory factor against immunosurveillance cells. J. Biol. Chem. 289, 5876–5888. doi: 10.1074/jbc.M113.544536
Ji, D., and Kim, Y. (2004). An entomopathogenic bacterium, Xenorhabdus nematophila, inhibits the expression of an antibacterial peptide, cecropin, of the beet armyworm, Spodoptera exigua. J. Insect physiol. 50, 489–496. doi: 10.1016/j.jinsphys.2004.03.005
Karp, N. A., Huber, W., Sadowski, P. G., Charles, P. D., Hester, S. V., and Lilley, K. S. (2010). Addressing accuracy and precision issues in iTRAQ quantitation. Mol. Cell. Proteom. 9, 1885–1897. doi: 10.1074/mcp.M900628-MCP200
Labropoulou, V., Douris, V., Stefanou, D., Magrioti, C., Swevers, L., and Iatrou, K. (2008). Endoparasitoid wasp bracovirus−mediated inhibition of hemolin function and lepidopteran host immunosuppression. Cell. Microbiol. 10, 2118–2128. doi: 10.1111/j.1462-5822.2008.01195.x
Mahlapuu, M., Håkansson, J., Ringstad, L., and Björn, C. (2016). Antimicrobial peptides: an emerging category of therapeutic agents. Front. Cell. Infect. Microbiol. 6:194. doi: 10.3389/fcimb.2016.00194
Mohan, M., Selvakumar, G., Sushil, S. N., Bhatt, J. C., and Gupta, H. S. (2011). Entomopathogenicity of endophytic Serratia marcescens strain SRM against larvae of Helicoverpa armigera (Noctuidae: Lepidoptera). World J. Microbiol. Biotechnol. 27, 2545–2551.
Moita, L. F., Vriend, G., Mahairaki, V., Louis, C., and Kafatos, F. C. (2006). Integrins of Anopheles gambiae and a putative role of a new β integrin, BINT2, in phagocytosis of E. coli. Insect Biochem. Mol. Biol. 36, 282–290. doi: 10.1016/j.ibmb.2006.01.004
Molina-Heredia, F. P., Houée-Levin, C., Berthomieu, C., Touati, D., Tremey, E., Favaudon, V., et al. (2006). Detoxification of superoxide without production of H2O2: antioxidant activity of superoxide reductase complexed with ferrocyanide. Proc. Natl. Acad. Sci. 103, 14750–14755. doi: 10.1073/pnas.0510828103
Moulton, J. K., Pepper, D. A., Jansson, R. K., and Dennehy, T. J. (2002). Pro-active management of beet armyworm (Lepidoptera: Noctuidae) resistance to tebufenozide and methoxyfenozide: baseline monitoring, risk assessment, and isolation of resistance. J. Econ. Entomol. 95, 414–424. doi: 10.1603/0022-0493-95.2.414
Munoz-Gomez, A., Corredor, M., Benitez-Paez, A., and Pelaez, C. (2014). Development of quantitative proteomics using iTRAQ based on the immunological response of Galleria mellonella larvae challenged with Fusarium oxysporum microconidia. PLoS One 9:e112179. doi: 10.1371/journal.pone.0112179
Nakamoto, M., Moy, R. H., Xu, J., Bambina, S., Yasunaga, A., Shelly, S. S., et al. (2012). Virus recognition by Toll-7 activates antiviral autophagy in Drosophila. Immunity 36, 658–667. doi: 10.1016/j.immuni.2012.03.003
Pascual, L., Jakubowska, A. K., Blanca, J. M., Cañizares, J., Ferré, J., Gloeckner, G., et al. (2012). The transcriptome of Spodoptera exigua larvae exposed to different types of microbes. Insect Biochem. Mol. Biol. 42, 557–570. doi: 10.1016/j.ibmb.2012.04.003
Pineda-Castellanos, M., Rodríguez-Segura, Z., Villalobos, F., Hernández, L., Lina, L., and Nuñez-Valdez, M. (2015). Pathogenicity of isolates of Serratia marcescens towards larvae of the scarab Phyllophaga blanchardi (Coleoptera). Pathogens 4, 210–228. doi: 10.3390/pathogens4020210
Ren, X.-L., Chen, R.-R., Zhang, Y., Ma, Y., Cui, J.-J., Han, Z.-J., et al. (2013). A Spodoptera exigua cadherin serves as a putative receptor for Bacillus thuringiensis Cry1Ca toxin and shows differential enhancement of Cry1Ca and Cry1Ac toxicity. Appl. Environ. Microbiol. 79, 5576–5583. doi: 10.1128/AEM.01519-13
Shen, D., Wang, L., Ji, J., Liu, Q., and An, C. (2018). Identification and characterization of C-type lectins in Ostrinia furnacalis (Lepidoptera: Pyralidae). J. Insect Sci. 18:24. doi: 10.1093/jisesa/iey011
Song, X.-Q., Zhao, Y., Weng, Q.-Y., Yuan, J.-C., Dong, Z.-P., Zhao, Z.-H., et al. (2018). Proteomic analysis of Zhangzagu3 (Setaria italica) and its parents based on iTRAQ technique. Biotechnol. Biotechnol. Equip. 32, 1407–1417.
Tamez-Guerra, P., Valadez-Lira, J., Alcocer-González, J., Oppert, B., Gomez-Flores, R., Tamez-Guerra, R., et al. (2008). Detection of genes encoding antimicrobial peptides in Mexican strains of Trichoplusia ni (Hübner) exposed to Bacillus thuringiensis. J. Invertebr. Pathol. 98, 218–227. doi: 10.1016/j.jip.2008.02.008
Thanigaivel, A., Chanthini, K. M.-P., Karthi, S., Vasantha-Srinivasan, P., Ponsankar, A., Sivanesh, H., et al. (2019). Toxic effect of essential oil and its compounds isolated from Sphaeranthus amaranthoides Burm. f. against dengue mosquito vector Aedes aegypti Linn. Pesticide Biochem. Physiol. 160, 163–170. doi: 10.1016/j.pestbp.2019.08.006
Van Munster, M., Prefontaine, G., Meunier, L., Elias, M., Mazza, A., Brousseau, R., et al. (2007). Altered gene expression in Choristoneura fumiferana and Manduca sexta in response to sublethal intoxication by Bacillus thuringiensis Cry1Ab toxin. Insect Mol. Biol. 16, 25–35. doi: 10.1111/j.1365-2583.2006.00692.x
Vonkavaara, M., Telepnev, M. V., Rydén, P., Sjöstedt, A., and Stöven, S. (2008). Drosophila melanogaster as a model for elucidating the pathogenicity of Francisella tularensis. Cell. Microbiol. 10, 1327–1338. doi: 10.1111/j.1462-5822.2008.01129.x
Wang, P., Zhuo, X.-R., Tang, L., Liu, X.-S., Wang, Y.-F., Wang, G.-X., et al. (2017). C-type lectin interacting with β-integrin enhances hemocytic encapsulation in the cotton bollworm, Helicoverpa armigera. Insect Biochem. Mol. Biol. 86, 29–40. doi: 10.1016/j.ibmb.2017.05.005
Wang, Y., Guan, H., Xie, D.-F., Xie, Y., Liu, X.-D., Wang, Q., et al. (2016). Proteomic analysis implicates dominant alterations of RNA metabolism and the proteasome pathway in the cellular response to carbon-ion irradiation. PLoS One 11:e0163896. doi: 10.1371/journal.pone.0163896
Wu, S., Zhang, X., He, Y., Shuai, J., Chen, X., and Ling, E. (2010). Expression of antimicrobial peptide genes in Bombyx mori gut modulated by oral bacterial infection and development. Dev. Comp. Immunol. 34, 1191–1198. doi: 10.1016/j.dci.2010.06.013
Xiu, W.-M., and Dong, S.-L. (2007). Molecular characterization of two pheromone binding proteins and quantitative analysis of their expression in the beet armyworm, Spodoptera exigua Hübner. J. Chem. Ecol. 33, 947–961. doi: 10.1007/s10886-007-9277-2
Yu, Q., Xiong, Y., Gao, H., Liu, J., Chen, Z., Wang, Q., et al. (2015). Comparative proteomics analysis of Spodoptera frugiperda cells during Autographa californica multiple nucleopolyhedrovirus infection. Virol. J. 12:115. doi: 10.1186/s12985-015-0346-9
Yuan, L.-L., Chen, X., Zong, Q., Zhao, T., Wang, J.-L., Zheng, Y., et al. (2015). Quantitative proteomic analyses of molecular mechanisms associated with cytoplasmic incompatibility in Drosophila melanogaster induced by Wolbachia. J. Proteom. Res. 14, 3835–3847. doi: 10.1021/acs.jproteome.5b00191
Keywords: Serratia marcescens, Spodoptera exigua, host immune defense, DEPs, iTRAQ
Citation: De Mandal S, Lin B, Shi M, Li Y, Xu X and Jin F (2020) iTRAQ-Based Comparative Proteomic Analysis of Larval Midgut From the Beet Armyworm, Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae) Challenged With the Entomopathogenic Bacteria Serratia marcescens. Front. Physiol. 11:442. doi: 10.3389/fphys.2020.00442
Received: 11 January 2020; Accepted: 08 April 2020;
Published: 08 May 2020.
Edited by:Senthil-Nathan Sengottayan, Manonmaniam Sundaranar University, India
Reviewed by:Sengodan Karthi, Manonmaniam Sundaranar University, India
You Ming Hou, Fujian Agriculture and Forestry University, China
Amrita Panda, Kumaun University, India
Copyright © 2020 De Mandal, Lin, Shi, Li, Xu and Jin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.