- Research article
- Open Access
Investigating knockdown resistance (kdr) mechanism against pyrethroids/DDT in the malaria vector Anopheles funestus across Africa
BMC Genetics volume 18, Article number: 76 (2017)
Understanding the molecular basis of insecticide resistance is key to improve the surveillance and monitoring of malaria vector populations under control. In the major malaria vector Anopheles funestus, little is currently known about the role of the knockdown resistance (kdr) mechanism. Here, we investigated the presence and contribution of knockdown resistance (kdr) to pyrethroids/DDT resistance observed in Anopheles funestus across Africa.
Pyrosequencing genotyping and sequencing of the voltage gated sodium channel (VGSC) gene did not detect the common L1014F mutation in field collected An. funestus across Africa. Amplification and cloning of the full-length of the sodium channel gene in pyrethroid resistant mosquitoes revealed evidences of alternative splicing events with three transcripts of 2092, 2061 and 2117 amino acids (93% average similarity to An. gambiae). Several amino acid changes were detected close to the domain II of the protein such as L928R, F938 W, I939S, L802S and T1008 M. However, all these mutations are found at low frequency and their role in pyrethroid resistance could not be established. The presence of the exclusive alternative splicing at exon 19 was not associated with resistance phenotype. Analysis of patterns of genetic diversity of the VGSC gene revealed a high polymorphism level of this gene across Africa with no evidence of directional selection suggesting a limited role for knockdown resistance in pyrethroid resistance in An. funestus. Patterns of genetic differentiation correlate with previous observations of the existence of barriers to gene flow Africa-wide with southern population significantly differentiated from other regions.
Despite an apparent limited role of knockdown resistance in An. funestus, it is necessary to continue to monitor the contribution of the mutations detected here as increasing selection from insecticide-based interventions may change the dynamic in field populations as previously observed in other vectors.
Malaria prevention relies heavily on insecticide-based vector control interventions such as Long Lasting Insecticide Nets (LLINs) and Indoor Residual Spraying (IRS) . However increasing reports of insecticide resistance is a serious cause of concern for the long-term effectiveness of these control tools. In the major malaria vector An. funestus s.s., resistance to pyrethroid insecticides (only class recommended for bed net impregnation) is spreading across Africa with reports of resistance in all regions of the continent. Resistance to both type I and II pyrethroids firstly detected in southern Africa [2, 3] has since spread across the region [4,5,6,7,8,9]. Similar reports of resistance have been made in East Africa [10,11,12], in West Africa [13,14,15,16,17] and Central Africa [18, 19]. To ensure a continued effectiveness of the current control tools it is imperative to implement successful resistance management strategies. However, this requires a thorough understanding of the molecular basis of resistance and also the design of field applicable diagnostic tools to easily detect and monitor the spread of resistance in the field.
To date, investigations of the underlying molecular basis of all these reports of pyrethroid resistance have indicated that metabolic resistance through over-expression of cytochrome P450 genes is the major resistance mechanism in An. funestus s.s. [20,21,22,23,24]. Synergist assay with Piperonyl Butoxide (PBO), a P450 inhibitor, have consistently showed a near full recovery of susceptibility to pyrethroids when mosquitoes were pre-exposure to PBO [2, 12] suggesting a major role played by cytochrome P450s. This was further supported by in-depth genetic and molecular characterisation of metabolic resistance in both laboratory and field population of An. funestus s.s. in southern Africa which revealed that the duplicated cytochrome P450 genes, CYP6P9a and CYP6P9b beside CYP6M7 are playing key role in pyrethroid resistance [21, 23, 25,26,27].
No knockdown resistance (kdr) mutation has been detected in An. funestus s.s. in stark contrast to other mosquito species including the major malaria vector An. gambiae for which kdr mutations, notably at the 1014 codon, have been widely reported [28,29,30]. Indeed, despite sequencing of the Exon 20 fragment, no kdr mutation has been detected in An. funestus [8, 12,13,14] although some amino acid changes (I877L, V881 L and A1007S in Cameroon, F1021C in Uganda) have been recently reported but with no association with resistance [12, 18]. However, a comparative analysis of the haplotype diversity of the voltage-gated sodium channel (VGSC) gene between permethrin resistant and susceptible mosquitoes from Malawi showed a possible correlation between resistance phenotype and haplotype distribution, suggesting that other kdr mutations different to the classical L1014F/S mutations could be associated with pyrethroid resistance in An. funestus s.s. [8, 9]. Indeed, in other insect species, different kdr resistance mutations than the L1014F have been detected such as in Aedes aegypti with mutations in 1011, 1016 and 1534 codon positions [31, 32] or in Musca domestica with 918 mutation  or recently in An. gambiae with the N1575Y mutation . It cannot be excluded that this could be the case in An. funestus s.s..
It is known that some resistance alleles could occur at a low frequency in field populations and gradually increased for many years without been detected, and then when a tipping point is reached, the frequency of that allele could significantly increase in the population in a short period of time. This phenomenon was observed in Ae. aegypti in Mexico  and also in An. coluzzii (previously An. gambiae M form) for which very low frequencies of the 1014F kdr allele (<2%) were detected before 2000 in Vallée du Kou in Burkina Faso  before a significant increase to 75% in ten years . Because it is important to detect resistance allele when it still at low frequency in vector populations in order to facilitate the implementation of suitable resistance management strategies, it is necessary to keep on monitoring for the presence of kdr mutations in An. funestus across the continent. The sequencing of the full-length of the VGSC gene will help to identify such resistance mutation and could lead to the design of suitable diagnostic tool to detect and monitor the evolution of this resistance in natural populations of this vector.
The present study aimed at establishing whether beside the predominant metabolic resistance mechanism, the target site resistance through knockdown resistance was also playing a role in pyrethroid resistance in An. funestus s.s.. To elucidate this question we cloned and sequenced the full length 6 kb of the VGSC gene in several populations of An. funestus across the continent, detected potential kdr resistance mutations and assessed their correlation with pyrethroid or DDT resistance. Furthermore, we assessed whether possible selection was acting on this gene and established pattern of genetic structure between African An. funestus s.s. populations based on VGSC diversity profiles.
Blood-fed female An. funestus adults resting indoor were collected in houses between 06.00 and 12.00 AM in several countries as previously reported. These locations include Mbinga (Katete District) (14° 11′ 0″ S, 31° 52′ 0″ E) in eastern Zambia in October 2010 ; Chikwawa District (0° 45′ N, 34° 5′E) in southern Malawi in July 2009 and April 2010 ; Tihuquine (Chokwe District) (24° 33′ 37″ S, 33° 1′ 20″ E) in southern Mozambique in January 2009 (Cuamba et al. 2010); Tororo in eastern Uganda (0°41′N, 34°10′E) ; Pahou (6° 23′ N, 2° 13′E) in southern Benin . The collection method and the rearing were done as described previously . Briefly, F1 adults were generated from field collected female mosquitoes using the forced-egg laying method  and were randomly mixed in cages for subsequent experiments.
All females used for individual oviposition were morphologically identified as belonging to the An. funestus group according to the key of Gillies and Coetzee . Genomic DNA was extracted using the Livak protocol . A PCR was carried out using the protocol of Koekemoer et al.  to confirm that all females that laid eggs were An. funestus s.s..
Genotyping of 1014 kdr coding position in field An. funestus
A pyrosequencing assay was used to genotype the VGSC 1014 coding position in a set of 100 field collected adult female An. funestus from southern Africa (Zambia, Mozambique, Malawi), from West Africa (Benin) and in East Africa (Uganda) to check the possible presence of the L1014F kdr mutation in this species. The pyrosequencing reactions were carried out according to Wondji et al. [40, 41] with a design allowing the simultaneous detection of both possible L1014F and L1014S kdr mutations reported in An. gambiae. Details of the pyrosequencing assay and the primers used (Additional file 1: Table S1) have been previously published .
Sequence analysis of exon19 in Uganda mosquitoes
A genomic fragment spanning exon 19 and the intron between exon 19 and 20 was amplified and directly sequenced for 9 resistant and 4 permethrin susceptible mosquitoes from Tororo. The genomic fragment was amplified using primers KdrEx19F and KdrEx19R (Additional file 1: Table S2). The PCR was carried out using 10 pmol of each primers and 30 ng of genomic DNA as template in 25 μl reactions containing 1X Kapa Taq buffer, 0.2 mM dNTPs, 1.5 mM MgCl2, 1 U Kapa Taq (Kapa Biosystems, Wilmington, MA, USA). The cycle parameters were: 1 cycle at 95 °C for 5 min; 35 cycles of 94 °C for 30s, 57 °C for 30s and elongation at 72 °C for 1 min; followed by 1 cycle at 72 °C for 10 min. Sequences were aligned using ClustalW  while haplotypes reconstruction and polymorphism analysis were done using DnaSP v5.10 . The phylogenetic Neighbour-joining trees were constructed using MEGA 6.0 .
cDNA preparation and full length amplification of voltage-gated sodium channel gene (VGSC)
RNA was extracted from three batches of ten 2- to 5-day-old An. funestus females alive after exposure to 0.75% permethrin (Resistant, R) and unexposed mosquitoes samples from Malawi, Mozambique, Zambia and Benin and from the susceptible FANG strain (Susceptible, S) using the Picopure RNA Isolation Kit (Thermo Fisher Scientific, Waltham, MA, USA). The quantity and quality of the extracted RNA were assessed using a NanoDrop ND1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and Bioanalyzer (Agilent, Santa Clara, CA, USA), respectively. 1 μg of total RNA from each of the three biological replicates from the resistant populations (R) and FANG (S) was used as the template for cDNA synthesis using Superscript III (Invitrogen, Carlsbad, CA, USA) with oligo-dT20 and RNase H according to the manufacturer’s instructions.
The full length of the Voltage Gated Sodium Channel gene (VGSC) was amplified using cDNA from permethrin resistant mosquito samples from Malawi, Mozambique, Zambia and Benin and from the susceptible FANG strain. Amplification was performed using the Phusion High-Fidelity DNA Polymerase (Thermo Fisher Scientific, Waltham, MA, USA) in a 25 μl reaction containing 10 pmol of each primers, 30 ng of cDNA as template, 1X Phusion® HF Buffer, 0.2 mM dNTPs, 1.5 mM MgCl2, 1 U Phusion Polymerase. The PCR conditions were as follows: 1 cycle at 95 °C for 5 min; 35 cycles of 94 °C for 20s, 57 °C for 30s and elongation at 72 °C for 2 min 30s; and 1 cycle at 72 °C for 5 min. The primers used are listed in Additional file 1: Table S1. The PCR products were purified using the QIAquick PCR Purification Kit (Qiagen, Valencia, CA, USA) and directly sequenced on both strands using a set of primers listed in Additional file 1: Table S1. To confirm mutations obtained after direct sequencing and to resolve the issue of overlapping sequencing peaks due to the presence of alternative splicing positions a portion of VGSC from exon 1 to exon 20 was cloned into the pJET1.2/blunt cloning vector using the CloneJET™ PCR Cloning Kit (Thermo Fisher Scientific, Waltham, MA, USA). Positive clones were purified by QIAprep® Miniprep (Qiagen, Valencia, CA, USA) and sequenced on both strands using a set of primers listed in Additional file 1: Table S2. Polymorphic positions were detected through manual analysis of sequence traces using BioEdit and as sequence differences in multiple alignments using ClustalW .
Analysis of alternative splicing of exon19 with intron amplification
An attempt was made to see if there was any association (difference in frequency) between the distribution of alternative splice version of exon 19 and susceptibility to pyrethroid as the FANG mosquitoes tended to only have the 19d form while the 19c and d were both present in resistant field samples. To assess this, a portion of genomic DNA spanning the exon 19 to 20 including the intron was cloned and sequenced using primers Exon18kdrF and KdrEx20R (Additional file 1: Table S2) between resistant and susceptible mosquitoes from Malawi, Benin and Uganda.
Pyrosequencing of candidate kdr mutations
Pyrosequencing assays were designed to genotype field populations for the candidate kdr mutations detected in respective countries in order to assess their possible involvement in pyrethroid resistance in Uganda, Malawi, Mozambique and Zambia. The pyrosequencing reactions were performed as previously described [36, 41, 45]. Briefly, three sequence specific primers (Additional file 1 Table S1) designed for each mutation using the software provided by Pyrosequencing AB (http://www.pyrosequencing.com) were used for amplification and genotyping of each mutation. Target DNA fragments for each mutation were first PCR-amplified in a reaction containing 10 pmol each of forward and biotinylated reverse primer with the same conditions as described previously .
Genetic differentiation of An. funestus across Africa based on a portion of the VGSC gene
A portion of the VGSC spanning intron 19 and exon 20 has previously been sequenced in Uganda, Kenya, Benin and Malawi in separate studies [8, 12, 13]. These sequences were used in the present study to assess the patterns of genetic differentiation of An. funestus across Africa based on the VGSC gene by estimating the levels of pair-wise genetic differentiation between these populations. Analyses were performed as implemented in dnaSP 5.10 using the K ST statistic . The significance of the K ST estimates was assessed by permutation of subpopulation identities and re-calculating K ST 10,000 times.
Phylogenetic tree of VGSC haplotypes
A maximum likelihood phylogenetic tree was constructed for the VGSC haplotypes from Uganda, Kenyan, Benin and Malawi samples using MEGA 6.0 . The best-fit substitution model was firstly assessed based on the Bayesian Information Criterion (BIC) and used to generate the maximum likelihood tree as implemented in MEGA 6.0  with 500 bootstrap replications to assess the robustness of the tree. A phylogenetic tree was also built with pairwise N ST genetic distances between samples using MEGA 6.0.
Pyrosequencing genotyping of L1014 across Africa
The pyrosequencing assay was successfully used to genotype the kdr 1014 coding position in a set of 100 field collected female adults An. funestus from southern Africa (Zambia, Mozambique, Malawi), from West Africa (Benin). Only the TTA codon for Leucine (L1014) (Additional file 2: Figure S1A) was detected indicating the absence of the 1014F or 1014S kdr mutations common in An. gambiae. However, in Uganda, two field mosquitoes were found to be heterozygotes TTA/T (Additional file 2: Figure S1B) suggesting the possible presence of the 1014F mutation in the An. funestus population of this country. Further steps were taken to validate this heterozygote status by amplifying a portion of exon 20 spanning the 1014 codon in these two mosquitoes. However, analysis of the sequencing traces confirm the absence of the 1014F mutation as both mosquitoes were TTA homozygotes. In addition, pyrosequencing genotyping of the 1014 position was further performed in 80 F1 mosquitoes which was permethrin resistant and 80 susceptible individuals from Tororo in Uganda; but no mutation was detected. Therefore, the two heterozygote TTA/T could be the result of PCR contamination from An. gambiae or induced by pyrosequencing errors.
Sequence analysis of exon19 in Uganda mosquitoes
Because kdr mutation could occur in other exons such as exon 19 where kdr mutations have previously been detected in other species , the Ugandan population was also assessed for the presence of polymorphism possibly associated with permethrin resistance in exon 19. Nine permethrin resistant and four susceptible mosquitoes from Tororo in Uganda were successfully sequenced for an 841 bp portion of the VGSC gene spanning the Exon 19 and part of intron 19. In total 24 substitution sites were detected with 13 haplotypes. Only 3 substitutions were detected on a 101 bp portion of exon19 with one, detected in a single resistant mosquito, leading to an amino acid change with the replacement of a leucine in codon 928 by an arginine (L928R) (Fig. 1a). Additionally, significant estimates (P < 0.05) of Tajima D (−1.96) and Fu and Li D’ (−2.89) were observed for the total sample indicating possible signature of positive selection on this fragment of VGSC in Tororo. However, the pyrosequencing genotyping of the L928R mutation in 50 permethrin resistant and 50 susceptible only successfully amplified the wild CTG genotype for L928 (Fig. 1b). Furthermore, construction of the maximum likelihood tree of the haplotypes did not detect a pattern of haplotype clustering according to phenotype (Fig. 1c).
Sequencing of full-length VGSC gene to detect potential kdr mutations
The detection of the L928R mutations and the previously reported association between VGSC haplotypes and pyrethroid resistance in Malawi  suggested that other potential kdr mutation could be identified in other exons of the VGSC gene. Therefore, a full-length cDNA of this gene was cloned and sequenced in field populations from Mozambique, Uganda and Benin. The full length of VGSC was successfully amplified in An. funestus with a size of 6276 bp corresponding to 2092 amino acids in a Mozambique clone (Acc No. KY499806), 6183 bp and 2061 amino acids in a clone from Uganda (Acc No. KY499807) and 6351 bp and 2117 amino acids in a clone from Benin (Acc No. KY499805). These complete sequences were compared to that of the FUMOZ strain obtained from the recently sequenced genome of An. funestus (AFUN000494-RA) of 6360 bp and 2120 amino acids obtained from VectorBase. These deduced amino acid sequences of An. funestus share 93% similarity with that of An. gambiae (AGAP004707-RA), 94% with that of An. minimus (AMIN007480-RA) and 91% with that of Aedes aegypti (AAEL006019-RA).
Alignment of the four An. funestus VGSC amino acid sequences revealed several mutations (Fig. 2) probably a result of alternative splicing events in this gene as previously reported in An. gambiae (Davies et al. 2007). These are characterised by insertion or deletion of portion of the gene such as the 33 bp in exon 2 detected in Mozambique encoding the “SDFGRKKKKKE” amino acid stretch. In exon 5, a 129 bp encoding domain I (S2-S3) is observed in the cDNA from Uganda (Fig. 2). Additionally, six amino acid changes were also observed with the majority being replacement between amino acids of same classes, such as I123V mutation in Mozambique. However, some amino acid changes such as the P625S (ANLG P/S RHSSY) observed only in FUMOZ could impact the activity of the protein.
Furthermore, to detect potential mutations associated with pyrethroid resistance, pools of cDNA from permethrin resistant mosquitoes from Mozambique, Malawi, Zambia, Benin, Uganda and from the susceptible FANG laboratory strain were directly sequenced. Analysis of these sequences confirmed the absence of the 1014 mutation in An. funestus. A main variation was observed at exon 19 with several polymorphisms observed in all samples but not in the susceptible FANG (Fig. 3a) reflecting the presence of both mutually exclusive exon 19c and 19d as previously reported in Culex quinquefasciatus . The FANG pooled cDNA sequence was equivalent to exon 19d while the field resistant samples exhibited both exon 19c and 19d. Among the many overlapping peaks were the possibility of two potential mutations which could lead to a L913F close to the I915M super kdr and also the L932F mutation (Fig. 3a) shown to confer pyrethroid resistance in Culex quinquefasciatus (Davies et al. 2007).
Because potential resistance mutations are reported in exon 19c and d in Culex quinquefasciatus , further sequencing of 26 individual clones of the fragment spanning exon 1 to 20 was carried out in Mozambique, Malawi, Zambia and Benin. This revealed a predominance of exon 19c (21 out of 26) over exon 19d (5 out of 26). The L913F and L932F mutations were not observed from the cloned samples and this was further validated by the pyrosequencing genotyping in 50 permethrin resistant and 50 susceptible individuals from Malawi and Zambia where only the wild type alleles were detected (Fig. 3b). Additionally, none of the other mutations observed in Cx quinquefasciatus (I915M in 19c and I936V in 19d) was detected. However, two amino acid changes were detected in Mozambique samples (each in a single clone), a T-to-C mutation at codon 802 leading to L802S substitution (TTA to TCA) in exon 18 (Fig. 4a and b) and a C-to-T at codon 1008 leading to T1008 M replacement (ACG to ATG) at exon 20 (Fig. 4c and d) which is located at the vicinity of the common L1014F mutation [LAT (1008) VVIGNL (1014) V]. Again, the pyrosequencing genotyping in 50 permethrin resistant and 50 susceptible from Malawi and Mozambique only detected the wild type alleles (Fig. 4e and f) suggesting that these mutations are probably just present at a very low frequency in the population and play no or only very minor role in the observed pyrethroid resistance The L802S and T1008 M mutation were also not found in the FUMOZ strain. One additional amino acid mutation, I818N, was also detected in Benin in a single clone but correlation with resistance was not assessed.
Analysis of alternative splicing of exon19 with intron amplification
Attempt was made to see if there was any association (difference in frequency) between the distribution of alternative splice version of exon 19 and susceptibility to pyrethroid as the FANG mosquitoes tended to only have one form while the 19c and d were both present in field resistant samples. To assess this, a portion of genomic DNA spanning the exon 18 to 20 including the intron was cloned and sequenced between resistant and susceptible mosquitoes from Malawi, Benin and Uganda. However, cloning of a genomic portion of exon18 to exon 20 between individual permethrin resistant and susceptible mosquitoes revealed that both alternative splice exons were equally present in both phenotypes with apparently no association with pyrethroid resistance.
Population structure of VGSC gene across African populations of An. funestus
Samples from 4 countries were used to assess patterns of genetic variability of the VGSC gene across African populations of An. funestus. These countries are Benin in West Africa, Uganda and Kenya in East Africa and Malawi in Southern Africa.
Haplotype distribution of VGSC across Africa
A total of 41 haplotypes were observed across Africa for the 875 bp fragment of VGSC spanning intron 19 and exon 20 with the lowest haplotype diversity observed in Malawi (Table 1) where a possible association between VGSC diversity and pyrethroid resistance has previously been observed . Overall, the VGSC exhibits a high polymorphism level for all the countries except in Malawi. This high polymorphism is shown first by the fact that a majority of haplotypes occurred as singletons (25/41), although a major haplotype (H4; 21/110) (Fig. 5a) is observed although only shared by the West and East African countries but not by Malawi. Second, the high number of mutational steps between haplotypes as shown on the TCS haplotype network further highlights the high polymorphism of the VGSC gene in these populations (Fig. 5b). Thirdly, there is no indication of selection acting on these populations as shown by the values of the selection test from Tajima (D) and Fu and Li (D*). However, a significant difference was observed between Malawi and other populations, as contrary to negative test values observed in other regions, suggesting an excess of low frequency polymorphisms (in line with high level of singleton haplotypes), Malawi exhibited positive and significant values (for D*) suggesting that this gene is under selection in this southern African populations in line with low haplotype diversity observed.
Population substructure at VGSC across Africa
Analysis of the genetic structure pattern at the VGSC fragment further supported the difference in genetic diversity observed between the Malawian populations and the West and East Africa populations. The construction of a maximum likelihood phylogenetic tree of VGSC sequences revealed that haplotypes from Malawi always cluster together separately from other countries (Fig. 5c). This difference is further supporter by the consistently high genetic differentiation estimates observed between Malawi and the other 3 countries (0.13 < KST < 0.17) whereas, the three West and East African populations exhibited low levels of genetic differentiation (0.035 < KST < 0.057) (Table 2) and cluster together on the Neighbor-joining tree of the genetic distances (Fig. 5d). This pattern of gene flow correlates with the genetic structure profiles of An. funestus populations across Africa also obtained with the DDT resistance gene GSTe2  or with microsatellite markers [50, 51].
Elucidating mechanisms of insecticide resistance is crucial in understanding the evolution and dynamic of resistance in field populations of vectors and to help design suitable resistance management strategies. This study investigated the potential role of the knockdown resistance mechanism in the increasingly reported pyrethroid and DDT resistance in the major malaria vector An. funestus across Africa.
Role of knockdown resistance is different to that observed in An. gambiae s.l.
This study has confirmed previous reports that the role of knockdown resistance in pyrethroid/DDT resistance in An. funestus is different from that observed in other major malaria vectors from An. gambiae complex, notably An. gambiae s.s., An. coluzzii and An. arabiensis. Indeed, despite an extensive effort to detect pyrethroids/DDT kdr mutations, the common L1014F mutation currently widespread in An. gambiae s.s. across Africa, was completely absent in An. funestus as previously observed [8, 12, 14, 18]. The absence of the 1014F mutation in An. funestus populations could also explain the lack of cross-resistance between pyrethroids and DDT in some regions such as in southern Africa where populations are highly resistant to pyrethroids but susceptible to DDT [5, 6] although moderate DDT resistance has recently been reported in this region . This difference between An. funestus and An. gambiae is also observed for the other target-site resistance mechanisms in the acetylcholinesterase gene where the common G119S mutation conferring carbamates and organophosphates resistance in An. gambiae and other mosquitoes such as Cx quinquefasciatus (Weill et al. 2004) is completely absent in An funestus [5, 8] but rather the N485I mutation is found in An. funestus conferring bendiocarb resistance . Furthermore, synergist assays performed with An. funestus across Africa (Riveron et al. 2014; Mulamba et al. 2014) also suggest that knockdown resistance is not a major mechanism contrary to what is observed in An. gambiae or Cx quinquesfaciatus and even Ae. aegypti highlighting the fact that one should not assume that mechanisms present in one Anopheles species apply to all. However, since most samples in this study are from 2009 and 2010, further studies should be performed with recent samples to assess whether the situation remains the same across Africa.
Metabolic resistance remains the predominant pyrethroid resistance mechanism in An. funestus across Africa
The absence of the 1014F/S mutation and other confirmed kdr mutations in this study supports the predominant role that metabolic resistance plays in the pyrethroid resistance in An. funestus as shown previously [23, 25, 26]. This is in contrast to other mosquito species such as An. gambiae for which knockdown resistance plays a significant role . However, despite the absence of the 1014F mutation in An. funestus, the detection of several amino acid replacements in the sodium channel gene in this species across Africa suggests that one cannot rule out a future role of knockdown resistance in pyrethroids and DDT resistance in this species. Indeed, the detection of the two amino acid changes L802S and T1008 M in Mozambique and that of I877L, V881 L and A1007S in Cameroon  or the F1021C in Uganda , albeit always at very low level suggests that such kdr resistance could occur in An. funestus but in other codons than 1014 as reported in Aedes aegypti with the mutations at the 1011, 1016 and 1534 codons [31, 32]. It is necessary to continue monitoring the frequency of the L802S and T1008 M and other mutations in field populations to assess their possible role in pyrethroid resistance. This is particularly important when taking into consideration the concept of a “tipping-point” in the evolution of a mutation in a population under selection as most Anopheles populations are currently with the scale up of insecticide-based control intervention throughout Africa. It is know that a resistance mutation can be present at low but gradually increasing frequency in a mosquito population for many years without being detected . However, once the tipping-point is reached after continuous selection pressure, the frequency of that resistance mutation may shoot up rapidly as seen with the rapid spread of the 1014F mutation in An. gambiae population across Africa  or for the kdr mutation in Ae. aegypti populations in Mexico . On the other hand, it is also possible that the observed amino acid changes have no association with knockdown resistance and just represent the natural polymorphism of the gene. Future monitoring will help establish their role.
Genetic diversity of VGSC supports a limited role of knockdown resistance and suggests restriction to gene flow
The high genetic diversity observed for the VGSC fragment across Africa further supports that knockdown resistance currently plays no role in the pyrethroid/DDT resistance in An. funestus. This lack of selection is supported by the high number of substitutions and high haplotypes number observed throughout Africa with most of the haplotypes been singletons. The absence of selection on this gene in An. funestus is in marked contrast to the extensive selective sweep reported around the VGSC in An. gambiae with significant reduced genetic and haplotype diversity . The restricted gene flow observed in this study between southern Africa and other regions is in line with previous analyses of the genetic diversity of different insecticide resistance genes such as the DDT resistance gene GSTe2  or the CYP6P9a and CYP6P9b pyrethroid resistance genes [22, 51, 54]. It is also in line with patterns of differentiation obtained with microsatellite markers [22, 50] further supporting a strong genetic differentiation between the southern African populations to those from East and West Africa. This further supports the presence of barriers to gene flow between An. funestus populations across Africa as previously suggested [22, 50]. Such restriction to gene flow could limit the spread of resistance mutations between populations of this species from different regions of Africa or lead to heterogeneous biological and genetic profile of this vectors across the continent.
Investigation of knockdown resistance in the malaria vector An. funestus did not establish a major role for this resistance mechanism in the pyrethroid or DDT resistance in this species across Africa despite the detection of some non-synonymous mutations. Nevertheless, these mutations in the VGSC gene should be constantly monitored to detect any increased frequency that could be associated with selection of this resistance mechanisms in An. funestus populations.
WHO: Global Plan for Insecticide Resistance Management (GPIRM). Geneva, Switzerland: World Health Organization; 2012.
Brooke BD, Kloke G, Hunt RH, Koekemoer LL, Temu EA, Taylor ME, Small G, Hemingway J, Coetzee M. Bioassay and biochemical analyses of insecticide resistance in southern African Anopheles Funestus (Diptera: Culicidae). Bull Entomol Res. 2001;91(4):265–72.
Hargreaves K, Koekemoer LL, Brooke BD, Hunt RH, Mthembu J, Coetzee M. Anopheles Funestus resistant to pyrethroid insecticides in South Africa. Med Vet Entomol. 2000;14(2):181–9.
Chanda E, Hemingway J, Kleinschmidt I, Rehman AM, Ramdeen V, Phiri FN, Coetzer S, Mthembu D, Shinondo CJ, Chizema-Kawesha E, et al. Insecticide resistance and the future of malaria control in zambia. PLoS One. 2011;6(9):e24336.
Cuamba N, Morgan JC, Irving H, Steven A, Wondji CS. High level of pyrethroid resistance in an Anopheles Funestus population of the Chokwe District in Mozambique. PLoS One. 2010;5(6):e11010.
Hunt R, Edwardes M, Coetzee M. Pyrethroid resistance in southern African Anopheles Funestus extends to Likoma Island in Lake Malawi. Parasit Vectors. 2010;3:122.
Vezenegho SB, Chiphwanya J, Hunt RH, Coetzee M, Bass C, Koekemoer LL. Characterization of the Anopheles Funestus group, including Anopheles Funestus-like, from northern Malawi. Trans R Soc Trop Med Hyg. 2013;107(12):753–62.
Wondji CS, Coleman M, Kleinschmidt I, Mzilahowa T, Irving H, Ndula M, Rehman A, Morgan J, Barnes KG, Hemingway J. Impact of pyrethroid resistance on operational malaria control in Malawi. Proc Natl Acad Sci U S A. 2012;109(47):19063–70.
Riveron JM, Chiumia M, Menze BD, Barnes KG, Irving H, Ibrahim SS, Weedall GD, Mzilahowa T, Wondji CS. Rise of multiple insecticide resistance in Anopheles Funestus in Malawi: a major concern for malaria vector control. Malar J. 2015;14:344.
Kawada H, Dida GO, Ohashi K, Komagata O, Kasai S, Tomita T, et al. Multimodal Pyrethroid Resistance in Malaria Vectors, Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus s.s. in Western Kenya. PLoS ONE 2011;6(8): e22574.
Morgan JC, Irving H, Okedi LM, Steven A, Wondji CS. Pyrethroid resistance in an Anopheles Funestus population from Uganda. PLoS One. 2010;5(7):e11872.
Mulamba C, Riveron JM, Ibrahim SS, Irving H, Barnes KG, Mukwaya LG, Birungi J, Wondji CS. Widespread pyrethroid and DDT resistance in the major malaria vector Anopheles Funestus in East Africa is driven by metabolic resistance mechanisms. PLoS One. 2014;9(10):e110058.
Djouaka R, Irving H, Tukur Z, Wondji CS. Exploring mechanisms of multiple insecticide resistance in a population of the malaria vector Anopheles Funestus in Benin. PLoS One. 2011;6(11):e27760.
Okoye PN, Brooke BD, Koekemoer LL, Hunt RH, Coetzee M. Characterisation of DDT, pyrethroid and carbamate resistance in Anopheles Funestus from Obuasi, Ghana. Trans R Soc Trop Med Hyg. 2008;102(6):591–8.
Djouaka R, Riveron JM, Yessoufou A, Tchigossou G, Akoton R, Irving H, Djegbe I, Moutairou K, Adeoti R, Tamo M, et al. Multiple insecticide resistance in an infected population of the malaria vector Anopheles Funestus in Benin. Parasit Vectors. 2016;9:453.
Djouaka RJ, Atoyebi SM, Tchigossou GM, Riveron JM, Irving H, Akoton R, Kusimo MO, Bakare AA, Wondji CS. Evidence of a multiple insecticide resistance in the malaria vector Anopheles Funestus in south West Nigeria. Malar J. 2016;15(1):565.
Samb B, Konate L, Irving H, Riveron JM, Dia I, Faye O, Wondji CS. Investigating molecular basis of lambda-cyhalothrin resistance in an Anopheles Funestus population from Senegal. Parasit Vectors. 2016;9(1):449.
Menze BD, Riveron JM, Ibrahim SS, Irving H, Antonio-Nkondjio C, Awono-Ambene PH, Wondji CS. Multiple insecticide resistance in the malaria vector Anopheles Funestus from northern Cameroon is mediated by metabolic resistance alongside potential target site insensitivity mutations. PLoS One. 2016;11(10):e0163261.
Wondji CS, Dabire RK, Tukur Z, Irving H, Djouaka R, Morgan JC. Identification and distribution of a GABA receptor mutation conferring dieldrin resistance in the malaria vector Anopheles Funestus in Africa. Insect Biochem Mol Biol. 2011;41(7):484–91.
Coetzee M, Koekemoer LL. Molecular systematics and insecticide resistance in the major African malaria vector Anopheles Funestus. Annu Rev Entomol. 2013;58:393–412.
Wondji CS, Irving H, Morgan J, Lobo NF, Collins FH, Hunt RH, Coetzee M, Hemingway J, Ranson H. Two duplicated P450 genes are associated with pyrethroid resistance in Anopheles Funestus, a major malaria vector. Genome Res. 2009;19(3):452–9.
Barnes KG, Irving H, Chiumia M, Mzilahowa T, Coleman M, Hemingway J, Wondji CS. Restriction to gene flow is associated with changes in the molecular basis of pyrethroid resistance in the malaria vector Anopheles Funestus. Proc Natl Acad Sci U S A. 2017;114(2):286–91.
Riveron JM, Irving H, Ndula M, Barnes KG, Ibrahim SS, Paine MJ, Wondji CS. Directionally selected cytochrome P450 alleles are driving the spread of pyrethroid resistance in the major malaria vector Anopheles Funestus. Proc Natl Acad Sci U S A. 2013;110(1):252–7.
Thomsen EK, Strode C, Hemmings K, Hughes AJ, Chanda E, Musapa M, Kamuliwo M, Phiri FN, Muzia L, Chanda J, et al. Underpinning sustainable vector control through informed insecticide resistance management. PLoS One. 2014;9(6):e99822.
Amenya DA, Naguran R, Lo TC, Ranson H, Spillings BL, Wood OR, Brooke BD, Coetzee M, Koekemoer LL. Over expression of a cytochrome P450 (CYP6P9) in a major African malaria vector, anopheles Funestus, resistant to pyrethroids. Insect Mol Biol. 2008;17(1):19–25.
Irving H, Riveron JM, Ibrahim SS, Lobo NF, Wondji CS. Positional cloning of rp2 QTL associates the P450 genes CYP6Z1, CYP6Z3 and CYP6M7 with pyrethroid resistance in the malaria vector Anopheles Funestus. Heredity (Edinb). 2012;109(6):383–92.
Riveron JM, Ibrahim SS, Chanda E, Mzilahowa T, Cuamba N, Irving H, Barnes KG, Ndula M, Wondji CS. The highly polymorphic CYP6M7 cytochrome P450 gene partners with the directionally selected CYP6P9a and CYP6P9b genes to expand the pyrethroid resistance front in the malaria vector Anopheles Funestus in Africa. BMC Genomics. 2014;15(1):817.
Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, Guillet P, Pasteur N, Pauron D. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles Gambiae s.S. Insect Mol Biol. 1998;7(2):179–84.
Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH. Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles Gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000;9(5):491–7.
Ranson H, N'Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V. Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol. 2011;27(2):91–8.
Brengues C, Hawkes NJ, Chandre F, McCarroll L, Duchon S, Guillet P, Manguin S, Morgan JC, Hemingway J. 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.
Harris AF, Rajatileka S, Ranson H. Pyrethroid resistance in Aedes Aegypti from grand Cayman. Am J Trop Med Hyg. 2010;83(2):277–84.
Soderlund DM, Knipple DC. The molecular biology of knockdown resistance to pyrethroid insecticides. Insect Biochem Mol Biol. 2003;33(6):563–77.
Jones CM, Liyanapathirana M, Agossa FR, Weetman D, Ranson H, Donnelly MJ, Wilding CS. Footprints of positive selection associated with a mutation (N1575Y) in the voltage-gated sodium channel of Anopheles Gambiae. Proc Natl Acad Sci U S A. 2012;109(17):6614–9.
Diabate A, Baldet T, Chandre F, Guiguemde RT, Brengues C, Guillet P, Hemingway J, Hougard JM. First report of the kdr mutation in Anopheles Gambiae M form from Burkina Faso, west Africa. Parassitologia. 2002;44(3–4):157–8.
Kwiatkowska RM, Platt N, Poupardin R, Irving H, Dabire RK, Mitchell S, Jones CM, Diabate A, Ranson H, Wondji CS. Dissecting the mechanisms responsible for the multiple insecticide resistance phenotype in Anopheles Gambiae s.S., M form, from Vallee du Kou, Burkina Faso. Gene. 2013;519(1):98–106.
Gillies MT, Coetzee M. A supplement to the Anophelinae of Africa south of the Sahara (Afrotropical region), vol. 55. South African Institute for medical research: Johannesburg; 1987.
Livak KJ. Organization and mapping of a sequence on the Drosophila Melanogaster X and Y chromosomes that is transcribed during spermatogenesis. Genetics. 1984;107(4):611–34.
Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles Funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66(6):804–11.
Wondji CS, Priyanka De Silva WA, Hemingway J, Ranson H, Parakrama Karunaratne SH. Characterization of knockdown resistance in DDT- and pyrethroid-resistant Culex Quinquefasciatus populations from Sri Lanka. Tropical Med Int Health. 2008;13(4):548–55.
Wondji CS, Morgan JC, Coetzee M, Hunt R, Steen K, Black WC, Hemingway J, Ranson H. Mapping a quantitative trait locus conferring pyrethroid resistance in the African malaria vector Anopheles Funestus. BMC Genomics. 2007;8:34.
Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22(22):4673–80.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
Witzig C, Parry M, Morgan JC, Irving H, Steven A, Cuamba N, Kerah-Hinzoumbe C, Ranson H, Wondji CS. Genetic mapping identifies a major locus spanning P450 clusters associated with pyrethroid resistance in kdr-free anopheles arabiensis from Chad. Heredity (Edinb). 2013;110(4):389–97.
Ibrahim SS, Manu YA, Tukur Z, Irving H, Wondji CS. High frequency of kdr L1014F is associated with pyrethroid resistance in anopheles coluzzii in Sudan savannah of northern Nigeria. BMC Infect Dis. 2014;14:441.
Hudson RR, Slatkin M, Maddison WP. Estimation of levels of gene flow from DNA sequence data. Genetics. 1992;132(2):583–9.
Davies TG, Field LM, Usherwood PN, Williamson MS. A comparative study of voltage-gated sodium channels in the Insecta: implications for pyrethroid resistance in Anopheline and other Neopteran species. Insect Mol Biol. 2007;16(3):361–75.
Riveron JM, Yunta C, Ibrahim SS, Djouaka R, Irving H, Menze BD, Ismail HM, Hemingway J, Ranson H, Albert A, et al. A single mutation in the GSTe2 gene allows tracking of metabolically-based insecticide resistance in a major malaria vector. Genome Biol. 2014;15(2):R27.
Michel AP, Ingrasci MJ, Schemerhorn BJ, Kern M, Le Goff G, Coetzee M, Elissa N, Fontenille D, Vulule J, Lehmann T, et al. Rangewide population genetic structure of the African malaria vector Anopheles Funestus. Mol Ecol. 2005;14(14):4235–48.
Barnes KG, Weedall GD, Ndula M, Irving H, Mzihalowa T, Hemingway J, Wondji CS. Genomic footprints of selective sweeps from metabolic resistance to Pyrethroids in African malaria vectors are driven by scale up of insecticide-based vector control. PLoS Genet. 2017;13(2):e1006539.
Ibrahim SS, Ndula M, Riveron JM, Irving H, Wondji CS. The P450 CYP6Z1 confers carbamate/pyrethroid cross-resistance in a major African malaria vector beside a novel carbamate-insensitive N485I acetylcholinesterase-1 mutation. Mol Ecol. 2016; In Press
Clarkson CS, Weetman D, Essandoh J, Yawson AE, Maslen G, Manske M, Field SG, Webster M, Antao T, MacInnis B, et al. Adaptive introgression between anopheles sibling species eliminates a major genomic island but not reproductive isolation. Nat Commun. 2014;5:4248.
Ibrahim SS, Riveron JM, Bibby J, Irving H, Yunta C, Paine MJ, Wondji CS. Allelic variation of cytochrome P450s drives resistance to Bednet insecticides in a major malaria vector. PLoS Genet. 2015;11(10):e1005618.
This work was supported by a Wellcome Trust Research Career Development Fellowship (083515/Z/07/Z) and a Wellcome Trust Senior Research Fellowship in Biomedical Sciences to CSW (101,893/Z/13/Z). The authors thanks Drs Jacob Riveron and Sulaiman Ibrahim for helpful comments.
Availability of data and materials
Most of the data supporting these findings are contained within the manuscript. However full sequencing data can be retrieved from Genbank and further details obtained by contacting the corresponding author.
The DNA sequences reported in this paper have been deposited in the GenBank database (accession numbers: GenBank KY499805- KY499833).
Ethics approval and consent to participate
Consent was obtained from village chiefs and from house owners before mosquito collection.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Pyrosequencing primers information for the Voltage-Gated Sodium channel gene mutations. Table S2. List of primers used to sequence full VGSC. (DOCX 15 kb)
Pyrosequencing genotyping of L1014 position in An. funestus. A) is the chromatogram showing the wild TTA genotype for L1014 whereas B) is the heterozygote TTA/T for L1014/1014F observed at 2 mosquitoes but not confirmed by sequencing. (TIFF 33 kb)