Larval midgut modifications associated with Bti resistance in the yellow fever mosquito using proteomic and transcriptomic approaches
- Guillaume Tetreau†1Email author,
- Krishnareddy Bayyareddy†2,
- Christopher M Jones3,
- Renaud Stalinski1,
- Muhammad A Riaz1,
- Margot Paris1,
- Jean-Philippe David1,
- Michael J Adang2, 4 and
- Laurence Després1
© Tetreau et al.; licensee BioMed Central Ltd. 2012
Received: 13 March 2012
Accepted: 25 May 2012
Published: 15 June 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.
KeywordsAedes aegypti Bacillus thuringiensis israelensis DIGE Microarray RT-qPCR Resistance Transcriptomics Proteomics Midgut Mosquito Candidate genes
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
Lethal concentrations and resistance ratio for the LiTOX and susceptible strains for Bti and Cry toxins
LC50in ng/mL (95% CI)
Bti Vectobac WG
Midgut transcriptome profiling
Midgut differential proteomics
Genes and proteins differentially expressed in the LiTOX strain
Protein identification of 30 spots with highest Mascot scores picked on deep purple stained 2D-gel
Vectorbase access number
Top ranking match
Predicted mass (kDa)
% sequence coverage
protease m1 zinc metalloprotease
protease m1 zinc metalloprotease
protease m1 zinc metalloprotease
protease m1 zinc metalloprotease
protease m1 zinc metalloprotease
ATP synthase beta subunit
ATP synthase beta subunit
V-ATP synthase subunit alpha
V-ATP synthase subunit alpha
V-ATP synthase subunit alpha
V-ATP synthase subunit beta
V-ATP synthase subunit beta
V-ATP synthase subunit beta
V-ATP synthase subunit beta
V-ATP synthase subunit beta
V-ATP synthase subunit beta
V-ATP synthase subunit E
V-ATP synthase subunit E
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
Total enzymatic activity and effect of protease inhibitors on the azocaseinolytic activity of midgut extract
Total enzymatic activity (OD at 440 nm)
Percentage of inhibition of total activity (%)
0.328 ± 0.009
92.98 ± 0.57
46.82 ± 1.02
23.42 ± 2.24
12.18 ± 1.74
0.356 ± 0.010
92.31 ± 0.23
42.01 ± 2.29
20.87 ± 1.64
−1.74 ± 4.55
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
Altered expression of known Bti Cry-binding proteins detected by transcriptomic and proteomic approaches
Alkaline phosphatase (ALP1)
Alkaline phosphatase (ALP2)
Alkaline phosphatase (ALP3)
Alkaline phosphatase (ALP4)
Alkaline phosphatase (ALP5)
+1.81; −2.74; −3.81
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
APN and ALP enrichments in final BBMV preparation relative to the initial midgut homogenate
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.
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.
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Weaver SC, Reisen WK: Present and future arboviral threats. Antivir Res. 2010, 85 (2): 328-345. 10.1016/j.antiviral.2009.10.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Kyle JL, Harris E: Global spread and persistence of dengue. Annu Rev Microbiol. 2008, 62: 71-92. 10.1146/annurev.micro.62.081307.163005.View ArticlePubMedGoogle Scholar
- Tomori O: Yellow fever: the recurring plague. Crit Rev Clin Lab Sci. 2004, 41 (4): 391-427. 10.1080/10408360490497474.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Griffitts JS, Aroian RV: Many roads to resistance: how invertebrates adapt to Bt toxins. BioEssays. 2005, 27 (6): 614-624. 10.1002/bies.20239.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Gygi SP, Rochon Y, Franza BR, Aebersold R: Correlation between protein and mRNA abundance in yeast. Mol Cell Biol. 1999, 19 (3): 1720-1730.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Feyereisen R: Insect CYP genes and P450 enzymes. Insect Molecular Biology and Biochemistry. Edited by: Gilbert LI. 2011, Academic Press , 236-316.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Girish V, Vijayalakshmi A: Affordable image analysis using NIH Image/ImageJ. Indian J Cancer. 2004, 41 (1): 47-PubMedGoogle Scholar
- WHO: Guidelines for Laboratory and Field Testing of Mosquito Larvicides. 2005, Geneva: World Health Organization, WHO/CDS/WHOPES/GCDPP/2005.13Google Scholar
- 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.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Pierleoni A, Martelli PL, Casadio R: PredGPI: a GPI-anchor predictor. BMC Bioinforma. 2008, 9:Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Terra WR, Ferreira C: Insect digestive enzymes—properties, compartmentalization and function. Comp Biochem Physiol B-Biochem Mol Biol. 1994, 109 (1): 1-62.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2007,http://www.R-project.org,Google Scholar