A cytochrome P450 allele confers pyrethroid resistance on a major African malaria vector, reducing insecticide-treated bednet efficacy

Distribution and insecticide resistance profile of the major malaria vector Anopheles funestus group across the African continent

There has been significant progress in malaria control in the last 2 decades, with a decline in mortality and morbidity. However, these gains are jeopardised by insecticide resistance, which negatively impacts the core interventions, such as insecticide-treated nets (ITN) and indoor residual spraying (IRS). While most malaria control and research efforts are still focused on Anopheles gambiae complex mosquitoes, Anopheles funestus remains an important vector in many countries and, in some cases, contributes to most of the local transmission. As countries move towards malaria elimination, it is important to ensure that all dominant vector species, including An. funestus, an important vector in some countries, are targeted. The objective of this review is to compile and discuss information related to A. funestus populations' resistance to insecticides and the mechanisms involved across Africa, emphasising the sibling species and their resistance profiles in relation to malaria elimination goals. Data on insecticide resistance in An. funestus malaria vectors in Africa were extracted from published studies. Online bibliographic databases, including Google Scholar and PubMed, were used to search for relevant studies. Articles published between 2000 and May 2023 reporting resistance of An. funestus to insecticides and associated mechanisms were included. Those reporting only bionomics were excluded. Spatial variation in species distribution and resistance to insecticides was recorded from 174 articles that met the selection criteria. It was found that An. funestus was increasingly resistant to the four classes of insecticides recommended by the World Health Organisation for malaria vector control; however, this varied by country. Insecticide resistance appears to reduce the effectiveness of vector control methods, particularly IRS and ITN. Biochemical resistance due to detoxification enzymes (P450s and glutathione-S-transferases [GSTs]) in An. funestus was widely recorded. However, An. funestus in Africa remains susceptible to other insecticide classes, such as organophosphates and neonicotinoids. This review highlights the increasing insecticide resistance of An. funestus mosquitoes, which are important malaria vectors in Africa, posing a significant challenge to malaria control efforts. While An. funestus has shown resistance to the recommended insecticide classes, notably pyrethroids and, in some cases, organochlorides and carbamates, it remains susceptible to other classes of insecticides such as organophosphates and neonicotinoids, providing potential alternative options for vector control strategies. The study underscores the need for targeted interventions that consider the population structure and geographical distribution of An. funestus, including its sibling species and their insecticide resistance profiles, to effectively achieve malaria elimination goals.

Genetic markers associated with the widespread insecticide resistance in malaria vector Anopheles funestus populations across Tanzania

BACKGROUND

Anopheles funestus is a leading vector of malaria in most parts of East and Southern Africa, yet its ecology and responses to vector control remain poorly understood compared with other vectors such as Anopheles gambiae and Anopheles arabiensis. This study presents the first large-scale survey of the genetic and phenotypic expression of insecticide resistance in An. funestus populations in Tanzania.

METHODS

We performed insecticide susceptibility bioassays on An. funestus mosquitoes in nine regions with moderate-to-high malaria prevalence in Tanzania, followed by genotyping for resistance-associated mutations (CYP6P9a, CYP6P9b, L119F-GSTe2) and structural variants (SV4.3 kb, SV6.5 kb). Generalized linear models were used to assess relationships between genetic markers and phenotypic resistance. An interactive R Shiny tool was created to visualize the data and support evidence-based interventions.

RESULTS

Pyrethroid resistance was universal but reversible by piperonyl-butoxide (PBO). However, carbamate resistance was observed in only five of the nine districts, and dichloro-diphenyl-trichloroethane (DDT) resistance was found only in the Kilombero valley, south-eastern Tanzania. Conversely, there was universal susceptibility to the organophosphate pirimiphos-methyl in all sites. Genetic markers of resistance had distinct geographical patterns, with CYP6P9a-R and CYP6P9b-R alleles, and the SV6.5 kb structural variant absent or undetectable in the north-west but prevalent in all other sites, while SV4.3 kb was prevalent in the north-western and western regions but absent elsewhere. Emergent L119F-GSTe2, associated with deltamethrin resistance, was detected in heterozygous form in districts bordering Mozambique, Malawi and the Democratic Republic of Congo. The resistance landscape was most complex in western Tanzania, in Tanganyika district, where all five genetic markers were detected. There was a notable south-to-north spread of resistance genes, especially CYP6P9a-R, though this appears to be interrupted, possibly by the Rift Valley.

CONCLUSIONS

This study underscores the need to expand resistance monitoring to include An. funestus alongside other vector species, and to screen for both the genetic and phenotypic signatures of resistance. The findings can be visualized online via an interactive user interface and could inform data-driven decision-making for resistance management and vector control. Since this was the first large-scale survey of resistance in Tanzania's An. funestus, we recommend regular updates with greater geographical and temporal coverage.

Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control

Insecticide resistance in malaria vector populations poses a major threat to malaria control, which relies largely on insecticidal interventions. Contemporary vector-control strategies focus on combatting resistance using multiple insecticides with differing modes of action within the mosquito. However, diverse genetic resistance mechanisms are present in vector populations, and continue to evolve. Knowledge of the spatial distribution of these genetic mechanisms, and how they impact the efficacy of different insecticidal products, is critical to inform intervention deployment decisions. We developed a catalogue of genetic-resistance mechanisms in African malaria vectors that could guide molecular surveillance. We highlight situations where intervention deployment has led to resistance evolution and spread, and identify challenges in understanding and mitigating the epidemiological impacts of resistance.

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

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

Adult mosquitoes of the sibling species Anopheles gambiae and Anopheles coluzzii exhibit contrasting patterns of susceptibility to four neonicotinoid insecticides along an urban-to-rural gradient in Yaoundé, Cameroon

BACKGROUND

Neonicotinoids are potential alternatives for controlling pyrethroid-resistant mosquitoes, but their efficacy against malaria vector populations of sub-Saharan Africa has yet to be investigated. The aim of the present study was to test the efficacy of four neonicotinoids against adult populations of the sibling species Anopheles gambiae and Anopheles coluzzii sampled along an urban-to-rural gradient.

METHODS

The lethal toxicity of three active ingredients for adults of two susceptible Anopheles strains was assessed using concentration-response assays, and their discriminating concentrations were calculated. The discriminating concentrations were then used to test the susceptibility of An. gambiae and An. coluzzii mosquitoes collected from urban, suburban and rural areas of Yaoundé, Cameroon, to acetamiprid, imidacloprid, clothianidin and thiamethoxam.

RESULTS

Lethal concentrations of neonicotinoids were relatively high suggesting that this class of insecticides has low toxicity against Anopheles mosquitoes. Reduced susceptibility to the four neonicotinoids tested was detected in An. gambiae populations collected from rural and suburban areas. By contrast, adults of An. coluzzii that occurred in urbanized settings were susceptible to neonicotinoids except acetamiprid for which 80% mortality was obtained within 72 h of insecticide exposure. The cytochrome inhibitor, piperonyl butoxide (PBO), significantly enhanced the activity of clothianidin and acetamiprid against An. gambiae mosquitoes.

CONCLUSIONS

These findings corroborate susceptibility profiles observed in larvae and highlight a significant variation in tolerance to neonicotinoids between An. gambiae and An. coluzzii populations from Yaoundé. Further studies are needed to disentangle the role of exposure to agricultural pesticides and of cross-resistance mechanisms in the development of neonicotinoid resistance in some Anopheles species.

Overexpression and nonsynonymous mutations of UDP-glycosyltransferases are potentially associated with pyrethroid resistance in Anopheles funestus

UDP-glycosyltransferases (UGTs) enzymes are pivotal in insecticide resistance by transforming hydrophobic substrates into more hydrophilic forms for efficient cell elimination. This study provides the first comprehensive investigation of Anopheles funestus UGT genes, their evolution, and their association with pyrethroid resistance. We employed a genome-wide association study using pooled sequencing (GWAS-PoolSeq) and transcriptomics on pyrethroid-resistant An. funestus, along with deep-targeted sequencing of UGTs in 80 mosquitoes Africa-wide. UGT310B2 was consistently overexpressed Africa-wide and significant gene-wise Fst differentiation was observed between resistant and susceptible populations: UGT301C2 and UGT302A3 in Malawi, and UGT306C2 in Uganda. Additionally, nonsynonymous mutations in UGT genes were identified. Gene-wise Tajima's D density curves provide insights into population structures within populations across these countries, supporting previous observations. These findings have important implications for current An. funestus control strategies facilitating the prediction of cross-resistance to other UGT-metabolised polar insecticides, thereby guiding more effective and targeted insecticide resistance management efforts.

Exploring the molecular mechanisms of increased intensity of pyrethroid resistance in Central African population of a major malaria vector Anopheles coluzzii

Molecular mechanisms driving the escalation of pyrethroid resistance in the major malaria mosquitoes of Central Africa remain largely uncharacterized, hindering effective management strategies. Here, resistance intensity and the molecular mechanisms driving it were investigated in a population of Anopheles coluzzii from northern Cameroon. High levels of pyrethroid and organochloride resistance were observed in An. coluzzii population, with no mortality for 1× permethrin; only 11% and 33% mortalities for 5× and 10× permethrin diagnostic concentrations, and <2% mortalities for deltamethrin and DDT, respectively. Moderate bendiocarb resistance (88% mortality) and full susceptibility to malathion were observed. Synergist bioassays with piperonyl butoxide recovered permethrin susceptibility, with mortalities increasing to 53.39%, and 87.30% for 5× and 10× permethrin, respectively, implicating P450 monooxygenases. Synergist bioassays with diethyl maleate (DEM) recovered permethrin and DDT susceptibilities (mortalities increasing to 34.75% and 14.88%, respectively), implicating glutathione S-transferases. RNA-seq-based genome-wide transcriptional analyses supported by quantitative PCR identified glutathione S-transferase, GSTe2 (RNA-seqFC = 2.93 and qRT-PCRFC = 8.4, p < 0.0043) and CYP450, CYP6Z2 (RNA-seqFC = 2.39 and qRT-PCRFC = 11.7, p < 0.0177) as the most overexpressed detoxification genes in the pyrethroid-resistant mosquitoes, compared to mosquitoes of the susceptible Ngousso colony. Other overexpressed genes include P450s, CYP6M2 (FC = 1.68, p < 0.0114), CYP4G16 (FC = 2.02, p < 0.0005), and CYP4G17 (FC = 1.86, p < 0.0276). While high frequency of the 1014F kdr mutation (50%) and low frequencies of 1014S (6.61%) and 1575Y (10.29%) were observed, no ace-1 mutation was detected in bendiocarb-resistant populations, suggesting the preeminent role of metabolic mechanism. Overexpression of metabolic resistance genes (including GSTe2 and CYP6Z2 known to confer resistance to multiple insecticides) in An. coluzzii from the Sudan Savannah of Cameroon highlights the need for alternative management strategies to reduce malaria burden in northern Cameroon.

Clothianidin-resistant Anopheles gambiae adult mosquitoes from Yaoundé, Cameroon, display reduced susceptibility to SumiShield® 50WG, a neonicotinoid formulation for indoor residual spraying

BACKGROUND

Chronic exposure of mosquito larvae to pesticide residues and cross-resistance mechanisms are major drivers of tolerance to insecticides used for vector control. This presents a concern for the efficacy of clothianidin, an agricultural neonicotinoid prequalified for Indoor Residual Spraying (IRS).

METHODS

Using standard bioassays, we tested if reduced susceptibility to clothianidin can affect the efficacy of SumiShield® 50WG, one of four new IRS formulations containing clothianidin. We simultaneously monitored susceptibility to clothianidin and to SumiShield 50WG, testing adults of Anopheles gambiae, An. coluzzii and Culex sp sampled from urban, suburban and agricultural areas of Yaoundé, Cameroon.

RESULTS

We found that in this geographic area, the level of susceptibility to the active ingredient predicted the efficacy of SumiShield 50WG. This formulation was very potent against populations that reached 100% mortality within 72 h of exposure to a discriminating concentration of clothianidin. By contrast, mortality leveled off at 75.4 ± 3.5% within 7 days of exposure to SumiShield 50WG in An. gambiae adults collected from a farm where the spraying of the two neonicotinoids acetamiprid and imidacloprid for crop protection is likely driving resistance to clothianidin.

CONCLUSIONS

Despite the relatively small geographic extend of the study, the findings suggest that cross-resistance may impact the efficacy of some new IRS formulations and that alternative compounds could be prioritized in areas where neonicotinoid resistance is emerging.

High genetic and haplotype diversity in vaccine candidate Pfceltos but not Pfrh5 among malaria-infected children in Ibadan, Nigeria

Malaria remains a global public health challenge. The disease has a great impact in sub-Saharan Africa among children under five years of age and pregnant women. Malaria control programs targeting the parasite and mosquitoes vectors with combinational therapy and insecticide-treated bednets are becoming obsolete due to the phenomenon of resistance, which is a challenge for reducing morbidity and mortality. Malaria vaccines would be effective alternative to the problem of parasite and insecticide resistance, but focal reports of polymorphisms in malaria candidate antigens have made it difficult to design an effective malaria vaccine. Therefore, studies geared towards elucidating the polymorphic pattern and how genes targeted for vaccine design evolve are imperative. We have carried out molecular and genetic analysis of two genes encoding vaccine candidates-the Plasmodium falciparum cell traversal ookinetes and sporozoites (Pfceltos) and P. falciparum reticulocyte binding protein 5 (Pfrh5) in parasite isolates from malaria-infected children in Ibadan, Nigeria to evaluate their genetic diversity, relatedness and pattern of molecular evolution. Pfceltos and Pfrh5 genes were amplified from P. falciparum positive samples. Amplified fragments were purified and sequenced using the chain termination method. Post-sequence edit of fragments and application of various population genetic analyses was done. We observed a higher number of segregating sites and haplotypes in the Pfceltos than in Pfrh5 gene, the former also presenting higher haplotype (0.942) and nucleotide diversity (θ = 0.01219 and π = 0.01148). In contrast, a lower haplotype (0.426) and nucleotide diversity (θ = 0.00125; π = 0.00095) was observed in the Pfrh5 gene. Neutrality tests do not show deviation from neutral expectations for Pfceltos, with the circulation of multiple low frequency haplotypes (Tajima's D = -0.21637; Fu and Li's D = -0.08164; Fu and Li's F = -0.14051). Strong linkage disequilibrium was observed between variable sites, in each of the genes studied. We postulate that the high diversity and circulation of multiple haplotypes has the potential of making a Pfceltos-subunit vaccine ineffective, while the low genetic diversity of Pfrh5 gene substantiates its evolutionary conservation and potential as a malaria vaccine candidate.

Overexpression of cytochrome P450 and esterase genes involved in permethrin resistance in larvae and adults of Culex quinquefasciatus

Mosquitoes are important vectors of several arthropod-borne diseases, which remain a priority for epidemiological research. Mosquito vector control strategies have traditionally relied on chemical insecticides such as synthetic pyrethroids. However, the indiscriminate use of pesticides has resulted in the development of resistance in many mosquito species. In insects, resistance evolves primarily through the overexpression of one or more gene products from the cytochrome P450, carboxylesterase, and glutathione superfamilies. The current study examined the expression of cytochrome P450 CYP6M2, CYP6AA7, CYP6Z2, CYP9J34, α-Esterase, Esterase B1, and neuroactin genes in larvae and adults of a permethrin-resistant (PerRes) and susceptible (Sus) Culex quinquefasciatus strains. The results showed that the CYP6AA7 gene was overexpressed (10-fold) in larvae and adults with PerRes (p < 0.01) followed by CYPJ34 (9.0-fold) and CYP6Z2 (5.0-fold) compared to the Sus, whereas fewer changes in CYP6M gene expression were observed in PerRes adults (p < 0.05), and no expression was found in larvae. The esterase gene was overexpressed in PerRes larvae (9.0-fold) followed by adults (2.5-fold) compared to the susceptible strain. Based on data, the present study suggests that cytochrome P450, CYP6AA7, CYP6Z2, CYP9J34, α-Esterase, Esterase B1, and neuroactin genes were involved in permethrin resistance in larval and adult Cx. quinquefasciatus.

Detection of insecticide resistance markers in Anopheles funestus from the Democratic Republic of the Congo using a targeted amplicon sequencing panel

Vector control strategies have been successful in reducing the number of malaria cases and deaths globally, but the spread of insecticide resistance represents a significant threat to disease control. Insecticide resistance has been reported across Anopheles (An.) vector populations, including species within the An. funestus group. These mosquitoes are responsible for intense malaria transmission across sub-Saharan Africa, including in the Democratic Republic of the Congo (DRC), a country contributing > 12% of global malaria infections and mortality events. To support the continuous efficacy of vector control strategies, it is essential to monitor insecticide resistance using molecular surveillance tools. In this study, we developed an amplicon sequencing ("Amp-seq") approach targeting An. funestus, and using multiplex PCR, dual index barcoding, and next-generation sequencing for high throughput and low-cost applications. Using our Amp-seq approach, we screened 80 An. funestus field isolates from the DRC across a panel of nine genes with mutations linked to insecticide resistance (ace-1, CYP6P4, CYP6P9a, GSTe2, vgsc, and rdl) and mosquito speciation (cox-1, mtND5, and ITS2). Amongst the 18 non-synonymous mutations detected, was N485I, in the ace-1 gene associated with carbamate resistance. Overall, our panel represents an extendable and much-needed method for the molecular surveillance of insecticide resistance in An. funestus populations.

Data-driven networking of global transcriptomics and male sexual development in the main malaria vector, Anopheles funestus

Deaths from malaria remain staggering despite global support that drives research into new territories. One major gap is our understanding of the sexual biological aspects of the male mosquito, which maintain the vector population solidity. Although Anopheles funestus s.s. is an extremely efficient African vector, little is known about the network between its sexual physiology and gene expression. The Culicidae male's sexual maturity involves a suite of physiological changes, such as genitalia rotation that is necessary for successful mating to occur. We show that mating success is guided by genes and physiological plasticity. Transcriptome analysis between newly emerged males (immature) versus males with rotating genitalia (maturing) provides insight into possible molecular mechanisms regulating male sexual behaviour. Putative transcripts that were associated with male sexual maturation were identified and validated. The discovery of the functions of these transcripts could lead to identifying potential targets for innovative vector control interventions, and mosquito population suppression.

Evidence of intensification of pyrethroid resistance in the major malaria vectors in Kinshasa, Democratic Republic of Congo

Assessing patterns and evolution of insecticide resistance in malaria vectors is a prerequisite to design suitable control strategies. Here, we characterised resistance profile in Anopheles gambiae and Anopheles funestus in Kinshasa and assess the level of aggravation by comparing to previous 2015 estimates. Both species collected in July 2021 were highly resistant to pyrethroids at 1×, 5× and 10× concentrations (mortality < 90%) and remain fully susceptible to bendiocarb and pirimiphos methyl. Compared to 2015, Partial recovery of susceptibility was observed in A. gambiae after PBO synergist assays for both permethrin and α-cypermethrin and total recovery of susceptibility was observed for deltamethrin in 2021. In addition, the efficacy of most bednets decreased significantly in 2021. Genotyping of resistance markers revealed a near fixation of the L1014-Kdr mutation (98.3%) in A. gambiae in 2021. The frequency of the 119F-GSTe2 resistant significantly increased between 2015 and 2021 (19.6% vs 33.3%; P = 0.02) in A. funestus. Transcriptomic analysis also revealed a significant increased expression (P < 0.001) of key cytochrome P450s in A. funestus notably CYP6P9a. The escalation of pyrethroid resistance observed in Anopheles populations from Kinshasa coupled with increased frequency/expression level of resistance genes highlights an urgent need to implement tools to improve malaria vector control.

High efficacy of chlorfenapyr-based net Interceptor® G2 against pyrethroid-resistant malaria vectors from Cameroon

BACKGROUND

The increasing reports of resistance to pyrethroid insecticides associated with reduced efficacy of pyrethroid-only interventions highlight the urgency of introducing new non-pyrethroid-only control tools. Here, we investigated the performance of piperonyl-butoxide (PBO)-pyrethroid [Permanet 3.0 (P3.0)] and dual active ingredients (AI) nets [Interceptor G2 (IG2): containing pyrethroids and chlorfenapyr and Royal Guard (RG): containing pyrethroids and pyriproxyfen] compared to pyrethroid-only net Royal Sentry (RS) against pyrethroid-resistant malaria vectors in Cameroon.

METHODS

The efficacy of these tools was firstly evaluated on Anopheles gambiae s.l. and Anopheles funestus s.l. from Gounougou, Mibellon, Mangoum, Nkolondom, and Elende using cone/tunnel assays. In addition, experimental hut trials (EHT) were performed to evaluate the performance of unwashed and 20 times washed nets in semi-field conditions. Furthermore, pyrethroid-resistant markers were genotyped in dead vs alive, blood-fed vs unfed mosquitoes after exposure to the nets to evaluate the impact of these markers on net performance. The XLSTAT software was used to calculate the various entomological outcomes and the Chi-square test was used to compare the efficacy of various nets. The odds ratio and Fisher exact test were then used to establish the statistical significance of any association between insecticide resistance markers and bed net efficacy.

RESULTS

Interceptor G2 was the most effective net against wild pyrethroid-resistant An. funestus followed by Permanet 3.0. In EHT, this net induced up to 87.8% mortality [95% confidence interval (CI): 83.5-92.1%) and 55.6% (95% CI: 48.5-62.7%) after 20 washes whilst unwashed pyrethroid-only net (Royal Sentry) killed just 18.2% (95% CI: 13.4-22.9%) of host-seeking An. funestus. The unwashed Permanet 3.0 killed up to 53.8% (95% CI: 44.3-63.4%) of field-resistant mosquitoes and 47.2% (95% CI: 37.7-56.7%) when washed 20 times, and the Royal Guard 13.2% (95% CI: 9.0-17.3%) for unwashed net and 8.5% (95% CI: 5.7-11.4%) for the 20 washed net. Interceptor G2, Permanet 3.0, and Royal Guard provided better personal protection (blood-feeding inhibition 66.2%, 77.8%, and 92.8%, respectively) compared to pyrethroid-only net Royal Sentry (8.4%). Interestingly, a negative association was found between kdrw and the chlorfenapyr-based net Interceptor G2 (χ2 = 138; P < 0.0001) with homozygote-resistant mosquitoes predominantly found in the dead ones.

CONCLUSIONS

The high mortality recorded with Interceptor G2 against pyrethroid-resistant malaria vectors in this study provides first semi-field evidence of high efficacy against these major malaria vectors in Cameroon encouraging the implementation of this novel net for malaria control in the country. However, the performance of this net should be established in other locations and on other major malaria vectors before implementation at a large scale.

Heritability of biting time behaviours in the major African malaria vector Anopheles arabiensis

BACKGROUND

The use of insecticide-treated nets for malaria control has been associated with shifts in mosquito vector feeding behaviour including earlier and outdoor biting on humans. The relative contribution of phenotypic plasticity and heritability to these behavioural shifts is unknown. Elucidation of the mechanisms behind these shifts is crucial for anticipating impacts on vector control.

METHODS

A novel portable semi-field system (PSFS) was used to experimentally measure heritability of biting time in the malaria vector Anopheles arabiensis in Tanzania. Wild An. arabiensis from hourly collections using the human landing catch (HLC) method were grouped into one of 3 categories based on their time of capture: early (18:00-21:00), mid (22:00-04:00), and late (05:00-07:00) biting, and placed in separate holding cages. Mosquitoes were then provided with a blood meal for egg production and formation of first filial generation (F1). The F1 generation of each biting time phenotype category was reared separately, and blood fed at the same time as their mothers were captured host-seeking. The resultant eggs were used to generate the F2 generation for use in heritability assays. Heritability was assessed by releasing F2 An. arabiensis into the PSFS, recording their biting time during a human landing catch and comparing it to that of their F0 grandmothers.

RESULTS

In PSFS assays, the biting time of F2 offspring (early: 18:00-21:00, mid: 22:00-04:00 or late: 05:00-07:00) was significantly positively associated with that of their wild-caught F0 grandmothers, corresponding to an estimated heritability of 0.110 (95% CI 0.003, 0.208). F2 from early-biting F0 were more likely to bite early than F2 from mid or late-biting F0. Similarly, the probability of biting late was higher in F2 derived from mid and late-biting F0 than from early-biting F0.

CONCLUSIONS

Despite modest heritability, our results suggest that some of the variation in biting time is attributable to additive genetic variation. Selection can, therefore, act efficiently on mosquito biting times, highlighting the need for control methods that target early and outdoor biting mosquitoes.

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

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

Species composition, infection rate and detection of resistant alleles in Anopheles funestus (Diptera: Culicidae) from Lare, a malaria hotspot district of Ethiopia

BACKGROUND

Anopheles funestus, which is considered as secondary vector of malaria in Ethiopia, is known to have several morphologically indistinguishable (sibling) species. Accurate identification of sibling species is crucial to understand their biology, behaviour and vector competence. In this study, molecular identification was conducted on the Ethiopian An. funestus populations. Moreover, insecticide resistance mechanism markers were detected, including ace N485I, kdr L1014F, L1014S, and CYP6P9a TaqMan qPCR was used to detect the infective stage of the parasite from field collected adult female An. funestus populations.

METHODS

Adult female mosquito collection was conducted from Lare, Gambella Regional State of Ethiopia between June 2018 to July 2020 using CDC light traps and HLC. Sub-samples of the morphologically identified An. funestus mosquitoes were molecularly identified using species-specific PCR, and the possible presence of insecticide resistance alleles was investigated using TaqMan qPCR (N485I-Ace-1), PCR-Sanger sequencing (L1014F-kdr), and PCR-RFLP (CYP6P9a resistance allele). Following head/thorax dissection, the TaqMan qPCR assay was used to investigate the presence of the infective stage Plasmodium parasite species.

RESULTS

A total of 1086 adult female An. funestus mosquitoes were collected during the study period. All sub-samples (N = 20) that were morphologically identified as An. funestus sensu lato (s.l.) were identified as An. funestus sensu stricto (s.s.) using species- specific PCR assay. The PCR-RFLP assay that detects the CYP6P9a resistance allele that confers pyrethroid resistance in An. funestus was applied in N = 30 randomly selected An. funestus s.l.

SPECIMENS

None of the specimens showed a digestion pattern consistent with the presence of the CYP6P9a resistance allele in contrast to what was observed in the positive control. Consequently, all samples were characterized as wild type. The qPCR TaqMan assay that detects the N485I acetylcholinesterase-1 mutation conferring resistance to organophosphates/carbamates in An. funestus was used in (N = 144) samples. All samples were characterized as wild type. The kdr L1014F and L1014S mutations in the VGSC gene that confer resistance to pyrethroids and DDT were analysed with direct Sanger sequencing after PCR and clean-up of the PCR products were also characterized as wild type. None of the samples (N = 169) were found positive for Plasmodium (P. falciparum/ovale/malariae/vivax) detection.

CONCLUSION

All An. funestus s.l. samples from Lare were molecularly identified as An. funestus s.s. No CYP6P9, N485I acetylcholinesterase 1, kdr L1014F or L1014S mutations were detected in the An. funestus samples. None of the An. funestus samples were positive for Plasmodium. Although the current study did not detect any insecticide resistant mechanism, it provides a reference for future vector monitoring programmes. Regular monitoring of resistance mechanisms covering wider geographical areas of Ethiopia where this vector is distributed is important for improving the efficacy of vector control programs.

Advances in the genetic characterization of the malaria vector, Anopheles funestus, and implications for improved surveillance and control

Anopheles mosquitoes present a major public health challenge in sub-Saharan Africa; notably, as vectors of malaria that kill over half a million people annually. In parts of the east and southern Africa region, one species in the Funestus group, Anopheles funestus, has established itself as an exceptionally dominant vector in some areas, it is responsible for more than 90% of all malaria transmission events. However, compared to other malaria vectors, the species is far less studied, partly due to difficulties in laboratory colonization and the unresolved aspects of its taxonomy and systematics. Control of An. funestus is also increasingly difficult because it has developed widespread resistance to public health insecticides. Fortunately, recent advances in molecular techniques are enabling greater insights into species identity, gene flow patterns, population structure, and the spread of resistance in mosquitoes. These advances and their potential applications are reviewed with a focus on four research themes relevant to the biology and control of An. funestus in Africa, namely: (i) the taxonomic characterization of different vector species within the Funestus group and their role in malaria transmission; (ii) insecticide resistance profile; (iii) population genetic diversity and gene flow, and (iv) applications of genetic technologies for surveillance and control. The research gaps and opportunities identified in this review will provide a basis for improving the surveillance and control of An. funestus and malaria transmission in Africa.

Effects of next-generation, dual-active-ingredient, long-lasting insecticidal net deployment on insecticide resistance in malaria vectors in Tanzania: an analysis of a 3-year, cluster-randomised controlled trial

BACKGROUND

Insecticide resistance among malaria-vector species is a pervasive problem that might jeopardise global disease-control efforts. Novel vector-control tools with different modes of action, including long-lasting insecticidal nets (LLINs) incorporating new active ingredients, are urgently needed to delay the evolution and spread of insecticide resistance. We aimed to measure phenotypic and genotypic insecticide-resistance profiles among wild Anopheles collected over 3 years to assess the longitudinal effects of dual-active-ingredient LLINs on insecticide resistance.

METHODS

For this analysis, data nested in a 3-year, four parallel-arm, superiority cluster-randomised controlled trial (cRCT) in Tanzania, collected from 84 clusters (39 307 households) formed of 72 villages in the Misungwi district, were used to measure insecticide-resistance profiles among female Anopheles mosquitoes via insecticide-resistance bioassays and quantitative RT-PCR of metabolic-resistance genes. Wild, blood-fed, indoor-resting mosquitoes were collected annually during the rainy seasons from house walls in clusters from all four trial groups. Mosquitoes were morphologically identified as An gambiae sensu lato (SL) or An funestus SL before separate bioassay testing. The primary outcomes were lethal-dose values for α-cypermethrin, permethrin, and piperonyl butoxide pre-exposure plus permethrin-resistance intensity bioassays, mortality 72 h after insecticidal exposure for chlorfenapyr bioassays, fertility reduction 72 h after insecticidal exposure for pyriproxyfen bioassays, and fold change in metabolic-enzyme expression relative to an insecticide-susceptible laboratory strain. All primary outcomes were measured in An funestus SL 1 year, 2 years, and 3 years after LLIN distribution. Primary outcomes were also assessed in An gambiae SL if enough mosquitoes were collected. The cRCT is registered with ClinicalTrials.gov (NCT03554616).

FINDINGS

Between May 24, 2019, and Oct 25, 2021, 47 224 female Anopheles were collected for resistance monitoring. In the pyrethroid (PY)-LLIN group, there were significant increases in α-cypermethrin-resistance intensity (year 1 LD50=9·52 vs year 2 76·20, p<0·0001) and permethrin-resistance intensity (year 1 13·27 vs year 2 35·83, p=0·0019) in An funestus SL. In the pyriproxyfen PY-LLIN group, there was similar increase in α-cypermethrin-resistance intensity (year 1 0·71 vs year 2 81·56, p<0·0001) and permethrin-resistance intensity (year 1 5·68 vs year 2 50·14, p<0·0001). In the piperonyl butoxide PY-LLIN group, α-cypermethrin-resistance intensity (year 1 33·26 vs year 3 70·22, p=0·0071) and permethrin-resistance intensity (year 1 47·09 vs year 3 2635·29, p<0·0001) also increased over time. In the chlorfenapyr PY-LLIN group, there were no effects on α-cypermethrin-resistance intensity (year 1 0·42 vs year 3 0·99, p=0·54) or permethrin-resistance intensity (data were not estimable due to nearly 100% mortality). There were also minimal reductions in chlorfenapyr susceptibility. However, in the chlorfenapyr PY-LLIN group, a significant decline in piperonyl-butoxide synergy was seen by year 3 (year 1 0·02 vs year 3 0·26, p=0·020). Highly over-expressed detoxification enzymes showed dynamic patterns of selection throughout the trial.

INTERPRETATION

Our phenotypic data supports trial epidemiological findings; chlorfenapyr PY-LLINs provided superior protection from malaria across multiple transmission seasons, with few effects on insecticide-resistance selection. Rapid pyrethroid-resistance intensification in the piperonyl butoxide PY-LLIN group and pre-existing tolerance of pyriproxyfen in vector populations might explain the poorer performance of these two interventions regarding malaria outcomes. Further work is required to elucidate the potential mechanisms driving cross-resistance between pyrethroids and novel active ingredients to better inform the design of pre-emptive resistance-management strategies.

FUNDING

UK Department for International Development; UK Medical Research Council; Wellcome Trust; UK Department of Health and Social Care; UK Foreign, Commonwealth and Development Office; and The Bill and Melinda Gates Foundation via the Innovative Vector Control Consortium.

Molecular drivers of insecticide resistance in the Sahelo-Sudanian populations of a major malaria vector Anopheles coluzzii

BACKGROUND

Information on common markers of metabolic resistance in malaria vectors from countries sharing similar eco-climatic characteristics can facilitate coordination of malaria control. Here, we characterized populations of the major malaria vector Anopheles coluzzii from Sahel region, spanning four sub-Saharan African countries: Nigeria, Niger, Chad and Cameroon.

RESULTS

Genome-wide transcriptional analysis identified major genes previously implicated in pyrethroid and/or cross-resistance to other insecticides, overexpressed across the Sahel, including CYP450s, glutathione S-transferases, carboxylesterases and cuticular proteins. Several, well-known markers of insecticide resistance were found in high frequencies-including in the voltage-gated sodium channel (V402L, I940T, L995F, I1527T and N1570Y), the acetylcholinesterase-1 gene (G280S) and the CYP4J5-L43F (which is fixed). High frequencies of the epidemiologically important chromosomal inversion polymorphisms, 2La, 2Rb and 2Rc, were observed (~80% for 2Rb and 2Rc). The 2La alternative arrangement is fixed across the Sahel. Low frequencies of these inversions (<10%) were observed in the fully insecticide susceptible laboratory colony of An. coluzzii (Ngoussou). Several of the most commonly overexpressed metabolic resistance genes sit in these three inversions. Two commonly overexpressed genes, GSTe2 and CYP6Z2, were functionally validated. Transgenic Drosophila melanogaster flies expressing GSTe2 exhibited extremely high DDT and permethrin resistance (mortalities <10% in 24h). Serial deletion of the 5' intergenic region, to identify putative nucleotide(s) associated with GSTe2 overexpression, revealed that simultaneous insertion of adenine nucleotide and a transition (T->C), between Forkhead box L1 and c-EST putative binding sites, were responsible for the high overexpression of GSTe2 in the resistant mosquitoes. Transgenic flies expressing CYP6Z2 exhibited marginal resistance towards 3-phenoxybenzylalcohol (a primary product of pyrethroid hydrolysis by carboxylesterases) and a type II pyrethroid, α-cypermethrin. However, significantly higher mortalities were observed in CYP6Z2 transgenic flies compared with controls, on exposure to the neonicotinoid, clothianidin. This suggests a possible bioactivation of clothianidin into a toxic intermediate, which may make it an ideal insecticide against populations of An. coluzzii overexpressing this P450.

CONCLUSIONS

These findings will facilitate regional collaborations within the Sahel region and refine implementation strategies through re-focusing interventions, improving evidence-based, cross-border policies towards local and regional malaria pre-elimination.

Limited Effect of Ground Floor Fogging on Mosquito Distribution in High-Rise Condominia

Fogging with insecticides is one of the main control measures for adult mosquito populations employed in countries that are affected by dengue. In many such countries, urban communities are increasingly characterised by high-density residence in high-rise condominia. Although fogging is typically applied at the ground level, its efficacy in three-dimensional urban environments is poorly understood. Here, we investigated the effect of fogging on vector mosquito distribution and abundance in high-rise condominia by conducting a before-after fogging survey. We showed that although mosquitoes were significantly concentrated at the lower levels in high-rise condominia, they were found throughout the three-dimensional environments. Fogging did not significantly alter this distribution or abundance pattern across any floor level. Thus, any fogging effect was short-lived as mosquito populations recovered within a few days before the subsequent scheduled treatment. In addition, increasing fogging frequency within practicable limits did not prolong the intended control effect. As urban mosquitoes are increasingly insusceptible to fogging due to insecticide resistance and vertical avoidance, this study demonstrates the need to implement other mosquito control strategies for high-rise condominia to manage mosquito populations.

A major 6 Mb superlocus is involved in pyrethroid resistance in the common bed bug Cimex lectularius

In the last few years, the bed bug Cimex lectularius has been an increasing problem worldwide, mainly due to the development of insecticide resistance to pyrethroids. The characterization of resistance alleles is a prerequisite to improve surveillance and resistance management. To identify genomic variants associated with pyrethroid resistance in Cimex lectularius, we compared the genetic composition of two recent and resistant populations with that of two ancient-susceptible strains using a genome-wide pool-seq design. We identified a large 6 Mb "superlocus" showing particularly high genetic differentiation and association with the resistance phenotype. This superlocus contained several clustered resistance genes and was also characterized by a high density of structural variants (inversions, duplications). The possibility that this superlocus constitutes a resistance "supergene" that evolved after the clustering of alleles adapted to insecticide and after reduction in recombination is discussed.

Resistance to clothianidin reduces the efficacy of SumiShield® 50WG, a neonicotinoid formulation for indoor residual spraying, against Anopheles gambiae

Chronic exposure of mosquito larvae to pesticide residues in agricultural areas is often associated with evolution of resistance to insecticides used for vector control. This presents a concern for the efficacy of clothianidin, an agricultural neonicotinoid qualified for Indoor Residual Spraying (IRS). Using standard bioassays, we tested if reduced susceptibility to clothianidin affects the efficacy of SumiShield® 50WG, one of the two newly approved formulations, which contains 50% clothianidin. We simultaneously monitored susceptibility to clothianidin and to SumiShield® 50WG, testing adults of Anopheles gambiae, An. coluzzii and Culex sp collected from urban, suburban and agricultural areas of Yaoundé. We found that the level of susceptibility to the active ingredient predicted the efficacy of SumiShield® 50WG. This formulation was very potent against populations that achieved 100% mortality within 72 h of exposure to a discriminating dose of clothianidin. By contrast, mortality leveled off at 75.4 ± 3.5% within 7 days of exposure to SumiShield® 50WG in An. gambiae adults collected from a farm where spraying of acetamiprid and imidacloprid is driving cross-resistance to clothianidin. These findings indicate that more potent formulations of clothianidin or different insecticides should be prioritized in areas where resistance is emerging.

Fine-scale spatial distribution of deltamethrin resistance and population structure of Anopheles funestus and Anopheles arabiensis populations in Southern Mozambique

BACKGROUND

Insecticide resistance in malaria vectors can be spatially highly heterogeneous, yet population structure analyses frequently find relatively high levels of gene flow among mosquito populations. Few studies have contemporaneously assessed phenotypic, genotypic and population structure analysis on mosquito populations and none at fine geographical scales. In this study, genetic diversity, population structure, and insecticide resistance profiles of Anopheles funestus and Anopheles arabiensis were examined across mosquito populations from and within neighbouring villages.

METHODS

Mosquitoes were collected from 11 towns in southern Mozambique, as well as from different neighbourhoods within the town of Palmeira, during the peak malaria transmission season in 2016. CDC bottle bioassay and PCR assays were performed with Anopheles mosquitoes at each site to determine phenotypic and molecular insecticide resistance profiles, respectively. Microsatellite analysis was conducted on a subsample of mosquitoes to estimate genetic diversity and population structure.

RESULTS

Phenotypic insecticide resistance to deltamethrin was observed in An. funestus sensu stricto (s.s.) throughout the area, though a high level of mortality variation was seen. However, 98% of An. funestus s.s. were CYP6P9a homozygous resistant. An. arabiensis was phenotypically susceptible to deltamethrin and 99% were kdr homozygous susceptible. Both Anopheles species exhibited high allelic richness and heterozygosity. Significant deviations from Hardy-Weinberg equilibrium were observed, and high linkage disequilibrium was seen for An. funestus s.s., supporting population subdivision. However, the FST values were low for both anophelines (- 0.00457 to 0.04213), Nm values were high (9.4-71.8 migrants per generation), AMOVA results showed almost 100% genetic variation among and within individuals, and Structure analysis showed no clustering of An. funestus s.s. and An. arabiensis populations. These results suggest high gene flow among mosquito populations.

CONCLUSION

Despite a relatively high level of phenotypic variation in the An. funestus population, molecular analysis shows the population is admixed. These data indicate that CYP6P9a resistance markers do not capture all phenotypic variation in the area, but also that resistance genes of high impact are likely to easily spread in the area. Conversely, other strategies, such as transgenic mosquito release programmes will likely not face challenges in this locality.

Contrasting Patterns of Asaia Association with Pyrethroid Resistance Escalation between the Malaria Vectors Anopheles funestus and Anopheles gambiae

Microbiome composition has been associated with insecticide resistance in malaria vectors. However, the contribution of major symbionts to the increasingly reported resistance escalation remains unclear. This study explores the possible association of a specific endosymbiont, Asaia spp., with elevated levels of pyrethroid resistance driven by cytochrome P450s enzymes and voltage-gated sodium channel mutations in Anopheles funestus and Anopheles gambiae. Molecular assays were used to detect the symbiont and resistance markers (CYP6P9a/b, 6.5 kb, L1014F, and N1575Y). Overall, genotyping of key mutations revealed an association with the resistance phenotype. The prevalence of Asaia spp. in the FUMOZ_X_FANG strain was associated with the resistance phenotype at a 5X dose of deltamethrin (OR = 25.7; p = 0.002). Mosquitoes with the resistant allele for the markers tested were significantly more infected with Asaia compared to those possessing the susceptible allele. Furthermore, the abundance correlated with the resistance phenotype at 1X concentration of deltamethrin (p = 0.02, Mann-Whitney test). However, for the MANGOUM_X_KISUMU strain, findings rather revealed an association between Asaia load and the susceptible phenotype (p = 0.04, Mann-Whitney test), demonstrating a negative link between the symbiont and permethrin resistance. These bacteria should be further investigated to establish its interactions with other resistance mechanisms and cross-resistance with other insecticide classes.

The duplicated P450s CYP6P9a/b drive carbamates and pyrethroids cross-resistance in the major African malaria vector Anopheles funestus

Cross-resistance to insecticides in multiple resistant malaria vectors is hampering resistance management. Understanding its underlying molecular basis is critical to implementation of suitable insecticide-based interventions. Here, we established that the tandemly duplicated cytochrome P450s, CYP6P9a/b are driving carbamate and pyrethroid cross-resistance in Southern African populations of the major malaria vector Anopheles funestus. Transcriptome sequencing revealed that cytochrome P450s are the most over-expressed genes in bendiocarb and permethrin-resistant An. funestus. The CYP6P9a and CYP6P9b genes are overexpressed in resistant An. funestus from Southern Africa (Malawi) versus susceptible An. funestus (Fold change (FC) is 53.4 and 17 respectively), while the CYP6P4a and CYP6P4b genes are overexpressed in resistant An. funestus in Ghana, West Africa, (FC is 41.1 and 17.2 respectively). Other up-regulated genes in resistant An. funestus include several additional cytochrome P450s (e.g. CYP9J5, CYP6P2, CYP6P5), glutathione-S transferases, ATP-binding cassette transporters, digestive enzymes, microRNA and transcription factors (FC<7). Targeted enrichment sequencing strongly linked a known major pyrethroid resistance locus (rp1) to carbamate resistance centering around CYP6P9a/b. In bendiocarb resistant An. funestus, this locus exhibits a reduced nucleotide diversity, significant p-values when comparing allele frequencies, and the most non-synonymous substitutions. Recombinant enzyme metabolism assays showed that both CYP6P9a/b metabolize carbamates. Transgenic expression of CYP6P9a/b in Drosophila melanogaster revealed that flies expressing both genes were significantly more resistant to carbamates than controls. Furthermore, a strong correlation was observed between carbamate resistance and CYP6P9a genotypes with homozygote resistant An. funestus (CYP6P9a and the 6.5kb enhancer structural variant) exhibiting a greater ability to withstand bendiocarb/propoxur exposure than homozygote CYP6P9a_susceptible (e.g Odds ratio = 20.8, P<0.0001 for bendiocarb) and heterozygotes (OR = 9.7, P<0.0001). Double homozygote resistant genotype (RR/RR) were even more able to survive than any other genotype combination showing an additive effect. This study highlights the risk that pyrethroid resistance escalation poses to the efficacy of other classes of insecticides. Available metabolic resistance DNA-based diagnostic assays should be used by control programs to monitor cross-resistance between insecticides before implementing new interventions.

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

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

Multi-centre discriminating concentration determination of broflanilide and potential for cross-resistance to other public health insecticides in Anopheles vector populations

Novel insecticides are urgently needed to control insecticide-resistant populations of Anopheles malaria vectors. Broflanilide acts as a non-competitive antagonist of the gamma-aminobutyric acid receptor and has shown prolonged effectiveness as an indoor residual spraying product (VECTRON T500) in experimental hut trials against pyrethroid-resistant vector populations. This multi-centre study expanded upon initial discriminating concentration testing of broflanilide, using six Anopheles insectary colonies (An. gambiae Kisumu KCMUCo, An. gambiae Kisumu NIMR, An. arabiensis KGB, An. arabiensis SENN, An. coluzzii N'Gousso and An. stephensi SK), representing major malaria vector species, to facilitate prospective susceptibility monitoring of this new insecticide; and investigated the potential for cross-resistance to broflanilide via the A296S mutation associated with dieldrin resistance (rdl). Across all vector species tested, the discriminating concentration for broflanilide ranged between LC₉₉ × 2 = 1.126-54.00 μg/ml or LC₉₅ × 3 = 0.7437-17.82 μg/ml. Lower concentrations of broflanilide were required to induce complete mortality of An. arabiensis SENN (dieldrin-resistant), compared to its susceptible counterpart, An. arabiensis KGB, and there was no association between the presence of the rdl mechanism of resistance and survival in broflanilide bioassays, demonstrating a lack of cross-resistance to broflanilide. Study findings provide a benchmark for broflanilide susceptibility monitoring as part of ongoing VECTRON T500 community trials in Tanzania and Benin.

Genetic Diversity of Cytochrome P450s CYP6M2 and CYP6P4 Associated with Pyrethroid Resistance in the Major Malaria Vectors Anopheles coluzzii and Anopheles gambiae from Yaoundé, Cameroon

Assessing the genetic diversity of metabolic resistance genes, such as cytochrome P450s, helps to understand the dynamics and evolution of resistance in the field. Here, we analyzed the polymorphisms of CYP6M2 and CYP6P4, associated with pyrethroid resistance in Anopheles coluzzii and Anopheles gambiae, to detect potential resistance markers. Field-caught resistant mosquitos and susceptible lab strains were crossed, and F4 was exposed to permethrin for 15 min and 90 min to discriminate highly susceptible (HS) and highly resistant (HR) mosquitos, respectively. Significant permethrin mortality reduction was observed after pre-exposure to PBO, suggesting the gene involvement of P450s. qPCR analysis revealed significant overexpression of CYP6M2 (FC = 19.57 [95% CI 13.96-25.18] for An. coluzzii; 10.16 [7.86-12.46] for An. gambiae) and CYP6P4 (FC = 6.73 [6.15-7.30] An. coluzzii; 23.62 [26.48-20.76] An. gambiae). Full-gene and ≈1 kb upstream were sequenced. For CYP6M2, the upstream region shows low diversity in HR and HS (overall Hd = 0.49, π = 0.018), whereas the full-gene shows allelic-variation but without evidence of ongoing selection. CYP6P4 upstream region showed a lower diversity in HR (Hd = 0.48) than HS (Hd = 0.86) of An. gambiae. These results highlighted that CYP6P4-associated resistance is potentially driven by modification in upstream region. However, further work is needed to determine the real causative variants that will help design rapid detection tools.

Marked Effects of Larval Salt Exposure on the Life History and Gut Microbiota of the Malaria Vector Anopheles merus (Diptera: Culicidae)

Anopheles merus can breed in a range of saltwater concentrations. The consequences of this ability on the life history of adult An. merus are poorly understood. This study examined the effects of exposure to 0, 2.1875, 4.375, 8.75, and 17.5 g/L of sodium chloride on An. merus. The effects on larval development, adult longevity, fertility, and fecundity, as well as deltamethrin tolerance were examined. The effect of larval salt exposure on the expression of defensin-1 in adults was examined by quantitative Real-Time PCR. Finally, the effect of the larval salt concentration on microbial dynamics was assessed by 16S Next Generation Sequencing. High concentrations of saltwater increased larval development time and number of eggs laid, as well as deltamethrin tolerance. Larval exposure to salt also reduced the expression of defensin-1. The exposure also had a significant effect on microbial diversity in larvae and adults. The diversity of larvae decreased once adults emerged. Salt-tolerant bacterial genera predominated in larvae but were absent in adults. High salt concentrations resulted in greater abundance of Plasmodium-protective genera in adults. Although this study was conducted on a laboratory strain of An. merus, these data suggest that osmoregulation has a significant effect on the life history of the species with potential epidemiological consequences.

Population genetic analysis of Plasmodium falciparum cell-traversal protein for ookinetes and sporozoite among malaria patients from southern Nigeria

Plasmodium falciparum immune escape mechanisms affect antigens being prioritized for vaccine design. As a result of the multiple surface antigens the parasite exhibits at different life cycle stages, designing a vaccine that would efficiently boost the immune system in clearing infections has been challenging. The P. falciparum cell-traversal protein for ookinetes and sporozoite (Pfceltos) is instrumental for ookinete traversal of the mosquito midgut and sporozoites invasion of the human liver cells. Pfceltos elicits both humoral and cellular immune response but has been reported with multiple single nucleotide polymorphisms in global isolates. A cross-sectional survey, conducted in southern Nigeria, between January-March 2021 recruited 283 individuals. Of this, 166 demonstrated P. falciparum infections (86 from Cross River and 80 from Edo), 48 (55.8%) while only 36 (45%) were amplified for Pfceltos gene from both sites respectively. Fifty amplified samples were sequenced and analysed for their diversity, polymorphisms and population structure of the gene. The number of segregating sites in Edo State was higher (34) than that of Cross River State. Though nucleotide diversity was higher for Edo compared to Cross River State (θw = 0.02505; π = 0.03993 versus θw = 0.00930; π = 0.01033 respectively), the reverse was the case for haplotype diversity (0.757 versus 0.890 for Edo and Cross River respectively). Of the twelve haplotypes observed from both states, only two (KASLPVEK and NAFLSFEK) were shared, with haplotype prevalence higher in Edo (16% and 36%) than Cross River (8% and 4%). The Tajima's D test was positive for both states, with Fst value showing a strong genetic differentiation (Fst = 0.25599), indicating the occurrence of balancing selection favoring haplotype circulation at a low frequency. The shared haplotypes, low Hst and Fst values presents a challenge to predict the extent of gene flow. High LD values present a grim public health consequence should a Pfceltos-conjugated vaccine be considered for prophylaxis in Nigeria.

Escalating pyrethroid resistance in two major malaria vectors Anopheles funestus and Anopheles gambiae (s.l.) in Atatam, Southern Ghana

BACKGROUND

Aggravation of insecticide resistance in malaria vectors is threatening the efforts to control malaria by reducing the efficacy of insecticide-based interventions hence needs to be closely monitored. This study investigated the intensity of insecticide resistance of two major malaria vectors An. funestus sensu stricto (s.s.) and An. gambiae sensu lato (s.l.) collected in southern Ghana and assessed the bio-efficacy of several long-lasting insecticidal nets (LLINs) against these mosquito populations.

METHODS

The insecticide susceptibility profiles of Anopheles funestus s.s. and Anopheles gambiae s.l. populations from Obuasi region (Atatam), southern Ghana were characterized and the bio-efficacy of some LLINs was assessed to determine the impact of insecticide resistance on the effectiveness of these tools. Furthermore, molecular markers associated with insecticide resistance in both species were characterized in the F0 and F1 populations using PCR and qPCR methods.

RESULTS

Anopheles funestus s.s. was the predominant species and was resistant to pyrethroids, organochlorine and carbamate insecticides, but fully susceptible to organophosphates. An. gambiae s.l. was resistant to all four insecticide classes. High intensity of resistance to 5 × and 10 × the discriminating concentration (DC) of pyrethroids was observed in both species inducing a considerable loss of efficacy of long-lasting insecticidal nets (LLINs). Temporal expression analysis revealed a massive 12-fold increase in expression of the CYP6P4a cytochrome P450 gene in An. funestus s.s., initially from a fold change of 41 (2014) to 500 (2021). For both species, the expression of candidate genes did not vary according to discriminating doses. An. gambiae s.l. exhibited high frequencies of target-site resistance including Vgsc-1014F (90%) and Ace-1 (50%) while these mutations were absent in An. funestus s.s.

CONCLUSIONS

The multiple and high intensity of resistance observed in both malaria vectors highlights the need to implement resistance management strategies and the introduction of new insecticide chemistries.

Characterization of the glutathione S-transferase genes in the sand flies Phlebotomus papatasi and Lutzomyia longipalpis shows expansion of the novel glutathione S-transferase xi (X) class

Leishmaniasis control often relies upon insecticidal control of phlebotomine sandfly vector populations. Such methods are vulnerable to the evolution of insecticide resistance via a range of molecular mechanisms. There is evidence that two major resistance mechanisms, target site insensitivity and metabolic resistance, have evolved in some sandfly populations and further genetic characterization of resistance would be useful to understand and combat it. To facilitate the study of the mechanisms of metabolic resistance, here we improved the annotation and characterized a major detoxification gene family, the glutathione-s-transferases (GST), in the genomes of two sand fly species: Phlebotomus papatasi and Lutzomyia longipalpis. The compositions of the GST gene family differ markedly from those of Aedes and Anopheles mosquitoes. Most strikingly, the xi (X) class of GSTs appears to have expanded in both sand fly genomes. Our results provide a basis for further studies of metabolic resistance mechanisms in these important disease vector species.

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.

Marked aggravation of pyrethroid resistance in major malaria vectors in Malawi between 2014 and 2021 is partly linked with increased expression of P450 alleles

BACKGROUND

Increased intensity of pyrethroid resistance is threatening the effectiveness of insecticide-based interventions to control malaria in Africa. Assessing the extent of this aggravation and its impact on the efficacy of these tools is vital to ensure the continued control of major vectors. Here we took advantage of 2009 and 2014 data from Malawi to establish the extent of the resistance escalation in 2021 and assessed its impact on various bed nets performance.

METHODS

Indoor blood-fed and wild female Anopheles (An) mosquitoes were collected with an electric aspirator in Chikwawa. Cocktail and SINE PCR were used to identify sibling species belonging to An. funestus group and An. gambiae complex. The susceptibility profile to the four classes of insecticides was assessed using the WHO tubes bioassays. Data were saved in an Excel file. Analysis was done using Vassarstats and figures by Graph Pad.

RESULTS

In this study, a high level of resistance was observed with pyrethroids (permethrin, deltamethrin and alpha-cypermethrin with mortality rate at 5x discriminating concentration (DC) < 50% and Mortality rate at 10x DC < 70%). A high level of resistance was also observed to carbamate (bendiocarb) with mortality rate at 5x DC < 25%). Aggravation of resistance was also noticed between 2009 and 2021. For pyrethroids, the mortality rate for permethrin reduced from 47.2% in 2009 to 13% in 2014 and 6.7% in 2021. For deltamethrin, the mortality rate reduced from 42.3% in 2009 to 1.75% in 2014 and 5.2% in 2021. For Bendiocarb, the mortality rate reduced from 60% in 2009 to 30.1% in 2014 and 12.2% in 2021. The high resistance observed is consistent with a drastic loss of pyrethroid-only bed nets efficacy although Piperonyl butoxide (PBO)-based nets remain effective. The resistance pattern observed was linked with high up-regulation of the P450 genes CYP6P9a, CYP6P9b and CYP6M7 in An. funestus s.s. mosquitoes surviving exposure to deltamethrin at 1x, 5x and 10x DC. A significant association was observed between the 6.5 kb structural variant and resistance escalation with homozygote resistant (SV+/SV+) more likely to survive exposure to 5x and 10x (OR = 4.1; P < 0.001) deltamethrin than heterozygotes. However, a significant proportion of mosquitoes survived the synergist assays with PBO suggesting that other mechanisms than P450s are present.

CONCLUSIONS

This resistance aggravation in An. funestus s.s. Malawian population highlights an urgent need to deploy novel control tools not relying on pyrethroids to improve the effectiveness of vector control.

Early Transcriptional Response to Monensin in Sensitive and Resistant Strains of Eimeria tenella

Eimeria parasites are the causative agents of coccidiosis, a common parasitic disease in poultry and livestock that causes significant economic losses to the animal husbandry industry. Ionophore coccidiostats, such as monensin and salinomycin, are widely used for prophylaxis of coccidiosis in poultry. Unfortunately, widespread drug resistance has compromised their efficacy. As a result, there is an increasing need to understand the targets and resistance mechanisms to anticoccidials. However, how Eimeria parasite genes respond to ionophores remains unclear. In this study, resistance to monensin was induced in E. tenella through serial generations of selection. Both sensitive and resistant E. tenella sporozoites were treated with 5 μg/ml monensin for 0, 2, and 4 h, respectively. Gene transcription profiles were then compared by high-throughput sequencing. The results showed that protein translation-related genes were significantly downregulated after drug induction. A total of 1,848 DEGs were detected in the sensitive strain after 2 h of exposure, whereas only 31 were detected in the resistant strain. Among these DEGs in the sensitive strain, genes associated with protein degradation were significantly upregulated, supporting the autophagy-like parasite killing theory. Then, 4 h of exposure resulted in additional 626 and 621 DEGs for sensitive and resistant strains, respectively. This result implies that the gene transcription in sensitive strain is more susceptible to monensin treatment. Our results provide gene expression landscapes of E. tenella following monensin treatment. These data will contribute to a better understanding of the mechanism of drug resistance to polyether ionophores in coccidia.

Reduced performance of community bednets against pyrethroid-resistant Anopheles funestus and Anopheles gambiae, major malaria vectors in Cameroon

BACKGROUND

Long-lasting insecticidal nets (LLINs) are a vital tool in the fight against malaria vectors. However, their efficacy in the field can be impacted by several factors, including patterns of usage, net age, mosquito resistance and the delayed mortality effect, all of which could influence malaria transmission. We have investigated the effectiveness of the various brands of LLINs available in markets and households in Cameroon on pyrethroid-resistant mosquitoes and assessed their post-exposure effect.

METHODS

Following quality control assessment on a susceptible laboratory mosquito strain, we evaluated the immediate and delayed mortality effects of exposure to LLINs (both newly bough LLINst and used ones collected from households in Elende village, Cameroon, in 2019) using standard WHO cone tests on Anopheles gambiae and Anopheles funestus populations collected from the Centre region of Cameroon. Alive female mosquitoes were genotyped for various resistance markers at different time points post-exposure to evaluate the impact of insecticide resistance on the efficacy of bednets.

RESULTS

The laboratory-susceptible strain experienced high mortality rates when exposed to all pyrethroid-only brands of purchased nets (Olyset® Net, Super Net, PermaNet® 2.0, Yorkool®, Royal Sentry®) (Mean±SEM: 68.66 ± 8.35% to 93.33 ± 2.90%). However, low mortality was observed among wild An. funestus mosquitoes exposed to the bednets (0 ± 0 to 28 ± 6.7%), indicating a reduced performance of these nets against field mosquitoes. Bednets collected from households also showed reduced efficacy on the laboratory strain (mortality: 19-66%), as well as displaying a significant loss of efficacy against the local wild strains (mortality: 0 ± 0% to 4 ± 2.6% for An. gambiae sensu lato and 0 ± 0% to 8 ± 3.2% for An. funestus). However, compared to the unexposed group, mosquitoes exposed to bednets showed a significantly reduced longevity, indicating that the efficacy of these nets was not completely lost. Mosquitoes with the CYP6P9a-RR and L119F-GSTe2 mutations conferring pyrethroid resistance showed greater longevity after exposure to the Olyset net than their susceptible counterparts, indicating the impact of resistance on bednet efficacy and delayed mortality.

CONCLUSION

These findings show that although standard bednets drastically lose their efficacy against pyrethroid-resistant field mosquitoes, they still are able to induce delayed mortality in exposed populations. The results of this study also provide evidence of the actual impact of resistance on the quality and efficacy of LLINs in use in the community, with mosquitoes carrying the CYP6P9a-RR and L119F-GSTe2 mutations conferring pyrethroid resistance living longer than their susceptible counterparts. These results highlight the need to use new-generation nets that do not rely solely on pyrethroids.

Gene Conversion Explains Elevated Diversity in the Immunity Modulating APL1 Gene of the Malaria Vector Anopheles funestus

Leucine-rich repeat proteins and antimicrobial peptides are the key components of the innate immune response to Plasmodium and other microbial pathogens in Anopheles mosquitoes. The APL1 gene of the malaria vector Anopheles funestus has exceptional levels of non-synonymous polymorphism across the range of An. funestus, with an average πn of 0.027 versus a genome-wide average of 0.002, and πn is consistently high in populations across Africa. Elevated APL1 diversity was consistent between the independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide resistance was observed. Although lacking the diversity of APL1, two further mosquito innate-immunity genes of the gambicin anti-microbial peptide family had πn/πs ratios greater than one, possibly driven by either positive or balancing selection. The cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, a result discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogues, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. In conclusion, we hypothesize that higher gene expression of APL1 than its paralogues is correlated with a more open chromatin formation, which enhances gene conversion and elevated diversity at this locus.

Experimental Hut Trials Reveal That CYP6P9a/b P450 Alleles Are Reducing the Efficacy of Pyrethroid-Only Olyset Net against the Malaria Vector Anopheles funestus but PBO-Based Olyset Plus Net Remains Effective

Malaria remains a major public health concern in Africa. Metabolic resistance in major malaria vectors such as An. funestus is jeopardizing the effectiveness of long-lasting insecticidal nets (LLINs) to control malaria. Here, we used experimental hut trials (EHTs) to investigate the impact of cytochrome P450-based resistance on the efficacy of PBO-based net (Olyset Plus) compared to a permethrin-only net (Olyset), revealing a greater loss of efficacy for the latter. EHT performed with progenies of F5 crossing between the An. funestus pyrethroid-resistant strain FUMOZ and the pyrethroid-susceptible strain FANG revealed that PBO-based nets (Olyset Plus) induced a significantly higher mortality rate (99.1%) than pyrethroid-only nets (Olyset) (56.7%) (p < 0.0001). The blood-feeding rate was higher in Olyset compared to Olyset Plus (11.6% vs. 5.6%; p = 0.013). Genotyping the CYP6P9a/b and the intergenic 6.5 kb structural variant (SV) resistance alleles showed that, for both nets, homozygote-resistant mosquitoes have a greater ability to blood-feed than the susceptible mosquitoes. Homozygote-resistant genotypes significantly survived more with Olyset after cone assays (e.g., CYP6P9a OR = 34.6; p < 0.0001) than homozygote-susceptible mosquitoes. A similar but lower correlation was seen with Olyset Plus (OR = 6.4; p < 0.001). Genotyping EHT samples confirmed that CYP6P9a/b and 6.5 kb_SV homozygote-resistant mosquitoes survive and blood-feed significantly better than homozygote-susceptible mosquitoes when exposed to Olyset. Our findings highlight the negative impact of P450-based resistance on pyrethroid-only nets, further supporting that PBO nets, such as Olyset Plus, are a better solution in areas of P450-mediated resistance to pyrethroids.

Multi-omics analysis identifies a CYP9K1 haplotype conferring pyrethroid resistance in the malaria vector Anopheles funestus in East Africa

Metabolic resistance to pyrethroids is a menace to the continued effectiveness of malaria vector controls. Its molecular basis is complex and varies geographically across Africa. Here, we used a multi‐omics approach, followed‐up with functional validation to show that a directionally selected haplotype of a cytochrome P450, CYP9K1 is a major driver of resistance in Anopheles funestus. A PoolSeq GWAS using mosquitoes alive and dead after permethrin exposure, from Malawi and Cameroon, detected candidate genomic regions, but lacked consistency across replicates. Targeted sequencing of candidate resistance genes detected several SNPs associated with known pyrethroid resistance QTLs. The most significant SNPs were in the cytochrome P450 CYP304B1 (Cameroon), CYP315A1 (Uganda) and the ABC transporter gene ABCG4 (Malawi). However, when comparing field resistant mosquitoes to laboratory susceptible, the pyrethroid resistance locus rp1 and SNPs around the ABC transporter ABCG4 were consistently significant, except for Uganda where SNPs in the P450 CYP9K1 was markedly significant. In vitro heterologous metabolism assays with recombinant CYP9K1 revealed that it metabolises type II pyrethroid (deltamethrin; 64% depletion) but not type I (permethrin; 0%), while moderately metabolising DDT (17%). CYP9K1 exhibited reduced genetic diversity in Uganda underlying an extensive selective sweep. Furthermore, a glycine to alanine (G454A) amino acid change in CYP9K1 was fixed in Ugandan mosquitoes but not in other An. funestus populations. This study sheds further light on the evolution of metabolic resistance in a major malaria vector by implicating more genes and variants that can be used to design field‐applicable markers to better track resistance Africa‐wide.

Pyrethroid resistance in the New World malaria vector Anopheles albimanus is mediated by cytochrome P450 CYP6P5

Pyrethroid resistance in the malaria vector Anopheles albimanus presents an obstacle to malaria elimination in the Americas. Here, An. albimanus CYP6P5 (the most overexpressed P450 in a Peruvian population) was functionally characterized. Recombinant CYP6P5 metabolized the type II pyrethroids, deltamethrin and α-cypermethrin with comparable affinities (K_M of 3.3 μM ± 0.4 and 3.6 μM ± 0.5, respectively), but exhibited a 2.7-fold higher catalytic rate for α-cypermethrin (k_cat of 6.02 min^-1 ± 0.2) versus deltamethrin (2.68 min^-1 ± 0.09). Time-course assays revealed progressive depletion of the above pyrethroids with production of four HPLC-detectable metabolites. Low depletion was obtained with type I pyrethroid, permethrin. Transgenic expression in Drosophila melanogaster demonstrated that overexpression of CYP6P5 alone conferred type II pyrethroid resistance, with only 16% and 55.3% mortalities in flies exposed to 0.25% α-cypermethrin and 0.15% deltamethrin, respectively. Synergist bioassays using P450 inhibitor piperonylbutoxide significantly recovered susceptibility (mortality = 73.6%, p < 0.001) in synergized flies exposed to 4% piperonylbutoxide, plus 0.25% α-cypermethrin, compared to non-synergized flies (mortality = 4.9%). Moderate resistance was also observed towards 4% DDT. These findings established the preeminent role of CYP6P5 in metabolic resistance in An. albimanus, highlighting challenges associated with deployment of insecticide-based control tools in the Americas.

CYP81A68 confers metabolic resistance to ALS and ACCase-inhibiting herbicides and its epigenetic regulation in Echinochloa crus-galli

Long-term and excessive herbicide use has led to some environmental concerns and especially, herbicide resistance evolution in weeds. Here, we confirmed acetolactate synthase (ALS) inhibiting herbicide penoxsulam resistance and cross resistance to acetyl-coenzyme carboxylase (ACCase) inhibiting herbicides (cyhalofop-butyl and metamifop) in a global weed Echinochloa crus-galli population resistant to these herbicides (R). Penoxsulam metabolism study indicated that degradation rate was significantly higher in R than susceptible E. crus-galli population (S). RNA-sequencing revealed that a cytochrome P450 (P450) gene, CYP81A68, expressed higher in R versus S. Rice seedlings overexpressing this CYP81A68 gene are resistant to penoxsulam, cyhalofop-butyl and metamifop, and penoxsulam resistance is due to enhanced metabolism via O-demethylation. Deletion analysis of the CYP81A68 gene promoter identified an efficient region, in which differential methylation of CpG islands occurred between R and S. Collectively, these results demonstrate that upregulation of E. crus-galli CYP81A68 gene endows generalist metabolic resistance to commonly used ALS- and ACCase-inhibiting herbicides in rice fields and epigenetic regulation may play a role in the resistance evolution. This research could contribute to strategies reducing herbicide environmental impacts by judicious selection of alternative herbicide and non-chemical control tactics.

A 6.5kb Intergenic Structural Variation Exacerbates the Fitness Cost of P450-Based Metabolic Resistance in the Major African Malaria Vector Anopheles funestus

Metabolic-based resistance to insecticides limit the control of medically important pests, and it is extremely detrimental in the ongoing struggle to control disease vectors. Elucidating the fitness cost of metabolic resistance in major malaria vectors is vital for successful resistance management. We established the fitness cost of the 6.5kb structural variant (6.5kb-sv) between the duplicated CYP6P9a/b P450s using the hybrid strain generated from the crossing between two An. funestus laboratory strains. Furthermore, we assessed the cumulative impact of this marker with the duplicated P450 genes. We established that individuals that were homozygote for the resistant structural variant (SV) presented reduced fecundity and slow development relative to those that were homozygote for the susceptible SV. Furthermore, we observed that 6.5kb act additively with CYP6P9a and CYP6P9b to exacerbate the reduced fecundity and the increased development time of resistant mosquitoes since double/triple homozygote susceptible (SS/SS/SS) significantly laid more eggs and developed faster than other genotypes. Moreover, a restoration of susceptibility was noted over 10 generations in the insecticide-free environment with an increased proportion of susceptible individuals. This study highlights the negative impact of multiple P450-based resistance on the key physiological traits of malaria vectors. Such high fitness costs suggest that in the absence of selection pressure, the resistant individuals will be outcompeted in the field. Therefore, this should encourage future strategies based on the rotation of insecticides to reduce selection pressure and to slow the spread of pyrethroid resistance.

Towards understanding transfluthrin efficacy in a pyrethroid-resistant strain of the malaria vector Anopheles funestus with special reference to cytochrome P450-mediated detoxification

Malaria vector control interventions rely heavily on the application of insecticides against anopheline mosquitoes, in particular the fast-acting pyrethroids that target insect voltage-gated sodium channels (VGSC). Frequent applications of pyrethroids have resulted in resistance development in the major malaria vectors including Anopheles funestus, where resistance is primarily metabolic and driven by the overexpression of microsomal cytochrome P450 monooxygenases (P450s). Here we examined the pattern of cross-resistance of the pyrethroid-resistant An. funestus strain FUMOZ-R towards transfluthrin and multi-halogenated benzyl derivatives, permethrin, cypermethrin and deltamethrin in comparison to the susceptible reference strain FANG. Transfluthrin and two multi-fluorinated derivatives exhibited micromolar potency - comparable to permethrin - to functionally expressed dipteran VGSC in a cell-based cation influx assay. The activity of transfluthrin and its derivatives on VGSC was strongly correlated with their contact efficacy against strain FUMOZ-R, although no such correlation was obtained for the other pyrethroids due to their rapid detoxification by the resistant strain. The low resistance levels for transfluthrin and derivatives in strain FUMOZ-R were only weakly synergized by known P450 inhibitors such as piperonyl butoxide (PBO), triflumizole and 1-aminobenzotriazole (1-ABT). In contrast, deltamethrin toxicity in FUMOZ-R was synergized > 100-fold by all three P450 inhibitors. The biochemical profiling of a range of fluorescent resorufin and coumarin compounds against FANG and FUMOZ-R microsomes identified 7-benzyloxymethoxy-4-trifluoromethylcoumarin (BOMFC) as a highly sensitive probe substrate for P450 activity. BOMFC was used to develop a fluorescence-based high-throughput screening assay to measure the P450 inhibitory action of potential synergists. Azole fungicides prochloraz and triflumizole were identified as extremely potent nanomolar inhibitors of microsomal P450s, strongly synergizing deltamethrin toxicity in An. funestus. Overall, the present study contributed to the understanding of transfluthrin efficacy at the molecular and organismal level and identified azole compounds with potential to synergize pyrethroid efficacy in malaria vectors.

•

Transfluthrin and derivatives lack cross-resistance in resistant Anopheles funestus.

•

Pyrethroid resistance in An. funestus is strongly synergized by azole fungicides.

•

BOMFC is a highly active fluorescent probe substrate for microsomal cytochrome P450 monooxygenases in An. funestus.

•

Azole fungicides are nanomolar inhibitors of microsomal cytochrome P450 monooxygenases in An. funestus.

Biochemical profiling of functionally expressed CYP6P9 variants of the malaria vector Anopheles funestus with special reference to cytochrome b5 and its role in pyrethroid and coumarin substrate metabolism

Cytochrome P450 monooxygenases (P450s) are well studied enzymes catalyzing the oxidative metabolism of xenobiotics in insects including mosquitoes. Their duplication and upregulation in agricultural and public health pests such as anopheline mosquitoes often leads to an enhanced metabolism of insecticides which confers resistance. In the laboratory strain Anopheles funestus FUMOZ-R the duplicated P450s CYP6P9a and CYP6P9b are highly upregulated and proven to confer pyrethroid resistance. Microsomal P450 activity is regulated by NADPH cytochrome P450 oxidoreductase (CPR) required for electron transfer, whereas the modulatory role of cytochrome b₅ (CYB5) on insect P450 activity is less clear. In previous studies CYP6P9a and CYP6P9b were recombinantly expressed in tandem with An. gambiae CPR using E. coli-expression systems and CYB5 added to the reaction mix to enhance activity. However, the precise role of CYB5 on substrate turn-over when combined with CYP6P9a and CYP6P9b remains poorly investigated, thus one objective of our study was to address this knowledge gap. In contrast to the CYP6P9 variants, the expression levels of both CYB5 and CPR were not upregulated in the pyrethroid resistant FUMOZ-R strain when compared to the susceptible FANG strain, suggesting no immediate regulatory role of these genes in pyrethroid resistance in FUMOZ-R. Here, for the first time we recombinantly expressed CYP6P9a and CYP6P9b from An. funestus in a baculovirus expression system using High-5 insect cells. Co-expression of each enzyme with CPR from either An. gambiae or An. funestus did not reveal noteworthy differences in catalytic capacity. Whereas the co-expression of An. funestus CYB5 - tested at different multiplicity of infection (MOI) ratios - resulted in a significantly higher metabolization of coumarin substrates as measured by fluorescence assays. This was confirmed by Michaelis-Menten kinetics using the most active substrate, 7-benzyloxymethoxy-4-trifluoromethylcoumarin (BOMFC). We observed a similar increase in coumarin substrate turnover by adding human CYB5 to the reaction mix. Finally, we compared by UPLC-MS/MS analysis the depletion rate of deltamethrin and the formation of 4'OH-deltamethrin by recombinantly expressed CYP6P9a and CYP6P9b with and without CYB5 and detected no difference in the extent of deltamethrin metabolism. Our results suggest that co-expression (or addition) of CYB5 with CYP6P9 variants, recombinantly expressed in insect cells, can significantly enhance their metabolic capacity to oxidize coumarins, but not deltamethrin.

Over-expression in cis of the midgut P450 CYP392A16 contributes to abamectin resistance in Tetranychus urticae

Cytochrome P450 mediated metabolism is a well-known mechanism of insecticide resistance. However, to what extent qualitative or quantitative changes are responsible for increased metabolism, is not well understood. Increased expression of P450 genes is most often reported, but the underlying regulatory mechanisms remain widely unclear. In this study, we investigate CYP392A16, a P450 from the polyphagous and major agricultural pest Tetranychus urticae. High expression levels of CYP392A16 and in vitro metabolism assays have previously associated this P450 with abamectin resistance. Here, we show that CYP392A16 is primarily localized in the midgut epithelial cells, as indicated by immunofluorescence analysis, a finding also supported by a comparison between feeding and contact toxicity bioassays. Silencing via RNAi of CYP392A16 in a highly resistant T. urticae population reduced insecticide resistance levels from 3400- to 1900- fold, compared to the susceptible reference strain. Marker-assisted backcrossing, using a single nucleotide polymorphism (SNP) found in the CYP392A16 allele from the resistant population, was subsequently performed to create congenic lines bearing this gene in a susceptible genetic background. Toxicity assays indicated that the allele derived from the resistant strain confers 3.6-fold abamectin resistance compared to the lines with susceptible genetic background. CYP392A16 is over-expressed at the same levels in these lines, pointing to cis-regulation of gene expression. In support of that, functional analysis of the putative promoter region from the resistant and susceptible parental strains revealed a higher reporter gene expression, confirming the presence of cis-acting regulatory mechanisms.

MAP4K4 controlled transcription factor POUM1 regulates PxABCG1 expression influencing Cry1Ac resistance in Plutella xylostella (L.)

Deciphering the molecular mechanisms of insect resistance to Bacillus thuringiensis (Bt) based biotechnology products including Bt sprays and Bt crops is critical for the long-term application of Bt technology. Previously, we established that down-regulation of the ABC transporter gene PxABCG1, trans-regulated by the MAPK signaling pathway, contributed to high-level resistance to Bt Cry1Ac toxin in diamondback moth, Plutella xylostella (L.). However, the underlying transcriptional regulatory mechanism was unknown. Herein, we identified putative binding sites (PBSs) of the transcription factor (TF) POUM1 in the PxABCG1 promoter and used a dual-luciferase reporter assay (DLRA) and yeast one-hybrid (Y1H) assay to reveal that POUM1 activates PxABCG1 via interaction with one of these sites. The expression of POUM1 was significantly decreased in the midgut tissue of Cry1Ac-resistant P. xylostella strains compared to a Cry1Ac-susceptible P. xylostella strain. Silencing of POUM1 expression resulted in reduced expression of the PxABCG1 gene and an increase in larval tolerance to Bt Cry1Ac toxin in the Cry1Ac-susceptible P. xylostella strain. Furthermore, silencing of PxMAP4K4 expression increased the expression of both POUM1 and PxABCG1 genes in the Cry1Ac-resistant P. xylostella strain. These results indicate that the POUM1 induces PxABCG1 expression, while the activated MAPK cascade represses PxABCG1 expression thus reducing Cry1Ac susceptibility in P. xylostella. This result deepens our understanding of the transcriptional regulatory mechanism of midgut Cry receptor genes and the molecular basis of the evolution of Bt resistance in insects.

Large-scale (Phase III) evaluation of broflanilide 50WP (VECTRON™ T500) for indoor residual spraying for malaria vector control in Northeast Tanzania: study protocol for a two-arm, non-inferiority, cluster-randomised community trial

BACKGROUND

Indoor residual spraying (IRS) is a major method of malaria vector control across sub-Saharan Africa. Effective control is being undermined by the rapid spread of insecticide resistance. There is major investment in development of new insecticides for IRS that possess novel modes of action, long residual activity, low mammalian toxicity and minimal cross-resistance. VECTRON™ T500, a new IRS product containing the active ingredient broflanilide as a 50% wettable powder (WP), has been shown to be efficacious against pyrethroid susceptible and resistant vector species on mud and concrete substrates in experimental hut (Phase II) trials.

METHODS

A two-arm non-inferiority cluster randomized controlled trial (Phase III) will be undertaken in Muheza District, Tanga Region, Tanzania. VECTRON™ T500 will be compared to the IRS product Fludora® Fusion (clothianidin 50% WP + deltamethrin 6.25% WP). The predominant malaria vectors in the study area are pyrethroid-resistant Anopheles gambiae s.s., An. arabiensis and An. funestus s.s. Sixteen village clusters will be pair-matched on baseline vector densities and allocated to reference and intervention arms. Consenting households in the intervention arm will be sprayed with VECTRON™ T500 and those in the reference arm will be sprayed with Fludora® Fusion. Each month, CDC light traps will collect mosquitoes to estimate changes in vector density, indoor biting, sporozoite and entomological inoculation rates (EIR). Susceptibility to IRS active ingredients will be assessed using World Health Organisation (WHO) bottle bioassays. Target site and metabolic resistance mechanisms will be characterised among Anopheles field populations from both trial arms. Residual efficacy of both IRS products will be monitored for 12 months post intervention. Questionnaire and focus group discussions will explore factors that influence adherence, adverse effects and benefits of IRS.

DISCUSSION

This protocol describes a large-scale non-inferiority evaluation of a novel IRS product to reduce the density and EIR of pyrethroid-resistant Anopheles vectors. If VECTRON™ T500 proves non-inferior to Fludora® Fusion, it will be considered as an additional vector control product for malaria prevention and insecticide resistance management.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT05150808, registered on 26 November 2021. Retrospectively registered.

Distribution of acetylcholinesterase (Ace-1R) target-site G119S mutation and resistance to carbamates and organophosphates in Anopheles gambiae sensu lato populations from Cameroon

BACKGROUND

Cameroon is considering the implementation of indoor residual spraying (IRS) as a complementary measure to control malaria in the context of high pyrethroid resistance in major malaria vectors. Non-pyrethroid insecticide classes such as organophosphates and carbamates may be utilized in IRS due to widespread pyrethroid resistance. However, the success of this strategy depends on good knowledge of the resistance status of malaria vectors to carbamates and organophosphates. Here, we assessed the susceptibility profile of Anopheles gambiae sensu lato with respect to carbamates and organophosphate and the distribution of the molecular mechanism underlying resistance to these insecticides.

METHODS

Anopheles gambiae s.l. mosquitoes were collected from nine settings across the country and bio-assayed with bendiocarb, propoxur and pirimiphos-methyl. The Ace-1 target-site G119S mutation was genotyped using a TaqMan assay. To investigate the polymorphism in the Ace-1 gene, a region of 924 base pairs in a sequence of the gene was amplified from both live and dead females of An. gambiae exposed to bendiocarb.

RESULTS

Pirimiphos-methyl induced full mortality in An. gambiae s.l. from all study sites, whereas for carbamates, resistance was observed in four localities, with the lowest mortality rate recorded in Mangoum (17.78 ± 5.02% for bendiocarb and 18.61 ± 3.86% for propoxur) in the southern part of Cameroon. Anopheles coluzzii was found to be the predominant species in the northern tropical part of the country where it is sympatric with Anopheles arabiensis. In the localities situated in southern equatorial regions, this species was predominant in urban settings, while An. gambiae was the most abundant species in rural areas. The G119S Ace-1 target-site mutation was detected only in An. gambiae and only in the sites located in southern Cameroon. Phylogenetic analyses showed a clustering according to the phenotype.

CONCLUSION

The occurrence of the Ace-1 target-site substitution G119S in An. gambiae s.l. populations highlights the challenge associated with the impending deployment of IRS in Cameroon using carbamates or organophosphates. It is therefore important to think about a resistance management plan including the use of other insecticide classes such as neonicotinoids or pyrrole to guarantee the implementation of IRS in Cameroon.

The Role of Cytochrome P450s in Insect Toxicology and Resistance

Insect cytochrome P450 monooxygenases (P450s) perform a variety of important physiological functions, but it is their role in the detoxification of xenobiotics, such as natural and synthetic insecticides, that is the topic of this review. Recent advances in insect genomics and postgenomic functional approaches have provided an unprecedented opportunity to understand the evolution of insect P450s and their role in insect toxicology. These approaches have also been harnessed to provide new insights into the genomic alterations that lead to insecticide resistance, the mechanisms by which P450s are regulated, and the functional determinants of P450-mediated insecticide resistance. In parallel, an emerging body of work on the role of P450s in defining the sensitivity of beneficial insects to insecticides has been developed. The knowledge gained from these studies has applications for the management of P450-mediated resistance in insect pests and can be leveraged to safeguard the health of important beneficial insects.