Evidence for gene duplication in the voltage-gated sodium channel gene of Aedes aegypti

Glycine to valine substitution in the short intracellular linkers of domain II enhances I1011M-mediated sodium channel resistance to Type I pyrethroids, but not Type II pyrethroids

Pyrethroids are widely used against agricultural pests and human disease vectors due to their broad insecticidal spectrum, fast action, and low mammalian toxicity. Unfortunately, overuse of pyrethroids has led to knockdown resistance (kdr) caused by mutations in voltage-gated sodium channels. Mutation I1011M was repeatedly detected in numerous pyrethroid-resistant Aedes aegypti populations from Latin American and Brazil. In addition, mutation G923V was first reported to coexist with I1011M in permethrin/DDT-resistant Ae. aegypti, whether G923V enhances the I1011M-mediated pyrethroid resistance in sodium channels remains unclear. In this study, we introduced mutations G923V and I1011M alone or in combination into the pyrethroid-sensitive sodium channel AaNav1-1 and examined the effects of these mutations on gating properties and pyrethroid sensitivity. We found mutations I1011M and G923V + I1011M shifted the voltage dependence of activation in the depolarizing direction, and none of mutations affect the voltage-dependence of inactivation. G923V and G923V + I1011M mutations reduced the channel sensitivity to both Type I and Type II pyrethroids. However, I1011M alone conferred resistance to Type I pyrethroids, not to Type II pyrethroids. Interestingly, significant synergism effects on Type I pyrethroids were observed between mutations G923V and I1011M. The effects of all mutations on channel sensitivity to DDT were identical with those to Type I pyrethroids. Our results confirm the molecular basis of resistance mediated by mutations G923V and I1011M and may contribute to develop molecular markers for monitoring pest resistance to pyrethroids.

Genetic structure and kdr mutations in Aedes aegypti populations along a road crossing the Amazon Forest in Amapá State, Brazil

Insecticide resistance in Aedes aegypti poses a significant threat to disease control. One form of resistance, caused by kdr mutations in the NaV gene, hinders vector control efforts in Brazil. Despite genetic differences typically accumulating among isolated populations, this mosquito can actively and passively disperse through human transportation. Our study investigated the genetic structure and spread of kdr mutations in Ae. aegypti populations across six localities in Amapá State, Brazil, within the Amazonian Forest. Using 12 microsatellite loci and qPCR methods, we assessed genetic structure and identified three common kdr mutations (V410L, V1016I, and F1534C). High prevalence of kdr alleles was observed in all localities, indicating widespread distribution in Amapá State. Microsatellite analysis revealed differentiation among mosquito populations, dividing them into two distinct clusters supported by Bayesian and DAPC analyses. Oiapoque, located along the northern border with French Guiana, exhibited the highest kdr frequencies and genetic differentiation compared to other localities. Our findings suggest genetic structure in Ae. aegypti populations in Amapá State, with some passive gene flow between clusters. The study underscores the need for continuous surveillance of Ae. aegypti populations to monitor the spread of insecticide resistance and inform effective vector control strategies.

The dynamics of deltamethrin resistance evolution in Aedes albopictus has an impact on fitness and dengue virus type-2 vectorial capacity

BACKGROUND

Worldwide invasion and expansion of Aedes albopictus, an important vector of dengue, chikungunya, and Zika viruses, has become a serious concern in global public health. Chemical insecticides are the primary means currently available to control the mosquito populations. However, long-term and large-scale use of insecticides has selected for resistance in the mosquito that is accompanied by a genetic load that impacts fitness.

RESULTS

A number of laboratory strains representing different resistance mechanisms were isolated and identified from laboratory-derived, deltamethrin-resistant Ae. albopictus recovered in previous work. Resistance levels and fitness costs of the strains were evaluated and compared to characterize the evolution of the resistance genotypes and phenotypes. The heterozygous F1534S mutation (1534F/S) in the voltage gated sodium channel (vgsc) gene product (VGSC), first detected in early stages of resistance evolution, not only confers high-level resistance, but also produces no significant fitness costs, leading to the rapid spread of resistance in the population. This is followed by the increase in frequency of homozygous F1534S (1534S/S) mosquitoes that have significant fitness disadvantages, prompting the emergence of an unlinked I1532T mutation with fewer side effects and a mating advantage better adapted to the selection and reproductive pressures imposed in the experiments. Metabolic resistance with no significant fitness cost and mediating a high-tolerance resistance phenotype may play a dominant role in the subsequent evolution of resistance. The different resistant strains had similar vector competence for dengue virus type-2 (DENV-2). Furthermore, a comparative analysis of vectorial capacity revealed that increased survival due to deltamethrin resistance balanced the negative fitness cost effects and contributed to the risk of dengue virus (DENV) transmission by resistant populations. The progressive evolution of resistance results in mosquitoes with both target-site insensitivity and metabolic resistance with lower fitness costs, which further leads to resistant populations with both high resistance levels and vectorial capacity.

CONCLUSIONS

This study reveals a possible mechanism for the evolution of deltamethrin resistance in Aedes albopictus. These findings will help guide practical strategies for insecticide use, resistance management and the prevention and control of mosquito-borne disease.

Detection of sex in adults and larvae of Leptinotarsa decemlineata on principle of copy number variation

The identification of sex in larvae of insects is usually challenging or even impossible, while in adults the sexual dimorphism is usually evident. Here, we used copy number analysis to develop a method of sex detection in Colorado potato beetle (Leptinotarsa decemlineata), which has an X0 sex determination system. The X linked gene LdVssc and autosomal gene LdUBE3B were identified as appropriate target and reference loci, respectively. The copy numbers (CNV) of LdVssc in males and females were estimated using standard droplet digital PCR (ddPCR) and real-time PCR (qPCR). With both methods, CNVs were bimodally distributed (BA_ddPCR = 0.709 and BA_qPCR = 0.683) with 100% ability to distinguish females from males. The use of qPCR-based sex detection in a broad collection of 448 random CPB adults showed a perfect association (Phi = 1.0, p < 0.05) with the true sexes of adults, with mean CNV in females of 2.032 (SD = 0.227) and 0.989 in males (SD = 0.147). In the collection of 50 random 4th instar larvae, 27 females and 23 males were identified, consistent with the expected 1:1 sex ratio (p = 0.689). The method is suitable for sexing in all stages of ontogenesis. The optimal cost-effective application of the method in large populations requires the DNA extraction using CTAB, the qPCR assay in one biological replicate and three technical replicates of each marker, and the use of one randomly chosen male per run to calibrate calculation of CNV.

Perspectives on gene copy number variation and pesticide resistance

Although the generation of evolutionary diversity by gene duplication has long been known, the implications for pesticide resistance are just now beginning to be appreciated. A few examples will be cited to illustrate the point that there are many variations on the theme that gene duplication does not follow a set pattern. Transposable elements may facilitate the process but the mechanistic details are obscure and unpredictable. New developments in DNA sequencing technology and genome assembly promise to reveal more examples, yet care must be taken in interpreting the results of transcriptome and genome assemblies and independent means of validation are important. Once a specific gene family is identified, special methods generally must be used to avoid underestimating population polymorphisms and being trapped in preconceptions about the simplicity of the process. © 2021 The Author. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Molecular Analysis of Targeted Insecticide Resistance Gene Mutations in Field-Caught Mosquitos of Medical Importance From Saudi Arabia

Gene mutations on target sites can be a valuable indicator of the status of insecticide resistance. Jeddah, a global commercial and major port-of-entry city, is bearing the brunt of dengue disease burden in Saudi Arabia. In the current study, six genotypes of three codon combinations (989, 1016, and 1534) were observed on voltage-gated sodium channel (VGSC) gene in Jeddah's Aedes aegypti population, with PGF/PGC as the dominant one. Two types of introns between exon 20 and 21 on VGSC have been identified for the first time in Ae. aegypti in Saudi Arabia. Statistical and phylogenetic analyses showed that the intron type was significantly associated with the 1016 allele and may reflect the history of insecticide treatment in different continents. In addition, fixation of the L1014F allele on VGSC and G119S on acetylcholinesterase 1 gene was detected in local Culex quinquefasciatus populations, with frequencies of 95.24 and 100%, respectively. To the best of our knowledge, this is the first report of resistant-associated mutations in field-caught Cx. quinquefasciatus in Saudi Arabia. The high prevalence of insecticide resistance gene mutations in local primary mosquito vector species highlights the urgent need to carry out comprehensive insecticide resistance surveillance in Saudi Arabia.

Comparative characterization of two putative duplicated sodium channel genes in the red flour beetle, Tribolium castaneum

The insect voltage-gated sodium channel is the primary target of pyrethroids and novel efficient insecticides such as indoxacarb and metaflumizone. In this study, we cloned and characterized two putative sodium channel genes, TcNav1 and TcNav2, in Tribolium castaneum. The composite TcNav1 and TcNav2 encode a protein of 2045 and 2037 amino acid residues, sharing 76.1% and 75.5% amino acid identity with Drosophila para, respectively. Comparative analysis of genomic organization showed that TcNav1 and TcNav2 contain 26 and 27 exons, respectively. Analysis of the expression patterns showed that the mRNA levels of TcNav1 and TcNav2 were predominantly expressed in the head. RNAi-mediated knockdown of both TcNav1 and TcNav2 adversely affected adult emergence and significantly decreased sensitivity to deltamethrin. Significantly reduced pupation rate and sensitivity to beta-cypermethrin were observed after injection of siRNA targeting TcNav1 but not TcNav2. Taken together, we provide evidence that sodium channel gene has undergone duplication in T. castaneum, resulting in diversified developmental and toxicological functions.

The Midgut Microbiota of Colombian Aedes aegypti Populations with Different Levels of Resistance to the Insecticide Lambda-cyhalothrin

Insecticide resistance in Aedes aegypti populations is a problem that hinders vector control and dengue prevention programs. In this study, we determined the susceptibility of Ae. aegypti populations from six Colombian regions to the pyrethroid lambda-cyhalothrin and evaluated the presence of the V1016I mutation in the sodium channel gene, which has been broadly involved in the resistance to this insecticide. The diversity of the gut microbiota of these mosquito populations was also analyzed. Only mosquitoes from Bello were susceptible to lambda-cyhalothrin and presented a lower allelic frequency of the V1016I mutation. Remarkably, there was not an important change in allelic frequencies among populations with different resistance ratios, indicating that other factors or mechanisms contributed to the resistant phenotype. Treatment of mosquitoes with antibiotics led us to hypothesize that the intestinal microbiota could contribute to the resistance to lambda-cyhalothrin. Beta diversity analysis showed significant differences in the species of bacteria present between susceptible and resistant populations. We identified 14 OTUs of bacteria that were unique in resistant mosquitoes. We propose that kdr mutations are important in the development of resistance to lambda-cyhalothrin at low insecticide concentrations but insect symbionts could play an essential role in the metabolization of pyrethroid insecticides at higher concentrations, contributing to the resistant phenotype in Ae. aegypti.

Kdr genotyping in Aedes aegypti from Brazil on a nation-wide scale from 2017 to 2018

Insecticide resistance is currently a threat to the control of Aedes agypti, the main vector of arboviruses in urban centers. Mutations in the voltage gated sodium channel (NaV), known as kdr (knockdown resistance), constitute an important selection mechanism for resistance against pyrethroids. In the present study, we investigated the kdr distribution for the Val1016Ile and Phe1534Cys alterations in Ae. aegypti from 123 Brazilian municipalities, based on SNP genotyping assays in over 5,500 mosquitoes. The alleles NaVS (1016Val+ + 1534Phe+), NaVR1 (1016Val+ + 1534Cyskdr) and NaVR2 (1016Ilekdr + 1534Cyskdr) were consistently observed, whereas kdr alleles have rapidly spread and increased in frequency. NaVS was the less frequent allele, mostly found in Northeastern populations. The highest allelic frequencies were observed for NaVR1, especially in the North, which was fixed in one Amazonian population. The double kdr NaVR2 was more prevalent in the Central-west and South-eastern populations. We introduce the 'kdr index', which revealed significant spatial patterns highlighting two to three distinct Brazilian regions. The 410L kdr mutation was additionally evaluated in 25 localities, evidencing that it generally occurs in the NaVR2 allele. This nationwide screening of a genetic mechanism for insecticide resistance is an important indication on how pyrethroid resistance in Ae. aegypti is evolving in Brazil.

A new knockdown resistance (kdr) mutation, F1534L, in the voltage-gated sodium channel of Aedes aegypti, co-occurring with F1534C, S989P and V1016G

BACKGROUND

Aedes aegypti is a primary vector of dengue, chikungunya and Zika infections in India. In the absence of specific drugs or safe and effective vaccines for these infections, their control relies mainly on vector control measures. The emergence of insecticide resistance in vectors, especially against pyrethroids, is a serious threat to the insecticide-based vector control programme. This study reports the presence of multiple knockdown resistance (kdr) mutations present in an Ae. aegypti population from Bengaluru (India), including a new mutation F1534L.

METHODS

Aedes aegypti collected from Bengaluru were subjected to insecticide susceptibility tests with DDT, deltamethrin and permethrin. The DNA sequencing of partial domain II, III and IV of the voltage-gated sodium channel (VGSC) was performed to screen kdr mutations present in the population and PCR-based assays were developed for their detection. Genotyping of kdr mutations was done using PCR-based assays, allelic frequencies were determined, and tests of genetic association of kdr mutations with the insecticide resistance phenotype were performed.

RESULTS

The Ae. aegypti population was resistant to DDT, deltamethrin and permethrin. The DNA sequencing of the VGSC revealed the presence of four kdr mutations, i.e. S989P and V1016G in domain II and two alternative kdr mutations F1534C and F1534L in domain III. Allele-specific PCR assays (ASPCR) were developed for the detection of kdr mutations S989P and V1016G and an existing PCR-RFLP based strategy was modified for the genotyping of all three known kdr mutations in domain III (F1534L, F1534C and T1520I). Genotyping of Ae. aegypti samples revealed a moderate frequency of S989P/V1016G (18.27%) and F1534L (17.48%), a relatively high frequency of F1534C (50.61%) and absence of T1520I in the population. Mutations S989P and V1016G were in complete linkage disequilibrium in this population while they were in linkage equilibrium with kdr mutations F1534C and F1534L. The alleles F1534C and F1534L are genetically associated with permethrin resistance.

CONCLUSIONS

A new kdr mutation, F1534L, was found in an Ae. aegypti population from Bengaluru (India), co-occurring with the other three mutations S989P, V1016G and F1534C. The findings of a new mutation have implications for insecticide resistance management.

Chronology of sodium channel mutations associated with pyrethroid resistance in Aedes aegypti

Aedes aegypti is the primary mosquito vector of dengue, yellow fever, Zika and chikungunya. Current strategies to control Ae. aegypti rely heavily on insecticide interventions. Pyrethroids are a major class of insecticides used for mosquito control because of their fast acting, highly insecticidal activities and low mammalian toxicity. However, Ae. aegypti populations around the world have begun to develop resistance to pyrethroids. So far, more than a dozen mutations in the sodium channel gene have been reported to be associated with pyrethroid resistance in Ae. aegypti. Co-occurrence of resistance-associated mutations is common in pyrethroid-resistant Ae. aegypti populations. As global use of pyrethroids in mosquito control continues, new pyrethroid-resistant mutations keep emerging. In this microreview, we compile pyrethroid resistance-associated mutations in Ae. aegypti in a chronological order, as they were reported, and summarize findings from functional evaluation of these mutations in an in vitro sodium channel expression system. We hope that the information will be useful for tracing possible evolution of pyrethroid resistance in this important human disease vector, in addition to the development of methods for global monitoring and management of pyrethroid resistance in Ae. aegypti.

Evaluation of resistance to pyrethroid and organophosphate adulticides and kdr genotyping in Aedes aegypti populations from Roraima, the northernmost Brazilian State

Background

Roraima, the northernmost State in Brazil, borders Venezuela and Guyana. Although mostly covered by the tropical forests, the urban centers of this state are highly infested with Ae. aegypti and are endemic for dengue, Zika and chikungunya. We accessed the insecticide resistance status of Ae. aegypti populations from the capital Boa Vista, two cities on international borders (Pacaraima and Bonfim) and Rorainópolis bordering Amazonas State, in order to evaluate the chemical control efficacy in these localities.

Methods

Tests with World Health Organization (WHO)-like tubes impregnated with the pyrethroid deltamethrin (0.05% and 0.12%) and the organophosphate malathion (0.7%) were conducted with Ae. aegypti from Boa Vista, Pacaraima, Bonfim and Rorainópolis, collected in 2016 and 2018. Genotyping of kdr mutations, related to resistance to pyrethroids, was performed for the SNP variations at sites 1016 and 1534 of the voltage gated sodium channel gene (Na_V) with a TaqMan qPCR approach.

Results

Aedes albopictus was absent in our collections, and therefore only Ae. aegypti was tested. All Ae. aegypti populations were susceptible to 0.7% malathion in 2016; however, mortality dropped to under 90% in Bonfim and Pacaraima populations in 2018. All populations were resistant to 0.05% deltamethrin in both years. The time that 50% of females suffered knockdown (KdT₅₀) under exposure to 0.05% deltamethrin was 3.3–5.9-fold longer in mosquitoes from the natural populations compared to the susceptible Rockefeller strain. Only the Pacaraima population (2018) remained resistant to 0.12% deltamethrin. Kdr genotyping revealed the absence of the wild-type Na_VS haplotype (1016Val + 1534Phe) in the populations from Roraima, indicating that all tested insects had a genetic background for pyrethroid resistance. The double kdr Na_VR2 haplotype (1016Ile + 15434Cys) was present in higher frequencies in all populations except for Rorainópolis, where this haplotype seems to have arrived recently.

Conclusions

These results are important for the knowledge about insecticide resistance status of Ae. aegypti populations from Roraima and will help improve vector control strategies that may be applied to diverse localities under similar geographical and urban conditions.

Evolution of kdr haplotypes in worldwide populations of Aedes aegypti: Independent origins of the F1534C kdr mutation

Aedes aegypti is the primary vector of dengue, chikungunya, Zika, and urban yellow fever. Insecticides are often the most effective tools to rapidly decrease the density of vector populations, especially during arbovirus disease outbreaks. However, the intense use of insecticides, particularly pyrethroids, has selected for resistant mosquito populations worldwide. Mutations in the voltage gated sodium channel (Na_V) are among the principal mechanisms of resistance to pyrethroids and DDT, also known as “knockdown resistance,” kdr. Here we report studies on the origin and dispersion of kdr haplotypes in samples of Ae. aegypti from its worldwide distribution. We amplified the IIS6 and IIIS6 Na_V segments from pools of Ae. aegypti populations from 15 countries, in South and North America, Africa, Asia, Pacific, and Australia. The amplicons were barcoded and sequenced using NGS Ion Torrent. Output data were filtered and analyzed using the bioinformatic pipeline Seekdeep to determine frequencies of the IIS6 and IIIS6 haplotypes per population. Phylogenetic relationships among the haplotypes were used to infer whether the kdr mutations have a single or multiple origin. We found 26 and 18 haplotypes, respectively for the IIS6 and IIIS6 segments, among which were the known kdr mutations 989P, 1011M, 1016I and 1016G (IIS6), 1520I, and 1534C (IIIS6). The highest diversity of haplotypes was found in African samples. Kdr mutations 1011M and 1016I were found only in American and African populations, 989P + 1016G and 1520I + 1534C in Asia, while 1534C was present in samples from all continents, except Australia. Based primarily on the intron sequence, IIS6 haplotypes were subdivided into two well-defined clades (A and B). Subsequent phasing of the IIS6 + IIIS6 haplotypes indicates two distinct origins for the 1534C kdr mutation. These results provide evidence of kdr mutations arising de novo at specific locations within the Ae. aegypti geographic distribution. In addition, our results suggest that the 1534C kdr mutation had at least two independent origins. We can thus conclude that insecticide selection pressure with DDT and more recently with pyrethroids is selecting for independent convergent mutations in Na_V.

Insecticide resistance is a global threat for the control of Aedes aegypti, the mosquito vector of aboviruses such as dengue, chikungunya and Zika. Mutations in the voltage gated sodium channel (Na_V), known as kdr, are one of the principal mechanisms related to resistance to pyrethroids, the class of insecticide most employed worldwide inside and around residences. We investigate whether the same kdr mutations found in Ae. aegypti populations from distinct regions of the world have a common origin and subsequently dispersed or if they emerged in unrelated populations at distinct moments. By evaluating the sequences of two fragments of the Na_V gene, obtained from DNA collections of Ae. aegypti from several countries, we found at least two independent origins for the F1534C kdr mutation in American, African and Asian populations. There was no evidence for multiple origins of the common kdr mutations V1016I and P989S + V1016G, which were exclusive to American and Asian populations. Our results increase our knowledge of insecticide resistance evolution in one of the main arboviral mosquito vectors of major global diseases.

High-throughput genotyping of a full voltage-gated sodium channel gene via genomic DNA using target capture sequencing and analytical pipeline MoNaS to discover novel insecticide resistance mutations

In insects, the voltage-gated sodium channel (VGSC) is the primary target site of pyrethroid insecticides. Various amino acid substitutions in the VGSC protein, which are selected under insecticide pressure, are known to confer insecticide resistance. In the genome, the VGSC gene consists of more than 30 exons sparsely distributed across a large genomic region, which often exceeds 100 kbp. Due to this complex genomic structure, it is often challenging to genotype full coding nucleotide sequences (CDSs) of VGSC from individual genomic DNA (gDNA). In this study, we designed biotinylated oligonucleotide probes from CDSs of VGSC of Asian tiger mosquito, Aedes albopictus. The probe set effectively concentrated (>80,000-fold) all targeted regions of gene VGSC from pooled barcoded Illumina libraries each constructed from individual A. albopictus gDNAs. The probe set also captured all orthologous VGSC CDSs, except some tiny exons, from the gDNA of other Culicinae mosquitos, A. aegypti and Culex pipiens complex, with comparable efficiency as a result of the high nucleotide-level conservation of VGSC. To improve efficiency of the downstream bioinformatic process, we developed an automated pipeline-MoNaS (Mosquito Na+ channel mutation Search)-which calls amino acid substitutions in the VGSC from NGS reads and compares those to known resistance mutations. The proposed method and our bioinformatic tool should facilitate the discovery of novel amino acid variants conferring insecticide resistance on VGSC and population genetic studies on resistance alleles (with respect to the origin, selection, and migration etc.) in both clinically and agriculturally important insect pests.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

Point mutations in the voltage-gated sodium channel gene associated with pyrethroid resistance in Iranian populations of the European red mite Panonychus ulmi

The European red mite Panonychus ulmi (Koch) is a major pest of apple trees worldwide and causes significant damage to apple orchards in Iran. Pyrethroid insecticides/acaricides, such as fenpropathrin and fenvalerate, are widely used to control P. ulmi, but their long-term use may lead to low efficacy. Earlier studies investigating pyrethroid resistance in closely related mites such as Tetranychus urticae revealed that pyrethroid resistance was associated with point mutations in the voltage-gated sodium channel gene (vgsc). The aim of this study was to investigate the biochemical and molecular mechanisms of fenpropathrin and fenvalerate resistance in Iranian populations of P. ulmi. Pyrethroid toxicity bioassays were carried out on different P. ulmi field populations. Marand (resistance ratio, RR = 149), Maraqeh (RR = 90) and Mianeh2 (RR = 71) populations exhibited high levels of resistance to fenpropathrin, compared to a susceptible field population (Shahin Dej). Resistance was also observed for fenvalerate with resistance ratio's ranging from 2- to 20-fold. Synergism experiments and enzyme activity assays predicted a minor role for classical detoxification enzymes. In contrast, two amino acid substitutions in the VGSC, L1024V and F1538I, that were previously shown to confer pyrethroid resistance, were detected in all three resistant P. ulmi populations and point towards target-site insensitivity as the most likely resistance mechanism. Furthermore, sequencing after cloning of vgsc fragments from single haploid males revealed the presence of multiple copies of vgsc in a highly resistant strain. The link between resistance mutations and vgsc copy number variation should be the subject of future study, as this might be used to develop molecular markers for monitoring pyrethroid resistance of P. ulmi in the field.

Voltage-gated sodium channel intron polymorphism and four mutations comprise six haplotypes in an Aedes aegypti population in Taiwan

Background

Knockdown resistance (kdr) to dichlorodiphenyltrichloroethane (DDT) and pyrethroids is known to link amino acid substitutions in the voltage-gated sodium channel (VGSC) in Aedes aegypti. Dengue fever primarily transmitted by Ae. aegypti is an annual public health issue in Taiwan. Accordingly, pyrethroid insecticides have been heavily used for decades to control mosquito populations in the summer and autumn. In Taiwan, an Ae. aegypti population with two VGSC mutations, V1016G and D1763Y, was described previously.

Methodology/Principal finding

Aedes aegypti (G0) were collected in Tainan and Kaohsiung in southern Taiwan. The VGSC gene polymorphisms of the kdr mutations and the intron flanked by exons 20 and 21 were verified. The first generation offspring (G1) were used to measure the resistance level to cypermethrin, a pyrethroid insecticide currently used in Taiwan. In addition to V1016G and D1763Y, we describe two new mutations, S989P and F1534C, which have not been reported in Taiwan. Moreover, we also identify two types (groups A and B) of introns between exons 20 and 21. Intriguingly, the kdr mutations S989P, V1016G and D1763Y are strictly located on the haplotype harboring the group A intron, whereas F1534C links to the group B intron. When those data were taken together, we proposed the following six haplotypes for VGSC genes in Taiwan today: (i)S989-intron A-V1016-F1534-D1763, (ii)S989-intron A-V1016G-F1534-D1763, (iii)S989P-intron A-V1016G-F1534-D1763, (iv)S989-intron A-V1016G-F1534-D1763Y, (v)S989-intron B-V1016-F1534-D1763 and (vi)S989-intron B-V1016-F1534C-D1763. Triple heterozygous mutations of either S989P/V1016G/F1534C or V1016G/F1534C/D1763Y can be found in one single Ae. aegypti mosquito. The proportions of the VGSC mutations were relevant to cypermethrin resistance. Notably, the presence of S989P and V1016G in the population could be a helpful reference to predict the resistance level to cypermethrin. This is the first study to demonstrate the coexistence of four kdr mutations in a population of Ae. aegypti.

Conclusions/Significance

Four kdr mutations (S989P, V1016G, F1534C and D1763Y) and two intron forms (Group A and B) were commonly found in local Ae. aegypti populations in Taiwan.

VGSC mutations of Aedes aegypti threaten vector control programs and have been brought to attention in dengue endemic areas. Taiwan has suffered dengue outbreaks, which usually begin from an imported case. Pyrethroid insecticides were used to kill infectious females and adults in the surrounding environment of each suspected case. In Taiwan, V1016G and F1763Y mutations of VGSC have been described previously. Here, we further describe two additional amino acid substitutions (S989P, F1534C) and two forms of the intron between exon 20 and 21. We also propose six haplotypes for VGSC genes in Taiwan today. In conclusion, four kdr mutations (S989P, V1016G, F1534C and D1763Y) and two intron forms (Group A and B) are commonly found in local Ae. aegypti populations in Taiwan.

Insecticide Resistance and Fitness: The Case of Four Aedes aegypti Populations from Different Brazilian Regions

Background

Chemical control is still a major strategy to constrain vector density and mitigate pathogen transmission. However, insecticide overuse poses a high selective pressure, favouring the spread of resistance alleles in natural populations. In an insecticide-free environment, a fitness cost is expected in resistant insects when compared to susceptible counterparts. This study investigates whether insecticide resistance to an organophosphate (temephos) and a pyrethroid (deltamethrin) is associated with fitness traits in four Aedes aegypti wild populations sampled every three months over one year.

Findings

We measured development time from larvae to adult, female survival, wing length, fecundity, and adult resistance to starvation in field insecticide resistant Ae. aegypti populations four times over a year. These results were confronted with resistance levels to temephos and deltamethrin and with potentially related mechanisms, including a kdr mutation in the pyrethroid target site. No differences in fitness cost were found after contrasting mosquitoes from the same population collected throughout a year, irrespective of differences in insecticide resistance levels. Additionally, significant differences were not observed among field populations. However, compared to the reference strain Rockefeller, field females survived significantly less. Moreover, larval development was equal or slower in three out of four field populations. In no case differences were evidenced in starvation tolerance, wing length, and fecundity.

Conclusions

Overall, field resistant mosquitoes seemed to have a slight fitness disadvantage when compared with the Rockefeller susceptible strain which might represent a potential fitness cost of insecticide resistance. However, after comparing Ae. aegypti from the same population but sampled at different moments, or from different field populations, mosquito life-history traits varied independently of resistance ratios. The metabolic deviations necessary to overcome the adverse effects of insecticides may cause an energy trade-off that affects energy allocation and ultimately basic demands of insect biology. The extent of fitness cost due to insecticide resistance is critical information to delay the evolution of resistance in wild vector populations.

Levels of Resistance to Pyrethroid among Distinct kdr Alleles in Aedes aegypti Laboratory Lines and Frequency of kdr Alleles in 27 Natural Populations from Rio de Janeiro, Brazil

Background

Several mutations in voltage gated sodium channel (Na_V) have been identified in Aedes aegypti populations worldwide. However, only few are related to knockdown resistance to pyrethroids, most of which with variations in the 1016 and 1534 Na_V sites. In Brazil, at least two Na_V alleles are known: Na_VR1, with a substitution in the 1534 (1016 Val⁺ + 1534 Ile^kdr) and Na_VR2, with substitutions in both 1016 and sites (1016Ile^kdr + 1534Cys^kdr). There is also the duplication in the Na_V gene, with one copy carrying the substitution Ile1011Met, although its effects on pyrethroid resistance remain to be clarified. Our goals in this study were (1) to determine the role of each kdr Na_V allele and the duplication on pyrethroid resistance and (2) to screen the frequency of the kdr alleles in 27 several natural Ae. aegypti populations from the metropolitan region of Rio de Janeiro.

Methods

Pyrethroid resistance was evaluated by a knockdown time (KdT) assay, an adaptation of the WHO test tubes with paper impregnated with deltamethrin. We used laboratory-selected Ae. aegypti lineages: R1R1 and R2R2 (homozygous for the kdr Na_VR1 and Na_VR2 alleles, respectively), Dup (with duplication in the Na_V gene), Rockefeller (the susceptibility reference control), and F1 hybrids among them. Genotyping of both 1016 and 1534 Na_V sites was performed in 811 Ae. aegypti sampled from 27 localities from Rio de Janeiro (17), Niterói (6) and Nova Iguaçu (4) cities, Rio de Janeiro State, Brazil, with a TaqMan real time PCR approach.

Results

The laboratory lineages R1R1, R2R2, and R1R2 were the only ones that needed more than 60 minutes to knock down all the insects exposed to the pyrethroid, being the KdT R2R2 > R1R2 > R1R1, corroborating the recessive nature of the kdr mutations. Frequency of kdr alleles Na_VR1 and Na_VR2 in field-caught mosquitoes varied from 0 to 52% and 43 to 86%, respectively, evidencing high levels of “resistant genotypes” (R1R1, R1R2, and R2R2), which together summed 60 to 100% in Ae. aegypti populations from Rio de Janeiro.

Conclusions

The Na_VR1 and Na_VR2 kdr alleles confer resistance to the pyrethroid deltamethrin in homozygotes and R1R2 heterozygotes, being the R2R2 most resistant genotype. The allele containing duplication in the Na_V gene, with a mutation in the 1011 site, did not confer resistance under the tested conditions. The frequencies of the “resistant genotypes” are elevated in Ae. aegypti natural populations from Rio de Janeiro.

Effect of Relaxation of Deltamethrin Pressure on Metabolic Resistance in a Pyrethroid-Resistant Aedes aegypti (Diptera: Culicidae) Strain Harboring Fixed P989P and G1016G kdr Alleles

Mutation of the voltage-gated sodium channel genes or knockdown resistance (kdr) and metabolic resistance in Aedes aegypti (L.) (Diptera: Culicidae) are important resistance mechanisms against pyrethroids. The present study investigated the effect of relaxation of deltamethrin selection pressure on the level of mixed-function oxidases (MFO), when the allele frequency of S989P+V1016G mutations is fixed in a resistant Ae. aegypti strain (UPK-R) from Chiang Mai, Thailand. The mosquitoes were divided into two groups, exposure and nonexposure groups, and maintained for 12 generations in an insectary room. Adults of the exposure group (F3 to F12) were treated with 0.05% deltamethrin-impregnated papers. The median lethal concentrations (LC50) of deltamethrin of larvae were determined by World Health Organization (WHO) bioassay. MFO activity was determined in F0 and F12. The results revealed that there was a decreasing trend of adult mortality rates in the exposure group over time. The larval LC50 values of the exposure group were gradually increased, whereas those of the nonexposure group were gradually decreased. The level of MFO activity in the nonexposure group (F12) was lower than the parent and exposure groups (F12) by 1.5 and 4-fold in the larvae, respectively, and 1.5 and 2.5-fold in the adult females, respectively. However, the frequency of P989+G1016 alleles in both groups was 100% up to F12 when the experiment ended. This study indicates that there was a significant but small reduction in the activity levels of MFOs when pyrethroid selection pressure is relaxed in this kdr strain of Ae. aegypti.

Copy number variation (CNV) and insecticide resistance in mosquitoes: evolving knowledge or an evolving problem?

Copy number variation (CNV) in insect genomes is a rich source of potentially adaptive polymorphism which may help overcome the constraints of purifying selection on conserved genes and/or permit elevated transcription. Classic studies of amplified esterases and acetylcholinesterase duplication in Culex pipiens quantified evolutionary dynamics of CNV driven by insecticidal selection. A more complex and potentially medically impactful form of CNV is found in Anopheles gambiae, with both heterogeneous duplications and homogeneous amplifications strongly linked with insecticide resistance. Metabolic gene amplification, revealed by shotgun sequencing, appears common in Aedes aegypti, but poorly understood in other mosquito species. Many methodologies have been used to detect CNV in mosquitoes, but relatively few can detect both duplications and amplifications, and contrasting methods should be combined. Genome scans for CNV have been rare to date in mosquitoes, but offer immense potential to determine the overall role of CNV as a component of resistance mechanisms.

The impact of insecticide applications on the dynamics of resistance: The case of four Aedes aegypti populations from different Brazilian regions

Background

In the tropics, the utilization of insecticides is still an important strategy for controlling Aedes aegypti, the principle vector of dengue, chikungunya and Zika viruses. However, increasing insecticide resistance in Ae. aegypti populations might hinder insecticide efficacy on a long-term basis. It will be important to understand the dynamics and evolution of insecticide resistance by assessing its frequency and the mechanisms by which it occurs.

Methodology/Principal findings

The insecticide resistance status of four Brazilian Ae. aegypti populations was monitored. Quantitative bioassays with the major insecticides employed in the country was performed: the adulticide deltamethrin (a pyrethroid—PY) and the larvicides, temephos (an organophosphate) and diflubenzuron (a chitin synthesis inhibitor). Temephos resistance was detected in all populations although exhibiting a slight decrease over time probably due to the interruption of field use. All vector populations were susceptible to diflubenzuron, recently introduced in the country to control Ae. aegypti. Resistance against deltamethrin was extremely high in three populations. Molecular assays investigated substitutions in the voltage gated sodium channel (Na_V), the PY target site, at positions 1011, 1016 and 1534. Elevated frequencies of substitutions Val1016Ile and Phe1534Cys related to high PY resistance levels were identified. Biochemical assays detected alterations in the activities of two detoxifying enzyme classes related to metabolic resistance, glutathion-S-transferases and esterases. The results obtained were evaluated in the context of both recent insecticide use and the records of dengue incidence in each locality.

Conclusions/Significance

The four Ae. aegypti populations evaluated were resistant to the neurotoxic insecticides, temephos and deltamethrin. However, they were still susceptible to diflubenzuron. A probable correlation between adult insect resistance to PY and the domestic application of insecticides is discussed, pointing to the need for awareness measures regarding the correct utilization by citizens. This work aims to contribute to the efficient and rational management of Ae. aegypti control of both larvae and adults.

Among the pathogens transmitted by Aedes aegypti, dengue virus is the most important due to the number of people affected or at risk and the high rate of mortality worldwide. The confirmation that Ae. aegypti is also the vector of Zika, chikungunya and urban yellow fever poses serious consequences for public health, pointing to the need of reevaluating current vector control strategies. Although there is growing recognition of the importance of social participation and community engagement to prevent high levels of infestation, insecticides are considered important vector control tools. Nevertheless, the massive and indiscriminate adoption of insecticides to control larvae and adults contributes to resistance spread. In particular, the domestic use of adulticides, especially in epidemic seasons, is assumed to induce high levels of resistance in Ae. aegypti populations. However, the consequences of insecticide interruption upon the resistance of field populations has been less investigated. We evaluated, in four Brazilian regions over one year, the dynamics of dengue vector population resistance to the principal insecticides used in the country. The main resistance mechanisms were also investigated. Data are discussed taking into account the potential relationship among dengue outbreaks, public and private chemical control and insecticide resistance.

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations

BACKGROUND

Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species.

METHODOLOGY/PRINCIPAL FINDINGS

Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci).

CONCLUSIONS/SIGNIFICANCE

Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex.

Detection and quantitation of copy number variation in the voltage-gated sodium channel gene of the mosquito Culex quinquefasciatus

Insecticide resistance is typically associated with alterations to the insecticidal target-site or with gene expression variation at loci involved in insecticide detoxification. In some species copy number variation (CNV) of target site loci (e.g. the Ace-1 target site of carbamate insecticides) or detoxification genes has been implicated in the resistance phenotype. We show that field-collected Ugandan Culex quinquefasciatus display CNV for the voltage-gated sodium channel gene (Vgsc), target-site of pyrethroid and organochlorine insecticides. In order to develop field-applicable diagnostics for Vgsc CN, and as a prelude to investigating the possible association of CN with insecticide resistance, three assays were compared for their accuracy in CN estimation in this species. The gold standard method is droplet digital PCR (ddPCR), however, the hardware is prohibitively expensive for widespread utility. Here, ddPCR was compared to quantitative PCR (qPCR) and pyrosequencing. Across all platforms, CNV was detected in ≈10% of mosquitoes, corresponding to three or four copies (per diploid genome). ddPCR and qPCR-Std-curve yielded similar predictions for Vgsc CN, indicating that the qPCR protocol developed here can be applied as a diagnostic assay, facilitating monitoring of Vgsc CN in wild populations and the elucidation of association between the Vgsc CN and insecticide resistance.

Spatial variation of insecticide resistance in the dengue vector Aedes aegypti presents unique vector control challenges

BACKGROUND

Dengue is a major public health problem in Mexico, where the use of chemical insecticides to control the principal dengue vector, Aedes aegypti, is widespread. Resistance to insecticides has been reported in multiple sites, and the frequency of kdr mutations associated with pyrethroid resistance has increased rapidly in recent years. In the present study, we characterized patterns of insecticide resistance in Ae. aegypti populations in five small towns surrounding the city of Merida, Mexico.

METHODS

A cross-sectional, entomological survey was performed between June and August 2013 in 250 houses in each of the five towns. Indoor resting adult mosquitoes were collected in all houses and four ovitraps were placed in each study block. CDC bottle bioassays were conducted using F0-F2 individuals reared from the ovitraps and kdr allele (Ile1016 and Cys1534) frequencies were determined.

RESULTS

High, but varying, levels of resistance to chorpyrifos-ethyl was detected in all study towns, complete susceptibility to bendiocarb in all except one town, and variations in resistance to deltamethrin between towns, ranging from 63-88% mortality. Significant associations were detected between deltamethrin resistance and the presence of both kdr alleles. Phenotypic resistance was highly predictive of the presence of both alleles, however, not all mosquitoes containing a mutant allele were phenotypically resistant. An analysis of genotypic differentiation (exact G test) between the five towns based on the adult female Ae. aegypti collected from inside houses showed highly significant differences (p < 0.0001) between genotypes for both loci. When this was further analyzed to look for fine scale differences at the block level within towns, genotypic differentiation was significant for both loci in San Lorenzo (Ile1016, p = 0.018 and Cys1534, p = 0.007) and for Ile1016 in Acanceh (p = 0.013) and Conkal (p = 0.031).

CONCLUSIONS

The results from this study suggest that 3 years after switching chemical groups, deltamethrin resistance and a high frequency of kdr alleles persisted in Ae. aegypti populations. The spatial variation that was detected in both resistance phenotypes and genotypes has practical implications, both for vector control operations as well as insecticide resistance management strategies.

Involvement of metabolic resistance and F1534C kdr mutation in the pyrethroid resistance mechanisms of Aedes aegypti in India

Pesticide resistance poses a serious problem for worldwide mosquito control programs. Resistance to insecticides can be caused by an increased metabolic detoxification of the insecticide and/or by target site insensitivity. In the present study, we estimated the tolerance of Indian Aedes aegypti populations using adult bioassays that revealed high resistance levels of the field populations to permethrin (RR-6, 5.8 and 5.1 folds) compared to our susceptible population. Enzymatic assays revealed increased activities of glutathione S-transferase and carboxylesterase enzymes in the field populations comparatively to the susceptible population. PBO synergist assays did not confirm that cytochrome P450 monooxygenase metabolic detoxification acted as a major cause of resistance. Hence the role of target site resistance was therefore investigated. A single substitution Phe1534Cys in the voltage gated sodium channel was found in domain III, segment 6 (III-S6) of the resistance populations (allele frequency=0.59, 0.51 and 0.47) suggesting its potential role in permethrin resistance in A. aegypti.

Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing

The capacity of mosquitoes to resist insecticides threatens the control of diseases such as dengue and malaria. Until alternative control tools are implemented, characterizing resistance mechanisms is crucial for managing resistance in natural populations. Insecticide biodegradation by detoxification enzymes is a common resistance mechanism; however, the genomic changes underlying this mechanism have rarely been identified, precluding individual resistance genotyping. In particular, the role of copy number variations (CNVs) and polymorphisms of detoxification enzymes have never been investigated at the genome level, although they can represent robust markers of metabolic resistance. In this context, we combined target enrichment with high-throughput sequencing for conducting the first comprehensive screening of gene amplifications and polymorphisms associated with insecticide resistance in mosquitoes. More than 760 candidate genes were captured and deep sequenced in several populations of the dengue mosquito Ae. aegypti displaying distinct genetic backgrounds and contrasted resistance levels to the insecticide deltamethrin. CNV analysis identified 41 gene amplifications associated with resistance, most affecting cytochrome P450s overtranscribed in resistant populations. Polymorphism analysis detected more than 30,000 variants and strong selection footprints in specific genomic regions. Combining Bayesian and allele frequency filtering approaches identified 55 nonsynonymous variants strongly associated with resistance. Both CNVs and polymorphisms were conserved within regions but differed across continents, confirming that genomic changes underlying metabolic resistance to insecticides are not universal. By identifying novel DNA markers of insecticide resistance, this study opens the way for tracking down metabolic changes developed by mosquitoes to resist insecticides within and among populations.

A sodium channel variant in Aedes aegypti as a candidate pathogen sensor for viral-associated molecular patterns

Recent work demonstrated that a splice variant of a human macrophage voltage-gated sodium channel expressed on endosomes acts as an intracellular sensor for dsRNA, a viral-associated molecular pattern. Here our goal was to identify a candidate gene in a clinically relevant invertebrate model with related cellular and pattern recognition properties. The para gene in drosophila and other insects encodes voltage-gated sodium channels with similar electrophysiological properties to those found in vertebrate excitable membranes. A database search revealed that the AAEL006019 gene in Aedes aegypti, the yellow fever mosquito, encodes a voltage-gated sodium channel that is distinct from genes that encode para-like sodium channels. As compared to para-like channels, the protein products from this gene have deletions in the N-terminus and in the DII-DIII linker region. When over-expressed in an Aedes aegypti cell line, CCL-125, the AAEL006019 channel demonstrated cytoplasmic expression on vesicular-like organelles. Electrophysiologic analysis revealed that the channel mediates small inward currents that are enhanced by synthetic mimics of viral-derived ssRNA, R848 and ORN02, but not the dsRNA mimic, poly I:C. R848 treatment of CCL-125 cells that express high levels of the channels led to increased expression of RelA and Ago2, two mediators of insect innate immunity. These results suggest that the AAEL006019 channel acts as an intracellular pathogen sensor for ssRNA molecular patterns.

Genotype to phenotype, the molecular and physiological dimensions of resistance in arthropods

The recent accumulation of molecular studies on mutations in insects, ticks and mites conferring resistance to insecticides, acaricides and biopesticides is reviewed. Resistance is traditionally classified by physiological and biochemical criteria, such as target-site insensitivity and metabolic resistance. However, mutations are discrete molecular changes that differ in their intrinsic frequency, effects on gene dosage and fitness consequences. These attributes in turn impact the population genetics of resistance and resistance management strategies, thus calling for a molecular genetic classification. Mutations in structural genes remain the most abundantly described, mostly in genes coding for target proteins. These provide the most compelling examples of parallel mutations in response to selection. Mutations causing upregulation and downregulation of genes, both in cis (in the gene itself) and in trans (in regulatory processes) remain difficult to characterize precisely. Gene duplications and gene disruption are increasingly reported. Gene disruption appears prevalent in the case of multiple, hetero-oligomeric or redundant targets.

Contrasting patterns of insecticide resistance and knockdown resistance (kdr) in the dengue vectors Aedes aegypti and Aedes albopictus from Malaysia

Background

Knowledge on the extent, distribution and mechanisms of insecticide resistance is essential for successful insecticide-based dengue control interventions. Here, we report an extensive resistance profiling of the dengue vectors Aedes aegypti and Aedes albopictus across Malaysia and establish the contribution of knockdown resistance mechanism revealing significant contrast between both species.

Methods

Aedes mosquitoes were collected from four states in Malaysia in 2010 using ovitraps and tested against six major insecticides using WHO bioassays. Knockdown resistance (kdr) was investigated in both species.

Results

A moderate resistance to temephos was detected from samples collected in 2010 in Penang, Kuala Lumpur, Johor Bharu and Kota Bharu (1.5 < RR < 3.3). A widespread and multiple resistances was observed in Ae. aegypti particularly against pyrethroids, DDT and bendiocarb. Mosquitoes from Kuala Lumpur consistently had the highest resistance levels and was the only population showing a moderate resistance to malathion (91% mortality). The resistance profile of Ae. albopictus contrasted to Ae. aegypti with full susceptibility to pyrethroids except in Kuala Lumpur where moderate resistance is observed. PBO synergist assays suggest metabolic resistance mechanisms play a major role in resistance in both species. Two kdr mutations, F1534C and V1016G, were detected in Ae. aegypti across Malaysia but neither of these mutations were found in Ae. albopictus. Additionally, signatures of selection were detected on the Voltage-gated sodium channel gene in Ae. aegypti but not in Ae. albopictus. The presence of the 1534C allele was significantly associated with pyrethroid resistance and an additive effect to pyrethroid resistance was observed in individuals containing both kdr alleles.

Conclusions

Findings from this study will help to design and implement successful insecticide-based interventions against Ae. aegypti and Ae. albopictus to improve dengue control across Malaysia.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0797-2) contains supplementary material, which is available to authorized users.

Novel mutations on the ace-1 gene of the malaria vector Anopheles albimanus provide evidence for balancing selection in an area of high insecticide resistance in Peru

Background

Resistance to multiple classes of insecticides has been detected in the malaria vector Anopheles albimanus in northwest Peru. Acetylcholinesterase (AChE) insensitivity has previously been associated with resistance to organophosphate (OP) and carbamate (CA) insecticides in arthropods. A single point mutation on the ace-1 gene (G119S) associated with resistance to OPs and CAs has been described previously in four anopheline species, but not in field-collected An. albimanus. The present study aimed to characterize the role of ace-1 in conferring resistance to both OPs and CAs in the An. albimanus population in Tumbes, Peru.

Methods

The frequency and intensity of resistance to OPs and CAs was quantified through bioassays of female An. albimanus collected between 2012 and 2014, and the presence of insensitive AChE was confirmed using biochemical assays. A portion of the ace-1 gene flanking codon 119 was amplified and sequenced from individuals used in the bioassays and biochemical assays, as well as from historical samples collected in 2008. Statistical analyses were conducted to determine: (1) associations between genotype and AChE insensitivity; and, (2) associations between genotype and resistance phenotype.

Results

After confirming high levels of resistance to fenitrothion, malathion, and bendiocarb through bioassays, two novel polymorphisms were identified at the first and second loci of codon 119, with all individuals from the 2012–2014 collections being heterozygous at the first base (G/T) and either heterozygous (G/C) or homozygous mutants (C/C) at the second base. Based on sequence data from historical samples, these mutations arose prior to 2008, but became fixed in the population between 2008 and 2012. Homozygotes at the second locus had significantly higher levels of AChE insensitivity than heterozygotes (p <0.05). Individuals phenotypically susceptible to OPs and CAs were more likely to be heterozygous at the second locus (p <0.01). Cloning identified four individuals each containing three distinct genotypes, suggesting that a duplication of the ace-1 gene may have occurred.

Conclusions

The occurrence of heterozygotes at two loci and the presence of three genotypes in four individuals suggest that balancing selection could be maintaining OP and CA resistance in this population, while minimizing associated fitness costs.

Distribution and dissemination of the Val1016Ile and Phe1534Cys Kdr mutations in Aedes aegypti Brazilian natural populations

BACKGROUND

The chemical control of the mosquito Aedes aegypti, the major vector of dengue, is being seriously threatened due to the development of pyrethroid resistance. Substitutions in the 1016 and 1534 sites of the voltage gated sodium channel (AaNaV), commonly known as kdr mutations, confer the mosquito with knockdown resistance. Our aim was to evaluate the allelic composition of natural populations of Brazilian Ae. aegypti at both kdr sites.

METHODS

The AaNaV IIIS6 region was cloned and sequenced from three Brazilian populations. Additionally, individual mosquitoes from 30 populations throughout the country were genotyped for 1016 and 1534 sites, based in allele-specific PCR. For individual genotypes both sites were considered as a single locus.

RESULTS

The 350 bp sequence spanning the IIIS6 region of the AaNaV gene revealed the occurrence of the kdr mutation Phe1534Cys in Brazil. Concerning the individual genotyping, beyond the susceptible wild-type (NaVS), two kdr alleles were identified: substitutions restricted to the 1534 position (NaVR1) or simultaneous substitutions in both 1016 and 1534 sites (NaVR2). A clear regional distribution pattern of these alleles was observed. The NaVR1kdr allele occurred in all localities, while NaVR2 was more frequent in the Central and Southeastern localities. Locations that were sampled multiple times in the course of a decade revealed an increase in frequency of the kdr mutations, mainly the double mutant allele NaVR2. Recent samples also indicate that NaVR2 is spreading towards the Northern region.

CONCLUSIONS

We have found that in addition to the previously reported Val1016Ile kdr mutation, the Phe1534Cys mutation also occurs in Brazil. Allelic composition at both sites was important to elucidate the actual distribution of kdr mutations throughout the country. Studies to determine gene flow and the fitness costs of these kdr alleles are underway and will be important to better understand the dynamics of Ae. aegypti pyrethroid resistance.