Skip to main content
Advertisement
  • Loading metrics

Laboratory selection of Aedes aegypti field populations with the organophosphate malathion: Negative impacts on resistance to deltamethrin and to the organophosphate temephos

  • Priscila Fernandes Viana-Medeiros,

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

    Affiliation Laboratório de Biologia Molecular de Flavivírus, Instituto Oswaldo Cruz (Fiocruz), Rio de Janeiro (RJ), Brasil

  • Diogo Fernandes Bellinato,

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

    Affiliation Laboratório de Biologia Molecular de Flavivírus, Instituto Oswaldo Cruz (Fiocruz), Rio de Janeiro (RJ), Brasil

  • Denise Valle

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

    dvalle@ioc.fiocruz.br

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

Abstract

Background

Resistance to pyrethroids and to the organophosphate temephos is widespread in Brazilian populations of the dengue vector, Aedes aegypti. Thereof, since 2009 Insect Growth Regulators are employed as larvicides, and malathion is used against adults.

Methodology/Principal findings

We performed laboratory selection with malathion of two A. aegypti field populations initially susceptible to this organophosphate but resistant to temephos and deltamethrin. A fixed malathion dose inducing at least 80% mortality in the first generation, was used throughout the selection process, interrupted after five generations, when the threshold of 20% mortality was reached. For each population, three experimental and two control groups, not exposed to insecticides, were kept independently. For both populations, quantitative bioassays revealed, in the selected groups, acquisition of resistance to malathion and negative impact of malathion selection on deltamethrin and temephos resistance levels. In the control groups resistance to all evaluated insecticides decreased except, unexpectedly, to deltamethrin. Analysis of the main resistance mechanisms employed routine methodologies: biochemical and molecular assays for, respectively, metabolic resistance and quantification of the NaV pyrethroid target main kdr mutations at positions 1016 and 1534. No diagnostic alteration could be specifically correlated with malathion selection, neither with the unusual deltamethrin increase in resistance levels observed in the control groups.

Conclusions/Significance

Our results confirm the multifactorial character of insecticide resistance and point to the need of high throughput methodologies and to the study of additional field vector populations in order to unravel resistance mechanisms.

Author summary

Dengue, Zika and chikungunya viruses affect millions of people worldwide. Due to the lack of specific antivirals or to the limited supply of vaccines, focus remains on the control of the main vector, Aedes aegypti. Although the importance of social participation in the elimination of A. aegypti breeding sites is increasingly recognized, chemical control is still an important component of vector control. The exaggerated use of insecticides results in the spread of resistance and, consequently, in the loss of their effectiveness. In Brazil, malathion is the last adulticide available to the control of A. aegypti, due to the widespread resistance to pyrethroids. In order to anticipate what could occur in the field, we exposed two vector populations to selection with malathion. Both malathion and temephos, a larvicide largely employed, are organophosphates; however, they are structurally distinct molecules and seem to elicit different resistance mechanisms. We confirmed this issue: selection with malathion had a negative impact on temephos resistance compared to groups reared without any insecticide. Indeed, the variety of responses of both vector populations to the various insecticides points to the participation of multiple resistance mechanisms and confirms previous assumptions regarding the difficulty of identifying diagnostic insecticide resistance mechanisms.

Introduction

Aedes aegypti (Linnaeus, 1762) is the main vector of dengue, Zika and chikungunya viruses, which are currently serious public health problems in Brazil [1]. Specific medications are not available for any of these viruses. One dengue vaccine was recently approved but it confers only partial protection [2]; other vaccines against dengue and Zika viruses are yet in the clinical trials phase [35], and efforts in this regard for the chikungunya virus are just beginning [6]. Moreover, since 2016, Brazil faces a public health emergency with respect of yellow fever [7]. Although its transmission is silvatic, by Haemagogus and Sabethes mosquitoes [8], with the low vaccination coverage of the population there is an imminent risk of reurbanization of the virus, and A. aegypti is its main urban vector [9,10].

Except for the yellow fever virus, in the current scenario prophylactic measures depend on control initiatives aimed at decreasing vector density [11,12]. Although the most recommended and effective measures are the mechanical control and society awareness, chemical insecticides are still widely used tools. However, the intense exposure of vector and pest populations to insecticides can select characteristics and mechanisms that confer resistance to these xenobiotics [13]. Vector populations resistant to insecticides are a threat to the success of control programs that prioritize chemical control [14,15]. Thus, monitoring the susceptibility status of natural populations to the main employed insecticides is one of the pillars of vector control programs. This is a way of guiding control actions and preventing resistance of reaching a threshold considered of risk [12,16].

The organophosphate (OP) temephos has been used in the control of A. aegypti larvae in Brazil since 1967. Resistance to temephos, monitored in vector populations throughout the country, started in 2000–2001 and increases since then. In 2009, due to the spread of resistance [17], temephos ceased to be the first-choice larvicide, being replaced by chitin synthesis inhibitors (CSI) [15,18].

Similar to the control of larvae, OP were also used against adult mosquitoes up to 2001, when pyrethroids (PY) replaced them [19]. However, shortly thereafter, reduction of PY susceptibility of various field populations started to be identified [20,21], a process that resulted in the decision of interrupting the use of this class of insecticides in the control of A. aegypti adults in 2009 [18,22].

The Brazilian Ministry of Health (MoH) adopts only insecticides recommended by WHO [23] for the control of A. aegypti. Taking into account resistance of adults to PY insecticides, the only viable alternative for ultra low volume applications was malathion [23]. Although malathion is also an OP, such as temephos, the structure of both compounds is distinct: while temephos is a closed-chain phosphorothionate, malathion, an open-chain compound, is a phosphorodithioate, and bears two sulfur atoms attached to the central phosphorus, rather than just one [24]. Regarding the perspective of vector control, such differences are important, since they contribute to elicit different resistance mechanisms, as is the case with Anopheles mosquitoes, for example [25,26]. Indeed, in Brazil, in contrast to the widespread A. aegypti resistance to temephos, registers of altered malathion susceptibility were scarce, and only derived from qualitative bioassays [27,28]. When PY were replaced by malathion in the control of adults, quantitative assays confirmed the susceptibility of Brazilian natural A. aegypti populations to this OP (Braga TA, personal communication), corroborating the pertinence of its use in the field.

Taking this scenario into account, the anticipation of a potential malathion resistance event, as well as the elucidation of the possible mechanisms involved, would be of great value for the definition of rational strategies that could prolong this OP as a viable alternative for the control of A. aegypti adult specimens in Brazil. This study describes the selection with malathion, in the laboratory, of two A. aegypti field populations resistant to temephos and to the PY deltamethrin. Alterations in the malathion susceptibility profile were investigated, as well as the consequences of malathion selection on the susceptibility profile of the resulting samples to the main insecticides used in the country against the vector: temephos, the CSI diflubenzuron and deltamethrin. Collection of A. aegypti population from the municipality of Aracaju, at Sergipe State (SE) was made in 2012 when the MoH applied both temephos and CSI against larvae, and only PY to control adults. Crato mosquitoes, from Ceará State (CE), the higher resistant population, was collected in 2013 when temephos applications had just been replaced by CSI, and control of adults employed both PY and malathion. More details regarding the history of insecticides use in both localities are shown elsewhere [29]. Our major aim was to investigate if selection with malathion for a limited number of generations would be enough to alter the susceptibility profile of the exposed populations. In addition, it would be of interest to evaluate to what extent different vector populations, bearing distinct backgrounds regarding challenges with insecticides and other xenobiotics, would be impacted by malathion selection.

Methods

Mosquito lines

Sampling of A. aegypti eggs representative of the evaluated municipalities was done with ovitraps [27,30]. In the laboratory, mosquito colonies were started with 349 positive ovitraps from Aracaju / SE installed in 2012, and 272 positive ones from Crato / CE, installed in 2013 [29]. Egg hatching and larval rearing were carried out as routinely [27]. The resulting A. aegypti adults were then identified and used to obtain F1 eggs. Although the initial objective was to employ F1 specimens to characterize the resistance and the resistance mechanisms of the original populations [29] and also to proceed with the selection, the amount of eggs obtained was not enough. Therefore, F2 generation specimens were used in the selection procedure.

Insecticides

In addition to malathion selection, the susceptibility status to the main insecticides applied by the Brazilian Dengue Control Program was analyzed. The study adopted technical grade compounds: the OP temephos (Pestanal; Sigma-Aldrich Brasil Ltda, São Paulo, SP, Brazil) and malathion (Cheminova Brasil Ltda, São Paulo, SP, Brazil), the PY deltamethrin (Pestanal; Sigma-Aldrich Brasil Ltda) and the CSI diflubenzuron (Pestanal; Sigma-Aldrich Brasil Ltda).

Malathion selection

A malathion dose capable of killing about 80% of specimens in the first generation was used. This same dose was kept fixed until the mortality was only 20%. Three biological replicates of the selection (S1, S2 and S3) were done per population, without exchange of specimens between the replicates. For each population two control groups (C1 and C2) were also included, kept simultaneously in the same conditions, but without any contact with insecticides.

Malathion concentrations were defined in previous trials: 0.06229 mg/L for Aracaju and 0.09876 mg/L for Crato. For the first generation of selection, at least 5,000 larvae of each population were used per replicate, and 3,000 specimens were used to initiate subsequent generations of each independent experimental group. Synchronized larvae were used in the whole selection process: hatching of A. aegypti eggs was induced for at most one hour in a BOD incubator at 28°C. Larvae were then kept at 26 ± 2°C, at the density of 500 larvae per basin, in 1 L dechlorinated water and 1 g of grounded cat food (Friskies, Purina/Camaquã/RS) was provided. Under these conditions, larvae collected 72 hours after hatching were essentially L3 ones, and approximately 5% of the specimens had moulted to L4.

Basically, for each strain and each replica, several groups of 100 larvae, 72 hours old, were placed in plastic cylindrical basins (10 X 12 cm) with 250 ml of dechlorinated water and malathion ethanolic solution (for Aracaju and Crato samples, respectively 346 and 548.6 μl of a 45 mg/l malathion solution were employed) or equivalent volume of ethanol, in the case of control replicas. Larvae were fed with 0.2 g of grounded cat food added each three days. Exposure to malathion was continuous up to pupation. Pupae were daily transferred to cylindrical cardboard cages (17 X 18 cm). Upon eclosion, adults were fed ad libitum (except before the blood meal) with a 10% sugar solution, replaced three times a week and kept in a temperature and humidity-controlled insectary (26 ± 1°C; 80 ± 10% rh).

To obtain eggs, females were weekly deprived of the sugar solution for 18–24 hours, and anaesthetized guinea pigs were then offered as blood source according to the “Formulary for laboratory animals” [31].

After five generations, all replicates exposed to malathion exhibited less than 20% mortality in the presence of the insecticide, and the experiment was interrupted.

Bioassays

Quantitative bioassays of larvae with temephos followed the parameters and procedures described by WHO [32]. Bioassays with diflubenzuron were made according to Martins et al. [33]. Adult bioassays with PY and malathion were done according to an adaptation of the insecticide impregnated papers methodology [34,35]. In all cases, three to four assays were performed for each population on different days. As an internal control, simultaneous trials were also carried out with the Rockefeller strain, a reference of insecticide susceptibility [36]. In each bioassay A. aegypti specimens were exposed to an insecticide spectrum of at least six concentrations. Larvae were exposed to temephos for 24 hours and to diflubenzuron until death or emergence of the last adult specimen of the control group, as previously defined by Braga et al. [37].

Adults were exposed to either deltamethrin or malathion for 1 hour and then transferred to a recovery chamber, with no insecticide, where they remained for another 24 hours, when mortality was recorded. For each assay with larvae, 320 to 640 specimens were used, totaling 960–1,920 specimens per insecticide tested by biological replicate (S1-S3, C1, C2). In the case of adults, 360–480 specimens were used per assay, corresponding to 1,080–1,440 specimens for each insecticide per biological replicate.

Genotyping assays

The TaqMan method was used to identify kdr mutations in the NaV gene in the post-selection samples, as previously described for pre-selection populations [29]. For each replica (C1, C2, S1, S2, S3), 30 individual males were used. Analysis of Aracaju samples was done with F7 specimens (offspring of the last generation of selection in the replicates S1-S3). In the case of Crato, analysis of C1 and C2 samples was done with F6, and S1-S3 replicates were investigated with F7 mosquitoes. Independent reactions were done for the substitutions Val1016Ile and Phe1534Cys. For each position, 1 μl, equivalent to 0.5% of the DNA content extracted from each specimen, was used in a 10 μl final reaction volume.

Biochemical assays

One-day-old adult females were evaluated according to methodology adapted from WHO and the US Centers for Disease Control and Prevention (CDC) [3840]. The same protocol, with some additional adaptations, was applied to larvae [41]. Activity of mixed function oxidases (MFO), esterases (EST) and glutathione S-transferases (GST) were quantified. Three substrates were used for ESTs: α-naphthyl, β- naphthyl and ρ-nitrophenyl acetates, related respectively to α-EST, β-EST and ρNPA-EST activities. For Acetylcholinesterase (Ace), the OP target site, both total activity (AChE) and the remaining activity after inhibition by propoxur (AChI) were evaluated. In order to calculate the specific enzymatic activities, total protein content of each sample was quantified using the Bio-Rad protein reagent (catalog number: 500–0006).

At least three assays were accomplished on different days for all enzymes. In each assay 40 individual specimens of each experimental replica were evaluated. Simultaneously, as an internal control of the assays, five specimens of the Rockefeller strain were also analyzed.

Data analysis

In the case of bioassays, adult emergence inhibition (EI) or lethal concentrations (LC) were calculated, respectively for diflubenzuron and neurotoxic insecticides, by probit analysis using the Polo-PC software (LeOra Software, Berkeley, CA) [42]. Results of the quantitative assays were expressed as the ratio of resistance between the LC (or EI) of the experimental group under test and the Rockefeller equivalent measure (RR). Another index, 'selection ratio' (SR), was also calculated, by comparison of LC (or EI) values of C and S groups with their corresponding parental (P) population ones [29]. Additionally, we assessed the 95% confidence intervals overlap range between the post- and pre-selection LC of each experimental sample.

The criterion used in Brazil to classify the resistance status of A. aegypti populations to temephos [43] was also applied to the other insecticides evaluated in this study. According to this criterion, populations with RR95 higher than 3.0 are classified as resistant, and there is a recommendation for replacement of the insecticide compound used in the field.

Individual NaV genotypes, as well as the allelic and genotypic frequencies of each population, were calculated based on the variations in positions 1016 and 1534, both on the same gene, as described elsewhere [44]. The 95% confidence intervals overlap range between the post- and pre-selection kdr allelic frequencies was also evaluated.

Enzymatic activities were classified according to a previously established criterion, based on the use of the 99th percentile of a reference strain as the cutoff point. In this case, references corresponded to the parental strains, from Aracaju and Crato [29]. After calculation of their 99th percentiles for each enzyme class activity, rates of corresponding specimens in the C and S groups above this value were estimated. The activities were classified as unaltered, altered or highly altered in cases where these rates were respectively less than 15, between 15 and 50, or above 50% [43,38]. In particular, in the case of the Ace inhibition assay with propoxur (AChI), the WHO criterion establishes that remaining activity exceeding 30% is indicative of resistance to OP insecticides [40]. All enzyme activity profiles were also compared using the Kruskal Wallis nonparametric test and the Dunn's multiple comparison post test (using the graphpad prism version 5.0).

Ethics statement

The blood feeding was done according to the Brazilian guidelines described in “The Brazilian legal framework on the scientific use of animals” [45], supported by a protocol approved by the "Ethics Committee in the Study of Animals" (CEUA/Fiocruz 2008), licenses L-011/09 and LW-20/14.

Results

Selection

In the first generation of selection with a fixed malathion dose, respectively 20.4 and 11.8% of Aracaju and Crato individuals survived. In the course of the selection process, the Aracaju survival percentage increased more evenly than that of Crato. Nevertheless, both populations reached the threshold of 80% survival (dotted line in Fig 1), defined as the time of selection interruption, in the 5th generation (F6). For the population of Crato, exposure to malathion was repeated with one additional generation, when 90.9% of survival was detected. After that, the potential impact of laboratory selection on the resistance status of malathion itself (in larvae and adults) and on other compounds used in the A. aegypti control in the country was assessed. It is ought to mention that mortality in the control groups did not exceed 1% during the whole process.

thumbnail
Fig 1. Survival rate of two Aedes aegypti populations subjected to selection with a fixed dose of malathion.

Exposure to the insecticide was continuous, between L3 and pupal phases. The x-axis indicates the generations of laboratory rearing.

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

Insecticide resistance: Bioassays

Table 1 summarizes the impacts of malathion selection and of adaptation to laboratory breeding on the RR profile to several insecticides employed in A. aegypti control. In order to facilitate a direct comparison, results of the original evaluations with the pre-selected populations [29] were included. Bold values in Table 1 point to resistance when the Brazilian MoH classification criteria, previously defined for temephos, is used (RR95 > 3.0, see Methods). Detailed results are presented in S1S5 Tables, including comparison of mortality values obtained for C and S groups with their corresponding parental (P) populations (the 'selection ratio', or SR analysis).

thumbnail
Table 1. Resistance status of A. aegypti larvae and adults from Aracaju/SE and Crato/CE to the insecticides used by the Brazilian MoH.

Samples before and after the malathion selection process are shown.

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

Both analytical approaches pointed out the effectiveness of the laboratory selection promoted in the study. Although SR data are slighter than RR ones: (a) both malathion RR and SR increased in the larvae of virtually all selected groups—malathion SR values increased in the range of 50% and more than twice for, respectively, Aracaju and Crato selected larvae; (b) in general, there was no overlap of CL values between C and S groups, pointing to a sound effect of malathion selection (S1 Fig); (c) for Crato, the selective pressure reverberated in the adult stage to the point that it could also be classified as resistant after selection, when the 'MoH RR95 criterion' was used. Malathion SR values increased in the range of 50% for Crato adults. Groups maintained without insecticide showed a slight reduction of malathion RR in the larval stage (which were already low), and a discrete increase in the adult stage, while SR revealed a slight increase in both stages (S1 and S2 Tables). In general, for each population, C replicas, as well as S ones, had an equivalent performance, although reared independently.

Laboratory rearing in the absence of insecticides (groups C1, C2) was enough to slightly reduce the initial temephos resistance levels exhibited by parental populations. Selection with malathion did not induce cross-resistance with temephos, although both are OP (Tables 1 and S3). In contrast, temephos resistance decrease was more noticeable in malathion selected groups than in the control ones–both RR and SR indexes suggested such a trend. In particular, temephos resistance reduction was more pronounced in Crato, the population with the highest initial RR.

When the present study started, diflubenzuron had been recently introduced in the country to control A. aegypti larvae. Resistance indexes for this CSI compares adult emergency inhibition rates (EI), and not lethal concentrations. When compared to the original populations, although some differences in RR and SR for DFB were observed between C and S groups, all values obtained were considered compatible with a susceptible status (Tables 1 and S4).

All samples were considered highly resistant to deltamethrin. Surprisingly, the simple laboratory rearing, in the absence of any insecticide (groups C1, C2), exacerbated PY resistance levels in both populations (Tables 1 and S5). In particular, the SR index confirmed the marked deltamethrin resistance increase in the control groups, mainly in Aracaju (S5 Table). This increase in deltamethrin resistance, however, was not observed in the S1-3 groups, selected with malathion, independently of the analytical method adopted (RR, SR and 95% confidence limits overlapping). In opposition, malathion selected groups, and particularly Crato mosquitoes, exhibited reduction of deltamethrin resistance levels.

In general, when compared to the parental samples and as judged by slope values, homogeneity of all evaluated groups remained unaltered or increased (S1S5 Tables). Only with deltamethrin bioassays there was a trend towards increased heterogeneity (except for Crato S groups). In this regard, the unusual and specific increase in deltamethrin resistance, noted in the control groups, should be considered.

Insecticide resistance mechanisms: kdr mutations genotyping

Table 2 shows the kdr genotypes at positions 1016 and 1534 of the NaV gene, the target site of PY insecticides. The 'wild type' susceptible allele, '1016 Val+ + 1534 Phe+', was called 'S'. The allele mutated only at position 1534 (1016 Val+ + 1534 Cyskdr) was named 'R1' and the other, with mutations at both positions (1016 Ilekdr + 1534 Cyskdr), 'R2'. Mutation exclusively at position 1016 was not found. The three alleles were detected in Aracaju and Crato. Since kdr mutations at positions 1016 and 1534 are recessive, PY resistance mediated by such changes is only expressed in homozygous individuals: R1R1, R1R2 and R2R2 [44] (highlighted in bold and summed in column ‘resistant genotypes’ in Table 2).

thumbnail
Table 2. Changes in the voltage-gated sodium channel (AaNaV), target of pyrethroids, in Aracaju and Crato, during the malathion selection process.

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

In Crato the wild-type (S) allele remained as the most frequent, followed by R1, a pattern kept in all sample groups from that locality. On the other hand, R2 was the most frequent allele in Aracaju, followed by the S allele but this arrangement tended to reverse after laboratory rearing. Regarding the original populations, in both cases, an increase in the frequency of the S allele was observed in almost all samples, selected with malathion or simply reared in the laboratory, without insecticides (note the IC 95% overlap range in S2 Fig).

The initial kdr genotype frequencies, exhibiting PY resistance, were 43% for Aracaju and 32% for Crato ('resistant genotypes' column in Table 2). In almost all cases, in both groups C and S, such frequencies decreased during the selection or rearing processes.

Insecticide resistance mechanisms: Metabolic resistance and OP target

Tables 3, 4 and S7 show the activity profiles of the major classes of detoxification enzymes in larvae and adults of Aracaju and Crato mosquito populations. Acetylcholinesterase activity, target site of OP, was also quantified; in this case the total activity (AChE) and its profile of inhibition by propoxur (AChI) are presented. Tables 3 and 4 depict, respectively, the percentage of Aracaju and Crato individuals with activity above the 99th percentile of each corresponding Parental strain; S7 Table compares median values of the same samples, using the Kruskal-Wallis nonparametric test.

thumbnail
Table 3. Quantification of detoxifying enzymes and acetylcholinesterase activity in Aracaju A. aegypti exposed to selection with malathion—Comparison with the 99th percentile of the parental population.

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

thumbnail
Table 4. Quantification of detoxifying enzymes and acetylcholinesterase activity in Crato A. aegypti exposed to selection with malathion—Comparison with the 99th percentile of the parental population.

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

For both populations more changes and higher activities were noted in adult mosquitoes, compared to the larval stage, in all C and S groups. The exceptions were α-EST and β-EST activities, which are higher in larvae (S7 Table).

In general, no differences were observed between control groups and those exposed to malathion which could be classified as diagnostic of the selected resistance.

In relation to the parental samples, Aracaju adult females tended to present higher MFO and GST activities, the later being revealed by both criteria in control samples (Tables 3 and S7) but detected only after comparison of median values in S ones (S7 Table). In Crato, marked alterations of all classes of enzymes were also observed in both C and S females, when compared to the Parental population; however, GST alterations were more attenuated in relation to the other enzyme classes, being detected only by medians comparison (Tables 4 and S7).

In the larval stage, no metabolic alterations were detected in Aracaju C and S samples, compared to the original population. In Crato larvae, α-EST EST activity was consistently enhanced and MFO alterations were noted in some samples when median values were compared, particularly in the malathion selected ones.

Total activity of the OP target site, AChE, remained equivalent to the original samples. We also investigated Ace using the WHO susceptibility criterion: inhibition of 70% or more of Ace's activity with propoxur (AChI); in all cases the results are indicative of a susceptible enzyme (Tables 3, 4 and S7). According to this criterion, Ace of both parental Aracaju and Crato populations has a susceptible profile.

Discussion

This study deals with the response to laboratory rearing and to selection with malathion of two A. aegypti field populations. Their resistance status to the main insecticides employed by the Brazilian MoH as well as the potentially associated resistance mechanisms were evaluated.

The outcome of malathion selection was considered consistent, since resistance to this OP increased in the three biological replicates of both Aracaju and Crato. This increase was seen in all the adopted analytical approaches—comparison with the Rockefeller susceptible strain or with the parental generation. Although comparison with the parental population (SR values) showed more slight increases than comparison with the susceptible reference strain (RR), both selection efficacy and the functional significance of this result are confirmed since there is no overlap of CL ranges between the selected groups and the control samples, nor with the parental ones. Considering also the initial and final mortalities of the malathion selection experiment, which ranged from 80 to less than 20% in the six independently selected groups (Fig 1), we can infer that malathion selection appears to have had a relevant biological impact.

The magnitude of such increases in resistance levels was also compatible with other records in the literature that subjected A. aegypti populations to different insecticides for a small number of generations [46,47].

Although selection has increased malathion resistance levels, resistance to temephos tended to decrease in both populations. Such reduction was more prominent in the samples selected with malathion than in the control groups. This situation was, to some extent, expected, considering the molecular structures of both insecticides and other reports in the literature [25,4850]. Collectively, such data reinforce the idea that these two organophosphates elicit different mechanisms of resistance [25,50,51]. This also occurs, for example, with species of the genera Anopheles and Culex: resistance to temephos has been found to be derived from the copy number amplification of esterases genes, culminating in a greater number of available enzyme molecules [51,52]; on the other hand, malathion metabolic resistance is related to alterations in the coding region of a particular esterase, called 'malaoxonase', resulting in a more efficient enzyme [26,53].

In fact, the two OPs differ in chemical structure, as stated in the Introduction. Such differences appear to result in the extremely slow hydrolysis of temephos, via esterases, a process that may last for days. In this case, it is even considered that temephos is sequestered [48,50]. As a consequence, esterase-mediated resistance to OP compounds, in general, would occur due to an increase in the amount of enzyme molecules, as is the case with Culex mosquitoes [51,54]. In our study, alterations in some specific resistance mechanisms were observed. However, such changes were not able to elucidate all the complexity observed in the responses of the different sample groups of Aracaju and Crato, related to the evaluated insecticides.

In the samples here evaluated, changes in the activity of enzymes related to metabolic resistance were noted, particularly in adult specimens of the vector. However, no difference was observed that was exclusive to the replicates submitted to selection, in both populations. This result suggests the participation of different molecular species and points to the need of high-performance methodological approaches in the identification of resistance mechanisms. However, this type of evaluation is certainly beyond the scope of the tests used in routine resistance monitoring [12,43,55,56].

Regarding Ace, target of OP insecticides, no significant changes were noted neither in total activity (AChE) nor in inhibition of activity (AChI). In general, AChI results were concordant for the two classification criteria employed, WHO [40] and Valle et al. [38]. Therefore, the study of the OP target revealed no alterations likely to explain the decrease in malathion susceptibility in the samples exposed to this insecticide.

Malathion selection appears to have had a negative impact on deltamethrin resistance, similar to what happened with temephos status of malathion selected mosquitoes. In both cases, temephos and deltamethrin, the RR reduction was more prominent in the malathion selected groups than in the control ones. One potential reason is the selection of malathion resistance mechanisms, in detriment of others, previously present in the original samples and related to resistance to temephos and deltamethrin. In other words, in the absence of changes in Ace, the target of OP, this scenario may result from the deviation of metabolic resources, in the scope of resistance mediated by detoxifying enzymes [57,58]. It is worth mentioning that decrease of PY resistance levels in the malathion selected groups was followed by the maintenance (in Crato) ou slight reduction (in Aracaju) of the kdr resistant homozygotes frequency.

Unexpectedly, a marked increase in deltamethrin resistance levels was observed in the control groups, maintained without any insecticide, relative to the original samples, in both populations. This occurred despite a slight decrease in the kdr frequencies, a known mechanism of PY resistance. This decrease in kdr mutations rate is somewhat expected since there are indications that, in the absence of PY, kdr mutations do not offer an adaptive advantage; instead, this feature is related to a significant evolutionary cost [59]. Of course, the possibility of contamination cannot be completely ruled out, although the whole rearing process has been tracked, as well all bioassays steps. However, taken together, such data suggest that (a) the kdr mutation frequency assessment should not be used as a single marker of PY resistance, at least in A. aegypti, and (b) additional mechanisms capable of conferring PY resistance were selected during the laboratory procedures.

In general, selection experiments performed under laboratory conditions tend to privilege multiple low-effect mechanisms, as would be expected in the field, according to assumptions of natural selection [60]. However, in the field, application of insecticides seems to be just one amongst the many challenges faced by mosquitoes. Aedes aegypti, as a generalist insect, can also deal with the presence of pollutants or secondary plant metabolites in its larval breeding sites. This condition could bestow a pre-adaptive advantage regarding the insecticide challenge and justify the early detection of resistance to newly introduced compounds, for instance, in the context of vector control programs [20,21]. In fact, some metabolites and pollutants were able to induce a significant increase in the expression of MFO enzymes in A. aegypti, corroborating such pre-adaptability for insecticide resistance [61].

As expected, higher changes in resistance levels were observed in the larval stage, submitted to selection. However, the metabolic resistance components evaluated were more altered in adults. The methodology used here for the quantification of resistance mechanisms is the approach of choice in the routine of the resistance monitoring of vector populations. However, this methodology has known limitations: it uses general substrates, which reveal classes of enzymes and not specific molecular entities. This means that eventual changes in some particular enzyme of a given class can be masked by the pool of enzymes acting on that substrate. In addition, it is possible that other mechanisms, not investigated by this approach, are in operation.

Evaluation of metabolic activities in adult Aracaju females indicated a potential participation of GST enzymes in the increased deltamethrin resistance: after comparison with both criteria, parental 99th percentile and the nonparametric statistics, higher GST activity was detected in Aracaju samples kept without insecticide—precisely those with a marked increase of PY resistance levels. Some studies suggest relation between GST and PY resistance of A. aegypti populations from different geographic origins [43,62,63]; others have even characterized this association [49,64]. However, it is worth noting that in adult females from Crato control groups, although also bearing an increased deltamethrin resistance, the change in GST was not prominent, being only detected through medians comparison. It is well known that insecticide resistance has a multifactorial character. There are several recent examples and there is growing understanding that different populations of the same vector species can select different resistance mechanisms against the same class of insecticides. This can derive from several parameters, intrinsic and extrinsic. Among the former are the genetic background and previously selected mechanisms, while insecticide selection pressure intensity and sources are among the extrinsic factors [15,46,63,65,66].

Although malathion is an adulticide, we opted to employ larvae in the selection procedure for operational reasons. On one hand, we were interested in carrying out this analysis in a timely manner to respond to the management of resistance in the country; on the other hand however, the technical difficulty of calibrating selection trials with adults is greater than with larvae. Selection with larvae was then a cost-effective way in which we decided to invest in order to ensure that all specimens would be exposed to malathion in a way as controlled as possible. Culicidae males are smaller and develop faster than females. In this case, the ideal selection with adults would require different experiments, with distinct malathion concentrations for males and females. Furthermore, in bioassays with adults, mainly with the available impregnated surfaces, it is difficult to control the level of exposure, since contact with insecticides occurs only during landing. Several studies report the use of larvae in selection experiments with adulticides and, in general, this apparent inconsistency is not even mentioned in the texts (examples in [6778]). Despite this, we are aware that the selection of larvae with malathion here shown may not represent the field situation reliably, and the results may have limited value in defining resistance management recommendations. For example, there is at least one report of the development of resistance to the adulticide deltamethrin that was more effective in larval selection, compared to adults [79]. In another study, selection of adult females with the same PY also resulted in only mild resistance status alteration of a Brazilian A. aegypti population [80].

It is possible that adaptation to laboratory rearing, in the absence of insecticides, has favored, or selected, one or more mechanisms relevant to the response to PY. Unlike changes in NaV (kdr), such mechanisms could have borrowed advantage to specimens carrying them. This might have occurred even in the absence of selective pressure and under optimal rearing conditions. If this is the case, malathion selection may have hampered the expression of such mechanisms in some way: lower levels of PY resistance were detected in the selected specimens compared to the control groups.

In general, it is expected that laboratory rearing, without selective pressure (as was the case with the control groups herein), will increase homogeneity regarding insecticide susceptibility [81]. However, our data indicate that this relationship, between absence of selective pressure and reduction of heterogeneity, although frequent, cannot be generalized. Slope is often a parameter put aside in the insecticide resistance monitoring of disease vectors. Few studies mention the heterogeneity of the evaluated samples, or often fail to inform details (reports we detected that mention slope results: [27,28,37,56,8288]). In the present study, homogeneity increase was noted in most groups, regardless the insecticide challenge, probably due to rearing under laboratory-controlled conditions. The main exception occurred with deltamethrin evaluation: control groups of both populations presented greater heterogeneity related to the PY. And it was precisely for deltamethrin that an unusual and marked resistance increase was detected in these control groups.

Laboratory selection of samples from two Brazilian A. aegypti natural populations with the OP malathion was attained. In both cases, alterations in the resistance profiles related to other insecticides employed in the country against A. aegypti were also observed. No cross-resistance was detected between malathion and temephos, also an OP, nor between malathion and deltamethrin, a PY. On one hand, malathion selection had a negative effect on resistance to other insecticides. On the other hand, maintenance of both populations in an insecticide-free environment resulted in an unexpected increase in resistance to deltamethrin. However, the methodological procedures routinely adopted by the control programs did not reveal the potential resistance mechanisms associated, in each case. This scenario suggests the pertinence of adopting high throughput investigation approaches which we are, in fact, providing with the material presented here.

Our data point to the feasibility of applying diflubenzuron and malathion in A. aegypti control in Brazil, respectively against larvae and adults. A larvicidal rotation scheme, as already recommended by the MoH, could certainly contribute to preserve diflubenzuron (or another CSI) in the field as a viable strategy. Unfortunately, a similar approach with malathion is not possible since all the alternative available adulticides are PY compounds [23], and A. aegypti resistance against this class of insecticides is widespread throughout the country [20,21,29,44,56]. Finally, considering the insecticides employed by control programs, our results reinforce the relevance of investigating the susceptibility status, and the associated mechanisms, of more field A. aegypti populations. This policy could subsidize the rational choice of the compounds and collaborate to keep the feasibility of the few available compounds.

Supporting information

S1 Fig. Malathion susceptibility levels (mg/l).

For each evaluated population and development stage, the 95% confidence limits of both lethal concentrations, LC50 and LC95, were plotted in order to estimate overlapping ranges among samples reared in the laboratory (C, S for control and malathion selected, respectively) and the parental (P) ones.

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

(TIF)

S2 Fig. kdr allelic frequencies.

The 95% confidence limit of each one of the alleles S (1016 Val + 1534Phe), R1 (1016 Val + 1534Cys) and R2 (1016Ile + 1534Cys) was plotted for both Aracaju and Crato mosquitoes S: green, R1: yellow, R2: pink. In the x-axis, P, C1-2 and S1-3 refer to, respectively, Parental, Control and Selected samples.

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

(TIF)

S1 Table. Details of malathion bioassays performed with Aedes aegypti larvae samples.

Results generated by probit analyses. Results of bioassays with the Rockefeller strain tested simultaneously to the experimental and control groups are also shown. SR: ‘Selection Rate’, ratio between the LC (or EI) values of C and S groups and their corresponding parental population (P). LC: lethal concentration. The 95% confidence interval is shown.

https://doi.org/10.1371/journal.pntd.0006734.s003

(PDF)

S2 Table. Details of malathion bioassays performed with Aedes aegypti adult samples.

Legend as in S1 Table.

https://doi.org/10.1371/journal.pntd.0006734.s004

(PDF)

S3 Table. Details of temephos bioassays performed with Aedes aegypti larvae samples.

Legend as in S1 Table.

https://doi.org/10.1371/journal.pntd.0006734.s005

(PDF)

S4 Table. Details of diflubenzuron bioassays performed with Aedes aegypti larvae samples.

Legend as in S1 Table. EI: Emergence inhibition.

https://doi.org/10.1371/journal.pntd.0006734.s006

(PDF)

S5 Table. Details of deltamethrin bioassays performed with Aedes aegypti adult samples.

Legend as in S1 Table.

https://doi.org/10.1371/journal.pntd.0006734.s007

(PDF)

S6 Table. Allelic and genotypic frequencies of A. aegypti NaV gene at positions 1016 and 1534, shown separately.

https://doi.org/10.1371/journal.pntd.0006734.s008

(PDF)

S7 Table. Quantification of detoxifying enzymes and acetylcholinesterase activity in Aracaju and Crato A. aegypti exposed to selection with malathion—Comparison of median values.

Samples' median values were compared with the corresponding medians of parental strains [29] using the Kruskal-Wallis test. (*): significantly different values (p<0.01). n1: number of individuals considered in the analysis. med2: median of enzymatic activities.

https://doi.org/10.1371/journal.pntd.0006734.s009

(PDF)

Acknowledgments

The authors thank all the team of Laboratório de Fisiologia e Controle de Artrópodes Vetores, Instituto Oswaldo Cruz (Fiocruz), Rio de Janeiro (RJ), Brasil, for logistical support, and especially Luana Carrara da Silva and Raquel Santos de Oliveira, for technical assistance. We also acknowledge Dr Ademir de Jesus Martins for assistance with the statistical assessments. Three anonymous reviewers are also appreciated for their pertinent comments.

References

  1. 1. Ferreira-de-Brito A, Ribeiro IP, Miranda RM, Fernandes RS, Campos SS, Silva KA et al. First detection of natural infection of Aedes aegypti with Zika virus in Brazil and throughout South America. Mem Inst Oswaldo Cruz. 2016; 111(10): 655–658. pmid:27706382
  2. 2. Sanofi Pasteur (2016). Comunicado de imprensa: a vacina contra dengue da Sanofi Pasteur é aprovada no Brasil. Available in: http://www.sanofi.com.br/l/br/medias/524EE272-23C5-4319-B849-70B1D9ECBFD0.pdf (Accessed on 25 January 2016).
  3. 3. Schwartz LM, Halloran ME, Durbin AP, Longini IM Jr. The dengue vaccine pipeline: implications for the future of dengue control. Vaccine. 2015; 33: 3293–3298. pmid:25989449
  4. 4. Khetarpal N, Khanna I. Dengue fever. causes, complications, and vaccine strategies. J Immunol Res. 2016; Article ID 6803098, 14 pages.
  5. 5. NIAID/NIH (2017). Zika virus vaccines. Available in: https://www.niaid.nih.gov/diseases-conditions/zika-vaccines (Acessed on 19 July 2017).
  6. 6. Erasmus JH, Auguste AJ, Kaelber JT, Luo H, Rossi SL, Fenton K et al. A chikungunya fever vaccine utilizing an insect-specific virus platform. Nat Med. 2016; 23(2): 192–199. pmid:27991917
  7. 7. Secretariat of Health Surveillance—Ministry of Health, 2017. Centro de operações de emergências em saúde pública sobre febre amarela: Informe 23/2017. http://portalarquivos.saude.gov.br/images/pdf/2017/fevereiro/24/COES-FEBRE-AMARELA-INFORME23-Atualizacao-23fev2017-13h.pdf (Acessed on 17 February 2017).
  8. 8. Goldani LZ. Yellow fever outbreak in Brazil, 2017. Braz J Infect Dis. 2017; 21(2): 123–124. pmid:28336123
  9. 9. Bonaldo MC, Gómez MM, dos Santos AA, Abreu FVS, Ferreira-de-Brito A, Miranda RM et al. Genome analysis of yellow fever virus of the ongoing outbreak in Brazil reveals polymorphisms. Mem Inst Oswaldo Cruz. 2017; 112(6): 447–451. pmid:28591405
  10. 10. Gómez MM, Abreu FVS, dos Santos AAC, Mello IS, Santos MP, Ribeiro IP et al. Genomic and structural features of the Yellow Fever virus from the 2016–2017 Brazilian outbreak. Preprint. Available from: bioRxiv 179481; https://doi.org/10.1101/179481.
  11. 11. Townson H, Nathan MB, Zaim M, Guillet P, Manga L, Bos Ret al. Exploiting the potential of vector control for disease prevention. Bull World Health Organ: Int J Public Health. 2005; 83(12): 942–947.
  12. 12. Ranson H, Burhani J, Lumjuan N, Black WC. Insecticide resistance in dengue vectors. TropIKA.net Journal. 2010; 1(1).
  13. 13. Beaty BJ, Maquardt WC. The biology of disease vectors. 1st Ed. Colorado: University Press of Colorado; 1996.
  14. 14. Poupardin R, Reynaud S, Strode C, Ranson H, Vontas J, David JP. Cross-induction of detoxification genes by environmental xenobiotics and insecticides in the mosquito Aedes aegypti: impact on larval tolerance to chemical insecticides. Insect Biochem Mol Biol. 2008; 38(5): 540–551. pmid:18405832
  15. 15. Valle D, Belinato TA, Martins AJ. Controle químico de Aedes aegypti, resistência a inseticidas e alternativas. In: Valle D, Pimenta DN, Cunha RV, Editors. Dengue: teorias e práticas. Rio de Janeiro: Editora Fiocruz; 2015. p. 93–126.
  16. 16. 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
  17. 17. Chediak M, Pimenta FG Jr, Coelho GE, Braga IA, Lima JBP, Cavalcante KRLJ et al. Spatial and temporal country-wide survey of temephos resistance in Brazilian populations of Aedes aegypti. Mem Inst Oswaldo Cruz. 2016; 111(5): 311–321. pmid:27143489
  18. 18. Secretariat of Health Surveillance—Ministry of Health, 2009. Manejo da resistência de Aedes aegypti a inseticidas. Nota Técnica n∘ 013/2009 CGPNCD/DIGES/SVS/MS. Ministério da Saúde, Brasília.
  19. 19. Braga IA, Valle D. Aedes aegypti: histórico do controle no Brasil. Epidemiol Serv Saúde. 2007; 16: 113–118.
  20. 20. 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. pmid:16113895
  21. 21. 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. pmid:19735371
  22. 22. Fiocruz (2011). Recomendação técnica sobre a interrupção do uso de inseticidas piretróides no controle do Aedes aegypti no Brasil. Nota Técnica 2/2011/IOC-Fiocruz/Diretoria, Instituto Oswaldo Cruz, Rio de Janeiro, Brazil. Available in: http://www.fiocruz.br/ioc/media/NotaTecnica%202%20IOC%20v1%20dengue%206%20set%202011.pdf
  23. 23. WHO (2012). World Health Organization: WHOPES recommended insecticides for space spraying. Geneva, Switzerland. http://www.who.int/whopes/Insecticides_for_space_spraying_Jul_2012.pdf (Accessed on 7 October 2015).
  24. 24. Bloomquist JR. Insecticides: Chemistries and Characteristics. In: University of Minnesota's electronic textbook of Integrated Pest Management; 2009. Available in: http://ipmworld.umn.edu/chapters/bloomq.htm
  25. 25. Hemingway J, Ranson H. Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000; 45: 371–391. pmid:10761582
  26. 26. Karunaratne SH, Hemingway J. Malathion resistance and prevalence of the malathion carboxylesterase mechanism in populations of mosquito vectors of disease in Sri Lanka. Bull World Health Org. 2001; 79(11): 1060–1064. pmid:11731814
  27. 27. Lima JB, da-Cunha MP, Da Silva RC, Galardo AK, Soares SS, Braga IA et al. Resistance of Aedes aegypti to organophosphates in several municipalities in the State of Rio de Janeiro and Espírito Santo, Brazil. Am J Trop Med Hyg. 2003; 68: 329–333. pmid:12685640
  28. 28. 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(8): 895–900. pmid:18209926
  29. 29. Viana-Medeiros PF, Bellinato DF, Martins AJ, Valle D. Insecticide resistance, associated mechanisms and fitness aspects in two Brazilian Stegomyia aegypti (= Aedes aegypti) populations. Med Vet Entomol. 2017; 31: 340–350 pmid:28752548
  30. 30. Braga IA, Gomes AC, Nelson M, Mello RCG, Bergamaschi DP, Souza JMP. Comparative study between larval surveys and ovitraps to monitor populations of Aedes aegypti. Rev Soc Bras Med Trop. 2000; 33: 347–353. pmid:10936947
  31. 31. Hawk CT, Leary SL. Formulary for laboratory animals. 2nd Ed. Iowa: University Press; 1999.
  32. 32. WHO. Guidelines for laboratory and field testing of mosquito larvicides. World Health Organization, 2005, Geneva, Switzerland (WHO/CDS/WHOPES/GCDPP/2005.13). Available: http://apps.who.int/iris/bitstream/10665/69101/1/WHO_CDS_WHOPES_GCDPP_2005.13.pdf. Accessed 17 May 2016.
  33. 33. Martins AJ, Belinato TA, Lima JBP, Valle D. Chitin synthesis inhibitor effect on Aedes aegypti populations susceptible and resistant to organophosphate temephos. Pest Manag Sci. 2008; 64: 676–680. pmid:18257034
  34. 34. World Health Organization (2006). World Health Organization: Guidelines for testing mosquito adulticides for indoor residual spraying and treatment of mosquito nets. Geneva, Switzerland (WHO/CDS/NTD/WHOPES/GCDPP/3).
  35. 35. Maciel-de-Freitas R, Avendanho FC, Santos R, Sylvestre G, Araújo SC, Lima JBP et al. Undesirable consequences of insecticide resistance following Aedes aegypti control activities due to a dengue outbreak. PLoS ONE. 2014; 9: e92424. pmid:24676277
  36. 36. 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: 957–971. pmid:21175042
  37. 37. Braga IA, Mello CB, Peixoto AA, Valle D. Evaluation of methoprene effect on Aedes aegypti (Diptera: culicidae) development in laboratory conditions. Mem Inst Oswaldo Cruz. 2005; 100(4): 435–440. pmid:16113894
  38. 38. Valle D, Montella IR, Viana-Medeiros PF, Ribeiro RA, Martins AJ, Lima JBP (2006). Quantification methodology for enzyme activity related to insecticide resistance in Aedes aegypti. Ministry of Health, Brasília. Available: http://bvsms.saude.gov.br/bvs/publicacoes/manual_novo_protocolo_dengue.pdf. Accessed 26 May 2018.
  39. 39. Centers for Disease Control and Prevention (1998) Insecticide Resistance Workshop: Microplate Assays. CDC, Atlanta, GA.
  40. 40. Hemingway J (1998). Techniques to detect insecticide resistance mechanisms (field and laboratory manual). (WHO/CDC/CPC/MAL/98.6). World Health Organization, Geneva.
  41. 41. Viana-Medeiros PF. Investigações sobre os mecanismos de resistência em larvas e adultos de Aedes aegypti, Linnaeus, 1762. MSc. Thesis, Fundação Oswaldo Cruz, Instituto Oswaldo Cruz. Rio de Janeiro, Brazil. 2011.
  42. 42. Raymond M. Presentation d’une programme d’analyse logprobit pour microordinateur. Cahiers Orstrom Sér Ent Med Parasitol. 1985; 22: 117–121.
  43. 43. 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–477. pmid:17827362
  44. 44. Linss JGB, 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. Parasite Vectors. 2014; 7: 25.
  45. 45. Filipecki AT, Machado CJ, Valle S, Teixeira MO. The Brazilian legal framework on the scientific use of animals. ILAR Journal. 2011; 52: E8–15. pmid:21447857
  46. 46. Saavedra-Rodriguez K, Suarez AF, Salas IF, Strode C, Ranson H et al. Transcription of detoxification genes after permethrin selection in the mosquito Aedes aegypti. Insect Mol Biol. 2012; 21(1): 61–77. pmid:22032702
  47. 47. Belinato TA, Valle D. The impact of selection with diflubenzuron, a chitin synthesis inhibitor, on the fitness of two Brazilian Aedes aegypti field populations. PLoS ONE. 2015; 10: e0130719. pmid:26107715
  48. 48. Hemingway J, Karunaratne SH. Mosquito carboxylesterases: a review of the molecular biology and biochemistry of a major insecticide resistance mechanism. Med Vet Entomol. 1998; 12(1): 1–12. pmid:9513933
  49. 49. Li X, Schuler MA, Berenbaum MR. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu Rev Entomol. 2007; 52: 231–253. pmid:16925478
  50. 50. 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–449. pmid:22666852
  51. 51. Hemingway J, Hawkes NJ, McCarroll L, Ranson H. The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004; 34: 653–665. pmid:15242706
  52. 52. Scott JA. The molecular genetics of resistance: resistance as a response to stress. Florid Entomol. 1995; 78(3): 399–414.
  53. 53. Whyard S, Downe AER, Walker VK. Characterization of a novel esterase conferring insecticide resistance in the mosquito Culex tarsalis. Arch Insect Biochem Physiol. 1995; 29: 329–342. pmid:7655057
  54. 54. Hawkes NJ, Hemingway J. Analysis of the promoters for the beta-esterase genes associated with insecticide resistance in the mosquito Culex quinquefasciatus. Biochem Biophys Acta. 2002; 1574(1): 51–62. pmid:11955613
  55. 55. Polson KA, Brogdon WG, Rawlins SC, Chadee DD. Characterization of insecticide resistance in Trinidadian strains of Aedes aegypti mosquitoes. Acta Trop. 2011; 117(1): 31–38. pmid:20858454
  56. 56. Bellinato DF, Viana-Medeiros PF, Araújo SC, Martins AJ, Lima JBP, Valle. Resistance status to the insecticides temephos, deltamethrin, and diflubenzuron in Brazilian Aedes aegypti populations. BioMed Res Int. 2016; Article ID 8603263, 12 pages.
  57. 57. Carrière Y, Deland J-P, Roff DA, Vincent C. Life-history costs associated with the evolution of insecticide resistance. Proc R Soc Lond B. 1994; 258: 35–40.
  58. 58. Rivero A, Vézilier J, Weill M, Read AF, Gandon S. Insecticide control of vector-borne diseases: when is insecticide resistance a problem? PLoS Pathogens. 2010; 6: e1001000. pmid:20700451
  59. 59. Brito LP, Linss JGB, Lima-Camara TN, Belinato TA, Peixoto AA, Lima JBP et al. Assessing the effects of Aedes aegypti kdr mutations on pyrethroid resistance and its fitness cost. PLoS ONE. 2013; 8: e60878. pmid:23593337
  60. 60. Mallet J. The evolution of insecticide resistance: Have the insects won? Trends Ecol Evol. 1989; 4(11): 336–340. pmid:21227375
  61. 61. Poupardin R, Reynaud S, Strode C, Ranson H, Vontas J, David JP. Cross-induction of detoxification genes by environmental xenobiotics and insecticides in the mosquito Aedes aegypti: impact on larval tolerance to chemical insecticides. Insect Biochem Mol Biol. 2008; 38(5): 540–551. pmid:18405832
  62. 62. Strode C, de Melo-Santos M, Magalhães T, Araújo A, Ayres C. Expression profile of genes during resistance reversal in a temephos selected strain of the dengue vector, Aedes aegypti. PLoS ONE. 2012; 7(8): e39439. pmid:22870187
  63. 63. Faucon F, Gaude T, Dusfour I, Navratil V, Corbel V, Juntarajumnong W et al. In the hunt for genomic markers of metabolic resistance to pyrethroids in the mosquito Aedes aegypti: an integrated next-generation sequencing approach. PLoS Negl Trop Dis. 2017; 11(4): e0005526. pmid:28379969
  64. 64. Vontas JG, Small GJ, Hemingway J. Glutathione S-transferases as antioxidant defence agents confer pyrethroid resistance in Nilaparvata lugens. Biochem J. 2001; 357: 65–72. pmid:11415437
  65. 65. 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–123. pmid:18070670
  66. 66. 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; 25(9): 1347–1359. pmid:26206155
  67. 67. González T, Bisset JA, Díaz C, Rodríguez MM, Diéguez L. The evolution of resistance in a Culex quinquefasciatus strain starting from selection with the pyrethroid insecticide lambdacyhalothrin. Rev Cubana Med Trop. 1996; 48(3): 218–223. pmid:9805056
  68. 68. Bisset JA, Rodríguez MM, Díaz C, Soca A. Resistance in a strain of Culex quinquefasciatus coming from Medellín, Colombia. Rev Cubana Med Trop. 1998; 50 (2):133–137. pmid:10349433
  69. 69. Hidayati H, Sofian MA, Nazni WA, Lee HL. Insecticide resistance development in Culex quinquefasciatus (Say), Aedes aegypti (L.) and Aedes albopictus (Skuse) larvae against malathion, permethrin and temephos. Trop Biomed. 2005; 22(1): 45–52. pmid:16880753
  70. 70. Rodríguez 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(4): 437–445. pmid:16506569
  71. 71. Selvi S, Edah MA, Nazni WA, Lee H-L, Tyagi BK, Sofian-Azirun M et al. Insecticide susceptibility and resistance development in malathion selected Aedes albopictus (Skuse). Trop Biomed. 2010; 27: 534–550. pmid:21399596
  72. 72. Bisset J, Rodríguez M, Soca A. Cross-resistance to malathion in Cuban Culex quinquefasciatus induced by larval selection with deltamethrin. Med Vet Entomol. 1998 Jan;12(1):109–12. pmid:9513948
  73. 73. Chandre F, Darriet F, Darder M, Cuany A, Doannio JM, Pasteur N, Guillet P. Pyrethroid resistance in Culex quinquefasciatus from west Africa. Med Vet Entomol. 1998 Oct;12(4):359–66. pmid:9824819
  74. 74. Hamdan H, Sofian-Azirun M, Nazni W, Lee HL. Insecticide resistance development in Culex quinquefasciatus (Say), Aedes aegypti (L.) and Aedes albopictus (Skuse) larvae against malathion, permethrin and temephos. Trop Biomed. 2005 Jun;22(1):45–52. pmid:16880753
  75. 75. Nazni WA, Lee HL, Sa'diyah I. Rate of resistance development in wild Culex quinquefasciatus (Say) selected by malathion and permethrin. Southeast Asian J Trop Med Public Health. 1998 Dec;29(4):849–55. pmid:10772575
  76. 76. Ramkumar G, Shivakumar MS. Laboratory development of permethrin resistance and cross-resistance pattern of Culex quinquefasciatus to other insecticides. Parasitol Res. 2015 Jul;114(7):2553–60. Epub 2015 Apr 10. pmid:25855351
  77. 77. Sarkar M, Bhattacharyya IK, Borkotoki A, Baruah I, Srivastava RB. Development of physiological resistance and its stage specificity in Culex quinquefasciatus after selection with deltamethrin in Assam, India. Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 104(5): 673–677, August 2009.
  78. 78. Selvi S, Edah MA, Nazni WA, Lee HL, Azahari AH. Resistance development and insecticide susceptibility in Culex quinquefasciatus against selection pressure of malathion and permethrin and its relationship to crossresistance towards propoxur Tropical Biomed 2005. 22(2): 103–113.
  79. 79. Kumar S, Thomas A, Sahgal A, Verma A, Samuel T, Pillai MK. Effect of the synergist, piperonyl butoxide, on the development of deltamethrin resistance in yellow fever mosquito, Aedes aegypti L. (Diptera: Culicidae). Arch Insect Biochem Physiol. 2002; 50(1): 1–8. pmid:11948970
  80. 80. Martins AJ, Ribeiro CD, Bellinato DF, Peixoto AA, Valle D, Lima JB. Effect of insecticide resistance on development, longevity and reproduction of field or laboratory selected Aedes aegypti populations. PLoS One. 2012; 7: e31889. pmid:22431967
  81. 81. Ferrari JA. Insecticide resistance. In: Beaty BJ, Marquardt WC, Editors. The biology of disease vectors. Niwot, CO: University Press of Colorado; 1996. p. 512–529.
  82. 82. Braga IA, Lima JB, Soares S da S, 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 (2): 199–203. pmid:15250476
  83. 83. 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. pmid:19857255
  84. 84. Alvarez LC, Ponce G, Oviedo M, Lopez B, Flores AE. Resistance to malathion and deltamethrin in Aedes aegypti (Diptera: Culicidae) from western Venezuela. J Med Entomol. 2013; 50(5): 1031–1039. pmid:24180108
  85. 85. Araújo AP, Araujo Diniz DF, Helvecio E, de Barros RA, de Oliveira CM, Ayres CF et al. The susceptibility of Aedes aegypti populations displaying temephos resistance to Bacillus thuringiensis israelensis: a basis for management. Par Vect. 2013; 6: 297.
  86. 86. Bisset JA, Marín R, Rodríguez MM, Severson DW, Ricardo Y, French L et al. Insecticide resistance in two Aedes aegypti (Diptera: Culicidae) strains from Costa Rica. J Med Entomol. 2013; 50: 352–361. pmid:23540124
  87. 87. Bisset JA, Rodríguez MM, French L, Severson DW, Gutiérrez G, Hurtado D et al. Insecticide resistance and metabolic mechanisms involved in larval and adult stages of Aedes aegypti insecticide-resistant reference strains from Cuba. J Am Mosq Control Assoc. 2014; 30(4): 298–304. pmid:25843136
  88. 88. Ishaaya I, Horowitz AR. Advances in insect control and resistance management: an overview. In: Horowitz AR, Ishaaya I, Editors. Advances in insect control and resistance management. Switzerland: Springer International Publishing; 2016. p. 1–8.