- Open Access
Bacterial microbiota of Aedes aegypti mosquito larvae is altered by intoxication with Bacillus thuringiensis israelensis
Parasites & Vectorsvolume 11, Article number: 121 (2018)
Insect microbiota is a dynamic microbial community that can actively participate in defense against pathogens. Bacillus thuringiensis (Bt) is a natural entomopathogen widely used as a bioinsecticide for pest control. Although Bt’s mode of action has been extensively studied, whether the presence of microbiota is mandatory for Bt to effectively kill the insect is still under debate. An association between a higher tolerance and a modified microbiota was already evidenced but a critical point remained to be solved: is the modified microbiota a cause or a consequence of a higher tolerance to Bt?
In this study we focused on the mosquito species Aedes aegypti, as no work has been performed on Diptera on this topic to date, and on B. thuringiensis israelensis (Bti), which is used worldwide for mosquito control. To avoid using antibiotics to cure bacterial microbiota, mosquito larvae were exposed to an hourly increasing dose of Bti during 25 hours to separate the most susceptible larvae dying quickly from more tolerant individuals, with longer survival.
Denaturing gradient gel electrophoresis (DGGE) fingerprinting revealed that mosquito larval bacterial microbiota was strongly affected by Bti infection after only a few hours of exposure. Bacterial microbiota from the most tolerant larvae showed the lowest diversity but the highest inter-individual differences. The proportion of Bti in the host tissue was reduced in the most tolerant larvae as compared to the most susceptible ones, suggesting an active control of Bti infection by the host.
Here we show that a modified microbiota is associated with a higher tolerance of mosquitoes to Bti, but that it is rather a consequence of Bti infection than the cause of the higher tolerance. This study paves the way to future investigations aiming at unraveling the role of host immunity, inter-species bacterial competition and kinetics of host colonization by Bti that could be at the basis of the phenotype observed in this study.
Bacillus thuringiensis (Bt) are entomopathogenic bacteria widely used as biological insecticides for the control of agricultural pests and disease-carrying insect populations . More than a hundred Bt subspecies have been described to date. Each subspecies produces a specific set of one or more toxins as a crystal during sporulation [2, 3]. After their ingestion by an insect together with bacterial spores, toxins perforate and disrupt the insect gut to allow Bt colonization of the hemocel, ultimately leading to host septicemia . Each Bt subspecies exhibits a high level of specificity toward an insect phylogenetic group [5, 6]. This adaptation to their host with the selection and conservation of key virulence factors makes Bt a pathogen rather than an opportunistic bacterium [7, 8]. Nevertheless, this view has been challenged during the last decade with contrasting reports about the role of the insect gut microbiota on the toxicity and infection capacity of Bt ; some authors suggested that host-associated microbiota might be mandatory for Bt to kill insects .
The microbiota of insects is a dynamic microbial community that shapes many host life history traits [10,11,12]. Interestingly, how the microbiota interacts with Bt toxicity and infection process is still far from being elucidated. Previous studies reported possible competition between Bt and microbiota. Bt has been shown to inhibit the growth of gut bacteria by producing bacteriocin  while microbiota can also inhibit Bt growth as well as degrade its toxins [14,15,16]. This protective role of gut microbial community against pathogens colonization, notably by niche occupation, nutrient competition or immune priming, is called “colonization resistance” [11, 17, 18]. Conversely, Bt and microbiota may also exhibit beneficial interactions. For instance, Bt protoxins can be activated into toxins by proteases produced by bacteria from insect gut . Moreover, some opportunistic bacteria, taking advantage of gut damages induced by Bt toxins, can spread across insect tissues and participate to host septicemia [20, 21].
In 2006, Broderick et al.  reported that gut microbiota was mandatory for Bt toxicity to gypsy moth larvae. Since then, several studies were conducted and led to contrasting results. Different research groups provided support for a key role of the microbiota on Bt toxicity based on an increased tolerance to Bt after gut bacterial curation by antibiotic treatments [9, 22, 23]. Interestingly, the involvement of gut microbiota seems to be species-dependent . In contrast, some scientists argued that the observed effect is mainly due to the residual antibiotic effect on Bt, resulting in a biased outcome of the experiments [7, 25, 26]. Additional studies performed by other groups found no experimental support for microbiota involvement in Bt toxicity, even on similar insect species [14, 21, 27, 28]. In light of these contrasted results, the relationship between microbiota and Bt toxicity remains unclear and it is therefore a burning topic to address . One might wonder whether these differences are exclusively due to the experimental design (antibiotics used, gut curation procedure) or if they are linked to the gut community itself, known to be highly influenced by the environment and to vary among species [10, 29] and among populations within species [30, 31]. The latter could explain the differences observed from one study to another at both intra- and inter-species levels.
Several studies have already shown that a higher tolerance to Bt is associated with a difference in bacterial community composition [32, 33]. Nevertheless, whether changes in microbiota content are the cause of an increased tolerance or a consequence of Bt infection remains unclear. Moreover, whether microbiota-Bt interaction effect on Bt toxicity is passive or active (i.e. exclusively opportunistic or evolutionary selected) is yet to be investigated. In the present study, we aimed at providing new insights to disentangle these two hypotheses. To do so, we used Aedes aegypti mosquito larvae, as the role of dipteran microbiota on Bt toxicity has not been studied to date, and B. thuringiensis israelensis (Bti), which is widely used as a biological insecticide to control mosquito populations. To avoid the bias associated with microbiota curation procedures, we studied the dynamics of the microbiota upon exposure to increasing Bt dose, which allowed us to conclude that microbiota modification is a rapid process occurring as a consequence of Bti infection.
Mosquito strain and rearing conditions
Although environment is known to influence microbial community, Aedes aegypti larvae are able to shape their microbiota that eventually significantly differ from water’s microbial community . The present work was performed on an A. aegypti laboratory strain (Bora-Bora), maintained for years in the laboratory but still exhibiting a high genetic variability . Adults were maintained in standard insectary conditions (27 °C, 14/10 h light/dark photoperiod, 80% relative humidity) in insect rearing cages and fed ad libitum with honey as previously described . After reproduction, adult females were blood-fed on mice twice a week and eggs were laid on Whatman paper disposed into crystallizing glass dishes (300 ml capacity) containing tap water. Papers were collected once or twice a week, left to dry and stored at room temperature for less than two months. For the experiment, papers containing eggs from different egg laying dates were put in water containing hay pellets to decrease the amount of oxygen in the water and promote egg hatching. Larvae were reared in tap water and fed with hay pellets in standard insectary conditions.
Bacillus thuringiensis israelensis (Bti) production
To avoid any effect of the insecticide formulation, a non-formulated Bacillus thuringiensis israelensis (Bti) toxins/spores suspension was produced in the laboratory using spores isolated from commercial Bti VectoBac WG. Suspensions of spores and crystals were produced on nutrient agar medium as previously described . Full sporulation and crystal production was verified under microscope. The quality of Bti production was assessed on SDS-PAGE. The concentration of the Bti suspensions was determined as the weight of dry pellet divided by the volume of water used for resuspension. The Bti suspension was conserved in water at -20 °C until use.
Phenotypic characterization of Bti tolerance
Considering that the mosquito strain used is highly susceptible to Bti, using a single dose of Bti would have not allowed discriminating precisely the most tolerant individuals from the most susceptible ones, as they would have all died within a very short time frame. Therefore, mosquito larvae were exposed to an hourly increasing dose of Bti to obtain a phenotypic distribution that maximize the differences between the most tolerant and the most susceptible larvae using the survival time as a proxy . A total of 560 early third-instar larvae were individually isolated in a plastic cup containing 20 ml of tap water. They were disposed in cups the day before the experiment without food to ensure that gut microbial community of all individuals were similar at T = 0. Thirty larvae, referred to as “Control” group, were unexposed to Bti, and were sampled at the beginning of the experiment. The 530 other larvae were exposed to Bti. At T = 0 and every hour afterwards, a dose of 27 μl of a 40 mg/l (1 μg) suspension of Bti spores/crystals was added within each plastic cup. Larval mortality was monitored every 15 min until the sixth hour of experiment and every 30 min onwards. Each dead larva was immediately collected and stored individually in 70% ethanol before performing DNA extraction. The frequent monitoring of larval mortality and the immediate sampling and storage of larvae limited the development of bacteria within larval tissues, which ensures that it did not bias microbiota analyses. The larval instar of each individual was determined during sampling to confirm they all remained as third-instars, and therefore ensured the maintenance of their microbiota. Each larva was labeled with a “T number” corresponding to the time point at which it died, and followed by a second number that gives the replicate number whenever several larvae died at the same time point (e.g. T4-5 is the fifth larva that died after four hours of exposure to Bti). Based on the phenotypic distribution, larvae exposed to Bti were grouped in three different categories based on their tolerance level: the 20 larvae that died during the first 6 h were qualified as “Susceptible”, the 30 larvae that survived more than 11 h were qualified as “Tolerant”, and the remaining ones (480 larvae in total) that survived between 6 and 11 h as “Intermediate”.
DNA extraction from larvae
A total of 15, 18, 15 and 20 larvae from the “Control”, “Susceptible”, “Intermediate” and “Tolerant” groups, respectively, were analyzed. Prior to DNA extraction, larvae were surface-disinfected to avoid any environmental contamination. Individuals were rinsed three times in sterile water, surface-disinfected for 5 min with 70% ethanol, and rinsed five times with sterile water. DNA extraction from larvae was slightly adapted from Minard et al.  who worked on adult mosquitoes. Considering that most, if not all, bacterial microbiota is contributed by the gut , gut larvae were not dissected. DNA extraction was performed on whole larvae to avoid potential bias induced by the dissection process (e.g. partial loss of gut content during extraction). Briefly, each larva was introduced into a 2 ml Eppendorf tube containing 5 mm diameter inox beads, plunged into liquid nitrogen for 10 s, and then crushed twice for 45 s using a Bioblock scientific MM 2000 mill (Retsch, Eragny sur Oise, France). After incubation with the extraction buffer, each sample was treated with 4 μl of RNase (100 mg/ml), and then kept at 37 °C for 5 min. Lipids and proteins were extracted with phenol-chloroform-isoamyl alcohol (25:24:1; v/v/v) and chloroform-isoamyl alcohol (24:1; v/v) steps, and DNA extracted with isopropyl alcohol. Samples were then centrifuged for 50 min at 13,200× rpm at 4 °C. DNA pellets were rinsed twice with 75% cold ethanol, air-dried under laminar flow, and re-suspended in 20 μl of TE buffer (10 mM Tris, 1 mM EDTA). The quantity of DNA was measured based on the absorbance at 260 nm and its quality was assessed by the absorbance ratio A260/280 and A260/230 (SAFAS UVmc2, Monaco).
PCR amplification of 16S rRNA fragments
A nested PCR approach was used to generate the PCR-DGGE profiles from previously extracted larval genomic DNA. The first PCR reaction was performed using the primer set pA (5'-AGA GTT TGA TCC TGG CTC AG-3') and pH (5'-AAG GAG GTG ATC CAG CCG CA-3') . PCR amplifications were carried out in 25 μl. Each reaction contained 23 μl of PCR mix which included: 200 μM of dNTPs, 0.5 μM of each primer (Thermo Fisher Scientific, Illkirch, France), 0.025 mg/ml of T4 gene 32 Protein (Roche, Boulogne-Billancourt, France), and 0.126 U/μl of Expand High Fidelity Enzyme Mix in 1× Expand High Fidelity Buffer containing MgCl2 (Roche), and nuclease-free water. Reactions were completed with 2 μl of DNA template (30 ng/μl). The amplification was carried out by performing an initial denaturation step at 95 °C for 3 min, followed by 35 cycles at 94 °C for 30 s, 55 °C for 40 s, and 72 °C for 90 s, and finalized by an extension step at 72 °C for 10 min.
The second PCR reaction targeted a fragment of approximately 200 bp of the bacterial conserved V3 region of 16S rRNA genes and was performed with a broad range bacterial primer set V3F-GC (5'-GCC GCC CGC CGC GCG CGG CGG GCG GGG CGG GGG CAC GGG GGG ACT CCT ACG GGA GGC AGC AG-3') and V3R (5'-ATT ACC GCG GCT GCT GG-3') [41, 42]. The nested PCR was carried out in 50 μl volume containing 3 μl of the first PCR run as template. Each PCR reaction contained 200 μM of dNTPs, 0.5 μM of each primer, 0.025 mg/ml of T4 gene 32 Protein, 0.04 U/μl of TaqDNA polymerase (Invitrogen) in 1× reaction PCR buffer without MgCl2, 1.5 mM MgCl2, and nuclease-free water. Amplification started by a denaturation step at 94 °C for 2 min, first followed by seven cycles at 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 1 min, second by 21 cycles at 92 °C for 30 s, 55 °C for 30 s, and 72 °C for 80 s, and finalized by an extension step at 72 °C for 10 min. All PCRs were run in a T Gradient thermocycler (Biometra, Göttingen, Germany). PCR product concentrations were assessed as previously described. Negative and positive controls were added to each PCR mix.
DGGE community fingerprinting
DGGE procedure was conducted using the Ingeny PhorU system (Apollo Intruments, Compiègne, France) as previously described . PCR products obtained from the V3 region amplification (3.5 μg per lane) were loaded onto 6% (w/v) polyacrylamide gel containing a 35–65% denaturant gradient of urea and formamide, and run in 1× TAE buffer at 60 °C for 16 h at 100 V. All gels were run with the same reference marker (1 kb Plus Ladder, Thermo Fisher Scientific) for normalization purpose in computer analyses. After the run, gels were stained using SYBR green (Thermo Fisher Scientific) for 30 min in the dark, rinsed with water, and photographed under UV-light. DGGE profiles were digitized using a camera and stored as TIFF files for computer analyses.
Sequencing of DGGE fragments and sequence analyses
The most abundant DGGE bands and some bands showing differences between conditions were excised from the gels with a sterile scalpel and rinsed individually three times in sterile ultra-pure water. Bands were incubated at 65 °C in the last washing solution and allowed to diffuse overnight at 4 °C in 100 μl of sterile water. Two microliters of eluate from individual bands were used to reamplify PCR products using the same bacterial primer set, except that primer V3F did not contain the GC clamp. Reaction conditions were the same as those described above. PCR products were sequenced at Biofidal-DTAMB (FR Bio-Environment and Health, Lyon, France). The sequences were analyzed with the BLASTN program. Sequences obtained were assigned to KC867313, KY124158, KM488465, KY608158, JQ58869, KY608117, FR821125, KY608117, KY124158 accession numbers in the GenBank database.
Analysis of DGGE patterns and statistical analyses
Gels were normalized using BioNumerics v.7.1 software (Applied Maths, Sint-Martens-Latem, Belgium). After normalization, a unique matrix containing the information related to the presence/absence of each band as well as its relative intensity was generated for all samples. Community composition and structure of each treatment group were compared by clustering lanes by the Jaccard similarity coefficient implemented in the software and using the unweighted-pair group method with arithmetic mean (UPGMA), rolling-disk background subtraction, and optimization at 0.5% [44,45,46]. Dice similarity coefficient was also tested and gave similar outputs (data not shown). This band-matching surface matrix was used as support to further DGGE profile analyses.
To investigate a potential negative correlation between intensity of bands corresponding to Bti and to Acinetobacter, Principal Components Analysis (PCA) was performed using princomp function in R 2.14.1 software . The band-matching surface matrix of the two bands of Bti (at 312 and 288 bp) and the two bands of Acinetobacter sp. (at 681 and 649 bp) were used as input and correlation circles were generated. For each band, the cosinus was calculated as a proxy of the linear correlation with the first three components. Cosinus2 were also calculated to indicate the quality of their representation by each component (from 0, bad representation; to 1, excellent representation).
The band-matching surface matrix including all bands was used to calculate Shannon and Simpson’s diversity indices and Pielou’s evenness values, and to analyze community structure by exporting it into PRIMER-E v.6.1.16 software (PRIMER-E Ltd, Plymouth, UK). The resemblance matrix was obtained using Bray-Curtis dissimilarity and nonmetric multidimensional scaling (NMDS) graphs were generated from our dataset previously modified by a fourth-root transformation . The values from Bray-Curtis index varies from 0 (similar) to 1 (completely dissimilar community composition). Euclidian distances were also tested to build the resemblance matrix and gave similar results (data not shown). An associated statistical analysis was performed based upon the analysis of similarity (ANOSIM, one-way analysis, 5000 permutations). The associated global R described the percentage of permutations related to the P-value, and the stress value indicated how faithful the relationships among samples are represented in the ordination plot.
In parallel to community structure analyses, the number of bands (taken as indicative for species richness), the Shannon (H’) and Simpson (1-λ') indices of diversity (with Simpson’s index being given more weight on dominant species compared to Shannon’s), and the Pielou’s index (a measurement of the community evenness) were used, which all together offer an insightful picture of the community α-diversity [49, 50]. Indices were all checked for normality and transformed when necessary. Only Simpson index values required log-transformation. Average differences between treatment groups were analyzed with linear models (package stats, ). Whenever ‘groups’ explained significant variation, post-hoc Tukey’s tests were used to compare pairs of groups (package multcomp, ).The normality of residuals was tested using Shapiro tests. Mean values of models and other averages were reported with their standard error. Statistical tests were based on threshold α = 0.05 and considered significant when P < 0.05. Analyses were conducted in R v.2.14.1 software .
Results and discussion
Larvae from the Bora-Bora Ae. aegypti mosquito strain were individually exposed to an hourly increasing dose of Bti (1 μg suspension of Bti spores/crystals applied every hour) in order to separate the most susceptible individuals from the most tolerant ones (Fig. 1). While this mosquito strain has been reported as highly susceptible to Bti , our phenotyping experiment revealed that individual survival upon Bti exposure was highly variable, with some larvae dying after 2.5 h of exposure while others survived 25 h exposure to a final dose of Bti 15 times greater. Larvae were attributed to three different groups based on their level of susceptibility (Fig. 1): the twenty most susceptible and the thirty most tolerant were attributed to “Susceptible” and “Tolerant” groups, respectively, while the 480 others (90.6% of the total number of exposed larvae) belonged to the “Intermediate” group. A fourth “Control” group was composed of unexposed larvae sampled at T = 0.
The stability of the microbiota of third instar larvae throughout the 24 h time frame of the experiment was verified by maintaining unexposed larvae during 0, 5, 10 and 24 h in the same laboratory conditions as the phenotyping experiment. As expected, bacterial communities did not change over time in unexposed larvae (Fig. 2, Additional file 1: Figure S1). It is well known that microbiota drastically change between developmental stages, and even between instars in the same stage [11, 53,54,55]. Nevertheless, it is quite stable within the same instar as long as the environment, which is known to be one of the major factor of microbiota shaping, is stable [11, 56]. The laboratory offers a fully controlled environment and the larvae were maintained overnight prior the experiment for ensuring that their bacterial community was stabilized. Therefore, the “Control” group is representative of the mosquito laboratory larval microbial community in the absence of insecticide and is stable throughout the 24 h time frame of the experiment.
To investigate the composition and diversity of bacterial communities of larvae from each group, PCR-DGGE fingerprints of the hypervariable V3 region of the rrs gene were produced from extracted genomic DNA larvae (Fig. 3). Some representative DGGE bands were excised from the gel, re-amplified and sequenced. Among the single sequences, BLAST analyses revealed that sequences were affiliated mostly with Bacillus, Acinetobacter, Staphylococcus and Delftia (Additional file 2: Figure S2).
Two bands were identified as Bti in the DGGE profiles (Fig. 3, Additional file 2: Figure S2). These bands were absent from profiles of most unexposed larvae, which is expected for the “Control” group (Fig. 3, Table 1). These bands were found in 89% and 100% of DGGE profiles of individuals from “Susceptible” and “Intermediate” groups, respectively. Their relative intensity ranged from 8.2 to 12.4% of total bands intensity of the profiles (Table 1). Interestingly, Bti bands were found in only 65% of “Tolerant” profiles. While the relative intensity of the higher band of Bti in “Tolerant” was 2.3 and 1.5-fold lower than the same one in “Susceptible” and “Intermediate” groups, respectively, the lower band had a similar relative intensity in all three exposed groups (11.8–12.7%, Table 1). Therefore, an overall decrease of the number of individuals with detectable Bti and of the proportion of Bti in the detectable bacterial community was observed in the most tolerant individuals compared to the most susceptible ones (Table 1). Considering that they were exposed to hourly increasing doses of Bti, the most tolerant larvae were exposed to higher doses for a longer time. If the larvae that survived longer (“Tolerant” group) did not immunologically respond to impede Bti infection, one could expect an increase in the proportion of Bti in the bacterial community, and we observed the opposite here. This suggests that, in spite of fresh Bti spores and crystals suspension added every hour to the water, the most tolerant larvae were able to control the infection, leading to a partial to complete clearance of Bti from their body. Nevertheless, at this point, these results do not allow determining whether Bti clearance is driven, at least partially, by the bacterial microbiota or exclusively under control of the larval immune system.
Interestingly, bands corresponding to Acinetobacter were negatively correlated to those of Bti (Fig. 4). While Acinetobacter was frequently associated with the “Control” group at a high intensity, the bands tended to decrease in intensity or disappear in the three other modalities raising questions about the possible antagonistic interaction between this bacterium/bacterial genus and Bti. Acinetobacter are ubiquitous bacteria that can be found in water, soil and living organisms . So far, only one study has reported the role played by Acinetobacter in the biology of a fly insect species Stomoxys calcitrans, for which these bacteria were required to ensure the complete development of larvae . Nonetheless, and in accordance with a potential antagonistic role, the Acinetobacter sp. strain KNF2022 was previously shown to produce an antiviral compound with inhibitory effects on the tobacco mosaic virus . Members of Acinetobacter are hosted by a number of mosquito species . However, most studies on mosquito-associated bacterial communities focused on adult stage and solely few studies have attempted to assess the infection status of Acinetobacter in larval stages. Recently, David et al.  showed that Acinetobacter species belonged to the core bacterial microbiota associated with Ae. aegypti as the mosquito stably harbors this bacterium throughout its lifespan.
Four diversity indices were calculated based on DGGE profiles to provide an overview of the larval bacterial community within each group (Fig. 5): species richness (total number of bands), species diversity (Shannon and Simpson indices), and species evenness (Pielou’s index). All groups significantly differed from each other, except those bearing “Susceptible” and “Intermediate” larvae (Fig. 5, Additional file 3: Table S1). The most tolerant larvae had the lowest species richness and diversity indices, and the lowest evenness index, indicating that few bacterial species became dominants (Additional file 4: Table S2). In addition, the resemblance matrix upon DGGE lanes showed that bacterial communities associated with tolerant larvae had the lowest similarity among each other (Additional file 5: Table S3). This indicates that the bacterial communities from tolerant larvae strongly differed from one individual to another.
Based on the resemblance matrix generated from DGGE profiles, pairwise comparisons from the ANOSIM analysis confirmed that all groups significantly differed from each other (Table 2). The clustering patterns observed in the non-metric multi-dimensional scaling (NMDS) plots revealed that bacterial community structure from each group did not overlap with each other (Fig. 6, Additional file 6: Figure S3). Strikingly, after only few hours of Bti exposure, the bacterial communities of the most susceptible larvae already harbored a bacterial microbiota significantly different from the control ones (Table 1, Fig. 6). The NMDS representation confirmed the higher dispersion among the “Tolerant” group compared to other groups (Fig. 6, Additional file 4: Table S2, Additional file 5: Table S3). Furthermore, this higher dispersion across microbiota of the most tolerant individual larvae seems to be associated with the survival time as there is a gradient from the larvae that survived 13–14 hours to the ones that survived 24.5 hours in the NMDS representation (Fig. 6). This cannot be dissociated from a potential dose-effect, considering that Bti concentration increased over time. Last but not least, the “Control” group is not embracing all the other groups but is rather separated with low dispersion in the NMDS (Fig. 6, Additional file 7: Figure S4). This observation excludes the selection in the more tolerant individuals of a specific bacterial community from the diversified microbiota present before exposure, otherwise “Control” individuals would have been spread in the NMDS analysis and be overlapping the other conditions. The pattern observed suggests that after exposure to Bti, the bacterial community is deeply affected and modified, comparatively to unexposed larvae. Therefore, bacterial microbiota composition associated with different levels of larval tolerance is here not the cause of the level of tolerance but rather a consequence of Bti infection.
It is noteworthy that control larvae were sampled alive while most exposed larvae were sampled dead, as a consequence of the phenotyping procedure. One could argue that the microbiota of the larvae begins to change when larvae approach to death and that this could explain the differences observed between the different larvae. This effect seems to be negligible considering that the most tolerant larvae still alive after 25 hours of exposure exhibited bacterial communities highly comparable to those dying only few hours before.
Based on our results, a scenario of the impact of the Bti infection on bacterial microbiota can be depicted. After the first exposure, microbiota began to change upon Bti ingestion while the most susceptible larvae died (Fig. 2). Only the most tolerant larvae survived by getting rid of, or significantly reducing, the quantity of Bti present in their body (Fig. 3). This process is likely to be unrelated to the microbiota, considering that no key bacteria taxa were found to be over-abundant in the “Tolerant” group in comparison with the “Control” one. In a recent similar experiment, exposure of larvae from the same laboratory Bora-Bora strain of A. aegypti to sublethal doses of Bti spores and crystals induced a significant 2.5- and 4.9-fold increased gene expression of two antimicrobial peptides (AMPs), a defensin and a cecropin, respectively (Table S3 in ). This suggests that larvae might directly control Bti by producing antimicrobial peptides (AMPs). While cecropin are generally more active against GRAM negative than GRAM positive bacteria such as Bt , anti-Bt activity of defensins has already been evidenced in lepidopterans  and coleopterans . Such an increased AMP gene expression has also been characterized in other insects exposed to different Bt subspecies [65,66,67,68]. In addition, the involvement of AMPs in tolerance to Bt has been functionally validated by RNAi experiments . Consistent with our results, this could indicate that insect larvae exposed to both toxins and spores might trigger an increased immune response to clear off bacteria out from host body while healing their gut damaged by the toxins . Increased AMPs expression can be at the basis of the strong decreased diversity and evenness of the bacterial communities observed in the “Tolerant” larvae (Fig. 5). The impact of the immune response of the most tolerant larvae on their microbiota might be time-dependent and/or dose-dependent considering that larvae from “Tolerant” group were collected during a wide range of time (from 11 to 25 h), therefore exposed to different amounts of Bti spores/crystals. Altogether, our results suggest that a modified microbiota is associated with a higher tolerance of mosquitoes to Bti but that it does not seem to be responsible for it. Higher tolerance is rather due to insect intrinsic, probably genetic, mechanisms yet to be further characterized .
Experiments linking kinetics of expression of candidate genes for immunity and overall microbial communities (e.g. by including fungi together with bacteria) are needed to better understand the association between microbiota diversity and mosquito tolerance. Although apparently not playing a direct role, bacteria may act afterwards during spore germination in larvae. Moreover, considering that immune priming induced by Bt exposure is transmissible to the next generation [71, 72] and that AMPs are driving gut colonization by symbiotic bacteria during the embryonic development , performing the same experiments on successive generations might allow identifying long term adaptation to Bt and uncover a potentially unexpected role of gut microbiota in such trans-generational Bt adaptation.
The role of microbiota in the toxicity of Bt to insects is still a burning question to be solved. It has been investigated in several species, mostly lepidopterans, by using antibiotics to cure insect gut from bacteria. Nevertheless, this approach is controversial, due to the direct effects of antibiotics on the insect physiology. Here, we provided an alternative approach to study this question, by exposing the larvae to increasing doses of Bti and comparing the microbiota of the most susceptible from that of the most tolerant individuals. Moreover, this allowed us to investigate whether a modified microbiota is the cause or a consequence of Bt exposure, which remained an open question. Our study was conducted on mosquito larvae with the israelensis subspecies of Bt, which has never been used previously on this topic. Our results revealed that a modified microbiota is associated with a higher tolerance of mosquitoes to Bti. Bacterial microbiota from the most tolerant larvae showed the lowest diversity but the highest inter-individual differences. The proportion of Bti in the host tissue was reduced in the most tolerant larvae as compared to the most susceptible ones, suggesting an active control of Bti infection by the host. Such modification of bacterial microbiota seems to be rather a consequence of Bti infection than the cause of the higher tolerance. This study paves the way to future investigations aiming at unraveling the role of host immunity, inter-species bacterial competition and kinetics of host colonization by Bti that could be at the basis of the phenotype observed in this study.
analysis of similarity
Bacillus thuringiensis israelensis
denaturing gradient gel electrophoresis
non-metric multi-dimensional scaling
principal components analysis
Bravo A, Likitvivatanavong S, Gill SS, Soberon M. Bacillus thuringiensis: A story of a successful bioinsecticide. Insect Biochem Molec. 2011;41(7):423–31.
Berry C, Crickmore N. Structural classification of insecticidal proteins - towards an in silico characterisation of novel toxins. J Invertebr Pathol. 2016;142:16–22.
Joung KB, Côté JC. A single phylogenetic analysis of Bacillus thuringiensis strains and bacilli species inferred from 16S rRNA gene restriction fragment length polymorphism is congruent with two independent phylogenetic analyses. J Appl Microbiol. 2002;93:1075–82.
Vachon V, Laprade R, Schwartz JL. Current models of the mode of action of Bacillus thuringiensis insecticidal crystal proteins: a critical review. J Invertebr Pathol. 2012;111(1):1–12.
de Maagd RA, Bravo A, Berry C, Crickmore N, Schnepf HE. Structure, diversity, and evolution of protein toxins from spore-forming entomopathogenic bacteria. Annu Rev Genet. 2003;37:409–33.
de Maagd RA, Bravo A, Crickmore N. How Bacillus thuringiensis has evolved specific toxins to colonize the insect world. Trends Genet. 2001;17(4):193–9.
Raymond B, Johnston PR, Nielsen-LeRoux C, Lereclus D, Crickmore N. Bacillus thuringiensis: an impotent pathogen? Trends Microbiol. 2010;18(5):189–94.
Lereclus D, Agaisse H, Grandvalet C, Salamitou S, Gominet M. Regulation of toxin and virulence gene transcription in Bacillus thuringiensis. Int J Med Microbiol. 2000;290(4–5):295–9.
Broderick NA, Raffa KF, Handelsman J. Midgut bacteria required for Bacillus thuringiensis insecticidal activity. Proc Natl Acad Sci USA. 2006;103(41):15196–9.
Blum JE, Fischer CN, Miles J, Handelsman J. Frequent replenishment sustains the beneficial microbiome of Drosophila melanogaster. mBio. 2013;4(6):e00860–13.
Engel P, Moran NA. The gut microbiota of insects - diversity in structure and function. FEMS Microbiol Rev. 2013;37(5):699–735.
Russell JA, Dubilier N, Rudger JA. Nature’s microbiome: introduction. Mol Ecol. 2014;23:1225–37.
Cherif A, Rezgui W, Raddadi N, Daffonchio D, Boudabous A. Characterization and partial purification of entomocin 110, a newly identified bacteriocin from Bacillus thuringiensis subsp. Entomocidus HD110. Microbiol Res. 2008;163(6):684–92.
Patil CD, Borase HP, Salunke BK, Patil SV. Alteration in Bacillus thuringiensis toxicity by curing gut flora: novel approach for mosquito resistance management. Parasitol Res. 2013;112(9):3283–8.
Takatsuka J, Kunimi Y. Intestinal bacteria affect growth of Bacillus thuringiensis in larvae of the oriental tea tortrix, Homona magnanima diakonoff (Lepidoptera: Tortricidae). J Invertebr Pathol. 2000;76(3):222–6.
Shan Y, Shu C, Crickmore N, Liu C, Xiang W, Song F, et al. Cultivable gut bacteria of scarabs (Coleoptera: Scarabaeidae) inhibit Bacillus thuringiensis multiplication. Environ Entomol. 2014;43(3):612–6.
Vollaard EJ, Clasener HAL. Colonization resistance. Antimicrob Agents Chemother. 1994;38(3):409–14.
Stecher B, Hardt WD. The role of microbiota in infectious disease. Trends Microbiol. 2008;16(3):107–14.
Regode V, Kuruba S, Mohammad AS, Sharma HC. Isolation and characterization of gut bacterial proteases involved in inducing pathogenicity of Bacillus thuringiensis toxin in cotton bollworm, Helicoverpa armigera. Front Microbiol. 2016;7:1567.
Mason KL, Stepien TA, Blum JE, Holt JF, Labbe NH, Rush JS, et al. From commensal to pathogen: translocation of Enterococcus faecalis from the midgut to the hemocoel of Manduca sexta. mBio. 2011;2(3):e00065–11.
van Frankenhuyzen K, Liu YH, Tonon A. Interactions between Bacillus thuringiensis subsp kurstaki HD-1 and midgut bacteria in larvae of gypsy moth and spruce budworm. J Invertebr Pathol. 2010;103(2):124–31.
Paramasiva I, Sharma HC, Krishnayya PV. Antibiotics influence the toxicity of the delta endotoxins of Bacillus thuringiensis towards the cotton bollworm, Helicoverpa armigera. BMC Microbiol. 2014;14:200.
Visweshwar R, Sharma HC, Akbar SM, Sreeramulu K. Elimination of gut microbes with antibiotics confers resistance to Bacillus thuringiensis toxin proteins in Helicoverpa armigera (Hubner). Appl Biochem Biotech. 2015;177(8):1621–37.
Broderick NA, Robinson CJ, McMahon MD, Holt J, Handelsman J, Raffa KF. Contributions of gut bacteria to Bacillus thuringiensis-induced mortality vary across a range of Lepidoptera. BMC Biol. 2009;7:11. https://doi.org/10.1186/1741-7007-7-11.
Beegle CC, Lewis LC, Lynch RE, Martinez AJ. Interaction of larval age and antibiotic on the susceptibility of three insect species to Bacillus thuringiensis. J Invertebr Pathol. 1981;37:143–53.
Ignoffo CM, Garcia C, Couch TL. Effect of antibiotics on the insecticidal activity of Bacillus thuringiensis. J Invertebr Pathol. 1977;30:277–8.
Johnston PR, Crickmore N. Gut bacteria are not required for the insecticidal activity of Bacillus thuringiensis toward the tobacco hornworm, Manduca sexta. Appl Environ Microbiol. 2009;75(15):5094–9.
Raymond B, Johnston PR, Wright DJ, Ellis RJ, Crickmore N, Bonsall MB. A mid-gut microbiota is not required for the pathogenicity of Bacillus thuringiensis to diamondback moth larvae. Environ Microbiol. 2009;11(10):2556–63.
Kikuchi Y, Hosokawa T, Fukatsu T. Insect-microbe mutualism without vertical transmission: a stinkbug acquires a beneficial gut symbiont from the environment every generation. Appl Environ Microbiol. 2007;73(13):4308–816.
Minard G, Mavingui P, Moro CV. Diversity and function of bacterial microbiota in the mosquito holobiont. Parasit Vectors. 2013;6:146.
Muturi EJ, Ramirez JL, Rooney AP, Kim CH. Comparative analysis of gut microbiota of mosquito communities in central Illinois. PLoS Negl Trop Dis. 2017;11(2):e0005377.
Caccia S, Di Lelio I, La Storia A, Marinelli A, Varricchio P, Franzetti E, et al. Midgut microbiota and host immunocompetence underlie Bacillus thuringiensis killing mechanism. Proc Natl Acad Sci USA. 2016;113(34):9486–91.
Paramasiva I, Shouche Y, Kulkarni GJ, Krishnayya PV, Akbar SM, Sharma HC. Diversity in gut microflora of Helicoverpa armigera populations from different regions in relation to biological activity of Bacillus thuringiensis delta-endotoxin Cry1Ac. Arch Insect Biochem. 2014;87(4):201–13.
Dada N, Jumas-Bilak E, Manguin S, Seidu R, Stenström TA, Overgaard HJ. Comparative assessment of the bacterial communities associated with Aedes aegypti larvae and water from domestic water storage containers. Parasit Vectors. 2014;7:391.
Paris M, Despres L. Identifying insecticide resistance genes in mosquito by combining AFLP genome scans and 454 pyrosequencing. Mol Ecol. 2012;21(7):1672–86.
Paris M, Tetreau G, Laurent F, Lelu M, Despres L, David J-P. Persistence of Bacillus thuringiensis israelensis (Bti) in the environment induces resistance to multiple Bti toxins in mosquitoes. Pest Manag Sci. 2011;67:122–8.
Tetreau G, Stalinski R, Kersusan D, Veyrenc S, David JP, Reynaud S, et al. Decreased toxicity of Bacillus thuringiensis subsp. israelensis to mosquito larvae after contact with leaf litter. Appl Environ Microbiol. 2012;78(15):5189–95.
Bonin A, Paris M, Frerot H, Bianco E, Tetreau G, Despres L. The genetic architecture of a complex trait: Resistance to multiple toxins produced by Bacillus thuringiensis israelensis in the dengue and yellow fever vector, the mosquito Aedes aegypti. Infect Genet Evol. 2015;35:204–13.
Minard G, Tran FH, Dubost A, Tran-Van V, Mavingui P, Moro CV. Pyrosequencing 16S rRNA genes of bacteria associated with wild tiger mosquito Aedes albopictus: a pilot study. Front Cell Infect Mi. 2014;4:59.
Edwards U, Rogall T, Blöcker H, Emde M, Böttger EC. Isolation and direct complete nucleotide determination of entire genes. Characterization of a gene coding for 16S ribosomal RNA. Nucl Acids Res. 1989;17(19):7843–53.
Muyzer G, de Waal EC, Uitterlinden AG. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl Environ Microbiol. 1993;59(3):695–700.
Zouache K, Raharimalala FN, Raquin V, Tran-Van V, Raveloson LH, Ravelonandro P, et al. Bacterial diversity of field-caught mosquitoes, Aedes albopictus and Aedes aegypti, from different geographic regions of Madagascar. FEMS Microbiol Ecol. 2011;75(3):377–89.
Zouache K, Voronin D, Tran-Van V, Mavingui P. Composition of bacterial communities associated with natural and laboratory populations of Asobara tabida infected with Wolbachia. Appl Environ Microbiol. 2009;75(11):3755–64.
Fromin N, Hamelin J, Tarnawski S, Roesti D, Jourdain-Miserez K, Forestier N, et al. Statistical analysis of denaturing gel electrophoresis (DGE) fingerprinting patterns. Environ Microbiol. 2002;4(11):634–43.
Kropf S, Läuter J, Eszlinger M, Krohn K, Paschke R. Nonparametric multiple test procedures with data-driven order of hypotheses and with weighted hypotheses. J Stat Plan Infer. 2004;125(1–2):31–47.
Rademaker JLW, Louws FJ, Rossbach U, Vinuesa P, de Bruijn FJ. Computer-assisted pattern analysis of molecular fingerprints and database construction. In: Akkermans ADL, van Elsas JD, de Bruijn FJ, editors. Molecular Microbial Ecology Manual. Dordrecht: Kluwer Academic Plublishers; 1999. p. 1–33.
R Development Core. Team: R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2011.
Grizard S, Dini-Andreote F, Tieleman BI, Salles JF. Dynamics of bacterial and fungal communities associated with eggshells during incubation. Ecol Evol. 2014;4(7):1140–57.
Jost L. The relation between evenness and diversity. Diversity. 2010;2(2):207–32.
Pielou EC. The measurement of diversity in different types of biological collections. J Theor Biol. 1966;13:131–44.
Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biometrical J. 2008;50(3):346–63.
Stalinski R, Tetreau G, Gaude T, Despres L. Pre-selecting resistance against individual Bti Cry toxins facilitates the development of resistance to the Bti toxins cocktail. J Invertebr Pathol. 2014;119:50–3.
Wang Y, Gilbreath TM III, Kukutla P, Yan G, Xu J. Dynamic gut microbiome across life history of the malaria mosquito Anopheles gambiae in Kenya. Plos One. 2011;6(9):e24767.
Johnston PR, Rolff J. Host and symbiont jointly control gut microbiota during complete metamorphosis. Plos Pathog. 2015;11(11):e1005246.
Moll RM, Romoser WS, Modrzakowski MC, Moncayo AC, Lerdthusnee K. Meconial peritrophic membranes and the fate of midgut bacteria during mosquito (Diptera: Culicidae) metamorphosis. J Med Entomol. 2001;38(1):29–32.
Yadav KK, Datta S, Naglot A, Bora A, Hmuaka V, Bhagyawant S, et al. Diversity of cultivable midgut microbiota at different stages of the asian tiger mosquito, Aedes albopictus from Tezpur, India. Plos One. 2016;11(12):e0167409.
Doughari HJ, Ndakidemi PA, Human IS, Benade S. The ecology, biology and pathogenesis of Acinetobacter spp.: an overview. Microbes Environ. 2011;26(2):101–12.
Lysyk TJ, Kalischuk-Tymensen L, Selinger LB, Lancaster RC, Wever L, Cheng KJ. Rearing stable fly larvae (Diptera: Muscidae) on an egg yolk medium. J Med Entomol. 1999;36(3):382–8.
Lee J-S. Acinetobacter antiviralis sp. nov., from tobacco plant roots. J Microbiol Biotechn. 2009;19(3):250–6.
David MR, Santos LM, Vicente AC, Maciel-de-Freitas R. Effects of environment, dietary regime and ageing on the dengue vector microbiota: evidence of a core microbiota throughout Aedes aegypti lifespan. Mem Inst Oswaldo Cruz. 2016;111(9):577–87.
Stalinski R, Laporte F, Despres L, Tetreau G. Alkaline phosphatases are involved in the response of Aedes aegypti larvae to intoxication with Bacillus thuringiensis subsp. israelensis Cry toxins. Environ Microbiol. 2016;18(3):1022–36.
Moore AJ, Beazly WD, Bibby MC, Devine DA. Antimicrobial activity of cecropins. J Antimicrob Chemoth. 1996;37:1077–89.
Xu XX, Zhang YQ, Freed S, Yu J, Gao YF, Wang S, et al. An anionic defensin from Plutella xylostella with potential activity against Bacillus thuringiensis. Bull Entomol Res. 2016;106(6):790–800.
Tonk M, Knorr E, Cabezas-Cruz A, Valdes JJ, Kollewe C, Vilcinskas A. Tribolium castaneum defensins are primarily active against Gram-positive bacteria. J Invertebr Pathol. 2015;132:208–15.
Crava CM, Jakubowska AK, Escriche B, Herrero S, Bel Y. Dissimilar regulation of antimicrobial proteins in the midgut of Spodoptera exigua larvae challenged with Bacillus thuringiensis toxins or baculovirus. Plos One. 2015;10(5):e0125991.
Contreras E, Benito-Jardon M, Lopez-Galiano MJ, Real MD, Rausell C. Tribolium castaneum immune defense genes are differentially expressed in response to Bacillus thuringiensis toxins sharing common receptor molecules and exhibiting disparate toxicity. Dev Comp Immunol. 2015;50(2):139–45.
Contreras E, Rausell C, Real MD. Tribolium castaneum apolipophorin-III acts as an immune response protein against Bacillus thuringiensis Cry3Ba toxic activity. J Invertebr Pathol. 2013;113(3):209–13.
Wang Q, Liu Y, He HJ, Zhao XF, Wang JX. Immune responses of Helicoverpa armigera to different kinds of pathogens. BMC Immunol. 2010;11:9.
Hwang J, Kim YRNA. interference of an antimicrobial peptide, gloverin, of the beet armyworm, Spodoptera exigua, enhances susceptibility to Bacillus thuringiensis. J Invertebr Pathol. 2011;108(3):194–200.
Castagnola A, Jurat-Fuentes JL. Intestinal regeneration as an insect resistance mechanism to entomopathogenic bacteria. Curr Opin Insect Sci. 2016;15:104–10.
Pigeault R, Garnier R, Rivero A, Gandon S. Evolution of transgenerational immunity in invertebrates. Proc Roy Soc B-Biol Sci. 2016;283:20161136.
Moret Y. “Trans-generational immune priming”: specific enhancement of the antimicrobial immune response in the mealworm beetle, Tenebrio molitor. Proc Roy Soc B-Biol Sci. 2006;273(1592):1399–405.
Login FH, Balmand S, Vallier A, Vincent-Monégat C, Vigneron A, Weiss-Gayet M, et al. Antimicrobial peptides keep insect endosymbionts under control. Science. 2011;334(6054):362–5.
This paper is dedicated to the memory of Florence-Hélène Tran who tragically left us in August 2016. We thank Sarah Ganet for helping in the use of the BioNumerics software.
This work was funded by ERA-NET BiodivERsA with the national funders ANR-13-EBID-0007-01, FWF I-1437 and DFG KL 2087/6-1 as part of the 2012–2013 BiodivERsA call for research proposals.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
Blood-feeding of adult mosquitoes was performed on mice. Mice were maintained in the animal house of the federative structure Environmental and Systems Biology (BEeSy) of Grenoble-Alpes University agreed by the French Ministry of animal welfare (agreement n° B 38 421 10 001) and used in accordance to European Union laws (directive 2010/63/UE). The use of animals for this study was approved by the ethic committee ComEth Grenoble-C2EA-12 mandated by the French Ministry of higher Education and Research (MENESR).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Picture of the DGGE gel containing unexposed larvae sampled at 0, 5, 10 and 24 h. (PDF 108 kb)
Figure S2. Composite picture of a DGGE gel containing all four groups of larvae. Bands excised and sequenced are indicated by a red rectangle. The corresponding species identified are indicated in the table. (PDF 242 kb)
Table S1. Tukey’s post-hoc test outputs from larval bacterial microbiota. Linear models were used to test the relationship between alpha-diversity metrics (diversity) and treatment groups (treatment). Tukey’s post-hoc tests were used to compare pair of treatment (Z-tests). P-values are then reported for each model and indicated in bold when significant (P < 0.05). (PDF 195 kb)
Table S2. Band counting numbers, Shannon index (diversity), Simpson index (diversity) and Pielou’s index (evenness). They are based on DGGE band-matching surfaces. Results are for larvae never exposed to Bacillus thuringiensis israelensis (Bti) (“Control”) and for larvae exposed to Bti, including larvae dead in less than 6 h of exposure to Bti (“Susceptible”), in between 6 and 11 h (“Intermediate”), and past 11 h (“Tolerant”). Results are given with mean (± standard deviation) (minimum-maximum values). (PDF 277 kb)
Table S3. Bacterial community similarity within and between treatment groups. Similarity values are based on DGGE band-matching surface matrix. Results are presented as the mean (± standard deviation) (minimum-maximum values). (PDF 9 kb)
Figure S3. Pictures of the four DGGE gels used in the analysis. (PDF 142 kb)
Figure S4. Non-metric multi-dimensional scaling plots of the bacterial communities of larvae based on DGGE lane analysis. 3D stress of the NMDS was 0.14. The NMDS representation is based on a Bray-Curtis dissimilarity matrix. Correspondence of sample names is indicated in the legend of Fig. 2. (PDF 65 kb)
About this article
- Bacillus thuringiensis
- Denaturing gradient gel electrophoresis (DGGE)
- Bacterial community fingerprinting
- Larval microbiota