Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates.

BACKGROUND: Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl). RESULTS: Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups. CONCLUSION: Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

Whole transcriptomic analysis reveals overexpression of salivary gland and cuticular proteins genes in insecticide-resistant Anopheles arabiensis from Western Kenya.

BACKGROUND: Effective vector control is key to malaria prevention. However, this is now compromised by increased insecticide resistance due to continued reliance on insecticide-based control interventions. In Kenya, we have observed heterogenous resistance to pyrethroids and organophosphates in Anopheles arabiensis which is one of the most widespread malaria vectors in the country. We investigated the gene expression profiles of insecticide resistant An. arabiensis populations from Migori and Siaya counties in Western Kenya using RNA-Sequencing. Centers for Disease Control and Prevention (CDC) bottle assays were conducted using deltamethrin (DELTA), alphacypermethrin (ACYP) and pirimiphos-methyl (PMM) to determine the resistance status in both sites. RESULTS: Mosquitoes from Migori had average mortalities of 91%, 92% and 58% while those from Siaya had 85%, 86%, and 30% when exposed to DELTA, ACYP and PMM, respectively. RNA-Seq analysis was done on pools of mosquitoes which survived exposure ('resistant'), mosquitoes that were not exposed, and the insecticide-susceptible An. arabiensis Dongola strain. Gene expression profiles of resistant mosquitoes from both Migori and Siaya showed an overexpression mainly of salivary gland proteins belonging to both the short and long form D7 genes, and cuticular proteins (including CPR9, CPR10, CPR15, CPR16). Additionally, the overexpression of detoxification genes including cytochrome P450s (CYP9M1, CYP325H1, CYP4C27, CYP9L1 and CYP307A1), 2 carboxylesterases and a glutathione-S-transferase (GSTE4) were also shared between DELTA, ACYP, and PMM survivors, pointing to potential contribution to cross resistance to both pyrethroid and organophosphate insecticides. CONCLUSION: This study provides novel insights into the molecular basis of insecticide resistance in An. arabiensis in Western Kenya and suggests that salivary gland proteins and cuticular proteins are associated with resistance to multiple classes of insecticides.

Key gene modules and hub genes associated with pyrethroid and organophosphate resistance in Anopheles mosquitoes: a systems biology approach.

Indoor residual spraying (IRS) and insecticide-treated nets (ITNs) are the main methods used to control mosquito populations for malaria prevention. The efficacy of these strategies is threatened by the spread of insecticide resistance (IR), limiting the success of malaria control. Studies of the genetic evolution leading to insecticide resistance could enable the identification of molecular markers that can be used for IR surveillance and an improved understanding of the molecular mechanisms associated with IR. This study used a weighted gene co-expression network analysis (WGCNA) algorithm, a systems biology approach, to identify genes with similar co-expression patterns (modules) and hub genes that are potential molecular markers for insecticide resistance surveillance in Kenya and Benin. A total of 20 and 26 gene co-expression modules were identified via average linkage hierarchical clustering from Anopheles arabiensis and An. gambiae, respectively, and hub genes (highly connected genes) were identified within each module. Three specific genes stood out: serine protease, E3 ubiquitin-protein ligase, and cuticular proteins, which were top hub genes in both species and could serve as potential markers and targets for monitoring IR in these malaria vectors. In addition to the identified markers, we explored molecular mechanisms using enrichment maps that revealed a complex process involving multiple steps, from odorant binding and neuronal signaling to cellular responses, immune modulation, cellular metabolism, and gene regulation. Incorporation of these dynamics into the development of new insecticides and the tracking of insecticide resistance could improve the sustainable and cost-effective deployment of interventions.

Entomological indicators of Plasmodium species transmission in Goma Tsé-Tsé and Madibou districts, in the Republic of Congo.

BACKGROUND: Malaria remains a major public health problem in the Republic of Congo, with Plasmodium falciparum being the deadliest species of Plasmodium in humans. Vector transmission of malaria is poorly studied in the country and no previous report compared rural and urban data. This study aimed to determine the Anopheles fauna and the entomological indices of malaria transmission in the rural and urban areas in the south of Brazzaville, and beyond. METHODS: Indoor household mosquitoes capture using electric aspirator was performed in rural and urban areas during raining and dry seasons in 2021. The identification of Anopheles species was done using binocular magnifier and nested-PCR. TaqMan and nested-PCR were used to detect the Plasmodium species in the head/thorax and abdomens of Anopheles. Some entomological indices including the sporozoite infection rate, the entomological inoculation rate and the man biting rate were estimated. RESULTS: A total of 699 Anopheles mosquitoes were collected: Anopheles gambiae sensu lato (s.l.) (90.7%), Anopheles funestus s.l. (6.9%), and Anopheles moucheti (2.4%). Three species of An. gambiae s.l. were identified including Anopheles gambiae sensu stricto (78.9%), Anopheles coluzzii (15.4%) and Anopheles arabiensis (5.7%). The overall sporozoite infection rate was 22.3% with a predominance of Plasmodium falciparum, followed by Plasmodium malariae and Plasmodium ovale. Anopheles aggressiveness rate was higher in households from rural area (1.1 bites/night) compared to that from urban area (0.8 ib/p/n). The overall entomological inoculation rate was 0.13 ib/p/n. This index was 0.17 ib/p/n and 0.092 ib/p/n in rural and in urban area, respectively, and was similar during the dry (0.18 ib/p/n) and rainy (0.14 ib/p/n) seasons. CONCLUSION: These findings highlight that malaria transmission remains high in rural and urban area in the south of Republic of Congo despite the ongoing control efforts, thereby indicating the need for more robust interventions.

Physico-chemical characterization of Anopheles gambiae s.l. breeding sites and kdr mutations in urban areas of Cotonou and Natitingou, Benin.

BACKGROUND: This study aimed to investigate the relationship between the physicochemical characteristics of An. gambiae s.s. and An. coluzzii breeding sites, the susceptibility profiles to commonly used insecticides in public health, and the underlying insecticide resistance mechanisms. METHODS: Anopheles breeding sites surveys were conducted in Cotonou and Natitingou in September 2020, January and August 2021. Physicochemical properties and bacterial loads were determined in individual breeding sites. The WHO susceptibility assays were carried out using the female of the emerging adult mosquitoes. Anopheles species were identified through PCR techniques. Kdr L1014F/S, N1575Y and G119S mutations were investigated using TaqMan genotyping assays. RESULTS: Molecular analysis showed that all mosquitoes analyzed in Cotonou were Anopheles coluzzii, while those of Natitingou were Anopheles gambiae s.s. Fecal coliforms were identified as playing a role in this distribution through their significant influence on the presence of An. coluzzii larvae. WHO susceptibility assay indicated a high level of resistance to deltamethrin in the two cities. The resistance levels to deltamethrin were higher in Cotonou (X2 = 31.689; DF = 1; P < 0.0001). There was a suspected resistance to bendiocarb in Cotonou, whereas the mosquito population in Natitingou was resistant. The kdr L1014F mutation was highly observed in both mosquito populations (frequence: 86-91%), while the Ace-1 mutation was found in a small proportion of mosquitoes. In Cotonou, salinity was the only recorded physicochemical parameter that significantly correlated with the resistance of Anopheles mosquitoes to deltamethrin (P < 0.05). In Natitingou, significant correlations were observed between the allelic frequencies of the kdr L1014F mutation and pH, conductivity, and TDS. CONCLUSION: These results indicate a high level of pyrethroid resistance in the anopheles populations of both Cotonou and Natitingou. Moreover, this study report the involvement of abiotic factors influencing Anopheles susceptibility profile.

Resistance to pirimiphos-methyl in West African Anopheles is spreading via duplication and introgression of the Ace1 locus.

Vector population control using insecticides is a key element of current strategies to prevent malaria transmission in Africa. The introduction of effective insecticides, such as the organophosphate pirimiphos-methyl, is essential to overcome the recurrent emergence of resistance driven by the highly diverse Anopheles genomes. Here, we use a population genomic approach to investigate the basis of pirimiphos-methyl resistance in the major malaria vectors Anopheles gambiae and A. coluzzii. A combination of copy number variation and a single non-synonymous substitution in the acetylcholinesterase gene, Ace1, provides the key resistance diagnostic in an A. coluzzii population from Côte d'Ivoire that we used for sequence-based association mapping, with replication in other West African populations. The Ace1 substitution and duplications occur on a unique resistance haplotype that evolved in A. gambiae and introgressed into A. coluzzii, and is now common in West Africa primarily due to selection imposed by other organophosphate or carbamate insecticides. Our findings highlight the predictive value of this complex resistance haplotype for phenotypic resistance and clarify its evolutionary history, providing tools to for molecular surveillance of the current and future effectiveness of pirimiphos-methyl based interventions.

Vectorial competence, insecticide resistance in Anopheles funestus and operational implications for malaria vector control strategies in Benin Republic.

The primary reason for the failure of malaria vector control across endemic regions is the widespread insecticide resistance observed in Anopheles vectors. The most dominant African vectors of malaria parasites are Anopheles gambiae and Anopheles funestus mosquitoes. These species often exhibit divergent behaviours and adaptive changes underscoring the importance of deploying active and effective measures in their control. Unlike An. gambiae, An. funestus mosquitoes are poorly studied in Benin Republic. However, recent reports indicated that An. funestus can adapt and colonize various ecological niches owing to its resistance against insecticides and adaptation to changing breeding habitats. Unfortunately, scientific investigations on the contribution of An. funestus to malaria transmission, their susceptibility to insecticide and resistance mechanism developed are currently insufficient for the design of better control strategies. In an attempt to gather valuable information on An. funestus, the present review examines the progress made on this malaria vector species in Benin Republic and highlights future research perspectives on insecticide resistance profiles and related mechanisms, as well as new potential control strategies against An. funestus. Literature analysis revealed that An. funestus is distributed all over the country, although present in low density compared to other dominant malaria vectors. Interestingly, An. funestus is being found in abundance during the dry seasons, suggesting an adaptation to desiccation. Among the An. funestus group, only An. funestus sensu stricto (s.s.) and Anopheles leesoni were found in the country with An. funestus s.s. being the most abundant species. Furthermore, An. funestus s.s. is the only one species in the group contributing to malaria transmission and have adapted biting times that allow them to bite at dawn. In addition, across the country, An. funestus were found resistant to pyrethroid insecticides used for bed nets impregnation and also resistant to bendiocarb which is currently being introduced in indoor residual spraying formulation in malaria endemic regions. All these findings highlight the challenges faced in controlling this malaria vector. Therefore, advancing the knowledge of vectorial competence of An. funestus, understanding the dynamics of insecticide resistance in this malaria vector, and exploring alternative vector control measures, are critical for sustainable malaria control efforts in Benin Republic.

Relative expression of key genes involved in nucleic acids methylation in Anopheles gambiae sensu stricto.

In vertebrates, enzymes responsible for DNA methylation, one of the epigenetic mechanisms, are encoded by genes falling into the cytosine methyltransferases genes family (Dnmt1, Dnmt3a,b and Dnmt3L). However, in Diptera, only the methyltransferase Dnmt2 was found, suggesting that DNA methylation might act differently for species in this order. Moreover, genes involved in epigenetic dynamics, such as Ten-eleven Translocation dioxygenases (TET) and Methyl-CpG-binding domain (MBDs), present in vertebrates, might play a role in insects. This work aimed at investigating nucleic acids methylation in the malaria vector Anopheles gambiae (Diptera: Culicidae) by analysing the expression of Dnmt2, TET2 and MBDs genes using quantitative real-time polymerase chain reaction (qRT-PCR) at pre-immature stages and in reproductive tissues of adult mosquitoes. In addition, the effect of two DNA methylation inhibitors on larval survival was evaluated. The qPCR results showed an overall low expression of Dnmt2 at all developmental stages and in adult reproductive tissues. In contrast, MBD and TET2 showed an overall higher expression. In adult mosquito reproductive tissues, the expression level of the three genes in males' testes was significantly higher than that in females' ovaries. The chemical treatments did not affect larval survival. The findings suggest that mechanisms other than DNA methylation underlie epigenetic regulation in An. gambiae.

Plasmodium falciparum msp1 and msp2 genetic diversity in parasites isolated from symptomatic and asymptomatic malaria subjects in the South of Benin.

Symptomatic and asymptomatic malaria patients are considered as the reservoirs of human Plasmodium. In the present study, we have evaluated the Plasmodium falciparum merozoite surface protein-1 (Pfmsp1) and protein-2 (Pfmsp2) genetic diversity among the symptomatic and asymptomatic malaria infection from health facilities in Cotonou, Benin Republic. A cross-sectional study recruited 158 individuals, including 77 from the asymptomatic and 81 from the symptomatic groups. The parasites were genotyped using Nested Polymerase Chain Reaction. Samples identified as Plasmodium falciparum were genotyped for their genetic diversity. No significant difference was observed in the overall multiplicity of infection (MOI) between the asymptomatic and symptomatic groups. In the symptomatic group, the overall frequency of K1, MAD20, and RO33 allelic family was more predominant (98.5%) followed by 3D7 (87.3%) and FC27 (83.1%). However, in asymptomatic group, the K1 alleles were the most prevalent (100%) followed by FC27 (89.9%), 3D7 (76.8%), MAD20 (60.5%), and RO33 (35.5%). The frequency of multiple allelic types (K1+MAD20+RO33) at the Pfmsp1 loci in the symptomatic infections was significantly higher when compared to that of the asymptomatic ones (97% vs. 34%, p < 0.05), whereas no difference was observed in the frequency of multiple allelic types (3D7 and FC27) at the Pfmsp2 loci between the two groups. The high presence of msp1 multiple infections in the symptomatic group compared to asymptomatic ones suggests an association between the genetic diversity and the onset of malaria symptoms. These data can provide valuable information in the development of a vaccine that could reduce the symptomatic disease.

Surveillance of Plasmodium malariae infection among inhabitants of rural areas in Ouidah-Kpomasse-Tori Bossito health district, Benin.

Among the Plasmodium species that infect humans, P. falciparum has been largely studied in malaria endemic areas. However, P. malariae infection is less documented among the human population. This study aimed to monitor the prevalence and distribution of P. malariae in Southern Benin. A cross-sectional survey was conducted in rural localities in the Ouidah-Kpomasse-Tori Bossito (OKT) health district in Southern Benin from June to October 2019. Socio-demographic data were collected using a questionnaire, while malaria infection data were obtained on the one hand by microscopy diagnosis and, on the other, by nested polymerase chain reaction (PCR). Based on microscopy, the prevalence of P. malariae mono-infection and coinfection of P. falciparum, P. malariae was respectively 2.3% and 1.2% in the OKT health district. This prevalence was higher (P < 0.01) than that reported by Damien et al. (2010) 10 years ago in the same study area with 0.7% and 0.3% of P. malariae and P. falciparum/P. malariae, respectively. Based on PCR analysis, P. malariae prevalence was 14.1%, including 5.2% of mono-infection and 8.9% of mixed infection with P. falciparum. Sub-microscopic Plasmodium infections were high (30.6%) and more pronounced in older participants (>20 years). The present study revealed that P. malariae increased in the OKT health district with a high prevalence of submicroscopic infection. Since our results provide valuable evidence of increasing P. malariae infection, the National Malaria Control Programs (NMCPs) must consider P. malariae when designing future measures for effective control and malaria treatment.

Experimental Transmission of Plasmodium malariae to Anopheles gambiae.

Our current knowledge of the clinical burden, biology, and transmission of Plasmodium malariae is extremely scarce. To start addressing some of those questions, we experimentally infected Anopheles gambiae mosquitoes with fresh P. malariae isolates obtained from asymptomatic individuals in Lambaréné, Gabon. The proportion of mosquitoes infected via direct membrane feeding assay with either P. malariae monoinfections (16% [19 of 121]) or coinfections (28% [31 of 112]) was higher after serum replacement than in parallel groups without serum replacement (4% [4 of 102] and 4% [2 of 45], respectively; P < .01). Our results show that isolates from asymptomatic carriers can be used for experimental studies of P. malariae transmission.

Putative pleiotropic effects of the knockdown resistance (L1014F) allele on the life-history traits of Anopheles gambiae.

BACKGROUND: Existing mechanisms of insecticide resistance are known to help the survival of mosquitoes following contact with chemical compounds, even though they could negatively affect the life-history traits of resistant malaria vectors. In West Africa, the knockdown resistance mechanism kdrR (L1014F) is the most common. However, little knowledge is available on its effects on mosquito life-history traits. The fitness effects associated with this knockdown resistance allele in Anopheles gambiae sensu stricto (s.s.) were investigated in an insecticide-free laboratory environment. METHODS: The life-history traits of Kisumu (susceptible) and KisKdr (kdr resistant) strains of An. gambiae s.s. were compared. Larval survivorship and pupation rate were assessed as well as fecundity and fertility of adult females. Female mosquitoes of both strains were directly blood fed through artificial membrane assays and then the blood-feeding success, blood volume and adult survivorship post-blood meal were assessed. RESULTS: The An. gambiae mosquitoes carrying the kdrR allele (KisKdr) laid a reduced number of eggs. The mean number of larvae in the susceptible strain Kisumu was three-fold overall higher than that seen in the KisKdr strain with a significant difference in hatching rates (81.89% in Kisumu vs 72.89% in KisKdr). The KisKdr larvae had a significant higher survivorship than that of Kisumu. The blood-feeding success was significantly higher in the resistant mosquitoes (84%) compared to the susceptible ones (34.75%). However, the mean blood volume was 1.36 µL/mg, 1.45 µL/mg and 1.68 µL/mg in Kisumu, homozygote and heterozygote KisKdr mosquitoes, respectively. After blood-feeding, the heterozygote KisKdr mosquitoes displayed highest survivorship when compared to that of Kisumu. CONCLUSIONS: The presence of the knockdown resistance allele appears to impact the life-history traits, such as fecundity, fertility, larval survivorship, and blood-feeding behaviour in An. gambiae. These data could help to guide the implementation of more reliable strategies for the control of malaria vectors.

Interplay Between Oxytetracycline and the Homozygote kdr (L1014F) Resistance Genotype on Fecundity in Anopheles gambiae (Diptera: Culicidae) Mosquitoes.

The insecticide resistance in Anopheles gambiae mosquitoes has remained the major threat for vector control programs but the fitness effects conferred by these mechanisms are poorly understood. To fill this knowledge gap, the present study aimed at testing the hypothesis that antibiotic oxytetracycline could have an interaction with insecticide resistance genotypes and consequently inhibit the fecundity in An. gambiae. Four strains of An. gambiae: Kisumu (susceptible), KisKdr (kdr (L1014F) resistant), AcerKis (ace-1 (G119S) resistant) and AcerKdrKis (both kdr (L1014F) and ace-1 (G119S) resistant) were used in this study. The different strains were allowed to bloodfeed on a rabbit previously treated with antibiotic oxytetracycline at a concentration of 39·10-5 M. Three days later, ovarian follicles were dissected from individual mosquito ovaries into physiological saline solution (0.9% NaCl) under a stereomicroscope and the eggs were counted. Fecundity was substantially lower in oxytetracycline-exposed KisKdr females when compared to that of the untreated individuals and oxytetracycline-exposed Kisumu females. The exposed AcerKis females displayed an increased fecundity compared to their nontreated counterparts whereas they had reduced fecundity compared to that of oxytetracycline-exposed Kisumu females. There was no substantial difference between the fecundity in the treated and untreated AcerKdrKis females. The oxytetracycline-exposed AcerKdrKis mosquitoes had an increased fecundity compared to that of the exposed Kisumu females. Our data indicate an indirect effect of oxytetracycline in reducing fecundity of An. gambiae mosquitoes carrying kdrR (L1014F) genotype. These findings could be useful for designing new integrated approaches for malaria vector control in endemic countries.

The ace-1 Locus Is Amplified in All Resistant Anopheles gambiae Mosquitoes: Fitness Consequences of Homogeneous and Heterogeneous Duplications.

Gene copy-number variations are widespread in natural populations, but investigating their phenotypic consequences requires contemporary duplications under selection. Such duplications have been found at the ace-1 locus (encoding the organophosphate and carbamate insecticides' target) in the mosquito Anopheles gambiae (the major malaria vector); recent studies have revealed their intriguing complexity, consistent with the involvement of various numbers and types (susceptible or resistant to insecticide) of copies. We used an integrative approach, from genome to phenotype level, to investigate the influence of duplication architecture and gene-dosage on mosquito fitness. We found that both heterogeneous (i.e., one susceptible and one resistant ace-1 copy) and homogeneous (i.e., identical resistant copies) duplications segregated in field populations. The number of copies in homogeneous duplications was variable and positively correlated with acetylcholinesterase activity and resistance level. Determining the genomic structure of the duplicated region revealed that, in both types of duplication, ace-1 and 11 other genes formed tandem 203kb amplicons. We developed a diagnostic test for duplications, which showed that ace-1 was amplified in all 173 resistant mosquitoes analyzed (field-collected in several African countries), in heterogeneous or homogeneous duplications. Each type was associated with different fitness trade-offs: heterogeneous duplications conferred an intermediate phenotype (lower resistance and fitness costs), whereas homogeneous duplications tended to increase both resistance and fitness cost, in a complex manner. The type of duplication selected seemed thus to depend on the intensity and distribution of selection pressures. This versatility of trade-offs available through gene duplication highlights the importance of large mutation events in adaptation to environmental variation. This impressive adaptability could have a major impact on vector control in Africa.

CYP6 P450 enzymes and ACE-1 duplication produce extreme and multiple insecticide resistance in the malaria mosquito Anopheles gambiae.

Malaria control relies heavily on pyrethroid insecticides, to which susceptibility is declining in Anopheles mosquitoes. To combat pyrethroid resistance, application of alternative insecticides is advocated for indoor residual spraying (IRS), and carbamates are increasingly important. Emergence of a very strong carbamate resistance phenotype in Anopheles gambiae from Tiassalé, Côte d'Ivoire, West Africa, is therefore a potentially major operational challenge, particularly because these malaria vectors now exhibit resistance to multiple insecticide classes. We investigated the genetic basis of resistance to the most commonly-applied carbamate, bendiocarb, in An. gambiae from Tiassalé. Geographically-replicated whole genome microarray experiments identified elevated P450 enzyme expression as associated with bendiocarb resistance, most notably genes from the CYP6 subfamily. P450s were further implicated in resistance phenotypes by induction of significantly elevated mortality to bendiocarb by the synergist piperonyl butoxide (PBO), which also enhanced the action of pyrethroids and an organophosphate. CYP6P3 and especially CYP6M2 produced bendiocarb resistance via transgenic expression in Drosophila in addition to pyrethroid resistance for both genes, and DDT resistance for CYP6M2 expression. CYP6M2 can thus cause resistance to three distinct classes of insecticide although the biochemical mechanism for carbamates is unclear because, in contrast to CYP6P3, recombinant CYP6M2 did not metabolise bendiocarb in vitro. Strongly bendiocarb resistant mosquitoes also displayed elevated expression of the acetylcholinesterase ACE-1 gene, arising at least in part from gene duplication, which confers a survival advantage to carriers of additional copies of resistant ACE-1 G119S alleles. Our results are alarming for vector-based malaria control. Extreme carbamate resistance in Tiassalé An. gambiae results from coupling of over-expressed target site allelic variants with heightened CYP6 P450 expression, which also provides resistance across contrasting insecticides. Mosquito populations displaying such a diverse basis of extreme and cross-resistance are likely to be unresponsive to standard insecticide resistance management practices.

Contemporary evolution of resistance at the major insecticide target site gene Ace-1 by mutation and copy number variation in the malaria mosquito Anopheles gambiae.

Functionally constrained genes are ideal insecticide targets because disruption is often fatal, and resistance mutations are typically costly. Synaptic acetylcholinesterase (AChE) is an essential neurotransmission enzyme targeted by insecticides used increasingly in malaria control. In Anopheles and Culex mosquitoes, a glycine-serine substitution at codon 119 of the Ace-1 gene confers both resistance and fitness costs, especially for 119S/S homozygotes. G119S in Anopheles gambiae from Accra (Ghana) is strongly associated with resistance, and, despite expectations of cost, resistant 119S alleles are increasing significantly in frequency. Sequencing of Accra females detected only a single Ace-1 119S haplotype, whereas 119G diversity was high overall but very low at non-synonymous sites, evidence of strong purifying selection driven by functional constraint. Flanking microsatellites showed reduced diversity, elevated linkage disequilibrium and high differentiation of 119S, relative to 119G homozygotes across up to two megabases of the genome. Yet these signals of selection were inconsistent and sometimes weak tens of kilobases from Ace-1. This unexpected finding is attributable to apparently ubiquitous amplification of 119S alleles as part of a large copy number variant (CNV) far exceeding the size of the Ace-1 gene, whereas 119G alleles were unduplicated. Ace-1 CNV was detectable in archived samples collected when the 119S allele was rare in Ghana. Multicopy amplification of resistant alleles has not been observed previously and is likely to underpin the recent increase in 119S frequency. The large CNV compromised localization of the strong selective sweep around Ace-1, emphasizing the need to integrate CNV analysis into genome scans for selection.

Estimation of allele-specific Ace-1 duplication in insecticide-resistant Anopheles mosquitoes from West Africa.

BACKGROUND: Identification of variation in Ace-1 copy number and G119S mutation genotype from samples of Anopheles gambiae and Anopheles coluzzii across West Africa are important diagnostics of carbamate and organophosphate resistance at population and individual levels. The most widespread and economical method, PCR-RFLP, suffers from an inability to discriminate true heterozygotes from heterozygotes with duplication. METHODS: In addition to PCR-RFLP, in this study three different molecular techniques were applied on the same mosquito specimens: TaqMan qPCR, qRTPCR and ddPCR. To group heterozygous individuals recorded from the PCR-RFLP analysis into different assumptive genotypes K-means clustering was applied on the Z-scores of data obtained from both the TaqMan and ddPCR methods. The qRTPCR analysis was used for absolute quantification of copy number variation. RESULTS: The results indicate that most heterozygotes are duplicated and that G119S mutation must now be regarded as a complex genotype ranging from primarily single-copy susceptible Glycine homozygotes to balanced and imbalanced heterozygotes, and multiply-amplified resistant Serine allele homozygotes. Whilst qRTPCR-based gene copy analysis suffers from some imprecision, it clearly illustrates differences in copy number among genotype groups identified by TaqMan or ddPCR. Based on TaqMan method properties, and by coupling TaqMan and ddPCR methods simultaneously on the same type of mosquito specimens, it demonstrated that the TaqMan genotype assays associated with the K-means clustering algorithm could provide a useful semi-quantitative estimate method to investigate the level of allele-specific duplication in mosquito populations. CONCLUSIONS: Ace-1 gene duplication is evidently far more complex in An. gambiae and An. coluzzii than the better-studied mosquito Culex quinquefasciatus, which consequently can no longer be considered an appropriate model for prediction of phenotypic consequences. These require urgent further evaluation in Anopheles. To maintain the sustained effectiveness carbamates and organophosphates as alternative products to pyrethroids for malaria vector control, monitoring of duplicated resistant alleles in natural populations is essential to guide the rational use of these insecticides.

Interplay between Plasmodium infection and resistance to insecticides in vector mosquitoes.

Despite its epidemiological importance, the impact of insecticide resistance on vector-parasite interactions and malaria transmission is poorly understood. Here, we explored the impact of Plasmodium infection on the level of insecticide resistance to dichlorodiphenyltrichloroethane (DDT) in field-caught Anopheles gambiae sensu stricto homozygous for the kdr mutation. Results showed that kdr homozygous mosquitoes that fed on infectious blood were more susceptible to DDT than mosquitoes that fed on noninfectious blood during both ookinete development (day 1 after the blood meal) and oocyst maturation (day 7 after the blood meal) but not during sporozoite invasion of the salivary glands. Plasmodium falciparum infection seemed to impose a fitness cost on mosquitoes by reducing the ability of kdr homozygous A. gambiae sensu stricto to survive exposure to DDT. These results suggest an interaction between Plasmodium infection and the insecticide susceptibility of mosquitoes carrying insecticide-resistant alleles. We discuss this finding in relation to vector control efficacy.

Insecticide exposure impacts vector-parasite interactions in insecticide-resistant malaria vectors.

Currently, there is a strong trend towards increasing insecticide-based vector control coverage in malaria endemic countries. The ecological consequence of insecticide applications has been mainly studied regarding the selection of resistance mechanisms; however, little is known about their impact on vector competence in mosquitoes responsible for malaria transmission. As they have limited toxicity to mosquitoes owing to the selection of resistance mechanisms, insecticides may also interact with pathogens developing in mosquitoes. In this study, we explored the impact of insecticide exposure on Plasmodium falciparum development in insecticide-resistant colonies of Anopheles gambiae s.s., homozygous for the ace-1 G119S mutation (Acerkis) or the kdr L1014F mutation (Kdrkis). Exposure to bendiocarb insecticide reduced the prevalence and intensity of P. falciparum oocysts developing in the infected midgut of the Acerkis strain, whereas exposure to dichlorodiphenyltrichloroethane reduced only the prevalence of P. falciparum infection in the Kdrkis strain. Thus, insecticide resistance leads to a selective pressure of insecticides on Plasmodium parasites, providing, to our knowledge, the first evidence of genotype by environment interactions on vector competence in a natural Anopheles-Plasmodium combination. Insecticide applications would affect the transmission of malaria in spite of resistance and would reduce to some degree the impact of insecticide resistance on malaria control interventions.

Contribution of Anopheles gambiae sensu lato mosquitoes to malaria transmission during the dry season in Djoumouna and Ntoula villages in the Republic of the Congo.

BACKGROUND: Mosquitoes belonging to the Anopheles gambiae sensu lato complex play a major role in malaria transmission across Africa. This study assessed the relative importance of members of An. gambiae s.l. in malaria transmission in two rural villages in the Republic of the Congo. METHODS: Adult mosquitoes were collected using electric aspirators from June to September 2022 in Djoumouna and Ntoula villages and were sorted by taxa based on their morphological features. Anopheles gambiae s.l. females were also molecularly identified. A TaqMan-based assay and a nested polymerase chain reaction (PCR) were performed to determine Plasmodium spp. in the mosquitoes. Entomological indexes were estimated, including man-biting rate, entomological inoculation rate (EIR), and diversity index. RESULTS: Among 176 mosquitoes collected, An. gambiae s.l. was predominant (85.8%), followed by Culex spp. (13.6%) and Aedes spp. (0.6%). Three members of the An. gambiae s.l. complex were collected in both villages, namely An. gambiae sensu stricto (74.3%), Anopheles coluzzii (22.9%) and Anopheles arabiensis (2.8%). Three Plasmodium species were detected in An. gambiae s.s. and An. coluzzii (Plasmodium falciparum, P. malariae and P. ovale), while only P. falciparum and P. malariae were found in An. arabiensis. In general, the Plasmodium infection rate was 35.1% (53/151) using the TaqMan-based assay, and nested PCR confirmed 77.4% (41/53) of those infections. The nightly EIR of An. gambiae s.l. was 0.125 infectious bites per person per night (ib/p/n) in Djoumouna and 0.08 ib/p/n in Ntoula. The EIR of An. gambiae s.s. in Djoumouna (0.11 ib/p/n) and Ntoula (0.04 ib/p/n) was higher than that of An. coluzzii (0.01 and 0.03 ib/p/n) and An. arabiensis (0.005 and 0.0 ib/p/n). CONCLUSIONS: This study provides baseline information on the dominant vectors and dynamics of malaria transmission in the rural areas of the Republic of the Congo during the dry season. In the two sampled villages, An. gambiae s.s. appears to play a predominant role in Plasmodium spp.