Genetic variation at the Cyp6m2 putative insecticide resistance locus in Anopheles gambiae and Anopheles coluzzii

Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control

Insecticide resistance in malaria vector populations poses a major threat to malaria control, which relies largely on insecticidal interventions. Contemporary vector-control strategies focus on combatting resistance using multiple insecticides with differing modes of action within the mosquito. However, diverse genetic resistance mechanisms are present in vector populations, and continue to evolve. Knowledge of the spatial distribution of these genetic mechanisms, and how they impact the efficacy of different insecticidal products, is critical to inform intervention deployment decisions. We developed a catalogue of genetic-resistance mechanisms in African malaria vectors that could guide molecular surveillance. We highlight situations where intervention deployment has led to resistance evolution and spread, and identify challenges in understanding and mitigating the epidemiological impacts of resistance.

The impact of agrochemical pollutant mixtures on the selection of insecticide resistance in the malaria vector Anopheles gambiae: insights from experimental evolution and transcriptomics

BACKGROUND

There are several indications that pesticides used in agriculture contribute to the emergence and spread of resistance of mosquitoes to vector control insecticides. However, the impact of such an indirect selection pressure has rarely been quantified and the molecular mechanisms involved are still poorly characterized. In this context, experimental selection with different agrochemical mixtures was conducted in Anopheles gambiae. The multi-generational impact of agrochemicals on insecticide resistance was evaluated by phenotypic and molecular approaches.

METHODS

Mosquito larvae were selected for 30 generations with three different agrochemical mixtures containing (i) insecticides, (ii) non-insecticides compounds, and (iii) both insecticide and non-insecticide compounds. Every five generations, the resistance of adults to deltamethrin and bendiocarb was monitored using bioassays. The frequencies of the kdr (L995F) and ace1 (G119S) target-site mutations were monitored every 10 generations. RNAseq was performed on all lines at generation 30 in order to identify gene transcription level variations and polymorphisms associated with each selection regime.

RESULTS

Larval selection with agrochemical mixtures did not affect bendiocarb resistance and did not select for ace1 mutation. Contrastingly, an increased deltamethrin resistance was observed in the three selected lines. Such increased resistance was not majorly associated with the presence of kdr L995F mutation in selected lines. RNA-seq identified 63 candidate resistance genes over-transcribed in at least one selected line. These include genes coding for detoxification enzymes or cuticular proteins previously associated with insecticide resistance, and other genes potentially associated with chemical stress response. Combining an allele frequency filtering with a Bayesian FST-based genome scan allowed to identify genes under selection across multiple genomic loci, supporting a multigenic adaptive response to agrochemical mixtures.

CONCLUSION

This study supports the role of agrochemical contaminants as a significant larval selection pressure favouring insecticide resistance in malaria vectors. Such selection pressures likely impact kdr mutations and detoxification enzymes, but also more generalist mechanisms such as cuticle resistance, which could potentially lead to cross-tolerance to unrelated insecticide compounds. Such indirect effect of global landscape pollution on mosquito resistance to public health insecticides deserves further attention since it can affect the nature and dynamics of resistance alleles circulating in malaria vectors and impact the efficacy of control vector strategies.

Genetic Diversity of Cytochrome P450s CYP6M2 and CYP6P4 Associated with Pyrethroid Resistance in the Major Malaria Vectors Anopheles coluzzii and Anopheles gambiae from Yaoundé, Cameroon

Assessing the genetic diversity of metabolic resistance genes, such as cytochrome P450s, helps to understand the dynamics and evolution of resistance in the field. Here, we analyzed the polymorphisms of CYP6M2 and CYP6P4, associated with pyrethroid resistance in Anopheles coluzzii and Anopheles gambiae, to detect potential resistance markers. Field-caught resistant mosquitos and susceptible lab strains were crossed, and F4 was exposed to permethrin for 15 min and 90 min to discriminate highly susceptible (HS) and highly resistant (HR) mosquitos, respectively. Significant permethrin mortality reduction was observed after pre-exposure to PBO, suggesting the gene involvement of P450s. qPCR analysis revealed significant overexpression of CYP6M2 (FC = 19.57 [95% CI 13.96-25.18] for An. coluzzii; 10.16 [7.86-12.46] for An. gambiae) and CYP6P4 (FC = 6.73 [6.15-7.30] An. coluzzii; 23.62 [26.48-20.76] An. gambiae). Full-gene and ≈1 kb upstream were sequenced. For CYP6M2, the upstream region shows low diversity in HR and HS (overall Hd = 0.49, π = 0.018), whereas the full-gene shows allelic-variation but without evidence of ongoing selection. CYP6P4 upstream region showed a lower diversity in HR (Hd = 0.48) than HS (Hd = 0.86) of An. gambiae. These results highlighted that CYP6P4-associated resistance is potentially driven by modification in upstream region. However, further work is needed to determine the real causative variants that will help design rapid detection tools.

The cytochrome P450 CYP325A is a major driver of pyrethroid resistance in the major malaria vector Anopheles funestus in Central Africa

The overexpression and overactivity of key cytochrome P450s (CYP450) genes are major drivers of metabolic resistance to insecticides in African malaria vectors such as Anopheles funestus s.s. Previous RNAseq-based transcription analyses revealed elevated expression of CYP325A specific to Central African populations but its role in conferring resistance has not previously been demonstrated. In this study, RT-qPCR consistently confirmed that CYP325A is highly over-expressed in pyrethroid-resistant An. funestus from Cameroon, compared with a control strain and insecticide-unexposed mosquitoes. A synergist bioassay with PBO significantly recovered susceptibility for permethrin and deltamethrin indicating P450-based metabolic resistance. Analyses of the coding sequence of CYP325A Africa-wide detected high-levels of polymorphism, but with no predominant alleles selected by pyrethroid resistance. Geographical amino acid changes were detected notably in Cameroon. In silico homology modelling and molecular docking simulations predicted that CYP325A binds and metabolises type I and type II pyrethroids. Heterologous expression of recombinant CYP325A and metabolic assays confirmed that the most-common Cameroonian haplotype metabolises both type I and type II pyrethroids with depletion rate twice that the of the DR Congo haplotype. Analysis of the 1 kb putative promoter of CYP325A revealed reduced diversity in resistant mosquitoes compared to susceptible ones, suggesting a potential selective sweep in this region. The establishment of CYP325A as a pyrethroid resistance metabolising gene further explains pyrethroid resistance in Central African populations of An. funestus. Our work will facilitate future efforts to detect the causative resistance markers in the promoter region of CYP325A to design field applicable DNA-based diagnostic tools.