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.

Windborne long-distance migration of malaria mosquitoes in the Sahel.

Over the past two decades efforts to control malaria have halved the number of cases globally, yet burdens remain high in much of Africa and the elimination of malaria has not been achieved even in areas where extreme reductions have been sustained, such as South Africa1,2. Studies seeking to understand the paradoxical persistence of malaria in areas in which surface water is absent for 3-8 months of the year have suggested that some species of Anopheles mosquito use long-distance migration3. Here we confirm this hypothesis through aerial sampling of mosquitoes at 40-290 m above ground level and provide-to our knowledge-the first evidence of windborne migration of African malaria vectors, and consequently of the pathogens that they transmit. Ten species, including the primary malaria vector Anopheles coluzzii, were identified among 235 anopheline mosquitoes that were captured during 617 nocturnal aerial collections in the Sahel of Mali. Notably, females accounted for more than 80% of all of the mosquitoes that we collected. Of these, 90% had taken a blood meal before their migration, which implies that pathogens are probably transported over long distances by migrating females. The likelihood of capturing Anopheles species increased with altitude (the height of the sampling panel above ground level) and during the wet seasons, but variation between years and localities was minimal. Simulated trajectories of mosquito flights indicated that there would be mean nightly displacements of up to 300 km for 9-h flight durations. Annually, the estimated numbers of mosquitoes at altitude that cross a 100-km line perpendicular to the prevailing wind direction included 81,000 Anopheles gambiae sensu stricto, 6 million A. coluzzii and 44 million Anopheles squamosus. These results provide compelling evidence that millions of malaria vectors that have previously fed on blood frequently migrate over hundreds of kilometres, and thus almost certainly spread malaria over these distances. The successful elimination of malaria may therefore depend on whether the sources of migrant vectors can be identified and controlled.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele-specific expression.

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb-resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and CRM sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele-specific expression (ASE) in hybrids of insecticide susceptible and resistant strains, suggesting cis-regulation is an important mechanism of gene expression regulation in A. gambiae. The genes showing ASE included a higher proportion of Anopheles-specific genes on average younger than genes with balanced allelic expression.

Comparison of different trapping methods to collect malaria vectors indoors and outdoors in western Kenya.

BACKGROUND: Vector surveillance is among the World Health Organization global vector control response (2017-2030) pillars. Human landing catches are a gold standard but difficult to implement and potentially expose collectors to malaria infection. Other methods like light traps, pyrethrum spray catches and aspiration are less expensive and less risky to collectors. METHODS: Three mosquito sampling methods (UV light traps, CDC light traps and Prokopack aspiration) were evaluated against human landing catches (HLC) in two villages of Rarieda sub-county, Siaya County, Kenya. UV-LTs, CDC-LTs and HLCs were conducted hourly between 17:00 and 07:00. Aspiration was done indoors and outdoors between 07:00 and 11:00 a.m. Analyses of mosquito densities, species abundance and sporozoite infectivity were performed across all sampling methods. Species identification PCR and ELISAs were done for Anopheles gambiae and Anopheles funestus complexes and data analysis was done in R. RESULTS: Anopheles mosquitoes sampled from 608 trapping efforts were 5,370 constituting 70.3% Anopheles funestus sensu lato (s.l.), 19.7% Anopheles coustani and 7.2% An. gambiae s.l. 93.8% of An. funestus s.l. were An. funestus sensu stricto (s.s.) and 97.8% of An. gambiae s.l. were Anopheles arabiensis. Only An. funestus were sporozoite positive with 3.1% infection prevalence. Indoors, aspiration captured higher An. funestus (mean = 6.74; RR = 8.83, P < 0.001) then UV-LT (mean = 3.70; RR = 3.97, P < 0.001) and CDC-LT (mean = 1.74; RR = 1.89, P = 0.03) compared to HLC. UV-LT and CDC-LT indoors captured averagely 0.18 An. arabiensis RR = 5.75, P = 0.028 and RR = 5.87, P = 0.028 respectively. Outdoors, UV-LT collected significantly higher Anopheles mosquitoes compared to HLC (An. funestus: RR = 5.18, P < 0.001; An. arabiensis: RR = 15.64, P = 0.009; An. coustani: RR = 11.65, P < 0.001). Anopheles funestus hourly biting indoors in UV-LT and CDC-LT indicated different peaks compared to HLC. CONCLUSIONS: Anopheles funestus remains the predominant mosquito species. More mosquitoes were collected using aspiration, CDC-LTs and UV-LTs indoors and UV-LTs and CD-LTs outdoors compared to HLCs. UV-LTs collected more mosquitoes than CDC-LTs. The varied trends observed at different times of the night suggest that these methods collect mosquitoes with diverse activities and care must be taken when interpreting the results.

LLIN Evaluation in Uganda Project (LLINEUP): modelling the impact of COVID-19-related disruptions on delivery of long-lasting insecticidal nets on malaria indicators in Uganda.

BACKGROUND: Disruptions in malaria control due to COVID-19 mitigation measures were predicted to increase malaria morbidity and mortality in Africa substantially. In Uganda, long-lasting insecticidal nets (LLINs) are distributed nationwide every 3-4 years, but the 2020-2021 campaign was altered because of COVID-19 restrictions so that the timing of delivery of new nets was different from the original plans made by the National Malaria Control Programme. METHODS: A transmission dynamics modelling exercise was conducted to explore how the altered delivery of LLINs in 2020-2021 impacted malaria burden in Uganda. Data were available on the planned LLIN distribution schedule for 2020-2021, and the actual delivery. The transmission model was used to simulate 100 health sub-districts, and parameterized to match understanding of local mosquito bionomics, net use estimates, and seasonal patterns based on data collected in 2017-2019 during a cluster-randomized trial (LLINEUP). Two scenarios were compared; simulated LLIN distributions matching the actual delivery schedule, and a comparable scenario simulating LLIN distributions as originally planned. Model parameters were otherwise matched between simulations. RESULTS: Approximately 70% of the study population received LLINs later than scheduled in 2020-2021, although some areas received LLINs earlier than planned. The model indicates that malaria incidence in 2020 was substantially higher in areas that received LLINs late. In some areas, early distribution of LLINs appeared less effective than the original distribution schedule, possibly due to attrition of LLINs prior to transmission peaks, and waning LLIN efficacy after distribution. On average, the model simulations predicted broadly similar overall mean malaria incidence in 2021 and 2022. After accounting for differences in cluster population size and LLIN distribution dates, no substantial increase in malaria burden was detected. CONCLUSIONS: The model results suggest that the disruptions in the 2020-2021 LLIN distribution campaign in Uganda did not substantially increase malaria burden in the study areas.

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.

Evaluation of community-based vector surveillance system for routine entomological monitoring under low malaria vector densities and high bed net coverage in western Kenya.

BACKGROUND: Entomological surveillance is traditionally conducted by supervised teams of trained technicians. However, it is expensive and limiting in the number of sites visited. Surveillance through community-based collectors (CBC) may be more cost-effective and sustainable for longitudinal entomological monitoring. This study evaluated the efficiency of CBCs in monitoring mosquito densities compared to quality-assured sampling conducted by experienced entomology technicians. METHODS: Entomological surveillance employing CBCs was conducted in eighteen clusters of villages in western Kenya using indoor and outdoor CDC light traps and indoor Prokopack aspiration. Sixty houses in each cluster were enrolled and sampled once every month. Collected mosquitoes were initially identified to the genus level by CBCs, preserved in 70% ethanol and transferred to the laboratory every 2 weeks. Parallel, collections by experienced entomology field technicians were conducted monthly by indoor and outdoor CDC light traps and indoor Prokopack aspiration and served as a quality assurance of the CBCs. RESULTS: Per collection, the CBCs collected 80% fewer Anopheles gambiae sensu lato (s.l.) [RR = 0.2; (95% CI 0.14-0.27)] and Anopheles coustani [RR = 0.2; (95% CI 0.06-0.53)] and 90% fewer Anopheles funestus [RR = 0.1; (95% CI 0.08-0.19)] by CDC light traps compared to the quality assured (QA) entomology teams. Significant positive correlations were however observed between the monthly collections by CBCs and QA teams for both An. gambiae and An. funestus. In paired identifications of pooled mosquitoes, the CBCs identified 4.3 times more Anopheles compared to experienced technicians. The cost per person-night was lower in the community-based sampling at $9.1 compared to $89.3 by QA per collection effort. CONCLUSION: Unsupervised community-based mosquito surveillance collected substantially fewer mosquitoes per trap-night compared to quality-assured collection by experienced field teams, while consistently overestimating the number of Anopheles mosquitoes during identification. However, the numbers collected were significantly correlated between the CBCs and the QA teams suggesting that trends observed by CBCs and QA teams were similar. Further studies are needed to evaluate whether adopting low-cost, devolved supervision with spot checks, coupled with remedial training of the CBCs, can improve community-based collections to be considered a cost-effective alternative to surveillance conducted by experienced entomological technicians.

Late morning biting behaviour of Anopheles funestus is a risk factor for transmission in schools in Siaya, western Kenya.

BACKGROUND: Children in Kenya spend a substantial amount of time at school, including at dawn and dusk when mosquitoes are active. With changing vector behaviour towards early morning biting, it is important to determine whether there is an additional risk of transmission in schools. This study sought to understand whether late morning biting by Anopheles funestus, previously documented in households in western Kenya, was replicated in schools. METHODS: From the 4th to the 6th of August 2023, human landing collections were conducted hourly in four schools in Alego Usonga sub-County, Siaya County. The collections were conducted in and outside five classrooms in each school and ran for 17 h, starting at 18:00 until 11:00 h the next morning. RESULTS: Anopheles funestus was the predominant species collected, forming 93.2% (N = 727) of the entire collection, with peak landing between 06:00 and 07:00 h and continuing until 11:00 h. More than half of the collected An. funestus were either fed or gravid, potentially indicative of multiple bloodmeals within each gonotrophic cycle, and had a sporozoite rate of 2.05%. CONCLUSION: School children spend up to 10 h of their daytime in schools, reporting between 06:00 and 07:00 h and staying in school until as late as 17:00 h, meaning that they receive potentially infectious mosquito bites during the morning hours in these settings. There is a need to consider vector control approaches targeting schools and other peridomestic spaces in the morning hours when An. funestus is active.

Modelling spatiotemporal trends in the frequency of genetic mutations conferring insecticide target-site resistance in African mosquito malaria vector species.

BACKGROUND: Resistance in malaria vectors to pyrethroids, the most widely used class of insecticides for malaria vector control, threatens the continued efficacy of vector control tools. Target-site resistance is an important genetic resistance mechanism caused by mutations in the voltage-gated sodium channel (Vgsc) gene that encodes the pyrethroid target-site. Understanding the geographic distribution of target-site resistance, and temporal trends across different vector species, can inform strategic deployment of vector control tools. RESULTS: We develop a Bayesian statistical spatiotemporal model to interpret species-specific trends in the frequency of the most common resistance mutations, Vgsc-995S and Vgsc-995F, in three major malaria vector species Anopheles gambiae, An. coluzzii, and An. arabiensis over the period 2005-2017. The models are informed by 2418 observations of the frequency of each mutation in field sampled mosquitoes collected from 27 countries spanning western and eastern regions of Africa. For nine selected countries, we develop annual predictive maps which reveal geographically structured patterns of spread of each mutation at regional and continental scales. The results show associations, as well as stark differences, in spread dynamics of the two mutations across the three vector species. The coverage of ITNs was an influential predictor of Vgsc allele frequencies, with modelled relationships between ITN coverage and allele frequencies varying across species and geographic regions. We found that our mapped Vgsc allele frequencies are a significant partial predictor of phenotypic resistance to the pyrethroid deltamethrin in An. gambiae complex populations. CONCLUSIONS: Our predictive maps show how spatiotemporal trends in insecticide target-site resistance mechanisms in African An. gambiae vary across individual vector species and geographic regions. Molecular surveillance of resistance mechanisms will help to predict resistance phenotypes and track their spread.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes.

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

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.

Studies of insecticide resistance provide insights into the capacity of populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs). Resistance in the major malaria vectors has increased in concert with the expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance is crucial for the development and deployment of resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border, with stepwise replacement of the double-mutant (Zanzibar-like TE and Cyp6p4-236 M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolize pyrethroids and is strongly predictive of resistance to pyrethroids most notably deltamethrin. Importantly, there was increased mortality in mosquitoes carrying the triple-mutation when exposed to nets cotreated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid cotreated LLINs across African countries.

Education and Socio-economic status are key factors influencing use of insecticides and malaria knowledge in rural farmers in Southern Côte d'Ivoire.

BACKGROUND: Insecticides play a key role in rural farming; however, their over- or misuse has been linked with a negative impact on malaria vector control policies. This study was conducted amongst agricultural communities in Southern Côte d'Ivoire to identify which insecticides are used by local farmers and how it relates to the perception of farmers on malaria. Understanding the use of insecticides may help in designing awareness programme on mosquito control and pesticides management. METHODS: A questionnaire was administered to 1399 farming households across ten villages. Farmers were interviewed on their education, farming practices (e.g. crops cultivated, insecticides use), perception of malaria, and the different domestic strategies of mosquito control they use. Based on some pre-defined household assets, the socioeconomic status (SES) of each household was estimated. Statistical associations were calculated between different variables, showing significant risk factors. RESULTS: The educational level of farmers was significantly associated with their SES (p < 0.0001). Most of the householders (88.82%) identified mosquitoes as the principal cause of malaria, with good knowledge of malaria resulting as positively related to high educational level (OR = 2.04; 95%CI: 1.35, 3.10). The use of indoor chemical compounds was strongly associated to the SES of the households, their education level, their use of ITNs and insecticide in agricultural (p < 0.0001). Indoor application of pyrethroid insecticides was found to be widespread among farmers as well as the use of such insecticide for crops protection. CONCLUSION: Our study shows that the education level remains the key factor influencing the use of insecticides by farmers and their awareness of malaria control. We suggest that better communication tailored to education level and including SES, controlled availability and access to chemical products, should be considered when designing campaigns on use of pesticides and vector borne disease control for local communities.

Evolution of the Insecticide Target Rdl in African Anopheles Is Driven by Interspecific and Interkaryotypic Introgression.

The evolution of insecticide resistance mechanisms in natural populations of Anopheles malaria vectors is a major public health concern across Africa. Using genome sequence data, we study the evolution of resistance mutations in the resistance to dieldrin locus (Rdl), a GABA receptor targeted by several insecticides, but most notably by the long-discontinued cyclodiene, dieldrin. The two Rdl resistance mutations (296G and 296S) spread across West and Central African Anopheles via two independent hard selective sweeps that included likely compensatory nearby mutations, and were followed by a rare combination of introgression across species (from A. gambiae and A. arabiensis to A. coluzzii) and across nonconcordant karyotypes of the 2La chromosomal inversion. Rdl resistance evolved in the 1950s as the first known adaptation to a large-scale insecticide-based intervention, but the evolutionary lessons from this system highlight contemporary and future dangers for management strategies designed to combat development of resistance in malaria vectors.

Effect of long-lasting insecticidal nets with and without piperonyl butoxide on malaria indicators in Uganda (LLINEUP): a pragmatic, cluster-randomised trial embedded in a national LLIN distribution campaign.

BACKGROUND: Long-lasting insecticidal nets (LLINs) are the primary malaria prevention tool, but their effectiveness is threatened by pyrethroid resistance. We embedded a pragmatic cluster-randomised trial into Uganda's national LLIN campaign to compare conventional LLINs with those containing piperonyl butoxide (PBO), a synergist that can partially restore pyrethroid susceptibility in mosquito vectors. METHODS: 104 health sub-districts, from 48 districts in Uganda, were randomly assigned to LLINs with PBO (PermaNet 3.0 and Olyset Plus) and conventional LLINs (PermaNet 2.0 and Olyset Net) by proportionate randomisation using an iterative process. At baseline 6, 12, and 18 months after LLIN distribution, cross-sectional surveys were done in 50 randomly selected households per cluster (5200 per survey); a subset of ten households per cluster (1040 per survey) were randomly selected for entomological surveys. The primary outcome was parasite prevalence by microscopy in children aged 2-10 years, assessed in the as-treated population at 6, 12, and 18 months. This trial is registered with ISRCTN, ISRCTN17516395. FINDINGS: LLINs were delivered to households from March 25, 2017, to March 18, 2018, 32 clusters were randomly assigned to PermaNet 3.0, 20 to Olyset Plus, 37 to PermaNet 2.0, and 15 to Olyset Net. In the as-treated analysis, three clusters were excluded because no dominant LLIN was received, and four clusters were reassigned, resulting in 49 PBO LLIN clusters (31 received PermaNet 3.0 and 18 received Olyset Plus) and 52 non-PBO LLIN clusters (39 received PermaNet 2.0 and 13 received Olyset Net). At 6 months, parasite prevalence was 11% (386/3614) in the PBO group compared with 15% (556/3844) in the non-PBO group (prevalence ratio [PR] adjusted for baseline values 0·74, 95% CI 0·62-0·87; p=0·0003). Parasite prevalence was similar at month 12 (11% vs 13%; PR 0·73, 95% CI 0·63-0·85; p=0·0001) and month 18 (12% vs 14%; PR 0·84, 95% CI 0·72-0·98; p=0·029). INTERPRETATION: In Uganda, where pyrethroid resistance is high, PBO LLINs reduced parasite prevalence more effectively than did conventional LLINs for up to 18 months. This study provides evidence needed to support WHO's final recommendation on use of PBO LLINs. FUNDING: The Against Malaria Foundation, UK Department for International Development, Innovative Vector Control Consortium, and Bill and Melinda Gates Foundation.

The genetic architecture of target-site resistance to pyrethroid insecticides in the African malaria vectors Anopheles gambiae and Anopheles coluzzii.

Resistance to pyrethroid insecticides is a major concern for malaria vector control. Pyrethroids target the voltage-gated sodium channel (VGSC), an essential component of the mosquito nervous system. Substitutions in the amino acid sequence can induce a resistance phenotype. We use whole-genome sequence data from phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) to provide a comprehensive account of genetic variation in the Vgsc gene across 13 African countries. In addition to known resistance alleles, we describe 20 other non-synonymous nucleotide substitutions at appreciable population frequency and map these variants onto a protein model to investigate the likelihood of pyrethroid resistance phenotypes. Thirteen of these novel alleles were found to occur almost exclusively on haplotypes carrying the known L995F kdr (knock-down resistance) allele and may enhance or compensate for the L995F resistance genotype. A novel mutation I1527T, adjacent to a predicted pyrethroid-binding site, was found in tight linkage with V402L substitutions, similar to allele combinations associated with resistance in other insect species. We also analysed genetic backgrounds carrying resistance alleles, to determine which alleles have experienced recent positive selection, and describe ten distinct haplotype groups carrying known kdr alleles. Five of these groups are observed in more than one country, in one case separated by over 3000 km, providing new information about the potential for the geographical spread of resistance. Our results demonstrate that the molecular basis of target-site pyrethroid resistance in malaria vectors is more complex than previously appreciated, and provide a foundation for the development of new genetic tools for insecticide resistance management.

Impact of seasonality and malaria control interventions on Anopheles density and species composition from three areas of Uganda with differing malaria endemicity.

BACKGROUND: Long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) are the malaria control interventions primarily responsible for reductions in transmission intensity across sub-Saharan Africa. These interventions, however, may have differential impact on Anopheles species composition and density. This study examined the changing pattern of Anopheles species in three areas of Uganda with markedly different transmission intensities and different levels of vector control. METHODS: From October 2011 to June 2016 mosquitoes were collected monthly using CDC light traps from 100 randomly selected households in three areas: Walukuba (low transmission), Kihihi (moderate transmission) and Nagongera (high transmission). LLINs were distributed in November 2013 in Walukuba and Nagongera and in June 2014 in Kihihi. IRS was implemented only in Nagongera, with three rounds of bendiocarb delivered between December 2014 and June 2015. Mosquito species were identified morphologically and by PCR (Polymerase Chain Reaction). RESULTS: In Walukuba, LLIN distribution was associated with a decline in Anopheles funestus vector density (0.07 vs 0.02 mosquitoes per house per night, density ratio [DR] 0.34, 95% CI: 0.18-0.65, p = 0.001), but not Anopheles gambiae sensu stricto (s.s.) nor Anopheles arabiensis. In Kihihi, over 98% of mosquitoes were An. gambiae s.s. and LLIN distribution was associated with a decline in An. gambiae s.s. vector density (4.00 vs 2.46, DR 0.68, 95% CI: 0.49-0.94, p = 0.02). In Nagongera, the combination of LLINs and multiple rounds of IRS was associated with almost complete elimination of An. gambiae s.s. (28.0 vs 0.17, DR 0.004, 95% CI: 0.002-0.009, p < 0.001), and An. funestus sensu lato (s.l.) (3.90 vs 0.006, DR 0.001, 95% CI: 0.0005-0.004, p < 0.001), with a less pronounced decline in An. arabiensis (9.18 vs 2.00, DR 0.15 95% CI: 0.07-0.33, p < 0.001). CONCLUSIONS: LLIN distribution was associated with reductions in An. funestus s.l. in the lowest transmission site and An. gambiae s.s. in the moderate transmission site. In the highest transmission site, a combination of LLINs and multiple rounds of IRS was associated with the near collapse of An. gambiae s.s. and An. funestus s.l. Following IRS, An. arabiensis, a behaviourally resilient vector, became the predominant species, which may have implications for malaria vector control activities. Development of interventions targeted at outdoor biting remains a priority.

LLIN Evaluation in Uganda Project (LLINEUP): factors associated with childhood parasitaemia and anaemia 3 years after a national long-lasting insecticidal net distribution campaign: a cross-sectional survey.

BACKGROUND: Recent reductions in malaria burden have been attributed largely to long-lasting insecticidal nets (LLINs). In March-June 2017, approximately 3 years after a national LLIN distribution campaign, a cross-sectional community survey was conducted to investigate factors associated with malaria parasitaemia and anaemia, in advance of Uganda's 2017-2018 LLIN campaign. METHODS: Households from 104 clusters in 48 districts were randomly selected using two-staged cluster sampling; 50 households were enrolled per cluster. Eligible children aged 2-10 years had blood obtained for a thick blood smear and those aged 2-4 years had haemoglobin measured. Associations between outcomes and variables of interest were assessed using log-binomial regression with generalized estimating equations to adjust for household clustering. RESULTS: In total, 5196 households, 8834 children with blood smear results, and 3753 with haemoglobin results were included. Only 16% of children lived in households with adequate LLIN coverage. Overall, parasite prevalence was 26.0%, ranging from 8.0% in the South West to 53.1% in East Central. Limiting data to children 2-4 years of age, parasite prevalence was 21.4%, up from 16.9% in 2014-2015 following the national LLIN campaign. In a multivariate analysis, factors associated with parasitaemia included region (East-Central vs South-Western; adjusted prevalence ratio [aPR] 6.45, 95% CI 5.55-7.50; p < 0.001), older age (8-10 vs 2-3 years; aPR 1.57, 95% CI 1.43-1.72; p < 0.001), living in a poorer household (poorest vs least poor tercile; aPR 2.32, 95% CI 2.05-2.63; p < 0.001), one constructed of traditional materials (aPR 1.13, 95% CI 1.03-1.24; p = 0.008), or without adequate LLIN coverage (aPR 1.30, 95% CI 1.14-1.48; p < 0.001). Overall, the prevalence of anaemia (haemoglobin < 10 g/dL) was 15.1% and varied geographically. In a multivariate analysis, factors associated with anaemia included region, younger age, living in a traditional house, and parasitaemia, which was the strongest predictor (aPR 2.50, 95% CI 2.12-2.95; p < 0.001). CONCLUSIONS: Three years after a national LLIN campaign, LLIN coverage was low and parasite prevalence had increased. Parasite prevalence varied widely across Uganda; older children, those living in poorer households, and those with inadequate LLIN coverage, were at highest risk of parasitaemia. LLINs may need to be distributed more frequently through mass campaigns or continuously through sustainable mechanisms. Targeting interventions to geographic areas and populations at highest risk should also be considered.

Implications of insecticide resistance for malaria vector control with long-lasting insecticidal nets: evidence from health facility data from Benin.

BACKGROUND: Insecticide-based interventions have averted more than 500 million malaria cases since 2000, but insecticide resistance in mosquitoes could bring about a rebound in disease and mortality. This study investigated whether insecticide resistance was associated with increased incidence of clinical malaria. METHODS: In an area of southern Benin with insecticide resistance and high use of insecticide-treated nets (ITNs), malaria morbidity and insecticide resistance were measured simultaneously in 30 clusters (villages or collections of villages) multiple times over the course of 2 years. Insecticide resistance frequencies were measured using the standard World Health Organization bioassay test. Malaria morbidity was measured by cases recorded at health facilities both in the whole population using routinely collected data and in a passively followed cohort of children under 5 years old. RESULTS: There was no evidence that incidence of malaria from routinely collected data was higher in clusters with resistance frequencies above the median, either in children aged under 5 (RR = 1.27 (95% CI 0.81-2.00) p = 0.276) or in individuals aged 5 or over (RR = 1.74 (95% CI 0.91-3.34) p = 0.093). There was also no evidence that incidence was higher in clusters with resistance frequencies above the median in the passively followed cohort (RR = 1.11 (0.52-2.35) p = 0.777). CONCLUSIONS: This study found no association between frequency of resistance and incidence of clinical malaria in an area where ITNs are the principal form of vector control. This may be because, as other studies have shown, ITNs continue to offer some protection from malaria even in the presence of insecticide resistance. Irrespective of resistance, nets provide only partial protection so the development of improved or supplementary vector control tools is required to reduce Africa's unacceptably high malaria burden.

Genome-wide transcriptional analyses in Anopheles mosquitoes reveal an unexpected association between salivary gland gene expression and insecticide resistance.

BACKGROUND: To combat malaria transmission, the Ugandan government has embarked upon an ambitious programme of indoor residual spraying (IRS) with a carbamate class insecticide, bendiocarb. In preparation for this campaign, we characterized bendiocarb resistance and associated transcriptional variation among Anopheles gambiae s.s. mosquitoes from two sites in Uganda. RESULTS: Gene expression in two mosquito populations displaying some resistance to bendiocarb (95% and 79% An. gambiae s.l. WHO tube bioassay mortality in Nagongera and Kihihi, respectively) was investigated using whole-genome microarrays. Significant overexpression of several genes encoding salivary gland proteins, including D7r2 and D7r4, was detected in mosquitoes from Nagongera. In Kihihi, D7r4, two detoxification-associated genes (Cyp6m2 and Gstd3) and an epithelial serine protease were among the genes most highly overexpressed in resistant mosquitoes. Following the first round of IRS in Nagongera, bendiocarb-resistant mosquitoes were collected, and real-time quantitative PCR analyses detected significant overexpression of D7r2 and D7r4 in resistant mosquitoes. A single nucleotide polymorphism located in a non-coding transcript downstream of the D7 genes was found at a significantly higher frequency in resistant individuals. In silico modelling of the interaction between D7r4 and bendiocarb demonstrated similarity between the insecticide and serotonin, a known ligand of D7 proteins. A meta-analysis of published microarray studies revealed a recurring association between D7 expression and insecticide resistance across Anopheles species and locations. CONCLUSIONS: A whole-genome microarray approach identified an association between novel insecticide resistance candidates and bendiocarb resistance in Uganda. In addition, a single nucleotide polymorphism associated with this resistance mechanism was discovered. The use of such impartial screening methods allows for discovery of resistance candidates that have no previously-ascribed function in insecticide binding or detoxification. Characterizing these novel candidates will broaden our understanding of resistance mechanisms and yield new strategies for combatting widespread insecticide resistance among malaria vectors.

Insecticide resistance in Anopheles gambiae from the northern Democratic Republic of Congo, with extreme knockdown resistance (kdr) mutation frequencies revealed by a new diagnostic assay.

BACKGROUND: Mutations in the voltage-gated sodium channel at codon 1014 confer knock-down resistance (kdr) to pyrethroids in a wide range of insects. Anopheles gambiae exhibits two mutant alleles at codon 1014, serine and phenylalanine; and both are now widespread across Africa. Existing screening methods only allow for one resistant allele to be detected per assay. A new locked nucleic acid (LNA) qPCR assay was developed for the simultaneous detection of both mutant alleles and the wild type allele in a single assay. This tri-allelic detection assay was assessed as part of a study of the insecticide resistance in An. gambiae sensu stricto (s.s.) in the previously un-sampled area of Nord Ubangi, Democratic Republic of the Congo. METHODS: Samples from three sites were tested for insecticide susceptibility using WHO bioassays, with and without the synergist PBO preceding pyrethroid exposures, and were subsequently analysed for frequency and resistance-association of the Vgsc-1014 and Vgsc-N1575Y mutations. Results from the LNA-kdr 1014 assay were compared to results from standard TaqMan-kdr assays. RESULTS: Anopheles gambiae sensu lato (s.l.) was by far the predominant vector captured (84%), with only low frequencies of Anopheles funestus s.l. (9%) detected in Nord Ubangi. Molecular identification found An. gambiae s.s. to be the principal vector (99%) although Anopheles coluzzii was detected at very low frequency. Anopheles gambiae were susceptible to the carbamate insecticide bendiocarb, but resistant to DDT and to the pyrethroids permethrin and deltamethrin. Susceptibility to both pyrethroids was partially restored with prior exposure to PBO suggesting likely involvement of metabolic resistance. Anopheles gambiae s.s. was homozygous for kdr resistant alleles with both the L1014F and L1014S mutations present, and the N1575Y polymorphism was present at low frequency. The LNA-kdr assay simultaneously detected both resistant alleles and gave results entirely consistent with those from the two TaqMan-kdr assays. CONCLUSION: This study provides rare data on insecticide resistance and mechanisms in Anopheles from the centre of Africa, with the first detection of N1575Y. Nord Ubangi populations of An. gambiae s.s. show insecticide resistance mediated by both metabolic mechanisms and Vgsc mutations. The LNA-kdr assay is particularly suitable for use in populations in which both 1014S and 1014F kdr alleles co-occur and provides robust results, with higher throughput and at a quarter of the cost of TaqMan assays.