Skip to main content

Advertisement

High frequencies of F1534C and V1016I kdr mutations and association with pyrethroid resistance in Aedes aegypti from Somgandé (Ouagadougou), Burkina Faso

Article metrics

Abstract

Background

Resistance to pyrethroid insecticides involving kdr mutations is widespread in Aedes aegypti (L.) (Diptera: Culicidae) and potentially could impact control efforts in endemic countries. Dengue cases had been sporadic in Burkina Faso for over a decade prior to the 2016–2017 outbreak that resulted in 15,074 suspected cases and 36 deaths, mainly in Ouagadougou. These outbreaks highlighted the lack of information on numerous aspects of the biology, behaviour and insecticide status of local dengue vector populations that are fundamental to vector control.

Results

We investigated the insecticide resistance profiles and the kdr mutations involved in pyrethroid resistance of Ae. aegypti from Somgandé, a district of Ouagadougou. WHO bioassays revealed that the local Ae. aegypti populations were highly resistant to pyrethroids with mortalities of 15% for permethrin and 37% for deltamethrin. Resistance to carbamates was also detected with mortalities of 55% for propoxur and 90% for bendiocarb, but high mortalities (> 97%) to organophosphates (malathion and fenitrothion) indicated susceptibility. Allele-specific PCR and voltage-gated sodium channel gene sequencing showed a very high frequency (97%) of the F1534C kdr allele whilst the V1016I kdr mutation frequency was 46%. Association of dual-locus kdr mutations was detected for permethrin resistance.

Conclusion

We conclude that in this locality of Burkina Faso, Ae. aegypti is resistant to pyrethroid and carbamate insecticides but remains susceptible to organophosphates, providing useful information for possible future control.

Background

Dengue occurs mainly in tropical and sub-tropical areas, especially in urban and semi-urban areas [1], and is a threat to approximately half of the world’s population [2]. In recent decades, dengue cases have increased significantly, largely due to rapidly growing urbanisation that creates vast areas within which the vector Aedes aegypti proliferates [3]. Dengue is becoming an acute public health problem [4] in low- and middle-income countries like Burkina Faso which experienced only sporadic cases or rare outbreaks [5, 6] until 2006–2008, when 683 confirmed cases of dengue fever were reported in Ouagadougou in the centre of Burkina Faso and Nouna in the north west [7]. Outbreaks followed with many cases in 2013, in 2016 and again in 2017, when 15,074 suspected cases and 36 deaths occurred mainly in Ouagadougou [8], within which the district of Somgandé has been highlighted as a dengue hotspot since 2013 [9].

Dengue control and prevention rely mainly on vector control, using primarily chemical insecticide classes such as organophosphates and pyrethroids for space spraying and pyrethroid-treated materials for personal protection [10, 11]. Pyrethroid insecticides act on the voltage-gated sodium channel (VGSC) causing depolarization, which disrupts electrical signalling in the nervous system [12, 13] resulting in rapid paralysis and death of mosquitoes known as “knock down”. Aedes aegypti has developed resistance to pyrethroids [14] as a result of non-synonymous kdr mutations in the VGSC, and via metabolic resistance [15]. To date, only two kdr mutations have been found in Africa: V1016I and F1534C [16]. The V1016I mutation in domain II of the VGSC is common in the Americas [15], but has been detected only at very low frequency in Ghana [17]. The F1534C mutation in domain IV in VGSC has a worldwide distribution [15] and was also recently detected in Ghana, but at much higher frequency than V1016I [17]. Metabolic resistance is likely to also have an important role in pyrethroid resistance in Ae. aegypti, and many genes have been shown to be overexpressed in transcriptomic studies [15]. Some have been demonstrated as being capable of metabolism via in vitro expression studies [18], though quantifying their contribution to resistant phenotypes is difficult and remains to be clarified.

To date, the two kdr mutations found in Africa involve nucleotide changes from T to G at codon 1534 resulting in the replacement of phenylalanine by cysteine [19], and from G to A at the codon 1016 resulting in the replacement of valine by isoleucine [20]. These mutants are associated in different ways with pyrethroid resistance in Ae. aegypti. The F1534C mutation often occurs alone, and its association with resistance to the type I pyrethroid permethrin has been established in vitro [21] and in many field studies [15]. However, association of the V1016I mutation alone with pyrethroid resistance has not been clearly established; in both Mexico and Brazil, co-occurrence of 1016I and 1534C is common and likely driven by selection for resistance [22, 23]. The main objective of this study was to assess the resistance status of Ae. aegypti to pyrethroids, organophosphates and carbamates and to assess the presence and contribution of kdr mutations to pyrethroid resistance in the dengue hotspot of Somgandé district in Ouagadougou (Burkina Faso).

Materials and methods

Larval collection

Larvae were collected during the 2016 and 2017 rainy seasons in Somgandé, a district situated in the eastern part of Ouagadougou where relatively high numbers of dengue cases have been recorded since 2013. Ae. aegypti larvae were collected from artificial outdoor breeding sites such as tires, drums, jars and bowls. Larvae were transported to the laboratory and reared to adults, which were fed with 10% sugar solution. All life stages were maintained in an insectary with 28 ± 1 °C temperature and 68 ± 4% relative humidity.

Insecticide bioassays

Insecticide bioassays were carried out according to WHO guidelines using the tube susceptibility test kit [24]. Six insecticides were tested: permethrin 0.75%, deltamethrin 0.05%, malathion 5%, fenitrothion 1%, bendiocarb 0.1% and propoxur 0.1%. It should be noted that whilst the diagnostic concentrations for fenitrothion, bendiocarb and propoxur are correct for Ae. aegypti and that for deltamethrin very close to a provisionally approved concentration, those for malathion and permethrin are higher than the recommended doses [24]. However, in both of these cases, we used the papers applied in the vast majority of other studies of Aedes worldwide [15]. Between 80 and 100 unfed 3- to 5-day-old female Ae. aegypti were exposed to insecticide-impregnated papers for 1 h. Bioassays were carried out at 28 ± 1 °C with 79 ± 7% relative humidity. After the exposure, mosquitoes were transferred to holding tubes and the number of dead mosquitoes was recorded after 24 h. Mortality rates were adjusted with Abbott’s correction when the mortality in the control tubes was between 5 and 20%. The susceptible Ae. aegypti Liverpool strain was used as control in bioassays to ensure appropriate performance of the tests (i.e. 100% mortality expected in this strain). At the end of the 24-h post-exposure period, live mosquitoes were killed by freezing and stored in silica gel-containing tubes at − 20 °C.

DNA extraction

Each mosquito was homogenised using a sterilised pestle in 200 μl of a buffer comprising of 0.1 M Tris (pH 9.0), 0.1 M ethylene diamine tetraacetic acid (EDTA), 1% sodium dodecyl sulfate (SDS) and 0.5% diethylpyrocarbonate (DEPC) in a 1.5-ml microtube. The homogenate was incubated at 70 °C for 30 min using a block incubator (ASTEC, B1-516). Following the addition of 44.8 μl of 5 M potassium acetate (KOAc), the mixture was vortexed and kept on ice for 30 min. The mixture was then centrifuged for 15 min at 20,000g at 4 °C. A volume of 180 μl of the supernatant was collected in a new 1.5-ml microtube, and 90 μl of isopropanol was added. After vortexing, the mixture was stored at − 25 °C for 30 min and centrifuged at 20,000×g at 4 °C for 20 min. The supernatant was discarded, and the DNA pellet was rinsed by adding 200 μl of 70% ethanol. Following centrifugation for 5 min at 20,000×g, the pellet was dried at room temperature and dissolved in 20 μl of TE (Tris-EDTA, pH 8.0) buffer.

Genotyping of F1534C and V1016I mutations

Allele-specific PCR (AS-PCR) followed published protocols for the V1016I (25) and the F1534C (26) mutations in each case with minor modifications.

For detection of the V1016I kdr mutation, a reaction volume of 12.5 μl contained 1 μl of target DNA; 2.5 μl of primer mixture comprising of Val1016f (0.25 μM), Iso1016f (0.25 μM) and the Iso1016r common reverse primer (0.125 μM) (Table 1); 2.75 μl of sterile water; and 6.25 μl of Taq polymerase mixture (AmpliTaq Gold® 360 Master Mix, Thermo Fisher Scientific). PCR cycling conditions were initial denaturation for 10 min at 95 °C, 30 cycles of extension at 95 °C for 30 s, 62 °C for 1 min and 72 °C for 45 s, and a final elongation at 72 °C for 5 min.

Table 1 List of primer sequences used for detecting allele-specific kdr mutations and direct sequencing of the sodium channel regions encompassing these mutations

Genotyping of the F1534C kdr mutation required two sets of PCR reactions: the first used a mix of primers M3-For, M3-Rev and M3-F and the second used a mix of M3-For, M3-Rev and M3-C for detection of the 1534F and 1534C alleles, respectively (Table 1). Each PCR reaction volume was 12.5 μl, containing 1 μl of target DNA, 2.5 μl of primer mixture (0.5 μM of final concentration for each primer), 2.75 μl of sterile water and 6.25 μl of Taq polymerase mixture (AmpliTaq Gold® 360 Master Mix, Thermo Fisher Scientific). PCR cycling conditions were initial denaturation for 10 min at 95 °C, 35 cycles of extension at 95 °C for 30 s, 60 °C for 1 min and 72 °C for 30 s, and final elongation at 72 °C for 7 min.

After amplification, the PCR products were mixed with 2.5 μl of 6× Loading Dye Purple Buffer (New England Biolabs), run on either a 3% (V1016I) or a 1.5% (F1534C) agarose gel in TAE buffer (Nacalai, Tesque, Inc., Kyoto, Japan) and stained with ethidium bromide solution for UV visualisation. For the reaction detecting V1016I, the sizes of amplified products were 98 bp for wild-type alleles and 78 bp for mutant alleles [25]. For the reaction detecting F1534C, product sizes for both wild type and mutant alleles were 284 bp [26].

To confirm the results from the allele-specific PCRs, some of the DNA samples were sequenced using primers described by Saingamsook et al. [27] (Table 1) for sequencing the VGSC region. Each PCR reaction mix contained 1 μl of DNA, 0.75 μl of 10 μM forward and reverse primers (IIPF and IIS6_R for V1016I and Ge-IIIS6_F and IIIS6_R for F1534C) (Table 1), 2.5 μl of dNTPs mix (2.5 mM), 1.5 μl of MgSO4 (25 mM), 2.5 μl of 10× PCR buffer and 0.5 μl of KOD-plus-neo DNA polymerase (1 U). The mixture was completed up to 25 μl with sterile water. PCR cycling conditions were denaturation at 94 °C for 1 min, 35 cycles of extension at 98 °C for 10 s, 68 °C for 1 min, and a final elongation at 68 °C for 3 min. PCR products were mixed with 5 μl of 6× Loading Dye Purple Buffer (New England Biolabs) and run on a 1.5% agarose gel (Nacalai, Tesque, Inc., Kyoto, Japan). Amplified band of either 581 bp for V1016I or 635 bp for F1534C were cut from the agarose gel for target DNA purification using the QIAEX® II Gel Extraction Kit according to the manufacturer’s protocol.

Sequencing reactions are comprised of 1 μl of purified DNA (5 ng/ul), 0.8 μl of 1 μM of forward (IIP_F or Ge-IIIS6_F) or reverse primer (IIS6_R or IIIS6_R), 2 μl of 5× Sequencing Buffer (Big Dye® Terminator V1.1, V3.1, Applied Biosystems), 0.5 μl of Cycle Sequencing RR-100 (Big Dye® Terminator V3.1, Applied Biosystems) and sterile water for a total reaction volume of 10 μl. Reaction conditions were denaturation for 1 min at 96 °C, 25 cycles of extension at 96 °C for 10 s, 50 °C for 5 s and 60 °C for 4 min. The product was run on an Applied Biosystems 3130xl Genetic Analyser, and sequencing chromatograms were analysed using Sequence Scanner 2 (Applied Biosystems).

Statistical analysis

Susceptibility bioassay results were interpreted according to WHO criteria [24] after 24 h: 98–100% mortality indicates susceptibility, 90–97% mortality suggests suspected resistance and mortality less than 90% indicates resistance. Chi-square tests of independence were performed to compare the mortality rates in 2016 and 2017. Kdr allele and genotype frequencies were calculated, and their 95% confidence interval was determined based on Agresti et al. [28]. A chi-square test for deviation from Hardy-Weinberg equilibrium was performed using R with Genetics package version 1.3.8.1. Fisher’s exact test was computed for association between genotype and the resistance/susceptible phenotype using R software version 3.4.3. Cohen’s kappa coefficient (κ) was calculated to measure the agreement between the AS-PCR and sequencing methods [29].

Results

Insecticide susceptibility

Bioassay results using the WHO method are summarised in Table 2. Low mortality was observed for pyrethroids with mortalities of 15% in 2016 and 26% in 2017 for permethrin and 20% in 2016 and 36% in 2017 for deltamethrin. Propoxur and bendiocarb were more toxic with mortalities of 44% in 2016 and 55% in 2017 for propoxur and 82% in 2016 and 92% in 2017 for bendiocarb. No difference in susceptibility was observed between 2016 and 2017 for all the insecticides tested (χ2 < 2, p > 0.16). High mortality was observed for the organophosphates fenitrothion and malathion insecticides (mortality between 97 and 100%). Comparatively, the exposure of the Liverpool susceptible Ae. aegypti strain to permethrin, the less effective insecticide tested, resulted in 100% mortality exposure, showing full susceptibility to this insecticide. These results indicate that the Ae. aegypti population from Somgandé is susceptible to organophosphates (or very nearly so) but resistant to permethrin, deltamethrin, bendiocarb and propoxur (Fig. 1).

Table 2 Ae. aegypti mortality rates after exposure to insecticides in 2016 and 2017 using WHO tubes standard bioassay guidelines. The chi-square is calculated, and the probability comparing 2016 and 2017 data and insecticides of the same family are shown
Fig. 1
figure1

Mortality (with standard errors) of Ae. aegypti following insecticide exposure in WHO tubes. The green dashed line indicates the threshold for susceptibility, and the black dashed line indicates the threshold for resistance

Allele and genotype frequencies of the V1016I and F1534C kdr mutations

In total, 578 and 554 mosquitoes were genotyped using the allele-specific PCR for the F1534C and V1016I mutations, respectively. At codon 1016, both wild-type (V/V) and mutant genotypes (V/I or I/I) were detected, whereas at codon 1534, wild-type genotypes (F/F) were absent. The majority (93%) of mosquitoes were homozygous 1534 C/C, and 21% homozygous 1016 I/I. Allele frequencies were 0.97 (95% CI 0.954–0.975) for 1534C and 0.46 (95% CI 0.429–0.488) for 1016I (Table 3).

Table 3 V1016I and F1534C genotype numbers and frequencies in parenthesis and the allelic frequency of the I and C mutations of Ae. aegypti. The chi-square and the probability comparing the genotypes is given

Only six from a possible nine genotypes across the two kdr mutations were identified from 552 Ae. aegypti mosquitoes from which genotypes at both loci were available (VV/FC, VV/CC, VI/FC, VI/CC, II/FC and II/CC) (Fig. 2). The genotype VI/CC was the most frequently observed (47%, n = 259), and 20% of individuals were homozygous for both kdr mutations (II/CC) (Fig. 2). Overall, the genotype frequency distribution was strongly deviant from Hardy-Weinberg equilibrium (χ2 < 1, p > 0.34) (Table 3). Forty-eight samples were sequenced for the sodium channel regions encompassing these mutations, and the sequencing data confirmed all six possible genotypes detected by AS-PCR.

Fig. 2
figure2

Frequency of the combined V1016I and F1534C genotypes where VV/FC, VV/CC, VI/FC, VI/CC, II/FC and II/CC mean respectively and 1016VV/1534CC, 1016VV/1534CC, 1016VI/1534FC, 1016VI/1534CC, 1016II/1534FC and 1016II/1534CC genotypes

Association of genotypes at V1016I/F1534C with pyrethroid resistance

As shown in Table 4, the frequency of the VV/CC genotype was relatively lower in live (deltamethrin, 11%; permethrin, 16%) than in dead mosquitoes (deltamethrin, 32%; permethrin, 29%). In contrast, the genotype frequency of VI/CC was higher in live mosquitoes (deltamethrin, 61%; permethrin, 52%), and the double mutant genotype II/CC was also detected more frequently in live mosquitoes (deltamethrin, 26%; permethrin, 28%). The association between the dual-locus kdr genotypes was significant for permethrin resistance (p < 0.001), but not deltamethrin resistance (p = 0.116) (Table 4). However, caution is required as the power of the test was limited by the low number of genotypes for dead mosquitoes.

Table 4 Co-occurrence of V1016I and F1534C mutations in pyrethroid-exposed mosquitoes. The number of dead and live mosquitoes and their frequencies in permethrin and deltamethrin exposure Ae. aegypti. The probability (p value) of genotype distribution in dead and alive mosquitoes is calculated

Discussion

This study investigated insecticide resistance of Ae. aegypti in Ouagadougou, the capital of Burkina Faso, and indicates that this important arbovirus vector is resistant to pyrethroid and carbamate insecticides but remains susceptible to organophosphates. Although the reasons for these patterns are not known, there are a number of events that may have contributed. Long-lasting insecticidal nets (LLINs) have been promoted by WHO [30] for malaria control throughout Africa with dramatic results in recent years [31]. Scaling up of the distribution of these pyrethroid-treated bednets in Burkina Faso started in 2010 and led to millions of nets being deployed throughout the country. Despite impacts on malaria mortality and morbidity, scaling up probably accelerated the development of pyrethroid resistance in populations of malaria vectors [32] and, given its peridomestic habits, may also have contributed to the development of pyrethroid resistance in Ae. aegypti. Selection pressure from widespread use of pyrethrum coils and insecticide aerosols in households is also likely [33, 34]. Resistance to carbamates is more difficult to explain because this class of insecticides is not used for vector control in Ouagadougou, though they are used for pest control in urban agriculture in Ouagadougou [35]. In Cameroon, deltamethrin and bendiocarb resistance was recorded in Ae. aegypti by Kamgang et al. [36] whilst 6 years earlier Ae. aegypti was found susceptible to deltamethrin and propoxur [37], suggesting a loss of susceptibility to these insecticides with time. However, in our study, no significant change in susceptibility is observed between 2016 and 2017 to all the insecticides tested, a short period indeed.

Resistance profiles from other African countries are quite variable with respect to carbamates, more common in pyrethroids, but remain very rare for organophosphates [16].

Allele-specific PCR genotyping of Ae. aegypti population from Somgandé revealed the occurrence of V1016I and F1534C mutations, the latter at very high frequency. The allelic frequencies of both mutations are higher than in Ghana, from where these two mutations were first reported in mainland Africa [17]. The near fixation of the F1534C allele in the Ae. aegypti population of Burkina resulted in an absence of genotypes including only 1016 kdr mutations, precluding independent analysis. When analysed together, a significant association with permethrin resistance of dual-kdr genotypes was detected in Ae. aegypti from Somgandé. This is consistent with results from Linss et al. [22] and also the model suggested by Vera-Maloof et al. [38] which suggests a rise in frequency of 1534C followed by a rise of 1016I, with concomitantly increased resistance. No association was detected in our study between kdr and deltamethrin resistance, and whilst power was limited because of the relatively small number of dead genotypes available, this was actually even more so for permethrin, suggesting that association between these kdr genotypes and deltamethrin might be weaker than for permethrin. However, the lack of any association demonstrated for deltamethrin should be taken with caution, and further work with larger sample sizes is required.

The V1016I kdr mutation is often found in association with the 1534C mutation in Southern and Central America [38,39,40] and in the Caribbean [41]. A similar profile of extreme 1534C frequency and moderate 1016I frequency has been detected in Madeira and is thought to have been imported recently [42]. Further investigation into the origin of the 1016I mutation in Burkina Faso is required as well as investigation of the role of metabolic resistance in Ae. aegypti [15].

Conclusion

The study has demonstrated that Ae. aegypti from a locality of Ougadougou are resistant to pyrethroids and carbamates but susceptible to organophosphates. The V1016I and especially F1534C kdr mutations were observed at high frequency and explain at least some of the observed pyrethroid resistance phenotypes. Further work is required to investigate spatial and longer-term temporal variations in the occurrence of these mechanisms and to identify other mechanisms involved in the resistance phenotypes. These results provide important information for the development of effective strategies for dengue control in Ouagadougou.

References

  1. 1.

    WHO. Comprehensive guidelines for prevention and control of dengue and dengue haemorrhagic fever. Regional Office for South-East Asia: World Health Organization; 2011. p. 196.

  2. 2.

    Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, et al. Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis. 2012;6:e1760.

  3. 3.

    Baruah S, Dutta P. Seasonal prevalence of Aedes aegypti in urban and industrial areas of Dibrugarh district. Assam Trop Biomed. 2013;30:434–43.

  4. 4.

    Murray NEA, Quam MB, Wilder-Smith A. Epidemiology of dengue: past, present and future prospects. Clin Epidemiol. 2013;5:299–309.

  5. 5.

    Were F. The dengue situation in Africa. Paediatr Int Child Health. 2012;32:18–21.

  6. 6.

    Messina JP, Brady OJ, Scott TW, Zou C, Pigott DM, Duda KA, et al. Global spread of dengue virus types: mapping the 70 year history. Trends Microbiol. 2014;22(3):138–46.

  7. 7.

    Collenberg E, Ouedraogo T, Ganame J, Fickenscher H, Kynast-Wolf G, Becher H, et al. Seroprevalence of six different viruses among pregnant women and blood donors in rural and urban Burkina Faso: a comparative analysis. J Med Virol. 2006;78:683–92.

  8. 8.

    WHO. Weekly bulletin and other emergencies. Week 44; 2017. p. 16p.

  9. 9.

    Ridde V, Agier I, Bonnet E, Carabali M, Dabiré R, Druetz T, et al. La Dengue Au Burkina Faso: Rech Interv Communaut pour l’Equite en sante au Burkina Faso 2014;1–4.

  10. 10.

    Esu E, Lenhart A, Smith L, Horstick O. Effectiveness of peridomestic space spraying with insecticide on dengue transmission; systematic review. Trop Med Int Heal. 2010;15(5):619–31.

  11. 11.

    WHO. Vector control operations framework for Zika virus. Geneva: World Health Organization; 2016. p. 10.

  12. 12.

    Catterall WA. Voltage-gated sodium channels at 60: structure, function and pathophysiology. J Physiol. 2012;590:2577–89.

  13. 13.

    Dong K, Du Y, Rinkevich F, Nomura Y, Xu P, Wang L, et al. Molecular biology of insect sodium channels and pyrethroid resistance. Insect Biochem Mol Biol. 2014;50:1–17.

  14. 14.

    Marcombe S, Darriet F, Tolosa M, Agnew P, Duchon S, Etienne M, et al. Pyrethroid resistance reduces the efficacy of space sprays for dengue control on the island of Martinique (Caribbean). PLoS Negl Trop Dis. 2011;5:e1202.

  15. 15.

    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:e0005625.

  16. 16.

    Weetman D, Kamgang B, Badolo A, Moyes C, Shearer F, Coulibaly M, et al. Aedes mosquitoes and Aedes-borne arboviruses in Africa: current and future threats. Int J Environ Res Public Health. 2018;15(2):220.

  17. 17.

    Kawada H, Higa Y, Futami K, Muranami Y, Kawashima E, Osei JHN, et al. Discovery of point mutations in the voltage-gated sodium channel from African Aedes aegypti populations: potential phylogenetic reasons for gene introgression. PLoS Negl Trop Dis. 2016;10:e0004780.

  18. 18.

    Stevenson BJ, Pignatelli P, Nikou D, Paine MJI. Pinpointing P450s associated with pyrethroid metabolism in the dengue vector , Aedes aegypti: developing new tools to combat insecticide resistance. PLoS Negl Trop Dis. 2012;6(3):e1595.

  19. 19.

    Yanola J, Somboon P, Walton C, Nachaiwieng W, Somwang P. 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. Trop Med Int Heal. 2011;16:501–9.

  20. 20.

    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.

  21. 21.

    Hirata K, Komagata O, Itokawa K, Yamamoto A, Tomita T. 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.

  22. 22.

    Linss JGB, Brito LP, Garcia GA, Araki AS, Bruno RV, Bento J, et al. Distribution and dissemination of the Val1016Ile and Phe1534Cys Kdr mutations in Aedes aegypti Brazilian natural populations. Parasit Vectors. 2014;7:25.

  23. 23.

    Vera-Maloof FZ, Saavedra-Rodriguez K, Elizondo-Quiroga AE, Lozano-Fuentes S, Black WC IV. 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.

  24. 24.

    WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva, Switzerland: WHO Document Production Services; 2016.

  25. 25.

    Martins AJ, Lima JBP, Peixoto AA, Valle D. Frequency of Val1016Ile mutation in the voltage-gated sodium channel gene of Aedes aegypti Brazilian populations. Trop Med Int Heal. 2009;14:1351–5.

  26. 26.

    Li C-X, Kaufman PE, Xue R, Zhao M, Wang G, Yan T, et al. Relationship between insecticide resistance and kdr mutations in the dengue vector Aedes aegypti in Southern China. Parasit Vectors. 2015;8:325.

  27. 27.

    Saingamsook J, Saeung A, Yanola J, Lumjuan N, Walton C, Somboon P. A multiplex PCR for detection of knockdown resistance mutations, V1016G and F1534C, in pyrethroid-resistant Aedes aegypti. Parasit Vectors. 2017;10:465.

  28. 28.

    Agresti A, Coull BA, Statistician TA, May N. Approximate is better than “exact” for interval estimation of binomial proportions. Am Stat. 1998;52:119–26.

  29. 29.

    Landis JR, Koch GG. The measurement of observer agreement for categorical data. Biometrics. 1977;33:159–74.

  30. 30.

    WHO. The Africa malaria report 2003. In: World Health Organization/UNICEF World Health Organization/UNICEF, vol. 2003; 2003. p. 1–118.

  31. 31.

    Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526(7572):207–11.

  32. 32.

    Toé KH, Jones CM, N’fale S, Ismai HM, Dabiré RK, Ranson H. Increased pyrethroid resistance in malaria vectors and decreased bed net effectiveness Burkina Faso. Emerg Infect Dis. 2014;20(10):1691–6.

  33. 33.

    Ogoma SB, Moore SJ, Maia MF. A systematic review of mosquito coils and passive emanators: defining recommendations for spatial repellency testing methodologies. Parasit Vectors. 2012;5:287.

  34. 34.

    Kudom AA, Mensah BA, Nunoo J. Assessment of anti mosquito measures in households and resistance status of Culex species in urban areas in southern Ghana: implications for the sustainability of ITN use. Asian Pac J Trop Med. 2013;6:859–64.

  35. 35.

    Tarnagda B, Tankoano A, François T, Sourabié Pane Bernadette AH, Ousman DAO, et al. Évaluation des pratiques agricoles des légumes feuilles: le cas des utilisations des pesticides et des intrants chimiques sur les sites maraîchers de Ouagadougou. J Appl Biosci. 2017;117:11658–68.

  36. 36.

    Kamgang B, Yougang AP, Tchoupo M, Riveron JM, Wondji C. Temporal distribution and insecticide resistance profile of two major arbovirus vectors Aedes aegypti and Aedes albopictus in Yaoundé, the capital city of Cameroon. Parasit Vectors. 2017;10:469.

  37. 37.

    Kamgang B, Marcombe S, Chandre F, Nchoutpouen E, Nwane P, Etang J, et al. Insecticide susceptibility of Aedes aegypti and Aedes albopictus in Central Africa. Parasit Vectors. 2011;4:79.

  38. 38.

    Vera-Maloof FZ, Saavedra-rodriguez K, Elizondo-quiroga AE, Lozano-fuentes S, Iv WCB. Mutations in the voltage gated sodium channel gene of Aedes aegypti in Mexico. PLoS Negl Trop Dis. 2015;9(12):1–22.

  39. 39.

    Collet ML, Frizzo C, Orlandin E, Rona LDP, Nascimento JC, Montano MAE, et al. Frequency of the Val1016Ile mutation on the kdr gene in Aedes aegypti (Diptera: Culicidae) in south Brazil. Genet Mol Res. 2016;15:gmr15048940.

  40. 40.

    Francis S, Saavedra-rodriguez K, Perera R, Paine M, Black C, Delgoda R. Insecticide resistance to permethrin and malathion and associated mechanisms in Aedes aegypti mosquitoes from St. Andrew Jamaica. PLoS One. 2017;12:1–13.

  41. 41.

    Goindin D, Delannay C, Gelasse A, Ramdini C, Gaude T, Faucon F, et al. Levels of insecticide resistance to deltamethrin, malathion, and temephos, and associated mechanisms in Aedes aegypti mosquitoes from the Guadeloupe and Saint Martin islands (French West Indies). Infect Dis Poverty. 2017;6:38.

  42. 42.

    Seixas G, Grigoraki L, Weetman D, Vicente JL, Silva AC, Pinto J, et al. Insecticide resistance is mediated by multiple mechanisms in recently introduced Aedes aegypti from Madeira Island (Portugal). PLoS Negl Trop Dis. 2017;11:e0005799.

Download references

Acknowledgements

We are also grateful to people of National Research Center for Protozoan Diseases, Obihiro, Japan, for supporting a stay in their Center.

Funding

This work was supported by a WHO/TDR grant (RCS-KM 2015 ID235974), and the International Collaborative Research Program for Tackling the NTDs Challenges in African countries from Japan Agency for Medical Research and Development, AMED (JP17jm0510002h0003). The study was also supported by JSPS KAKENHI (Grant Numbers JP26304038 to SF).

Availability of data and materials

The datasets including the bioassay raw data and PCR data collected during the current study are available from the corresponding author on request.

Author information

AB, PJM, DW, HK and SF conceived the project. AS, ES, SF, TS and HK designed the experiments which were performed by AS, ES, FY, TS, TS and SF. AS, AB and DW wrote the manuscript. All authors edited the draft and approved the final manuscript.

Correspondence to Athanase Badolo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sombié, A., Saiki, E., Yaméogo, F. et al. High frequencies of F1534C and V1016I kdr mutations and association with pyrethroid resistance in Aedes aegypti from Somgandé (Ouagadougou), Burkina Faso. Trop Med Health 47, 2 (2019) doi:10.1186/s41182-018-0134-5

Download citation

Keywords

  • Aedes aegypti
  • Insecticide
  • Resistance
  • kdr mutation
  • AS-PCR
  • Sequencing
  • Burkina Faso
  • F1534C
  • V1016I