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.

BackgroundThe capacity of Aedes mosquitoes to resist chemical insecticides threatens the control of major arbovirus diseases worldwide. Until alternative control tools are widely deployed, monitoring insecticide resistance levels and identifying resistance mechanisms in field mosquito populations is crucial for implementing appropriate management strategies. Metabolic resistance to pyrethroids is common in Aedes aegypti but the monitoring of the dynamics of resistant alleles is impeded by the lack of robust genomic markers.Methodology/Principal findingsIn an attempt to identify the genomic bases of metabolic resistance to deltamethrin, multiple resistant and susceptible populations originating from various continents were compared using both RNA-seq and a targeted DNA-seq approach focused on the upstream regions of detoxification genes. Multiple detoxification enzymes were over transcribed in resistant populations, frequently associated with an increase in their gene copy number. Targeted sequencing identified potential promoter variations associated with their over transcription. Non-synonymous variations affecting detoxification enzymes were also identified in resistant populations.Conclusion /SignificanceThis study not only confirmed the role of gene copy number variations as a frequent cause of the over expression of detoxification enzymes associated with insecticide resistance in Aedes aegypti but also identified novel genomic resistance markers potentially associated with their cis-regulation and modifications of their protein structure conformation. As for gene transcription data, polymorphism patterns were frequently conserved within regions but differed among continents confirming the selection of different resistance factors worldwide. Overall, this study paves the way of the identification of a comprehensive set of genomic markers for monitoring the spatio-temporal dynamics of the variety of insecticide resistance mechanisms in Aedes aegypti.

BackgroundIn the Guadeloupe and Saint Martin islands, Aedes aegypti mosquitoes are the only recognized vectors of dengue, chikungunya, and Zika viruses. For around 40 years, malathion was used as a mosquito adulticide and temephos as a larvicide. Since the European Union banned the use of these two insecticide molecules in the first decade of the 21st century, deltamethrin and Bacillus thuringiensis var. israelensis are the remaining adulticide and larvicide, respectively, used in Guadeloupe. In order to improve the management of vector control activities in Guadeloupe and Saint Martin, we investigated Ae. aegypti resistance to and mechanisms associated with deltamethrin, malathion, and temephos.Methods Ae. aegypti mosquitoes were collected from six different localities of Guadeloupe and Saint Martin. Larvae were used for malathion and temephos bioassays, and adult mosquitoes for deltamethrin bioassays, following World Health Organization recommendations. Knockdown resistance (Kdr) genotyping for V1016I and F1534C mutations, and expression levels of eight enzymes involved in detoxification mechanisms were examined in comparison with the susceptible reference Bora Bora strain.ResultsResistance ratios (RR50) calculated for Ae. aegypti larvae showed high resistance levels to temephos (from 8.9 to 33.1-fold) and low resistance levels to malathion (from 1.7 to 4.4-fold). Adult females displayed moderate resistance levels to deltamethrin regarding the time necessary to affect 50% of individuals, varying from 8.0 to 28.1-fold. Molecular investigations on adult mosquitoes showed high resistant allele frequencies for V1016I and F1534C (from 85 to 96% and from 90 to 98%, respectively), as well as an overexpression of the glutathione S-transferase gene, GSTe2, the carboxylesterase CCEae3a, and the cytochrome genes 014614, CYP6BB2, CYP6M11, and CYP9J23.Conclusions Ae. aegypti populations from Guadeloupe and Saint Martin exhibit multiple resistance to organophosphates (temephos and malathion), and pyrethroids (deltamethrin). The mechanisms associated with these resistance patterns show strong frequencies of F1534C and V1016I Kdr mutations, and an over-expression of CCEae3a, GSTe2, and four cytochrome P450 genes (014614, CYP9J23, CYP6M11, CYP6BB2). These results will form the baseline for a deeper understanding of the insecticide resistance levels and associated mechanisms of Ae. aegypti populations and will be used to improve vector control strategies in Guadeloupe and Saint Martin.Electronic supplementary materialThe online version of this article (doi:10.1186/s40249-017-0254-x) contains supplementary material, which is available to authorized users.