METHODS: Anopheles gambiae (s.l.) mosquitoes were identified to species level using PCR techniques. Standard WHO insecticide susceptibility bioassays were carried out to detect resistance to deltamethrin (0.05%), DDT (4%) and bendiocarb (0.1%). TaqMan assays were performed on random samples of deltamethrin-resistant phenotyped and pyrethrum spray collected individuals to determine Vgsc-1014 knockdown resistance mutations.
RESULTS: Anopheles arabiensis accounted for 99.9% of any anopheline species collected across all sites. Bioassay screening indicated that mosquitoes remained susceptible to bendiocarb but were resistance to deltamethrin and DDT in all areas. There were significant increases in deltamethrin resistance over the four years, with overall mean percent mortality to deltamethrin declining from 81.0% (95% CI: 77.6-84.3%) in 2011 to 47.7% (95% CI: 43.5-51.8%) in 2014. The rate of increase in phenotypic deltamethrin-resistance was significantly slower in the LLIN + IRS arm than in the LLIN-only arm (Odds ratio 1.34; 95% CI: 1.02-1.77). The frequency of Vgsc-1014F mutation varied spatiotemporally with highest frequencies in Galabat (range 0.375-0.616) and New Halfa (range 0.241-0.447). Deltamethrin phenotypic-resistance correlated with Vgsc-1014F frequency.
CONCLUSION: Combining LLIN and IRS, with different classes of insecticide, may delay pyrethroid resistance development, but the speed at which resistance develops may be area-specific. Continued monitoring is vital to ensure optimal management and control.
METHODOLOGY/PRINCIPAL FINDINGS: Genome-wide microarray-based transcription analysis was carried out to detect the genes associated with metabolic resistance in these populations. Comparisons of the susceptible New Orleans strain to three non-exposed multiple insecticide resistant field strains; Penang, Kuala Lumpur and Kota Bharu detected 2605, 1480 and 425 differentially expressed transcripts respectively (fold-change>2 and p-value ≤ 0.05). 204 genes were commonly over-expressed with monooxygenase P450 genes (CYP9J27, CYP6CB1, CYP9J26 and CYP9M4) consistently the most up-regulated detoxification genes in all populations, indicating that they possibly play an important role in the resistance. In addition, glutathione S-transferases, carboxylesterases and other gene families commonly associated with insecticide resistance were also over-expressed. Gene Ontology (GO) enrichment analysis indicated an over-representation of GO terms linked to resistance such as monooxygenases, carboxylesterases, glutathione S-transferases and heme-binding. Polymorphism analysis of CYP9J27 sequences revealed a high level of polymorphism (except in Joho Bharu), suggesting a limited directional selection on this gene. In silico analysis of CYP9J27 activity through modelling and docking simulations suggested that this gene is involved in the multiple resistance in Malaysian populations as it is predicted to metabolise pyrethroids, DDT and bendiocarb.
CONCLUSION/SIGNIFICANCE: The predominant over-expression of cytochrome P450s suggests that synergist-based (PBO) control tools could be utilised to improve control of this major dengue vector across Malaysia.