Abstract
Aedes aegypti is a major vector of Zika, dengue, and other arboviruses. Permethrin adulticidal spraying, which targets the voltage-gated sodium channel (VGSC), is commonly done to reduce local mosquito populations and protect humans from exposure to arbovirus pathogens transmitted by this dangerous pest. Permethrin resistance, however, is a growing problem and understanding its underlying molecular basis may identify avenues to combat it. We identified a single G:C polymorphism in pre-miR-33 that was genetically associated with permethrin resistance; resulting isoforms had structural differences that may affect DICER-1/pre-miRNA processing rates. We then assessed the effects of overexpression of pre-miR-33 isoforms on permethrin toxicological phenotypes, VGSC transcript abundance and protein levels for two genetically related mosquito strains. One strain had its naturally high permethrin resistance levels maintained by periodic treatment, and the other was released from selection. VGSC protein levels were lower in the permethrin resistant strain than in the related permethrin-susceptible strain. Overexpression of the G-pre-miR-33 isoform reduced VGSC expression levels in both strains. To further elucidate changes in gene expression associated with permethrin resistance, exome-capture gDNA deep sequencing, genetic association mapping and subsequent gene set enrichment analysis revealed that transport genes, in particular, were selected in resistant versus susceptible mosquitoes. Collectively, these data indicate that miR-33 regulates VGSC expression as part of a nuanced system of neuronal regulation that contributes to a network of heritable features determining permethrin resistance.