Transcriptomic analysis of insecticide resistance in the lymphatic filariasis vector Culex quinquefasciatus.

Culex quinquefasciatus plays an important role in transmission of vector-borne diseases of public health importance, including lymphatic filariasis (LF), as well as many arboviral diseases. Currently, efforts to tackle C. quinquefasciatus vectored diseases are based on either mass drug administration (MDA) for LF, or insecticide-based interventions. Widespread and intensive insecticide usage has resulted in increased resistance in mosquito vectors, including C. quinquefasciatus. Herein, the transcriptome profile of Ugandan bendiocarb-resistant C. quinquefasciatus was explored to identify candidate genes associated with insecticide resistance. High levels of insecticide resistance were observed for five out of six insecticides tested, with the lowest mortality (0.97%) reported to permethrin, while for DDT, lambdacyhalothrin, bendiocarb and deltamethrin the mortality rate ranged from 1.63-3.29%. Resistance to bendiocarb in exposed mosquitoes was marked, with 2.04% mortality following 1 h exposure and 58.02% after 4 h. Genotyping of the G119S Ace-1 target site mutation detected a highly significant association (p < 0.0001; OR = 25) between resistance and Ace1-119S. However, synergist assays using the P450 inhibitor PBO, or the esterase inhibitor TPP resulted in markedly increased mortality (to ≈80%), suggesting a role of metabolic resistance in the resistance phenotype. Using a novel, custom 60 K whole-transcriptome microarray 16 genes significantly overexpressed in resistant mosquitoes were detected, with the P450 Cyp6z18 showing the highest differential gene expression (>8-fold increase vs unexposed controls). These results provide evidence that bendiocarb resistance in Ugandan C. quinquefasciatus is mediated by both target-site mechanisms and over-expression of detoxification enzymes.

Development and application of a tri-allelic PCR assay for screening Vgsc-L1014F kdr mutations associated with pyrethroid and organochlorine resistance in the mosquito Culex quinquefasciatus.

BACKGROUND: Culex quinquefasciatus has a widespread distribution across tropical and sub-tropical regions, and plays an important role in the transmission of vector-borne diseases of public health importance, including lymphatic filariasis (LF) and multiple arboviruses. Increased resistance to insecticides threatens the efficacy and sustainability of insecticide-based anti-vector interventions which mitigate the burden of mosquito transmitted diseases in endemic regions. In C. quinquefasciatus two non-synonymous voltage gated sodium channel (Vgsc) variants, both resulting in a leucine to phenylalanine change at codon 1014, are associated with resistance to pyrethroids and DDT. This tri-allelic variation has compromised the ability to perform high-throughput single-assay screening. To facilitate the detection and monitoring of the Vgsc-1014 locus in field-caught mosquitoes, an Engineered-Tail Allele-Specific-PCR (ETAS-PCR) diagnostic assay was developed and applied to wild mosquitoes from Brazil, Tanzania and Uganda. RESULTS: This new cost-effective, single-tube assay was compared to two, well-established, genotyping approaches, pyrosequencing and TaqMan. The ETAS-PCR assay showed high specificity for discriminating the three alleles at Vgsc-L1014F, with genotyping results strongly correlated with pyrosequencing and TaqMan results (98.64% and 100% respectively). CONCLUSIONS: Our results support the utility of the ETAS-PCR/Vgsc-1014 diagnostic assay, which stands as an effective alternative for genotyping tri-allelic variants.

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.

Insecticide resistance in Aedes aegypti populations from Ceará, Brazil.

BACKGROUND: Organophosphates and pyrethroids are used widely in Brazil to control Aedes aegypti, the main vector of dengue viruses, under the auspices of the National Programme for Dengue Control. Resistance to these insecticides is widespread throughout Brazil. In Ceará the vector is present in 98% of districts and resistance to temephos has been reported previously. Here we measure resistance to temephos and the pyrethroid cypermethrin in three populations from Ceará and use biochemical and molecular assays to characterise resistance mechanisms. RESULTS: Resistance to temephos varied widely across the three studied populations, with resistance ratios (RR(95)) of 7.2, 30 and 192.7 in Juazeiro do Norte, Barbalha and Crato respectively. The high levels of resistance detected in Barbalha and Crato (RR(95) ≥ 30) imply a reduction of temephos efficacy, and indeed in simulated field tests reduced effectiveness was observed for the Barbalha population. Two populations (Crato and Barbalha) were also resistant to cypermethrin, whilst Juazeiro do Norte showed only an altered susceptibility. The Ile1011Met kdr mutation was detected in all three populations and Val1016Ile in Crato and Juazeiro do Norte. 1011Met was significantly associated with resistance to cypermethrin in the Crato population. Biochemical tests showed that only the activity of esterases and GSTs, among the tested detoxification enzymes, was altered in these populations when compared with the Rockefeller strain. CONCLUSIONS: Our results demonstrate that two A. aegypti populations from Ceará are under strong selection pressure by temephos, compromising the field effectiveness of this organophosphate. Our results also provide evidence that the process of reducing resistance to this larvicide in the field is difficult and slow and may require more than seven years for reversal. In addition, we show resistance to cypermethrin in two of the three populations studied, and for the first time the presence of the allele 1016Ile in mosquito populations from northeastern Brazil. A significant association between 1011Met and resistance was observed in one of the populations. Target-site mechanisms seem not to be implicated in temephos resistance, reinforcing the idea that for the studied populations, detoxification enzymes most likely play a major role in the resistance to this insecticide.