Karla Saavedra-Rodriguez

A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti

Pyrethroids are commonly used as mosquito adulticides and evolution of resistance to these compounds is a major threat to public health. ‘Knockdown resistance’ to pyrethroids (kdr) is frequently caused by nonsynonymous mutations in the voltage‐gated sodium channel transmembrane protein (para) that reduce pyrethroid binding. Early detection of kdr is critical to the development of resistance management strategies in mosquitoes including Aedes aegypti, the most prevalent vector of dengue and yellow fever viruses. Brengues et al. described seven novel mutations in hydrophobic segment 6 of domain II of para in Ae. aegypti. Assays on larvae from strains bearing these mutations indicated reduced nerve sensitivity to permethrin inhibition. Two of these occurred in codons Iso1011 and Val1016 in exons 20 and 21 respectively. A transition in the third position of Iso1011 encoded a Met1011 replacement and a transversion in the second position of Val1016 encoded a Gly1016 replacement. We have screened this same region in 1318 mosquitoes in 32 additional strains; 30 from throughout Latin America. While the Gly1016 allele was never detected in Latin America, we found two new mutations in these same codons. A transition in the first position of codon 1011 encodes a Val replacement while a transition in the first position of codon 1016 encodes an Iso replacement. We developed PCR assays for these four mutations that can be read either on an agarose gel or as a melting curve. Selection experiments, one with deltamethrin on a field strain from Santiago de Cuba and another with permethrin on a strain from Isla Mujeres, Mexico rapidly increased the frequency of the Iso1016 allele. Bioassays of F3 offspring arising from permethrin susceptible Val1016 homozygous parents and permethrin resistant Iso1016 homozygous parents show that Iso1016 segregates as a recessive allele in conferring kdr. Analysis of segregation between alleles at the 1011 and 1016 codons in the F3 showed a high rate of recombination even though the two codons are only separated by a ~250 bp intron. The tools and information presented provide a means for early detection and characterization of kdr that is critical to the development of strategies for resistance management.

Recent Rapid Rise of a Permethrin Knock Down Resistance Allele in Aedes aegypti in México

Background 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. Methods and Findings 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 1990s 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 (R2 = 0.98). Conclusions 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.

Quantitative Trait Loci Mapping of Genome Regions Controlling Permethrin Resistance in the Mosquito Aedes aegypti

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.

Transcription of detoxification genes after permethrin selection in the mosquito Aedes aegypti.

Changes in gene expression before, during and after five generations of permethrin laboratory selection were monitored in six strains of Aedes aegypti: five F2–F3 collections from the Yucatán Peninsula of Mexico and one F2 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 (LC50, KC50) 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 LC50, KC50 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 P450 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.

Coevolution of the Ile1,016 and Cys1,534 Mutations in the Voltage Gated Sodium Channel Gene of Aedes aegypti in Mexico.

Background 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. Methodology/Principal Findings 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. Conclusion/Significance 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.

Rhythms and synchronization patterns in gene expression in the Aedes aegypti mosquito

BackgroundAedes aegypti is arguably the most studied of all mosquito species in the laboratory and is the primary vector of both Dengue and Yellow Fever flaviviruses in the field. A large number of transcriptional studies have been made in the species and these usually report transcript quantities observed at a certain age or stage of development. However, circadian oscillation is an important characteristic of gene expression in many animals and plants, modulating both their physiology and behavior. Circadian gene expression in mosquito species has been previously reported but for only a few genes directly involved in the function of the molecular clock.ResultsHerein we analyze the transcription profiles of 21,494 messenger RNAs using an Ae. aegypti Agilent® microarray. Transcripts were quantified in adult female heads at 24 hours and then again at 72 hours and eight subsequent time points spaced four hours apart. We document circadian rhythms in multiple molecular pathways essential for growth, development, immune response, detoxification/pesticide resistance. Circadian rhythms were also noted in ribosomal protein genes used for normalization in reverse transcribed PCR (RT-PCR) to determine transcript abundance. We report pervasive oscillations and intricate synchronization patterns relevant to all known biological pathways.ConclusionThese results argue strongly that transcriptional analyses either need to be made over time periods rather than confining analyses to a single time point or development stage or exceptional care needs to be made to synchronize all mosquitoes to be analyzed and compared among treatment groups.

Differential transcription profiles in Aedes aegypti detoxification genes after temephos selection.

The mosquito Aedes aegypti is the main vector of Dengue and Yellow Fever flaviviruses. The organophosphate insecticide temephos is a larvicide that is used globally to control Ae. aegypti populations; many of which have in turn evolved resistance. Target site alteration in the acetylcholine esterase of this species has not being identified. Instead, we tracked changes in transcription of metabolic detoxification genes using the Ae. aegypti ‘Detox Chip’ microarray during five generations of temephos selection. We selected for temephos resistance in three replicates in each of six collections, five from Mexico, and one from Peru. The response to selection was tracked in terms of lethal concentrations. Uniform upregulation was seen in the epsilon class glutathione‐S‐transferase (eGST) genes in strains from Mexico prior to laboratory selection, while eGSTs in the Iquitos Peru strain became upregulated after five generations of temephos selection. While expression of many carboxyl/cholinesterase esterase (CCE) genes increased with selection, no single esterase was consistently upregulated and this same pattern was noted in the cytochrome P450 monooxygenase (CYP) genes and in other genes involved in reduction or oxidation of xenobiotics. Bioassays using glutathione‐S‐transferase (GST), CCE and CYP inhibitors suggest that various CCEs instead of GSTs are the main metabolic mechanism conferring resistance to temephos. We show that temephos‐selected strains show no cross resistance to permethrin and that genes associated with temephos selection are largely independent of those selected with permethrin in a previous study.

Frequency of V1016I and F1534C mutations in the voltage-gated sodium channel gene in Aedes aegypti in Venezuela.

BACKGROUND The V1016I and F1534C mutations in the voltage-gated sodium channel gene have been associated with resistance to pyrethroids and DDT in Aedes aegypti mosquitoes. A study was carried out to determine the frequency of I1016 and C1534 by real-time PCR in five natural populations of Ae. aegypti in Venezuela during 2008, 2010 and 2012, as well as in a strain selected with 0.14 µg of deltamethrin for 15 generations. RESULTS In natural populations, frequencies of I1016 varied between 0.01 and 0.37, and frequencies of C1534 between 0.35 and 1.0. In the Pampanito strain, the frequency of I1016 increased from 0.02 in F1 up to 0.5 in F15 and from 0.35 up to fixation for C1534 after selection with deltamethrin. CONCLUSION The results showed that C1534 frequencies are higher than I1016 frequencies in natural populations of Ae. aegypti in Venezuela, and that deltamethrin selected the C1534 more rapidly than I1016.

Differential transcription profiles inAedes aegyptidetoxification genes after temephos selection

The mosquito Aedes aegypti is the main vector of Dengue and Yellow Fever flaviviruses. The organophosphate insecticide temephos is a larvicide that is used globally to control Ae. aegypti populations; many of which have in turn evolved resistance. Target site alteration in the acetylcholine esterase of this species has not being identified. Instead, we tracked changes in transcription of metabolic detoxification genes using the Ae. aegypti ‘Detox Chip’ microarray during five generations of temephos selection. We selected for temephos resistance in three replicates in each of six collections, five from Mexico, and one from Peru. The response to selection was tracked in terms of lethal concentrations. Uniform upregulation was seen in the epsilon class glutathione‐S‐transferase (eGST) genes in strains from Mexico prior to laboratory selection, while eGSTs in the Iquitos Peru strain became upregulated after five generations of temephos selection. While expression of many carboxyl/cholinesterase esterase (CCE) genes increased with selection, no single esterase was consistently upregulated and this same pattern was noted in the cytochrome P450 monooxygenase (CYP) genes and in other genes involved in reduction or oxidation of xenobiotics. Bioassays using glutathione‐S‐transferase (GST), CCE and CYP inhibitors suggest that various CCEs instead of GSTs are the main metabolic mechanism conferring resistance to temephos. We show that temephos‐selected strains show no cross resistance to permethrin and that genes associated with temephos selection are largely independent of those selected with permethrin in a previous study.

Towards a Casa Segura: A Consumer Product Study of the Effect of Insecticide-Treated Curtains on Aedes aegypti and Dengue Virus Infections in the Home

The home, or domicile, is the principal environment for transmission of dengue virus (DENV) between humans and mosquito vectors. Community-wide distribution of insecticide-treated curtains (ITCs), mimicking vector control program-driven interventions, has shown promise to reduce DENV infections. We conducted a Casa Segura consumer product intervention study in Mérida, Mexico to determine the potential to reduce intradomicillary DENV transmission through ITC use in individual homes. Dengue virus infections in mosquitoes and in humans were reduced in homes with ITCs in one of two study subareas. Overall, ITCs reduced intradomicillary DENV transmission; ITC homes were significantly less likely to experience multiple DENV infections in humans than NTC homes. Dengue virus-infected Aedes aegypti females were reduced within the ITC homes where curtain use was highest. Some homes yielded up to nine infected Ae. aegypti females. This study provides insights regarding best practices for Casa Segura interventions to protect homes from intradomicillary DENV transmission.