Background Resistance to major public health insecticides in Côte d’Ivoire has intensified and now threatens the long-term effectiveness of malaria vector control interventions. Methods This study evaluated the bioefficacy of conventional and next-generation long-lasting insecticidal nets (LLINs), determined resistance profiles, and characterized molecular and metabolic mechanisms in wild Anopheles coluzzii from South-East Côte d’Ivoire in 2019. Results Phenotypic resistance was intense: more than 25% of mosquitoes survived exposure to ten times the doses of pyrethroids required to kill susceptible populations. Similarly, 24-hour mortality to deltamethrin-only LLINs was very low and not significantly different to an untreated net. Sub-lethal pyrethroid exposure did not induce significant delayed vector mortality 72 hours later. In contrast, LLINs containing the synergist piperonyl butoxide (PBO), or new insecticides, clothianidin and chlorfenapyr, were highly toxic to An. coluzzii. Pyrethroid-susceptible An. coluzzii were significantly more likely to be infected with malaria, compared to those that survived insecticidal exposure. Pyrethroid resistance was associated with significant over-expression of CYP6P4, CPY6Z1 and CYP6P3. Conclusions Study findings raise concerns regarding the operational failure of standard LLINs and support the urgent deployment of vector control interventions incorporating PBO, chlorfenapyr or clothianidin in areas of high resistance intensity in Côte d’Ivoire.

Insecticide resistance among mosquito species is now a pervasive phenomenon that threatens to jeopardize global malaria vector control efforts. Evidence of links between the mosquito microbiota and insecticide resistance is emerging, with significant enrichment of insecticide degrading bacteria and enzymes in resistant populations. Using 16S rRNA amplicon sequencing, we characterized and compared the microbiota of Anopheles coluzzii in relation to their deltamethrin resistance and exposure profiles. Comparisons between 2- and 3-day-old deltamethrin-resistant and -susceptible mosquitoes demonstrated significant differences in microbiota diversity. Ochrobactrum, Lysinibacillus, and Stenotrophomonas genera, each of which comprised insecticide-degrading species, were significantly enriched in resistant mosquitoes. Susceptible mosquitoes had a significant reduction in alpha diversity compared to resistant individuals, with Asaia and Serratia dominating microbial profiles. There was no significant difference in deltamethrin-exposed and -unexposed 5- to 6-day-old individuals, suggesting that insecticide exposure had minimal impact on microbial composition. Serratia and Asaia were also dominant in 5- to 6-day-old mosquitoes, which had reduced microbial diversity compared to 2- to 3-day-old mosquitoes. Our findings revealed significant alterations of Anopheles coluzzii microbiota associated with deltamethrin resistance, highlighting the potential for identification of novel microbial markers for insecticide resistance surveillance. qPCR detection of Serratia and Asaia was consistent with 16S rRNA sequencing, suggesting that population-level field screening of bacterial microbiota may be feasibly integrated into wider resistance monitoring, if reliable and reproducible markers associated with phenotype can be identified. IMPORTANCE Control of insecticide-resistant vector populations remains a significant challenge to global malaria control and while substantial progress has been made elucidating key target site mutations, overexpressed detoxification enzymes and alternate gene families, the contribution of the mosquito microbiota to phenotypic insecticide resistance has been largely overlooked. We focused on determining the effects of deltamethrin resistance intensity on Anopheles coluzzii microbiota and identifying any microbial taxa associated with phenotype. We demonstrated a significant reduction in microbial diversity between deltamethrin-resistant and -susceptible mosquitoes. Insecticide degrading bacterial species belonging to Ochrobactrum, Lysinibacillus, and Stenotrophomonas genera were significantly enriched in resistant mosquitoes, while Asaia and Serratia dominated microbial profiles of susceptible individuals. Our results revealed significant alterations of Anopheles coluzzii microbiota associated with deltamethrin resistance, highlighting the potential for identification of novel microbial markers for surveillance and opportunities for designing innovative control techniques to prevent the further evolution and spread of insecticide resistance.

Surveillance of malaria vector species and the monitoring of insecticide resistance are essential to inform malaria control strategies and support the reduction of infections and disease. Genetic barcoding of mosquitoes is a useful tool to assist the high-throughput surveillance of insecticide resistance, discriminate between sibling species and to detect the presence of Plasmodium infections. In this study, we combined multiplex PCR, custom designed dual indexing, and Illumina next generation sequencing for high throughput single nucleotide polymorphism (SNP)-profiling of four species from the Anopheles (An.) gambiae complex ( An. gambiae sensu stricto, An. coluzzii, An. arabiensis and An. melas ). By amplifying and sequencing only 14 genetic fragments (500 bp each), we were able to simultaneously detect Plasmodium infection; insecticide resistance-conferring SNPs in ace1 , gste 2, vgsc and rdl genes; the partial sequences of nuclear ribosomal internal transcribed spacers (ITS1 and ITS2) and intergenic spacers (IGS), Short INterspersed Elements (SINE), as well as mitochondrial genes ( cox1 and nd4 ) for species identification and genetic diversity. Using this amplicon sequencing approach with the four selected An. gambiae complex species, we identified a total of 15 non-synonymous mutations in the insecticide target genes, including previously described mutations associated with resistance and two new mutations (F1525L in vgsc and D148E in gste2 ). Overall, we present a reliable and cost-effective high-throughput panel for surveillance of An. gambiae complex mosquitoes in malaria endemic regions.