Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment

Phenotypic evidence of deltamethrin resistance and identification of selective sweeps in Anopheles mosquitoes on the Bijagós Archipelago, Guinea-Bissau

Vector control in the Bijagós Archipelago of Guinea-Bissau currently relies on pyrethroid insecticide-treated nets. However, data on insecticide resistance in Guinea-Bissau is limited. This study identified deltamethrin resistance in the Anopheles gambiae sensu lato complex on Bubaque island using WHO tube tests in November 2022. Whole genome sequencing of An. gambiae sensu stricto mosquitoes identified six single nucleotide polymorphisms (SNPs) previously associated with, or putatively associated with, insecticide resistance: T791M, L995F, N1570Y, A1746S and P1874L in the vgsc gene, and L119V in the gste2 gene. Twenty additional non-synonymous SNPs were identified in insecticide-resistance associated genes. Four of these SNPs were present at frequencies over 5% in the population: T154S, I126F and G26S in the vgsc gene and A65S in ace1. Genome wide selection scans using Garud's H12 statistic identified two selective sweeps: one in chromosome X and one in chromosome 2R. Both selective sweeps overlap with metabolic genes previously associated with insecticide resistance, including cyp9k1 and the cyp6aa/cyp6p gene cluster. This study presents the first phenotypic testing for deltamethrin resistance and the first whole genome sequence data for Anophelesgambiae mosquitoes from the Bijagós, contributing data of significance for vector control policy in this region.

Cryptic population structure and insecticide resistance in Anopheles gambiae from the southern Democratic Republic of Congo

The Democratic Republic of Congo (DRC) suffers from one of the highest malaria burdens worldwide, but information on its Anopheles vector populations is relatively limited. Preventative malaria control in DRC is reliant on pyrethroid-treated nets, raising concerns over the potential impacts of insecticide resistance. We sampled Anopheles gambiae from three geographically distinct populations (Kimpese, Kapolowe and Mikalayi) in southern DRC, collecting from three sub-sites per population and characterising mosquito collections from each for resistance to pyrethroids using WHO tube bioassays. Resistance to each of three different pyrethroids was generally high in An. gambiae with < 92% mortality in all tests, but varied between collections, with mosquitoes from Kimpese being the most resistant. Whole genome sequencing of 165 An. gambiae revealed evidence for genetic differentiation between Kimpese and Kapolowe/Mikalayi, but not between the latter two sample sites despite separation of approximately 800 km. Surprisingly, there was evidence of population structure at a small spatial scale between collection subsites in Kimpese, despite separation of just tens of kilometres. Intra-population (H12) and inter-population (FST) genome scans identified multiple peaks corresponding to genes associated with insecticide resistance such as the voltage gated sodium channel (Vgsc) target site on chromosome 2L, a Cyp6 cytochrome P450 cluster on chromosome arm 2R, and the Cyp9k1 P450 gene on chromosome X. In addition, in the Kimpese subsites, the P450 redox partner gene Cpr showed evidence for contemporary selection (H12) and population differentiation (FST) meriting further exploration as a potential resistance associated marker.

Parallel Evolution in Mosquito Vectors-A Duplicated Esterase Locus is Associated With Resistance to Pirimiphos-methyl in Anopheles gambiae

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programs. The organophosphate (OP), pirimiphos-methyl, is a relatively new chemical in the vector control armory but is now widely used in indoor-residual spray campaigns. While generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in A. gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to OPs in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in A. gambiae, A. coluzzii and A. arabiensis. As in C. pipiens, copy number variants have arisen at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in A. gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programs.

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.

Association of a rapidly selected 4.3kb transposon-containing structural variation with a P450-based resistance to pyrethroids in the African malaria vector Anopheles funestus

Deciphering the evolutionary forces controlling insecticide resistance in malaria vectors remains a prerequisite to designing molecular tools to detect and assess resistance impact on control tools. Here, we demonstrate that a 4.3kb transposon-containing structural variation is associated with pyrethroid resistance in central/eastern African populations of the malaria vector Anopheles funestus. In this study, we analysed Pooled template sequencing data and direct sequencing to identify an insertion of 4.3kb containing a putative retro-transposon in the intergenic region of two P450s CYP6P5-CYP6P9b in mosquitoes of the malaria vector Anopheles funestus from Uganda. We then designed a PCR assay to track its spread temporally and regionally and decipher its role in insecticide resistance. The insertion originates in or near Uganda in East Africa, where it is fixed and has spread to high frequencies in the Central African nation of Cameroon but is still at low frequency in West Africa and absent in Southern Africa. A marked and rapid selection was observed with the 4.3kb-SV frequency increasing from 3% in 2014 to 98% in 2021 in Cameroon. A strong association was established between this SV and pyrethroid resistance in field populations and is reducing pyrethroid-only nets' efficacy. Genetic crosses and qRT-PCR revealed that this SV enhances the expression of CYP6P9a/b but not CYP6P5. Within this structural variant (SV), we identified putative binding sites for transcription factors associated with the regulation of detoxification genes. An inverse correlation was observed between the 4.3kb SV and malaria parasite infection, indicating that mosquitoes lacking the 4.3kb SV were more frequently infected compared to those possessing it. Our findings highlight the underexplored role and rapid spread of SVs in the evolution of insecticide resistance and provide additional tools for molecular surveillance of insecticide resistance.

LLIN Evaluation in Uganda Project (LLINEUP)-effects of a vector control trial on Plasmodium infection prevalence and genotypic markers of insecticide resistance in Anopheles vectors from 48 districts of Uganda

Pyrethroid bednets treated with the synergist piperonyl butoxide (PBO) offer the possibility of improved vector control in mosquito populations with metabolic resistance. In 2017-2019, we conducted a large-scale, cluster-randomised trial (LLINEUP) to evaluate long-lasting insecticidal nets (LLINs) treated with a pyrethroid insecticide plus PBO (PBO LLINs), as compared to conventional, pyrethroid-only LLINs across 104 health sub-districts (HSDs) in Uganda. In LLINEUP, and similar trials in Tanzania, PBO LLINs were found to provide greater protection against malaria than conventional LLINs, reducing parasitaemia and vector density. In the LLINEUP trial, we conducted cross-sectional household entomological surveys at baseline and then every 6 months for two years, which we use here to investigate longitudinal changes in mosquito infection rate and genetic markers of resistance. Overall, 5395 female Anopheles mosquitoes were collected from 5046 households. The proportion of mosquitoes infected (PCR-positive) with Plasmodium falciparum did not change significantly over time, while infection with non-falciparum malaria decreased in An. gambiae s.s., but not An. funestus. The frequency of genetic markers associated with pyrethroid resistance increased significantly over time, but the rate of change was not different between the two LLIN types. The knock-down resistance (kdr) mutation Vgsc-995S declined over time as Vgsc-995F, the alternative resistance mutation at this codon, increased. Vgsc-995F appears to be spreading into Uganda. Distribution of LLINs in Uganda was previously found to be associated with reductions in parasite prevalence and vector density, but here we show that the proportion of infective mosquitoes remained stable across both PBO and non-PBO LLINs, suggesting that the potential for transmission persisted. The increased frequency of markers of pyrethroid resistance indicates that LLIN distribution favoured the evolution of resistance within local vectors and highlights the potential benefits of resistance management strategies.Trial registration: This study is registered with ISRCTN, ISRCTN17516395. Registered 14 February 2017, http://www.isrctn.com/ISRCTN17516395 .

Pyrethroid resistance and gene expression profile of a new resistant An. gambiae colony from Uganda reveals multiple resistance mechanisms and overexpression of Glutathione-S-Transferases linked to survival of PBO-pyrethroid combination

Background

The effectiveness of long-lasting insecticidal nets (LLINs) are being threatened by growing resistance to pyrethroids. To restore their efficacy, a synergist, piperonyl butoxide (PBO) which inhibits cytochrome P450s has been incorporated into pyrethroid treated nets. A trial of PBO-LLINs was conducted in Uganda from 2017 and we attempted to characterize mechanisms of resistance that could impact intervention efficacy.

Methods

We established an Anopheles gambiae s.s colony in 2018 using female mosquitoes collected from Busia district in eastern Uganda. We first assessed the phenotypic resistance profile of this colony using WHO tube and net assays using a deltamethrin dose-response approach. The Busia colony was screened for known resistance markers and RT-qPCR targeting 15 genes previously associated with insecticide resistance was performed.

Results

The Busia colony had very high resistance to deltamethrin, permethrin and DDT. In addition, the colony had moderate resistance to alpha-cypermethrin and lambda-cyhalothrin but were fully susceptible to bendiocarb and fenitrothion. Exposure to PBO in combination with permethrin and deltamethrin resulted in higher mortality rates in both net and tube assays, with a higher mortality observed in net assays than tube assays. The kdr marker, Vgsc-995S was at very high frequency (91.7-98.9%) whilst the metabolic markers Coeae1d and Cyp4j5-L43F were at very low (1.3% - 11.5%) and moderate (39.5% - 44.7%) frequencies respectively. Our analysis showed that gene expression pattern in mosquitoes exposed to deltamethrin, permethrin or DDT only were similar in comparison to the susceptible strain and there was significant overexpression of cytochrome P450s, glutathione-s-transferases (GSTs) and carboxyl esterases (COEs). However, mosquitoes exposed to both PBO and pyrethroid strikingly and significantly only overexpressed closely related GSTs compared to unexposed mosquitoes while major cytochrome P450s were underexpressed.

Conclusions

The high levels of pyrethroid resistance observed in Busia appears associated with a wide range of metabolic gene families.

Copy number variants underlie the major selective sweeps in insecticide resistance genes in Anopheles arabiensis from Tanzania

To keep ahead of the evolution of resistance to insecticides in mosquitoes, national malaria control programmes must make use of a range of insecticides, both old and new, while monitoring resistance mechanisms. Knowledge of the mechanisms of resistance remains limited in Anopheles arabiensis, which in many parts of Africa is of increasing importance because it is apparently less susceptible to many indoor control interventions. Furthermore, comparatively little is known in general about resistance to non-pyrethroid insecticides such as pirimiphos-methyl (PM), which are crucial for effective control in the context of resistance to pyrethroids. We performed a genome-wide association study to determine the molecular mechanisms of resistance to deltamethrin (commonly used in bednets) and PM, in An. arabiensis from two regions in Tanzania. Genomic regions of positive selection in these populations were largely driven by copy number variants (CNVs) in gene families involved in resistance to these two insecticides. We found evidence of a new gene cluster involved in resistance to PM, identifying a strong selective sweep tied to a CNV in the Coeae2g-Coeae6g cluster of carboxylesterase genes. Using complementary data from An. coluzzii in Ghana, we show that copy number at this locus is significantly associated with PM resistance. Similarly, for deltamethrin, resistance was strongly associated with a novel CNV allele in the Cyp6aa / Cyp6p cluster. Against this background of metabolic resistance, target site resistance was very rare or absent for both insecticides. Mutations in the pyrethroid target site Vgsc were at very low frequency in Tanzania, yet combining these samples with three An. arabiensis individuals from West Africa revealed a startling diversity of evolutionary origins of target site resistance, with up to 5 independent origins of Vgsc-995 mutations found within just 8 haplotypes. Thus, despite having been first recorded over 10 years ago, Vgsc resistance mutations in Tanzanian An. arabiensis have remained at stable low frequencies. Overall, our results provide a new copy number marker for monitoring resistance to PM in malaria mosquitoes, and reveal the complex picture of resistance patterns in An. arabiensis.

Genomic Profiling of Insecticide Resistance in Malaria Vectors: Insights into Molecular Mechanisms

Malaria control faces challenges from widespread insecticide resistance in major Anopheles species. This study, employing a cross-species approach, integrates RNA-Sequencing, whole-genome sequencing, and microarray data to elucidate drivers of insecticide resistance in Anopheles gambiae complex and An. funestus. Findings show an inverse relationship between genetic diversity and gene expression, with highly expressed genes experiencing stronger purifying selection. These genes cluster physically in the genome, revealing potential coordinated regulation. We identified known and novel candidate insecticide resistance genes, enriched in metabolic, cuticular, and behavioural functions. We also present AnoExpress, a Python package, and an online interface for user-friendly exploration of resistance candidate expression. Despite millions of years of speciation, convergent gene expression responses to insecticidal selection pressures are observed across Anopheles species, providing crucial insights for malaria vector control. This study culminates in a rich dataset that allows us to understand molecular mechanisms, better enabling us to combat insecticide resistance effectively.

The dynamics of deltamethrin resistance evolution in Aedes albopictus has an impact on fitness and dengue virus type-2 vectorial capacity

BACKGROUND

Worldwide invasion and expansion of Aedes albopictus, an important vector of dengue, chikungunya, and Zika viruses, has become a serious concern in global public health. Chemical insecticides are the primary means currently available to control the mosquito populations. However, long-term and large-scale use of insecticides has selected for resistance in the mosquito that is accompanied by a genetic load that impacts fitness.

RESULTS

A number of laboratory strains representing different resistance mechanisms were isolated and identified from laboratory-derived, deltamethrin-resistant Ae. albopictus recovered in previous work. Resistance levels and fitness costs of the strains were evaluated and compared to characterize the evolution of the resistance genotypes and phenotypes. The heterozygous F1534S mutation (1534F/S) in the voltage gated sodium channel (vgsc) gene product (VGSC), first detected in early stages of resistance evolution, not only confers high-level resistance, but also produces no significant fitness costs, leading to the rapid spread of resistance in the population. This is followed by the increase in frequency of homozygous F1534S (1534S/S) mosquitoes that have significant fitness disadvantages, prompting the emergence of an unlinked I1532T mutation with fewer side effects and a mating advantage better adapted to the selection and reproductive pressures imposed in the experiments. Metabolic resistance with no significant fitness cost and mediating a high-tolerance resistance phenotype may play a dominant role in the subsequent evolution of resistance. The different resistant strains had similar vector competence for dengue virus type-2 (DENV-2). Furthermore, a comparative analysis of vectorial capacity revealed that increased survival due to deltamethrin resistance balanced the negative fitness cost effects and contributed to the risk of dengue virus (DENV) transmission by resistant populations. The progressive evolution of resistance results in mosquitoes with both target-site insensitivity and metabolic resistance with lower fitness costs, which further leads to resistant populations with both high resistance levels and vectorial capacity.

CONCLUSIONS

This study reveals a possible mechanism for the evolution of deltamethrin resistance in Aedes albopictus. These findings will help guide practical strategies for insecticide use, resistance management and the prevention and control of mosquito-borne disease.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae and Anopheles coluzzii

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of malaria control, but the genetics of resistance are only partially understood. We performed a large scale multi-country genome-wide association study of resistance to two widely used insecticides: deltamethrin and pirimiphos-methyl, using sequencing data from An. gambiae and An. coluzzii from ten locations in West Africa. Resistance was highly multi-genic, multi-allelic and variable between populations. While the strongest and most consistent association with deltamethrin resistance came from Cyp6aa1, this was based on several independent copy number variants (CNVs) in An. coluzzii, and on a non-CNV haplotype in An. gambiae. For pirimiphos-methyl, signals included Ace1, cytochrome P450s, glutathione S-transferases and the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes showed evidence of cross-resistance to both insecticides. These locally-varying, multi-allelic patterns highlight the challenges involved in genomic monitoring of resistance, and may form the basis for improved surveillance methods.

LLIN Evaluation in Uganda Project (LLINEUP) - Plasmodium infection prevalence and genotypic markers of insecticide resistance in Anopheles vectors from 48 districts of Uganda

Background

In 2017-2019, we conducted a large-scale, cluster-randomised trial (LLINEUP) to evaluate long-lasting insecticidal nets (LLINs) treated with a pyrethroid insecticide plus the synergist piperonyl butoxide (PBO LLINs), as compared to conventional, pyrethroid-only LLINs across 104 health sub-districts (HSDs) in Uganda. In LLINEUP, and similar trials in Tanzania, PBO LLINs were found to provide greater protection against malaria than conventional LLINs, reducing parasitaemia and vector density. In the LLINEUP trial, cross-sectional entomological surveys were carried out at baseline and then every 6 months for two years. In each survey, ten households per HSD were randomly selected for indoor household entomological collections.

Results

Overall, 5395 female Anopheles mosquitoes were collected from 5046 households. The proportion of mosquitoes infected with Plasmodium falciparum did not change significantly over time, while infection with non-falciparum malaria decreased in An. gambiae s.s, but not An. funestus. The frequency of genetic markers associated with pyrethroid resistance increased significantly over time, but the rate of change was not different between the two LLIN types. The knock-down resistance (kdr) mutation Vgsc-995S declined over time as Vgsc-995F, the alternative resistance mutation at this codon, increased. Vgsc-995F appears to be spreading into Uganda.

Conclusions

Distribution of LLINs in Uganda was associated with reductions in parasite prevalence and vector density, but the proportion of infective mosquitoes remained stable, suggesting that the potential for transmission persisted. The increased frequency of markers of pyrethroid resistance indicates that LLIN distribution favoured the evolution of resistance within local vectors and highlights the potential benefits of resistance management strategies.Trial registration:: This study is registered with ISRCTN, ISRCTN17516395. Registered 14 February 2017, http://www.isrctn.com/ISRCTN17516395.

RNA-Seq-Pop: Exploiting the sequence in RNA sequencing-A Snakemake workflow reveals patterns of insecticide resistance in the malaria vector Anopheles gambiae

We provide a reproducible and scalable Snakemake workflow, called RNA-Seq-Pop, which provides end-to-end analysis of RNA sequencing data sets. The workflow allows the user to perform quality control, perform differential expression analyses and call genomic variants. Additional options include the calculation of allele frequencies of variants of interest, summaries of genetic variation and population structure, and genome-wide selection scans, together with clear visualizations. RNA-Seq-Pop is applicable to any organism, and we demonstrate the utility of the workflow by investigating pyrethroid resistance in selected strains of the major malaria mosquito, Anopheles gambiae. The workflow provides additional modules specifically for An. gambiae, including estimating recent ancestry and determining the karyotype of common chromosomal inversions. The Busia laboratory colony used for selections was collected in Busia, Uganda, in November 2018. We performed a comparative analysis of three groups: a parental G24 Busia strain; its deltamethrin-selected G28 offspring; and the susceptible reference strain Kisumu. Measures of genetic diversity reveal patterns consistent with that of laboratory colonization and selection, with the parental Busia strain exhibiting the highest nucleotide diversity, followed by the selected Busia offspring, and finally, Kisumu. Differential expression and variant analyses reveal that the selected Busia colony exhibits a number of distinct mechanisms of pyrethroid resistance, including the Vgsc-995S target-site mutation, upregulation of SAP genes, P450s and a cluster of carboxylesterases. During deltamethrin selections, the 2La chromosomal inversion rose in frequency (from 33% to 86%), supporting a previous link with pyrethroid resistance. RNA-Seq-Pop is hosted at: github.com/sanjaynagi/rna-seq-pop. We anticipate that the workflow will provide a useful tool to facilitate reproducible, transcriptomic studies in An. gambiae and other taxa.

Widespread occurrence of copy number variants and fixation of pyrethroid target site resistance in Anopheles gambiae (s.l.) from southern Côte d'Ivoire

Resistance to pyrethroid and organophosphate insecticides in the malaria vector Anopheles gambiae (s.l.) is conferred by a variety of genetic mutations, including single nucleotide polymorphisms (SNPs) and copy number variants (CNVs). Knowledge of the distribution of these mutations in mosquito populations is a prerequisite for establishing better strategies for their management. In this study, a total of 755 Anopheles gambiae (s.l.) from southern Côte d'Ivoire were exposed to deltamethrin or pirimiphos-methyl insecticides and were screened to assess the distribution of SNPs and CNVs known or believed to confer resistance to one or other of the insecticide classes. Most individuals from the An. gambiae (s.l.) complex were identified by molecular tests as Anopheles coluzzii. Survival to deltamethrin (from 94% to 97%) was higher than to pirimiphos-methyl (from 10% to 49%). In An. gambiae (s.s.), the SNP in the Voltage Gated Sodium Channel (Vgsc) at the 995F locus (Vgsc-995F) was fixed, while other target site mutations were rare or absent (Vgsc-402L: 0%; Vgsc-1570Y: 0%, Acetylcholinesterase Acel-280S: 14%). In An. coluzzii, Vgsc-995F was the target site SNP found at highest frequency (65%) followed by other target site mutations (Vgsc-402L: 36%; Vgsc-1570Y: 0.33%; Acel-280S: 45%). The Vgsc-995S SNP was not present. The presence of the Ace1-280S SNP was found to be significantly linked to the presence of the Ace1-CNV, Ace1_AgDup. Significant association was found between the presence of the Ace1_AgDup and pirimiphos-methyl resistance in An. gambiae (s.s.) but not in An. coluzzii. The deletion Ace1_Del97 was found in one specimen of An. gambiae (s.s.). Four CNVs in the Cyp6aa/Cyp6p gene cluster, which contains genes of known importance for resistance, were detected in An. coluzzii, the most frequent being Dup 7 (42%) and Dup 14 (26%). While none of these individual CNV alleles were significantly associated with resistance, copy number in the Cyp6aa gene region in general was associated with increased resistance to deltamethrin. Elevated expression of Cyp6p3 was nearly associated with deltamethrin resistance, although there was no association of resistance with copy number. Use of alternative insecticides and control methods to arrest resistance spread in An. coluzzii populations is merited.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae s.l

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of the most widespread tools currently used to control malaria. The genetic underpinnings of resistance are still only partially understood, with much of the variance in resistance phenotype left unexplained. We performed a multi-country large scale genome-wide association study of resistance to two insecticides widely used in malaria control: deltamethrin and pirimiphos-methyl. Using a bioassay methodology designed to maximise the phenotypic difference between resistant and susceptible samples, we sequenced 969 phenotyped female An. gambiae and An. coluzzii from ten locations across four countries in West Africa (Benin, Côte d'Ivoire, Ghana and Togo), identifying single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) segregating in the populations. The patterns of resistance association were highly multiallelic and variable between populations, with different genomic regions contributing to resistance, as well as different mutations within a given region. While the strongest and most consistent association with deltamethrin resistance came from the region around Cyp6aa1 , this resistance was based on a combination of several independent CNVs in An. coluzzii , and on a non-CNV bearing haplotype in An. gambiae . Further signals involved a range of cytochrome P450, mitochondrial, and immunity genes. Similarly, for pirimiphos-methyl, while the strongest signal came from the region of Ace1 , more widespread signals included cytochrome P450s, glutathione S-transferases, and a subunit of the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes were associated with resistance to both insecticide classes, suggesting possible cross-resistance mechanisms. These locally-varying, multigenic and multiallelic patterns highlight the challenges involved in genomic monitoring and surveillance of resistance, and form the basis for improvement of methods used to detect and predict resistance. Based on simulations of resistance variants, we recommend that yet larger scale studies, exceeding 500 phenotyped samples per population, are required to better identify associated genomic regions.

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.