Joint spatial modelling of malaria incidence and vector's abundance shows heterogeneity in malaria‐vector geographical relationships

Limited attention from the modelling community has been given to ecological approaches which aim to predict geographical patterns of malaria by accounting for the joint effects of different vectors and environmental drivers. A hierarchical multivariate joint spatial Gaussian generalised linear model was developed to provide joint parameters inference and mapping of counts of Anopheles gambiae, An. funestus, An. nili and malaria incidence collected in an area of Cote d'Ivoire. Variable‐selection methods were applied to select important predictors for each mosquito species and malaria incidence. The proposed joint model led to a general reduction of the variance in the estimates compared to independent modelling. There was high variability in the composition of Anopheles mosquito species in the villages with each species suitability only partly overlapping geographically. Abundances of An. gambiae, An. funestus and An. nili were primarily determined by temperature. None of the species were found as a significant predictor for the others. Anopheles gambiae was the predominant species and only An. gambiae female abundance was an important variable (linear predictor) for malaria incidence. However, the geographic correlation analyses show that the rest of Anopheles species are likely playing a role in malaria suitability. Residuals from the models of mosquito abundance and malaria cases are also correlated with each other and overlapping but in geographic patches, meaning that local drivers of vector‐malaria suitability are still present and not represented by the predictors used in the model. Synthesis and applications: Joint modelling improve predictive estimation compared to individual modelling. The accurate predictions highlighted high diversity in the association between malaria and vector species, with most of the area having more than one species suitability correlated with malaria suitability. These zones are unlikely to benefit from species‐specific interventions. Areas with correlated malaria and vector species suitability residuals contain local information, not included in the model, that requires further investigation. This will identify additional communal malaria and vectors factors that need to be considered for optimal malaria control and elimination strategies since these factors are expected to be linked to the local malaria transmission

Small-scale field evaluation of PermaNet ® Dual (a long-lasting net coated with a mixture of chlorfenapyr and deltamethrin) against pyrethroid-resistant Anopheles gambiae mosquitoes from Tiassalé, Côte d’Ivoire

Background: Due to the rapid expansion of pyrethroid-resistance in malaria vectors in Africa, Global Plan for Insecticide Resistance Management (GPIRM) has recommended the development of long-lasting insecticidal nets (LLINs), containing insecticide mixtures of active ingredients with different modes of action to mitigate resistance and improve LLIN efficacy. This good laboratory practice (GLP) study evaluated the efficacy of the chlorfenapyr and deltamethrin-coated PermaNet® Dual, in comparison with the deltamethrin and synergist piperonyl butoxide (PBO)-treated PermaNet® 3.0 and the deltamethrin-coated PermaNet® 2.0, against wild free-flying pyrethroid-resistant Anopheles gambiae sensu lato (s.l.), in experimental huts in Tiassalé, Côte d’Ivoire (West Africa). Methods: PermaNet® Dual, PermaNet® 3.0 and PermaNet® 2.0, unwashed and washed (20 washes), were tested against free-flying pyrethroid-resistant An. gambiae s.l. in the experimental huts in Tiassalé, Côte d’Ivoire from March to August 2020. Complementary laboratory cone bioassays (daytime and 3-min exposure) and tunnel tests (nightly and 15-h exposure) were performed against pyrethroid-susceptible An. gambiae sensu stricto (s.s.) (Kisumu strain) and pyrethroid-resistant An. gambiae s.l. (Tiassalé strain). Results: PermaNet® Dual demonstrated significantly improved efficacy, compared to PermaNet® 3.0 and PermaNet® 2.0, against the pyrethroid-resistant An. gambiae s.l. Indeed, the experimental hut trial data showed that the mortality and blood-feeding inhibition in the wild pyrethroid-resistant An. gambiae s.l. were overall significantly higher with PermaNet® Dual compared with PermaNet® 3.0 and PermaNet® 2.0, for both unwashed and washed samples. The mortality with unwashed and washed samples were 93.6 ± 0.2% and 83.2 ± 0.9% for PermaNet® Dual, 37.5 ± 2.9% and 14.4 ± 3.9% for PermaNet® 3.0, and 7.4 ± 5.1% and 11.7 ± 3.4% for PermaNet® 2.0, respectively. Moreover, unwashed and washed samples produced the respective percentage blood-feeding inhibition of 41.4 ± 6.9% and 43.7 ± 4.8% with PermaNet® Dual, 51.0 ± 5.7% and 9.8 ± 3.6% with PermaNet® 3.0, and 12.8 ± 4.3% and − 13.0 ± 3.6% with PermaNet® 2.0. Overall, PermaNet® Dual also induced higher or similar deterrence, exophily and personal protection when compared with the standard PermaNet® 3.0 and PermaNet® 2.0 reference nets, with both unwashed and washed net samples. In contrast to cone bioassays, tunnel tests predicted the efficacy of PermaNet® Dual seen in the current experimental hut trial. Conclusion: The deltamethrin-chlorfenapyr-coated PermaNet® Dual induced a high efficacy and performed better than the deltamethrin-PBO PermaNet® 3.0 and the deltamethrin-only PermaNet® 2.0, testing both unwashed and 20 times washed samples against the pyrethroid-susceptible and resistant strains of An. gambiae s.l. The inclusion of chlorfenapyr with deltamethrin in PermaNet® Dual net greatly improved protection and control of pyrethroid-resistant An. gambiae populations. PermaNet® Dual thus represents a promising tool, with a high potential to reduce malaria transmission and provide community protection in areas compromised by mosquito vector resistance to pyrethroids

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

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 programmes. The organophosphate, pirimiphos-methyl, is a relatively new chemical in the vector control armoury but is now widely used in indoor residual spray campaigns. Whilst 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 Anopheles 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 organophosphates 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 An. gambiae, An. coluzzii and An. arabiensis. As in Cx. 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 Anopheles 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 programmes

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

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, multiallelic 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 form the basis for improved surveillance methods

Joint spatial modelling of malaria incidence and vector's abundance shows heterogeneity in malaria-vector geographical relationships

Limited attention from the modelling community has been given to ecological approaches which aim to predict geographical patterns of malaria by accounting for the joint effects of different vectors and environmental drivers. A hierarchical multivariate joint spatial Gaussian generalised linear model was developed to provide joint parameters inference and mapping of counts of Anopheles gambiae, An. funestus, An. nili and malaria incidence collected in an area of Cote d'Ivoire. Variable‐selection methods were applied to select important predictors for each mosquito species and malaria incidence. The proposed joint model led to a general reduction of the variance in the estimates compared to independent modelling. There was high variability in the composition of Anopheles mosquito species in the villages with each species suitability only partly overlapping geographically. Abundances of An. gambiae, An. funestus and An. nili were primarily determined by temperature. None of the species were found as a significant predictor for the others. Anopheles gambiae was the predominant species and only An. gambiae female abundance was an important variable (linear predictor) for malaria incidence. However, the geographic correlation analyses show that the rest of Anopheles species are likely playing a role in malaria suitability. Residuals from the models of mosquito abundance and malaria cases are also correlated with each other and overlapping but in geographic patches, meaning that local drivers of vector‐malaria suitability are still present and not represented by the predictors used in the model. Synthesis and applications: Joint modelling improve predictive estimation compared to individual modelling. The accurate predictions highlighted high diversity in the association between malaria and vector species, with most of the area having more than one species suitability correlated with malaria suitability. These zones are unlikely to benefit from species‐specific interventions. Areas with correlated malaria and vector species suitability residuals contain local information, not included in the model, that requires further investigation. This will identify additional communal malaria and vectors factors that need to be considered for optimal malaria control and elimination strategies since these factors are expected to be linked to the local malaria transmission

Small-scale field evaluation of PermaNet® Dual (a long-lasting net coated with a mixture of chlorfenapyr and deltamethrin) against pyrethroid-resistant Anopheles gambiae mosquitoes from Tiassalé, Côte d’Ivoire

Abstract Background Due to the rapid expansion of pyrethroid-resistance in malaria vectors in Africa, Global Plan for Insecticide Resistance Management (GPIRM) has recommended the development of long-lasting insecticidal nets (LLINs), containing insecticide mixtures of active ingredients with different modes of action to mitigate resistance and improve LLIN efficacy. This good laboratory practice (GLP) study evaluated the efficacy of the chlorfenapyr and deltamethrin-coated PermaNet® Dual, in comparison with the deltamethrin and synergist piperonyl butoxide (PBO)-treated PermaNet® 3.0 and the deltamethrin-coated PermaNet® 2.0, against wild free-flying pyrethroid-resistant Anopheles gambiae sensu lato (s.l.), in experimental huts in Tiassalé, Côte d’Ivoire (West Africa). Methods PermaNet® Dual, PermaNet® 3.0 and PermaNet® 2.0, unwashed and washed (20 washes), were tested against free-flying pyrethroid-resistant An. gambiae s.l. in the experimental huts in Tiassalé, Côte d’Ivoire from March to August 2020. Complementary laboratory cone bioassays (daytime and 3-min exposure) and tunnel tests (nightly and 15-h exposure) were performed against pyrethroid-susceptible An. gambiae sensu stricto (s.s.) (Kisumu strain) and pyrethroid-resistant An. gambiae s.l. (Tiassalé strain). Results PermaNet® Dual demonstrated significantly improved efficacy, compared to PermaNet® 3.0 and PermaNet® 2.0, against the pyrethroid-resistant An. gambiae s.l. Indeed, the experimental hut trial data showed that the mortality and blood-feeding inhibition in the wild pyrethroid-resistant An. gambiae s.l. were overall significantly higher with PermaNet® Dual compared with PermaNet® 3.0 and PermaNet® 2.0, for both unwashed and washed samples. The mortality with unwashed and washed samples were 93.6 ± 0.2% and 83.2 ± 0.9% for PermaNet® Dual, 37.5 ± 2.9% and 14.4 ± 3.9% for PermaNet® 3.0, and 7.4 ± 5.1% and 11.7 ± 3.4% for PermaNet® 2.0, respectively. Moreover, unwashed and washed samples produced the respective percentage blood-feeding inhibition of 41.4 ± 6.9% and 43.7 ± 4.8% with PermaNet® Dual, 51.0 ± 5.7% and 9.8 ± 3.6% with PermaNet® 3.0, and 12.8 ± 4.3% and − 13.0 ± 3.6% with PermaNet® 2.0. Overall, PermaNet® Dual also induced higher or similar deterrence, exophily and personal protection when compared with the standard PermaNet® 3.0 and PermaNet® 2.0 reference nets, with both unwashed and washed net samples. In contrast to cone bioassays, tunnel tests predicted the efficacy of PermaNet® Dual seen in the current experimental hut trial. Conclusion The deltamethrin-chlorfenapyr-coated PermaNet® Dual induced a high efficacy and performed better than the deltamethrin-PBO PermaNet® 3.0 and the deltamethrin-only PermaNet® 2.0, testing both unwashed and 20 times washed samples against the pyrethroid-susceptible and resistant strains of An. gambiae s.l. The inclusion of chlorfenapyr with deltamethrin in PermaNet® Dual net greatly improved protection and control of pyrethroid-resistant An. gambiae populations. PermaNet® Dual thus represents a promising tool, with a high potential to reduce malaria transmission and provide community protection in areas compromised by mosquito vector resistance to pyrethroids

Genes significantly overexpressed (relative to susceptible samples) in (A) Tiassalé bendiocarb resistant samples in Exp1, and (B). both Tiassalé and Kovié samples.

<p>Plots show: (A). Log₂-transformed fold-changes (FC) plotted against -log₁₀ transformed q-values (multiple-testing-corrected probabilities) for bendiocarb-selected Tiassalé samples versus the average of the three susceptible populations; (B) Comparison of Kovié FC against Tiassalé-selected FC for probes significant in both experiments. For genes represented by multiple probes, numbers in parentheses indicate the number of probes significant/total.</p