Transcriptional profiling analysis of Spodoptera litura larvae challenged with Vip3Aa toxin and possible involvement of trypsin in the toxin activation

Comparative transcriptome analysis reveals disruption of Plutella xylostella immune system by fungal peptide cyclosporin C

The diamondback moth (Plutella xylostella), a major threat to crucifers across the globe, has developed resistance against the majority of insecticides enhancing the need for alternate control measures against this pest. Recently cyclosporin C, a secondary metabolite produced by the insect pathogenic fungus Purpeocillium lilacinum, has been reported to induce lethal and sub-lethal effects against P. xylostella. To date, little is known about the molecular mechanisms of interaction between cyclosporin C and P. xylostella immune systems. This study reports the transcriptome-based immune response of P. xylostella to cyclosprin C treatment. Our results showed differential expression of 322, 97, and 504 differentially expressed genes (DEGS) in P. xylostella treated with cyclosporin C compared to control 24, 48, and 72 h post-treatment, respectively. Thirteen DEGs were commonly expressed at different time intervals in P. xylostella larvae treated with cyclosporin C compared to control. Cyclosporin C treatment induced the down-regulated expression of majority of immune-related genes related to pattern recognition responses, signal modulation, Toll and IMD pathways, antimicrobial peptides and antioxidant responses confirming the ability to suppress immune response of P. xylostella. These results will further improve our knowledge of the infection mechanism and complex biochemical processes involved in interaction between cyclosporin C and insect immune systems.

N-Terminal α-Helices in Domain I of Bacillus thuringiensis Vip3Aa Play Crucial Roles in Disruption of Liposomal Membrane

Bacillus thuringiensis Vip3 toxins form a tetrameric structure crucial for their insecticidal activity. Each Vip3Aa monomer comprises five domains. Interaction of the first four α-helices in domain I with the target cellular membrane was proposed to be a key step before pore formation. In this study, four N-terminal α-helix-deleted truncations of Vip3Aa were produced and, it was found that they lost both liposome permeability and insecticidal activity against Spodoptera litura. To further probe the role of domain I in membrane permeation, the full-length domain I and the fragments of N-terminal α-helix-truncated domain I were fused to green fluorescent protein (GFP), respectively. Only the fusion carrying the full-length domain I exhibited permeability against artificial liposomes. In addition, seven Vip3Aa-Cry1Ac fusions were also constructed by combination of α-helices from Vip3Aa domains I and II with the domains II and III of Cry1Ac. Five of the seven combinations were determined to show membrane permeability in artificial liposomes. However, none of the Vip3Aa-Cry1Ac combinations exhibited insecticidal activity due to the significant reduction in proteolytic stability. These results indicated that the N-terminal helix α1 in the Vip3Aa domain I is essential for both insecticidal activity and liposome permeability and that domain I of Vip3Aa preserved a high liposome permeability independently from domains II-V.

Cry9A and Vip3A protein-induced transcriptional changes correspond to their synergistic damage to the midgut of Chilo suppressalis

Cry and Vip3 proteins are both pore-forming toxins produced by Bacillus thuringiensis that show synergistic insecticidal activity against different insect pests. However, the synergistic effect of Cry and Vip3 proteins on the midgut in target insects is still unclear. In this study, faster and more serious damage was observed after treatment with both Cry9A and Vip3A proteins in the Chilo suppressalis midgut compared to single-protein treatment. Through RNA sequencing, midgut transcriptomic comparison was performed between dual- and single-protein treatments according to midgut injury. After 6 h, 609 differentially expressed genes were found with the combined Cry9A and Vip3A treatments, which was much more than that in the single treatment, corresponding to faster and more serious damage. These genes were mainly enriched in similar pathways, such as lipid metabolic, oxidation-reduction and carbohydrate metabolic process, peptide secretion and cell-cell adhesion; however, the number and expression level of differentially expressed genes are increased. For specific genes significantly regulated by induction of Cry9A and Vip3A, lipases, phospholipid scramblase, probable tape measure protein and arylsulfatase J were significantly downregulated after 6 h treatment. In addition, regular genes related to the activation and receptor binding of B. thuringiensis toxins were differentially regulated, such as ATP-binding cassette subfamily G member 1 and serine protease. Validation with RT-qPCR showed agreement with the sequencing results. Overall, our results support that stronger and faster midgut responses at the cellular and transcriptional levels are induced by the synergistic toxicity of Cry9A and Vip3A in C. suppressalis.

Sequential and Simultaneous Interactions of Plant Allelochemical Flavone, Bt Toxin Vip3A, and Insecticide Emamectin Benzoate in Spodoptera frugiperda

Target pests of genetically engineered crops producing both defensive allelochemicals and Bacillus thuringiensis (Bt) toxins often sequentially or simultaneously uptake allelochemicals, Bt toxins, and/or insecticides. How the three types of toxins interact to kill pests remains underexplored. Here we investigated the interactions of Bt toxin Vip3A, plant allelochemical flavone, and insecticide emamectin benzoate in Spodoptera frugiperda. Simultaneous administration of flavone LC25 + Vip3A LC25, emamectin benzoate LC25 + Vip3A LC25, and flavone LC15 + emamectin benzoate LC15 + Vip3A LC15 but not flavone LC25 + emamectin LC25 yielded a mortality significantly higher than their expected additive mortality (EAM). One-day pre-exposure to one toxin at LC5 followed by six-day exposure to the same toxin at LC5 plus another toxin at LC50 showed that the mortality of flavone LC5 + Vip3A LC50, emamectin benzoate LC5 + Vip3A LC50, and Vip3A LC5 + emamectin benzoate LC50 were significantly higher than their EAM, while that of flavone LC5 + emamectin benzoate LC50 was significantly lower than their EAM. No significant difference existed among the mortalities of Vip3A LC5 + flavone LC50, emamectin benzoate LC5 + flavone LC50, and their EAMs. The results suggest that the interactions of the three toxins are largely synergistic (inductive) or additive, depending on their combinations and doses.

Midgut serine proteinases participate in dietary adaptations of the castor (Eri) silkworm Samia ricini Anderson transferred from Ricinus communis to an ancestral host, Ailanthus excelsa Roxb

Dietary change influenced the life-history traits, nutritional utilization, and midgut serine proteinases in the larvae of the domesticated polyphagous S. ricini, transferred from R. communis (common name: castor; family Euphorbiaceae; the host plant implicated in its domestication) to A. excelsa (common name: Indian tree of heaven; family Simaroubaceae; an ancestral host of wild Samia species). Significantly higher values for fecundity and body weight were observed in larvae feeding on R. communis (Scr diet), and they took less time to reach pupation than insects feeding on A. excelsa (Scai diet). Nevertheless, the nutritional index for efficiency of conversion of digested matter (ECD) was similar for larvae feeding on the two plant species, suggesting the physiological adaptation of S. ricini (especially older instars) to an A. excelsa diet. In vitro protease assays and gelatinolytic zymograms using diagnostic substrates and protease inhibitors revealed significantly elevated levels (p ≤ 0.05) of digestive trypsins, which may be associated with the metabolic costs influencing slow growth in larvae feeding on A. excelsa. RT-PCR with semidegenerate serine proteinase gene-specific primers, and cloning and sequencing of 3' cDNA ends identified a large gene family comprising at least two groups of putative chymotrypsins (i.e., Sr I and Sr II) resembling invertebrate brachyurins/collagenases with wide substrate specificities, and five groups of putative trypsins (i.e., Sr III, Sr IV, Sr V, Sr VII, and Sr VIII). Quantitative RT-PCR indicated that transcripts belonging to the Sr I, Sr III, Sr IV, and Sr V groups, especially the Sr IV group (resembling achelase I from Lonomia achelous), were expressed differentially in the midguts of fourth instars reared on the two plant species. Sequence similarity indicated shared lineages with lepidopteran orthologs associated with expression in the gut, protein digestion, and phytophagy. The results obtained are discussed in the context of larval serine proteinases in dietary adaptations, domestication, and exploration of new host plant species for commercial rearing of S. ricini.

Involvement of an Enhanced Immunity Mechanism in the Resistance to Bacillus thuringiensis in Lepidopteran Pests

Bacillus thuringiensis (Bt) is the safest, economically successful entomopathogen to date. It is extensively produced in transgenic crops or used in spray formulations to control Lepidopteran pests. The most serious threat to the sustainable usage of Bt is insect resistance. The resistance mechanisms to Bt toxins depend not only on alterations in insect receptors, but also on the enhancement of insect immune responses. In this work, we review the current knowledge of the immune response and resistance of insects to Bt formulations and Bt proteins, mainly in Lepidopteran pests. We discuss the pattern recognition proteins for recognizing Bt, antimicrobial peptides (AMPs) and their synthetic signaling pathways, the prophenoloxidase system, reactive oxygen species (ROS) generation, nodulation, encapsulation, phagocytosis, and cell-free aggregates, which are involved in immune response reactions or resistance to Bt. This review also analyzes immune priming, which contributes to the evolution of insect resistance to Bt, and puts forward strategies to improve the insecticidal activity of Bt formulations and manage insect resistance, targeting the insect immune responses and resistance.

Dynamics of Zea mays transcriptome in response to a polyphagous herbivore, Spodoptera litura

Zea mays defense response is well-crafted according to the physical and chemical weapons utilized by their invaders during the coevolutionary period. Maize plants employ diversified defense strategies and alter the spatiotemporal distribution of several classes of defensive compounds to affect insect herbivore performance. However, only little knowledge is available about the defense orchestration of maize in response to Spodoptera litura, a voracious Noctuidae pest. In order to decipher the defense status of Zea mays (African tall variety) against S. litura, a comparative feeding bioassay was executed, which revealed reduced performance of the herbivore on maize. In order to understand the molecular mechanism behind maize tolerance against S. litura, a microarray-based genome-wide expression analysis was performed. The comparative analysis displayed 792 differentially expressed genes (DEGs), wherein 357 genes were upregulated and 435 genes were downregulated at fold change ≥ 2 and p value ≤ 0.05. The upregulated genes were identified and categorized as defense-related, oxidative stress-related, transcription regulatory genes, protein synthesis genes, phytohormone-related, and primary and secondary metabolism-related. In contrast, downregulated genes were mainly associated with plant growth and development, indicating a balance of growth and defense response and utilization of a highly evolved C-diversion response were noticed. Maize plants showed better tolerance against herbivory and maintained its fitness using a combinatorial strategy. This peculiar response of Zea mays against S. litura offers an excellent possibility of managing polyphagous pests by spicing up the plant's defensive response with tolerance mechanism.

Vegetative Insecticidal Protein (Vip): A Potential Contender From Bacillus thuringiensis for Efficient Management of Various Detrimental Agricultural Pests

Bacillus thuringiensis (Bt) bacterium is found in various ecological habitats, and has natural entomo-pesticidal properties, due to the production of crystalline and soluble proteins during different growth phases. In addition to Cry and Cyt proteins, this bacterium also produces Vegetative insecticidal protein (Vip) during its vegetative growth phase, which is considered an excellent toxic candidate because of the difference in sequence homology and receptor sites from Cry proteins. Vip proteins are referred as second-generation insecticidal proteins, which can be used either alone or in complementarity with Cry proteins for the management of various detrimental pests. Among these Vip proteins, Vip1 and Vip2 act as binary toxins and have toxicity toward pests belonging to Hemiptera and Coleoptera orders, whereas the most important Vip3 proteins have insecticidal activity against Lepidopteran pests. These Vip3 proteins are similar to Cry proteins in terms of toxicity potential against susceptible insects. They are reported to be toxic toward pests, which can’t be controlled with Cry proteins. The Vip3 proteins have been successfully pyramided along with Cry proteins in transgenic rice, corn, and cotton to combat resistant pest populations. This review provides detailed information about the history and importance of Vip proteins, their types, structure, newly identified specific receptors, and action mechanism of this specific class of proteins. Various studies conducted on Vip proteins all over the world and the current status have been discussed. This review will give insights into the significance of Vip proteins as alternative promising candidate toxic proteins from Bt for the management of pests in most sustainable manner.

Response Mechanisms of Invertebrates to Bacillus thuringiensis and Its Pesticidal Proteins

Extensive use of chemical insecticides adversely affects both environment and human health. One of the most popular biological pest control alternatives is bioinsecticides based on Bacillus thuringiensis This entomopathogenic bacterium produces different protein types which are toxic to several insect, mite, and nematode species. Currently, insecticidal proteins belonging to the Cry and Vip3 groups are widely used to control insect pests both in formulated sprays and in transgenic crops. However, the benefits of B. thuringiensis-based products are threatened by insect resistance evolution. Numerous studies have highlighted that mutations in genes coding for surrogate receptors are responsible for conferring resistance to B. thuringiensis Nevertheless, other mechanisms may also contribute to the reduction of the effectiveness of B. thuringiensis-based products for managing insect pests and even to the acquisition of resistance. Here, we review the relevant literature reporting how invertebrates (mainly insects and Caenorhabditis elegans) respond to exposure to B. thuringiensis as either whole bacteria, spores, and/or its pesticidal proteins.

Dynamic changes of transcriptome of fifth-instar spodoptera litura larvae in response to insecticide

Spodoptera litura is a major insect with a cosmopolitan distribution and strong resistance to multiple insecticides. Determining the molecular basis and key candidate genes of the insecticide resistance of S. litura may help in managing this insect. In this study, fifth-instar S. litura larvae were subjected to transcriptome analysis at 6, 12, 24, 48, and 72 h after feeding on an LC20 dose of avermectin. The result showed that genes responding to avermectin changed dynamically with different gene counts and resistance mechanisms at the fifth instar based on a metabolic pathway map. These responses included degrading the insecticide by a series of P450 and glutathione-S-transferase enzymes starting at the 12 h time point, with subsequent increases in the number of genes involved and shifts to TOLL and immune deficiency (IMD) pathways at 48 h after feeding the insecticide. Weighted correlation network analysis (WGCNA) determined a co-expression module related to the avermectin response at 12 and 24 h (r = 0.403, p = 0.0371; r = 0.436, p = 0.023), in which a hub gene (LOC111358940) related to metalloproteinase activity was identified. In addition, Analysis of the genes in the co-expression module further revealed that eight genes encoding UDP-glucuronosyltransferases were directly associated with insecticide response in S. litura. These results provide better understanding of the avermectin response mechanism of S. litura and may be useful in developing improved control strategies for this species.

SUPPLEMENTARY INFORMATION

The online version of this article (10.1007/s13205-021-02651-9) contains supplementary material, which is available to authorized users.

Autophagy induced by Vip3Aa has a pro-survival role in Spodoptera frugiperda Sf9 cells

Vip3Aa is an insecticidal protein that can effectively control certain lepidopteran pests and has been used widely in biological control. However, the mechanism of action of Vip3Aa is unclear. In the present study, we showed that Vip3Aa could cause autophagy in Sf9 cells, which was confirmed by the increased numbers of GFP-Atg8 puncta, the appearance of autophagic vacuoles, and an elevated Atg8-II protein level. Moreover, we found that the AMPK-mTOR-ULK1 pathway is involved in Vip3Aa-induced autophagy, which might be associated with the destruction of ATP homeostasis in Vip3Aa-treated cells. Both the elevated p62 level and the increased numbers of GFP-RFP-Atg8 yellow fluorescent spots demonstrated that autophagy in Sf9 cells was inhibited at 24 h after Vip3Aa treatment. With the prolongation of Vip3Aa treatment time, this inhibition became more serious and led to autophagosome accumulation. Genetic knockdown of ATG5 or the use of the autophagy inhibitor 3-MA further increased the sensitivity of Sf9 cells to Vip3Aa. Overexpression of ATG5 reduced the cell mortality of Vip3Aa-treated cells. In summary, the results revealed that autophagy induced by Vip3Aa has a pro-survival role, which might be related to the development of insect resistance.

In silico identification of effector proteins from generalist herbivore Spodoptera litura

BACKGROUND

The common cutworm, Spodoptera litura Fabricius is a leaf and fruit feeding generalist insect of the order Lepidoptera and a destructive agriculture pest. The broad host range of the herbivore is due to its ability to downregulate plant defense across different plants. The identity of Spodoptera litura released effectors that downregulate plant defense are largely unknown. The current study aims to identify genes encoding effector proteins from salivary glands of S. litura (Fab.).

RESULTS

Head and salivary glands of Spodoptera litura were used for de-novo transcriptome analysis and effector prediction. Eight hundred ninety-nine proteins from the head and 330 from salivary gland were identified as secretory proteins. Eight hundred eight proteins from the head and 267 from salivary gland proteins were predicted to be potential effector proteins.

CONCLUSIONS

This study is the first report on identification of potential effectors from Spodoptera litura salivary glands.

Current Insights on Vegetative Insecticidal Proteins (Vip) as Next Generation Pest Killers

Bacillus thuringiensis (Bt) is a Gram negative soil bacterium. This bacterium secretes various proteins during different growth phases with an insecticidal potential against many economically important crop pests. One of the important families of Bt proteins is vegetative insecticidal proteins (Vip), which are secreted into the growth medium during vegetative growth. There are three subfamilies of Vip proteins. Vip1 and Vip2 heterodimer toxins have an insecticidal activity against many Coleopteran and Hemipteran pests. Vip3, the most extensively studied family of Vip toxins, is effective against Lepidopteron. Vip proteins do not share homology in sequence and binding sites with Cry proteins, but share similarities at some points in their mechanism of action. Vip3 proteins are expressed as pyramids alongside Cry proteins in crops like maize and cotton, so as to control resistant pests and delay the evolution of resistance. Biotechnological- and in silico-based analyses are promising for the generation of mutant Vip proteins with an enhanced insecticidal activity and broader spectrum of target insects.

The Tripartite Interaction of Host Immunity-Bacillus thuringiensis Infection-Gut Microbiota

Bacillus thuringiensis (Bt) is an important cosmopolitan bacterial entomopathogen, which produces various protein toxins that have been expressed in transgenic crops. The evolved molecular interaction between the insect immune system and gut microbiota is changed during the Bt infection process. The host immune response, such as the expression of induced antimicrobial peptides (AMPs), the melanization response, and the production of reactive oxygen species (ROS), varies with different doses of Bt infection. Moreover, B. thuringiensis infection changes the abundance and structural composition of the intestinal bacteria community. The activated immune response, together with dysbiosis of the gut microbiota, also has an important effect on Bt pathogenicity and insect resistance to Bt. In this review, we attempt to clarify this tripartite interaction of host immunity, Bt infection, and gut microbiota, especially the important role of key immune regulators and symbiotic bacteria in the Bt killing activity. Increasing the effectiveness of biocontrol agents by interfering with insect resistance and controlling symbiotic bacteria can be important steps for the successful application of microbial biopesticides.

Immune responses to Bacillus thuringiensis in the midgut of the diamondback moth, Plutella xylostella

The diamondback moth, Plutella xylostella, is the first insect to develop resistance to Bacillus thuringiensis (Bt) in the field. To date, little is known about the molecular mechanism of the interaction between Bt and midgut immunity in P. xylostella. Here, we report immune responses in the P. xylostella midgut to Bt strain Bt8010 using a combined approach of transcriptomics and quantitative proteomics. Many genes in the Toll, IMD, JNK and JAK-STAT pathways and antimicrobial peptide genes were activated at 18 h post-infection. In the prophenoloxidase (PPO) cascade, four serpin genes were activated, and the PPO1 gene was suppressed by Bt8010. Inhibition of the two PPO proteins was observed at 18 h post-infection. Feeding Bt8010-infected larvae recombinant PPOs enhanced their survival. These results revealed that the Toll, IMD, JNK and JAK-STAT pathways were triggered and participated in the immune defence of the midgut against Bt8010, while the PPO cascade was inhibited and played an important role in this process.

Activity of Vip3Aa1 against Periplaneta Americana

Bacillus thuringiensis (Bt) is a well-known entomopathogen. In this study, we cloned the vip3Aa1 gene from Bt strain GIM1.147 and investigated the insecticidal activity of Bt Vip3Aa1 protein produced by Escherichia coli against Periplaneta americana and Blattella germanica. The results showed that purified Vip3Aa1 exhibited an LC50 at 24 h against P. americana and B. germanica of 0.182 mg·ml-1 and 0.276 mg·ml-1, respectively. Investigations of its mode of action showed that Vip3Aa1 could be proteolyzed into a 62-kDa toxic protein by P. americana gut-soluble proteases. In addition, Vip3Aa1 caused severe damage to the columnar colon and the midgut, as observed through hematoxylin-eosin staining and scanning electron microscopy. The 62-kDa activated Vip3Aa1 protein could form ion channels in the colon and the midgut in vitro. Based on protease activity analysis, Vip3Aa1 at concentrations of 0.125 mg·ml-1 and 0.031 mg·ml-1 could downregulate the activities of glutathione S-transferase, α-NA esterase, trypsin, and chymotrypsin. This report provides the first description of the activity of Vip3Aa1 toxins toward P. americana and B. germanica and demonstrates that the mechanism through which Vip3Aa1 kills P. americana and B. germanica differs from that involved in the killing of lepidopteran insects.

Midgut de novo transcriptome analysis and gene expression profiling of Spodoptera exigua larvae exposed with sublethal concentrations of Cry1Ca protein

Spodoptera exigua (Hübner) is a polyphagous pest on agricultural crops, whose control is based mainly on the application of chemical insecticides. Bacillus thuringiensis (Bt) is one of the most important biological agents that have been successfully applied as a biological control, and Cry1Ca protein is considered to be active against S. exigua. Therefore, to understand the response of S. exigua to Cry1Ca protein, high-throughput sequencing was used to analyse the S. exigua larval midgut after treatment with sublethal concentrations of Cry1Ca protein. Transcriptome data showed that a total of 98,571 unigenes with an N50 value of 1135 bp and a mean length of 653 bp were obtained. Furthermore, 2962 differentially expressed genes (DEGs) were identified after Cry1Ca challenge, including 1508 up-regulated and 1454 down-regulated unigenes. Among these DEGs, detoxification (CYP, CarE, and GST) and Bt resistance (ALP, APN, and ABC transporter)-related genes were differentially expressed in the midgut of S. exigua after Cry1Ca treatment. However, most DEGs of protective enzymes were down-regulated, while most DEGs related with serine protease and REPAT were up-regulated. Furthermore, almost all DEGs related to the immune signaling pathway, antimicrobial protein, and lysozyme were up-regulated by Cry1Ca treatment. These results indicated that the detoxification enzyme, protective enzymes, Bt resistance-related genes, serine protease, REPAT, and the immune response might have been involved in the response of S. exigua to Cry1Ca protein. In summary, analysis of the transcriptomal expression of genes involved in Cry1Ca protein against S. exigua provided potential clues for elucidating the host response processes and defensive mechanisms underlying Cry1Ca toxicity.

Identification of transcriptome and fluralaner responsive genes in the common cutworm Spodoptera litura Fabricius, based on RNA-seq

Background

Fluralaner is a novel isoxazoline insecticide with a unique action site on the γ-aminobutyric acid receptor (GABAR), shows excellent activity on agricultural pests including the common cutworm Spodoptera litura, and significantly influences the development and fecundity of S. litura at either lethal or sublethal doses. Herein, Illumina HiSeq Xten (IHX) platform was used to explore the transcriptome of S. litura and to identify genes responding to fluralaner exposure.

Results

A total of 16,572 genes, including 451 newly identified genes, were observed in the S. litura transcriptome and annotated according to the COG, GO, KEGG and NR databases. These genes included 156 detoxification enzyme genes [107 cytochrome P450 enzymes (P450s), 30 glutathione S-transferases (GSTs) and 19 carboxylesterases (CarEs)] and 24 insecticide-targeted genes [5 ionotropic GABARs, 1 glutamate-gated chloride channel (GluCl), 2 voltage-gated sodium channels (VGSCs), 13 nicotinic acetylcholine receptors (nAChRs), 2 acetylcholinesterases (AChEs) and 1 ryanodine receptor (RyR)]. There were 3275 and 2491 differentially expressed genes (DEGs) in S. litura treated with LC₃₀ or LC₅₀ concentrations of fluralaner, respectively. Among the DEGs, 20 related to detoxification [16 P450s, 1 GST and 3 CarEs] and 5 were growth-related genes (1 chitin and 4 juvenile hormone synthesis genes). For 26 randomly selected DEGs, real-time quantitative PCR (RT-qPCR) results showed that the relative expression levels of genes encoding several P450s, GSTs, heat shock protein (HSP) 68, vacuolar protein sorting-associated protein 13 (VPSAP13), sodium-coupled monocarboxylate transporter 1 (SCMT1), pupal cuticle protein (PCP), protein takeout (PT) and low density lipoprotein receptor adapter protein 1-B (LDLRAP1-B) were significantly up-regulated. Conversely, genes encoding esterase, sulfotransferase 1C4, proton-coupled folate transporter, chitinase 10, gelsolin-related protein of 125 kDa (GRP), fibroin heavy chain (FHC), fatty acid synthase and some P450s were significantly down-regulated in response to fluralaner.

Conclusions

The transcriptome in this study provides more effective resources for the further study of S. litura whilst the DEGs identified sheds further light on the molecular response to fluralaner.

HT-SuperSAGE of the gut tissue of a Vip3Aa-resistant Heliothis virescens (Lepidoptera: Noctuidae) strain provides insights into the basis of resistance

Multitoxin Bt-crops expressing insecticidal toxins with different modes of action, for example, Cry and Vip, are expected to improve resistance management in target pests. While Cry1A resistance has been relatively well characterized in some insect species, this is not the case for Vip3A, for which no mechanism of resistance has yet been identified. Here we applied HT-SuperSAGE to analyze the transcriptome of the gut tissue of tobacco budworm Heliothis virescens (F.) laboratory-selected for Vip3Aa resistance. From a total of 1 324 252 sequence reads, 5 895 126-bp tags were obtained representing 17 751 nonsingleton unique transcripts (UniTags) from genetically similar Vip3Aa-resistant (Vip-Sel) and susceptible control (Vip-Unsel) strains. Differential expression was significant (≥2.5 fold or ≤0.4; P < 0.05) for 1989 sequences (11.2% of total UniTags), where 420 represented overexpressed (OE) and 1569 underexpressed (UE) genes in Vip-Sel. BLASTN searches mapped 419 UniTags to H. virescens sequence contigs, of which, 416 (106 OE and 310 UE) were unambiguously annotated to proteins in NCBI nonredundant protein databases. Gene Ontology distributed 345 of annotated UniTags in 14 functional categories with metabolism (including serine-type hydrolases) and translation/ribosome biogenesis being the most prevalent. A UniTag homologous to a particular member of the REsponse to PAThogen (REPAT) family was found among most overexpressed, while UniTags related to the putative Vip3Aa-binding ribosomal protein S2 (RpS2) were underexpressed. qRT-PCR of a subset of UniTags validated the HT-SuperSAGE data. This study is the first providing lepidopteran gut transcriptome associated with Vip3Aa resistance and a foundation for future attempts to elucidate the resistance mechanism.

Bacillus thuringiensis Suppresses the Humoral Immune System to Overcome Defense Mechanism of Plutella xylostella

Background: Plutella xylostella has become a notorious pest of cruciferous crops all over the world. Delta-endotoxins of Bacillus thuringiensis are widely used insecticidal proteins for controlling P. xylostella. However, the interaction mechanism of B. thuringiensis with the immune system of P. xylostella, at the genomic level, is still unclear. This study explored the immune response of P. xylostella to B. thuringiensis, at different time intervals, 6 h, 12 h, 18 h, 24 h, and 36 h, by using RNA-Sequencing (RNA-Seq) and RT-qPCR. Results: In total, 167 immunity-related genes were identified and placed into different families, including pattern recognition receptors (PRRs), signal modulators, immune pathways (Toll, IMD, and JAK/STAT), and immune effectors. It is worth mentioning that the analyses of the differentially expressed immunity-related genes revealed that most of the differentially expressed genes (DEGs) (87, 56, 76, 67, and 73 genes) were downregulated in P. xylostella following B. thuringiensis oral infection at 6 h, 12 h, 18 h, 24 h, and 36 h. Interestingly, our RNA-Seq analysis also revealed reduced expression of antimicrobial peptides, that play a vital role in the humoral immune system of P. xylostella. Conclusion: This study demonstrates that B. thuringiensis plays a novel role in controlling P. xylostella, by suppressing the immune system.

Artefactual band patterns by SDS-PAGE of the Vip3Af protein in the presence of proteases mask the extremely high stability of this protein

Vip3 proteins are secretable proteins from Bacillus thuringiensis with important characteristics for the microbiological control of agricultural pests. The exact details of their mode of action are yet to be disclosed and the crystallographic structure is still unknown. Vip3 proteins are expressed as protoxins that have to be activated by the insect gut proteases. A previous study on the peptidase processing of Vip3Aa revealed that the protoxin produced artefactual band patterns by SDS-PAGE due to the differential stability of this protein and the peptidases to SDS and heating (Bel et al., 2017 Toxins 9:131). To determine whether this phenomenon also applies to other Vip3A proteins, here we chose a different Vip3A protein (Vip3Af) and subjected it to commercial trypsin and midgut juice from a target insect species (Spodoptera frugiperda). The misleading degradation patterns were also observed with Vip3Af, both with trypsin and midgut juice. However, gel filtration chromatography showed that, under native conditions, Vip3Af is found as a tetramer and that peptidases only act upon primary cleavage sites. The proteolytic cleavage renders two fragments of approximately 20 kDa and 65 kDa which remain together in the tretameric structure and that are no further processed even at high peptidase concentrations.

Midgut Transcriptional Variation of Chilo suppressalis Larvae Induced by Feeding on the Dead-End Trap Plant, Vetiveria zizanioides

Chilo supprressalis is one of the most important rice pests that causes serious damage to production in the rice growth area of Asia. Vetiver grass (Vetiveria zizanioides) was previously found to effectively attract female adults of C. suppressalis laying eggs on vetiver leaves, while the larvae cannot complete their life cycles by feeding on vetiver, indicating a potential means of controlling this pest. In the present study, the transcriptomes of midguts of rice-fed and vetiver-fed C. suppressalis larvae were profiled, which aimed to clarify the molecular mechanism of vetiver as a dead-end trap plant preliminarily. We found that ingestion of vetiver provoked a robust transcriptional response in the larval midguts, and a total of 1,849 differentially expressed UniGenes were identified. We focused on 12 digestion-related genes, four immune-related genes and three detoxification-related genes. Most of these genes were significantly down regulated in the larval midguts at 6, 8, and 10 days after feeding on vetiver compared to on rice. Transcriptional dynamics suggested that these genes might be involved in toxicity responses following exposure to vetiver. Taken together, this study provides an initial molecular framework for developing biological control strategies for C. suppressalis in an effort to protect economically important rice crops.

Draft genome and reference transcriptomic resources for the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

The pine processionary moth Thaumetopoea pityocampa (Lepidoptera: Notodontidae) is the main pine defoliator in the Mediterranean region. Its urticating larvae cause severe human and animal health concerns in the invaded areas. This species shows a high phenotypic variability for various traits, such as phenology, fecundity and tolerance to extreme temperatures. This study presents the construction and analysis of extensive genomic and transcriptomic resources, which are an obligate prerequisite to understand their underlying genetic architecture. Using a well-studied population from Portugal with peculiar phenological characteristics, the karyotype was first determined and a first draft genome of 537 Mb total length was assembled into 68,292 scaffolds (N50 = 164 kb). From this genome assembly, 29,415 coding genes were predicted. To circumvent some limitations for fine-scale physical mapping of genomic regions of interest, a 3X coverage BAC library was also developed. In particular, 11 BACs from this library were individually sequenced to assess the assembly quality. Additionally, de novo transcriptomic resources were generated from various developmental stages sequenced with HiSeq and MiSeq Illumina technologies. The reads were de novo assembled into 62,376 and 63,175 transcripts, respectively. Then, a robust subset of the genome-predicted coding genes, the de novo transcriptome assemblies and previously published 454/Sanger data were clustered to obtain a high-quality and comprehensive reference transcriptome consisting of 29,701 bona fide unigenes. These sequences covered 99% of the cegma and 88% of the busco highly conserved eukaryotic genes and 84% of the busco arthropod gene set. Moreover, 90% of these transcripts could be localized on the draft genome. The described information is available via a genome annotation portal (http://bipaa.genouest.org/sp/thaumetopoea_pityocampa/).

Changes in gene expression and apoptotic response in Spodoptera exigua larvae exposed to sublethal concentrations of Vip3 insecticidal proteins

The insecticidal Vip3 proteins from Bacillus thuringiensis (Bt), along with the classical Bt Cry proteins, are currently used in Bt-crops to control insect pests, since they do not share the same mode of action. Here we characterized the response of Spodoptera exigua larvae after Vip3 challenge. The expression profile of 47 genes was analyzed in larvae challenged with three concentrations of Vip3Ca. Results showed that the up-regulated genes were mainly involved in immune response, whereas the down-regulated genes were mainly involved in the digestion process. Other mechanisms of cellular response to the damage such as apoptosis were analyzed. For this analysis, sections from the midguts were examined by terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) staining. The nuclei of the midgut epithelial cells were stained at the highest concentration of the Vip3Ca protein and at lower concentrations of Vip3Aa in agreement with the different potency of the two proteins. In addition, apoptosis was also examined by the analysis of the expression of five caspase genes. The present study shows that exposure of S. exigua larvae to sublethal concentrations of Vip3 proteins activates different insect response pathways which trigger the regulation of some genes, APN shedding, and apoptotic cell death.

Insights into the Structure of the Vip3Aa Insecticidal Protein by Protease Digestion Analysis

Vip3 proteins are secretable proteins from Bacillus thuringiensis whose mode of action is still poorly understood. In this study, the activation process for Vip3 proteins was closely examined in order to better understand the Vip3Aa protein stability and to shed light on its structure. The Vip3Aa protoxin (of 89 kDa) was treated with trypsin at concentrations from 1:100 to 120:100 (trypsin:Vip3A, w:w). If the action of trypsin was not properly neutralized, the results of SDS-PAGE analysis (as well as those with Agrotis ipsilon midgut juice) equivocally indicated that the protoxin could be completely processed. However, when the proteolytic reaction was efficiently stopped, it was revealed that the protoxin was only cleaved at a primary cleavage site, regardless of the amount of trypsin used. The 66 kDa and the 19 kDa peptides generated by the proteases co-eluted after gel filtration chromatography, indicating that they remain together after cleavage. The 66 kDa fragment was found to be extremely resistant to proteases. The trypsin treatment of the protoxin in the presence of SDS revealed the presence of secondary cleavage sites at S-509, and presumably at T-466 and V-372, rendering C-terminal fragments of approximately 29, 32, and 42 kDa, respectively. The fact that the predicted secondary structure of the Vip3Aa protein shows a cluster of beta sheets in the C-terminal region of the protein might be the reason behind the higher stability to proteases compared to the rest of the protein, which is mainly composed of alpha helices.

Knockdown of the MAPK p38 pathway increases the susceptibility of Chilo suppressalis larvae to Bacillus thuringiensis Cry1Ca toxin

The bacterium Bacillus thuringiensis (Bt) produces a wide range of toxins that are effective against a number of insect pests. Identifying the mechanisms responsible for resistance to Bt toxin will improve both our ability to control important insect pests and our understanding of bacterial toxicology. In this study, we investigated the role of MAPK pathways in resistance against Cry1Ca toxin in Chilo suppressalis, an important lepidopteran pest of rice crops. We first cloned the full-length of C. suppressalis mitogen-activated protein kinase (MAPK) p38, ERK1, and ERK2, and a partial sequence of JNK (hereafter Csp38, CsERK1, CsERK2 and CsJNK). We could then measure the up-regulation of these MAPK genes in larvae at different times after ingestion of Cry1Ca toxin. Using RNA interference to knockdown Csp38, CsJNK, CsERK1 and CsERK2 showed that only knockdown of Csp38 significantly increased the mortality of larvae to Cry1Ca toxin ingested in either an artificial diet, or after feeding on transgenic rice expressed Cry1Ca. These results suggest that MAPK p38 is responsible for the resistance of C. suppressalis larvae to Bt Cry1Ca toxin.