- Open Access
Insecticide resistance in Aedes aegypti populations from Ceará, Brazil
Parasites & Vectorsvolume 4, Article number: 5 (2011)
Organophosphates and pyrethroids are used widely in Brazil to control Aedes aegypti, the main vector of dengue viruses, under the auspices of the National Programme for Dengue Control. Resistance to these insecticides is widespread throughout Brazil. In Ceará the vector is present in 98% of districts and resistance to temephos has been reported previously. Here we measure resistance to temephos and the pyrethroid cypermethrin in three populations from Ceará and use biochemical and molecular assays to characterise resistance mechanisms.
Resistance to temephos varied widely across the three studied populations, with resistance ratios (RR95) of 7.2, 30 and 192.7 in Juazeiro do Norte, Barbalha and Crato respectively. The high levels of resistance detected in Barbalha and Crato (RR95 ≥ 30) imply a reduction of temephos efficacy, and indeed in simulated field tests reduced effectiveness was observed for the Barbalha population. Two populations (Crato and Barbalha) were also resistant to cypermethrin, whilst Juazeiro do Norte showed only an altered susceptibility. The Ile1011Met kdr mutation was detected in all three populations and Val1016Ile in Crato and Juazeiro do Norte. 1011Met was significantly associated with resistance to cypermethrin in the Crato population. Biochemical tests showed that only the activity of esterases and GSTs, among the tested detoxification enzymes, was altered in these populations when compared with the Rockefeller strain.
Our results demonstrate that two A. aegypti populations from Ceará are under strong selection pressure by temephos, compromising the field effectiveness of this organophosphate. Our results also provide evidence that the process of reducing resistance to this larvicide in the field is difficult and slow and may require more than seven years for reversal. In addition, we show resistance to cypermethrin in two of the three populations studied, and for the first time the presence of the allele 1016Ile in mosquito populations from northeastern Brazil. A significant association between 1011M et and resistance was observed in one of the populations. Target-site mechanisms seem not to be implicated in temephos resistance, reinforcing the idea that for the studied populations, detoxification enzymes most likely play a major role in the resistance to this insecticide.
In Brazil, the wide distribution of Aedes aegypti - vector of the three circulating dengue virus serotypes (DENV-1, -2 and -3) - and the lack of a specific treatment or vaccine for dengue, are considered serious public health issues. This situation is exacerbated by the recent re-introduction of serotype 4 (DENV-4) in the State of Roraima, in the North of Brazil .
Aedes aegypti is present in all states and in most cities of Brazil. In Ceará State, this mosquito has been detected in approximately 98% of cities . In the last 23 years dengue has been endemic in the state, with four major epidemic periods registered in 1987, 1994, 2001 and 2008. In the latter year, Ceará was the most affected Brazilian state with respect to dengue cases, with 44,244 registered cases - an incidence of 530.77 cases/100,000 inhabitants. In 2009, although case numbers were lower than in the previous year, the mortality rate for dengue hemorrhagic fever was 34.6% and among severe dengue cases was 45% , significantly higher than the average mortality elsewhere in Brazil .
For the past 13 years, insecticides of the organophosphate (OP) and pyrethroid (PY) classes have been utilized intensively as the main strategy to control A. aegypti, run by the National Programme for Dengue Control (Programa Nacional de Controle da Dengue - PNCD) [3, 4]. Monitoring of the susceptibility status of A. aegypti populations to these insecticides has been undertaken through the National Network for Monitoring Resistance to Insecticides in Aedes aegypti (Rede Nacional de Monitoramento da Resistência de Aedes aegypti a Inseticidas - MoReNAa), since 1999 . The first signs of incipient resistance to temephos were registered in 1999 in mosquito populations from the State of São Paulo , and in other states since 2001 . Currently, it is known that temephos resistance is widespread in A. aegypti populations throughout the country [5, 7–12]. Consequently, the MoReNAa network has also been working on assaying the activity of detoxification enzymes in A. aegypti populations since 2001, aiming at identifying alterations that may be related to temephos resistance [12, 13].
In Ceará, resistance to temephos was detected in surveys conducted from 2000-2002 in A. aegypti populations from Fortaleza, Caucaia  and Juazeiro do Norte . In an effort to manage resistance in the field, temephos was replaced by the biolarvicide Bacillus thuringiensis israelensis (Bti) in these cities, and others, as a recommendation of the MoReNAa network . However, in most of the municipalities of Ceará (164) temephos continued to be utilized until the beginning of 2010, including Crato and Barbalha, where resistance to this insecticide had been first detected in A. aegypti in 2003 . Although the dengue control program at city level prioritizes the use of insecticides, a continuous evaluation of the effectiveness of temephos-based products in the field is not performed in most of them , despite PNCD recommendations to do so.
Data released by the MoReNAa network obtained from sentinel cities, revealed that A. aegypti populations were also resistant to malathion (another OP used as an adulticide) in 1999/2000, when its use was interrupted and replaced by pyrethroids as a resistance management action [5, 10]. However, this change led to cypermethrin/deltamethrin resistance three years later . After detecting resistance to all insecticides used by PNCD, this network carried out studies to characterize resistance mechanisms in Brazilian field populations.
Generally, two main mechanisms are commonly associated with resistance to chemical insecticides: metabolic, via increased activity of detoxification enzymes such as esterases, mixed function oxidases (cytochrome P450s) and glutathione S-transferases (GSTs); and structural modifications in insecticide binding sites, such as acetylcholinesterase and the voltage-gated sodium (NaV) channel. Recent studies suggest that in addition to metabolic resistance, mutations in the sodium channel (the target of pyrethroids as well as DDT) may be playing a role in the resistance to cypermethrin and deltamethrin in mosquito populations from Brazil .
PY compounds act on the insect nervous system, targeting the NaV channel. This channel is composed of four domains (I-IV) and each domain comprises six transmembrane helices (S1-S6) [18, 19]. Pesticides such as PY and DDT retard the activation and inactivation potential of NaV channels, triggering a series of repetitive discharges in motor and sensory axons, and resulting in paralysis ("knock-down") and death . However, several insect species, including Anopheles gambiae, Culex pipiens, Culex quinquefasciatus and A. aegypti, may present a resistance phenotype to chemicals which target the NaV, commonly called knock-down resistance [17, 21–27]. The genetic basis of knock-down resistance was first elucidated in the housefly Musca domestica. The knock-down resistance trait (named kdr) and another kdr-related trait (super-kdr), which confers greatly elevated resistance in combination with kdr, were mapped to chromosome 3 . Both traits have been associated with a lower electrophysiological sensitivity of elements from the nervous system and a reduced function of the NaV channel. Many studies have focused on finding mutations in NaV channel sequences from knock-down resistant populations. Characterization of sequences from A. gambiae and C. quinquefasciatus pyrethroid-resistant strains showed that the most common mutation is a leucine to phenylalanine substitution in the S6 hydrophobic segment of domain II [21, 22], although a leucine to serine mutation has also been reported at the same 1014 site [24, 29]. Brengues et al.  have failed to find the same mutation in A. aegypti, however, other studies have found different mutations correlated with kdr in this vector, such as: Gly923Val, Leu982Trp, Ile1011Val, Ile1011Met, Val1016Ile and Val1016Gly[17, 25, 30]. A few reports have shown that kdr genotyping is a good predictor of susceptibility to pyrethroid and DDT, and, at the moment, it is considered the best tool for predicting the efficacy of these compounds in the field .
The aim of this study was to measure resistance levels to two insecticides and to characterize resistance mechanisms at both the molecular and biochemical levels in A. aegypti populations from Ceará State.
Characterization of study area
mosquito populations - this study was undertaken in three municipalities, Juazeiro do Norte, Crato and Barbalha, located in the south region of Ceará State, Brazil, (Figure 1), 10 Km equidistant from each other. The local climate is tropical semi-arid, with average temperatures of 24-26°C and rainfall of 925 to 1,156 mm per year .
Juazeiro do Norte (7° 12' 47" S, 39° 18' 55" W) is 248.55 Km2 in area, and has a population of 249,831 with 50,000 properties . This city is the main commercial centre of the south region and is also one of the main religious routes of Brazil, attracting thousands of visitors each year from throughout the country, but particularly from the Northeast region. Crato (7° 14' 03" S, 39° 24' 34" W) has an area of 1,009.20 km2, 116,759 residents and 20,185 properties. This city is on the slope of the Araripe plateau. Barbalha (7° 18' 40" S, 39° 18' 15" W) has an area of 479.184 km², 52,496 residents and 7,032 properties.
Mosquito samples were collected in the field in November and December 2009 with the use of 100 oviposition traps (ovitraps), adapted from , distributed homogeneously in each area, based on the recommendations of the MoReNAa network [12, 35]. More than 50% of these ovitraps were positive for Aedes spp eggs. Field derived colonies from each selected population were established with approximately 2,000 mosquitoes (males and females) from these eggs. Populations were kept under controlled conditions in the insectary of the Federal University of Ceará (Cariri Campus), in order to obtain filial generation (F1), which was utilized in the biological, biochemical and molecular experiments described below. The Rockefeller strain, the standard strain for insecticide susceptibility tests, and used in all experiments as a susceptible reference strain, was obtained from a sub-strain kindly provided by Núcleo de Pesquisa/Sucen/Marília-SP (member of the MoReNAa/MS network, and initially provided by CDC, Puerto Rico).
temephos, an organophosphate larvicide, was used under two formulations: a standard powder (Pestanal, 97.5%, Sigma lot 6333X), used to make a stock solution in ethanol for in vivo laboratory bioassays, and as sand granules with 1% active ingredient (Fersol, lot 197), utilized in the simulated field trials. The pyrethroid adulticide cypermethrin, as a standard powder (98%, Sigma lot 127K1099), was dissolved in acetone for the impregnated bottles assay.
In vivo bioassays with larvae
Dose-response bioassays were undertaken according to the methodology proposed by the World Health Organization to evaluate larval susceptibility to temephos . In these experiments, third-instar larvae (L3) were exposed to 10 different concentrations of the larvicide determined according to previous exploratory bioassays. For each concentration, and for the control, three replicates of 20 larvae were tested. Larval mortality was checked 24 h after exposure. All tests were repeated at least three times. Mortality data were utilized to calculate the lethal concentrations for 95% of exposed individuals (LC95) through linear regression log probit analysis, using SPSS software (version 8.0 for Windows). These values were compared to those obtained for the Rockefeller strain in order to estimate the resistance ratio (RR) for each population. The following criteria proposed by Mazzari and Georghiou  was adopted to classify the resistance level of populations: low (RR < 5), moderate (RR > 5 < 10) or high (RR >10).
Simulated field trial
the aim of this trial was to verify the effectiveness of temephos in its sand granules formulation to control larvae in the field presenting distinct resistance levels, using the concentration recommended by the PNCD (1 ppm). The trials were performed in the external area of the Department of Entomology of the School of Medicine in Cariri/Barbalha, based on recommendations by the MoReNAa network . Plastic reservoirs filled with 20 L of tap water, treated with 2 g of temephos sand granules and colonized individually with 20 L3 larvae, were placed in a covered area protected from the sun and rain. One third of the water in the containers was renewed three times per week. For the Rockefeller and each field population, one control and three treated replicates were included. Mortality was verified 24 h after insecticide exposure and live larvae were removed before the weekly re-colonization of containers. The evaluation was continued for 60 days, a period of time similar to the interval between insecticide applications in the field performed by the national program. In these tests, only populations from Barbalha and Juazeiro do Norte were used, since samples from Crato were scarce.
In vivo bioassays with adults
For each population, three tests were performed using bottles impregnated with insecticides, according to the methodology described by Brogdon and McAllister , and Melo-Santos . For these, around 130 unfed one day-old females were utilized. Mosquitoes were exposed to a diagnostic dose of 8 μg/bottle of cypermethrin. Mortality was checked every 15 minutes over a period of 2 h, based on a lethal time to kill 100% of the population (LT100) of 15 mins at this concentration for the susceptible Rockefeller strain. All live and dead (those that did not move throughout the observation period) mosquitoes were transferred to another container free of insecticides and observed 24 h later, to determine the final mortality rate and the resistance status of populations. According to , populations were classified as susceptible when mortality was > 98%, in need of further evaluation when > 80% < 98%, or resistant when < 80%. Differences between resistant and susceptible individuals were evaluated among populations through Fisher's exact test. Following bioassays, females were stored according to phenotype, at -80°C for analysis of kdr mutations, which could possibly be associated with pyrethroid resistance.
Individual mosquitoes were separated according to the resistance phenotype, established at the final mortality (24 h) as resistant (survivors) and susceptible (dead mosquitoes) and homogenized in DNAzol® (Invitrogen, Carlsbad, CA, USA), for DNA extraction according to the manufacturer's instructions. DNA samples were later processed in order to analyze Ile1011Met and Val1016Ile Aa NaV mutations. To detect the Aa NaVIle1011Met mutation, primers were designed based on the IIS6 region, as described by Saavedra-Rodriguez et al  and Martins et al [17, 25]. PCR reactions were performed using 0.2 units of Pfu polymerase (Promega), 0.5 μM dNTPs, 2 μM of each primer, 1.5 μM MgCl2 and 10 ng of DNA. PCR conditions were: 5 min at 94°C, followed by 30 cycles of 1 min at 94°C, 1 min at 58°C and 2 min at 72°C, with a final extension at 72°C for 10 min. PCR products were purified using GFX™ PCR DNA and Gel Band Purification Kit (Amersham Pharmacia Biotech) and sequenced (ABI 3000, Applied Biosystems). Sequences were manually edited. Detection of Aa NaVVal1016Ile mutation was carried out using the following allele-specific primers, designed by Saavedra-Rodriguez et al: reverse primer (5'-GGATGAACCGAAATTGGACAAAAGC-3') and allele-specific forward primers (Val 5'-GCGGGCAGGGCGGCGGGGGCGGGGCCACAAATTGTTTCCCACCCGCACCGG-3'and Ile 5'-GCGGGCACAATTGTTTCCCACCCGCACTGA-3). PCR reactions and conditions were identical to the ones mentioned above. PCR products were run in 5% high-resolution agarose (Sigma-Aldrich) gels to determine genotypes (Val - 78 bp and Ile - 98 bp). To test for association between kdr mutations and resistance we used Haploview v4 .
approximately 100 females, randomly sampled from each field populations, were assayed by biochemical tests in order to verify the activity of the following insecticide detoxification enzymes: esterases (evaluated through the use of α and β naphthyl for α-Est and β-Est, respectively, and p-nitrophenyl PNPA), mixed function oxidases (MFO), glutathione S-transferases (GSTs), and acetylcholinesterase (AChE; the target of temephos). Microplate assays were performed on one-day-old adult females kept at -80°C . The Rockefeller strain was used as the susceptible reference. The proportion of individuals from each population presenting enzyme activity higher than that observed for the 99 percentile of Rockefeller strain was used for status classification as follows: unaltered <15%, altered >15% <50%, and highly altered >50% (for more details see [12, 41]).
Status of susceptibility to temephos
all tested populations presented some resistance level to temephos, with individuals from Crato and Barbalha presenting the highest levels and Juazeiro do Norte a moderate level (LC95 = 3.279 mg/ml, 0.510 mg/ml and 0.123 mg/ml respectively. c.f. 0.017 mg/ml for Rockefeller) (Table 1). Temephos has been used extensively over the period 2002-2009 in both Crato and Barbalha and in these municipalities resistance ratios have increased over this period with, for example, an increase in RR to temephos from 9.0 to 192.7 in Crato (21 fold increase) and from 7.5 to 30 in Barbalha (4 fold increase). In Juazeiro do Norte where temephos was substituted by Bti for larval control the RR to temephos decreased by 30% over this period (10.4 to 7.2) (Table 2).
Effectiveness of temephos in the field
in the simulated field trial, temephos killed 100% of the larvae in the two tested populations up to the 4th week after treatment. The larvicidal effect on the Juazeiro do Norte population varied between 90% and 98% from the 4th to the 8th week, while the mortality of larvae from the Barbalha population decreased to 66.6% in the 5th week, indicating a lower mortality rate than that recommended for validating the residual effect of temephos, according to the MoReNAa network. As expected, the insecticide killed all larvae of the control group (Rockefeller) throughout the experiment.
Status of susceptibility to cypermethrin
two of the three tested populations, Crato and Barbalha, were resistant to cypermethrin, with both populations presenting average mortality rates lower than 80% after 15 min and 24 h (Table 3). Crato was considered the most resistant population to this insecticide, since 51% of the exposed mosquitoes recovered 24 h after being removed to an insecticide-free recovery bottle. The mortality rate of individuals from Juazeiro do Norte indicated reduced susceptibility to this pyrethroid (Table 3). The differences in phenotypic distribution among populations were significant (p = 1.054 × 10-11).
Sequencing results showed that the Ile1011Met kdr mutation of the Nav channel was found in all three populations. The 1016Ile allele was detected in Crato and Juazeiro do Norte (in heterozygotes only) but not in Barbalha. Significant differences were observed in genotype frequencies between populations: for kdr1011; Crato vs. Barbalha (p = 0.0218), and for kdr1016; Crato vs. Juazeiro do Norte (p = 3.81 × 10-5) and Crato vs. Barbalha (7.0844 × 10-7). Phenotype frequencies of each population classified by genotype are shown in Figure 2 and Tables 4 and 5.
Barbalha and Juazeiro do Norte populations were found to be in Hardy-Weinberg equilibrium. However, the Crato population showed a significant deviation from Hardy-Weinberg expectation, due to a heterozygote deficit for Ile1011Met (p = 0.0197), and heterozygote excess for Val1016Ile (p = 0.015). The frequency of the 1011Met allele was higher in resistant individuals from Crato, and significantly associated with resistance to cypermethrin (p = 0.0242). In Barbalha and Juazeiro do Norte, the resistance did not show any association with this allele (Table 4). Although the resistant allele 1016Ile was observed in both resistant and susceptible A. aegypti from Crato, no association was found between this allele and the resistance status. For kdr 1016 in samples from Barbalha all individuals were homozygous for the 1011Val allele, thus, no association test was calculated. Since all individuals from Juazeiro do Norte were all phenotyped as susceptible, no test was performed for either putative kdr mutation (Table 5).
these assays showed an alteration in the activity of GST and α-esterase in all populations, in particular in individuals from Barbalha (Table 6, Figure 3). Esterases that metabolize PNPA were also altered only in the Barbalha and Juazeiro do Norte populations. Activity of other enzymes was not significantly different between populations when compared to the Rockefeller, including acetylcholinesterase (Table 6).
Over the last 10 years, resistance to organophosphate and pyrethroid insecticides employed for control of the dengue vector Aedes aegypti has been detected in all Brazilian regions, including Rio de Janeiro [4, 12], Espírito Santo  and São Paulo [7, 9–12] representing the Southeast, in all states from the Northeast region, in Distrito Federal and Goiás representing the Midwest region, and in Pará and Amazonas representing the North . Here we now show that two populations of A. aegypti from Crato and Barbalha (Ceará State, Northeast Brazil) are resistant to both temephos and cypermethrin. The temephos resistance level of the A. aegypti population from Crato (RR = 192.7) was the second highest ever recorded in Brazil, surpassed only by Araripina, Pernambuco, Northeast Brazil (RR = 240) .
Over recent years resistance to temephos has risen markedly in Crato (2003 RR = 9.0; 2009 RR = 192.7) and Barbalha (2003 RR = 7.5; 2009 RR = 30). In this study the resistance level of a neighbouring population (Juazeiro do Norte) was much lower (RR = 7.2), which is only marginally below the 2003 level (RR = 10.2). Here, temephos has not been used for Aedes control for at least seven years, having been substituted by Bti. The increased resistance levels to temephos observed for populations from Crato and Barbalha, compared to Juazeiro do Norte where increase in resistance has been halted, indicate that the resistance management strategy utilized in the field in Juazeiro do Norte should have been extended to these neighboring cities. Our results also demonstrate that even though temephos has been substituted in Juazeiro do Norte, recovery of susceptibility has been slow. Melo-Santos et al. , in a study that simulated different field situations for reversion of temephos resistance in A. aegypti, suggested that besides the probable fitness cost associated with resistance, interruption of the selective pressure alone is insufficient to completely revert resistance and may prolong the state of intermediate or low resistance unless susceptible individuals are introduced. However, this species has a low dispersal rate in densely populated environments  suggesting gene flow between A. aegypti populations is low, despite their close proximity, as demonstrated previously by Ayres et al.  in their study of genetic variability of populations from Northeast Brazil.
The high resistance levels detected for Crato and Barbalha suggest that temephos will have low efficacy for Aedes control. Through simulated field trials we have confirmed the low effectiveness of temephos-based products, used at the field dosage, for the Barbalha population (RR = 30), corroborating Montella et al. , who showed that an RR > 10 can compromise insecticide effectiveness in control programs. Our results show clearly that the period for temephos larvicidal effectiveness (30 days) was much lower than the interval (60 days) recommended by PNCD between temephos application cycles in A. aegypti breeding sites. Consequently, there are serious faults in the vector control activities as currently undertaken in Barbalha and, probably in Crato (where RR = 192.7), although we did not undertake simulated field-trial evaluation of this population. Recent entomological indexes based on larval surveys conducted in these cities provide further evidence that the control program has not been successful in maintaining mosquito populations at low densities, since A. aegypti was detected in over 1% of houses . The continued use of temephos in such areas will certainly decrease the efficacy of the insecticide within the scope of PNCD.
Resistance to temephos can arise through alterations in the target site of the insecticide (the acetylcholinesterase gene ACE-1) or through elevated levels or differential efficacy of metabolic genes. We utilized biochemical assays to characterize the mechanisms involved in insecticide resistance. No evidence was found for target-site resistance (AChE test). However, increased activity levels of GSTs, α-esterases and PNPA-esterases were detected indicating that members of these enzyme classes may play roles in detoxification and deserve further attention. Interestingly, A. aegypti from Crato, which presented the highest RR to temephos, has the lowest activity profile for these enzyme classes, indicating that the association between the increase of activity and level of temephos resistance is not simple. Thus, metabolic mechanisms do not seem to explain fully the elevated resistance levels to temephos, unless very few genes are involved whose activities could not be detected by our broad biochemical assays or unless other mechanisms, such as reduced insecticide penetration, are involved in temephos resistance. Future studies with synergists are warranted in order to provide additional information on metabolic mediated resistance mechanisms.
In Brazil, pyrethroids were first utilized by PNCD in 2000 [3, 17], when resistance to temephos was spreading throughout the country . Surveillance of cypermethrin susceptibility in A. aegypti populations from Brazil has been performed since 2001 when da Cunha et al.  demonstrated that populations from seven out of sixteen Brazilian cities presented an incipient alteration to this insecticide. Resistance increased over the following two years with 80% of the analyzed populations classified as resistant in 2003 . Our results demonstrate that A. aegypti populations from Juazeiro do Norte present only a slight alteration in cypermethrin susceptibility, whilst for the Crato and Barbalha populations the mortality rates to PY is among the lowest registered for Brazilian populations [10, 16]. Cypermethrin, alpha-cypermethrin or deltamethrin are sprayed monthly for mosquito control in strategic points (tire repair centres, cemeteries and junkyards). Additionally, ultra-low volume (ULV) spraying is recommended by the PNCD to block viral transmission as a complementary strategy when a high number of dengue cases are reported in specific areas  or in places with high mosquito infestation levels (>5%) in order to prevent outbreaks. In Juazeiro do Norte, the registered number of dengue cases was lower than Crato and Barbalha in the last year, implying that this population should have been subject to a lower insecticide pressure. In addition, since temephos is no longer used in this location, potential pyrethroid cross-resistance mechanisms would be minimized. As seen in other studies, resistance to PY seems to be rapidly evolving in populations from Ceará, where this insecticide started being utilized less than 10 years ago [10, 16, 45]. Additional factors that may have contributed to this rapid rise in resistance are: 1) the mechanisms involved in temephos resistance are also involved in PY resistance [46, 47]; 2) a concomitant use of PY insecticides in other control programs, such as those for Leishmaniasis, malaria and Chagas disease; 3) large-scale PY usage to control urban pests (e.g. cockroaches, mites, ants and scorpions), performed by private companies acting independently of, and unsupervised by, the Secretaries of Health and Environment; and 4) domestic use of spray insecticides.
In recent years, various studies have shown an association between mutations in the voltage-gated sodium (Nav) channel and resistance to pyrethroids in various mosquito species [21, 24, 30, 48–51]. In A. aegypti, a number of Nav mutations in segment 6 of domain II have been associated with kdr-like resistance, such as Gly923Val, Leu982Trp, Ile1011Met, Ile1011Val, Val1016Ile and Val1016Gly[17, 25, 30, 52], Asp1794Tyr between segments 5 and 6 of domain II , and, more recently, Phe1534Cys in segment 6 of domain III . Brengues et al. detected Ile1011Met within a population of A. aegypti from Belém (North region) resistant to pyrethroids and DDT. In electrophysiological studies, this mutation was associated with significantly different levels of nerve sensitivity to both permethrin and lambda cyhalothrin . However, the first countrywide screening for this mutation and its correlation with resistance did not occur until 2009 [17, 23]. Martins et al.  observed that in A.aegypti from Natal (Northeast Brazil) the frequency of the resistant allele (1011Met = 0.529) was higher in individuals resistant to cypermethrin than in susceptible individuals (0.425), although this difference was not significant. In the present study, we observed a significant association between the 1011Met frequency and resistance to cypermethrin in the population from Crato (p = 0.0242). Martins et al.  showed that the 1016Ile allele was found almost exclusively in the Midwest region of Brazil, being absent in the Southeast and Northeast, although no test for association with resistance was performed. In the present study, the 1016Ile allele was detected for the first time in A. aegypti populations from the Northeast of Brazil (Crato and Juazeiro do Norte). We did not detect any association between the 1016Ile variant and resistance.
The detection of significant association with the resistance phenotype for a kdr allele (1011Met) shown previously to have a functional effect on nerve response following insecticide exposure , and hence with potential importance for the knockdown phenotype, has important implications for control programmes since the continuous use of PY adulticides may promote a drastic increase in frequency in local mosquito populations potentially leading to fixation, as demonstrated by García et al.  in a study on A. aegypti from Mexico. Similarly, Lynd et al. , observed a rapid increase in 1014Phe frequency in An. gambiae from southern Ghana over a 5-year period. In addition to mutations in the target site of PY, resistance may also have a metabolic basis. Our study detected no significant alterations for MFOs (cytochrome P450s) which have been demonstrated previously to be capable of metabolizing pyrethroids in vitro. However, microarray studies of either field samples or laboratory colonies phenotyped for insecticide resistance have previously implicated particular cytochrome P450s in the resistance phenotype [57–59]. Alterations in activity of one, or a few, P450s may not be detectable by the broad MFO assay. However, we did detect significant increase in GST activity and previous studies have shown an association between an increase of GSTs activity and resistance to PY, although some of the populations studied, including those studied here, have also been exposed to OP [12, 17, 49, 60], complicating interpretation.
The status of temephos and cypermethrin resistance in these A. aegypti populations reinforces the need for a constant surveillance of mosquito susceptibility against insecticides used in control programs, as well as their effectiveness in the field. This must begin before, or as soon as the insecticide starts to be utilized, such that initial resistance levels are determined, so facilitating resistance management. In general, resistance to chemical insecticides is a multi-factorial trait that may be affected by environmental (availability and types of breeding sites), operational (insecticide application frequency and amount, period of exposure) and genetic (metabolic genes and target-site alteration) factors. Additionally, A. aegypti populations from Brazil possess different genetic backgrounds, especially those from the Northeast region , a fact that could be behind the diverse resistance mechanisms to the same insecticide developed in different populations. This reinforces the suggestion of Rawlins et al. , that control strategies must be adjusted to each place according to local peculiarities. A broad program such as PNCD should take this into consideration.
Fundação Nacional de Saúde: Nota técnica n° 110, de 6 de agosto de 2010. 2010, Brasília: Ministério da Saúde
Secretaria de Saúde do Ceará: Informe Semanal Dengue - Semana Epidemiológica 01 a 21. 2010, Fortaleza
Braga IA, Valle D: Aedes aegypti: insecticides, mechanisms of action and resistance. Epidemiologia e Serviços de Saúde. 2007, 16: 295-302.
Braga IA, Lima JB, Soares SS, Valle D: Aedes aegypti resistance to temephos during 2001 in several municipalities in the states of Rio de Janeiro, Sergipe, and Alagoas, Brazil. Mem Inst Oswaldo Cruz. 2004, 99: 199-203. 10.1590/S0074-02762004000200015.
Fundação Nacional de Saúde: Reunião Técnica Para Discutir Status de Resistência de Aedes aegypti e Definir Estratégias a Serem Implantadas Para Monitoramento da Resistência no Brasil. 1999, Brasília: Ministério da Saúde
Macoris MLG, Andrighetti MT, Takaku L, Glasser CM, Garbeloto VC, Cirino VC: Changes in susceptibility of Aedes aegypti to organophosphates in municipalities in the state of São Paulo, Brazil. Rev Saúde Pública. 1999, 33: 521-522. 10.1590/S0034-89101999000500013.
Beserra EB, Fernandes CR, Queiroga MF, de Castro FP: Resistance of Aedes aegypti (L.) (Diptera: Culicidae) populations to organophosphate temephos in the Paraíba State, Brazil. Neotrop Entomol. 2007, 36: 303-307. 10.1590/S1519-566X2007000200019.
Carvalho MS, Caldas ED, Degallier N, Vilarinhos PT, Souza LC, Yoshizawa MA, Knox MB, Oliveira C: Susceptibility of Aedes aegypti larvae to the insecticide temephos in the Federal District, Brazil. Rev Saude Publica. 2004, 38: 623-629. 10.1590/S0034-89102004000500002.
Lima EP, Filho AMO, Lima JWO, Junior ANR, Cavalcanti LPG, Pontes RJ: Aedes aegypti resistance to temefos in counties of Ceará State. Rev Soc Bras Med Trop. 2006, 39: 259-263. 10.1590/S0037-86822006000300006.
Macoris ML, Andrighetti MT, Otrera VC, Carvalho LR, Caldas Junior AL, Brogdon WG: Association of insecticide use and alteration on Aedes aegypti susceptibility status. Mem Inst Oswaldo Cruz. 2007, 102: 895-900. 10.1590/S0074-02762007000800001.
Melo-Santos MAV, Varjal-Melo JJ, Araújo AP, Gomes TCS, Paiva MH, Regis LN, Furtado AF, Magalhaes T, Macoris ML, Andrighetti MT, Ayres CFJ: Resistance to the organophosphate temephos: mechanisms, evolution and reversion in an Aedes aegypti laboratory strain from Brazil. Acta Trop. 2010, 113: 180-189. 10.1016/j.actatropica.2009.10.015.
Montella IR, Martins AJ, Viana-Medeiros PF, Lima JB, Braga IA, Valle D: Insecticide resistance mechanisms of Brazilian Aedes aegypti populations from 2001 to 2004. Am J Trop Med Hyg. 2007, 77: 467-477.
Macoris ML, Andrighetti MT, Takaku L, Glasser CM, Garbeloto VC, Bracco JE: Resistance of Aedes aegypti from the state of Sao Paulo, Brazil, to organophosphates insecticides. Mem Inst Oswaldo Cruz. 2003, 98: 703-708. 10.1590/S0074-02762003000500020.
Filho AMO, Melo MTV, Santos CE, Griffo HG, Lustosa ELB, Ribeiro ZM: Resistência a Temefós em linhagens de Aedes aegypti provenientes do Ceará. In: Resumos do XXXVII Congresso da Sociedade Brasileira de Medicina Tropical, Salvador. 2001, 291-
Secretaria de Vigilância em Saúde/Ministério da Saúde: Reunião técnica para discutir e avaliar os resultados do monitoramento de resistência do Aedes aegypti a inseticidas. 2004, Brasília-DF
da-Cunha MP, Lima JB, Brogdon WG, Moya GE, Valle D: Monitoring of resistance to the pyrethroid cypermethrin in Brazilian Aedes aegypti (Diptera: Culicidae) populations collected between 2001 and 2003. Mem Inst Oswaldo Cruz. 2005, 100: 441-444. 10.1590/S0074-02762005000400017.
Martins AJ, Lins RM, Linss JG, Peixoto AA, Valle D: Voltage-gated sodium channel polymorphism and metabolic resistance in pyrethroid-resistant Aedes aegypti from Brazil. Am J Trop Med Hyg. 2009, 81: 108-115.
Hemingway J, Hawkes NJ, McCarroll L, Ranson H: The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004, 34: 653-665. 10.1016/j.ibmb.2004.03.018.
Soderlund DM, Knipple DC: The molecular biology of knockdown resistance to pyrethroid insecticides. Insect Biochem Mol Biol. 2003, 33: 563-577. 10.1016/S0965-1748(03)00023-7.
Vais H, Williamson MS, Devonshire AL, Usherwood PN: The molecular interactions of pyrethroid insecticides with insect and mammalian sodium channels. Pest Manag Sci. 2001, 57: 877-888. 10.1002/ps.392.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, Guillet P, Pasteur N, Pauron D: Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998, 7: 179-184. 10.1046/j.1365-2583.1998.72062.x.
Martinez-Torres D, Foster SP, Field LM, Devonshire AL, Williamson MS: A sodium channel point mutation is associated with resistance to DDT and pyrethroid insecticides in the peach-potato aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae). Insect Mol Biol. 1999, 8: 339-346. 10.1046/j.1365-2583.1999.83121.x.
Martins AJ, Lima JB, Peixoto AA, Valle D: Frequency of Val1016Ile mutation in the voltage-gated sodium channel gene of Aedes aegypti Brazilian populations. Trop Med Int Health. 2009, 14: 1351-1355. 10.1111/j.1365-3156.2009.02378.x.
Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH: Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000, 9: 491-497. 10.1046/j.1365-2583.2000.00209.x.
Saavedra-Rodriguez K, Urdaneta-Marquez L, Rajatileka S, Moulton M, Flores AE, Fernandez-Salas I, Bisset J, Rodriguez M, McCall PJ, Donnelly MJ, Ranson H, Hemingway J, Black 4th WC: A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti. Insect Mol Biol. 2007, 16: 785-798. 10.1111/j.1365-2583.2007.00774.x.
Wondji CS, Priyanka De Silva WA, Hemingway J, Ranson H, Parakrama Karunaratne SH: Characterization of knockdown resistance in DDT- and pyrethroid-resistant Culex quinquefasciatus populations from Sri Lanka. Trop Med Int Health. 2008, 13: 548-555. 10.1111/j.1365-3156.2008.02033.x.
Xu Q, Wang H, Zhang L, Liu N: Kdr allelic variation in pyrethroid resistant mosquitoes, Culex quinquefasciatus (S.). Biochem Biophys Res Commun. 2006, 345: 774-780. 10.1016/j.bbrc.2006.04.155.
Busvine JR: Mechanism of resistance to insecticide in houseflies. Nature. 1951, 168: 193-195. 10.1038/168193a0.
Wondji CS, Hemingway J, Ranson H: Identification and analysis of single nucleotide polymorphisms (SNPs) in the mosquito Anopheles funestus, malaria vector. BMC Genomics. 2007, 8: 5-10.1186/1471-2164-8-5.
Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S, Guillet P, Manguin S, Morgan JC, Hemingway J: Pyrethroid and DDT cross-resistance in Aedes aegypti is correlated with novel mutations in the voltage-gated sodium channel gene. Med Vet Entomol. 2003, 17: 87-94. 10.1046/j.1365-2915.2003.00412.x.
Donnelly MJ, Corbel V, Weetman D, Wilding CS, Williamson MS, Black 4th WCt: Does kdr genotype predict insecticide-resistance phenotype in mosquitoes?. Trends Parasitol. 2009, 25: 213-219. 10.1016/j.pt.2009.02.007.
Instituto de Pesquisa e Estratégia Econômica do Ceará: Perfil básico municipal de Juazeiro do Norte. 2007
Instituto Brasileiro de Geografia e Estatística: Pesquisa de Informações Básicas Municipais. 2010, Juazeiro do Norte
Fay RW, Eliason DA: A preferred oviposition site as a surveillance method for Aedes aegypti. Mosq News. 1966, 26: 531-535.
Lima JB, Da-Cunha MP, Da Silva RC, Galardo AK, Soares SS, Braga IA, Ramos RP, Valle D: Resistance of Aedes aegypti to organophosphates in several municipalities in the State of Rio de Janeiro and Espirito Santo, Brazil. Am J Trop Med Hyg. 2003, 68: 329-333.
World Health Organization: Instructions for determining the susceptibility or resistance of mosquito larvae to insecticides. 1981, Geneve
Mazzarri MB, Georghiou GP: Characterization of resistance to organophosphate, carbamate, and pyrethroid insecticides in field populations of Aedes aegypti from Venezuela. J Am Mosq Control Assoc. 1995, 11: 315-322.
Brogdon WG, McAllister JC: Simplification of adult mosquito bioassays through use of time-mortality determinations in glass bottles. J Am Mosq Control Assoc. 1998, 14: 159-164.
World Health Organization: Test procedures for insecticide resistance montoring in malaria vectors, bio-efficacy and persistence of insecticides on treated surfaces. 1998, Geneva, Switzerland
Barrett JC FB, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
Ministério da Saúde: Metodologia para Quantificacão de Atividade de Enzimas Relacionadas com a Resistência a Inseticidas em Aedes aegypti. 2006, Brasília-DF
Harrington LC, Scott TW, Lerdthusnee K, Coleman RC, Costero A, Clark GG, Jones JJ, Kitthawee S, Kittayapong P, Sithiprasasna R, Edman JD: Dispersal of the dengue vector Aedes aegypti within and between rural communities. Am J Trop Med Hyg. 2005, 72: 209-220.
Ayres CFJ, Melo-Santos MAV, Prota JR, Solé-Cava AM, Regis LN, Furtado AF: Genetic structure of natural populations of Aedes aegypti at the micro- and macrogeographic levels in Brazil. J Am Mosq Control Assoc. 2004, 20: 350-356.
Fundação Nacional de Saúde: Dengue - Instruções para pessoal de combate ao vetor. 2001, Brasília: Ministério da Saúde
Luna JE, Martins MF, Anjos AF, Kuwabara EF, Navarro-Silva MA: Susceptibility of Aedes aegypti to temephos and cypermethrin insecticides, Brazil. Rev Saude Publica. 2004, 38: 842-843. 10.1590/S0034-89102004000600013.
Rodriguez MM, Bisset JA, Mila LH, Calvo E, Diaz C, Soca LA: Levels of insecticide resistance and its mechanisms in a strain of Aedes aegypti of Santiago de Cuba. Rev Cubana Med Trop. 1999, 51: 83-88.
Wirth MC, Georghiou GP: Selection and characterization of temephos resistance in a population of Aedes aegypti from Tortola, British Virgin Islands. J Am Mosq Control Assoc. 1999, 15: 315-320.
Diabate A, Brengues C, Baldet T, Dabire KR, Hougard JM, Akogbeto M, Kengne P, Simard F, Guillet P, Hemingway J, Chandre F: The spread of the Leu-Phe kdr mutation through Anopheles gambiae complex in Burkina Faso: genetic introgression and de novo phenomena. Trop Med Int Health. 2004, 9: 1267-1273. 10.1111/j.1365-3156.2004.01336.x.
Etang J, Vicente JL, Nwane P, Chouaibou M, Morlais I, Do Rosario VE, Simard F, Awono-Ambene P, Toto JC, Pinto J: Polymorphism of intron-1 in the voltage-gated sodium channel gene of Anopheles gambiae s.s. populations from Cameroon with emphasis on insecticide knockdown resistance mutations. Mol Ecol. 2009, 18: 3076-3086. 10.1111/j.1365-294X.2009.04256.x.
Lynd A, Ranson H, McCall PJ, Randle NP, Black 4th WC, Walker ED, Donnelly MJ: A simplified high-throughput method for pyrethroid knock-down resistance (kdr) detection in Anopheles gambiae. Malar J. 2005, 4: 16-10.1186/1475-2875-4-16.
Morgan JC, Irving H, Okedi LM, Steven A, Wondji CS: Pyrethroid Resistance in an Anopheles funestus Population from Uganda. PLoS One. 2010, 5: e11872-10.1371/journal.pone.0011872.
Rajatileka S, Black 4th WC, Saavedra-Rodriguez K, Trongtokit Y, Apiwathnasorn C, McCall PJ, Ranson H: Development and application of a simple colorimetric assay reveals widespread distribution of sodium channel mutations in Thai populations of Aedes aegypti. Acta Trop. 2008, 108: 54-57. 10.1016/j.actatropica.2008.08.004.
Chang C, Shen WK, Wang TT, Lin YH, Hsu EL, Dai SM: A novel amino acid substitution in a voltage-gated sodium channel is associated with knockdown resistance to permethrin in Aedes aegypti. Insect Biochem Mol Biol. 2009, 39: 272-278. 10.1016/j.ibmb.2009.01.001.
Yanola J, Somboon P, Walton C, Nachaiwieng W, Prapanthadara L: A novel F1552/C1552 point mutation in the Aedes aegypti voltage-gated sodium channel gene associated with permethrin resistance. Pesticide Biochemistry and Physiology. 2010, 96: 127-131. 10.1016/j.pestbp.2009.10.005.
Garcia GP, Flores AE, Fernandez-Salas I, Saavedra-Rodriguez K, Reyes-Solis G, Lozano-Fuentes S, Guillermo Bond J, Casas-Martinez M, Ramsey JM, Garcia-Rejon J, Domínguez-Galera M, Ranson H, Hemingway J, Eisen L, Black 4th WC: Recent rapid rise of a permethrin knock down resistance allele in Aedes aegypti in Mexico. PLoS Negl Trop Dis. 2009, 3: e531-10.1371/journal.pntd.0000531.
Lynd A, Weetman D, Barbosa S, Egyir Yawson A, Mitchell S, Pinto J, Hastings I, Donnelly MJ: Field, genetic, and modeling approaches show strong positive selection acting upon an insecticide resistance mutation in Anopheles gambiae s.s. Mol Biol Evol. 2010, 27: 1117-1125. 10.1093/molbev/msq002.
Muller P, Warr E, Stevenson BJ, Pignatelli PM, Morgan JC, Steven A, Yawson AE, Mitchell SN, Ranson H, Hemingway J, Paine MJ, Donnelly MJ: Field-caught permethrin-resistant Anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids. PLoS Genet. 2008, 4: 1000286-10.1371/journal.pgen.1000286.
Djouaka RF, Bakare AA, Coulibaly ON, Akogbeto MC, Ranson H, Hemingway J, Strode C: Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s.s. from Southern Benin and Nigeria. BMC Genomics. 2008, 9: 538-10.1186/1471-2164-9-538.
Muller P, Donnelly MJ, Ranson H: Transcription profiling of a recently colonised pyrethroid resistant Anopheles gambiae strain from Ghana. BMC Genomics. 2007, 8: 36-10.1186/1471-2164-8-36.
Rodriguez MM, Bisset JA, De Armas Y, Ramos F: Pyrethroid insecticide-resistant strain of Aedes aegypti from Cuba induced by deltamethrin selection. J Am Mosq Control Assoc. 2005, 21: 437-445. 10.2987/8756-971X(2006)21[437:PISOAA]2.0.CO;2.
Rawlins SC: Spatial distribution of insecticide resistance in Caribbean populations of Aedes aegypti and its significance. Rev Panam Salud Publica. 1998, 4: 243-251. 10.1590/S1020-49891998001000004.
We are indebted to all the health staff (health officers and operational vector-control staff) of the Secretaria de Saúde do estado do Ceará and the Núcleo de Entomologia Médica do Cariri, who participated in the mosquito collection in the municipalities of Ceará State. We also thank the Program for Technological development in Tools for Health- PDTIS-FIOCRUZ for providing sequencing assistance.
The work was funded by National Council for Scientific and Technological Development (CNPq, Grants ID: 479825/2009-2 and 550116/2010-9), CAPES/PROCAD/NF 23, Rede Nordeste de Biotecnologia (RENORBIO) and Fundação de Amparo à Ciência e Tecnologia do Estado de Pernambuco (FACEPE, Grant ID: APQ 1361-4.00/08).
The authors declare that they have no competing interests.
EPL and MHSP contributed equally to the study. EPL, AES and MOFG conceived the idea of the study; MAVMS designed and coordinated the bioassays, biochemical and semi field tests; UMS, LNO, CNB and CCPN have collected the samples and performed the field-simulated experiments; EPL, APA and EVGS carried out the bioassay and biochemical tests and its analysis; CFJA designed and coordinated the molecular study; EPL and MHSP carried out the molecular assays; CSW and MHSP performed the sequence alignment and statistical analysis; EPL, MHSP, MOFG, MAVMS, CSW and CFJA helped to draft the manuscript and all authors read and approved the final version. MAVMS is guarantor of the paper.