Skip to main content
Advertisement
  • Loading metrics

Pyrethroid resistance persists after ten years without usage against Aedes aegypti in governmental campaigns: Lessons from São Paulo State, Brazil

  • Maria de Lourdes Macoris ,

    Contributed equally to this work with: Maria de Lourdes Macoris, Ademir Jesus Martins

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Writing – original draft, Writing – review & editing

    Affiliation Laboratório de Entomologia Aplicada, Sucen, Marília, SP. Brazil

  • Ademir Jesus Martins ,

    Contributed equally to this work with: Maria de Lourdes Macoris, Ademir Jesus Martins

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing

    ademirjr@ioc.fiocruz.br (AJM); dvalle@ioc.fiocruz.br (DV)

    Affiliations Laboratório de Fisiologia e Controle de Artrópodes Vetores, Instituto Oswaldo Cruz, Fiocruz. Rio de Janeiro, RJ. Brazil, Instituto Nacional de Ciência e Tecnologia em Entomologia Molecular (INCT-EM), Rio de Janeiro, Brazil

  • Maria Teresa Macoris Andrighetti,

    Roles Data curation, Funding acquisition, Project administration, Resources, Writing – review & editing

    Affiliation Laboratório de Entomologia Aplicada, Sucen, Marília, SP. Brazil

  • José Bento Pereira Lima,

    Roles Methodology, Resources, Writing – review & editing

    Affiliation Laboratório de Fisiologia e Controle de Artrópodes Vetores, Instituto Oswaldo Cruz, Fiocruz. Rio de Janeiro, RJ. Brazil

  • Denise Valle

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    ademirjr@ioc.fiocruz.br (AJM); dvalle@ioc.fiocruz.br (DV)

    Affiliations Instituto Nacional de Ciência e Tecnologia em Entomologia Molecular (INCT-EM), Rio de Janeiro, Brazil, Laboratório de Biologia Molecular de Flavivírus, Instituto Oswaldo Cruz, Fiocruz. Rio de Janeiro, RJ. Brazil

Abstract

Background

Aedes aegypti, vector of dengue, chikungunya and Zika viruses, is found at high densities in tropical urban areas. The dissemination of this vector is partially the consequence of failures in current vector control methods, still mainly relying upon insecticides. In the State of São Paulo (SP), Brazil, public health managers employed pyrethroids against Ae. aegypti adults from 1989 to 2000, when a robust insecticide resistance monitoring system detected resistance to pyrethroids in several Ae. aegypti populations. However, pyrethroids are also the preferred compounds engaged in household applications due to their rapid knockdown effect, lower toxicity to mammals and less irritating smell.

Methodology/Principal findings

We evaluated pyrethroid resistance in Ae. aegypti populations over the course of a decade, from 2004 to 2015, after interruption of pyrethroid public applications in SP. Qualitative bioassays with papers impregnated with a deltamethrin diagnostic dose (DD) performed with insects from seven SP municipalities and evaluated yearly from 2006 to 2014, detected resistance in most of the cases. Quantitative bioassays were also carried out with four populations in 2011, suggesting a positive correlation between resistance level and survivorship in the DD bioassays. Biochemical tests conducted with seven insect populations in 2006 and 2015, detected increasing metabolic alterations of all major classes of detoxifying enzymes, mostly of mixed function oxidases. Genotyping of the voltage-gated sodium channel (AaNaV, the pyrethroid target-site) with a TaqMan real time PCR based technique was performed from 2004 to 2014 in all seven localities. The two kdr mutations, Val1016Ile and Phe1534Cys, known to be spread throughout Brazil, were always present with a severe decrease of the susceptible allele over time.

Conclusions/Significance

These results are discussed in the context of public and domestic insecticide use, the necessity of implementation of a strong integrated vector control strategy and the conceptual misunderstanding between 'vector control' and 'chemical control of vectors'.

Author summary

Insecticides will continue to exert an essential role in the reduction of vector density and personal protection for a long time, especially during epidemics. In Brazil, these chemicals are extensively incorporated in national campaigns against Ae. aegypti, favoring the rise and dispersion of resistant populations. The State of São Paulo (SP), the capital (also named São Paulo) being the 5th most populated city in the world, interrupted pyrethroid application in 2000 given the high resistance levels registered. Nevertheless, SP Ae. aegypti populations remain resistant to pyrethroids with an ever increasing frequency of kdr mutations, one of the main mechanisms responsible for resistance to these chemicals. Although no longer employed by public authorities, pyrethroids are readily available on the commercial market and are extensively used by households as well as private companies, which may partially account for the high resistance ratios still in effect. All in all, it is essential to invest in basic infrastructural sanitation and community educational methods to consciously avoid larval breading sites as well as encourage the correct use of insecticides. Without this, we are fated to forfeit the efficacy of all insecticide classes as well as undermine any other new methods of Aedes control.

Introduction

Latin America has been the epicenter not only of a severe increase in the number of dengue cases, but also of the more recently emerging chikungunya and Zika viruses [1, 2]. The recent cluster of microcephaly in newborns together with other neurological disorders reported in Brazil following a similar scenario in French Polynesia in 2014, has led the World Health Organization (WHO) to declare Zika as a Public Health Emergency of International Concern [3, 4]. The lack of specific drugs and effective vaccines against these infections reinforce the importance of controlling their main vector, the Aedes aegypti mosquito. Elimination or sealing of artificial containers and water accumulation in domestic and peridomestic surroundings is increasingly considered the most effective way to avoid proliferation of potential larval breeding sites. However, insecticides still play a relevant role in vector control, either as adulticides through space spraying or as larvicides applied in man-made water collections [5]. Unfortunately, the intense use of these insecticides has resulted in the selection of resistant mosquito populations to compounds of several classes [6], hampering Ae. aegypti control and in consequence, the pathogens transmitted by this vector.

The current goal of public health policy makers is to keep Ae. aegypti infestation levels below considered risk thresholds [7, 8]. However in the 1950's, Ae. aegypti erradication was already a challenge with the intense use of organochlorine DDT selecting resistant populations virtually at a worldwide scale [9]. Similarly, resistance to pyrethroids is now widespread throughout the globe [10], mainly due to a serious conceptual misunderstanding between 'vector control' and 'chemical control of vectors' [11]. Ideally insecticide use should be an activity complementary to the mechanical control, the latter based on the removal or elimination of potential breeding sites. Nevertheless, pyrethroids compared to other classes of insecticides are safer for man, bear high insecticidal potency and possess a rapid knock-down effect [12], which contribute to their intense use by both public managers and private initiative. Synthetic pyrethroid resistance dates from the 1990s for Ae. aegypti populations from US, Caribbean and Asia [9] and has been registered in Brazil since 2000 [1315]. More recent studies point to increasing resistance levels to pyrethroids in consequence of their use against dengue in endemic areas, especially in Asia, Caribbean and South America [1619].

Throughout the XX century, Brazil was considered free of Ae. aegypti twice, in 1958 and in 1973, but this mosquito was detected again in 1967 and later in 1976. However, it was only since 1985 that dengue outbreaks began uninterruptedly in the country [8,20]. In São Paulo (SP), the most urbanized Brazilian State, a novel Ae. aegypti infestation was detected in 1985, followed by the first local dengue case in its western region, in 1986 [21]. At the beginning of the 1990’s, dengue transmission occurred in the northern SP region quickly expanding to its northwest and coastal regions with high incidence around the seaside city of Santos. High incidence (number of cases/100,000 inhabitants) of dengue was registered in SP between 2006 (150.2) and 2007 (275.9) with transmission expanding toward the central west of the State, reaching 503.0 in 2010. The Brazilian Ministry of Health considers dengue incidence rates above 300 indicative of an epidemic scenario [22]. In 2014 the incidence of dengue was 515.2 in SP disseminating throughout the whole State with distinct incidence patterns [23]. Currently in Brazil besides dengue, chikungunya and Zika viruses are also autochthonously transmitted by Ae. aegypti [24]. In 2016, SP alone registered 155,972, 3,857 and 232 confirmed autochthonous cases of dengue, Zika and chikungunya, respectively [23].

The Superintendência de Controle de Endemias (Sucen) is the SP division of the Health Secretariat responsible for coordinating insect vector surveillance and control activities in the whole State. Space spraying has been adopted since 1985, applications generally restricted to the summer and autumn in conjunction with the seasonality of both mosquito infestation levels and dengue cases [25]. Several insecticides have been employed against adults in SP, in particular the carbamate propoxur (1986–1989), the organophosphate malathion (1985–1992) and the pyrethroid cypermethrin (1989 to 2000) [26]. In practice, SP State started and stopped pyrethroid application against Ae. aegypti adults long before the rest of the country. In 1996, Sucen initiated an insecticide resistance surveillance program of Ae. aegypti populations in SP [27]. This program monitored municipalities subjected to an intense use of insecticides conducted by public health managers due to high dengue incidence [28]. A Brazilian network of Ae. aegypti insecticide resistance monitoring covering mosquito populations throughout the country was also implemented, in 1999 [29]. In 2000, cypermethrin resistance detection in several Ae. aegypti populations effected the interruption of pyrethroids and the resumption of malathion in SP. However, the uncontrolled domestic and private use of pyrethroids prevailed by means of ordinary commercial spraying and fogging hired from private agencies.

DDT and pyrethroids induce a characteristic intoxication known as the knockdown effect in which the insect undergoes fast and repetitive muscle spasms followed by paralysis and eventually death [30]. One of the main resistance mechanisms selected against these classes of insecticides is a change in their molecular target, the voltage gated sodium channel (NaV). NaV mutations are related to knockdown resistance in several insects, therefore denominated kdr mutations, most of which are substitutions of Leu, in the 1014 codon, for Phe, Ser or Hys, as witnessed in Anopheles and Culex mosquitoes [31, 32]. Due to differences in the codon usage, kdr mutations at position 1014 have not been detected in Ae. aegypti. However, at least seven other punctual mutations have been evidenced worldwide in this mosquito, among which substitutions in the 1016 and 1534 codon sites exhibit the strongest correlation with knockdown resistance [16, 3335]. In Brazil, there are at least two kdr alleles, NaVR1 and NaVR2. The former has a mutation in the 1534 site of the channel (Phe1534Cys) and the other presents an additional substitution in the 1016 site (Val1016Ile + 1534Cys) [31, 33].

Metabolic resistance is another major physiological mechanism and refers to an increase in the synthesis of detoxifying enzymes or in their specificity to metabolize the insecticide, both resulting in an enhancement of the insect detoxifying capacity [36]. Glutathione S-transferases (GST), carboxylesterases (EST) and multi-function oxidases (MFO, also known as P450) are the key classes of enzymes generally enrolled in this process, all belonging to families composed of several genes [11, 36]. Some of these enzymes potentially act over a broad nature of xenobiotics, resulting in cross-resistance among different classes of insecticides [37, 38]. In addition to these generic activities, changes in some genes can be selected to detoxify a specific compound [39]. Among these detoxifying enzymes, the MFO class deserves attention in resistance of Ae. aegypti to pyrethroids, as revealed by high-throughput assays, by comparing the overall profile at genomic and transcriptome levels between resistant and susceptible populations [4043].

Insecticide bioassays are the primary methods for monitoring the resistance status of natural populations [44]. Further investigation of selected mechanisms may anticipate the presence of resistance alleles or specific activities before resistance reaches critical levels in a population. Based on these data, surveillance is likely to guide control actions more effectively and locally oriented. We show the temporal analysis of the pyrethroid susceptibility status of several Ae. aegypti populations from SP State along with the profiles of the main potentially related mechanisms, target site and metabolic resistance. Our aim is to describe this resistance scenario since the history of insecticide utilization in SP is different from the rest of the country, as stated above. If on one hand dengue incidence depends on many factors such as mosquito abundance, immune population frequency and vector competence of mosquitoes, on the other hand, dengue outbreaks can often lead to an uncontrolled increase in the domestic use of insecticides, available in the retail marked. The intensity of insecticide applications by governmental campaigns as well as the dynamics of Ae. aegypti resistance status in SP have been well monitored [15, 26, 45]. We analysed data related to pyrethroid resistance status, enzymatic profiles of detoxifying enzymes and the frequency of kdr alleles in Ae. aegypti populations from strategic areas of SP under a time series perspective. Although based on a regional scale, conclusions of this study may contribute to unfold the phenomena of insecticide resistance and help to design sustainable vector control strategies.

Materials and methods

Sampling

In SP, the insecticide resistance of Ae. aegypti populations from 16 municipalities is regularly monitored. In this study we focused the analysis on seven of them with different profiles of dengue incidence, an indirect indicator of insecticide use intensity. As representatives of their respective regions, Ae. aegypti samples of discriminated municipalities are listed in descending order of dengue cumulative incidence from 1998 to 2014 (cases per 100,000 inhabitants): Sao José do Rio Preto (20,203), Araçatuba (14,779), Santos (14,010), Ribeirão Preto (9,346), Campinas (5,800), Presidente Prudente (2,446) and Marília (2,084). Field sampling was performed yearly with ovitraps as described elsewhere [46] from November to December, which corresponds to the pre-epidemic season. The Rockefeller strain (Rock), continually maintained in the laboratory, was always referred to in parallel both as an insecticide susceptibility control and an internal control of laboratory assays [47].

Bioassays

Although cypermethrin was employed in SP from 1989 to 2001, in 2000 deltamethirn was implemented in the whole country, except SP. However, as the insecticide resistance monitoring evaluations were conducted at a national level, deltamethrin was the compound of choice for evaluating adult resistance to pyrethroids. The status of susceptibility/resistance to the pyrethroid deltamethrin was evaluated by dose-diagnostic (DD) mortality bioassays with impregnated papers and the WHO test tube kits according to the WHO guidelines [48, 49]. Papers were impregnated in the laboratory with the recommended diagnostic dosage of 0.05% (18 mg/a.i/m2). Around 20 non-blood fed 2–5 day old adult females were exposed during 1 hour to deltamethrin and then transferred to the resting tube free of insecticide where mortality was recorded 24 hours later. The assays for each population were performed in triplicate and were repeated at least three times with F1 or F2 females, average mortality below 80% indicating resistance, according to the WHO criteria in force at the time of the assays [48, 49]. In addition to DD tests, a dose-response (DR) assay was also performed with F1 or F2 generations of four of the populations collected in 2011. DR procedures were similar to the DD test and incorporated four to seven different deltamethrin concentrations. Data were submitted to Probit transformation followed by linear regression analysis with the help of the software Polo PC [50] in order to calculate the lethal doses (LD).

Biochemical assays

The activity of esterases (with both substrates α– and β–naphthyl acetates), multi-function oxidases (MFO) and glutathione-S-transferases (GST) enzymes was quantified in samples of around 30 individual adult females one-day post-emergence. These tests were carried out according to adaptations to the instructions from the Centers for Disease Control (CDC) and WHO [51, 52].

kdr genotyping

DNA was extracted from each of the abdomen deprived females which were tittered in 200 μL of “squishing buffer” (10 mM Tris-HCl pH 8.2, 2 mM EDTA and 0.2% Triton X-100) according to Jowett (1986), followed by an incubation with 0.2 mg/L proteinase K (Promega) at 56°C overnight and an inactivation final step at 95°C for 5 minutes. Both 1016 (Val+ and Ilekdr) and 1534 (Phe+ and Cyskdr) sites of the Ae. aegypti NaV were genotyped with a customized TaqMan genotyping assay (Thermo Fischer Scientific) independently for each site. The sequences of primers and probes are available in Table 1. The reactions were conducted in a 10 μL mix containing 1 μL of DNA (~10 ng), 1 X TaqMan Genotyping Master-mix and the TaqMan assay combined primers and probes 1X for 1016 assay or 0.5 X for 1534 assay. Around 30 specimens of each population were individually evaluated in a 96 well microplate together with the positive controls: SS (Rock strain), RR (Rock-kdr strain, [53]) and RS (an equimolar mix or Rock and Rock-kdr). The thermocycling program was in accordance with the manual instructions (TaqMan genotyping assay, Thermo Fischer Scientific) in StepOne Plus or 7500 real time machines (Thermo Fischer Scientific). As previously described [33], allelic and genotypic frequencies took into account that the 1016 and 1534 SNPs are placed in a single locus, constituted by the alleles: NaVS (1016 Val+ + 1534 Phe+), NaVR1 (1016 Val+ + 1534 Cyskdr) and NaVR2 (1016 Ilekdr + 1534 Cyskdr). Table 2 evidences how each genotype was obtained based on the result of each SNP. Hardy-Weinberg equilibrium was assessed by the classical equation [54], the null hypothesis of equilibrium checked by a chi-square test with one or three degrees of freedom when three or six genotypes, respectively, were evidenced.

thumbnail
Table 1. Sequences of primers and probes used in the genotyping reactions for 1016 and 1534 NaV sites of Aedes aegypti.

https://doi.org/10.1371/journal.pntd.0006390.t001

thumbnail
Table 2. Possible genotypes for Latin American Aedes aegypti populations, considering the single nucleotide polymorphisms in the 1016 and 1534 NaV sites.

https://doi.org/10.1371/journal.pntd.0006390.t002

Results

Bioassays

Aedes aegypti samples of seven SP municipalities, field-collected every year from 2006 to 2014, were exposed to DD bioassays with the pyrethroid deltamethrin (Table 3). In this period the criteria for interpreting bioassays with diagnostic dose classified all populations as resistant, except Campinas in 2013 and Marília in 2009 as their mortality rates exceeded 80% [48]. According to the current criterion, however, all the analysed populations would have been considered resistant to deltamethrin [55, 56]. In principle, DD tests are qualitative, i.e. useful only for classification of a population as resistant or susceptible. In addition, quantitative DR deltamethrin assays were performed with four populations collected in 2011 (Table 4), confirming the resistance status previously detected by the DD assays. Notwithstanding, we found a significant correlation between DD mortality levels and DR derived resistance ratios (RR) in the four populations submitted to both assays. In other words, the higher the RR of a population, the more specimens survived the diagnostic dose. This correlation served for both RR50 (R2 = 0.9485, p = 0.0261) and RR95 (R2 = 0.9399, p = 0.0305), suggesting that some quantitative character could be attributed in this case, comparatively, to the DD assays.

thumbnail
Table 3. Qualitative (DD) deltamethrin bioassays with F1 or F2 Aedes aegypti female adults derived from São Paulo State populations collected from 2006 to 2014.

https://doi.org/10.1371/journal.pntd.0006390.t003

thumbnail
Table 4. Quantitative (DR) deltamethrin bioassays with Aedes aegypti female adults derived from São Paulo State populations collected in 2011.

https://doi.org/10.1371/journal.pntd.0006390.t004

Due to the explosive aspect of dengue epidemics, it is possible to observe a significant intensification of the domestic use of insecticides, usually pyrethroids, during outbreaks. We investigated, indirectly, whether this indiscriminate chemical control could impact resistance to pyrethroids: dengue incidence rates were plotted against the results obtained with the bioassays with the diagnostic dose of deltamethrin. We found a positive correlation between the cumulative dengue incidence (number of cases/ 100,000 inhabitants) from 1995 to 2014 and survival to the deltamethrin DD in 2014 (R2 = 0.4196, p = 0.1157). This correlation increases (R2 = 0.7372, p = 0.0286) if we exclude Araçatuba from the analysis (Fig 1), considering that the mortality index in DD tests for this locality in 2014 substantially differed from its time series and other localities that year (Table 3).

thumbnail
Fig 1. Does the intensification of domestic insecticide use due to dengue outbreaks may have an impact on Ae. aegypti resistance to pyrethroids?

Correlation analysis between the cumulative incidence of dengue cases (1995–2014) in six localities of São Paulo State and Aedes aegypti survival after exposure to the deltamethrin qualitative dose-diagnostic bioassay of samples collected in 2014. Araçatuba was plotted in the figure but not considered in the correlation analysis (see text for details).

https://doi.org/10.1371/journal.pntd.0006390.g001

Kdr genotyping

Kdr genotyping was carried out with around 30 specimens of the seven SP localities exposed to DD, as mentioned above. In all cases, three samples were employed, collected from 2003 to 2014. The allelic composition considered the result of both genotyping reactions (for 1016 and 1534 NaV sites) for each insect, herein referred to as NaVS (1016 Val+ + 1534 Phe+), NaVR1 (1016 Val+ + 1534 Cyskdr) and NaVR2 (1016 Ilekdr + 1534 Cyskdr). The theoretically possible allele NaVR3 (1016 Ilekdr + 1534 Phe+) was not detected. Fig 2 displays the variations in the frequencies of each NaV allele and the pooled “Resistant” genotypes, indicating their respective localities in the map. Overall, there was an increase in the frequencies of the resistant kdr alleles. In this scenario, considering the difference between final and initial frequencies, the minimal observed decline of the NaVS allele in the period was 45% in Araçatuba, the greatest difference, 76%, noted in Ribeirão Preto. One must be aware that the NaVS allele was not detected in Santos since its first sampling, in 2006. Although both NaVR1 and NaVR2 kdr alleles increased in the period with differing dynamics, the latter exhibited a greater rise, lowest and highest increase being 32% (Campinas) and 58% (Ribeirão Preto), respectively. In contrast, the lowest and highest increase of NaVR1 were 5% and 30%, respectively, in Araçatuba and Santos (S1 Table).

thumbnail
Fig 2. Frequencies of susceptible (NaVS) or kdr (NaVR1 and NaVR2) alleles and resistance genotypes (R1R1, R1/R2 and R2R2) in Aedes aegypti populations from São Paulo State (SP) between 2003 and 2014.

The composition of each allele considered both 1016 and 1534 NaV sites for each genotyped individual: NaVS (1016 Val+ + 1534 Phe+), NaVR1 (1016 Val+ + 1534 Cyskdr) and NaVR2 (1016 Ilekdr + 1534 Cyskdr). The municipalities are placed inside the highlighted map of SP State.

https://doi.org/10.1371/journal.pntd.0006390.g002

Considering that knockdown resistance based on the kdr mutations is a recessive trait, individuals with the genotypes R1R1, R1R2 and R2R2 should potentially represent the pyrethroid resistant portion of the population for this target-site mechanism at least. In this case, for the most recent sampling (2014), Marília presented the lowest frequency of resistant genotypes (22%), in contrast with 100% observed in Santos and Ribeirão Preto (Table 5). All individuals from Santos were genotyped as”resistant” in the three sampling times. Meanwhile, Ribeirão Preto mosquitoes exhibited the highest rise in the resistant genotypes, from 0 in 2004 to 100% a decade later. The decrease of the S allele in all populations between 2004 and 2014 can be seen in S1 Fig.

thumbnail
Table 5. Kdr genotypes in Aedes aegypti populations from São Paulo State from 2004 to 2014.

https://doi.org/10.1371/journal.pntd.0006390.t005

Biochemical assays

Alterations in the activity of the main detoxifying enzyme classes were investigated in Ae. aegypti from all SP municipalities evaluated, collected both in 2006 and 2015 (Fig 3), as an indication for the selection of metabolic resistance mechanisms. In general, an increase in Esterases, MFO and GSTs activities was evidenced during the study time. The activities were considered either ‘altered’ or ‘highly altered’ if more than 15% or 50%, respectively, of the population presented values beyond the Rockefeller’s 99 percentile [51]. In this sense with the exception of Campinas and Marília, the remaining populations were classified as ‘altered’ for Esterases, evaluated with both αNA and βNA substrates. Araçatuba, Marília and Santos were ‘highly altered’ for MFO, presenting a substantial increase compared to the previous evaluation. Excluding Campinas, GST was ‘altered’ in all assessments and Santos (2015) ‘highly altered’.

thumbnail
Fig 3. Evaluation of enzymes related to metabolic resistance in Aedes aegypti populations from SP State in two moments from 2006 to 2015.

The bars represent the enzymatic activity in % beyond the 99 percentile of Rockefeller lineage, run in parallel (white bars = 2006 and grey bars = 2015). The values beyond 15 and 50%, indicated by the dotted lines, were respectively classified as altered or highly altered. The “X” mark indicates samples not evaluated.

https://doi.org/10.1371/journal.pntd.0006390.g003

Discussion

The susceptibility status of Aedes aegypti to insecticides has been annually monitored in SP State since 1996 by the Superintendência de Controle de Endemias (SUCEN) and accompanied as an integrant of the national network (MoReNAa Network) since 1999 [29]. Organophosphate compounds have been adopted in Brazil against larvae and adults of Ae. aegypti since the first contemporary dengue epidemic outbreak in the 1980’s. São Paulo State started employing pyrethroids in governmental campaigns targeting Ae. aegypti adults in 1989, at least ten years prior to the rest of the country. Ultralow volume (ULV) assays with Ae. aegypti populations from SP performed in 2000/2001, and afterwards between 2007 and 2009, however, revealed lack of efficacy of the pyrethroids but not of the organophosphate malathion. Based on this scenario, malathion was resumed in the State, replacing the pyrethroids [26]. Again, this change in the rest of the country started later, in 2009 [57], given the broader evidence of disseminated of pyrethroid resistance in Brazilian Ae. aegypti populations [14, 33, 58]. In this sense, laboratory bioassays have been very useful to the surveillance programs since resistance detected in these tests are correlated with failure of field control operations [18, 26]. Additional semi-field trials of residual treatment indicated lack of efficiency of cypermethrin, malathion, bendiocarb and fenitrothion in Santos but not in Marilia (2002–2003). Afterwards in 2011, this same kind of test displayed ineffectiveness of fenitrothion, deltamethrin and bendiocarb for Marilia, Araçatuba and Santos vector populations [26].

Although pyrethroids have no longer been used for several years in governmental campaigns against Aedes in SP, all vector populations here evaluated with bioassays remained resistant to deltamethrin. Araçatuba sampled in 2014 exhibited the lowest mortality level after exposure to this pyrethroid, suggesting that mosquitoes from that locality could substantially differ from the others. However, a previous population genetics study with Ae. aegypti populations from SP did not provide any evidence of an unusual genetic difference of Araçatuba insects [59]. In that study, Araçatuba, Presidente Prudente and Marília were grouped in a clade. On the other hand, insecticide residual spraying (IRS) with pyrethroids occurs since 1999 in Araçatuba, due to cases of urban visceral leishmaniosis [60]. This intensified use of pyrethroids probably contributed to the special low mortality rates of Araçatuba mosquitoes in deltamethrin bioassays, when compared to other SP localities and reinforce the need for integrated vector control methods. The increase in the frequencies of kdr alleles in all other localities suggests that a selection pressure with pyrethroids still persists. Although hard to precisely quantify, an increase in the use of household sprays, all of them formulated with pyrethroids, was noted. A tendency of intensification of metabolic resistance was also observed. However, this could also be attributed to organophosphate pressure [10]. Further bioassays with organophosphates will help us to better address this hypothesis. The lack of a well-structured integrated vector control strategy with emphasis on mechanical control and community engagement towards the prevention and elimination of breeding sites, is needed. In this sense, chemical control should be envisaged as a complementary approach regarding vector control activities, and, ideally, insecticide adoption should take into account the potential overlap of different synantropic vector species.

Herein, we followed the dynamics of the NaV alleles during a decade in seven Ae. aegypti populations of SP. The three alleles, NaVS, NaVR1 and NaVR2 previously identified in Brazilian populations [33, 61], were confirmed. As kdr mutations are physiologically expressed as a recessive trait, only the kdr homozygous insects are potentially resistant to pyrethroids [32]. Mutant alleles can remain at low levels for a long period and therefore, resistance based on this mechanism may delay to be established in the field. This happens because initially most of the mutant alleles are carried by heterozygotes, which are not favourably selected by insecticide exposure. On the other hand under similar conditions, as kdr homozygotes appear, these alleles can spread very rapidly [62]. For instance, in the present study at Ribeirão Preto, no NaVR1 or NaVR2 kdr homozygote was found in 2004. However 48% of the population presented one of these alleles as heterozygotes, meaning that none of them would be potentially resistant to pyrethroids through this target site mechanism. In 2014, 100% of genotyped Ae. aegypti from Ribeirão Preto presented a genotype with the potential for resistance (R1R1, R1R2 and R2R2).

Besides genotyping of known molecular markers for resistance, additional approaches are unraveling other factors enrolled in the physiology of metabolic resistance mechanisms. This is the case of genomic and transcriptomic analysis, based on microarrays developed to specifically assess detox genes, and high-throughput sequencing with resistant populations of Ae. aegypti. Beyond mutations on protein coding sequences, copy gene amplifications have been identified in detox genes as polymorphisms associated with pyrethroid resistance in Ae. aegypti [40, 63, 64]. These mechanisms are so dynamic and variable among populations that identifying a specific gene or a mutation as a diagnostic marker for resistance is a difficult task [65].

The so called 'biochemical tests' in turn only indicate general activity alterations of broad classes of enzymes from a test population compared to a reference lineage [66]. Based on this method, an increase in the activity of Esterases and GSTs was evident in Ae. aegypti populations from Brazil between 2001 and 2004 [52]. Herein, we also observed a great change in the MFO profile. Accordingly, the key participation of Cyp450 genes is detected by genomic and transcriptomic studies. For instance, three unrelated deltamethrin resistant Ae. aegypti populations from French Guiana (South America), Guadaloupe islands (Lesser Antilles) and New Caledonia (Pacific Ocean) shared alterations in four out of five common up-regulated Cyp450 genes [16]. Although it is not easy to trace linear correlations between enzymatic activity and mortality values, some general relationships can be considered. For instance, Santos and Ribeirão Preto presented the most altered metabolic resistance profiles in 2015 related to 2006 (Fig 3). These were also the localities where 100% of screened insects presented a kdr ‘R genotypes’ in 2014 (Table 5).

Failure in controlling vector population densities ultimately results in high incidence of the diseases derived from the arboviruses/parasites they transmit. There are several operational reasons for this flaw, insecticide resistance being an important issue when insecticides are the major tool towards vector control [10, 67]. It is not uncommon that failure of vector control with insecticides triggers the use of even more insecticides, leading to resistance increase and dissemination. When an epidemic situation is established this vicious cycle can reach extremely high levels. This is a consequence, among others, of the great increase in the domestic use of insecticides, motivated by the population's panic. This situation has been recently reported in different Brazilian municipalities [68]. The overall positive correlation found in SP State municipalities between cumulative dengue case incidence and survival levels of Ae. aegypti in the pyrethroid bioassays is an indirect evidence of this phenomenon: the greater the incidence of dengue, the more intense 'individual' attempts to control its vector; in this case with the use of insecticides for which the resistance is already installed. The consequence is evident, that is, more resistance increase. This indirect correlation, together with data demonstrating lack of insecticide efficacy in the field when laboratory bioassays had indicated loss of susceptibility [26], corroborates the importance of insecticide resistance monitoring programs, especially when the results are considered for vector control strategies.

Currently the Brazilian Dengue Control Program indicates the use of Insect Growth Regulators against larvae and the organophosphate malathion for adults in Ae. aegypti control campaigns [57], where no resistance has yet been reported in the country, except for some very specific localities. The continuous susceptibility status surveillance of field populations based on bioassays is important to evaluate the efficiency of regularly employed insecticides as well as the potential alternatives. Dose-response assays were performed in 2011, aiming to evaluate whether any quantitative character in terms of resistance level could be inferred from the dose-diagnostic tests, in principle a qualitative assay, that classifies populations as susceptible or resistant. We witnessed a significant correlation among resistance ratios (obtained from dose-response assays) and mortality values (from dose-diagnostic assays) based on the four populations assessed in both tests. Field response was also in line with the diagnostic-dose assay, where pyrethroid ULV field test applications resulted in failure for populations classified as resistant through laboratory bioassays [26]. The São Paulo State monitoring program uses only diagnostic-dose testing since it is more rapid and requires fewer insects and lower amounts of impregnated papers compared to the quantitative assays. Additionally, in the case of kdr mutations, genotyping of molecular markers for resistance to pyrethroids can predict the dynamics of the genetic background of a given population before resistance is established. However, when the frequency of the kdr allele exceeds the “tipping-point” and selection pressure persists in the environment, the resistant status is very rapidly achieved in the population [69,70].

Conclusions

Although pyrethroids have no longer been administrated by public health campaigns against Ae. aegypti in SP State since 2001, the mosquito is still resistant to this class of insecticide. Both the frequency of kdr mutations and the altered activity of detoxifying enzymes have been increasing. High dengue incidence rates in the successive epidemics throughout the country, along with a strong bias towards chemical control, are probably the main reasons motivating people to utilize excessive amounts of domestic insecticides which are mostly pyrethroids. Therefore, the selection pressure with pyrethroids continues to be relevant. Also, the insecticide based strategies for controlling other vector species, such as Culex, Anopheles, triatomines and sandflies, should be integrated considering insecticide resistance peculiarities and management of each species. Measures taken by the public health managers to avoid increase and spread of insecticide resistance and cross-resistance should also consider domestic use. The most efficient manner for vector control of dengue, chikungunya and Zika is undoubtedly avoiding accumulated water in domestic and urban containers. Nevertheless, chemical control still inspires general overconfidence in society. Community awareness regarding the consequences and limitations of the chemical strategy for controlling Ae. aegypti and towards the benefits of vector prevention measures is therefore crucial. An integrated vector control strategy is imperative for the optimal use of resources and significant amelioration in the sustainability of chemical control, not only of Aedes but also other insect vectors.

Supporting information

S1 Table. Variation between initial and final frequencies of Aedes aegypti NaV alleles from São Paulo populations.

https://doi.org/10.1371/journal.pntd.0006390.s001

(DOCX)

S1 Fig. Frequencies of the susceptible NaV allele (S: 1016 Val+ + 1534 Phe+) of Aedes aegypti populations from different municipalities of São Paulo State between 2004 and 2014.

Beyond the data points, regression lines are also exhibited.

https://doi.org/10.1371/journal.pntd.0006390.s002

(TIFF)

Acknowledgments

We thank Luana Carrara for technical assistance with molecular assays, Heloisa Diniz (IOC/Fiocruz) for taking care of the figures in this paper and Rafi Ur Rahman for critical language review. English review and revision was further conducted by Mitchell Raymond Lishon, native of Chicago, Illinois, USA–UCLA 1969.

References

  1. 1. PAHO. Annual Report 2014. Epidemiological Alerts and Updates Washington D.C.: Pan American Health Organization; 2015.
  2. 2. Faria NR, Azevedo Rdo S, Kraemer MU, Souza R, Cunha MS, Hill SC, et al. Zika virus in the Americas: Early epidemiological and genetic findings. Science. 2016;352(6283):345–9. pmid:27013429
  3. 3. WHO Director-General summarizes the outcome of the Emergency Committee regarding clusters of microcephaly and Guillain-Barré syndrome [Internet]. World Health Organization. 2016 [cited June 2016]. Available from: http://www.who.int/mediacentre/news/statements/2016/emergency-committee-zika-microcephaly/en/.
  4. 4. Heymann DL, Hodgson A, Sall AA, Freedman DO, Staples JE, Althabe F, et al. Zika virus and microcephaly: why is this situation a PHEIC? Lancet. 2016;387(10020):719–21. pmid:26876373
  5. 5. WHOPES. Pesticides and their application—For the control of vectors and pests of public health importance. Sixth ed: World Health Organization; 2006.
  6. 6. Ranson H, Burhani J, Lumjuan N, Black IV WC. Insecticide resistance in dengue vectors. TropIKAnet. 2010;1(1).
  7. 7. Brasil. Diretrizes nacionais para prevenção e controle de epidemias de dengue. Brasília: Ministério da Saúde, Secretaria de Vigilância em Saúde, Departamento de Vigilância Epidemiológica; 2009.
  8. 8. Braga IA, San Martin JL. Histórico do Controle de Aedes aegypti. In: Fiocruz , editor. Dengue: Teorias e práticas. Rio de Janeiro: Editora Fiocruz; 2015.
  9. 9. WHO. Vector resistance to pesticides. Control WECoVBa, editor. Geneva: World Health Organization; 1992. 62 p.
  10. 10. Moyes CL, Vontas J, Martins AJ, Ng LC, Koou SY, Dusfour I, et al. Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans. PLoS Negl Trop Dis. 2017;11(7):e0005625. pmid:28727779
  11. 11. Valle D, Belinato TA, Martins AJ. Controle Químico de Aedes aegypti, resistência a inseticidas e alternativas In: Fiocruz , editor. Dengue: Teorias e práticas. Rio de Janeiro: Editora Fiocruz; 2015. p. 93–126.
  12. 12. WHO. Safety of pyrethroids for public health use: World Health Organization; 2005.
  13. 13. Luna JE, Martins MF, Anjos AF, Kuwabara EF, Navarro-Silva MA. [Susceptibility of Aedes aegypti to temephos and cypermethrin insecticides, Brazil]. Revista de saude publica. 2004;38(6):842–3. pmid:15608903
  14. 14. 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(4):441–4. pmid:16113895
  15. 15. Macoris Mde L, 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(8):895–900. pmid:18209926
  16. 16. Dusfour I, Zorrilla P, Guidez A, Issaly J, Girod R, Guillaumot L, et al. Deltamethrin Resistance Mechanisms in Aedes aegypti Populations from Three French Overseas Territories Worldwide. PLoS Negl Trop Dis. 2015;9(11):e0004226. pmid:26588076
  17. 17. Kushwah RB, Dykes CL, Kapoor N, Adak T, Singh OP. Pyrethroid-resistance and presence of two knockdown resistance (kdr) mutations, F1534C and a novel mutation T1520I, in Indian Aedes aegypti. PLoS Negl Trop Dis. 2015;9(1):e3332. pmid:25569164
  18. 18. Marcombe S, Darriet F, Tolosa M, Agnew P, Duchon S, Etienne M, et al. Pyrethroid resistance reduces the efficacy of space sprays for dengue control on the island of Martinique (Caribbean). PLoS Negl Trop Dis. 2011;5(6):e1202. Epub 2011/06/30. pmid:21713017
  19. 19. Chuaycharoensuk T, Juntarajumnong W, Boonyuan W, Bangs MJ, Akratanakul P, Thammapalo S, et al. Frequency of pyrethroid resistance in Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in Thailand. J Vector Ecol. 2011;36(1):204–12. pmid:21635659
  20. 20. Braga IA, Martins JLS. Histórico do controle de Aedes aegypti. In: Valle D, Pimenta DN, Cunha RV, editors. Dengue—teorias e práticas Rio de Janeiro: Ed Fiocruz; 2015. p. 61–73.
  21. 21. Glasser CM, Gomes AdC. Infestação do Estado de São Paulo por Aedes aegypti e Aedes albopictus. Revista de saude publica. 2000;34:570–7. pmid:11175600
  22. 22. Teixeira MG, Costa MCN, Barreto ML, Barreto FR. Epidemiologia da Dengue. In: Valle D, Pimenta DN, Cunha RV, editors. Dengue—teorias e práticas. Rio de Janeiro: Ed Fiocruz; 2015. p. 293–315.
  23. 23. São Paulo. Secretaria de Estado de Saúde. Dados estatísticos dengue, Zika e chikungunya SP: Governo do Estado de São Paulo; 2017. Available from: http://portal.saude.sp.gov.br/cve-centro-de-vigilancia-epidemiologica-prof.-alexandre-vranjac/areas-de-vigilancia/doencas-de-transmissao-por-vetores-e-zoonoses/agravos/dengue/dados-estatisticos.
  24. 24. Zanluca C, Melo VC, Mosimann AL, Santos GI, Santos CN, Luz K. First report of autochthonous transmission of Zika virus in Brazil. Mem Inst Oswaldo Cruz. 2015;110(4):569–72. pmid:26061233
  25. 25. Braga IA, Valle D. Aedes aegypti: histórico do controle no Brasil. Epidemiol Serv Saúde. 2007;16(2):113–8. http://dx.doi.org/10.5123/S1679-49742007000200006
  26. 26. Macoris MGL, Andrighetti MTM, Wanderley DMV, Ribolla PEM. Impact of insecticide resistance on the field control of Aedes aegypti in the State of São Paulo. Rev Soc Bras Med Tro. 2014;47(5):573–8.
  27. 27. Andrighetti MT, Macoris MLG, Nalcon KCR, Otrera VCG. Monitoramento da suscetibilidade de Aedes aegypti aos inseticidas utilizados para seu controle no Estado de São Paulo. Boletim Epidemiol Paulista (BEPA) 2006;3(suppl I):50–5.
  28. 28. Macoris MdLdG, Andrighetti MTM, Otrera VCG, Carvalho LRd, Caldas Júnior AL, Brogdon WG. Association of insecticide use and alteration on Aedes aegypti susceptibility status. Mem I Oswaldo Cruz. 2007;102:895–900.
  29. 29. Braga IA, Valle D. Aedes aegypti: vigilância, monitoramento da resistência e alternativas de controle no Brasil. Epidemiol Serv Saúde. 2007;16:295–302.
  30. 30. Martins AJ, Valle D. The Pyrethroid Knockdown Resistance. In: Soloneski S, Larramendy M, editors. Insecticides—Basic and Other Applications Rijeka: InTech; 2012. p. 17–38.
  31. 31. Seixas G, Salgueiro P, Silva AC, Campos M, Spenassatto C, Reyes-Lugo M, et al. Aedes aegypti on Madeira Island (Portugal): genetic variation of a recently introduced dengue vector. Mem Inst Oswaldo Cruz. 2013;108 Suppl 1:3–10. pmid:24473797
  32. 32. Dong K, Du Y, Rinkevich F, Nomura Y, Xu P, Wang L, et al. Molecular biology of insect sodium channels and pyrethroid resistance. Insect Biochem Mol Biol. 2014;50:1–17. pmid:24704279
  33. 33. Linss JG, Brito LP, Garcia GA, Araki AS, Bruno RV, Lima JB, et al. Distribution and dissemination of the Val1016Ile and Phe1534Cys Kdr mutations in Aedes aegypti Brazilian natural populations. Parasit Vectors. 2014;7(1):25. Epub 2014/01/17. pmid:24428880
  34. 34. Saavedra-Rodriguez K, Urdaneta-Marquez L, Rajatileka S, Moulton M, Flores AE, Fernandez-Salas I, et al. A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti. Insect Mol Biol. 2007;16(6):785–98. Epub 2007/12/21. pmid:18093007
  35. 35. Harris AF, Rajatileka S, Ranson H. Pyrethroid resistance in Aedes aegypti from Grand Cayman. Am J Trop Med Hyg. 2010;83(2):277–84. Epub 2010/08/05. pmid:20682868
  36. 36. Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000;45:371–91. Epub 2000/04/13. pmid:10761582
  37. 37. Hemingway J, Hawkes NJ, McCarroll L, Ranson H. The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004;34(7):653–65. Epub 2004/07/10. pmid:15242706
  38. 38. Montella IR, Schama R, Valle D. The classification of esterases: an important gene family involved in insecticide resistance—A review. Mem Inst Oswaldo Cruz. 2012;107(4):437–49. Epub 2012/06/06. pmid:22666852
  39. 39. Ranson H, Claudianos C, Ortelli F, Abgrall C, Hemingway J, Sharakhova MV, et al. Evolution of supergene families associated with insecticide resistance. Science. 2002;298(5591):179–81. pmid:12364796
  40. 40. Strode C, Wondji CS, David JP, Hawkes NJ, Lumjuan N, Nelson DR, et al. Genomic analysis of detoxification genes in the mosquito Aedes aegypti. Insect Biochem Mol Biol. 2008;38(1):113–23. pmid:18070670.
  41. 41. Marcombe S, Poupardin R, Darriet F, Reynaud S, Bonnet J, Strode C, et al. Exploring the molecular basis of insecticide resistance in the dengue vector Aedes aegypti: a case study in Martinique Island (French West Indies). Bmc Genomics. 2009;10:494. Epub 2009/10/28. pmid:19857255; PubMed Central PMCID: PMC2770535.
  42. 42. Bariami V, Jones CM, Poupardin R, Vontas J, Ranson H. Gene amplification, ABC transporters and cytochrome P450s: unraveling the molecular basis of pyrethroid resistance in the dengue vector, Aedes aegypti. PLoS Negl Trop Dis. 2012;6(6):e1692. pmid:22720108
  43. 43. Faucon F, Dusfour I, Gaude T, Navratil V, Boyer F, Chandre F, et al. Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing. Genome Res. 2015;25(9):1347–59. pmid:26206155
  44. 44. WHO. Monitoring and managing insecticide resistance in Aedes mosquito populations—Interim guidance for entomologists. Geneva: World Health Organization; 2016.
  45. 45. Macoris M, Andrighetti MT, Takaku L, Glasser CM, Garbeloto VC, Cirino VC. [Changes in susceptibility of Aedes aegypti to organophosphates in municipalities in the state of Sao Paulo, Brazil]. Revista de saude publica. 1999;33(5):521–2. pmid:10576756
  46. 46. Lima JB, Da-Cunha MP, Da Silva RC, Galardo AK, Soares Sda S, Braga IA, et al. 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(3):329–33. pmid:12685640
  47. 47. Kuno G. Early history of laboratory breeding of Aedes aegypti (Diptera: Culicidae) focusing on the origins and use of selected strains. J Med Entomol. 2010;47(6):957–71. pmid:21175042
  48. 48. WHO. Instructions for determining the suscetibility or resistance of adult mosquitoes to organochlorine, organophosphate and carbamate insecticides. Stablishment of the baseline. Geneva: World Health Organization; 1981. p. 4.
  49. 49. WHO. Criteria and Meaning of Tests for Determining the susceptibility or Resistance of Insects to insecticides. Geneva: World Health Organization; 1981. p. 4.
  50. 50. Raymond M. Presentation d’une programme d’analyse logprobit pour microordinateur cahiers Orstrom. Sér Ent Med Parasitol. 1985;22:117–21.
  51. 51. Valle D, Montella IR, Medeiros PFV, Ribeiro RA, Martins AJ, Lima JBP. Quantification metodology for enzyme activity related to insecticide resistance in Aedes aegypti. Brasília: Ministério da Saúde/ Brasil; 2006.
  52. 52. 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(3):467–77. pmid:17827362
  53. 53. Brito LP, Linss JG, Lima-Camara TN, Belinato TA, Peixoto AA, Lima JB, et al. Assessing the effects of Aedes aegypti kdr mutations on pyrethroid resistance and its fitness cost. PloS one. 2013;8(4):e60878. pmid:23593337
  54. 54. Shorrocks B. The Genesis of Diversity. London: Hodder and Stoughton; 1978.
  55. 55. WHO (2016a). Test procedures for insecticide resistance monitoring in malaria vector mosquitoes - 2nd ed. Edited by Organization WH. Geneva, Switzerland: World Health Organization; 2016: 48.
  56. 56. WHO (2016b). Monitoring and managing insecticide resistance in Aedes mosquito populations: Interim guidance for entomologists. In: WHO/ZIKV/VC/161: World Health Organization; 2016.
  57. 57. Brasil. [Informa mudança de uso de inseticidas larvicidas e adulticidas na rotina do Programa Nacional de Controle da Dengue]. Brasília: Ministério da Saúde, Secretaria de Vigilância em Saúde, 2009.
  58. 58. Bellinato DF, Viana-Medeiros PF, Araujo SC, Martins AJ, Lima JB, Valle D. Resistance Status to the Insecticides Temephos, Deltamethrin, and Diflubenzuron in Brazilian Aedes aegypti Populations. BioMed research international. 2016;2016:8603263. pmid:27419140
  59. 59. Santos VMd Macoris MdLdG, Andrighetti MTM Avila PE, Kirchgatter K. Analysis of genetic relatedness between populations of Aedes aegypti from different geographic regions of São Paulo state, Brazil. Revista do Instituto de Medicina Tropical de Sao Paulo. 2003;45:99–101. pmid:12754576
  60. 60. Cardim MF, Rodas LA, Dibo MR, Guirado MM, Oliveira AM, Chiaravalloti-Neto F. Introduction and expansion of human American visceral leishmaniasis in the state of Sao Paulo, Brazil, 1999–2011. Revista de saude publica. 2013;47(4):691–700. pmid:24346660
  61. 61. Chapadense FG, Fernandes EK, Lima JB, Martins AJ, Silva LC, Rocha WT, et al. Phenotypic and genotypic profile of pyrethroid resistance in populations of the mosquito Aedes aegypti from Goiania, Central West Brazil. Rev Soc Bras Med Trop. 2015;48(5):607–9. pmid:26516973
  62. 62. Barbosa S, Black WCt, Hastings I. Challenges in estimating insecticide selection pressures from mosquito field data. PLoS Negl Trop Dis. 2011;5(11):e1387. pmid:22069506
  63. 63. David JP, Coissac E, Melodelima C, Poupardin R, Riaz MA, Chandor-Proust A, et al. Transcriptome response to pollutants and insecticides in the dengue vector Aedes aegypti using next-generation sequencing technology. Bmc Genomics. 2010;11:216. Epub 2010/04/02. pmid:20356352
  64. 64. Faucon F, Dusfour I, Gaude T, Navratil V, Boyer F, Chandre F, et al. Unravelling genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing. Genome Res. 2015.
  65. 65. Saavedra-Rodriguez K, Suarez AF, Salas IF, Strode C, Ranson H, Hemingway J, et al. Transcription of detoxification genes after permethrin selection in the mosquito Aedes aegypti. Insect Mol Biol. 2012;21(1):61–77. pmid:22032702
  66. 66. Hemingway J. Techniques to Detect Insecticide Resistance Mechanisms (field and laboratory manual). Organization WH, editor: World Health Organization; 1998. 35 p.
  67. 67. Mnzava AP, Knox TB, Temu EA, Trett A, Fornadel C, Hemingway J, et al. Implementation of the global plan for insecticide resistance management in malaria vectors: progress, challenges and the way forward. Malaria journal. 2015;14:173. pmid:25899397
  68. 68. Garcia GA, David MR, Martins AJ, Freitas RM, Linss JGB, Araujo SC, et al. (2018) The impact of insecticide applications on the dynamics of resistance: The case of four Aedes aegypti populations from different Brazilian regions. PLoS Negl Trop Dis 12(2): e0006227. https://doi.org/10.1371/journal.pntd.0006227 pmid:29432488
  69. 69. WHO (2012). Global Plan for Insecticide Resistance Management in Malaria Vectors (GPIRM). Geneva, Switzerland: World Health Organization; 2012: 132.
  70. 70. Du Y, Nomura Y, Zhorov BS, Dong K. Sodium Channel Mutations and Pyrethroid Resistance in Aedes aegypti. Insects. 2016;7(4). pmid:27809228; PubMed Central PMCID: PMCPMC5198208.
  71. 71. Vera-Maloof FZ, Saavedra-Rodriguez K, Elizondo-Quiroga AE, Lozano-Fuentes S, Black Iv WC. Coevolution of the Ile1,016 and Cys1,534 Mutations in the Voltage Gated Sodium Channel Gene of Aedes aegypti in Mexico. PLoS Negl Trop Dis. 2015;9(12):e0004263. pmid:26658798