Molecular Ecology of Pyrethroid Knockdown Resistance in Culex pipiens pallens Mosquitoes

Pyrethroid insecticides have been extensively used in China and worldwide for public health pest control. Accurate resistance monitoring is essential to guide the rational use of insecticides and resistance management. Here we examined the nucleotide diversity of the para-sodium channel gene, which confers knockdown resistance (kdr) in Culex pipiens pallens mosquitoes in China. The sequence analysis of the para-sodium channel gene identified L1014F and L1014S mutations. We developed and validated allele-specific PCR and the real-time TaqMan methods for resistance diagnosis. The real-time TaqMan method is more superior to the allele-specific PCR method as evidenced by higher amplification rate and better sensitivity and specificity. Significant positive correlation between kdr allele frequency and bioassay-based resistance phenotype demonstrates that the frequency of L1014F and L1014S mutations in the kdr gene can be used as a molecular marker for deltamethrin resistance monitoring in natural Cx. pipiens pallens populations in the East China region. The laboratory selection experiment found that L1014F mutation frequency, but not L1014S mutation, responded to deltamethrin selection, suggesting that the L1014F mutation is the key mutation conferring resistance to deltamethrin. High L1014F mutation frequency detected in six populations of Cx. pipens pallens suggests high prevalence of pyrethroid resistance in Eastern China, calling for further surveys to map the resistance in China and for investigating alternative mosquito control strategies.     

First report of L1014F kdr mutation in Culex quinquefasciatus in Mexico

The L1014F mutation in the voltage‐sodium channel gene has been associated with resistance to DDT and pyrethroids in various arthropod species including mosquitoes. We determined the frequency of the L1014F kdr mutation in 16 field populations of Culex quinquefasciatus from Northeastern Mexico collected between 2008 and 2013. The L1014F was present in all populations analyzed with the lowest frequency (3.33%) corresponding to the population from Monclova collected in 2012, and the highest frequency (63.63%) from the Monterrey population collected in 2012. The presence of a kdr mutation in populations of Cx. quinquefasciatus from northeastern Mexico provides evidence of pyrethroid resistance. This requires a special attention, considering that pyrethroid‐based insecticides are commonly used in vector‐control campaigns, especially against Aedes aegypti (L.).     

Widespread pyrethroid resistance in Australian diamondback moth, Plutella xylostella (L.), is related to multiple mutations in the para sodium channel gene

Populations of Plutella xylostella, extending over 3800 km in southern Australia, show no genetic structure as assessed by microsatellite markers; yet outbreaks of pyrethroid resistance occur sporadically in cropping areas. Since mutations in the para voltage-gated sodium channel gene have been implicated in pyrethroid resistance, we looked for DNA sequence variation at this target among Australian moths. We found two resistance mutations previously reported for this species (L1014F and T929I), as well as a novel substitution (F1020S). Of the eight possible haplotypes formed by combinations of these three biallelic polymorphisms, only four were found in Australian populations: the wild-type allele (w), the kdr mutation allele (kdr) with only L1014F, the super-kdr-like combination of L1014F and T929I (skdrl), and the crashdown allele with only F1020S (cdr). Comparison of genotype frequencies among survivors of permethrin assays with those from untreated controls identified three resistant genotypes: skdrl homozygotes, cdr homozygotes and the corresponding heterozygote, cdr/skrdl - the heterozygote being at least as resistant as either homozygote. Spatial heterogeneity of allele frequencies was conspicuous, both across the continent and among local collections, consistent with reported spatial heterogeneity of pyrethroid resistance. Further, high resistance samples were sometimes associated with high frequency of cdr, sometimes high frequency of skdrl, or sometimes with a high combined cdr+skdrl frequency. The skdrl and cdr alleles explain a high proportion of the Australia-wide resistance variation. These data add to evidence that nerve insensitivity by mutations in the para-sodium channel gene is a common pyrethroid resistance mechanism in P. xylostella.     

Identification of mutations associated with pyrethroid resistance in the para-type sodium channel of the cat flea, Ctenocephalides felis

Knockdown resistance (kdr) to pyrethroid insecticides is caused by point mutations in the pyrethroid target site, the para-type sodium channel of nerve membranes. This most commonly involves alterations within the domain II (S4–S6) region of the channel protein where five different mutation sites have been identified across a range of insect species. To investigate the incidence of this mechanism in cat fleas, we have cloned and sequenced the IIS4–IIS6 region of the para sodium channel gene from seven laboratory flea strains. Analysis of these sequences revealed two amino acid replacements at residues previously implicated in pyrethroid resistance. One is the ‘common’ kdr mutation, a leucine to phenylalanine substitution (equivalent to L1014F of housefly) reported previously in several other insects. The other is a threonine to valine substitution (equivalent to T929V) and is a novel variant of the T929I mutation first identified in diamondback moth. The L1014F mutation was found at varying frequency in all of the laboratory flea strains, whereas the T929V mutation was found only in the highly resistant Cottontail strain. We have developed rapid PCR-based diagnostic assays for the detection of these mutations in individual cat fleas and used them to show that both L1014F and T929V are common in UK and US flea populations. This survey revealed a significant number of fleas that carry only the V929 allele indicating that co-expression with the F1014 allele is not necessary for flea viability.     

Insecticide resistance in Culex quinquefasciatus mosquitoes after the introduction of insecticide-treated bed nets in Macha, Zambia

Culex quinquefasciatus, an arboviral and filarial vector, is present in high numbers throughout sub-Saharan Africa, and insecticide-resistant populations have been reported worldwide. In order to determine the insecticide resistance status of Cx. quinquefasciatus in Macha, Zambia, adult mosquitoes reared from eggs collected from oviposition traps were tested by bioassay. High levels of resistance to DDT, pyrethroids, malathion, and deltamethrin-treated net material were detected, and molecular assays revealed that the knockdown resistance (kdr) allele was frequent in the Cx. quinquefasciatus population, with 7.0% homozygous for the kdr L1014 allele and 38.5% heterozygous (0.263 kdr frequency). The kdr frequency was significantly higher in mosquitoes that had successfully fed on human hosts, and screening archived specimens revealed that kdr was present at lower frequency prior to the introduction of ITNs, indicating that ITNs might be a selective force in this population. Additionally, metabolic detoxification enzyme activity assays showed upregulated glutathione S-transferases, α-esterases, and β-esterases. Continued monitoring and assessment of the Cx. quinquefasciatus population is necessary to determine levels of resistance.     

High incidence of ace-1 duplicated haplotypes in resistant Culex pipiens mosquitoes from Algeria

The status of genes conferring resistance to organophosphate and carbamate insecticides has been examined in Culex pipiens pipiens mosquitoes sampled in Algeria. Presence of overproduced esterases was sporadic, but acetylcholinesterase-1 resistant alleles were observed in almost all samples. We focused our study on the AChE1 G119S substitution characterized in almost all samples, mostly at the heterozygous state. A genetic test revealed the presence of ace-1 duplication associating a susceptible and a resistant ace-1 copy. Molecular characterization showed a high occurrence of ace-1 duplication with six distinct duplicated alleles out of four samples. The inferred frequency of duplicated allele suggests that it is replacing the single resistant G119S allele. Finally, we discuss the mechanism at the origin of these duplicated haplotypes and their consequences on the management of insecticide resistance.     

Preliminary report of knockdown resistance in Culex pipiens pallens and Aedes koreicus from Korea

Pyrethroid insecticides have been effective and powerful for controlling mosquitoes. However, abuse of these insecticides increases the number of resistant mosquitoes. In this study, Culex pipiens pallens and Aedes koreicus were collected from an artificial reservoir in the vicinity of a populated area in Korea, which is also a migratory bird catchment area. To monitor resistance to pyrethroid insecticides in mosquitoes, genomic DNA from the collected mosquitoes was sequenced for the kdr mutation in the voltage‐gated sodium channel (VGSC) gene. As a result, three samples with homozygous resistance (17.6%) and one with heterozygous resistance (5.9%) were found among 17 Cx. pipiens pallens specimens. One of the samples had a unique sequence at the amplified VGSC region. Of the 15 Ae. koreicus, no insecticide resistant individuals were found. In Korea, this is the first report of kdr genetic traits in Ae. koreicus and Cx. pipiens pallens and of a unique VGSC allele in Cx. pipiens pallens. Further investigation is needed to monitor the kdr resistance of these species in Korea and to determine how the unique sequence found in Cx. pipiens pallens is related to insecticide resistance.     

Multiple origins of knockdown resistance mutations in the Afrotropical mosquito vector Anopheles gambiae

How often insecticide resistance mutations arise in natural insect populations is a fundamental question for understanding the evolution of resistance and also for modeling its spread. Moreover, the development of resistance is regarded as a favored model to study the molecular evolution of adaptive traits. In the malaria vector Anopheles gambiae two point mutations (L1014F and L1014S) in the voltage-gated sodium channel gene, that confer knockdown resistance (kdr) to DDT and pyrethroid insecticides, have been described. In order to determine whether resistance alleles result from single or multiple mutation events, genotyping of the kdr locus and partial sequencing of the upstream intron-1 was performed on a total of 288 A. gambiae S-form collected from 28 localities in 15 countries. Knockdown resistance alleles were found to be widespread in West Africa with co-occurrence of both 1014S and 1014F in West-Central localities. Differences in intron-1 haplotype composition suggest that kdr alleles may have arisen from at least four independent mutation events. Neutrality tests provided evidence for a selective sweep acting on this genomic region, particularly in West Africa. The frequency and distribution of these kdr haplotypes varied geographically, being influenced by an interplay between different mutational occurrences, gene flow and local selection. This has important practical implications for the management and sustainability of malaria vector control programs.     

Inheritance of L1014F and M918T sodium channel mutations associated with pyrethroid resistance in Myzus persicae

Two amino acid substitutions (L1014F and M918T) in the voltage-gated sodium channel confer target-site resistance to pyrethroid insecticides in the peach potato aphid, Myzus persicae. Pyrethroid-resistant and -susceptible M. persicae clones with various combinations of these mutations were crossed under laboratory conditions, and the genotypes of aphid progeny were analysed by direct DNA sequencing of the IIS4-S6 region of the sodium channel gene. Segregation patterns showed that in aphids heterozygous for both L1014F and M918T, both mutations were present in the same resistance allele. Despite these mutations appearing largely recessive in other pest species, such aphids exhibited strong resistance to pyrethroids in leaf-dip bioassays. These results have important implications for the spread and management of pyrethroid resistance in field populations.     

Multiple duplications of the rare ace-1 mutation F290V in Culex pipiens natural populations

Two amino acid substitutions in acetylcholinesterase 1 (AChE1), G119S and F290V, are responsible for resistance to organophosphate and carbamate insecticides in Culex pipiens mosquitoes. These mutations generate very different levels of insensitivity to insecticide inhibitors. We described here a biochemical method that rapidly identifies AChE1 variants (susceptible, G119S and F290V, named S, R and V, respectively) present in individual mosquitoes. We investigated the frequency of AChE1 phenotypes in 41 field samples collected around the Mediterranean Sea. F290V substitution was found only in 15 samples and at low frequency, whereas G119S was highly spread in all samples. However, seven V distinct alleles were identified whereas only one R allele was present. The [V] enzymatic phenotype was never observed alone, and the V allele was always found associated with the susceptible and/or G119S AChE1 ([VS], [VR] or [VRS] phenotypes). Furthermore, we showed the presence of duplicated alleles, associating a susceptible and a V copy of the ace-1 gene, in most individuals analyzed for its presence. Evolutionary forces driving the large number of F290V ace-1 alleles and their low frequency in Mediterranean countries are discussed.     

我国蚊虫体内感染的Wolbachia的wsp基因序列测定与分析

测定了我国尖音库蚊复合组和白纹伊蚊蚊虫体内感染的Wolbachia株的wsp基因序列。核苷酸和氨基酸的同源性及系统关系分析表明,我国尖音库蚊复合组和白纹伊蚊中Wolbachia株的wsp基因序列与Pip组其它株的核苷酸及氨基酸同源性分别为98%～100%和97%～100%, 属B大组Wolbachia中的Pip组。     

Distribution and Frequency of kdr Mutations within Anopheles gambiae s.l. Populations and First Report of the Ace.1G119S Mutation in Anopheles arabiensis from Burkina Faso (West Africa)

An entomological survey was carried out at 15 sites dispersed throughout the three eco-climatic regions of Burkina Faso (West Africa) in order to assess the current distribution and frequency of mutations that confer resistance to insecticides in An. gambiae s.l. populations in the country. Both knockdown (kdr) resistance mutation variants (L1014F and L1014S), that confer resistance to pyrethroid insecticides, were identified concomitant with the ace-1 G119S mutation confirming the presence of multiple resistance mechanisms in the An. gambiae complex in Burkina Faso. Compared to the last survey, the frequency of the L1014F kdr mutation appears to have remained largely stable and relatively high in all species. In contrast, the distribution and frequency of the L1014S mutation has increased significantly in An. gambiae s.l. across much of the country. Furthermore we report, for the first time, the identification of the ace.1 G116S mutation in An. arabiensis populations collected at 8 sites. This mutation, which confers resistance to organophosphate and carbamate insecticides, has been reported previously only in the An. gambiae S and M molecular forms. This finding is significant as organophosphates and carbamates are used in indoor residual sprays (IRS) to control malaria vectors as complementary strategies to the use of pyrethroid impregnated bednets. The occurrence of the three target-site resistance mutations in both An. gambiae molecular forms and now An. arabiensis has significant implications for the control of malaria vector populations in Burkina Faso and for resistance management strategies based on the rotation of insecticides with different modes of action.     

库蚊野生种群扩增等位基因重组现象分析

羧酸酯酶 (carboxylesterases)的过量产生是库蚊Culexpipiens对有机磷杀虫剂 (OP)产生抗性的主要机制。由est 3和est 2组成的酯酶超级基因座 (estersuper locus)的基因扩增是引起酯酶基因扩增的主要遗传学基础。通过淀粉电泳研究了采自广州、佛山、郑州的库蚊野生蚊虫种群 ,发现在这些种群中存在着扩增等位基因重组现象。该现象可能是蚊虫受到杀虫药剂的选择压力、等位基因多样性和等位基因型频率的影响。这将提供一个研究抗性进化的自然模型。     

Three years of insecticide resistance monitoring in Anopheles gambiae in Burkina Faso: resistance on the rise?

ABSTRACT: Background and methods A longitudinal Anopheles gambiae s.l. insecticide-resistance monitoring programme was established in four sentinel sites in Burkina Faso. For three years, between 2008 and 2010, WHO diagnostic dose assays were used to measure the prevalence of resistance to all the major classes of insecticides at the beginning and end of the malaria transmission season. Species identification and genotyping for target site mutations was also performed and the sporozoite rate in adults determined. RESULTS: At the onset of the study, resistance to DDT and pyrethroids was already prevalent in An. gambiae s.l. from the south-west of the country but mosquitoes from the two sites in central Burkina Faso were largely susceptible. Within three years, DDT and permethrin resistance was established in all four sites. Carbamate and organophosphate resistance remains relatively rare and largely confined to the south-western areas although a small number of bendiocarb survivors were found in all sites by the final round of monitoring. The ace-1R target site resistance allele was present in all localities and its frequency exceeded 20% in 2010 in two of the sites. The frequency of the 1014 F kdr mutation increased throughout the three years and by 2010, the frequency of 1014 F in all sites combined was 0.02 in Anopheles arabiensis, 0.56 in An. gambiae M form and 0.96 in An. gambiae S form. This frequency did not differ significantly between the sites. The 1014 S kdr allele was only found in An. arabiensis but its frequency increased significantly throughout the study (P = 0.0003) and in 2010 the 1014 S allele frequency was 0.08 in An. arabiensis. Maximum sporozoite rates (12%) were observed in Soumousso in 2009 and the difference between sites is significant for each year. CONCLUSION: Pyrethroid and DDT resistance is now established in An. gambiae s.l. throughout Burkina Faso. Results from diagnostic dose assays are highly variable within and between rounds of testing, and hence it is important that resistance monitoring is carried out on more than one occasion before decisions on insecticide procurement for vector control are made. The presence of 1014 S in An. gambiae s.l., in addition to 1014 F, is not unexpected given the recent report of 1014 S in Benin but highlights the importance of monitoring for both mutations throughout the continent. Future research must now focus on the impact that this resistance is having on malaria control in Burkina Faso.     

CRISPR/Cas9 modified An. gambiae carrying kdr mutation L1014F functionally validate its contribution in insecticide resistance and combined effect with metabolic enzymes

Insecticide resistance in Anopheles mosquitoes is a major obstacle in maintaining the momentum in reducing the malaria burden; mitigating strategies require improved understanding of the underlying mechanisms. Mutations in the target site of insecticides (the voltage gated sodium channel for the most widely used pyrethroid class) and over-expression of detoxification enzymes are commonly reported, but their relative contribution to phenotypic resistance remain poorly understood. Here we present a genome editing pipeline to introduce single nucleotide polymorphisms in An. gambiae which we have used to study the effect of the classical kdr mutation L1014F (L995F based on An. gambiae numbering), one of the most widely distributed resistance alleles. Introduction of 1014F in an otherwise fully susceptible genetic background increased levels of resistance to all tested pyrethroids and DDT ranging from 9.9-fold for permethrin to >24-fold for DDT. The introduction of the 1014F allele was sufficient to reduce mortality of mosquitoes after exposure to deltamethrin treated bednets, even as the only resistance mechanism present. When 1014F was combined with over-expression of glutathione transferase Gste2, resistance to permethrin increased further demonstrating the critical combined effect between target site resistance and detoxification enzymes in vivo. We also show that mosquitoes carrying the 1014F allele in homozygosity showed fitness disadvantages including increased mortality at the larval stage and a reduction in fecundity and adult longevity, which can have consequences for the strength of selection that will apply to this allele in the field.     

First report on co-occurrence knockdown resistance mutations and susceptibility to beta-cypermethrin in Anopheles sinensis from Jiangsu Province, China

The increasing prevalence of insecticide resistance in Anopheles sinensis, a major vector of malaria in Jiangsu province in eastern China, threatens to compromise the successful use of insecticides in malaria control strategies. It is therefore vital to understand the insecticide resistance status of An. sinensis in the region. This study examined the nucleotide diversity of the para-sodium channel and knockdown resistance (kdr) in five field populations of adult An. sinensis mosquitoes collected in Jiangsu province, identifying the L1014F and L1014C substitutions for the first time. Competitive polymerase chain reaction (PCR) amplification of specific allele (cPASA) and polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) for resistance diagnosis were developed and validated. Comparing the results with direct sequencing revealed that the PCR-RFLP method was more sensitive and specific whereas the cPASA method was more convenient and suitable. The significant positive correlation between kdr allele frequency and bioassay-based resistance phenotype demonstrates that the frequency of L1014F and L1014C substitutions in the kdr gene provides a useful molecular marker for monitoring beta-cypermethrin resistance in natural populations of An. sinensis. Our results point to the L1014F substitution as the key mutation associated with beta-cypermethrin resistance. The high resistance and mutation frequency detected in the five populations also suggest cross-resistance with other pyrethroids may occur in An. sinensis, highlighting the need for further surveys to map insecticide resistance in China and the adoption of a rational management of insecticide application for resistance management and mosquito vector control.     

Whole-genome sequencing to detect mutations associated with resistance to insecticides and Bt proteins in Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a major pest native to the Americas that has recently invaded the Old World. Point mutations in the target-site proteins acetylcholinesterase-1 (ace-1), voltage-gated sodium channel (VGSC) and ryanodine receptor (RyR) have been identified in S. frugiperda as major resistance mechanisms to organophosphate, pyrethroid and diamide insecticides respectively. Mutations in the adenosine triphosphate-binding cassette transporter C2 gene (ABCC2) have also been identified to confer resistance to Cry1F protein. In this study, we applied a whole-genome sequencing (WGS) approach to identify point mutations in the target-site genes in 150 FAW individuals collected from China, Malawi, Uganda and Brazil. This approach revealed three amino acid substitutions (A201S, G227A and F290V) of S. frugiperda ace-1, which are known to be associated with organophosphate resistance. The Brazilian population had all three ace-1 point mutations and the 227A allele (mean frequency = 0.54) was the most common. Populations from China, Malawi and Uganda harbored two of the three ace-1 point mutations (A201S and F290V) with the 290V allele (0.47–0.58) as the dominant allele. Point mutations in VGSC (T929I, L932F and L1014F) and RyR (I4790M and G4946E) were not detected in any of the 150 individuals. A novel 12-bp insertion mutation in exon 15 of the ABCC2 gene was identified in some of the Brazilian individuals but absent in the invasive populations. Our results not only demonstrate robustness of the WGS-based genomic approach for detection of resistance mutations, but also provide insights for improvement of resistance management tactics in S. frugiperda.     

Crossing type variability associated with cytoplasmic incompatibility in Australian populations of the mosquito Culex quinquefasciatus Say

An analysis of cytoplasmic crossing type variation in Australian populations of Culex quinquefasciatus, a member of the Culex pipiens complex of mosquitoes, revealed high levels of variability causing partial incompatibility between natural populations. Segregating crossing types were commonly found together within sampled sites. No correlation was evident between similarity of crossing type and environmental parameters of the sites, nor distance between sites. The nature of the observed variation did not support the hypothesis of paternally expressed nuclear 'restorer' genes. Such high levels of crossing type variation would be likely to impede attempts to control populations of the Culex pipiens complex using cytoplasmic incompatibility.     

Preliminary studies on levels of sensitivity to phosphoric esters in the Italian population of the Culex pipiens complex

A preliminary survey of the susceptibility level to organophosphate (OP) insecticides in six populations of Culex pipiens from Central Italy has been carried out. The correlation between OP resistance and highly active esterase electromorphs has also been investigated. The bioassay and electrophoresis results show an association between OP resistance and the highly active esterase allele Est-3A and confirm what it was found in populations of C. pipiens from Southern France. Est-3A appears to be a reliable indicator of OP resistance in mosquitoes of italian C. pipiens.     

Mutations in the house fly Vssc1 sodium channel gene associated with super-kdr resistance abolish the pyrethroid sensitivity of Vssc1/tipE sodium channels expressed in Xenopus oocytes

The super-kdr insecticide resistance trait of the house fly confers resistance to pyrethroids and DDT by reducing the sensitivity of the fly nervous system. The super-kdr genetic locus is tightly linked to the Vssc1 gene, which encodes a voltage-sensitive sodium channel alpha subunit that is the principal site of pyrethroid action. DNA sequence analysis of Vssc1 alleles from several independent super-kdr fly strains identified two amino acid substitutions associated with the super-kdr trait: replacement of leucine at position 1014 with phenylalanine (L1014F), which has been shown to cause the kdr resistance trait in this species, and replacement of methionine at position 918 with threonine (M918T). We examined the functional significance of these mutations by expressing house fly sodium channels containing them in Xenopus laevis oocytes and by characterizing the biophysical properties and pyrethroid sensitivities of the expressed channels using two-electrode voltage clamp. House fly sodium channels that were specifically modified by site-directed mutagenesis to contain the M918T/L1014F double mutation gave reduced levels of sodium current expression in oocytes but otherwise exhibited functional properties similar to those of wildtype channels and channels containing the L1014F substitution. However, M918T/L1014F channels were completely insensitive to high concentrations of the pyrethroids cismethrin and cypermethrin. House fly sodium channels specifically modified to contain the M918T single mutation, which is not known to exist in nature except in association with the L1014F mutation, gave very small sodium currents in oocytes. Assays of these currents in the presence of high concentrations of cismethrin suggest that this mutation alone is sufficient to abolish the pyrethroid sensitivity of house fly sodium channels. These results define the functional significance of the Vssc1 mutations associated with the super-kdr trait of the house fly and are consistent with the hypothesis that the super-kdr trait arose by selection of a second-site mutation (M918T) that confers to flies possessing it even greater resistance than the kdr allele containing the L1014F mutation.