Research | Open | Published:
Relationship between insecticide resistance and kdr mutations in the dengue vector Aedes aegypti in Southern China
Parasites & Vectorsvolume 8, Article number: 325 (2015)
Aedes aegypti is an important vector for dengue virus and thus has been targeted with pyrethroid insecticides in many areas of the world. As such, resistance has been detected to several of these insecticides, including in China, but the mechanisms of the resistance are not well understood in this country.
Using the World Health Organization larval mosquito bioassay, five field populations of Aedes aegypti from Southern China were characterized for their resistance to cypermethrin and cyhalothrin. RNA extraction with PCR amplification, cloning and sequencing of the sodium channel gene was followed by comparisons of susceptible and wild mosquito strains Additionally, genomic DNA was used for Allele-specific PCR (AS-PCR) genotyping of the sodium channel genes to detect S989P, V1016G and F1534C mutations and allow for correlation analysis of resistance expression for the different mutations.
All wild strains expressed resistance to cypermethrin and cyhalothrin and the resistance expression between the two insecticides was highly correlated suggesting cross-resistance between these two pyrethroids. The AS-PCR technique effectively distinguished individual genotypes for all three mutations. Among the five wild strains tested, two strains carried all three mutations. Although the S989P and V1016G mutations were positively correlated to resistance expression of both pyrethroids, the F1534C mutation was negatively correlated.
Our methodology proved highly reliable and will aid future detection of kdr mutations. The three sodium channel mutations were common in the Ae. aegypti strains sampled from Southern China. The V1016G mutation appears to be the most important kdr mutation in Ae. aegypti strains in Southern China.
The mosquito Aedes aegypti (L.), a primary vector of dengue fever and dengue hemorrhagic fever, mainly occurs in tropical areas. Because there is currently no dengue vaccine, mosquito control has become the most important method of preventing transmission of these diseases  and insecticides are expected to remain a key component of dengue control for the foreseeable future . Pyrethroids, characterized by their high toxicity to target insects and relatively few adverse effects on mammals, have been used widely in mosquito control. However, many mosquito populations have developed resistance to these compounds following repeated exposure to them [3, 4].
Both pyrethroids and DDT target the voltage-gated sodium channel of insect neurons and single amino acid substitutions in the sodium channel have been associated with resistance to both these insecticides . This form of resistance, known as knockdown resistance (kdr), has been observed in several insect species including the mosquitoes Anopheles gambiae Giles , An. stephensi Liston , Culex quinquefasciatus Say  and Ae. aegypti .
In many insect species, mutations related to both pyrethroid and DDT resistance largely have been located in the IIS6 region. Of these, the replacement of a Leu for a Phe in the 1014 site (Leu1014Phe), first described in a DDT resistant strain of Musca domestica L. , was the most common. Other substitutions in the same homologous site have been described in a number of species, such as Leu1014Ser in the mosquitoes An. gambiae  and Culex pipiens L.  and Leu1014His in the tobacco budworm Heliothis virescens (Fabricius) . In Ae. aegypti, substitutions at site 1014 are unlikely, as two independent changes in the same codon would be necessary . Instead, mutations in different positions have been observed in Latin American and Southeast Asian populations of Ae. aegypti. So far, a number of kdr mutations have been identified in Ae. aegypti (S989P, I1011M/V, V1016G/I, F1269C, F1534C) [9, 14, 15] and at least two of these are known to be related to pyrethroid resistance; 1016 (Val to Ile or Gly) and 1534 (Phe to Cys) in the IIS6 and IIIS6 segments of the Vssc gene, respectively [9, 14, 16, 17].
Ae. aegypti is present in many residential areas of Southern China. In recent years, regular monitoring has revealed that many Chinese populations of Ae. aegypti have become resistant to pyrethroid insecticides [18, 19], but the underlying mechanism responsible for this resistance is unknown. The aims of this study were to detect and identify mutations in the sodium channel gene and to determine whether any of these are linked to pyrethroid resistance in wild Ae. aegypti populations in China. We used the WHO susceptibility test to evaluate the resistance of Ae. aegypti larvae to commonly used pyrethroid insecticides and molecular assays to identify specific point mutations in the sodium channel gene of Ae. aegypti collected from various regions in Southern China. The relationship between pyrethroid knockdown resistance and the frequencies of different kdr alleles is analysed and discussed. A novel method has also been developed to detect the mutation found in Chinese Ae. aegypti field strains and to determine the role of the mutations in cypermethrin and cyhalothrin resistance in China. Determining the distribution of one of the major mechanisms underpinning kdr mutations in wild mosquito populations is important to improving the effectiveness of dengue control operations.
Larvae of Ae. aegypti were collected from five sites in Southern China in 2013: Qishui (20°45′46.3″N, 109°45′35.0″E), Dongfang (19°6′3.8″N, 108°35′8.6″E), Jinghan (24°01′51.5″N, 97°53′15.6″E), Yundang (24°01′25.6″N, 97°52′58.9″E) and Munao (24°01′08.6″N, 97°51′47.8″E). Larvae were collected from standing water sources such as jars, bowls, tanks, tin cans and drums. No specific permissions were required to collect larvae as their collection did not involve, or affect, any endangered or protected species. Collected larvae were transported to the laboratory and colonized. Bioassays for pesticide resistance were conducted on the first (F1) generation of larvae produced in captivity and the results compared to those of a pyrethroid-susceptible Ae. aegypti strain. The susceptible Ae. aegypti strain originated from Hainan province, China, and has been isolated in a laboratory for more than 10 years without exposure to insecticides. Larvae were flash-frozen at −70 °C after pesticide exposure trials.
The resistance of larvae to two pyrethroids, cypermethrin (≥98 %, Sigma, USA) and cyhalothrin (>99 %, Sigma, USA), was evaluated. Thirty late 3rd and early 4th instar larvae were placed in a glass container that held 199 ml water and 1 mL of an insecticide concentration using the methods outlined by the WHO . All analytical grade insecticides were diluted five to seven concentrations with acetone. The serial concentrations were 0.0002, 0.0005, 0.001, 0.003, 0.005, 0.008 and 0.01 ppm for the susceptible strain, Qishui strain and Dongfang strain. The serial concentrations were 0.005, 0.01, 0.03, 0.05, 0.08 ppm for Jinghan strain and Yundang strain. The serial concentrations were 0.025, 0.05, 0.1, 0.15 and 0.2 ppm for Munao strain. A control with no insecticide was done in order to control for the natural mortality of the strain tested. All experiments were repeated three times and larval mortality recorded 24 h after treatment. Larval mortality was determined by dividing the number of dead larvae by the total number tested. No food was provided to larvae during bioassays, which were conducted under a 14 L: 10D photoperiod, 75 % relative humidity and temperature of 26 ± 1 °C.
LC50 values were calculated using Schoofs and Willhite’s  probit analysis program. As a measure of resistance we calculated the resistance ratio (RR) , which was the ratio of the estimated LC50 of the F1 captive generation to that of larvae of the susceptible strain. In order to examine whether resistance was correlated among sampling sites, relationships between estimated RR values of various insecticides were analyzed by Pearson correlation analysis. P-values <0.05 were considered statistically significant.
Extraction of mosquito RNA and synthesis of mosquito cDNA
Total RNA of 30 female mosquitoes from the susceptible laboratory strain and the five wild strains were extracted with TRIzol reagent according to the manufacturer’s protocol (Invitrogen, U.S.A.). These samples came from the same generation and age classes (late 3rd and early 4th instar) used for the resistance bioassay. The yield and purity of the extracted RNA were assessed by determining its absorbance (Abs) at 260 and 280 nm. RNA was used only if its Abs260 nm/Abs280 nm ratio was >1.8. Extracted total RNA was stored at −70 °C for later use.
Two microliters of total RNA from each sample was used to synthesize the first-strand cDNA using a oligo(dT) 12–18 primer in a 20-μl reaction following the Superscript TM RT-PCR System (Invitrogen, U.S.A.). The first-strand cDNA reaction was stored at −20 °C until required.
PCR amplification of the Ae. aegypti sodium channel gene
Gene-specific primers based on the published sequence of the Ae. aegypti para-sodium channel gene (GenBank Accession No.: EU399179.1) were designed in NCBI Primer BLAST and used to amplify the sodium channel gene of specimens from each strain. The sodium channel gene is 6554 bp in length and was amplified as 11 fragments. The primers used are shown in Table 1. The PCR reaction consisted of one cycle of 94 °C for 3 min and 35 cycles at 94 °C for 30 s, 46 ~ 63 °C for 30 s (different Tm for 11 fragments’ amplification), 72 °C for 1 min and a final extension step at 72 °C for 7 min.
Cloning and sequencing of PCR products
PCR products were run on a 1.5 % agarose gel and purified using Wizard PCR Preps DNA Purification System (Promega, U.S.A.). Purified PCR products were ligated into pMD 18-T easy vector (Takara, Japan) and plasmids cloned into E. coli XL1-blue competent cells. Plasmids were purified using a miniBEST plasmid purification kit (Takara, Japan). At least 30 clones of each PCR product were selected. Sequencing reactions were carried out on recombinant plasmids using M13 forward and reverse primers. Sequence alignment was performed with the Clustal W algorithm .
Comparison of sodium channel genes from susceptible and wild strains
The amino acid and cDNA sequences of susceptible and wild strains collected in 2013 were analysed with DNAstar software . Blast research was performed at http://www.ncbi.nlm.nih.gov/blast/. Sequences were aligned with ClustalX software  and viewed with GeneDoc software .
Genomic DNA isolation
Genomic DNA was isolated from 32 ~ 36 female mosquitoes from each strain using a Universal Genomic DNA Extraction Kit according to the manufacturer’s protocol (Takara, Japan). DNA was solubilized in a final volume of 100 μl ddH2O.
Allele-specific PCR (AS-PCR) genotyping of sodium channel genes
AS-PCR was used to detect S989P, V1016G and F1534C mutations in each of the five wild Ae. aegypti strains. A test using two PCR reactions for each individual was developed to diagnose each mutation (Fig. 1). The two reactions were exactly the same except that one contained a susceptible-specific primer and the other contained a mutant-specific primer. The outer primers and allele-specific inner primers for the three target mutations are listed in Table 2. A PCR diagnostic test was performed in accordance with the standard procedure with a total reaction volume of 25 μl, consisting of 80–100 ng genomic DNA, 14 mM Tris–HCl, pH 8.3, 70 mM KCl, 4.5 mM MgCl2, 0.15 mM each dNTP and 0.67 U TakaRa rTaq. PCR conditions were one cycle of 94 °C for 3 min, then 35 cycles of 94 °C for 30s, 60 or 62 °C (60 °C for V1016G and F1534C, 62 °C for S989P) for 30 s and 72 °C for 1 min, followed by one cycle of 72 °C for 7 min. PCR products were checked by electrophoresis on 1.5 % agarose gel in TAE buffer. Bands were visualized by ethidium bromide staining. The size of the PCR products for the detection of kdr alleles were 240 bp (S989P), 348 bp (V1016G) and 284 bp (F1534C), whereas the size of the products used as allele-nonspecific outer primers were 594 bp (S989P), 592 bp (V1016G) and 517 bp (F1534C).
Correlation of resistance with the frequencies of different mutations
In order to assess whether the frequency of different kdr alleles in Ae. aegypti strains was related to resistance, the relationship between allele frequency and the LC50 values for cypermethrin and cyhalothrin was analysed with SAS Software 6.22 (SAS Institute Inc.). Coefficients of correlation were calculated between the frequency of different mutations in different strains and the corresponding LC50 values for each pyrethroid. To determine whether the correlation was significant, the null hypothesis of r = 0, i.e. no correlation, was tested with a Student’s t-test. Significance levels of these tests were adjusted using the Bonferroni Correction  to compensate for the increased likelihood of obtaining significant results by chance.
Correlation of resistance to different insecticides
All wild strains had various levels of resistance to the insecticides tested (Table 3). The resistance ratio (RR) of each strain for cypermethrin and cyhalothrin ranged from 6.82 to 88.82-fold and 4.82 to 88.87-fold, respectively. The Munao strain had the highest resistance to these two pyrethroids. Resistance to cypermethrin and cyhalothrin was significantly correlated, i.e. strains with high resistance to cypermethrin tended to also be highly resistant to cyhalothrin (R2 > 0.9, P < 0.05) (Fig. 2) which indicates cross-resistance to these two insecticides.
Identification of sodium channel gene mutations
Three mutations were identified from the cDNA of a pooled sample of F1 specimens of wild Ae. aegypti: S989P (TCC-CCC), V1016G (GTA-GGA), and F1534C (TTC-TGC), compared to the susceptible strain.
AS-PCR genotyping of mosquitoes carrying kdr mutations
After optimizing experimental conditions, AS-PCR could effectively distinguish individual mosquitoes that were homozygous or heterozygous for the S989P, V1016G and F1534C mutations (Fig. 3). All five wild strains were tested and the mutation frequencies are summarized in Table 4. The Qishui and Dongfang strains had no mutations at sites S989 and V1016, however, three strains from Ruili (Jinghan, Yundang and Munao) had high frequencies of these mutations which reached 100 % in the Munao strain. At the F1534 locus, the opposite situation was apparent; mutation frequency was 100 % in the Qishui strain but very low in the three Ruili strains, indeed, there were no mutations at this locus in the Munao strain. The Jinghan and Yundang strains carried all three mutations.
Correlation between insecticide resistance and mutation frequency
A significant positive correlation between the frequency of the S989P and V1016G point mutations and resistance to cypermethrin and cyhalothrin was observed when comparing LC50 values and the frequency of these sodium channel gene mutations (Fig. 4). However, the frequency of the F1534C point mutation was negatively correlated with resistance to these two pyrethroids (Fig. 4).
Correlation between the S989P and V1016G mutations
The results show a strong linear relationship between the frequency of the S989P and V1016G mutations (Fig. 5). Co-occurrence of the S989P and V1016G mutations in five wild strains are presented (Table 5). The Qishui and Dongfang strains were all SS homozygotes and the Munao strain were all RR homozygotes, but the Jinghan and Yundang strains were either RS heterozygotes (13.9 % and 9.1 %) or RR homozygotes (66.7 %). In the remaining small proportion of individual mosquitoes, if the S989P mutation was detected as SS, the V1016G mutation was expressed as RS (2.8 % in Jinghan strain). Whereas, if the S989P mutation was expressed as RS, the V1016G mutation was expressed as RR (16.7 % in Jinghan strain, 24.2 % in Yundang strain).
Ae. aegypti is an important vector of dengue fever in China. It was previously thought that this species was only found in China south of latitude 22°N. Climate change, the rapid development of tourism and transportation and increasing urbanization have resulted in an increase in Ae. aegypti breeding sites and consequently a marked change in the species’ range in China. In recent years, the range of Ae. aegypti has expanded to approximately latitude 25°N .
The research presented in this paper is the first study of kdr resistance mechanisms in wild Ae. aegypti strains in China. We used the bioassay recommended by the WHO to determine the level of pyrethroid resistance of five wild strains of Ae. aegypti. The results show that the three wild strains obtained from Ruili in Yunnan (the Jinghan, Yundang and Munao strains) generally had high pyrethroid resistance. Highest and lowest levels of resistance to cypermethrin were 88.82 and 27.51 times, respectively, that of the susceptible strain, and for cyhalothrin, 88.87 and 30.7 times that of the susceptible strain. These high levels of resistance are probably related to extensive pesticide use in the Ruili region. Although Ae. aegypti was only first detected breeding in Yunnan in 2006 , pyrethroid insecticides have been extensively used since its discovery for mosquito control. The resultant selection pressure appears to have produced high levels of resistance to pyrethroids; our results show that cypermethrin and cyhalothrin are already unsuitable for controlling Ae. aegypti in the Ruili region.
Ae. aegypti collected from the Qishui and Dongfang were 7.13 and 6.82 times as resistant to cypermethrin, and 6.24 and 4.82 times as resistant to cyhalothrin, as the susceptible strain, which were much lower as compared to those of the three Ruili strains. The relatively low pyrethroid resistance in these Ae. aegypti strains may be the result of limited pyrethroid use in this economically-indigent area.
Our results indicate a significant positive correlation between the frequency of the S989 and V1016 mutations and resistance to cypermethrin and cyhalothrin in the five Ae. aegypti strains examined. As such, the higher the frequency of these mutations, the higher the level of resistance to cypermethrin and cyhalothrin one would expect. Lin et al.  found a correlation between the frequency of the V1016G kdr mutation and resistance to transfluthrin, d-allethrin, metofluthrin, esbiothrin and prallethrin. Moreover, the S989P mutation often co-exists with V1016G in pyrethroid resistant strains of Ae. aegypti . However, V1016G/S989P individuals are no more resistant to pyrethroids than those carrying just one of the mutations . This suggests that S989P does not enhance the resistance conferred by the V1016G mutation. Du et al.  speculated that S989P could compensate for any reduction in fitness caused by the V1016G mutation. Hirata et al.  states that S989P alone did not alter deltamethrin sensitivity, but when combined with V1016G, mosquitoes were ten times more sensitive to deltamethrin. Kawada et al.  states that S989P has always been found with V1016G, but V1016G has been found alone. S989P has never been found alone. From our results, the frequency of occurrence of S989P and V1016G heterozygotes and homozygotes, we speculate that the V1016G occurred before the S989P mutation.
Hu et al.  previously confirmed by the insertion of the 1534C equivalent mutation into the sodium channel gene of the cockroach that the F1534C mutation confers sodium channel resistance to type I, but not type II, pyrethroids. The F to C mutation drastically reduced sodium channel sensitivity to three type I pyrethroids; permethrin, NRDC 157 (a deltamethrin analogue lacking the a-cyano group) and bioresemthrin, but not to three type II pyrethroids (cypermethrin, deltamethrin and cyhalothrin). We found that the F1534C mutation was negatively correlated with resistance to cypermethrin and cyhalothrin, which suggests that the higher the frequency of this mutation in mosquito strains, the lower their resistance to these two insecticides; thus, the presence of this mutation not only does not enhance resistance to cypermethrin and cyhalothrin but in fact reduces it.
Research on Ae. aegypti in Thailand found, regardless of insecticide exposure status, that no homozygous 1016G mutants expressed the homozygous form of the 1534C mutation . Mosquitoes with the homozygous 1534C mutation are more susceptible to deltamethrin than those heterozygous and homozygous for 1534 F . This is a very interesting result which suggests that the F1534C mutation could potentially compensate for any reduction in fitness caused by the S989P and V1016G mutations.
The methods we established proved to be highly reliable and should aid future studies aimed at further determining the extent of kdr mutations worldwide. The frequency of the S989P, V1016G and F1534C mutations were correlated with resistance to cypermethrin and cyhalothrin, and all three mutations were found to be common in Ae. aegypti strains throughout Southern China. Our results also provide the first evidence of the S989P, V1016G and F1534C mutations in Ae. aegypti in Southern China, thereby increasing the knowledge of these mutations. The V1016G mutation appears to be the most important mutation in Ae. aegypti strains in Southern China at the present time, however, there exists the possibility that other mechanisms were also present and acting to confer resistance.
Gubler DJ. Dengue and dengue hemorrhagic fever. Clin Microbiol Rev. 1998;11:480–96.
Townson H, Nathan MB, Zaim M, Guillet P, Manga L, Bos R, et al. Exploiting the potential of vector control for disease prevention. Bull World Health Organ. 2005;83:942–7.
Rodriguez MM, Bisset J, de Fernandez DM, Lauza’ n L, Soca A. Detection of insecticide resistance in Aedes aegypti (Diptera: Culicidae) from Cuba and Venezuela. J Med Entomol. 2001;38:623–8.
Lin YH, Teng HJ, Ho CM, Pai HH, Wu SC, Hsu EL. Insecticide resistance in Aedes aegypti during dengue epidemics in Taiwan. Formosan Entomol. 2003;23:263–74.
Narahashi T. Neuronal ion channels as the target sites of insecticides. Parmacol Toxicol. 1996;78:1–14.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.
Enayati AA, Vatandoost H, Ladonni H, Townson H, Hemingway J. Molecular evidence for a kdr-like pyrethroid resistance mechanism in the malaria vector mosquito Anopheles stephensi. Med Vet Entomol. 2003;17:138–44.
Martinez-Torres D, Chevillon C, Brun-Barale A, Berge JB, Pasteur N, Pauron D. Voltage-dependent Na + channels in pyrethroid-resistant Culex pipiens L mosquitoes. Pest Sci. 1999;55:1012–20.
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:87–94.
Ingles PJ, Adams PM, Knipple DC, Soderlund DM. Characterization of voltage-sensitive sodium channel gene coding sequences from insecticide susceptible and knockdown-resistant house fly strains. Insect Biochem Mol Biol. 1996;26:319–26.
Pinto J, Lynd A, Elissa N, Donnelly MJ, Costa C, Gentile G, et al. Co-occurrence of East and West African kdr mutations suggests high levels of resistance to pyrethroid insecticides in Anopheles gambiae from Libreville, Gabon. Med Vet Entomol. 2006;20:27–32.
Chen L, Zhong D, Zhang D, Shi L, Zhou G, Gong M, et al. Molecular ecology of pyrethroid knockdown resistance in Culex pipiens pallens mosquitoes. PLoS ONE. 2010;5:e11681.
Park Y, Taylor MFJ. A novel mutation L1029H in sodium channel gene hscp associated with pyrethroid resistance for Heliothis virescens (Lepidoptera: Noctuidae). Insect Biochem Mol Biol. 1997;27:9–13.
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:785–98.
Kawada H, Higa Y, Komagata O, Kasai S, Tomita T, Thi Yen N, et al. Widespread distribution of a newly found point mutation in voltage-gated sodium channel in pyrethroid-resistant Aedes aegypti populations in Vietnam. PLoS Negl Trop Dis. 2009;3:e527.
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–5.
Harris AF, Rajatileka S, Ranson H. Pyrethroid resistance in Aedes aegypti from Grand Cayman. Am Trop Med Hyg. 2010;83:277–84.
Cui F, Raymond M, Qiao CL. Insecticide resistance in vector mosquitoes in China. Pest Manag Sci. 2006;62:1013–22.
Zeng LH, Sun DW, Zhao W, Li SG, Yang X. Determination of the susceptibility of Aedes aegypti and Aedes albopictus to commonly used insecticides in Hainan province. Chin J Vector Biol Control. 2010;21:148–9.
WHO [World Health Organization]. Report of the WHO Informal Consultation on the evaluation and testing of insecticides. Geneva, Switzerland: World Health Organization; 1996. CTD/WHO PES/IC/96.1.
Schoofs GM, Willhite CC. A probit analysis program for the personal computer. J Appl Toxicol. 1984;4:141–4.
Orshan L, Kelbert M, Pener H. Patterns of insecticide resistance in larval Culex pipiens populations in Israel: dynamics and trends. J Vector Ecol. 2005;30:289–94.
Aiyar A. The use of CLUSTAL W and CLUSTAL X for multiple sequence alignment. Methods Mol Biol. 2000;132:221–41.
Tompson JD, Higgins DG, Gibson TJ. Clustal-W-improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalities and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Rice WR. Analysing tables of statistical tests. Evolution. 1989;43:223–5.
Xie H, Zhou HN, Yang YM. Advances in the research on the primary dengue vector Aedes aegypti in China. Chin J Vector Biol Control. 2011;22:194–7.
Wang PY, Zhou HN, Wu C, Gu YA, Dong LM, Zhu GJ, et al. Investigation of dengue vector Aedes aegypti in Yunnan province. Chin J Vector Biol Control. 2006;17:507–8.
Lin YH, Tsen WL, Tien NY, Luo YP. Biochemical and molecular analyses to determine pyrethroid resistance in Aedes aegypti. Pestic Biochem Physiol. 2013;107:266–76.
Rinkevich FD, Du Y, Dong K. Diversity and convergence of sodium channel mutations involved in resistance to pyrethroids. Pestic Biochem Physiol. 2013;106:93–100.
Du Y, Nomura Y, Satar G, Hu Z, Nauen R, He SY, et al. Molecular evidence for dual pyrethroid-receptor sites on a mosquito sodium channel. Proc Natl Acad Sci U S A. 2013;110:11785–90.
Hirata K, Komagata O, Itokawa K, Yamamoto A, Tomita T, Kasai S. A single crossing-over event in voltage-sensitive Na + channel genes may cause critical failure of Dengue mosquito control by insecticides. PLoS Negl Trop Dis. 2014;8:e3085.
Kawada H, Oo SZM, Thaung S, Kawashima E, Maung YNM, Thu HM, et al. Co-occurrence of point mutations in the voltage-gated sodium channel of pyrethroid-resistant Aedes aegypti populations in Myanmar. PLoS Negl Trop Dis. 2014;8:e3032.
Hu Z, Du Y, Nomura Y, Dong K. A sodium channel mutation identified in Aedes aegypti selectively reduces cockroach sodium channel sensitivity to type I, but not type II pyrethroids. Insect Biochem Mol Biol. 2011;41:9–13.
Stenhouse SA, Plernsub S, Yanola J, Lumjuan N, Dantrakool A, Choochote W, et al. Detection of the V1016G mutation in the voltage-gated sodium channel gene of Aedes aegypti (Diptera: Culicidae) by allele-specific PCR assay, and its distribution and effect on deltamethrin resistance in Thailand. Parasit Vectors. 2013;6:253.
We would like to thank Lois Wood, Mark Halvorsen and Chris Holderman for assistance with laboratory procedures conducted in the University of Florida, U.S.A. This work was funded by grants from the Infective Diseases Prevention and Cure Project of National Ministry of Public Health of China (No. 2012ZX10004219).
The authors declare that they have no competing interests.
CXL and TYZ jointly designed and coordinated the study, with the contribution of PEK and RDX. CXL drafted the article. MHZ and GW collected samples from different sites. TY, XXG, YMZ, YDD, DX and HDZ carried out the laboratory work and performed the statistical analysis. All authors read and approved the final manuscript submitted.