Insecticide resistance patterns in Uganda and the effect of indoor residual spraying with bendiocarb on kdr L1014S frequencies in Anopheles gambiae s.s

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

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

Insights and challenges of insecticide resistance modelling in malaria vectors: a review

BACKGROUND

Malaria is one of the most devastating tropical diseases, resulting in loss of lives each year, especially in children under the age of 5 years. Malaria burden, related deaths and stall in the progress against malaria transmission is evident, particularly in countries that have moderate or high malaria transmission. Hence, mitigating malaria spread requires information on the distribution of vectors and the drivers of insecticide resistance (IR). However, owing to the impracticality in establishing the critical need for real-world information at every location, modelling provides an informed best guess for such information. Therefore, this review examines the various methodologies used to model spatial, temporal and spatio-temporal patterns of IR within populations of malaria vectors, incorporating pest-biology parameters, adopted ecological principles, and the associated modelling challenges.

METHODS

The review focused on the period ending March 2023 without imposing restrictions on the initial year of publication, and included articles sourced from PubMed, Web of Science, and Scopus. It was also limited to publications that deal with modelling of IR distribution across spatial and temporal dimensions and excluded articles solely focusing on insecticide susceptibility tests or articles not published in English. After rigorous selection, 33 articles met the review's elibility criteria and were subjected to full-text screening.

RESULTS

Results show the popularity of Bayesian geostatistical approaches, and logistic and static models, with limited adoption of dynamic modelling approaches for spatial and temporal IR modelling. Furthermore, our review identifies the availability of surveillance data and scarcity of comprehensive information on the potential drivers of IR as major impediments to developing holistic models of IR evolution.

CONCLUSIONS

The review notes that incorporating pest-biology parameters, and ecological principles into IR models, in tandem with fundamental ecological concepts, potentially offers crucial insights into the evolution of IR. The results extend our knowledge of IR models that provide potentially accurate results, which can be translated into policy recommendations to combat the challenge of IR in malaria control.

Sibling Species Composition and Susceptibility Status of Anopheles gambiae s.l. to Insecticides Used for Indoor Residual Spraying in Eastern Uganda

Background

Malaria remains one of the most critical disease causing morbidity and mortality in Uganda. Indoor residual spraying (IRS) and the use of insecticide-treated bed nets are currently the predominant malaria vector control interventions. However, the emergence and spread of insecticide resistance among malaria vectors threaten the continued effectiveness of these interventions to control the disease, particularly in high transmission areas. To inform decisions on vector control, the current study evaluated the Anopheles malaria vector species and their susceptibility levels to 0.1% bendiocarb and 0.25% pirimiphos-methyl insecticides used in IRS intervention program in Namutumba district, Eastern Uganda.

Methods

Anopheles larvae were collected between March and May 2017 from different breeding sites in the parishes of Nsinze and Nawaikona in Nsinze sub-county and reared to adults to assess the susceptibility status of populations in the study area. Mosquitoes were identified using morphological keys and species-specific polymerase chain reaction (PCR) assays. Susceptibility tests were conducted on 2- to 5-day-old non-blood-fed adult female Anopheles that emerged using insecticide-impregnated papers with 0.1% bendiocarb and 0.25% pirimiphos-methyl following standard World Health Organization (WHO) insecticide susceptibility bioassays. A Log-probit regression model was used to derive the knock-down rates for 50% and 95% of exposed mosquitoes.

Results

A total of 700 mosquito larvae were collected from different breeding sites. Morphological identification showed that 500 individuals that emerged belonged to Anopheles gambiae sensu lato (s.l.), the main malaria vector. The PCR results showed that the dominant sibling species under the A. gambiae complex was Anopheles arabiensis 99.5% (395/397). WHO bioassay tests revealed that the population of mosquitoes exhibited high levels of susceptibility (24-hour post-exposure mortality 98-100%) to both insecticides tested. The median knock-down time, KDT₅₀, ranged from 6.6 to 81.4 minutes, while the KDT₉₅ ranged from 21.6 to 118.9 minutes for 0.25% pirimiphos-methyl. The KDT₅₀ for 0.1% bendiocarb ranged from 2.8 to 62.9 minutes, whereas the KDT₉₅ ranged from 36.0 to 88.5 minutes.

Conclusions

These findings indicate that bendiocarb and pirimiphos-methyl are still effective against the major malaria vector, A. arabiensis in Nsinze sub-county, Namutumba district, Uganda and can be effectively used for IRS. The study has provided baseline information on the insecticide susceptibility status on malaria vectors in the study area. However, routine continuous monitoring program of insecticide susceptibility and malaria vector composition is required so as to guide future decisions on insecticide use for IRS intervention toward malaria elimination and to track future changes in vector population.

In Vitro and In Vivo Validation of CYP6A14 and CYP6N6 Participation in Deltamethrin Metabolic Resistance in Aedes albopictus

The extensive use of chemical insecticides for public health and agricultural purposes has increased the occurrence and development of insecticide resistance. This study used transcriptome sequencing to screen 10 upregulated metabolic detoxification enzyme genes from Aedes albopictus resistant strains. Of these, CYP6A14 and CYP6N6 were found to be substantially overexpressed in the deltamethrin-induced expression test, indicating their role in deltamethrin resistance in Ae. albopictus. Furthermore, the corresponding 60-kDa recombinant proteins, CYP6A14 and CYP6N6, were successfully expressed using the Escherichia coli expression system. Enzyme activity studies revealed that CYP6A14 (5.84 U/L) and CYP6N6 (6.3 U/L) have cytochrome P450 (CYP450) enzyme activity. In vitro, the metabolic analysis revealed that the recombinant proteins degraded deltamethrin into 1-oleoyl-sn-glycero-3-phosphoethanolamine and 2',2'-dibromo-2'-deoxyguanosine. Subsequently, the CYP450 genes in larvae of Ae. albopictus were silenced by RNA interference technology to study deltamethrin resistance in vivo. The silencing of CYP6A14 and CYP6N6 increased the mortality rate of mosquitoes without affecting their survival time, spawning quantity, hatching rate, and other normal life activities. Altogether, CYP6A14 and CYP6N6 belong to the CYP6 family and mutually increase deltamethrin resistance in Ae. albopictus.

Colonization and Authentication of the Pyrethroid-Resistant Anopheles gambiae s.s. Muleba-Kis Strain; an Important Test System for Laboratory Screening of New Insecticides

Simple Summary

Malaria control and prevention have traditionally relied on the use of insecticides in the form of treated bed nets or residual spraying in households. However, scaling up of these interventions—based on few available insecticide classes—resulted in the development and spread of insecticide resistance in malaria-transmitting mosquitoes. There is therefore an urgent need for introducing and applying new insecticides that are effective against these mosquitoes. Laboratories tasked with evaluating the efficacy of novel insecticides need to establish a large colony of resistant mosquitoes. In this study, we report the procedures used and challenges faced during the establishment and maintenance of a resistant mosquito strain in the laboratory which reflects the characteristics of the wild-resistant mosquito populations found in East Africa.

Abstract

Background: The emergence and spread of insecticide resistance in malaria vectors to major classes of insecticides call for urgent innovation and application of insecticides with novel modes of action. When evaluating new insecticides for public health, potential candidates need to be screened against both susceptible and resistant mosquitoes to determine efficacy and to identify potential cross-resistance to insecticides currently used for mosquito control. The challenges and lessons learned from establishing, maintaining, and authenticating the pyrethroid-resistant An. gambiae s.s. Muleba-Kis strain at the KCMUCo-PAMVERC Test Facility are described in this paper. Methods: Male mosquitoes from the F₁ generation of wild-pyrethroid resistant mosquitoes were cross-bred with susceptible female An. gambiae s.s. Kisumu laboratory strain followed by larval selection using a pyrethroid insecticide solution. Periodic screening for phenotypic and genotypic resistance was done. WHO susceptibility tests and bottle bioassays were used to assess the phenotypic resistance, while Taqman™ assays were used to screen for known target-site resistance alleles (kdr and ace-1). Additionally, the strains were periodically assessed for quality control by monitoring adult weight and wing length. Results: By out-crossing the wild mosquitoes with an established lab strain, a successful resistant insectary colony was established. Intermittent selection pressure using alphacypermethrin has maintained high kdr mutation (leucine-serine) frequencies in the selected colony. There was consistency in the wing length and weight measurements from the year 2016 to 2020, with the exception that one out of four years was significantly different. Mean annual wing length varied between 0.0142–0.0028 mm compared to values obtained in 2016, except in 2019 where it varied by 0.0901 mm. Weight only varied by approximately 0.001 g across four years, except in 2017 where it differed by 0.005 g. Routine phenotypic characterization on Muleba-Kis against pyrethroids using the WHO susceptibility test indicated high susceptibility when type I pyrethroids were used compared to type II pyrethroids. Dynamics on susceptibility status also depended on the lapse time when the selection was last done. Conclusions: This study described the procedure for introducing, colonizing, and maintaining a resistant An. gambiae s.s. strain in the laboratory with leucine to serine substitution kdr allele which reflects the features of the wild-resistant population in East Africa. Challenges in colonizing a wild-resistant mosquito strain were overcome by out-crossing between mosquito strains of desired traits followed by intermittent insecticide selection at the larval stage to select for the resistant phenotype.

Anopheline Mosquito Species Composition, Kdr Mutation Frequency, and Parasite Infectivity Status in Northern Tanzania

The scaling-up of malaria control interventions in northern Tanzania has resulted in a decline in malaria prevalence and vector species composition. Despite this achievement, residual malaria transmission remains a concern in the area. The main aim of this study was to investigate malaria vector species composition, parasite infectivity rates, and the presence of insecticide knockdown resistance (kdr) mutations in three sites that have experienced a significant decline in malaria in northern Tanzania. Adult mosquitoes were sampled using light traps in houses and hand-aspirators in cowsheds, whereas the standard dipping method was used for sampling mosquito larvae. Adult mosquitoes identified as Anopheles gambiae s.l. and An. funestus s.l. and larval stages III and IV of An. gambiae s.l. were stored in absolute ethanol for further laboratory molecular identification. The identified species in the An. gambiae complex were An. gambiae s.s., An. merus, An. quadriannulatus, and An. arabiensis, whereas the An. funestus group comprised An. funestus s.s., An. rivulorum, and An. leesoni. For An. gambiae s.s. analyzed from Zeneth, 47.6% were kdr-East homozygous susceptible, 35.7% kdr-East heterozygous resistant, 9.6% kdr-East homozygous resistant, and 7.1% undefined, whereas specimens from Kwakibuyu were 45.5% kdr-East homozygous susceptible, 32.7% kdr-East heterozygous resistant, 16.3% kdr-East homozygous resistant, and 5.5% undefined. There were no kdr-West alleles identified from any specimen. The overall malaria parasite infectivity rate was 0.75%. No infections were found in Moshi. The findings indicate that populations of the major malaria vector mosquitoes are still present in the study area, with An. funestus taking a lead in malaria transmission.

Transcriptomic analysis of insecticide resistance in the lymphatic filariasis vector Culex quinquefasciatus

Culex quinquefasciatus plays an important role in transmission of vector-borne diseases of public health importance, including lymphatic filariasis (LF), as well as many arboviral diseases. Currently, efforts to tackle C. quinquefasciatus vectored diseases are based on either mass drug administration (MDA) for LF, or insecticide-based interventions. Widespread and intensive insecticide usage has resulted in increased resistance in mosquito vectors, including C. quinquefasciatus. Herein, the transcriptome profile of Ugandan bendiocarb-resistant C. quinquefasciatus was explored to identify candidate genes associated with insecticide resistance. High levels of insecticide resistance were observed for five out of six insecticides tested, with the lowest mortality (0.97%) reported to permethrin, while for DDT, lambdacyhalothrin, bendiocarb and deltamethrin the mortality rate ranged from 1.63-3.29%. Resistance to bendiocarb in exposed mosquitoes was marked, with 2.04% mortality following 1 h exposure and 58.02% after 4 h. Genotyping of the G119S Ace-1 target site mutation detected a highly significant association (p < 0.0001; OR = 25) between resistance and Ace1-119S. However, synergist assays using the P450 inhibitor PBO, or the esterase inhibitor TPP resulted in markedly increased mortality (to ≈80%), suggesting a role of metabolic resistance in the resistance phenotype. Using a novel, custom 60 K whole-transcriptome microarray 16 genes significantly overexpressed in resistant mosquitoes were detected, with the P450 Cyp6z18 showing the highest differential gene expression (>8-fold increase vs unexposed controls). These results provide evidence that bendiocarb resistance in Ugandan C. quinquefasciatus is mediated by both target-site mechanisms and over-expression of detoxification enzymes.

LLIN Evaluation in Uganda Project (LLINEUP): a cross-sectional survey of species diversity and insecticide resistance in 48 districts of Uganda

BACKGROUND

Long-lasting insecticidal nets (LLINs) are the principal tool for malaria control in Africa and are presently treated with a single class of insecticide; however, increasing levels of insecticide resistance threaten their success. In response to this threat nets have been developed that incorporate the synergist, piperonyl butoxide (PBO), which inhibits the activity of cytochrome P450s which is one main mechanisms of insecticide resistance, allowing resistance to pyrethroids to be reversed. However, data on the value and cost effectiveness of these nets is lacking. A large-scale cluster randomised trial of conventional LLINs and PBO-LLINs was conducted in Uganda in 104 health sub-districts (HSDs) in 2017-2019. Prior to the mass distribution of LLINs, a baseline entomological survey was carried out, the results of which are reported herein. Ten households from each HSD were randomly selected for entomological surveillance at baseline which included household mosquito collections.

RESULTS

Prior to LLIN distribution entomological collections were carried out in 1029 houses across the 104 HSDs. Anopheles gambiae (s.l.) was the principal vector in all but 9 of the 71 HSDs that yielded vector species. Molecular analysis found An. gambiae (s.s.) to be the predominant vector collected. Plasmodium falciparum was detected in 5.5% of An. gambiae (s.s.) and in 4.0% of An. funestus (s.s.) examined. Infection rates of other plasmodium species (P. vivax, P. ovale and P. malariae) were lower with infection rates of 1.2% and 1.7% for An. gambiae (s.s.) and An. funestus (s.s.), respectively. The knockdown resistance (kdr) mutation Vgsc-L1014S was found at very high frequency in An. gambiae (s.s.) with the Vgsc-L1014F mutation at low frequency and the wild-type allele virtually absent. In An. arabiensis the wild-type allele was predominant. The resistance-associated alleles, Cyp4j5-L43F and Coeae1d were found at moderate frequencies which varied across the study site. Vgsc-N1575Y mutation was not found in any samples examined.

CONCLUSIONS

No significant differences between planned intervention arms was observed in vector densities, sporozoite infection rate or insecticide resistance marker frequency across the study site prior to the distribution of LLINs. Very high levels of kdr resistance were observed in all areas; however, the resistance-associated markers Cyp4j5-L43F and Coeae1d were found at varying frequencies across the study site which may have implications for the effectiveness of standard LLINs. Trial registration This study is registered with ISRCTN, ISRCTN17516395. Registered 14 February 2017, http://www.isrctn.com/ISRCTN17516395.

The 2015-2016 malaria epidemic in Northern Uganda; What are the implications for malaria control interventions?

Vector control and effective case management are currently the backbone strategies of malaria control. Kitgum district, an area of perennial holoendemic malaria transmission intensity in Northern Uganda, appears to have experienced a malaria epidemic in 2015. This study aimed to describe the malaria trends in Kitgum General Hospital from 2011 to 2017 in relation to climatic factors and the application of population-based malaria control interventions. Hospital records were examined retrospectively to calculate malaria normal channels, malaria cases per 1000 population, test positivity rates (TPR) and to enumerate pregnancy malaria, hospitalizations and deaths. Climatic factors (humidity, temperature and rainfall) and population-based malaria control interventions that had been applied during this period were described. Kitgum district experienced an epidemic between the years 2015 and 2016; the malaria burden rose above the established normal channels. At its peak the number of malaria cases attending KGH was over 20 times above the normal channels. The total number of cases per 1000 population increased from 7 in 2014 to 113 in 2015 and 114 in 2016 (p value for trend <0.0001). Similarly, TPR increased from 10.5% to 54.6% between 2014 and 2016 (p value for trend <0.0001). This trend was also observed for malaria attributable hospitalizations, and malaria in pregnancy. There were no significant changes in any of the climatic factors assessed (p value = 0.92, 0.99, 0.52 for relative humidity, max temperature, and rainfall, respectively). The malaria upsurge occurred in conjunction with a general decline in the use and application of malaria control interventions. Specifically, indoor residual spraying was interrupted in 2014. In response to the epidemic, IRS was reapplied together with mass distribution of long-lasting insecticide treated nets (LLINs) in 2017. Subsequently, there was a decline in all malaria indicators. The epidemic in Kitgum occurred in association with the interruption of IRS and appears to have abated following its re-introduction alongside LLINs. The study suggests that to enable malaria elimination in areas of high malaria transmission intensity, effective control measures may need to be sustained for the long-term.

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations

BACKGROUND

Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species.

METHODOLOGY/PRINCIPAL FINDINGS

Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci).

CONCLUSIONS/SIGNIFICANCE

Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex.