Skip to main content
Advertisement
  • Loading metrics

Quantification of permethrin resistance and kdr alleles in Florida strains of Aedes aegypti (L.) and Aedes albopictus (Skuse)

  • Alden S. Estep ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Project administration, Supervision, Writing – original draft

    alden.estep@ars.usda.gov

    Affiliation CMAVE Detachment, Navy Entomology Center of Excellence, Gainesville, FL, United States of America

  • Neil D. Sanscrainte,

    Roles Data curation, Formal analysis, Project administration, Writing – original draft

    Affiliation Mosquito and Fly Research Unit, United States Department of Agriculture, Agricultural Research Service, Center for Medical, Agricultural, and Veterinary Entomology, Gainesville, FL, United States of America

  • Christy M. Waits,

    Roles Data curation, Investigation, Visualization, Writing – original draft

    Affiliation CMAVE Detachment, Navy Entomology Center of Excellence, Gainesville, FL, United States of America

  • Sarah J. Bernard,

    Roles Data curation, Visualization, Writing – original draft

    Affiliation CMAVE Detachment, Navy Entomology Center of Excellence, Gainesville, FL, United States of America

  • Aaron M. Lloyd,

    Roles Investigation, Resources, Writing – original draft

    Affiliation Pasco County Mosquito Control District, Odessa, FL, United States of America

  • Keira J. Lucas,

    Roles Formal analysis, Investigation, Resources, Writing – original draft

    Affiliation Collier Mosquito Control District, Naples, FL, United States of America

  • Eva A. Buckner,

    Roles Investigation, Resources, Writing – original draft

    Current address: Florida Medical Entomology Laboratory, Vero Beach, FL, United States of America

    Affiliation Manatee Mosquito Control District, Palmetto, FL, United States of America

  • Rajeev Vaidyanathan,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Clarke Inc., Saint Charles, IL, United States of America

  • Rachel Morreale,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Lee County Mosquito Control, Lehigh Acres, FL, United States of America

  • Lisa A. Conti,

    Roles Project administration, Writing – review & editing

    Affiliation Florida Department of Agriculture and Consumer Services, Tallahassee, FL, United States of America

  • James J. Becnel

    Roles Formal analysis, Funding acquisition, Supervision, Writing – original draft

    Affiliation Mosquito and Fly Research Unit, United States Department of Agriculture, Agricultural Research Service, Center for Medical, Agricultural, and Veterinary Entomology, Gainesville, FL, United States of America

Abstract

Recent outbreaks of locally transmitted dengue and Zika viruses in Florida have placed more emphasis on integrated vector management plans for Aedes aegypti (L.) and Aedes albopictus Skuse. Adulticiding, primarily with pyrethroids, is often employed for the immediate control of potentially arbovirus-infected mosquitoes during outbreak situations. While pyrethroid resistance is common in Ae. aegypti worldwide and testing is recommended by CDC and WHO, resistance to this class of products has not been widely examined or quantified in Florida. To address this information gap, we performed the first study to quantify both pyrethroid resistance and genetic markers of pyrethroid resistance in Ae. aegypti and Ae. albopictus strains in Florida. Using direct topical application to measure intrinsic toxicity, we examined 21 Ae. aegypti strains from 9 counties and found permethrin resistance (resistance ratio (RR) = 6-61-fold) in all strains when compared to the susceptible ORL1952 control strain. Permethrin resistance in five strains of Ae. albopictus was very low (RR<1.6) even when collected from the same containers producing resistant Ae. aegypti. Characterization of two sodium channel kdr alleles associated with pyrethroid-resistance showed widespread distribution in 62 strains of Ae. aegypti. The 1534 phenylalanine to cysteine (F1534C) single nucleotide polymorphism SNP was fixed or nearly fixed in all strains regardless of RR. We observed much more variation in the 1016 valine to isoleucine (V1016I) allele and observed that an increasing frequency of the homozygous V1016I allele correlates strongly with increased RR (Pearson corr = 0.905). In agreement with previous studies, we observed a very low frequency of three kdr genotypes, IIFF, VIFF, and IIFC. In this study, we provide a statewide examination of pyrethroid resistance, and demonstrate that permethrin resistance and the genetic markers for resistance are widely present in FL Ae. aegypti. Resistance testing should be included in an effective management program.

Author summary

Aedes aegypti (Yellow-fever mosquito) and Aedes albopictus (Asian Tiger mosquito) can vector a variety of arboviruses that cause diseases and are thus a public health concern. Pyrethroid insecticide resistance is common in Ae. aegypti in many locations worldwide and can adversely affect vector control operations. However, the resistance status of these vectors in Florida is largely unreported and recent local transmission of dengue and Zika viruses has made this information critical for effective control operations. In this study, we showed that permethrin resistance and two common SNPs of the voltage gated sodium channel (V1016I and F1534C) previously associated with pyrethroid resistance were widely present in Florida Ae. aegypti strains. We also observed a strong correlation between the dilocus knock down response (kdr) genotype and resistance ratio (RR) as determined by topical application, which suggests, as have others, that kdr frequency may be a useful indicator of resistance in Aedes aegypti.

Introduction

Local vector control programs play a part in public health. In many countries, these programs serve as the primary defense against the spread of several mosquito-borne diseases. Effective Integrated Vector Management (IVM) programs rely on surveillance information coupled with multiple vector control strategies such as chemical adulticiding to reduce vector populations and arbovirus transmission. Limited recent transmission of locally acquired dengue and Zika viruses in the southeastern continental US, primarily Florida, has brought renewed attention to the importance of IVM programs where the potential vectors, Aedes aegypti (L.) and Ae. albopictus Skuse, have long been present. However, for IVM programs in Florida to effectively control Aedes vectors and to reduce dengue and Zika virus transmission during an outbreak, it is essential to know which adulticide products are effective against local Ae. aegypti and Ae. albopictus strains.

Organophosphates and pyrethroids are the only two classes of insecticides available to public health agencies for control of adult mosquitoes in the US. Compared to organophosphate insecticides, pyrethroids have higher public acceptance, rapid knockdown, relatively low costs, and are generally the product class of choice when adulticiding is required [1]. Unfortunately, years of pyrethroid insecticide use and previous DDT usage has increased the frequency of genetic resistance and enhanced enzymatic detoxification in insects like mosquitoes. Genetic target site changes of the sodium channel, known as knockdown resistance (kdr) mutations, are a relatively common insect response to selective pressure by pyrethroids [2,3]. Although a variety of other single nucleotide polymorphisms (SNPs) have been noted in Ae. aegypti, the primary two SNPs assessed to determine pyrethroid resistance are at codons 1016 and 1534 (positions according to standard M. domestica notation) [4,5]. One allele of the 1016 mutation is geographically distinct to the Western hemisphere and results in the replacement of the normal valine with an isoleucine (V1016I), while the 1534 mutation results in a phenylalanine to cysteine change (F1534C). There is currently debate about the individual toxicological effect of these two mutations, but they are consistently present in resistant strains [4]. Heterozygous and homozygous combinations of these two SNPs could result in nine possible genotypes, but strong linkage between the SNPs has been noted with only six of the genotypes observed in a study of Mexican Ae. aegypti [6].

More recent work indicates that the kdr genotype in Ae. aegypti may be more complicated than the combination of 1016I and 1534C SNPs. An additional SNP at position 410 has been identified and this valine to leucine (V410L) change, in addition to the other two SNPs, may be the strongly resistant phenotype [7]. The noted strong linkage between 410L and 1016I indicates that assessment of 1016 may still be adequate to assess the strongly resistant phenotype [7]. As expected based on the dilocus genotype distribution, several of the trilocus genotypes were also not frequently present in the wild [8].

Pyrethroid resistance and the distribution of kdr alleles have been well documented in Ae. aegypti strains from Central America, South America and the Caribbean [913]. Recent testing as part of the Zika emergency response in Puerto Rico has shown that isolated, early reports of resistance were indicative of widespread resistance on the island [1316]. Studies have also shown the pyrethroid resistance is not specific to a particular pyrethroid but is generally class-wide to type I, type II, and non-ester pyrethroids [13, 17, 18]. Kdr alleles and pyrethroid resistance are widely distributed in Mexico, including Nuevo Laredo which lies just south of the US-Mexico border [11].

In contrast, little has been published about pyrethroid resistance in Aedes strains within the continental US. In the Garcia et al. [11] study mentioned above, the authors did not find kdr alleles in a Houston, TX collection from 1999. Recent reviews of the Aedes resistance literature listed no reports of resistance in continental US Ae. aegypti [4, 19], but Cornel et al. [20] recently demonstrated toxicological resistance and sodium channel mutations in invasive populations of Ae. aegypti in California. Two recent studies using CDC bottle bioassays do indicate resistance in strains from the southern US, but these studies do not provide any quantification of the strength of the resistance nor did they examine the presence of kdr alleles [21, 22]. In contrast, Ae. albopictus does not appear to frequently develop kdr resistance, and most US strains tested thus far have only shown very minimal pyrethroid resistance [4, 23, 24, 25].

Resistance surveillance is recommended by the CDC and statewide initiatives to map pyrethroid resistance have begun in Florida and California [22,26]. Resistance information is critically important for operational decision making as part of an effective IVM program. In this study, a collaborative group of government, academic, industry and vector control district stakeholders collected Ae. aegypti and Ae. albopictus adults, eggs or larvae from more than 200 locations throughout Florida to assess the extent and intensity of pyrethroid resistance. The goal of this study was to improve vector control operations by producing a resistance map and apply this information to make more effective control decisions. We determined permethrin resistance ratios (RRs) relative to the susceptible ORL1952 strain for 21 wild type strains of Ae. aegypti and 5 strains of Ae. albopictus using direct topical application, the WHO gold standard assay for determination of intrinsic toxicity [27]. While the CDC bottle bioassay is a more commonly used assay, it is a threshold assay and does not quantitate resistance levels which was one of our primary goals. Nearly 5,000 Ae. aegypti from numerous locations were genotyped by allele-specific PCR to assess the frequency of V1016I and F1534C alleles and rapidly visualize the pattern of resistance throughout the state of Florida.

Methods

Mosquito strains and rearing

The toxicological profiles and rearing procedures for the susceptible (ORL1952) and resistant (Puerto Rico–BEIResources, NR-48830) strains of Ae. aegypti used in this study have been described previously [17]. Strain information, including specific location information, collectors, and dates of collection are noted in Table 1. Field strains for toxicology testing were often collected as mixed (Ae. albopictus and Ae. aegypti) eggs laid on seed germination paper (Anchor Paper, St. Paul, MN) placed in a variety of oviposition containers including black plastic stadium cups, plastic cemetery vases, and glass containers. Field collected eggs were hatched by soaking papers for 48–72 hours at 27 oC in rearing trays of untreated well water or deionized water (diH2O). Papers were carefully removed, and larvae were reared through adult emergence with a reduced feeding regimen compared to the standard rearing protocol for the ORL1952 susceptible strain [28] due to sensitivity to overfeeding. Strains from Monroe, Seminole, Orange, Hernando and Sarasota counties were collected as larvae from a variety of manmade and natural containers (tires, plant pots, bottles, buckets, etc.), rinsed with diH2O and then placed into the standard larval rearing procedure described above. Pupae were collected from rearing trays and placed in emergence chambers or 12” x 12” screen cages (BioQuip Models 1425 & 1450B, Rancho Dominquez, CA). After emergence, wildtype mosquitoes were briefly chilled to 4 oC and sorted to species. Strains for toxicology testing were produced from locations that had more than 30 wildtype founders. If multiple nearby locations were combined to produce a strain, the GPS location in Table 1 represents the GPS centroid of the sites that contributed the mosquitoes. Individual oviposition cup locations and life stages tested are provided in S1 File. Eggs were produced using standard rearing methods [17]. Colony strains were provided 2–10 bloodmeals on weekly intervals to collect F1 eggs. Warmed bovine blood was provided to produce 2,000 or more eggs per strain. If feeding was poor, warmed blood was spiked with 1 mM ATP as a phagostimulant. Bovine blood for mosquito feeding was purchased by the USDA under contract from a local, licensed abattoir. Blood was collected during normal operations of the abattoir from the waste stream after animal slaughter. Under CFR9, Parts 1–3, tissues, including blood, collected from dead livestock intended as food are exempt from IACUC regulation.

thumbnail
Table 1. Collection information for Florida Aedes aegypti and Ae. albopictus used in this study.

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

Aedes aegypti strains from St. Augustine, Clearwater, and Vero Beach, FL were provided as F1 eggs produced at Anastasia Mosquito Control District, Pinellas County Mosquito Control, and the Florida Medical Entomology Laboratory, respectively.

To produce mosquitoes for bioassay testing, eggs from the F1 or F2 generation of field collected strains, the lab susceptible strains (ORL1952) and the pyrethroid-resistant strain were hatched and reared as described above in covered trays at a density of approximately 1,000 mosquito larvae/tray. Adult mosquitoes emerged into 12” x 12” screen cages and provided with cotton saturated with 10% sucrose in diH2O. Females used for bioassays were 3–7 days post-emergence.

Determination of resistance ratio by topical application

The adult topical bioassay has been described previously in detail [28, 29]. For these studies, the permethrin stock solution in DMSO (Product #N-12848-250MG, Chemservice, Westchester, PA) and all dilutions in acetone were prepared gravimetrically. An initial 10-fold dilution series was prepared over a range of relevant concentrations [17]. Sub-dilutions were prepared from the 10-fold dilution series as necessary to determine the critical region of the dose curve. Three assays (n = 3) were performed for all strains with listed LD50s unless limited numbers of F1 test mosquitoes allowed only two replicates. Average mass/female was calculated for each strain by weighing a cohort of 50–100 females before each replicate.

LD50s, standard errors, and goodness of fit were determined from dose-mortality curves using SigmaPlot v13 (Systat software Inc., San Jose, CA) with data fit to a four-parameter logistic model. To provide a comparative metric between strains that may have different body sizes, doses applied were divided by the average mass of the mosquitoes before curve fitting. This results in an LD50 of ng of active ingredient per mg of mosquito. Resistance ratios for strains were calculated using the LD50 of the field strain divided by the LD50 of the susceptible ORL1952 strain included in the same assay. In this study, we use the WHO scale to define the levels of resistance [30]. When RR is less than 5, the field population is considered susceptible. When the RR is between 5 and10, the mosquitoes are considered to have moderate resistance. A RR greater than 10 indicates that mosquitoes are highly resistant [31].

Strain kdr genotyping

Genotypes for individual mosquitoes or eggs were determined using melt curve analysis with previously described allele-specific primers for the 1016 and 1534 SNPs [32, 33]. Preliminary assays were conducted using eggs and adults from the same generation to verify that the resulting frequencies were similar (within 10%) between life stages rather than showing agenotype bias due to varied hatch rates or differing larval mortality (S2 File). Assays were performed in 96 well plates on a StepOnePlus (Applied Biosystems) or QuantStudio5 (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA) using SYBR green chemistry. Plates were loaded with 8 μl of PCR master mix containing SYBR Select Master Mix (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA), nuclease free water (NFW), and three primers (Table 2). Due to differences in efficiency from the added GC tails, initial primer titering was necessary to ensure that the melt curves of homozygous controls would be accurately called by the analysis software. Individual adult females were homogenized for 60 seconds at max speed in 100 μl of NFW on a bead beater (BioSpec, Bartlesville, OK). Individual eggs, randomly sampled from all available egg papers of the same generation, were similarly treated but homogenized in 50 μl of NFW. Immediately after homogenization, samples were centrifuged at 10,000 relative centrifugal force (rcf) for 60 seconds to pellet solids. Two microliters of each supernatant were added to 8 μl of PCR mix for each primer set and then subjected to standard cycling conditions (3 min @ 95 °C; 40 cycles @ 95 °C for 3 sec, 60 °C for 15 sec). Melt curve analysis followed cycling with acquisition of fluorescence data every 0.3 °C as the temperature was ramped from 60 °C to 95 °C. Characteristic melting temperature (Tm) peaks in the derivative fluorescence data indicate the presence of specific alleles [32, 33]. For the 1016 mutation, a codon for the susceptible valine has a Tm peak at 86±0.3 °C while the isoleucine codon has a Tm of 77.3±0.3 °C. The 1534 phenylalanine has a Tm of 79.8±0.3 °C while the mutant cysteine produces a peak at 84.7±0.3 °C. Homozygotes for either allele produce one peak while heterozygotes produce peaks at both Tms. All assay plates included two susceptible ORL1952 strain and two resistant PR strain samples as negative and positive controls, respectively. Most plates also contained artificial heterozygotes created by including ORL1952 and PR homogenates in the same sample. Well positions of individual mosquitoes were maintained in both plates (one plate for each locus) to allow genotyping of an individual for both SNPs. Frequencies for each of the nine genotypes (VVFF, VVFC, VVCC, VIFF, VIFC, VICC, IIFF, IIFC, IICC) were calculated by dividing the specific genotype by the total number tested from each area.

Mapping

Base maps were created using ArcGIS software by Esri. ArcGIS® is the intellectual property of ESRI and is used herein under license. Copyright © ESRI. GIS data sources were ESRI and Tele Atlas. All rights reserved. (For more information about Esri® software, please visit www.esri.com). Permission to publish this content was verified from the ESRI Redistribution Rights Matrix at https://www.esri.com/~/media/Files/Pdfs/legal/pdfs/redist_rights_106.pdf?la=en

Maps were exported to GIMP 2.8 and additional layers with the kdr or resistance ratio data was added. Pie chart representations of kdr allele frequencies were created with Microsoft Excel and exported to layers added to the basemaps.

Results

Topical bioassay

All Florida strains of Ae. aegypti were resistant to permethrin when compared to the ORL1952 strain, which has been in a continuous laboratory colony for nearly seventy years (Fig 1 and Table 3). The field strains showed varied levels of resistance, from 6-fold to 61-fold compared to the ORL1952 strain. The two least resistant Ae. aegypti strains were collected from flowerpots on Big Coppitt Key (Monroe County) and from oviposition cups in Cortez (Manatee County) at 6.0 and 6.8-fold, respectively. Ae. aegypti mosquitoes collected from a tire facility in Orlando (Orange County) were about 17.2-fold more resistant than Ae. aegypti collected from the same area generations earlier and used to produce the ORL1952 strain. While most of the FL strains were 15 to 35-fold resistant compared to the lab strains, several strains with higher RR were identified. The strain from New Port Richey in Pasco County had the highest RR at 61.3-fold, which is like the RR of the Puerto Rico-resistant reference strain [17]. Strains from Miami Beach and Fort Myers had RRs above 50.

thumbnail
Fig 1. Permethrin resistance ratios of Ae. aegypti and Ae. albopictus (in blue) strains compared to the susceptible ORL1952 laboratory strain.

The Ae. albopictus RR for St John’s County (marked with asterisks) are taken from our previous study [24]. Base maps were sourced from ESRI and Tele Atlas data through ArcGIS Online under an enterprise license with USDA. Additional layers were added to the base map using GIMP 2.8.

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

thumbnail
Table 3. Permethrin resistance ratios of select Florida strains of Aedes aegypti and Ae. albopictus calculated from direct topical application.

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

The variability we observed in RR throughout the state was also seen at finer resolution within several, but not all, counties. In Miami-Dade County, the Miami Beach and East Wynwood strains were relatively resistant (57-fold and 33-fold), while nearby locations like South Wynwood and central Little River were much less resistant (Fig 1). Manatee County also had a range of RRs. The Anna Maria Island strain, collected from a densely populated barrier island, had a higher RR than nearby strains from Palmetto or Cortez (Fig 1 and Table 3). We also observed this same trend in Lee County but did not observe universal variation in RR.

Aedes albopictus competes with Ae. aegypti for oviposition sites, and in several locations Ae. albopictus eggs were collected in conjunction with Ae. aegypti eggs. Ae. albopictus strains were also subjected to topical application along with the Ae. aegypti from the same locations (Fig 1, in blue). In the Ae. albopictus strains we tested, only slight resistance to permethrin (<2-fold) was observed when compared to the ORL1952 strain. In Miami-Dade County, we observed large differences in RR between Ae. albopictus and Ae. aegypti even when collected from the same sites in areas such as south Wynwood and central Little River. St. Johns (St. Augustine), Orange (Orlando), and Lee (Alva) counties also had resistant Ae. aegypti and low RR Ae. albopictus.

Kdr allele distribution

Examination of kdr alleles in Ae. aegypti strains from 62 locations showed a range of genotypes (Figs 2 and 3). We observed that most strains were fixed or nearly fixed for the F1534C SNP. For many strains, more than 95% of the tested mosquitoes were homozygous for the 1534C (1534CC) and the remainder of the strain was made up of a few 1534 heterozygotes (1534FC). Only strains from Big Coppitt Key (13% FC), Cortez in Manatee County (38% FC), Central Wynwood (22% FC) and 2 strains from Little River (20% and 14% FC) had more than 10% of mosquitoes that were heterozygous at position 1534. Most strains had no mosquitoes without at least one copy of the 1534C allele but we did find two strains from central Wynwood (8% FF) and Cortez (11% FF) that still had appreciable numbers of susceptible alleles.

thumbnail
Fig 2. Frequency of kdr genotypes in 40 strains of FL Aedes aegypti.

Genotype frequencies were determined using the methods of [28, 29] as described in the methods section. Specific collection locations that are included in each tested population are noted in S1 File. Strains included a minimum of 25 individual organisms. Base maps were sourced from ESRI and Tele Atlas data through ArcGIS Online under an enterprise license with USDA. Additional layers were added to the base map using GIMP 2.8. Graphical representation of kdr frequencies was produced using Microsoft Excel and included as an additional layer.

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

thumbnail
Fig 3. Kdr genotype frequencies in 20 populations of Aedes aegypti from Miami-Dade County, FL.

Genotype frequencies were determined using the methods of [32, 33] as described in the methods section. Specific collection locations that are included in each tested population are noted in S1 File. All data is based on a minimum of 76 tested individuals. Base maps were sourced from ESRI and Tele Atlas data through ArcGIS Online under an enterprise license with USDA. Additional layers were added to the base map using GIMP 2.8. Graphical representation of kdr frequencies was produced using Microsoft Excel and included as an additional layer.

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

There was much more variation throughout the state at position 1016. Strains from Big Coppitt Key, Longboat Key, Orlando and Clearwater had the lowest percentages of the homozygous 1016II at 11.0, 14.1, 14.9, and 17.8%, respectively. We did observe strains with high levels of 1016II. Six of eight strains examined from Pasco County had 1016II frequencies above 75%, including two strains from New Port Richey at 100% 1016II. Due to the 2016 Zika outbreak and resulting massive public health response, we heavily sampled strains from southeast Florida for allele frequencies. Two strains from Broward and several from Miami-Dade County had 1016II frequencies above 75% including a Miami Beach strain at 91% 1016II (Fig 3). Select strains from Lee and Collier counties in southwest Florida were also above 75% 1016II (Fig 2).

As with RR, in several counties we observed large differences in allele frequency from one part of the county to another. In Lee County (Fig 2), an inland strain from Alva had a relatively low level 1016II frequency (20%) while four strains from near the Gulf Coast all had 1016II frequencies of greater than 55%. This variation was also observed at much closer scale between strains in neighboring cities. The Manatee County Anna Maria Island strain had a much higher 1016II allele frequency than nearby strains from Palmetto or Longboat Key (Fig 2). While overall IICC frequencies were high in Pasco County, we did find variation at the neighborhood level in south Pasco County, where a strain with 75% 1016II was found just blocks from a strain with 30% 1016II (Fig 2). In Miami-Dade, we were able to test four strains from Wynwood and three from Little River collected during the 2016 Zika response (Fig 3). Again, we observed variations in the levels of 1016II within each neighborhood. North Wynwood was slightly more than 25% IICC while south Wynwood was nearly 50%. In Little River, the four strains ranged from 20 to 60%. Just across the water from Wynwood, the Miami Beach strain was greater than 90% IICC.

Correlating kdr genotype with resistance ratio

We regularly observed only six genotypes in the field (Table 4). Only a combined seven mosquitoes of the 4,810 analyzed in this study had genotypes of IIFF, IIFC, or VIFF, which would require an allele coding for a resistant isoleucine at 1016 (1016I) and a susceptible phenylalanine at 1534 (1534F) to be contributed from at least one parent. We did commonly observe the reverse where homozygous or heterozygous susceptible 1016V alleles paired with the homozygous resistant 1534C allele (VVCC or VICC). Plotting the combined dataset for strains with both kdr genotype frequencies and resistance ratios from topical application of permethrin showed a strong correlation between increasing RR and increasing 1016II frequency (Fig 4, Pearson correlation coefficient = 0.905, P<0.00001, n = 20). Strains from Miami Beach, Fort Myers, and New Port Richey, with high frequencies of 1016II, which, based on the rarity of IIFF and IIFC mosquitoes implies the genotype IICC, were the strains with the highest resistance ratios. In contrast, the strain from Big Coppitt Key had the lowest 1016II frequency and the lowest RR even though the 1534CC frequency was relatively high.

thumbnail
Fig 4. Regression of genotype IICC frequency to resistance ratio.

Plot of RR versus IICC frequency indicates a strong correlation between the two factors (Pearson correlation coefficient ρ = 0.905). Comparison of RR to other genotypes did not indicate a strong correlation.

https://doi.org/10.1371/journal.pntd.0006544.g004

thumbnail
Table 4. Combined genotype frequencies for 1016 and 1534 alleles in 4,810 Florida Ae. aegypti.

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

Discussion

Recent local transmission of Zika virus during 2016 and small outbreaks of dengue virus in 2010 and 2011 demonstrate that an effective IVM plan that attacks multiple life stages to reduce mosquito numbers is a necessity. However, when disease transmission is active, chemical adulticiding can be the only means to immediately reduce the population of potentially infected mosquitoes. While there is some debate as to overall efficacy of adulticiding against Ae. aegypti and Ae. albopictus, there is little question that it is a necessary part of the response that allows other slower methods of control like larviciding and source reduction to gain a foothold. Pyrethroids are the major class of chemical adulticides that Florida mosquito control programs use operationally. Thus, the efficacy of pyrethroids on Aedes vectors is a critically important part of a control program.

Direct topical application of permethrin clearly demonstrated that resistance was widespread in Ae. aegypti strains throughout Florida. Resistance ratios ranged from about 6-fold in strains from Manatee County and Big Coppitt Key to approximately 60-fold in Miami Beach, Cape Coral, and New Port Richey. Genotyping thousands of individuals indicated that common kdr alleles were also widely distributed in the state. While these findings are noteworthy as they represent the first published report of widespread pyrethroid resistance and kdr alleles in the southeast US, this result is not surprising and could be predicted based on the results from other locations including Puerto Rico, Mexico and, more recently, in invasive CA Ae. aegypti strains [11, 13, 20].

The dataset described in this study does reveal wide variations in both RR and kdr alleles within small geographic areas. Examining Miami-Dade County, the strain collected from Miami Beach was highly resistant and had high levels of kdr. However, strains from inland Wynwood and Little River were much less resistant than the Miami Beach strain. Even strains from the opposite ends of neighborhoods differed. We saw this disparate pattern in south Miami-Dade, Manatee, and Lee counties. In Pasco County, we saw very different genotypes in strains geographically close to one another, separated only by a major highway. This wide variation in resistance alleles has been observed in Mexico and has been proposed to predominate over natural gene flow [34]. Considering this variability, along with the relatively short flight ranges and limited immigration in Ae. aegypti [35], the very different allele frequencies we observed in geographically close strains would support performing testing in numerous areas of a control district to get an accurate resistance picture.

In contrast to the resistance in the Ae. aegypti strains, we observed very little permethrin resistance in Ae. albopictus statewide. This was true whether Ae. aegypti were present or absent from the same collections. Waits et al. [24] showed very low levels of resistance in St. Johns County strains collected from areas without Ae. aegypti. Miami-Dade, Orange, and Lee counties all had the same pattern of resistant Ae. aegypti and much less resistant Ae. albopictus. It has been proposed that the development of pyrethroid resistance is much more difficult in Ae. albopictus due to sequence differences in the voltage-gated sodium channel (NaV) that make kdr much less likely, although recent reports indicate it may be possible [36, 37]. Our laboratory efforts to induce permethrin resistance in wildtype Florida Ae. albopictus by gentle pressuring have failed. At this time, pyrethroid resistance in Florida Ae. albopictus does not appear to be an issue that could lead to adulticide failure.

An important observation made due to this work is the correlation between increasing RR and the frequency of the IICC genotype in Ae. aegypti, which has been anecdotally observed in other studies using the CDC bottle bioassay [6, 11, 13, 38]. The linkage between kdr mutations and a strong correlation with resistance ratios has also been observed in other dipterans [39, 40]. While more research must be done to validate the correlation, this dataset adds another 20 strains with both resistance and kdr data that support using kdr genotype as a surrogate to estimate pyrethroid resistance levels in mosquitoes [9, 41]. The use of allele frequencies has several potential benefits. Genotype data are relatively easy to collect, the results are produced in hours, and dead mosquitoes collected during standard surveillance activities can be used to provide information on resistance. With limited budgets and personnel at the operational vector control district level, predictive estimation of resistance levels could produce useful operational data from activities already being done without requiring additional efforts to collect or produce mosquitoes for bioassay testing.

Nearly a decade after Donnelly et al. [41] asked whether kdr genotypes are predictive in mosquitoes there is, to our knowledge, no published study that shows a pyrethroid-resistant strain of Ae. aegypti without also showing the presence of kdr alleles. We suggest that the major benefits to be gained from use of allele frequencies as an estimator of resistance would be improvements in coverage area (this study, 62 strains with allele frequencies vs. 26 by direct topical application) and better operational decision-making as vector control programs could access more timely information on area specific resistance levels. These challenges of getting wide coverage and providing this information on an operationally useful timeline have been present in recent response efforts to Zika virus. The efforts of CDC and vector control units in Puerto Rico and the efforts of the authors and others in Florida to develop wide-area pyrethroid resistance maps by relying strictly on bioassay data show that it is currently a slow, labor intensive process. In Florida at least, the epidemic had passed by the time more than a few strains had been tested by bioassay. Until there is reliable, published evidence to argue against the use of kdr frequencies as a predictor, it is at present the only rapid way to assess strains across a large area.

This study shows that permethrin resistance is widely present in variable intensity in Ae. aegypti throughout Florida. This variability points to the need to include resistance testing as part of an IVM plan as well as examine resistance in more than one location. But how do we use this resistance information to improve vector control? Clearly, the strain of Ae. aegypti in Miami Beach is very different from nearby Wynwood and would likely call for different treatment strategies. A treatment with permethrin would likely have much less effect in Miami Beach in comparison to other less resistant locations.

Our study also points to the value of a collaborative approach from motivated stakeholders to develop resistance information. The state of Florida began regular conferences to bring together vector control districts, researchers, and public health resources months before the first locally transmitted case of Zika was reported in 2016. Like the CDC response in Puerto Rico, development of resistance information was an early and ongoing part of this process. The dataset in this study represents the result of thousands of hours of effort from vector control districts, vector control contractors, the state of Florida, and the federal government to produce operationally useful resistance information to protect public health and improve the efficacy of control operations. However, even with these efforts, this study is very limited in scope. We examined only permethrin resistance and although the literature shows the patterns we observed would likely be applicable to other pyrethroids [13, 17], work to define statewide patterns of resistance to synergized products or organophosphates still needs to be addressed.

Supporting information

S1 File. Collection information and kdr data for strains used in this study.

This excel file contains detailed location information on collection sites, life stages tested, numbers of organisms tested and raw kdr allele data.

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

(XLS)

S2 File. Comparison of 1016 and 1534 kdr allele frequencies in adult and egg life stages of several strains of Ae. aegypti.

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

(DOCX)

Acknowledgments

The Authors would like to acknowledge the contributions of Rachel Bales, Lettie Cronin, J. Cotter, Lucy Li, and Jessica Louton for assistance with bioassays and genotyping.

Disclaimer: This work represents the opinions of the authors and does not represent an official position of any local, state or the US government or corporation.

References

  1. 1. Komalamisra N, Srisawat R, Apiwathanasorn C, Samung Y, Kaisri P. Residual effect of 10% bifenthrin wp on mosquitoes, and community acceptance, in eastern Thailand. Southeast Asian J Trop Med Public Health. 2009;40(6):1221. pmid:20578456
  2. 2. Soderlund DM, Bloomquist JR. Molecular mechanisms of insecticide resistance. In: Roush RT, Tabashnik BE, editors. Pesticide resistance in arthropods. Boston: Springer;1990. p. 58–96.
  3. 3. Georghiou GP. Overview of insecticide resistance. In: Green MB, LeBaron HM, Moberg WK, editors. ACS Symposium series—Managing Resistance to Agrochemicals. American Chemical Society (USA);1990. p. 18–41.
  4. 4. Smith LB, Kasai S, Scott JG. Pyrethroid resistance in Aedes aegypti and Aedes albopictus: Important mosquito vectors of human diseases. Pestic Biochem Physiol. 2016;133:1–12. pmid:27742355
  5. 5. Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S, Guillet P, et al. Pyrethroid and DDT cross—resistance in Aedes aegypti is correlated with novel mutations in the voltage—gated sodium channel gene. Med Vet Entomol. 2003;17(1):87–94. pmid:12680930
  6. 6. Vera-Maloof FZ, Saavedra-Rodriguez K, Elizondo-Quiroga AE, Lozano-Fuentes S, Black IV WC. Coevolution of the Ile1, 016 and Cys1, 534 mutations in the voltage gated sodium channel gene of Aedes aegypti in Mexico. PLoS Negl Trop Dis. 2015;9(12):e0004263. pmid:26658798
  7. 7. Haddi K, Tomé HV, Du Y, Valbon WR, Nomura Y, Martins GF, et al. Detection of a new pyrethroid resistance mutation (V410L) in the sodium channel of Aedes aegypti: a potential challenge for mosquito control. Sci Rep. 2017;7:46549. pmid:28422157
  8. 8. Saavedra-Rodriguez K, Maloof FV, Campbell CL, Garcia-Rejon J, Lenhart A, Penilla P, et al. Parallel evolution of vgsc mutations at domains IS6, IIS6 and IIIS6 in pyrethroid resistant Aedes aegypti from Mexico. Sci Rep. 2018;8:6747. pmid:29712956
  9. 9. Rodríguez MM, Bisset JA, Fernández D. Levels of insecticide resistance and resistance mechanisms in Aedes aegypti from some Latin American countries. J Am Mosq Control Assoc. 2007;23(4):420–429. pmid:18240518
  10. 10. Harris AF, Rajatileka S, Ranson H. Pyrethroid resistance in Aedes aegypti from Grand Cayman. Am J Trop Med Hyg. 2010;83(2):277–284. pmid:20682868
  11. 11. García GP, Flores AE, Fernández-Salas I, Saavedra-Rodríguez K, Reyes-Solis G, Lozano-Fuentes S, et al. Recent rapid rise of a permethrin knock down resistance allele in Aedes aegypti in Mexico. PLoS Negl Trop Dis. 2009;3(10):e531. pmid:19829709
  12. 12. Linss JG, Brito LP, Garcia GA, Araki AS, Bruno RV, Lima JB, et al. Distribution and dissemination of the Val1016Ile and Phe1534Cys Kdr mutations in Aedes aegypti Brazilian natural populations. Parasit Vectors. 2014;7(1):25.
  13. 13. Ponce-García G, Del Río-Galvan S, Barrera R, Saavedra-Rodriguez K, Villanueva-Segura K, Felix G, et al. Knockdown Resistance Mutations in Aedes aegypti (Diptera: Culicidae) From Puerto Rico. J Med Entomol. 2016;53(6):1410–1414. pmid:27493252
  14. 14. Centers for Disease Control and Prevention [Internet]. Atlanta: CDC; 2017 [cited 2018 Aug 06]. Available from: https://www.cdc.gov/zika/vector/testing-puertorico.html
  15. 15. Hemingway J, Boddington RG, Harris J, Dunbar SJ. Mechanisms of insecticide resistance in Aedes aegypti (L.)(Diptera: Culicidae) from Puerto Rico. Bull Entomol Res. 1989;79(1):123–130.
  16. 16. Reid WR, Thornton A, Pridgeon JW, Becnel JJ, Tang F, Estep A, et al. Transcriptional analysis of four family 4 P450s in a Puerto Rico strain of Aedes aegypti (Diptera: Culicidae) compared with an Orlando strain and their possible functional roles in permethrin resistance. J Med Entomol. 2014;51(3):605–615. pmid:24897853
  17. 17. Estep AS, Sanscrainte ND, Waits CM, Louton JE, Becnel JJ. Resistance status and resistance mechanisms in a strain of Aedes aegypti (Diptera: Culicidae) from Puerto Rico. J Med Entomol. 2017;54(6):1643–1648. pmid:28981681
  18. 18. Agramonte NM, Bloomquist JR, Bernier UR. Pyrethroid resistance alters the blood-feeding behavior in Puerto Rican Aedes aegypti mosquitoes exposed to treated fabric. PLoS Negl Trop Dis. 2017;11(9):e0005954. pmid:28931018
  19. 19. Vontas J, Kioulos E, Pavlidi N, Morou E, Della Torre A, Ranson H. Insecticide resistance in the major dengue vectors Aedes albopictus and Aedes aegypti. Pestic Biochem Physiol. 2012;104(2):126–131.
  20. 20. Cornel AJ, Holeman J, Nieman CC, Lee Y, Smith C, Amorino M, et al. Surveillance, insecticide resistance and control of an invasive Aedes aegypti (Diptera: Culicidae) population in California. F1000Res. 2016;5.
  21. 21. Richards SL, Balanay JA, Fields M, Vandock K. Baseline insecticide susceptibility screening against six active ingredients for Culex and Aedes (Diptera: Culicidae) mosquitoes in the United States. J Med Entomol. 2017;54(3):682–695. pmid:28399272
  22. 22. Connelly CR, and Parker CN. 2018. Statewide results of insecticide susceptibility testing of using the CDC Bottle Bioassay–what is next for Florida? In: Pesticide Resistance Management–from Theory to Practice, Symposium II. Paper 186. American Mosquito Control Association 84th Annual Meeting, Kansas City, MO. 1 March 2018.
  23. 23. Marcombe S, Farajollahi A, Healy SP, Clark GG, Fonseca DM. Insecticide resistance status of United States populations of Aedes albopictus and mechanisms involved. PLoS One. 2014;9(7):e101992. pmid:25013910
  24. 24. Waits CM, Fulcher A, Louton JE, Richardson AG, Becnel JJ, Xue RD, et al. A comparative analysis of resistance testing methods in Aedes albopictus (Diptera: Culicidae) from St. Johns County, Florida. Fla Entomol. 2017;100(3):571–577.
  25. 25. 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
  26. 26. Centers for Disease Control and Prevention [Internet]. Atlanta: CDC; 2016 [cited 2018 Aug 06]. Available from: https://www.cdc.gov/zika/vector/vector-control.html
  27. 27. World Health Organization. Guidelines for efficacy testing of insecticides for indoor and outdoor ground-applied space spray applications. 2009. Available from: http://apps.who.int/iris/bitstream/handle/10665/70070/WHO_HTM_NTD_WHOPES_2009.2_eng.pdf?sequence=1
  28. 28. Pridgeon JW, Pereira RM, Becnel JJ, Allan SA, Clark GG, Linthicum KJ. Susceptibility of Aedes aegypti, Culex quinquefasciatus Say, and Anopheles quadrimaculatus Say to 19 pesticides with different modes of action. J Med Entomol. 2008;45(1):82–87. pmid:18283946
  29. 29. Chang F, Dutta S, Becnel JJ, Estep AS, Mascal M. Synthesis of the insecticide prothrin and its analogues from biomass-derived 5-(chloromethyl) furfural. J Agric Food Chem. 2014;62(2):476–80. pmid:24354505
  30. 30. World Health Organization. Monitoring and managing insecticide resistance in Aedes mosquito populations: Interim guide for entomologists. 2016. Available from: http://apps.who.int/iris/bitstream/10665/204588/2/WHO_ZIKV_VC_16.1_eng.pdf
  31. 31. Chouaibou MS, Chabi J, Bingham GV, Knox TB, N’Dri L, Kesse NB, et al. Increase in susceptibility to insecticides with aging of wild Anopheles gambiae mosquitoes from Côte d’Ivoire. BMC Infect Dis. 2012;12(1):214.
  32. 32. Saavedra—Rodriguez K, Urdaneta—Marquez L, Rajatileka S, Moulton M, Flores AE, Fernandez—Salas I, et al. A mutation in the voltage—gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti. Insect Mol Biol. 2007;16(6):785–798. pmid:18093007
  33. 33. Yanola J, Somboon P, Walton C, Nachaiwieng W, Somwang P, Prapanthadara LA. High—throughput assays for detection of the F1534C mutation in the voltage—gated sodium channel gene in permethrin—resistant Aedes aegypti and the distribution of this mutation throughout Thailand. Trop Med Int Health. 2011;16(4):501–509. pmid:21342372
  34. 34. Saavedra-Rodriguez K, Beaty M, Lozano-Fuentes S, Denham S, Garcia-Rejon J, Reyes-Solis G, et al. Local evolution of pyrethroid resistance offsets gene flow among Aedes aegypti collections in Yucatan State, Mexico. Am J Trop Med Hyg. 2015;92(1):201–209. pmid:25371186
  35. 35. Harrington LC, Scott TW, Lerdthusnee K, Coleman RC, Costero A, Clark GG, et al. Dispersal of the dengue vector Aedes aegypti within and between rural communities. Am J Trop Med Hyg. 2005;72(2):209–220. pmid:15741559
  36. 36. Kasai S, Ng LC, Lam-Phua SG, Tang CS, Itokawa K, Komagata O, et al. First detection of a putative knockdown resistance gene in major mosquito vector, Aedes albopictus. Jpn J Infect Dis. 2011;64(3):217–221. pmid:21617306
  37. 37. Ishak IH, Jaal Z, Ranson H, Wondji CS. Contrasting patterns of insecticide resistance and knockdown resistance (kdr) in the dengue vectors Aedes aegypti and Aedes albopictus from Malaysia. Parasit Vectors. 2015;8(1):181.
  38. 38. Harris AF, Rajatileka S, Ranson H. Pyrethroid resistance in Aedes aegypti from Grand Cayman. Am J Trop Med Hyg. 2010;83(2):277–284. pmid:20682868
  39. 39. Jamroz RC, Guerrero FD, Kammlah DM, Kunz SE. Role of the kdr and super-kdr sodium channel mutations in pyrethroid resistance: correlation of allelic frequency to resistance level in wild and laboratory populations of horn flies (Haematobia irritans). Insect Biochem Mol Biol. 1998;28(12):1031–1037. pmid:9887519
  40. 40. Williamson MS, Martinez-Torres D, Hick CA, Devonshire AL. Identification of mutations in the houseflypara-type sodium channel gene associated with knockdown resistance (kdr) to pyrethroid insecticides. Mol Gen Genet. 1996;252(1–2):51–60. pmid:8804403
  41. 41. Donnelly MJ, Corbel V, Weetman D, Wilding CS, Williamson MS, Black IV WC. Does kdr genotype predict insecticide-resistance phenotype in mosquitoes?. Trends Parasitol. 2009;25(5):213–219. pmid:19369117