Larval midgut modifications associated with Bti resistance in the yellow fever mosquito using proteomic and transcriptomic approaches
BMC Genomics volume 13, Article number: 248 (2012)
Bacillus thuringiensis var. israelensis (Bti) is a natural larval mosquito pathogen producing pore-forming toxins targeting the midgut of Diptera larvae. It is used worldwide for mosquito control. Resistance mechanisms of an Aedes aegypti laboratory strain selected for 30 generations with field-collected leaf litter containing Bti toxins were investigated in larval midguts at two levels: 1. gene transcription using DNA microarray and RT-qPCR and 2. differential expression of brush border membrane proteins using DIGE (Differential In Gel Electrophoresis).
Several Bti Cry toxin receptors including alkaline phosphatases and N-aminopeptidases and toxin-binding V-ATPases exhibited altered expression levels in the resistant strain. The under-expression of putative Bti- receptors is consistent with Bt-resistance mechanisms previously described in Lepidoptera. Four soluble metalloproteinases were found under-transcribed together with a drastic decrease of metalloproteinases activity in the resistant strain, suggesting a role in resistance by decreasing the amount of activated Cry toxins in the larval midgut.
By combining transcriptomic and proteomic approaches, we detected expression changes at nearly each step of the ingestion-to-infection process, providing a short list of genes and proteins potentially involved in Bti-resistance whose implication needs to be validated. Collectively, these results open the way to further functional analyses to better characterize Bti- resistance mechanisms in mosquitoes.
Mosquito control represents a major public health concern as mosquitoes transmit many pathogens causing fatal human diseases including malaria, filariasis, dengue, yellow fever, and Chikungunya . Vector borne diseases represent a major health threat and economic burden in disease-endemic countries and are currently expanding worldwide [2, 3]. As no specific treatment exists for most of these diseases, the most effective way of reducing the incidence of these diseases is to control the vector mosquitoes [4, 5]. Chemical insecticides still used in endemic countries have shown their limits as resistance has evolved in all target species together with environmental concerns due to their high persistence and toxicity for non-target organisms, including humans .
The bacterium Bacillus thuringiensis var. israelensis (Bti) produces a mosquitocidal toxic crystal during sporulation and represents the best alternative to chemical insecticides for mosquito larval control due to its high potency and specificity . The action of Bti begins when larvae ingest Bti spores and toxic crystals. In susceptible larvae, the toxic crystal is dissolved in the alkaline pH of the midgut, protoxins are then activated by digestive proteases to activated-toxins that bind to specific membrane receptors, form pores, disrupt the midgut epithelium, allowing spore penetration and bacterial proliferation in the host tissues [7, 8]. The receptors for mosquitocidal Bti Cry toxins are similar to the lepidopteran-active Cry toxins which utilize N-aminopeptidase, alkaline phosphatase and cadherin proteins as midgut receptors .
In contrast to Bacillus thuringiensis subspecies active against lepidopteran and coleopteran species where cases of insect resistance in the field have been reported [10–13], only one study reported Bti resistance in field mosquitoes . However, subsequent confirmations of this case have not been reported. The delay in the evolution of resistance to Bti is believed to be due to its composite toxic crystal containing four major toxins (Cry4Aa, Cry4Ba, Cry11Aa and Cyt1Aa) . Cyt toxins are known to largely enhance Cry toxins activity due to synergic effects and to drastically decrease resistance development [8, 15]. Although Bti is known to have a low persistence in the environment, recent studies suggest that it can persist and possibly proliferate in specific conditions [16–18]. In the French Rhône-Alpes region, decaying leaf litters collected in mosquito breeding sites several months after a Bti treatment revealed a high toxicity against mosquito larvae due to the presence of large amounts of Bti. This toxic leaf litter was used to select an Aedes aegypti strain in laboratory conditions. After 18 generations, the selected strain (named LiTOX) was only moderately resistant to the whole Bti toxins mixture, but up to 30 fold resistant to individual Cry toxins . Although resistance to Bti has already been selected in laboratory conditions [20, 21], this is the only reported case of resistance obtained by using field-collected material containing residual Bti toxins. Therefore, this Bti-resistant LiTOX strain provides a unique opportunity to better understand the mode of action of Bti toxins and to elucidate the mechanisms of resistance developed by mosquitoes exposed to field residual Bti toxins.
To identify the resistance mechanisms developed by the LiTOX strain, a genome scan and a transcriptome scan were previously performed on whole larvae twelve generations ago [22, 23]. The main bias of these whole-larvae approaches is that many genes are identified that may not be directly related to Bti resistance. Indeed, selection was shown to have induced many changes in the LiTOX strain, including decreased egg survival to desiccation, longer larval development time and decreased female fecundity , reflecting the evolution of resistance costs that are not directly involved in resistance to Bti toxins. Because insect midgut is the primary target site for Bti toxins our aim in the present work is to focus on constitutive expression changes in midgut proteins of resistant versus susceptible larvae. For that purpose, we combine a comparative analysis of brush border membrane proteins using 2D-DIGE (2-Dimensional Differential in Gel Electrophoresis) with a midgut transcriptome profiling using DNA microarrays. In addition, altered gene expression of known Bti Cry toxins receptors (i.e. alkaline phosphatases, cadherins, N-aminopeptidases) between the two strains were investigated using RT-qPCR. Finally, because the DiGE didn’t allow detecting proteins with high molecular size such as cadherins, we performed Western blots with anti-cadherins antibodies.
Resistance levels to Bti toxins in the LiTOX strain
After 30 generations of selection with leaf litter containing Bti, bioassays indicated that the LiTOX strain exhibited a moderate 3.5-fold resistance to commercial Bti mixture Vectobac® WG compared to the susceptible strain at the larval stage (Table 1). When Bti Cry toxins were tested separately, the LiTOX strain showed an increased resistance of 68-fold, 9-fold and 9-fold to Cry4Aa, Cry4Ba and Cry11Aa protoxins respectively. The relatively important variability observed for the LC50 for Cry4Aa toxin of the LiTOX strain is mainly due to a higher variability in larval mortality in the replicates than for the susceptible strain and for the other toxins. As resistance is not fixed yet in the LiTOX strain , this variability between replicates might reflect a large range of different combinations of Cry4A resistance alleles between individuals.
Midgut transcriptome profiling
Comparative transcriptome profiling between total mRNAs extracted from midguts of larvae from the LiTOX and the susceptible strains was performed using a DNA microarray representing 14204 of the more than 17000 Ae. aegypti transcripts identified in Vectorbase. A total of 3512 transcripts were detected in at least 5 hybridizations out of 6 [ArrayExpress: E-MTAB-1094] (Additional file 1). Among them, 24 and 46 genes were significantly over- and under-transcribed respectively in the LiTOX strain (≥3-fold and corrected P-value <0.01) (Additional file 2). Distribution of transcription ratios was well balanced between over and under-transcription ranging from 20.9-fold under-expression to 18.9-fold over-expression (Figure 1). RT-qPCR validation of transcription ratios for 15 selected genes revealed a good correspondence between the two techniques, supporting the reliability of microarray data (Additional file 3).
Differentially transcribed genes were further analyzed according to their biological function by classifying them into 13 different categories (Figure 2). Genes of unknown function represented 34% of detected transcripts while genes not assigned to any category (other functions) represented 17%. Enzymes represented 30% of detected transcripts and were strongly over-represented among under- and over-transcribed genes (55% and 60% respectively) (Figure 2B & C). Proteases were equally represented in over- and under-expressed genes, while detoxification enzymes were more often under- than over-transcribed (7 under-transcribed versus 2 over-transcribed genes). Transaminases, represented by only 11 genes in the Ae. aegypti genome, were over-represented in under- and over-transcribed genes while dehydrogenases were strongly over-represented only in over-transcribed genes (23% of enzymes compared to 10% overall).
Midgut differential proteomics
Midgut membrane proteins were compared between larvae of the LiTOX and susceptible strains using 2D-DIGE (Figure 3). Dye-swapping for each biological sample showed no dye-dependent spot changes on the gels (Additional file 4). Spot locations were reproducible between the biological replicates, but the signal intensity was higher for the second replicate, revealing additional spots differing between the two strains (Additional file 4). A total of 56 distinct protein spots differently expressed between the two strains were processed and 35 unique proteins were identified (Figure 3B, Additional file 5). The MS/MS analyses gave the same protein identifications between biological replicates for spots 2, 8, 20, 21 and 24 with high Mascot scores (from 110 to 249) while spots 14, 42 and 49, showing Mascot scores lower than 100, were assigned to different proteins (Additional file 5). Indeed, none of the spots with low Mascot scores were considered for further analyses. Different spots yielded the same identified protein for 10 proteins, with a maximum of six spots for AAEL005798 (V-ATP synthase subunit beta).
Genes and proteins differentially expressed in the LiTOX strain
Proteome analysis identified two N-aminopeptidase proteins (APN, annotated as ‘protease m1 zinc metalloprotease’) differentially expressed in the LiTOX strain (Table 2): two spots matching APN AAEL012774 were up-regulated and two of the three spots matching APN AAEL012776 were down-regulated in the LiTOX strain. Transcriptomic approach detected thirteen APN (including AAEL012774 and AAEL012776) with transcription level ranging from −1.82 to +1.96 fold changes (Additional file 1) but none was significant.
Two alkaline phosphatases (ALP) proteins, matching AAEL003313 and AAEL003298, were under-expressed in the LiTOX strain while transcriptomics identified two other ALP genes (AAEL011175 and AAEL015070) significantly over-transcribed in the LiTOX strain with transcription ratios of +4.63 and +3.95 fold respectively. All ALPs but AAEL003298 have predicted glycosylphosphatidylinositol (GPI)-anchor domains allowing them to tether to the epithelial membrane and be potential membrane-bound Cry toxin receptors (Additional file 6).
Four proteins matching ATP synthase subunits alpha, beta and epsilon, with two to seven different spots for the same protein, had levels from −1.51 to −2.39 lower in the LiTOX strain. This tendency was consistent with microarray data for AAEL008787-RA (−1.19 fold), AAEL005798 (−1.37 fold) and AAEL012035 (−1.19 fold) although P values were not significant.
A unique calreticulin protein was picked and identified from DIGE gels. Initially detected as only one big spot with +10.74 fold change, the second biological replicate allowed clearly identifying four different spots respectively −2.56, +2.75, +3.01 and +3.14 fold differentially expressed in the LiTOX strain. In microarray experiment, no significant differential transcription of this gene was found in the LiTOX strain.
Using BLASTP software, we managed to identify putative functions for the 15 transcripts of unknown function differentially transcribed in the LiTOX strain with protein identities ranging from 25 to 99% (Additional file 2). Among them, two were strongly over-transcribed in the resistant strain (AAEL013584 19-fold and AAEL010435 9.6-fold) and matched to a putative G12 protein in Ae. aegypti (77% protein identity, Additional file 2).
Among the five cytochrome P450 monooxygenase transcripts identified by microarray analyses, CYP4D24 was over-transcribed while the others (CYP6N9, CYP6Z7, CYP6Z8 and CYP9M9) were under-transcribed in the LiTOX strain. DIGE experiments identified one protein matching to a cytochrome P450 3-fold over-expressed in the LiTOX strain.
Transcriptomic data detected four metalloproteinases significantly under-transcribed from −3.16 to −5.29 fold (Additional file 2). The presence of conserved domains of soluble astacin-like metalloproteinases together with the absence of detected GPI-anchor domain (Additional file 6) suggests that these four metalloproteinases are probably secreted extracellular enzymes, explaining why they were not identified in the BBMV by the DIGE analysis.
Global and specific proteolytic activities
To determine if the modifications in protease transcription levels observed in the resistant strain result in changes in gut proteolytic activities, we compared the protease activities of secreted proteins from larval midgut of each strain using azocasein as substrate. Total proteolytic activity was 8.5% higher in the resistant strain compared to the susceptible strain (Table 3). The use of specific protease inhibitors revealed that more than 90% of the total proteolytic activity is due to serine protease for both strains. Among them, chymotrypsins and trypsins represented respectively more than 40% and 20% of the total activity in both strains. The use of the metalloproteinase inhibitor EDTA showed that 12% of the total proteolytic activity was due to metalloproteinases in the susceptible strain whereas no inhibition was measured in the resistant strain (Wilcoxon test; P-value <0.05), suggesting a strong reduction of metalloproteinase activity in the LiTOX strain.
Resistance levels to Bti toxins in the LiTOX strain
After 30 generations of selection, resistance to Cry4Aa in the LiTOX strain has more than doubled as compared to twelve generations ago, while resistance ratios did not change for Cry4Ba and Cry11Aa [19, 22]. Resistance to Bti is moderate (3.5 fold) but higher than at generation 18 (2-fold), indicating that resistance alleles are not all fixed yet. These results are consistent with previous attempts to select Ae. aegypti Culex pipiens and Cx. quinquefasciatus with Bti which obtained moderate resistance (2 to 3 fold) after 20 to 30 generations [20, 21, 25, 26]. The increased Bti resistance observed may be due to the increase in Cry4Aa resistance, and most changes observed in the present study may be related to Cry4Aa resistance. The discrepancy between Bti and Cry4Aa increased resistances is likely to be due to the presence of Cyt toxin in Bti, known to overcome Cry resistance in insects .
Midgut transcriptome and proteome analyses
Our comparison of midgut transcripts and brush border proteins between the susceptible and LiTOX strains revealed an overlapping but distinct set of transcripts/proteins differentially expressed. Transcriptome profiling with a microarray representing more than 81% of known Ae. aegypti transcripts lead to the detection of 3512 transcripts of which 70 were differentially transcribed in the LiTOX strain. This relatively low number of transcripts detected (about 25%) is probably due to the low transcription level or absence of transcription of several genes in this particular organ —the larval midgut. Little overlap was observed between the previous transcriptomic analysis, performed on whole larvae 12 generations ago using a DGETP approach , and the present study, focusing on midgut gene expression using microarrays. This is possibly due to the technical differences between the two studies and to the fact that no resistance gene is fixed yet in the LiTOX strain, indicated by the still increasing resistance to Bti and to Cry toxins . Moreover, as genes conferring resistance to Bti toxins are likely to be expressed in larval midgut, focusing on midguts rather than on whole larvae allows to considerably reduce the candidate gene dataset, and to consider only genes likely to be directly involved in resistance, rather than those only indirectly affected by selection side-effects (genetic drift) or compensatory mechanisms (resistance costs) . The 2D-DIGE analysis resolved about 400 distinct proteins in larval BBMV fractions, 56 spots were picked of which 50 differed by more than 2-fold between the two strains. The difference between the number of spots picked (56) and the number of unique identified proteins (35) is due to different spots for the same protein, as for example up to six spots observed for one V-ATPase. The multiple spots for the same protein are most likely due to post-translational modifications (glycosylation, phosphorylation) that cause shift in protein mobility. Eight spots common to the replicated DIGE experiments were picked and identified twice. Among them all but one, calreticulin (AAEL01005), showed similar levels of differential expression supporting the consistency of biological replicates (Additional file 4). Both transcriptomic and proteomic data identified more under than over-expressed genes/proteins in the LiTOX strain, which is congruent with a previous transcriptome analysis performed on whole larvae 12 generations ago . Such asymmetry is not surprising considering that mechanisms of resistance to Bt can involve a decreased activation of protoxins or a decreased toxin-binding to the epithelium membrane .
Little overlap was found between data obtained by transcriptomic and proteomic approaches. This could be explained by both biological processes and technical limitations inherent to each method. Regarding DIGE, BBMV were used, which are enriched for proteins attached to apical brush border midgut membrane via scaffolding and proteins attached to the inner membrane leaflet via acylation. Therefore, except few soluble proteins trapped in re-folded membranes, soluble intracellular proteins and proteins excreted inside the gut lumen are typically absent in BBMV preparations [28, 29]. In contrast, mRNAs extracted from whole larval midguts should be representative of all transcripts present in midgut cells. Another factor limiting overlapping data may be the consequence of the relatively stringent filtering of the microarrays dataset (3-fold threshold). Several studies also showed that mRNA transcription profiles fit poorly with protein levels because of numerous post-transcriptional regulatory activities and post-translational events [30–32]. Such events generate a high diversity of proteins while gene expression remains unchanged, and this source of variation is so far under-explored in studies on fast adaptive changes like the evolution of insecticide resistance. It is likely that the two complementary approaches used in the present work detected distinct mechanisms of resistance acting at different steps in the mode of action of Bti (i.e. crystal solubilization, toxin activation and binding to receptors).
Altered expression and activities of proteases from the LiTOX strain
Four soluble astacin-like metalloproteinases were found significantly under-transcribed in the LiTOX strain. This observation was correlated with a strong decrease of metalloproteinases activity among the enzymes secreted in the midgut lumen of larvae from the resistant strain. To our knowledge, this is the first time astacin-like metalloproteinases are associated to Bt resistance. The observed decrease in metalloproteinases in the resistant strain might reflect an alteration in Bti Cry toxins activation in the gut lumen of LiTOX larvae. Further experiments based on measuring proteolytic activities and performing bioassays with activated toxins will clarify the potential role that alteration of protoxins processing, notably for Cry4Aa, could play in the resistance phenotype.
Altered expression of known Bti- binding proteins in the LiTOX strain
To validate the expression alteration of putative Bti-receptors observed in microarrays and DIGE approaches, RT-qPCR analyses were performed on five N-aminopeptidases (APN1 to 5), two cadherins (Cad1 and Cad2) and three alkaline phosphatases (ALP1 to 3) previously described as binding proteins for Cry4Ba  or Cry11Aa [33–36] (Table 4).
The cadherin described as a Cry11Aa-receptor in Ae. aegypti (AAEL007488)  was found 1.47 fold under-transcribed in both microarrays and RT-qPCR experiments. However, no cadherin was detected by DIGE approach. The inability to detect cadherin in the DIGE analysis is not surprising as they are large proteins (>170 kDa) present in low amounts in insect brush border membranes . Blotting of BBMV using two anti-cadherin antibodies showed that most of the cadherins were degraded, even in a freshly prepared UGAL Aedes BBMV preparation, confirming that cadherins in BBMV are very unstable (Additional file 7). Western blots showed that cadherin(s), notably a ~32 kDa fragment, is strongly over-represented in the LiTOX strain compared to the susceptible strain (Additional file 7). Further analyses of the toxins-binding properties of the detected cadherin(s) are needed to better understand the role they could play in the resistance phenotype.
Alkaline phosphatases (ALPs), typically anchored by GPI-moieties, are known to be Cry toxin receptors in Lepidoptera [38, 39] and mosquitoes [36, 40]. Recently, a decrease in ALP amounts and activities were linked to Cry-resistance of larvae from three lepidopteran genera . The ALPs detected as over-transcribed in the LiTOX strain by the two transcriptomic approaches were either not identified as differentially expressed (AAEL000931, AAEL009077 and AAEL015070) or identified as under-expressed (AAEL003313) by DIGE approach. These results suggest that the lower ALP protein abundance in the epithelium membrane might rather be due to post-translational events than under-expression. Indeed, our DIGE analyses identified three ALPs showing a decreased expression in the LiTOX strain (AAEL003313 and two spots of AAEL003298). Moreover, the ALP AAEL003298 was also detected as under-transcribed in microarrays, RT-qPCR and in a previous transcriptomic study . These two ALP have already been described as Cry4Ba-binding proteins . The reduction of potential Cry-receptor ALPs proteins on the brush border of LiTOX larvae is consistent with the resistant phenotype.
N-Aminopeptidases (APNs) are a third major class of Cry toxin receptors in Lepidoptera  and mosquitoes [43, 44] and their alteration correlates with Cry1A-resistance in Helicoverpa armigera and Trichoplusia ni. DIGE experiments revealed two spots of APN AAEL012776 being under-expressed, congruent with transcriptomic data, while another spot was over-expressed in the LiTOX strain. Two spots matching APN AAEL012774 were over-expressed in the LiTOX strain, as also found by transcriptomic approaches. These two ALP proteins have been previously described as potential receptors for Cry11Aa in Ae. aegypti but it is still unclear how their altered expression could lead to a higher Bti-resistance.
In general, Bt resistance involves changes in the Cry receptors structure rather than in their expression [27, 47], although some cases of differential expression of cadherin and aminopeptidase have been reported in resistant strains [48–50]. These changes in expression can be the result of diverse genetic mechanisms including mutations in regulatory regions or even genome rearrangements that can drive rapid adaptation to new environmental pressures such as an insecticide treatment. Moreover, in the case of Bti, the presence of Cyt toxins, known to act as Cry receptors , might contribute to overcome receptor alterations in the LiTOX strain. Further analysis of the binding capacities of Cry toxins to the putative receptors found differentially expressed here will contribute to evaluate their relative roles in Bti resistance. Only few studies have focused on Cry4Aa toxin binding to our knowledge [52, 53] and nothing is known about its potential membrane receptors. Such experiments will determine if its receptors are highly specific, explaining the high differences in the resistance ratio between Cry4Aa and the other Cry toxins in the LiTOX strain, or if Cry4Aa shares all or a part of its receptors with Cry4Ba and Cry11Aa that could lead to cross-resistance.
Other mechanisms potentially involved in the resistance
All the spots of the four ATP synthases detected by our proteomic approach showed an under-expression pattern in the LiTOX strain. Vacuolar H+-ATPases (V-ATPase) subunits B to E are known to bind Cry4Ba  and subunits A and B have been described to bind for Cry1Ac in Heliothis virescens and Helicoverpa armigera. Moreover, V-ATPases are localized in the posterior midgut of mosquito larvae , where Cry4Aa, Cry4Ba and Cry11Aa toxins exhibit the highest affinity to the epithelium membrane [52, 53, 57, 58]. Nevertheless, their role as Bti toxins receptors has not been demonstrated yet. V-ATPases are strongly implicated in the alkalinization of the midgut pH by establishing a proton motive force by transporting proton across membranes leading to a pH gradient and transmembrane voltage [59–61]. Onken et al. (2008) inhibited all the proteins implicated in the alkalinization process in the midgut of Ae. aegypti larvae and they showed that only the inhibition of V-ATPases induced a strong acidification of the midgut pH . As pH affects numerous aspects of toxin action like Bt crystal solubility , toxin conformation [63, 64], gut enzymes activity  and pore formation [66, 67], an alteration of gut pH could have a general effect on reducing Bti toxicity. Comparing internal larval midgut pH between resistant and susceptible strains will allow to confirm/infirm if the observed ATPases decreased expression induce an acidification of the gut lumen.
Multiple detoxification enzymes were found under-transcribed in the resistant strain. Such enzymes are often involved in the degradation of small chemicals such as insecticides and plant allele-chemicals [68, 69], but they are unlikely to process large proteins such as Bti toxins. Synthesis of detoxification enzymes represents an important energetic cost for the insect . Moreover, several detoxification genes found under-transcribed in our dataset, were found over-transcribed in Ae. aegypti larvae submitted to a chemical challenge . Although the over-expression of particular detoxification genes in the resistant strain can be linked to larval response to tannins contained in the toxic leaf litter , the frequent under-expression of these enzymes in the resistant strain may reflect compensatory mechanisms.
Bti has evolved to infect Diptera such as mosquitoes and blackflies through a sequential mechanism. The multi-step mode of action of Bti and its toxins from ingestion to spore germination and proliferation offers many resistance ways for mosquito larvae. By combining transcriptomic and proteomic approaches, we detected expression alteration at nearly each step of the ingestion-to-infection process. Our study paves the way to further functional studies to characterize resistance mechanisms to this bioinsecticide. This information will be of extreme value as this environmentally safe bioinsecticide is increasingly used for vector control worldwide with virtually no knowledge and no suspicion so far about how mosquitoes can develop resistance in the field.
The Ae. aegypti laboratory strain Bora-Bora, susceptible to all insecticides, was used for selection with field-collected leaf litter containing Bti spores and toxins . This material, highly toxic after ingestion by mosquito larvae, was used for laboratory selection during 30 generations to obtain the LiTOX strain. Selection consisted in exposing 6000 third instar larvae to toxic leaf litter to obtain about 70% of larval mortality after 48 h exposure . Both susceptible and resistant strains were reared in standard insectary conditions (27°C, 14/10 h light/dark period, 80% relative humidity). Larvae were reared in tap water and fed with standard amount of larval food [19, 73].
Production of individual Bti Cry toxins
To produce Bti Cry toxins separately, we used a crystal negative strain of Bacillus thuringiensis var. israelensis (4Q2-81) transformed with the plasmids pHT606, pHT618 or pWF53 producing respectively Cry4A, Cry4B and Cry11 toxins obtained from the Pasteur Institute (Paris, France) or from Prof. B. Federici (University of Riverside, USA). Transformed Bti bacteria were grown on Nutrient Agar solid medium (Sigma Aldrich) supplemented with erythromycin antibiotic (25 μg/mL). Spores and crystals were recovered using cell scrapers (BD Falcon) after 7 days at 30°C and purified as previously described . This protocol ensures producing large amount of high quality toxin. Toxins were corun on SDS-PAGE with BSA at five concentrations (from 20, 40, 60, 80, 100 μg/mL). Intensity of each band was estimated and toxin concentration was calculated using BSA as standard using Imagej software v.1.41o .
Comparative bioassays between the LiTOX and the susceptible Bora-Bora strains were conducted after 30 generations of laboratory selection. Larvae from each strain were exposed to 6 concentrations of Cry4Aa, Cry4Ba, Cry11Aa and commercial Bti (Vectobac WG, 3500 ITU/mg) for 24 h to obtain 5% to 95% mortality. Bioassays were performed in triplicate on 20 third-instar larvae in 50 mL of insecticide solution or tap water (control) according to the standard bioassay procedure described by the World Health Organisation . LC50 (lethal concentration for 50% individuals) were calculated for each strain and each toxin using a probit statistical model with the module ‘dose’ of XLSTAT v.2009.4.06 (Addinsoft). For each toxin, resistance ratios (RR50) were calculated by dividing LC50 of the LiTOX strain by LC50 of the susceptible strain.
Larval midgut RNA extraction
For each strain, three biological replicates of 150 dissected midguts from early fourth instar larvae were prepared and conserved overnight at 4°C in RNAlater® (Ambion). After a brief centrifugation, supernatant was discarded and total RNA was extracted using RNAqueous®-4PCR kit (Ambion) following manufacturer’s instructions. Quantity and quality of RNA were assessed by spectrophotometry (Nanodrop ND-1000 spectrophotometer). To digest remaining genomic DNA, RNA samples were treated with DNAseI (Ambion) following manufacturer’s instructions. RNA were then concentrated using ammonium acetate and linear acrylamide to obtain at least 70 ng/μL of total RNA for each sample. Because the LiTOX strain was selected from the Bora-Bora strain, they share the same genetic background, and both were bred together in the same insectarium standard conditions, so that any constitutive change in gene expression between these strains is likely to result from Bti selection.
Larval midguts transcriptome profiling by DNA microarray
Low Input Quick Amp Labeling Kit, two-color (Agilent), containing Cy5 and Cy3 fluorescent dyes, was used to amplify and label messenger RNA. Labeled RNAs were then purified using Absolutely RNA® Nanoprep Kit (Stratagene) following manufacturer’s instructions with two elution steps in a final volume of 25 μL. Quantity and quality of RNA and labeling efficiency were assessed using Nanodrop spectrophotometer and Bioanalyzer® (Agilent).
Microarray hybridizations were performed with the 15 K Agilent ‘Aedes detox chip plus’ DNA microarray (ArrayExpress accession number A-MEXP-1966), containing eight replicated arrays of 60-mers oligo-probes representing 14204 different Ae. aegypti transcripts and several control probes. For each biological replicate, two hybridizations were performed in which the Cy3 and Cy5 labels were swapped between samples for a total of six hybridizations. For each hybridization, 300 ng of labeled mRNA were used. After 17 h hybridization, non-specific probes were washed off according to manufacturer’s instructions. Slides were scanned with an Agilent G2205B microarray scanner. Spot finding and signal quantification were performed using the Agilent Feature Extraction software (Agilent Technologies).
Data were analyzed using GeneSpring GX v9.0 software (Agilent). Only transcripts present in at least 5 hybridizations out of 6 were kept for further analyses. Transcripts exhibiting more than 3-fold transcription and a Benjamini-Hochberg  corrected P-value <0.01 were considered significantly differentially transcribed between the LiTOX and the susceptible strain. Midgut transcripts detected by microarrays were then classified into thirteen different categories based on their putative biological functions: receptors, transport, DNA interaction, cytoskeleton, ribosomal proteins, proteases, detoxification enzymes, kinases-phosphatases, transaminases, dehydrogenases, other enzymes, other function and unknown function. For genes of ‘unknown function’, the putative function was further investigated using BLASTP software, but they were not considered for functional analysis.
Real-time quantitative PCR (RT-qPCR) validation of microarray data
Transcription levels of 15 genes detected differentially transcribed with the microarray approach were validated by RT-qPCR using the same RNA extracts used in microarrays. In addition, transcription levels of two more genes (ALP1 and Cad1) encoding known Bti Cry toxins binding proteins were also compared between both strains by RT-qPCR. Three technical replicates were performed for each of the three biological replicates. Specific primers were designed for each gene using Beacon Designer v.5.10 software (Premier Biosoft International) (Additional file 8). Their specificity to the target gene was verified by BLAST analysis against Ae. aegypti genome. First-strand cDNA synthesis was obtained from 4 μg RNA by incubating them at 50°C for 1 h with SuperScript III (Invitrogen) reverse transcriptase, oligo-dT20 primers (2.5 μM), dNTPs (0.5 mM each), DTT (5 mM) and RNase Out (40 U, Invitrogen). Real-time quantitative PCR reaction was performed in 25 μL total reaction volume with specific primers (0.3 μM each), 12.5 μL iQ SYBR Green supermix (Bio-Rad) and 5 μL diluted cDNA on an iQ5 system (Bio-Rad). After an initial denaturing step at 95°C for 3 min, 40 cycles were performed each consisting in a denaturing step 15 s at 95°C and an annealing step 30 s at the optimal temperature of each primers couple (Additional file 8) [22, 71]. Specificity of DNA amplification was assessed by performing a melt curve analysis and verifying PCR product Tm. To check for any contamination, “no template controls” (NTC) were added in each PCR plate.
For each gene analyzed, a serial dilution of pooled cDNA from both strains was used to estimate PCR efficiency. Genes encoding ribosomal proteins RPL8 and RPS7 (housekeeping genes) were used for gene expression normalization taking into account PCR efficiency using ΔΔCT method, calculated using the iQ5 software (Bio-Rad) [77, 78]. Mean transcription ratios are expressed for the resistant strain relative to the susceptible strain.
GPI-anchor domain detection
To see whether some midgut enzymes detected with DNA microarrays were membrane-bound, we looked for glycosylphosphatidylinositol (GPI)-anchor domains using four complementary GPI domains predictors: big-PI Predictor v.3.0 (http://mendel.imp.ac.at/gpi/gpi_server.html) , PredGPI (http://gpcr.biocomp.unibo.it/predgpi/pred.htm) , FragAnchor (http://navet.ics.hawaii.edu/~fraganchor/NNHMM/NNHMM.html)  and GPI-SOM (http://gpi.unibe.ch/) .
Brush border membrane vesicles (BBMV) preparation
For each strain, two independent biological replicates were prepared. The day before midgut dissection, water was changed and food discarded. Early fourth instar larvae were chilled on ice for at least 20 min. Larvae were then dried on a clean paper. Midguts were dissected and mixed together in MET buffer (300 mM Mannitol, 5 mM EGTA, 17 mM TrisHCl, pH 7.5) with Complete Protease Inhibitor (Roche) to be conserved at −80°C until use. About 1500 larvae were dissected for each larval strain and biological replicate. 500 μg of midguts were centrifuged 5 min at 12,000 g to discard the old buffer, resuspended in ice-cold fresh MET buffer containing 1 mM PMSF and homogenized with 30 strokes of a glass-teflon homogenizer. BBMV were prepared following magnesium precipitation method as previously described . BBMV protein concentration was determined by a Bradford assay using BSA as standard . About 600 dissected guts yielded 500 μg of BBMV based on protein amount. Quality of BBMV was assessed by measuring the enrichment of two brush border enzymes: alkaline phosphatases (ALP) and aminopeptidases (APN). ALP and APN activities were measured using 4-nitrophenyl phosphate disodium and L-leucine-p-nitroanilide as substrates, respectively [85, 86]. APN and ALP enrichments are obtained by dividing the activity in the final BBMV preparation by the activity in the initial midgut homogenate (Table 5).
150 μg of BBMV proteins from each strain were used for each 2D-DIGE experiment and were purified using 2D-clean up kit (Amersham Bioscience) as described by the manufacturer. 100 μg of proteins were labeled with either Cy3 or Cy5 and the remaining 50 μg of proteins from the each strain were pooled and labeled with Cy2 as an internal standard. A dye swap was performed to be sure that the observed differences between the two strains were not due to different efficiencies of the dyes to label different proteins. The CyDye minimal labeling of the purified proteins was performed following manufacturer’s instructions (GE Healthcare). Labeled proteins were then mixed together and diluted to a final volume of 340 μL with rehydration buffer (2 M Thiourea, 7 M Urea, 3% CHAPS, 1% SB3-10, 13 mM DTT, 1% Immobilized pH Gradient (IPG) buffer pH 4–7, 0.002% Bromophenol blue (w/v)) and loaded on an IPG strip (pH 4–7 nonlinear, 18 cm) overlaid with 2 mL of plus-one IPG strip cover fluid (GE Healthcare). After 17 h of passive rehydration, the first dimension was run on a Multiphor-II flatbed system (GE Healthcare) at 20°C with the following program: 15 min at 300 V, 15 min at 500 V, and 9 h at 3500 V. This step allows proteins to migrate on the strip till a region in which pH is equal to their pI (isoelectric point).
After the Isoelectric Focusing, strips were reduced in equilibration buffer (6 M Urea, 75 mM Tris pH 8.8, 2%SDS, 29.3% Glycerol (v/v), 0.002% Bromophenol blue (w/v)) containing 1% of DTT (w/v) for 15 min and then alkylated in equilibration buffer with 2.5% of Iodoacetamide for 15 min. The IPG strip was then transferred on a pre-casted 12.5% SDS-PAGE gel (GE Healthcare) and second dimensional electrophoresis, separating proteins in function of their molecular size, was run at 22°C for 1 h at 2.5 W/gel followed by 5 h at 17 W/gel on an Ettan DALTsix vertical electrophoresis system (GE Healthcare). DIGE Gels were scanned using a Typhoon 9400 imager (GE Healthcare). As CyDye labeling induces a size modification of 1–2% of the amount of all the proteins that could bias the protein identification, non-labeled proteins were also prepared in parallel following the same protocol (except CyDye labeling) for mass spectra analyses and regular gels were co-run with DIGE gels to avoid modification in spot patterns due to different migrations. Regular gels were stained with Deep Purple stain (GE Healthcare), scanned using 532/610 nm excitation/emission wavelengths and used for spot picking.
2D-DIGE gels were analyzed using Decyder v7.0 software (GE Healthcare). The Decyder detection algorithm 5.0 was used to generate a list of spots with their coordinates and level of expression in the resistant strain relative to the susceptible strain. Only spots showing at least 1.5 fold differences between the two strains were considered for further analyses. 29 spots were picked from the first biological replicate and 35 from the second one, with 8 spots shared between them. Excised spots were digested with trypsin before subjecting peptides fragments to MALDI-ToF/ToF (time-of-flight) . To increase the likelihood of protein identification, each protein was identified by searching MS/MS data against an Ae. aegypti local database or other dipteran database when no significant match was obtained. To ensure accurate protein identification, we compared observed and expected pI values, molecular size, percentage of amino acid coverage and Mascot scores for Mascot search engine (http://www.matrixscience.com/search_intro.html) or z-scores for ProFound (http://prowl.rockefeller.edu).
Cadherin detection by immunoblotting
20 μg of proteins from BBMV prepared from the susceptible Bora-Bora strain, the LiTOX strain and the UGAL strain were separated by SDS-PAGE on 4–20% gradient TGX gels (Biorad). BBMV from the UGAL strain were prepared a few days before the experiments to compare the cadherin conservation in those fresh BBMV to the previously prepared BBMV from the two other strains. Proteins were either stained with coomassie blue to control that equal amount of proteins were stained from all the strains, or electroblotted to polyvinylidene fluoride (PVDF) filters for immunoblotting. Filters were blocked with 3% bovine serum albumin (BSA) in PBST (PBS + 0.1% Tween20) for 1 h at room temperature and then probed with α-AgCad1 antibodies (1:5000 dilution) , α-AgCad2 antibodies (1:500 dilution) or with pre-immune serum from the AgCad2 rabbit in PBST-0.1% BSA for 2 h. α-AgCad2 antibodies were prepared against an E. coli expressed cadherin peptide AgCad2 (Hua et al., unpublished work). Filters were then washed and detected by an anti-rabbit IgG-peroxidase conjugate (1:25000 dilution) in PBST-0.1% BSA for 1 h at room temperature. Filters were developed with an ECL kit (GE Healthcare) and chemiluminescence was detected with a ChemiImager (Alpha Innotech). All the Western blots were performed in duplicate.
Larval midgut proteolytic activities
For each strain, three biological replicates of soluble protein extracted from midgut juice were prepared. 20 midguts of early fourth instars were extracted and placed into 50 μL of distilled water and homogenized using a vortex for 30 s. Sample were centrifuged at 12,000 g for 10 min at 4°C. All the supernatants from larvae of the same biological replicate were mixed together, protein concentration was quantified by a Bradford assay using BSA as standard  and aliquots were conserved at −20°C until use. Total protease activity was measured using azocasein as substrate (Sigma Aldrich) as described in . All activities were normalized according to the amount of total protein from each replicate. For each biological replicate, six technical replicates were performed and absorbance was measured at 440 nm. Percentages of protease activity due to serine proteases, chymotrypsins, trypsins and metallo-enzymes were measured using respectively PMSF (30 mM), TPCK (1.5 mM), TLCK (1.5 mM) and EDTA (1 mM) (Sigma Aldrich) . Statistical differences between the two strains were measured by a Wilcoxon test performed with R 2.8.1 software .
J. P.D., M.J.A. and L.D. designed research; G.T. and K.B. performed DiGE experiments; K.B. did Western blot experiments; G.T., K.B. and M.J.A. analyzed proteomic data; G.T. and M.A.R. prepared mRNA samples for microarrays and did gene expression analyses; C.M.J. performed microarray experiment; G.T. performed RT-qPCR experiments; G.T. and J.P.D analyzed transcriptomic data; G.T. and R.S. did enzymatic experiments; G.T. and M.P. performed the selection and rearing of the LiTOX strain; G.T. wrote the paper; K.B., R.S., J.P.D., M.J.A. and L.D. reviewed and helped improving the manuscript. All authors read and approved the manuscript.
G.T. performed this study during his Ph.D. at the Laboratoire d’Ecologie Alpine (LECA), University of Grenoble. His research interests goes from the understanding of the fate of pesticides in the field to the adaptive strategies of target insects. K.B. worked on this study during his Ph.D. tenure in the Dept. of Entomology, University of Georgia. Currently, he is a Postdoctoral Research Associate in the Dept. of Microbiology, UGA. His research interest is Oxidative protein damage and protein repair in Helicobacter pylori. C.M.J. is a post-doctoral researcher in the Vector Group of the LSTM where his current interests lie in vector biology and genetics with a strong emphasis on insecticide resistance. R.S. participated to this study during his Master internship at the LECA. He is interested in molecular ecology and agricultural science. M.P. participated to this study during her Ph.D. She is now a post-doc in Zurich, and is interested in population genetics, molecular adaptation and evolution. J.P.D. is a senior CNRS researcher having a strong experience in vector control and adaptive mechanisms developed by insects to insecticides, pollutants and toxins. M.J.A., a professor at the University of Georgia and CSO of InsectiGen, investigates Bt toxins and their action in pest insects. L.D. is a Professor at University of Grenoble. Her main research interests are in adaptive patterns in natural populations and underlying evolutionary processes.
Brush border membrane vesicles
Differential in gel electrophoresis
Tandem mass spectrometry
Nα-tosyl-L-lysine chloromethyl ketone hydrochloride
N-p-Tosyl-L-phenylalanine chloromethyl ketone.
Hemingway J, Ranson H: Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000, 45: 371-391. 10.1146/annurev.ento.45.1.371.
Thai KTD, Anders KL: The role of climate variability and change in the transmission dynamics and geographic distribution of dengue. Exp Biol Med. 2011, 236 (8): 944-954. 10.1258/ebm.2011.010402.
Weaver SC, Reisen WK: Present and future arboviral threats. Antivir Res. 2010, 85 (2): 328-345. 10.1016/j.antiviral.2009.10.008.
Kyle JL, Harris E: Global spread and persistence of dengue. Annu Rev Microbiol. 2008, 62: 71-92. 10.1146/annurev.micro.62.081307.163005.
Tomori O: Yellow fever: the recurring plague. Crit Rev Clin Lab Sci. 2004, 41 (4): 391-427. 10.1080/10408360490497474.
van den Berg H: Global status of DDT and its alternatives for use in vector control to prevent disease. Cienc Saude Coletiva. 2011, 16 (2): 575-590. 10.1590/S1413-81232011000200021.
Lacey LA: Bacillus thuringiensis serovariety israelensis and Bacillus sphaericus for mosquito control. J Am Mosq Control Assoc. 2007, 23 (2): 133-163. 10.2987/8756-971X(2007)23[133:BTSIAB]2.0.CO;2.
Bravo A, Gill SS, Soberon M: Mode of action of Bacillus thuringiensis Cry and Cyt toxins and their potential for insect control. Toxicon. 2007, 49 (4): 423-435. 10.1016/j.toxicon.2006.11.022.
Bravo A, Likitvivatanavong S, Gill SS, Soberon M: Bacillus thuringiensis: a story of a successful bioinsecticide. Insect Biochem Mol Biol. 2011, 41 (7): 423-431. 10.1016/j.ibmb.2011.02.006.
Gassmann AJ, Petzold-Maxwell JL, Keweshan RS, Dunbar MW: Field-evolved resistance to Bt maize by Western corn rootworm. PLoS One. 2011, 6 (7): e22629-10.1371/journal.pone.0022629.
Dhurua S, Gujar GT: Field-evolved resistance to Bt toxin Cry1Ac in the pink bollworm, Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae), from India. Pest Manag Sci. 2011, 67 (8): 898-903. 10.1002/ps.2127.
Tabashnik BE, Van Rensburg JBJ, Carriere Y: Field-evolved insect resistance to Bt crops: definition, theory, and data. J Econ Entomol. 2009, 102 (6): 2011-2025. 10.1603/029.102.0601.
Storer NP, Babcock JM, Schlenz M, Meade T, Thompson GD, Bing JW, Huckaba RM: Discovery and characterization of field resistance to Bt maize: Spodoptera frugiperda (Lepidoptera: Noctuidae) in Puerto Rico. J Econ Entomol. 2010, 103 (4): 1031-1038. 10.1603/EC10040.
Paul A, Harrington LC, Zhang L, Scott JG: Insecticide resistance in Culex pipiens from New York. J Am Mosq Control Assoc. 2005, 21 (3): 305-309. 10.2987/8756-971X(2005)21[305:IRICPF]2.0.CO;2.
Wirth MC, Park HW, Walton WE, Federici BA: Cyt1A of Bacillus thuringiensis delays evolution of resistance to Cry11A in the mosquito Culex quinquefasciatus. Appl Environ Microbiol. 2005, 71 (1): 185-189. 10.1128/AEM.71.1.185-189.2005.
Tilquin M, Paris M, Reynaud S, Despres L, Ravanel P, Geremia RA, Gury J: Long lasting persistence of Bacillus thuringiensis Subsp israelensis (Bti) in mosquito natural habitats. PLoS One. 2008, 3 (10): e3432-10.1371/journal.pone.0003432.
de Melo-Santos MAV, de Araujo AP, Rios EMM, Regis L: Long lasting persistence of Bacillus thuringiensis serovar. israelensis larvicidal activity in Aedes aegypti (Diptera: Culicidae) breeding places is associated to bacteria recycling. Biol Control. 2009, 49 (2): 186-191. 10.1016/j.biocontrol.2009.01.011.
Shaheen R, Svensson B, Andersson MA, Christiansson A, Salkinoja-Salonen M: Persistence strategies of Bacillus cereus spores isolated from dairy silo tanks. Food Microbiol. 2010, 27 (3): 347-355. 10.1016/j.fm.2009.11.004.
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-128. 10.1002/ps.2046.
Goldman IF, Arnold J, Carlton BC: Selection for resistance to Bacillus thuringiensis subspecies israelensis in field and laboratory populations of the mosquito Aedes aegypti. J Invertebr Pathol. 1986, 47 (3): 317-324. 10.1016/0022-2011(86)90102-3.
Mittal PK, Adak T, Subbarao SK: Laboratory selection to investigate the development of resistance to Bacillus thuringiensis var. israelensis H-14 in Culex quinquefasciatus Say (Diptera : Culicidae). Natl Acad Sci Lett-India. 2005, 28 (7-8): 281-283.
Bonin A, Paris M, Tetreau G, David JP, Despres L: Candidate genes revealed by a genome scan for mosquito resistance to a bacterial insecticide: sequence and gene expression variations. BMC Genomics. 2009, 10: 551-10.1186/1471-2164-10-551.
Paris M, Melodelima C, Coissac E, Tetreau G, Reynaud S, David JP, Despres L: Transcription profiling of resistance to Bti toxins in the mosquito Aedes aegypti using next-generation sequencing. J Invertebr Pathol. 2012, 109 (2): 201-208. 10.1016/j.jip.2011.11.004.
Paris M, David JP, Despres L: Fitness costs of resistance to Bti toxins in the dengue vector Aedes aegypti. Ecotoxicology. 2011, 20: 1184-1194. 10.1007/s10646-011-0663-8.
Georghiou GP, Wirth MC: Influence of exposure to single versus multiple toxins of Bacillus thuringiensis subsp. israelensis on development of resistance in the mosquito Culex quinquefasciatus (Diptera: Culicidae). Appl Environ Microbiol. 1997, 63 (3): 1095-1101.
Saleh MS, El-Meniawi FA, Kelada NL, Zahran HM: Resistance development in mosquito larvae Culex pipiens to the bacterial agent Bacillus thuringiensis var. israelensis. J Appl Entomol-Zeitschrift Fur Angew Entomol. 2003, 127 (1): 29-32.
Griffitts JS, Aroian RV: Many roads to resistance: how invertebrates adapt to Bt toxins. BioEssays. 2005, 27 (6): 614-624. 10.1002/bies.20239.
Popova-Butler A, Dean DH: Proteomic analysis of the mosquito Aedes aegypti midgut brush border membrane vesicles. J Insect Physiol. 2009, 55 (3): 264-272. 10.1016/j.jinsphys.2008.12.008.
Bayyareddy K, Andacht TM, Abdullah MA, Adang MJ: Proteomic identification of Bacillus thuringiensis subsp israelensis toxin Cry4Ba binding proteins in midgut membranes from Aedes (Stegomyia) aegypti Linnaeus (Diptera, Culicidae) larvae. Insect Biochem Mol Biol. 2009, 39 (4): 279-286. 10.1016/j.ibmb.2009.01.002.
Gygi SP, Rochon Y, Franza BR, Aebersold R: Correlation between protein and mRNA abundance in yeast. Mol Cell Biol. 1999, 19 (3): 1720-1730.
De Wit M, Keil D, van der Ven K, Vandamme S, Witters E, De Coen W: An integrated transcriptomic and proteomic approach characterizing estrogenic and metabolic effects of 17 alpha-ethinylestradiol in zebrafish (Danio rerio). Gen Comp Endocrinol. 2010, 167 (2): 190-201. 10.1016/j.ygcen.2010.03.003.
Foss EJ, Radulovic D, Shaffer SA, Goodlett DR, Kruglyak L, Bedalov A: Genetic variation shapes protein networks mainly through non-transcriptional mechanisms. PLoS Biol. 2011, 9 (9): e1001144-10.1371/journal.pbio.1001144.
Likitvivatanavong S, Chen JW, Evans AM, Bravo A, Soberon M, Gill SS: Multiple receptors as targets of cry toxins in mosquitoes. J Agric Food Chem. 2011, 59 (7): 2829-2838. 10.1021/jf1036189.
Chen JW, Aimanova KG, Fernandez LE, Bravo A, Soberon M, Gill SS: Aedes aegypti cadherin serves as a putative receptor of the Cry11Aa toxin from Bacillus thuringiensis subsp israelensis. Biochem J. 2009, 424: 191-200. 10.1042/BJ20090730.
Chen JW, Aimanova KG, Pan SQ, Gill SS: Identification and characterization of Aedes aegypti aminopeptidase N as a putative receptor of Bacillus thuringiensis Cry11A toxin. Insect Biochem Mol Biol. 2009, 39 (10): 688-696. 10.1016/j.ibmb.2009.08.003.
Fernandez LE, Aimanova KG, Gill SS, Bravo A, Soberon M: A GPI-anchored alkaline phosphatase is a functional midgut receptor of Cry11Aa toxin in Aedes aegypti larvae. Biochem J. 2006, 394: 77-84. 10.1042/BJ20051517.
McNall RJ, Adang MJ: Identification of novel Bacillus thuringiensis Cry1Ac binding proteins in Manduca sexta midgut through proteomic analysis. Insect Biochem Mol Biol. 2003, 33 (10): 999-1010. 10.1016/S0965-1748(03)00114-0.
Arenas I, Bravo A, Soberon M, Gomez I: Role of alkaline phosphatase from Manduca sexta in the mechanism of action of Bacillus thuringiensis Cry1Ab toxin. J Biol Chem. 2010, 285 (17): 12497-12503. 10.1074/jbc.M109.085266.
Martins ES, Monnerat RG, Queiroz PR, Dumas VF, Braz SV, Aguiar RWD, Gomes A, Sanchez J, Bravo A, Ribeiro BM: Midgut GPI-anchored proteins with alkaline phosphatase activity from the cotton boll weevil (Anthonomus grandis) are putative receptors for the Cry1B protein of Bacillus thuringiensis. Insect Biochem Mol Biol. 2010, 40 (2): 138-145. 10.1016/j.ibmb.2010.01.005.
Hua G, Zhang R, Bayyareddy K, Adang MJ: Anopheles gambiae alkaline phosphatase is a functional receptor of Bacillus thuringiensis jegathesan Cry11Ba toxin. Biochemistry. 2009, 48 (41): 9785-9793. 10.1021/bi9014538.
Jurat-Fuentes JL, Karumbaiah L, Jakka SRK, Ning CM, Liu CX, Wu KM, Jackson J, Gould F, Blanco C, Portilla M: Reduced levels of membrane-bound alkaline phosphatase are common to Lepidopteran strains resistant to cry toxins from Bacillus thuringiensis. PLoS One. 2011, 6 (3): e17606-10.1371/journal.pone.0017606.
Pigott CR, Ellar DJ: Role of receptors in Bacillus thuringiensis crystal toxin activity. Microbiol Mol Biol Rev. 2007, 71 (2): 255-10.1128/MMBR.00034-06.
Abdullah MAF, Moussa S, Taylor MD, Adang MJ: Manduca sexta (Lepidoptera: Sphingidae) cadherin fragments function as synergists for Cry1A and Cry1C Bacillus thuringiensis toxins against noctuid moths Helicoverpa zea, Agrotis ipsilon and Spodoptera exigua. Pest Manag Sci. 2009, 65 (10): 1097-1103. 10.1002/ps.1798.
Saengwiman S, Aroonkesorn A, Dedvisitsakul P, Sakdee S, Leetachewa S, Angsuthanasombat C, Pootanakit K: In vivo identification of Bacillus thuringiensis Cry4Ba toxin receptors by RNA interference knockdown of glycosylphosphatidylinositol-linked aminopeptidase N transcripts in Aedes aegypti larvae. Biochem Biophys Res Commun. 2011, 407 (4): 708-713. 10.1016/j.bbrc.2011.03.085.
Zhang SP, Cheng HM, Gao YL, Wang GR, Liang GM, Wu KM: Mutation of an aminopeptidase N gene is associated with Helicoverpa armigera resistance to Bacillus thuringiensis Cry1Ac toxin. Insect Biochem Mol Biol. 2009, 39 (7): 421-429. 10.1016/j.ibmb.2009.04.003.
Tiewsiri K, Wang P: Differential alteration of two aminopeptidases N associated with resistance to Bacillus thuringiensis toxin Cry1Ac in cabbage looper. Proc Natl Acad Sci U S A. 2011, 108 (34): 14037-14042. 10.1073/pnas.1102555108.
Ferre J, Escriche B, Bel Y, Vanrie J: Biochemistry and genetics of insect resistance to Bacillus-thuringiensis insecticidal crystal proteins. FEMS Microbiol Lett. 1995, 132 (1–2): 1-7.
Yang YL, Zhu YC, Ottea J, Husseneder C, Leonard BR, Abel C, Luttrell R, Huang FN: Down regulation of a gene for cadherin, but not alkaline phosphatase, associated with Cry1Ab resistance in the sugarcane borer Diatraea saccharalis. PLoS One. 2011, 6 (10): e25783-10.1371/journal.pone.0025783.
Herrero S, Gechev T, Bakker PL, Moar WJ, de Maagd RA: Bacillus thuringiensis Cry1Ca-resistant Spodoptera exigua lacks expression of one of four Aminopeptidase N genes. BMC Genomics. 2005, 6 (96):
Jurat-Fuentes JL, Adang MJ: Characterization of a Cry1Ac-receptor alkaline phosphatase in susceptible and resistant Heliothis virescens larvae. Eur J Biochem. 2004, 271 (15): 3127-3135. 10.1111/j.1432-1033.2004.04238.x.
Perez C, Fernandez LE, Sun JG, Folch JL, Gill SS, Soberon M, Bravo A: Bacillus thuringiensis subsp israelensis Cyt1Aa synergizes Cry11Aa toxin by functioning as a membrane-bound receptor. Proc Natl Acad Sci U S A. 2005, 102 (51): 18303-18308. 10.1073/pnas.0505494102.
Ravoahangimalala O, Charles JF: In-vitro binding of Bacillus-thuringiensis var israelensis individual toxins to midgut cells of Anopheles-gambiae larvae (Diptera, Culicidae). FEBS Lett. 1995, 362 (2): 111-115. 10.1016/0014-5793(95)00220-4.
Yamagiwa M, Kamauchi S, Okegawa T, Esaki M, Otake K, Amachi T, Komano T, Sakai H: Binding properties of Bacillus thuringiensis Cry4A toxin to the apical microvilli of larval midgut of Culex pipiens. Biosci Biotechnol Biochem. 2001, 65 (11): 2419-2427. 10.1271/bbb.65.2419.
Krishnamoorthy M, Jurat-Fuentes JL, McNall RJ, Andacht T, Adang MJ: Identification of novel CrylAc binding proteins in midgut membranes from Heliothis virescens using proteomic analyses. Insect Biochem Mol Biol. 2007, 37 (3): 189-201. 10.1016/j.ibmb.2006.10.004.
Chen LZ, Liang GM, Zhang J, Wu KM, Guo YY, Rector BG: Proteomic analysis of novel Cry1Ac binding proteins in Helicoverpa armigera (Hubner). Arch Insect Biochem Physiol. 2010, 73 (2): 61-73. 10.1002/arch.20340.
Harvey WR, Okech BA, Linser PJ, Becnel JJ, Ahearn GA, Sterling KM: H+ V-ATPase-energized transporters in brush border membrane vesicles from whole larvae of Aedes aegypti. J Insect Physiol. 2010, 56 (10): 1377-1389. 10.1016/j.jinsphys.2010.04.017.
Ravoahangimalala O, Charles JF, Schoellerraccaud J: Immunological localization of Bacillus-thuringiensis serovar israelensis toxins in midgut cells of intoxicated Anopheles-gambiae larvae (Diptera, Culicidae). Res Microbiol. 1993, 144 (4): 271-278. 10.1016/0923-2508(93)90011-P.
Ruiz LM, Segura C, Trujillo J, Orduz S: In vivo binding of the Cry11Bb toxin of Bacillus thuringiensis subsp medellin to the midgut of mosquito larvae (Diptera : Culicidae). Memorias Do Inst Oswaldo Cruz. 2004, 99 (1): 73-79. 10.1590/S0074-02762004000100013.
Onken H, Moffett DF: Revisiting the cellular mechanisms of strong luminal alkalinization in the anterior midgut of larval mosquitoes. J Exp Biol. 2009, 212 (3): 373-377. 10.1242/jeb.023580.
Boudko DY, Moroz LL, Linser PJ, Trimarchi JR, Smith PJS, Harvey WR: In situ analysis of pH gradients in mosquito larvae using noninvasive, self-referencing, pH-sensitive microelectrodes. J Exp Biol. 2001, 204 (4): 691-699.
Beyenbach KW, Wieczorek H: The V-type H+ ATPase: molecular structure and function, physiological roles and regulation. J Exp Biol. 2006, 209 (4): 577-589. 10.1242/jeb.02014.
Onken H, Moffett SB, Moffett DF: Alkalinization in the isolated and perfused anterior midgut of the larval mosquito, Aedes aegypti. J Insect Sci (Online). 2008, 8: 1-20.
Feng Q, Becktel WJ: pH-induced conformational transitions of Cry IA(a), Cry IA(c), and Cry IIIA delta-endotoxins in Bacillus-thuringiensis. Biochemistry. 1994, 33 (28): 8521-8526. 10.1021/bi00194a017.
Venugopal MG, Wolfersberger MG, Wallace BA: Effects of pH on conformational properties related to the toxicity of Bacillus-thuringiensis delta-endotoxin. Biochim Biophys Acta. 1992, 1159 (2): 185-192. 10.1016/0167-4838(92)90024-8.
Karumbaiah L, Oppert B, Jurat-Fuentes JL, Adang MJ: Analysis of midgut proteinases from Bacillus thuringiensis-susceptible and -resistant Heliothis virescens (Lepidoptera : Noctuidae). Comp Biochem Physiol B-Biochem Mol Biol. 2007, 146 (1): 139-146. 10.1016/j.cbpb.2006.10.104.
Tran LB, Vachon V, Schwartz JL, Laprade R: Differential effects of pH on the pore-forming properties of Bacillus thuringiensis insecticidal crystal toxins. Appl Environ Microbiol. 2001, 67 (10): 4488-4494. 10.1128/AEM.67.10.4488-4494.2001.
Fortier M, Vachon V, Kirouac M, Schwartz JL, Laprade R: Differential effects of ionic strength, divalent cations and pH on the pore-forming activity of Bacillus thuringiensis insecticidal toxins. J Membr Biol. 2005, 208 (1): 77-87. 10.1007/s00232-005-0820-1.
Despres L, David JP, Gallet C: The evolutionary ecology of insect resistance to plant chemicals. Trends Ecol Evol. 2007, 22 (6): 298-307. 10.1016/j.tree.2007.02.010.
Feyereisen R: Insect CYP genes and P450 enzymes. Insect Molecular Biology and Biochemistry. Edited by: Gilbert LI. 2011, Academic Press , 236-316.
Zeng RS, Wen Z, Niu G, Schuler MA, Berenbaum MR: Enhanced toxicity and induction of cytochrome P450s suggest a cost of “eavesdropping” in a multitrophic interaction. J Chem Ecol. 2009, 35 (5): 526-532. 10.1007/s10886-009-9640-6.
Poupardin R, Reynaud S, Strode C, Ranson H, Vontas J, David JP: Cross-induction of detoxification genes by environmental xenobiotics and insecticides in the mosquito Aedes aegypti: impact on larval tolerance to chemical insecticides. Insect Biochem Mol Biol. 2008, 38 (5): 540-551. 10.1016/j.ibmb.2008.01.004.
David JP, Boyer S, Mesneau A, Ball A, Ranson H, Dauphin-Villemant C: Involvement of cytochrome P450 monooxygenases in the response of mosquito larvae to dietary plant xenobiotics. Insect Biochem Mol Biol. 2006, 36 (5): 410-420. 10.1016/j.ibmb.2006.02.004.
Riaz MA, Poupardin R, Reynaud S, Strode C, Ranson H, David JP: Impact of glyphosate and benzo[a]pyrene on the tolerance of mosquito larvae to chemical insecticides. Role of detoxification genes in response to xenobiotics. Aquat Toxicol. 2009, 93 (1): 61-69. 10.1016/j.aquatox.2009.03.005.
Girish V, Vijayalakshmi A: Affordable image analysis using NIH Image/ImageJ. Indian J Cancer. 2004, 41 (1): 47-
WHO: Guidelines for Laboratory and Field Testing of Mosquito Larvicides. 2005, Geneva: World Health Organization, WHO/CDS/WHOPES/GCDPP/2005.13
Benjamini Y, Hochberg Y: Controlling the false discovery rate—a practical and powerful approach to multiple testing. J R Stat Soc Ser B-Methodol. 1995, 57 (1): 289-300.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29 (9): e45-10.1093/nar/29.9.e45.
Poupardin R, Riaz MA, Vontas J, David JP, Reynaud S: Transcription profiling of eleven cytochrome P450s potentially involved in xenobiotic metabolism in the mosquito Aedes aegypti. Insect Mol Biol. 2010, 19 (2): 185-193. 10.1111/j.1365-2583.2009.00967.x.
Eisenhaber B, Bork P, Yuan YP, Loffler G, Eisenhaber F: Automated annotation of GPI anchor sites: case study C-elegans. Trends in Biochem Sci. 2000, 25 (7): 340-341. 10.1016/S0968-0004(00)01601-7.
Pierleoni A, Martelli PL, Casadio R: PredGPI: a GPI-anchor predictor. BMC Bioinforma. 2008, 9:
Poisson G, Chauve C, Chen X, Bergeron A: FragAnchor: a large-scale predictor of glycosylphosphatidylinositol anchors in eukaryote protein sequences by qualitative scoring. Genomics, Proteomics bioinforma / Beijing Genomics Inst. 2007, 5 (2): 121-130.
Fankhauser N, Maser P: Identification of GPI anchor attachment signals by a Kohonen self-organizing map. Bioinformatics. 2005, 21 (9): 1846-1852. 10.1093/bioinformatics/bti299.
Abdul-Rauf M, Ellar DJ: Isolation and characterization of brush border membrane vesicles from whole Aedes aegypti larvae. J Invertebr Pathol. 1999, 73 (1): 45-51. 10.1006/jipa.1998.4792.
Bradford MM: Rapid and sensitive method for quantitation of microgram quantities of protein utilizing principle of protein-dye binding. Anal Biochem. 1976, 72 (1–2): 248-254.
Terra WR, Ferreira C: Insect digestive enzymes—properties, compartmentalization and function. Comp Biochem Physiol B-Biochem Mol Biol. 1994, 109 (1): 1-62.
Hernandez-Martinez P, Navarro-Cerrillo G, Caccia S, de Maagd RA, Moar WJ, Ferre J, Escriche B, Herrero S: Constitutive activation of the midgut response to Bacillus thuringiensis in Bt-resistant Spodoptera exigua. PLoS One. 2010, 5 (9): e12795-10.1371/journal.pone.0012795.
Hua G, Zhang R, Abdullah MAF, Adang MJ: Anopheles gambiae cadherin AgCad1 binds the Cry4Ba toxin of Bacillus thuringiensis israelensis and a fragment of AgCad1 synergizes toxicity. Biochemistry. 2008, 47 (18): 5101-5110. 10.1021/bi7023578.
Hosseininaveh V, Bandani A, Azmayeshfard P, Hosseinkhani S, Kazzazi M: Digestive proteolytic and amylolytic activities in Trogoderma granarium Everts (Dermestidae : Coleoptera). J Stored Prod Res. 2007, 43: 515-522. 10.1016/j.jspr.2007.02.003.
R Development Core Team: R: A Language and Environment for Statistical Computing. 2007,http://www.R-project.org,
The authors would like to thank B. Federici for providing recombinant Bt strains. They also thank E. Bianco, A. Bonin and A. Proust for their help for larval midgut dissections. They thank S. Veyrenc and T. Gaude for technical assistance with mosquito breeding and insectary maintenance. They thank Dr. B. MacCallum for useful help with microarray data deposition in Vectorbase. This study was founded by the French National Research Agency (ANR, project ANR-08-CES-006-01 DIBBECO). Guillaume Tetreau was supported by the French Ministry of Research.
The authors declare that no competing interests exist.
Guillaume Tetreau, Krishnareddy Bayyareddy contributed equally to this work.
Electronic supplementary material
Additional file 1: All the 3512 transcripts detected by microarrays experiments in at least 5 hybridizations out of 6. For each transcript, accession number, corrected p-value, expression level changes, Vectorbase annotation and functional category are indicated. (XLS 570 KB)
Additional file 2: 70 transcripts significantly (corrected P-val<0.01) more than 3-fold differentially transcribed in the LiTOX strain. Transcripts are classified according to their putative function using the 13 functional categories. For each transcript, accession number, corrected P-value, expression level changes, Vectorbase annotation and supercontig are indicated. For transcripts of ‘unknown functions’, their putative function with corresponding score, ID, accession number and species of the best hit found using BLASTP software are indicated. (XLS 156 KB)
Additional file 3: Validation of microarray data by RT-qPCR on fifteen selected genes. Both experiments were performed on the same mRNA extracted from dissected larval midguts. ALP2, Alkaline phosphatase AAEL003298; ALP3, AAEL003313; ALP5, AAEL015070; ALP6, AAEL011175; APN1, N-Aminopeptidase AAEL012774; APN2, AAEL012776; APN3, AAEL012778; Cad2, Cadherin AAEL007488; HP1, Conserved hypothetical protein AAEL010435; HP2, AAEL013584; SE1, Serine-type endopeptidase AAEL007938; SE2, Serine-type endopeptidase AAEL011917; Cytochrome P450: CYP6Z7, AAEL009130; CYP6Z8, AAEL009131 and CYP4D24, AAEL007815. (TIFF 219 KB)
Additional file 4: 2D-DIGE gels from the two biological replicates and dye-swapping. BBMV prepared from first (A and B) and second (C and D) biological replicate are separated in function of their size (kDa) and their isoelectric point (pI). BBMV from Bti resistant strain are labeled with Cy3 and susceptible strain with Cy5 (A and C) or resistant strain with Cy5 and susceptible with Cy3 (B and D). (TIFF 2 MB)
Additional file 5: Protein identification of the 56 spots picked on deep purple stained 2D-gel. When different spots pointed to the same protein, they were differentiated using arbitrary letters after the access number. For each identification, the predicted pI, the predicted mass in kilodaltons, the percentage of sequence coverage, their functional category, and the species and database matched are indicated. (XLS 42 KB)
Additional file 6: Glycosylphosphatidylinositol (GPI)-anchor domains detection by four predictive computational programs. For each gene and protein, their accession number, the transcript and protein sizes are indicated. Results from the big-GPI and GPI-SOM softwares are indicated as ‘YES’ when they found a potential GPI-domain and ‘NO’ when no GPI-domain was determined. For PredGPI, presence is indicated by ‘Highly probable’, ‘Weakly probable’ or ‘Probable’ and absence by ‘NO’. For FragAnchor, presence of GPI domain is indicated by ‘Highly probable’ or ‘Probable’, absence by ‘NO’ and when prediction is uncertain by ‘Potential false positive’. (XLS 20 KB)
Additional file 7: Cadherin detection by immunoblotting. BBMV proteins from the susceptible Bora-Bora strain (lane 1), LiTOX strain (lane 2) and the UGAL Aedes strain (lane 3) were separated in SDS-PAGE and stained with coomassie blue (panel A) or probed with α-AgCad1 antibodies (panel B), α-AgCad2 antibodies (panel C) or with pre-immune serum from α-AgCad2 rabbit (panel D). (TIFF 3 MB)
Additional file 8: Primer pairs used for RT-qPCR analyses. For each primer pair, sequence, corresponding gene name and accession number, product length, Tm and optimal annealing temperature used in PCR program are indicated. PCR efficiency and different parameters of the calibration curves (R2, slope and y-intercept) are also indicated. Specificity of each primer pair was first assessed by BLAST analysis against Ae. aegypti genome and then verified by performing a melt curve analysis. A high specificity is indicated as “YES” when the primer pair matched to a unique position in the Ae. aegypti genome and when PCR product Tm was correct. (XLS 148 KB)
About this article
Cite this article
Tetreau, G., Bayyareddy, K., Jones, C.M. et al. Larval midgut modifications associated with Bti resistance in the yellow fever mosquito using proteomic and transcriptomic approaches. BMC Genomics 13, 248 (2012). https://doi.org/10.1186/1471-2164-13-248