Changes in gene expression in the larval gut of Ostrinia nubilalis in Response to Bacillus thuringiensis Cry1Ab protoxin ingestion

Culex quinquefasciatus Resistant to the Binary Toxin from Lysinibacillus sphaericus Displays a Consistent Downregulation of Pantetheinase Transcripts

Culex quinquefasciatus resistance to the binary (Bin) toxin, the major larvicidal component from Lysinibacillus sphaericus, is associated with mutations in the cqm1 gene, encoding the Bin-toxin receptor. Downregulation of the cqm1 transcript was found in the transcriptome of larvae resistant to the L. sphaericus IAB59 strain, which produces both the Bin toxin and a second binary toxin, Cry48Aa/Cry49Aa. Here, we investigated the transcription profiles of two other mosquito colonies having Bin resistance only. These confirmed the cqm1 downregulation and identified transcripts encoding the enzyme pantetheinase as the most downregulated mRNAs in both resistant colonies. Further quantification of these transcripts reinforced their strong downregulation in Bin-resistant larvae. Multiple genes were found encoding this enzyme in Cx. quinquefasciatus and a recombinant pantetheinase was then expressed in Escherichia coli and Sf9 cells, with its presence assessed in the midgut brush border membrane of susceptible larvae. The pantetheinase was expressed as a ~70 kDa protein, potentially membrane-bound, which does not seem to be significantly targeted by glycosylation. This is the first pantetheinase characterization in mosquitoes, and its remarkable downregulation might reflect features impacted by co-selection with the Bin-resistant phenotype or potential roles in the Bin-toxin mode of action that deserve to be investigated.

Transcriptome profiling analysis of the intoxication response in midgut tissue of Agrotis ipsilon larvae to Bacillus thuringiensis Vip3Aa protoxin

Vip insecticidal proteins are produced by Bacillus thuringiensis (Bt) during its vegetative growth phase. In the present study, Vip3Aa11 and Vip3Aa39 proteins were investigated. These two proteins present 39 amino acid differential sites and they shared 95.06% amino acid sequence similarity. They are effective against some Lepidoptera insect larvae. In a previous study, using artificial diet bioassays, we estimated the LC₅₀ of Vip3Aa11 and Vip3Aa39 strains against Agrotis ipsilon larvae were 73.41 μg/mL (with 95% confidence interval of 2.34-11.19) and 5.43 μg/mL (with 95% confidence interval of 43.20-115.03), respectively. To investigate the response of Agrotis ipsilon transcriptome in defending against Vip3Aa11 and Vip3Aa39 toxins, we performed high-throughput RNA-sequencing on cDNA generated from the midguts of Agrotis ipsilon larvae that consumed a control diet (CK-M-A), Vip3Aa11 (Vip3Aa11-M-A) and Vip3Aa39 (Vip3Aa39-M-A) proteins. We generated about 98.87 Gb bases in total on BGISEQ-500 sequencing platform. After assembling all samples together and filtering the abundance, we got 51,887 unigenes, the total length, average length, N50 and GC content of unigenes are 64,523,651 bp, 1243 bp, 2330 bp and 41.81% respectively. We revealed 558 midgut genes differential expressed in Vip3Aa11-M-A and 65 midgut genes differentially expressed in Vip3Aa39-M-A. The differentially expressed genes were enriched for serine proteases and potential Bt Vip toxin midgut receptor genes. Eleven serine proteases related genes and 13 Bt toxin potential receptor genes with differential expression were found. Based on transcriptome profiling, we focused on validation the sensitivity of these two Vip3Aa proteins to trypsin and their binding properties to Agrotis ipsilon midgut BBMV (Brush Border Membrane Vesicles). The results show that the sensitivity of the two proteins to trypsin is similar. Binding experiments revealed that both proteins can bind to Agrotis ipsilon midgut BBMV, and there is a competitive binding between them. This transcriptome dataset provided a comprehensive sequence resource of Agrotis ipsilon and provides a foundation for comparative studies with other species of insects.

Role of the peritrophic matrix in insect-pathogen interactions

The peritrophic matrix (PM) is an acellular chitin and glycoprotein layer that lines the invertebrate midgut. The PM has long been considered a physical as well as a biochemical barrier, protecting the midgut epithelium from abrasive food particles, digestive enzymes and pathogens infectious per os. This short review will focus on the latter function, as a barrier to pathogens infectious per os. We focus on the evidence confirming the role of the PM as protective barrier against pathogenic microorganisms of insects, mainly bacteria and viruses, as well as the evolution of a variety of mechanisms used by pathogens to overcome the PM barrier.

Impact of Antibiotics on Efficacy of Cry Toxins Produced in Two Different Genetically Modified Bt Maize Varieties in Two Lepidopteran Herbivore Species, Ostrinia nubilalis and Spodoptera littoralis

The insecticidal crystal proteins from Bacillus thuringiensis (Bt) are widely-used biopesticides that are used both as Bt spore-crystal preparations in sprayable formulations and as activated toxins in genetically modified (GM) plants. Models for their modes of action have been proposed but many issues remain unresolved. Among those is the role of commensal gut bacteria in target insect death: previous studies showed that antibiotics attenuate the toxicity of Bt sprays. We tested whether antibiotics interfere with the effects of GM plant-produced Bt toxins in larvae of two Lepidopteran species, the European corn borer Ostrinia nubilalis and the cotton leafworm Spodoptera littoralis. The larvae were reared on artificial diet with or without antibiotics and, thereafter, fed two varieties of Bt GM maize in comparison to conventional non-Bt maize leaves sprayed with antibiotic solution and/or with a Bt formulation. Antibiotics significantly reduced or delayed the toxicity of Cry toxins, although to a lesser extent than previously reported for Bt-sprays. This supports the hypothesis that Cry toxins induce mortality by themselves in the absence of Bt bacteria and spores, and of commensal gut bacteria. However, larvae that were not treated with antibiotics died faster and at a higher rate which was further compounded by plant variety and species sensitivity. These findings support a hypothesis that toxicemia alone can inflict significant mortality. However, in the absence of antibiotics, the gut bacteria likely enhance the Cry toxin effect by inflicting, additionally, bacterial septicemia. This has important implications in field situations where antibiotic substances are present—e.g., from manure of animals from conventional production systems—and for ecotoxicological testing schemes of Bt toxins and nontarget organisms that are often using artificial diets enriched with high concentrations of antibiotics.

Carotenoids moderate the effectiveness of a Bt gene against the European corn borer, Ostrinia nubilalis

We assessed the effectiveness of a biofortified maize line (4BtxHC) which accumulates high levels of antioxidant carotenoids that also expressed the insecticidal Cry1Ac Bacillus thuringiensis (Bt) gene against the European corn borer Ostrinia nubilalis. This line had been previously engineered to accumulate carotenoids specifically in the seed endosperm, whereas the Bt gene was expressed constitutively. The concentrations of Bt toxin (Cry 1Ac) in the leaves of the 4Bt and 4BtxHC lines were not significantly different at 47±6 μg/g of fresh weight (FW); neither were they in the kernels of both lines (35±3 μg/g FW). The kernels and leaves were toxic to the larvae of O. nubilalis. However, the insecticidal activity was substantially lower (ca. 20%) than that of lines that expressed only Bt in spite that the two lines showed a quantity of toxin not significantly different in kernels or in leaves. Although the reduced effectiveness of Cry1Ac in kernels may not be entirely surprising, the observation of the same phenomenon in vegetative tissues was unexpected. When semi-artificial diets containing kernels from 4Bt supplemented with different levels of β-carotene were used in insect bioassays, the β-carotene moderated the effectiveness of the Bt similarly to the plant material with carotenoid enrichment. To elucidate the biochemical basis of the reduced effectiveness of Bt toxin in the carotenoid-enriched plants, we measured the activity of three enzymes known to be implicated in the detoxification defence, namely, catalase, superoxide dismutase and glutathione S-transferase. Whereas Cry1Ac expression significantly increased SOD and CAT enzymatic activity in the absence of carotenoids, carotenoids, either in 4BtxHC or in artificial diets enriched with β-carotene, significantly lowered CAT activity. Carotenoids can therefore moderate the susceptibility of the maize borer O. nubilalis to Cry1Ac, and we hypothesize that their role as antioxidants could explain this phenomenon via their scavenging of reactive oxygen species produced during Cry1Ac detoxification in the larvae. The involvement of this mechanism in the decreased mortality caused by Cry1Ac when carotenoids are present in the diet is discussed.

Comparison of gene expression profiles in the aquatic midge (Chironomus tentans) larvae exposed to two major agricultural pesticides

We developed a high-resolution expression microarray based on 2456 unique transcripts from a cDNA library of the aquatic midge (Chironomus tentans). By using the microarray, we detected that 146, 434 and 243 genes were differentially expressed after C. tentans larvae were exposed to chlorpyrifos (organophosphate insecticide) at 0.1 and 0.5 μg/L, and atrazine (triazine herbicide) at 1000 μg/L, respectively, for 48 h. The number of differentially expressed genes in the larvae exposed to chlorpyrifos at 0.5 μg/L was three times of that in the larvae exposed to chlorpyrifos at 0.1 μg/L. Among the differentially expressed genes in response to chlorpyrifos exposures, 76 genes showed significant Blast hits, and among them 42 were in common between the chlorpyrifos and atrazine exposures. In 19 differentially expressed xenobiotic detoxification genes, 16 were significantly up-regulated in the larvae exposed to chlorpyrifos and/or atrazine. Two cytochrome P450 genes (CtCYP6EV1 and CtCYP4DG2) were specifically up-regulated by chlorpyrifos, whereas three cytochrome P450 genes (CtCYP4DG1, CtCYP6EX3 and CtCYP6EV3) were specifically up-regulated by atrazine. Our results showed that chlorpyrifos exposures even at low concentrations can lead to significant changes in gene expression. The significant transcriptional responses are likely attributed to larval intoxication by the insecticide. These results not only support our previous studies in which candidate gene approaches were used, but also can potentially help develop specific molecular markers for monitoring pesticide exposures in non-target organisms in aquatic systems.

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.

FOXA transcriptional factor modulates insect susceptibility to Bacillus thuringiensis Cry1Ac toxin by regulating the expression of toxin-receptor ABCC2 and ABCC3 genes

Cry toxins produced by Bacillus thuringiensis (Bt) are insecticidal proteins widely used in insect control. Recently, it was shown that ATP-binding cassette transporter proteins (ABC) such as ABCC2, ABCC3, ABCG1 and ABCA2 are implicated in the insecticidal action of Cry toxins as putative receptors. However, the transcriptional regulators involved in the expression of ABC transporter genes remain unknown. Sequence analysis of promoter regions of ABCC2 gene from Helicoverpa armigera and ABCC3 gene from Spodoptera litura Sl-HP cultured cells, revealed the potential participation of Forkhead box protein A (FOXA), a transcription factor that regulates the expression of genes through remodeling chromatin. To determine if FOXA was involved in regulating expression of ABCC2 and ABCC3 genes, the expression of FOXA, ABCC2 and ABCC3 was compared in Sl-HP cells that are sensitive to Cry1Ac toxin with those in S. frugiperda Sf9 cells that are not sensitive to the toxin. Expression levels of those genes were significantly higher in Sl-HP than in Sf9 cells. Transient expression of FOXA in Sf9 cells activated ABCC2 and ABCC3 transcription, which directly correlated with enhanced Cry1Ac-susceptibility in these cells. Silencing of FOXA gene expression by RNAi in H. armigera larvae resulted in a decreased expression of ABCC2 and ABCC3 without affecting expression of other Cry toxin receptor genes such as alkaline phosphatase, aminopeptidase or cadherin. Silencing of FOXA gene expression also resulted in a Cry1Ac-tolerant phenotype since lower mortality and higher pupation rate were observed in diet containing Cry1Ac protoxin in comparison with the control group. These results demonstrate that FOXA up-regulates expression of the Cry1Ac-toxin receptor ABCC2 and ABCC3 genes, and that lower FOXA expression correlates with tolerance to Cry toxin in cell lines and in lepidopteran larvae.

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.

Comparisons of Transcriptional Profiles of Gut Genes between Cry1Ab-Resistant and Susceptible Strains of Ostrinia nubilalis Revealed Genes Possibly Related to the Adaptation of Resistant Larvae to Transgenic Cry1Ab Corn

A microarray developed on the basis of 2895 unique transcripts from larval gut was used to compare gut gene expression profiles between a laboratory-selected Cry1Ab-resistant (R) strain and its isoline susceptible (S) strain of the European corn borer (Ostrinia nubilalis) after the larvae were fed the leaves of transgenic corn (MON810) expressing Cry1Ab or its non-transgenic isoline for 6 h. We revealed 398 gut genes differentially expressed (i.e., either up- or down-regulated genes with expression ratio ≥2.0) in S-strain, but only 264 gut genes differentially expressed in R-strain after being fed transgenic corn leaves. Although the percentages of down-regulated genes among the total number of differentially expressed genes (50% in S-strain and 45% in R-strain) were similar between the R- and S-strains, the expression ratios of down-regulated genes were much higher in S-strain than in R-strain. We revealed that 17 and 9 significantly up- or down-regulated gut genes from S and R-strain, respectively, including serine proteases and aminopeptidases. These genes may be associated with Cry1Ab toxicity by degradation, binding, and cellular defense. Overall, our study suggests enhanced adaptation of Cry1Ab-resistant larvae on transgenic Cry1Ab corn as revealed by lower number and lower ratios of differentially expressed genes in R-strain than in S-strain of O. nubilalis.

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

Vip proteins, a new group of insecticidal toxins produced by Bacillus thuringiensis, are effective against specific pests including Spodoptera litura. Here, we report construction of a transcriptome database of S. litura by de novo assembly along with detection of the transcriptional response of S. litura larvae to Vip3Aa toxin. In total, 56,498 unigenes with an N50 value of 1,853 bp were obtained. Results of transcriptome abundance showed that Vip3Aa toxin provoked a wide transcriptional response of the S. litura midgut. The differentially expressed genes were enriched for immunity-related, metabolic-related and Bt-related genes. Twenty-nine immunity-related genes, 102 metabolic-related genes and 62 Bt-related genes with differential expression were found. On the basis of transcriptional profiling analysis, we focus on the functional validation of trypsin which potentially participated in the activation of Vip3Aa protoxin. Zymogram analysis indicated that the presence of many proteases, including trypsin, in S. litura larvae midgut. Results of enzymolysis in vitro of Vip3Aa by trypsin, and bioassay and histopathology of the trypsin-digested Vip3Aa toxin showed that trypsin was possibly involved in the Vip3Aa activation. This study provides a transcriptome foundation for the identification and functional validation of the differentially expressed genes in an agricultural important pest, S. litura.

Genome-wide identification and expression profiling of serine proteases and homologs in the diamondback moth, Plutella xylostella (L.)

Background

Serine proteases (SPs) are crucial proteolytic enzymes responsible for digestion and other processes including signal transduction and immune responses in insects. Serine protease homologs (SPHs) lack catalytic activity but are involved in innate immunity. This study presents a genome-wide investigation of SPs and SPHs in the diamondback moth, Plutella xylostella (L.), a globally-distributed destructive pest of cruciferous crops.

Results

A total of 120 putative SPs and 101 putative SPHs were identified in the P. xylostella genome by bioinformatics analysis. Based on the features of trypsin, 38 SPs were putatively designated as trypsin genes. The distribution, transcription orientation, exon-intron structure and sequence alignments suggested that the majority of trypsin genes evolved from tandem duplications. Among the 221 SP/SPH genes, ten SP and three SPH genes with one or more clip domains were predicted and designated as PxCLIPs. Phylogenetic analysis of CLIPs in P. xylostella, two other Lepidoptera species (Bombyx mori and Manduca sexta), and two more distantly related insects (Drosophila melanogaster and Apis mellifera) showed that seven of the 13 PxCLIPs were clustered with homologs of the Lepidoptera rather than other species. Expression profiling of the P. xylostella SP and SPH genes in different developmental stages and tissues showed diverse expression patterns, suggesting high functional diversity with roles in digestion and development.

Conclusions

This is the first genome-wide investigation on the SP and SPH genes in P. xylostella. The characterized features and profiled expression patterns of the P. xylostella SPs and SPHs suggest their involvement in digestion, development and immunity of this species. Our findings provide a foundation for further research on the functions of this gene family in P. xylostella, and a better understanding of its capacity to rapidly adapt to a wide range of environmental variables including host plants and insecticides.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2243-4) contains supplementary material, which is available to authorized users.

Novel genetic factors involved in resistance to Bacillus thuringiensis in Plutella xylostella

The widespread and sustainable exploitation of the entomopathogen Bacillus thuringiensis (Bt) in pest control is threatened by the evolution of resistance. Although resistance is often associated with loss of binding of the Bt toxins to the insect midgut cells, other factors have been implicated. Here we used suppressive subtractive hybridization and gene expression suppression to identify additional molecular components involved in Bt-resistance in Plutella xylostella. We isolated transcripts from genes that were differentially expressed in the midgut of larvae from a resistant population, following ingestion of a Bt kurstaki HD1 strain-based commercial formulation (DiPel), and compared with a genetically similar susceptible population. Quantitative real-time polymerase-chain reaction (RT-PCR) analysis confirmed the differential basal expression of a subset of these genes. Gene expression suppression of three of these genes (P. xylostella cyclin-dependent kinase 5 regulatory subunit associated protein 1-like 1, stromal cell-derived factor 2-like 1 and hatching enzyme-like 1) significantly increased the pathogenicity of HD1 to the resistant population. In an attempt to link the multitude of factors reportedly influencing resistance to Bt with the well-characterized loss of toxin binding, we also considered Bt-resistance models in P. xylostella and other insects.

Cloning and Tissue-Specific Expression of a Chitin Deacetylase Gene from Helicoverpa armigera (Lepidoptera: Noctuidae) and Its Response to Bacillus thuringiensis

Chitin deacetylases (CDAs) convert chitin into chitosan, the N-deacetylated form of chitin, which influences the mechanical and permeability properties of structures such as the cuticle and peritrophic matrices. In this article, a new CDA encoding gene, Hacda2, was cloned by reverse transcription-polymerase chain reaction method in Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae), with an open reading frame of 1,611 bp. The deduced protein composed of 536 amino acid residues with a signal peptide, a chitin-binding domain, a low-density lipoprotein receptor class A domain, and a polysaccharide deacetylase-like catalytic domain. The highest expression level of Hacda2 was detected in fat body among tissues tested in the fifth-instar larvae using real-time quantitative polymerase chain reaction method. Feeding of Bacillus thuringiensis (Bt) (Bacillales: Bacillaceae) diet changed the expression level of Hacda1, Hacda2, Hacda5a, and Hacda5b significantly and differentially in the third-instar larvae. Hacda5a and Hacda5b expression were initially down-regulated and then up-regulated, whereas, the expression level of Hacda1 and Hacda2 was suppressed constantly postfeeding on Bt diet. These results suggested that HaCDAs may be involved in the response of H. armigera larvae to Bt and may be helpful to elucidate the roles of HaCDAs in the action of Bt cry toxin. The potential of HaCDAs to be used as synergists of Bt insecticidal protein needs to be further tested.

Reproductive workers show queenlike gene expression in an intermediately eusocial insect, the buff-tailed bumble bee Bombus terrestris

Bumble bees represent a taxon with an intermediate level of eusociality within Hymenoptera. The clear division of reproduction between a single founding queen and the largely sterile workers is characteristic for highly eusocial species, whereas the morphological similarity between the bumble bee queen and the workers is typical for more primitively eusocial hymenopterans. Also, unlike other highly eusocial hymenopterans, division of labour among worker subcastes is plastic and not predetermined by morphology or age. We conducted a differential expression analysis based on RNA-seq data from 11 combinations of developmental stage and caste to investigate how a single genome can produce the distinct castes of queens, workers and males in the buff-tailed bumble bee Bombus terrestris. Based on expression patterns, we found males to be the most distinct of all adult castes (2411 transcripts differentially expressed compared to nonreproductive workers). However, only relatively few transcripts were differentially expressed between males and workers during development (larvae: 71 and pupae: 162). This indicates the need for more distinct expression patterns to control behaviour and physiology in adults compared to those required to create different morphologies. Among female castes, reproductive workers and their nonreproductive sisters displayed differential expression in over ten times more transcripts compared to the differential expression found between reproductive workers and their mother queen. This suggests a strong shift towards a more queenlike behaviour and physiology when a worker becomes fertile. This contrasts with eusocial species where reproductive workers are more similar to nonreproductive workers than the queen.

Dissimilar Regulation of Antimicrobial Proteins in the Midgut of Spodoptera exigua Larvae Challenged with Bacillus thuringiensis Toxins or Baculovirus

Antimicrobial peptides (AMPs) and lysozymes are the main effectors of the insect immune system, and they are involved in both local and systemic responses. Among local responses, midgut immune reaction plays an important role in fighting pathogens that reach the insect body through the oral route, as do many microorganisms used in pest control. Under this point of view, understanding how insects defend themselves locally during the first phases of infections caused by food-borne pathogens is important to further improve microbial control strategies. In the present study, we analyzed the transcriptional response of AMPs and lysozymes in the midgut of Spodoptera exigua (Lepidoptera: Noctuidae), a polyphagous pest that is commonly controlled by products based on Bacillus thuringiensis (Bt) or baculovirus. First, we comprehensively characterized the transcripts encoding AMPs and lysozymes expressed in S. exigua larval midgut, identifying 35 transcripts that represent the S. exigua arsenal against microbial infection. Secondly, we analyzed their expression in the midgut after ingestion of sub-lethal doses of two different pore-forming B. thuringiensis toxins, Cry1Ca and Vip3Aa, and the S. exigua nucleopolyhedrovirus (SeMNPV). We observed that both Bt toxins triggered a similar, wide and in some cases high transcriptional activation of genes encoding AMPs and lysozymes, which was not reflected in the activation of the classical systemic immune-marker phenoloxidase in hemolymph. Baculovirus ingestion resulted in the opposed reaction: Almost all transcripts coding for AMPs and lysozymes were down-regulated or not induced 96 hours post infection. Our results shed light on midgut response to different virulence factors or pathogens used nowadays as microbial control agents and point out the importance of the midgut immune response contribution to the larval immunity.