Worldwide the mosquito Aedes aegypti (L.) is the principal urban vector of dengue viruses. Currently 2.5 billion people are at risk for infection and reduction of Ae. aegypti populations is the most effective means to reduce the risk of transmission. Pyrethroids are used extensively for adult mosquito control, especially during dengue outbreaks. Pyrethroids promote activation and prolong the activation of the voltage gated sodium channel protein (VGSC) by interacting with two distinct pyrethroid receptor sites [1], formed by the interfaces of the transmembrane helix subunit 6 (S6) of domains II and III. Mutations of S6 in domains II and III synergize so that double mutants have higher pyrethroid resistance than mutants in either domain alone. Computer models predict an allosteric interaction between mutations in the two domains. In Ae. aegypti, a Ile1,016 mutation in the S6 of domain II was discovered in 2006 and found to be associated with pyrethroid resistance in field populations in Mexico. In 2010 a second mutation, Cys1,534 in the S6 of domain III was discovered and also found to be associated with pyrethroid resistance and correlated with the frequency of Ile1,016. A linkage disequilibrium analysis was performed on Ile1,016 and Cys1,534 in Ae. aegypti collected in Mexico from 2000–2012 to test for statistical associations between S6 in domains II and III in natural populations. We estimated the frequency of the four dilocus haplotypes in 1,016 and 1,534: Val1,016/Phe1,534 (susceptible), Val1,016/Cys1,534, Ile1,016/Phe1,534, and Ile1,016/Cys1,534 (resistant). The susceptible Val1,016/Phe1,534 haplotype went from near fixation to extinction and the resistant Ile1,016/Cys1,534 haplotype increased in all collections from a frequency close to zero to frequencies ranging from 0.5–0.9. The Val1,016/Cys1,534 haplotype increased in all collections until 2008 after which it began to decline as Ile1,016/Cys1,534 increased. However, the Ile1,016/Phe1,534 haplotype was rarely detected; it reached a frequency of only 0.09 in one collection and subsequently declined. Pyrethroid resistance in the vgsc gene requires the sequential evolution of two mutations. The Ile1,016/Phe1,534 haplotype appears to have low fitness suggesting that Ile1,016 was unlikely to have evolved independently. Instead the Cys1,534 mutation evolved first but conferred only a low level of resistance. Ile1,016 in S6 of domain II then arose from the Val1,016/Cys1,534 haplotype and was rapidly selected because double mutants confer higher pyrethroid resistance. This pattern suggests that knowledge of the frequencies of mutations in both S6 in domains II and III are important to predict the potential of a population to evolve kdr. Susceptible populations with high Val1,016/Cys1,534 frequencies are at high risk for kdr evolution, whereas susceptible populations without either mutation are less likely to evolve high levels of kdr, at least over a 10 year period. Constant use of pyrethroid insecticides has driven mosquito populations to develop resistance. In Aedes aegypti, the primary mosquito vector of dengue, yellow Fever, and chikungunya viruses, pyrethroid resistance is primarily associated with mutations in the voltage-gated sodium channel protein. One mutation occurs in codon 1,016 and involves a replacement of valine with isoleucine (Ile1, 016), and a second located in subunit 6 of domain III in codon 1,534, replaces phenylalanine with cysteine (Cys1,534). In Mexico, we found that Cys1,534 was present in the same mosquito collections that were previously analyzed for Ile1,016. In this study, we performed a linkage disequilibrium analysis on both Ile1,016 and Cys1,534 in Mexican collections from 2000–2012. Our analysis suggests that pyrethroid resistance requires the sequential evolution of the two mutations and that Cys1,534 must occur first and appears to enable the Ile1,016 mutation to survive.

Aedes aegypti, the ‘yellow fever mosquito’, is the primary vector to humans of dengue and yellow fever flaviviruses (DENV, YFV), and is a known vector of the chikungunya alphavirus (CV). Because vaccines are not yet available for DENV or CV or are inadequately distributed in developing countries (YFV), management of Ae. aegypti remains the primary option to prevent and control outbreaks of the diseases caused by these arboviruses. Permethrin is one of the most widely used active ingredients in insecticides for suppression of adult Ae. aegypti. In 2007, we documented a replacement mutation in codon 1,016 of the voltage-gated sodium channel gene (para) of Ae. aegypti that encodes an isoleucine rather than a valine and confers resistance to permethrin. Ile1,016 segregates as a recessive allele conferring knockdown resistance to homozygous mosquitoes at 5–10 µg of permethrin in bottle bioassays. A total of 81 field collections containing 3,951 Ae. aegypti were made throughout México from 1996 to 2009. These mosquitoes were analyzed for the frequency of the Ile1,016 mutation using a melting-curve PCR assay. Dramatic increases in frequencies of Ile1,016 were recorded from the late 1990's to 2006–2009 in several states including Nuevo León in the north, Veracruz on the central Atlantic coast, and Yucatán, Quintana Roo and Chiapas in the south. From 1996 to 2000, the overall frequency of Ile1,016 was 0.04% (95% confidence interval (CI95) = 0.12%; n = 1,359 mosquitoes examined). The earliest detection of Ile1,016 was in Nuevo Laredo on the U.S. border in 1997. By 2003–2004 the overall frequency of Ile1,016 had increased ∼100-fold to 2.7% (±0.80% CI95; n = 808). When checked again in 2006, the frequency had increased slightly to 3.9% (±1.15% CI95; n = 473). This was followed in 2007–2009 by a sudden jump in Ile1,016 frequency to 33.2% (±1.99% CI95; n = 1,074 mosquitoes). There was spatial heterogeneity in Ile1,016 frequencies among 2007–2008 collections, which ranged from 45.7% (±2.00% CI95) in the state of Veracruz to 51.2% (±4.36% CI95) in the Yucatán peninsula and 14.5% (±2.23% CI95) in and around Tapachula in the state of Chiapas. Spatial heterogeneity was also evident at smaller geographic scales. For example within the city of Chetumal, Quintana Roo, Ile1,016 frequencies varied from 38.3%–88.3%. A linear regression analysis based on seven collections from 2007 revealed that the frequency of Ile1,016 homozygotes accurately predicted knockdown rate for mosquitoes exposed to permethrin in a bioassay (R² = 0.98). We have recorded a dramatic increase in the frequency of the Ile1,016 mutation in the voltage-gated sodium channel gene of Ae. aegypti in México from 1996 to 2009. This may be related to heavy use of permethrin-based insecticides in mosquito control programs. Spatial heterogeneity in Ile1,016 frequencies in 2007 and 2008 collections may reflect differences in selection pressure or in the initial frequency of Ile1,016. The rapid recent increase in Ile1,016 is predicted by a simple model of positive directional selection on a recessive allele. Unfortunately this model also predicts rapid fixation of Ile1,016 unless there is negative fitness associated with Ile1,016 in the absence of permethrin. If so, then spatial refugia of susceptible Ae. aegypti or rotational schedules of different classes of adulticides could be established to slow or prevent fixation of Ile1,016. Pyrethroid insecticides prolong the opening of voltage-dependent sodium channels in insect nerves to produce instant paralysis and “knock-down.” Many insects have evolved knock-down resistance through nonsynonymous mutations that reduce pyrethroid binding in the channels. In 2006 we discovered one such mutation in the arbovirus mosquito vector Aedes aegypti, called Ile1,016, that confers very high knockdown resistance to the pyrethroid insecticide permethrin in mosquitoes homozygous for this mutation. We examined collections of Ae. aegypti from México during 1996–2009 and found that the overall Ile1,016 frequency increased from <0.1% in 1996–2000, to 2%–5% in 2003–2006, to 38.3%–88.3% in 2007–2009 depending upon collection location. We also demonstrate a strong linear relationship between the frequency of Ile1,016 homozygotes and knockdown rate in bioassays and speculate that widespread use of permethrin-based insecticides in México may be impacting the frequency of Ile1,016. Such a rapid increase is predicted by a simple model of positive directional selection acting on a recessive allele. Unfortunately this model also predicts rapid fixation of Ile1,016 unless there is a negative fitness associated with Ile1,016 in the absence of permethrin and if insecticidal pressure can be reduced.

Changes in gene expression before, during and after five generations of permethrin laboratory selection were monitored in six strains of Aedes aegypti : five F ₂ – F ₃ collections from the Yucatán Peninsula of Mexico and one F ₂ from Iquitos, Peru. Three biological replicate lines were generated for each strain. The response to selection was measured as changes in the lethal and knockdown permethrin concentrations (LC ₅₀ , KC ₅₀ ) and in the frequency of the Ile1,016 substitution in the voltage‐gated sodium channel ( para ) gene. Changes in expression of 290 metabolic detoxification genes were measured using the ‘ Aedes Detox ’ microarray. Selection simultaneously increased the LC ₅₀ , KC ₅₀ and Ile1,016 frequency. There was an inverse relationship between Ile1,016 frequency and the numbers of differentially transcribed genes. The Iquitos strain lacked the Ile1,016 allele and 51 genes were differentially transcribed after selection as compared with 10–18 genes in the Mexican strains. Very few of the same genes were differentially transcribed among field strains but 10 cytochrome P ₄₅₀ genes were upregulated in more than one strain. Laboratory adaptation to permethrin in Ae. aegypti is genetically complex and largely conditioned by geographic origin and pre‐existing target site insensitivity in the para gene. The lack of uniformity in the genes that responded to artificial selection as well as differences in the direction of their responses challenges the assumption that one or a few genes control permethrin metabolic resistance. Attempts to identify one or a few metabolic genes that are predictably associated with permethrin adaptation may be futile.

The mosquito Aedes aegypti is the principal vector of dengue and yellow fever flaviviruses. Permethrin is an insecticide used to suppress Ae. aegypti adult populations but metabolic and target site resistance to pyrethroids has evolved in many locations worldwide. Quantitative trait loci (QTL) controlling permethrin survival in Ae. aegypti were mapped in an F3 advanced intercross line. Parents came from a collection of mosquitoes from Isla Mujeres, México, that had been selected for permethrin resistance for two generations and a reference permethrin-susceptible strain originally from New Orleans. Following a 1-hr permethrin exposure, 439 F3 adult mosquitoes were phenotyped as knockdown resistant, knocked down/recovered, or dead. For QTL mapping, single nucleotide polymorphisms (SNPs) were identified at 22 loci with potential antixenobiotic activity including genes encoding cytochrome P450s (CYP), esterases (EST), or glutathione transferases (GST) and at 12 previously mapped loci. Seven antixenobiotic genes mapped to chromosome I, six to chromosome II, and nine to chromosome III. Two QTL of major effect were detected on chromosome III. One corresponds with a SNP previously associated with permethrin resistance in the para sodium channel gene and the second with the CCEunk7o esterase marker. Additional QTL but of relatively minor effect were also found. These included two sex-linked QTL on chromosome I affecting knockdown and recovery and a QTL affecting survival and recovery. On chromosome II, one QTL affecting survival and a second affecting recovery were detected. The patterns confirm that mutations in the para gene cause target-site insensitivity and are the major source of permethrin resistance but that other genes dispersed throughout the genome contribute to recovery and survival of mosquitoes following permethrin exposure.