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

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

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

Genomic Signatures of Microgeographic Adaptation in Anopheles coluzzii Along an Anthropogenic Gradient in Gabon

Species distributed across heterogeneous environments often evolve locally adapted populations, but understanding how these persist in the presence of homogenizing gene flow remains puzzling. In Gabon, Anopheles coluzzii, a major African malaria mosquito is found along an ecological gradient, including a sylvatic population, away of any human presence. This study identifies into the genomic signatures of local adaptation in populations from distinct environments including the urban area of Libreville, and two proximate sites 10km apart in the La Lopé National Park (LLP), a village and its sylvatic neighborhood. Whole genome re-sequencing of 96 mosquitoes unveiled ∼ 5.7millions high-quality single nucleotide polymorphisms. Coalescent-based demographic analyses suggest an ∼ 8,000-year-old divergence between Libreville and La Lopé populations, followed by a secondary contact ( ∼ 4,000 ybp) resulting in asymmetric effective gene flow. The urban population displayed reduced effective size, evidence of inbreeding, and strong selection pressures for adaptation to urban settings, as suggested by the hard selective sweeps associated with genes involved in detoxification and insecticide resistance. In contrast, the two geographically proximate LLP populations showed larger effective sizes, and distinctive genomic differences in selective signals, notably soft-selective sweeps on the standing genetic variation. Although neutral loci and chromosomal inversions failed to discriminate between LLP populations, our findings support that microgeographic adaptation can swiftly emerge through selection on standing genetic variation despite high gene flow. This study contributes to the growing understanding of evolution of populations in heterogeneous environments amid ongoing gene flow and how major malaria mosquitoes adapt to human.

Significance

Anopheles coluzzii , a major African malaria vector, thrives from humid rainforests to dry savannahs and coastal areas. This ecological success is linked to its close association with domestic settings, with human playing significant roles in driving the recent urban evolution of this mosquito. Our research explores the assumption that these mosquitoes are strictly dependent on human habitats, by conducting whole-genome sequencing on An. coluzzii specimens from urban, rural, and sylvatic sites in Gabon. We found that urban mosquitoes show de novo genetic signatures of human-driven vector control, while rural and sylvatic mosquitoes exhibit distinctive genetic evidence of local adaptations derived from standing genetic variation. Understanding adaptation mechanisms of this mosquito is therefore crucial to predict evolution of vector control strategies.

Pyrethroid resistance and gene expression profile of a new resistant An. gambiae colony from Uganda reveals multiple resistance mechanisms and overexpression of Glutathione-S-Transferases linked to survival of PBO-pyrethroid combination

Background

The effectiveness of long-lasting insecticidal nets (LLINs) are being threatened by growing resistance to pyrethroids. To restore their efficacy, a synergist, piperonyl butoxide (PBO) which inhibits cytochrome P450s has been incorporated into pyrethroid treated nets. A trial of PBO-LLINs was conducted in Uganda from 2017 and we attempted to characterize mechanisms of resistance that could impact intervention efficacy.

Methods

We established an Anopheles gambiae s.s colony in 2018 using female mosquitoes collected from Busia district in eastern Uganda. We first assessed the phenotypic resistance profile of this colony using WHO tube and net assays using a deltamethrin dose-response approach. The Busia colony was screened for known resistance markers and RT-qPCR targeting 15 genes previously associated with insecticide resistance was performed.

Results

The Busia colony had very high resistance to deltamethrin, permethrin and DDT. In addition, the colony had moderate resistance to alpha-cypermethrin and lambda-cyhalothrin but were fully susceptible to bendiocarb and fenitrothion. Exposure to PBO in combination with permethrin and deltamethrin resulted in higher mortality rates in both net and tube assays, with a higher mortality observed in net assays than tube assays. The kdr marker, Vgsc-995S was at very high frequency (91.7-98.9%) whilst the metabolic markers Coeae1d and Cyp4j5-L43F were at very low (1.3% - 11.5%) and moderate (39.5% - 44.7%) frequencies respectively. Our analysis showed that gene expression pattern in mosquitoes exposed to deltamethrin, permethrin or DDT only were similar in comparison to the susceptible strain and there was significant overexpression of cytochrome P450s, glutathione-s-transferases (GSTs) and carboxyl esterases (COEs). However, mosquitoes exposed to both PBO and pyrethroid strikingly and significantly only overexpressed closely related GSTs compared to unexposed mosquitoes while major cytochrome P450s were underexpressed.

Conclusions

The high levels of pyrethroid resistance observed in Busia appears associated with a wide range of metabolic gene families.

The impact of agrochemical pollutant mixtures on the selection of insecticide resistance in the malaria vector Anopheles gambiae: insights from experimental evolution and transcriptomics

BACKGROUND

There are several indications that pesticides used in agriculture contribute to the emergence and spread of resistance of mosquitoes to vector control insecticides. However, the impact of such an indirect selection pressure has rarely been quantified and the molecular mechanisms involved are still poorly characterized. In this context, experimental selection with different agrochemical mixtures was conducted in Anopheles gambiae. The multi-generational impact of agrochemicals on insecticide resistance was evaluated by phenotypic and molecular approaches.

METHODS

Mosquito larvae were selected for 30 generations with three different agrochemical mixtures containing (i) insecticides, (ii) non-insecticides compounds, and (iii) both insecticide and non-insecticide compounds. Every five generations, the resistance of adults to deltamethrin and bendiocarb was monitored using bioassays. The frequencies of the kdr (L995F) and ace1 (G119S) target-site mutations were monitored every 10 generations. RNAseq was performed on all lines at generation 30 in order to identify gene transcription level variations and polymorphisms associated with each selection regime.

RESULTS

Larval selection with agrochemical mixtures did not affect bendiocarb resistance and did not select for ace1 mutation. Contrastingly, an increased deltamethrin resistance was observed in the three selected lines. Such increased resistance was not majorly associated with the presence of kdr L995F mutation in selected lines. RNA-seq identified 63 candidate resistance genes over-transcribed in at least one selected line. These include genes coding for detoxification enzymes or cuticular proteins previously associated with insecticide resistance, and other genes potentially associated with chemical stress response. Combining an allele frequency filtering with a Bayesian FST-based genome scan allowed to identify genes under selection across multiple genomic loci, supporting a multigenic adaptive response to agrochemical mixtures.

CONCLUSION

This study supports the role of agrochemical contaminants as a significant larval selection pressure favouring insecticide resistance in malaria vectors. Such selection pressures likely impact kdr mutations and detoxification enzymes, but also more generalist mechanisms such as cuticle resistance, which could potentially lead to cross-tolerance to unrelated insecticide compounds. Such indirect effect of global landscape pollution on mosquito resistance to public health insecticides deserves further attention since it can affect the nature and dynamics of resistance alleles circulating in malaria vectors and impact the efficacy of control vector strategies.

LLIN Evaluation in Uganda Project (LLINEUP) - Plasmodium infection prevalence and genotypic markers of insecticide resistance in Anopheles vectors from 48 districts of Uganda

Background

In 2017-2019, we conducted a large-scale, cluster-randomised trial (LLINEUP) to evaluate long-lasting insecticidal nets (LLINs) treated with a pyrethroid insecticide plus the synergist piperonyl butoxide (PBO LLINs), as compared to conventional, pyrethroid-only LLINs across 104 health sub-districts (HSDs) in Uganda. In LLINEUP, and similar trials in Tanzania, PBO LLINs were found to provide greater protection against malaria than conventional LLINs, reducing parasitaemia and vector density. In the LLINEUP trial, cross-sectional entomological surveys were carried out at baseline and then every 6 months for two years. In each survey, ten households per HSD were randomly selected for indoor household entomological collections.

Results

Overall, 5395 female Anopheles mosquitoes were collected from 5046 households. The proportion of mosquitoes infected with Plasmodium falciparum did not change significantly over time, while infection with non-falciparum malaria decreased in An. gambiae s.s, but not An. funestus. The frequency of genetic markers associated with pyrethroid resistance increased significantly over time, but the rate of change was not different between the two LLIN types. The knock-down resistance (kdr) mutation Vgsc-995S declined over time as Vgsc-995F, the alternative resistance mutation at this codon, increased. Vgsc-995F appears to be spreading into Uganda.

Conclusions

Distribution of LLINs in Uganda was associated with reductions in parasite prevalence and vector density, but the proportion of infective mosquitoes remained stable, suggesting that the potential for transmission persisted. The increased frequency of markers of pyrethroid resistance indicates that LLIN distribution favoured the evolution of resistance within local vectors and highlights the potential benefits of resistance management strategies.Trial registration:: This study is registered with ISRCTN, ISRCTN17516395. Registered 14 February 2017, http://www.isrctn.com/ISRCTN17516395.

Widespread occurrence of copy number variants and fixation of pyrethroid target site resistance in Anopheles gambiae (s.l.) from southern Côte d'Ivoire

Resistance to pyrethroid and organophosphate insecticides in the malaria vector Anopheles gambiae (s.l.) is conferred by a variety of genetic mutations, including single nucleotide polymorphisms (SNPs) and copy number variants (CNVs). Knowledge of the distribution of these mutations in mosquito populations is a prerequisite for establishing better strategies for their management. In this study, a total of 755 Anopheles gambiae (s.l.) from southern Côte d'Ivoire were exposed to deltamethrin or pirimiphos-methyl insecticides and were screened to assess the distribution of SNPs and CNVs known or believed to confer resistance to one or other of the insecticide classes. Most individuals from the An. gambiae (s.l.) complex were identified by molecular tests as Anopheles coluzzii. Survival to deltamethrin (from 94% to 97%) was higher than to pirimiphos-methyl (from 10% to 49%). In An. gambiae (s.s.), the SNP in the Voltage Gated Sodium Channel (Vgsc) at the 995F locus (Vgsc-995F) was fixed, while other target site mutations were rare or absent (Vgsc-402L: 0%; Vgsc-1570Y: 0%, Acetylcholinesterase Acel-280S: 14%). In An. coluzzii, Vgsc-995F was the target site SNP found at highest frequency (65%) followed by other target site mutations (Vgsc-402L: 36%; Vgsc-1570Y: 0.33%; Acel-280S: 45%). The Vgsc-995S SNP was not present. The presence of the Ace1-280S SNP was found to be significantly linked to the presence of the Ace1-CNV, Ace1_AgDup. Significant association was found between the presence of the Ace1_AgDup and pirimiphos-methyl resistance in An. gambiae (s.s.) but not in An. coluzzii. The deletion Ace1_Del97 was found in one specimen of An. gambiae (s.s.). Four CNVs in the Cyp6aa/Cyp6p gene cluster, which contains genes of known importance for resistance, were detected in An. coluzzii, the most frequent being Dup 7 (42%) and Dup 14 (26%). While none of these individual CNV alleles were significantly associated with resistance, copy number in the Cyp6aa gene region in general was associated with increased resistance to deltamethrin. Elevated expression of Cyp6p3 was nearly associated with deltamethrin resistance, although there was no association of resistance with copy number. Use of alternative insecticides and control methods to arrest resistance spread in An. coluzzii populations is merited.

Prevalence and predictors of adverse events following exposure to long-lasting insecticidal nets used for malaria prevention: a community based cross-sectional study in the Democratic Republic of the Congo

BACKGROUND

Malaria morbidity and mortality increase in the Democratic Republic of the Congo (DRC) may be the consequence of the low utilization rate of long-lasting insecticidal nets (LLINs) resulting from poor compliance due to adverse events (AEs). This study aimed at determining the prevalence and predictors of AEs following the mass distribution of LLINs in the Kisantu Health Zone (KHZ), a high malaria-endemic region in the DRC.

METHODS

A community-based cross-sectional study embedded was conducted within a randomized controlled trial (RCT) after the mass distribution of LLINs in 30 villages located in DRC KHZ. A three-stage sampling method was used without replacement to select 1790 children. Data was collected on adverse events (AEs) using a reporting form and information on demographics, nutritional status, and house characteristics. This was done using a structured questionnaire administered to household heads. Logistic regression models were used to identify predictors of AEs following the mass distribution of LLINs.

RESULT

In a total of 1790 children enrolled, 17.8% (95% CI 16.1-19.7) experienced AEs. The most common AEs were respiratory-related (61%). Around 60% of AEs occurred within 24 h of use, and 51% were resolved without treatment. Sleeping under deltamethrin LLINs (Adjusted OR, 95% CI 5.5 [3.8-8.0]) and zinc roofing (Adjusted OR, 95% CI 1.98 [1.1-3.57]) were associated with the risk of reporting an AE following the mass distribution of LLINs.

CONCLUSION

Approximately 1 out of 5 children had an AE within 24 h following LLIN use. These adverse events were often respiratory-related. LLINs and roofing types were associated with a higher risk of reporting AEs. However, further research using a robust study design is needed to confirm these findings. Future studies should design and implement interventions aiming to reduce AEs and improve compliance with LLINs.

Malaria among children under 10 years in 4 endemic health areas in Kisantu Health Zone: epidemiology and transmission

BACKGROUND

The Democratic Republic of the Congo (DRC) is the second most malaria-affected country in the world with 21,608,681 cases reported in 2019. The Kongo Central (KC) Province has a malaria annual incidence of 163 cases/per 1000 inhabitants which are close to the national average of 153.4/1000. However, the malaria prevalence varies both between and within health zones in this province. The main objective of this study was to describe the epidemiology and transmission of malaria among children aged 0 to 10 years in the 4 highest endemic health areas in Kisantu Health Zone (HZ) of KC in DRC.

METHODS

A community-based cross-sectional study was conducted from October to November 2017 using multi-stage sampling. A total of 30 villages in 4 health areas in Kisantu HZ were randomly selected. The prevalence of malaria was measured using a thick blood smear (TBS) and known predictors and associated outcomes were assessed. Data are described and association determinants of malaria infection were analysed.

RESULTS

A total of 1790 children between 0 and 10 years were included in 30 villages in 4 health areas of Kisantu HZ. The overall prevalence in the study area according to the TBS was 14.8% (95% CI: 13.8-16.6; range: 0-53). The mean sporozoite rate in the study area was 4.3% (95% CI: 2.6-6.6). The determination of kdr-west resistance alleles showed the presence of both L1014S and L1014F with 14.6% heterozygous L1014S/L1014F, 84.4% homozygous 1014F, and 1% homozygous 1014S. The risk factors associated with malaria infection were ground or wooden floors aOR: 15.8 (95% CI: 8.6-29.2), a moderate or severe underweight: 1.5 (1.1-2.3) and to be overweight: 1.9 (95% CI: 1.3-2.7).

CONCLUSION

Malaria prevalence differed between villages and health areas within the same health zone. The control strategy activities must be oriented by the variety in the prevalence and transmission of malaria in different areas. The policy against malaria regarding long-lasting insecticidal nets should be based on the evidence of metabolic resistance.

Characterizing pyrethroid resistance and mechanisms in Anopheles gambiae (s.s.) and Anopheles arabiensis from 11 districts in Uganda

Insecticide resistance threatens recent progress on malaria control in Africa. To characterize pyrethroid resistance in Uganda, Anopheles gambiae (s.s.) and Anopheles arabiensis were analyzed from 11 sites with varied vector control strategies. Mosquito larvae were collected between May 2018 and December 2020. Sites were categorized as receiving no indoor-residual spraying ('no IRS', n ​= ​3); where IRS was delivered from 2009 to 2014 and in 2017 and then discontinued ('IRS stopped', n ​= ​4); and where IRS had been sustained since 2014 ('IRS active', n ​= ​4). IRS included bendiocarb, pirimiphos methyl and clothianidin. All sites received long-lasting insecticidal nets (LLINs) in 2017. Adult mosquitoes were exposed to pyrethroids; with or without piperonyl butoxide (PBO). Anopheles gambiae (s.s.) and An. arabiensis were identified using PCR. Anopheles gambiae (s.s.) were genotyped for Vgsc-995S/F, Cyp6aa1, Cyp6p4-I236M, ZZB-TE, Cyp4j5-L43F and Coeae1d, while An. arabiensis were examined for Vgsc-1014S/F. Overall, 2753 An. gambiae (s.l.), including 1105 An. gambiae (s.s.) and 1648 An. arabiensis were evaluated. Species composition varied by site; only nine An. gambiae (s.s.) were collected from 'IRS active' sites, precluding species-specific comparisons. Overall, mortality following exposure to permethrin and deltamethrin was 18.8% (148/788) in An. gambiae (s.s.) and 74.6% (912/1222) in An. arabiensis. Mortality was significantly lower in An. gambiae (s.s.) than in An. arabiensis in 'no IRS' sites (permethrin: 16.1 vs 67.7%, P ​< ​0.001; deltamethrin: 24.6 vs 83.7%, P ​< ​0.001) and in 'IRS stopped' sites (permethrin: 11.3 vs 63.6%, P ​< ​0.001; deltamethrin: 25.6 vs 88.9%, P ​< ​0.001). When PBO was added, mortality increased for An. gambiae (s.s.) and An. arabiensis. Most An. gambiae (s.s.) had the Vgsc-995S/F mutation (95% frequency) and the Cyp6p4-I236M resistance allele (87%), while the frequency of Cyp4j5 and Coeae1d were lower (52% and 55%, respectively). Resistance to pyrethroids was widespread and higher in An. gambiae (s.s.). Where IRS was active, An. arabiensis dominated. Addition of PBO to pyrethroids increased mortality, supporting deployment of PBO LLINs. Further surveillance of insecticide resistance and assessment of associations between genotypic markers and phenotypic outcomes are needed to better understand mechanisms of pyrethroid resistance and to guide vector control.

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.

Molecular Drivers of Multiple and Elevated Resistance to Insecticides in a Population of the Malaria Vector Anopheles gambiae in Agriculture Hotspot of West Cameroon

(1) Background: Malaria remains a global public health problem. Unfortunately, the resistance of malaria vectors to commonly used insecticides threatens disease control and elimination efforts. Field mosquitoes have been shown to survive upon exposure to high insecticide concentrations. The molecular mechanisms driving this pronounced resistance remain poorly understood. Here, we elucidated the pattern of resistance escalation in the main malaria vector Anopheles gambiae in a pesticide-driven agricultural hotspot in Cameroon and its impact on vector control tools; (2) Methods: Larval stages and indoor blood-fed female mosquitoes (F0) were collected in Mangoum in May and November and forced to lay eggs; the emerged mosquitoes were used for WHO tube, synergist and cone tests. Molecular identification was performed using SINE PCR, whereas TaqMan-based PCR was used for genotyping of L1014F/S and N1575Y kdr and the G119S-ACE1 resistance markers. The transcription profile of candidate resistance genes was performed using qRT-PCR methods. Characterization of the breeding water and soil from Mangoum was achieved using the HPLC technique; (3) Results: An. gambiae s.s. was the only species in Mangoum with 4.10% infection with Plasmodium. These mosquitoes were resistant to all the four classes of insecticides with mortality rates <7% for pyrethroids and DDT and <54% for carbamates and organophophates. This population also exhibited high resistance intensity to pyrethroids (permethrin, alpha-cypermethrin and deltamethrin) after exposure to 5× and 10× discriminating doses. Synergist assays with PBO revealed only a partial recovery of susceptibility to permethrin, alpha-cypermethrin and deltamethrin. Only PBO-based nets (Olyset plus and permaNet 3.0) and Royal Guard showed an optimal efficacy. A high amount of alpha-cypermethrin was detected in breeding sites (5.16-fold LOD) suggesting ongoing selection from agricultural pesticides. The 1014F-kdr allele was fixed (100%) whereas the 1575Y-kdr (37.5%) and the 119S Ace-1R (51.1%) were moderately present. Elevated expression of P450s, respectively, in permethrin and deltamethrin resistant mosquitoes [CYP6M2 (10 and 34-fold), CYP6Z1(17 and 29-fold), CYP6Z2 (13 and 65-fold), CYP9K1 (13 and 87-fold)] supports their role in the observed resistance besides other mechanisms including chemosensory genes as SAP1 (28 and 13-fold), SAP2 (5 and 5-fold), SAP3 (24 and 8-fold) and cuticular genes as CYP4G16 (6 and 8-fold) and CYP4G17 (5 and 27-fold). However, these candidate genes were not associated with resistance escalation as the expression levels did not differ significantly between 1×, 5× and 10× surviving mosquitoes; (4) Conclusions: Intensive and multiple resistance is being selected in malaria vectors from a pesticide-based agricultural hotspot of Cameroon leading to loss in the efficacy of pyrethroid-only nets. Further studies are needed to decipher the molecular basis underlying such resistance escalation to better assess its impact on control interventions.

Organophosphate Insecticide Exposure Impacts Reproductive Success in Insensitive Acetylcholinesterase Anopheles gambiae Mosquitoes

Extensive use of insecticides has led to the selection of resistance alleles in malaria vectors threatening the control programs. Even if mosquitoes are not killed directly in the contact of insecticide-treated bed nets, their capacity to transmit malaria parasite could be decreased because of the consequences on their life-history traits after repeated exposure. The current work investigated the effects of organochlorine, carbamate, organophosphate, and pyrethroid insecticide exposure on the reproductive success in Anopheles gambiae s.s. Two Anopheles gambiae strains, AcerKis, KisKdr, were used. According to WHO recommendations, female mosquitoes of these resistant strains were exposed to discriminant doses of DDT, chlorpyriphos-methyl, bendiocarb, and permethrin insecticides. Surviving mosquitoes were then fed and allowed to lay eggs. Fecundity was assessed by examining the number of eggs per mosquito, the number of larvae per egg batch and larval hatching rates were used to evaluate the fertility. The data showed that AcerKis females surviving chlorpyriphos-methyl exposure significantly laid few eggs. No significant difference in the hatching rate was noticed in AcerKis females exposed to bendiocarb compared to their control. No significant effect on the fecundity and fertility was observed in KisKdr females exposed to permethrin. Our finding showed that organophosphate insecticides represented here by chlorpyriphos-methyl could hamper egg-laying in insensitive acetylcholinesterase An. gambiae female mosquitoes. This knowledge could help design alternative vector control strategies targeting fecundity and fertility in resistant malaria vectors.

Persistent, bioaccumulative, and toxic chemicals in insects: Current state of research and where to from here?

The ongoing decline in the biomass, abundance, and species number of insects is an established fact. Persistent, bioaccumulative, and toxic chemicals (PBTs) - persistent organic pollutants (POPs) and, in the case of our study, mercury (Hg) - play an important role, but their effect on insect populations is insufficiently investigated. Here, the current state of research on PBTs related to insects is examined with a systematic literature study using Web of Science™. We investigate time trends of research intensity compared with other organisms, insect orders and chemicals analyzed, chemicals' effects on insects, and geographical aspects. We show that research intensity increased in the early 1990s, but studies on PBTs in insects are still underrepresented compared with other organisms. The taxonomic focus lies strongly on dipterans. The predominance of studies on DDT suggests its relevance in the context of disease-vector management. Phenotypic and acute effects on insects were more often investigated than genotypic and chronic effects. Laboratory-bred insects and wild-bred insects were examined equally often, pollutant exposure and analysis were conducted predominantly in the laboratory. Mostly habitats with a medium or high human impact were studied, and natural and near-natural habitats are understudied. The sources of the substances are often unknown. Most studies were carried out in economically rich continents, including North America, Europe, and Australia. The numbers of publications dealing with Asia, South America, and Africa are comparatively low, although the control of vector-borne diseases with POPs is still intensively practiced there. We identify gaps in the research - among others, refined analytical methods for biomarkers and for the examination of chronic effects, combinations of field and laboratory experiments to analyze the same problem, and a global approach for the monitoring of PBTs will be needed for accelerating the dearly needed progress in the research of PBTs in insects.

Detection and population genetic analysis of kdr L1014F variant in eastern Ethiopian Anopheles stephensi

Anopheles stephensi is a malaria vector that has been recently introduced into East Africa, where it threatens to increase malaria disease burden. The use of insecticides, especially pyrethroids, is still one of the primary malaria vector control strategies worldwide. The knockdown resistance (kdr) mutation in the IIS6 transmembrane segment of the voltage-gated sodium channel (vgsc) is one of the main molecular mechanisms of pyrethroid resistance in Anopheles. Extensive pyrethroid resistance in An. stephensi has been previously reported in Ethiopia. Thus, it is important to determine whether or not the kdr mutation is present in An. stephensi populations in Ethiopia to inform vector control strategies. In the present study, the kdr locus was analyzed in An. stephensi collected from ten urban sites (Awash Sebat Kilo, Bati, Dire Dawa, Degehabur, Erer Gota, Godey, Gewane, Jigjiga, Semera, and Kebridehar) situated in Somali, Afar, and Amhara regions, and Dire Dawa Administrative City, to evaluate the frequency and evolution of kdr mutations and the association of the mutation with permethrin resistance phenotypes. Permethrin is one of the pyrethroid insecticides used for vector control in eastern Ethiopia. DNA extractions were performed on adult mosquitoes from CDC light trap collections and those raised from larval and pupal collections. PCR and targeted sequencing were used to analyze the IIS6 transmembrane segment of the vgsc gene. Of 159 An. stephensi specimens analyzed from the population survey, nine (5.7%) carried the kdr mutation (L1014F). An. stephensi with kdr mutations were only observed from Bati, Degehabur, Dire Dawa, Gewane, and Semera. We further selected randomly twenty resistant and twenty susceptible An. stephensi mosquitoes from Dire Dawa post-exposure to permethrin and investigated the role of kdr in pyrethroid resistance by comparing the vgsc gene in the two populations. We found no kdr mutations in the permethrin-resistant mosquitoes. Population genetic analysis of the sequences, including neighboring introns, revealed limited evidence of non-neutral evolution (e.g., selection) at this locus. The low kdr mutation frequency detected and the lack of kdr mutation in the permethrin-resistant mosquitoes suggest the existence of other molecular mechanisms of pyrethroid resistance in eastern Ethiopian An. stephensi.

In vivo functional validation of the V402L voltage gated sodium channel mutation in the malaria vector An. gambiae

BACKGROUND

Pyrethroids are the most widely used insecticides for the control of malaria transmitting Anopheles gambiae mosquitoes and rapid increase in resistance to this insecticide class is of major concern. Pyrethroids target the Voltage Gated Sodium Channels (VGSCs), that have a key role in the normal function of the mosquitoes' nervous system. VGSC mutations L995F and L995S have long been associated with pyrethroid resistance and screening for their presence is routine in insecticide resistance management programs. Recently, a VGSC haplotype containing two amino acid substitutions associated with resistance in other species, V402L and I1527T, was identified. These two VGSC mutations are found in tight linkage and are mutually exclusive to the classical L995F/S mutations.

RESULTS

We identify the presence of the V402L-I1527T haplotype in resistant An. coluzzii colonized strains and in field populations from Burkina Faso, at frequencies higher than previously reported; in some cases almost reaching fixation. Functional validation of V402L in insecticide resistance using a CRISPR/Cas9 genome modified line showed that it confers reduced mortality after exposure to all tested pyrethroids and DDT, but at lower levels compared to L995F. In contrast to L995F however, no fitness costs were identified for mosquitoes carrying V402L under laboratory conditions.

CONCLUSION

The V402L substitution confers pyrethroid resistance in An. gambiae in the absence of any other VGSC substitution and/or alternative resistance mechanisms. The lower fitness cost associated with this kdr mutation may provide a selective advantage over the classical kdr in some settings and genotyping at this locus should be added in the list of resistant alleles for routine screening.

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

Background

Existing mechanisms of insecticide resistance are known to help the survival of mosquitoes following contact with chemical compounds, even though they could negatively affect the life-history traits of resistant malaria vectors. In West Africa, the knockdown resistance mechanism kdr^R (L1014F) is the most common. However, little knowledge is available on its effects on mosquito life-history traits. The fitness effects associated with this knockdown resistance allele in Anopheles gambiae sensu stricto (s.s.) were investigated in an insecticide-free laboratory environment.

Methods

The life-history traits of Kisumu (susceptible) and KisKdr (kdr resistant) strains of An. gambiae s.s. were compared. Larval survivorship and pupation rate were assessed as well as fecundity and fertility of adult females. Female mosquitoes of both strains were directly blood fed through artificial membrane assays and then the blood-feeding success, blood volume and adult survivorship post-blood meal were assessed.

Results

The An. gambiae mosquitoes carrying the kdr^R allele (KisKdr) laid a reduced number of eggs. The mean number of larvae in the susceptible strain Kisumu was three-fold overall higher than that seen in the KisKdr strain with a significant difference in hatching rates (81.89% in Kisumu vs 72.89% in KisKdr). The KisKdr larvae had a significant higher survivorship than that of Kisumu. The blood-feeding success was significantly higher in the resistant mosquitoes (84%) compared to the susceptible ones (34.75%). However, the mean blood volume was 1.36 µL/mg, 1.45 µL/mg and 1.68 µL/mg in Kisumu, homozygote and heterozygote KisKdr mosquitoes, respectively. After blood-feeding, the heterozygote KisKdr mosquitoes displayed highest survivorship when compared to that of Kisumu.

Conclusions

The presence of the knockdown resistance allele appears to impact the life-history traits, such as fecundity, fertility, larval survivorship, and blood-feeding behaviour in An. gambiae. These data could help to guide the implementation of more reliable strategies for the control of malaria vectors.

Investigating molecular mechanisms of insecticide resistance in the Eastern Democratic Republic of the Congo

BACKGROUND

Malaria vector control in the Democratic Republic of the Congo is plagued by several major challenges, including inadequate infrastructure, lack of access to health care systems and preventative measures, and more recently the widespread emergence of insecticide resistance among Anopheles mosquitoes. Across 26 provinces, insecticide resistance has been reported from multiple sentinel sites. However, to date, investigation of molecular resistance mechanisms among Anopheles vector populations in DRC has been more limited.

METHODS

Adult Anopheles gambiae sensu lato (s.l.) and Anopheles funestus s.l. were collected from two sites in Sud-Kivu province and one site in Haut-Uélé province and PCR-screened for the presence of 11 resistance mutations, to provide additional information on frequency of resistance mechanisms in the eastern DRC, and to critically evaluate the utility of these markers for prospective country-wide resistance monitoring.

RESULTS

L1014F-kdr and L1014S-kdr were present in 75.9% and 56.7% of An. gambiae s.l. screened, respectively, with some individuals harbouring both resistant alleles. Across the three study sites, L43F-CYP4J5 allele frequency ranged from 0.42 to 0.52, with evidence for ongoing selection. G119S-ace1 was also identified in all sites but at lower levels. A triple mutant haplotype (comprising the point mutation CYP6P4-I236M, the insertion of a partial Zanzibar-like transposable element and duplication of CYP6AA1) was present at high frequencies. In An. funestus s.l. cis-regulatory polymorphisms in CYP6P9a and CYP6P9b were detected, with allele frequencies ranging from 0.82 to 0.98 and 0.65 to 0.83, respectively.

CONCLUSIONS

This study screened the most up-to-date panel of DNA-based resistance markers in An. gambiae s.l. and An. funestus s.l. from the eastern DRC, where resistance data is lacking. Several new candidate markers (CYP4J5, G119S-ace1, the triple mutant, CYP6P9a and CYP6P9b) were identified, which are diagnostic of resistance to major insecticide classes, and warrant future, larger-scale monitoring in the DRC to inform vector control decisions by the National Malaria Control Programme.

A whole transcriptomic approach provides novel insights into the molecular basis of organophosphate and pyrethroid resistance in Anopheles arabiensis from Ethiopia

The development of insecticide resistance in malaria vectors is of increasing concern in Ethiopia because of its potential implications for vector control failure. To better elucidate the specificity of resistance mechanisms and to facilitate the design of control strategies that minimize the likelihood of selecting for cross-resistance, a whole transcriptomic approach was used to explore gene expression patterns in a multi-insecticide resistant population of Anopheles arabiensis from Oromia Region, Ethiopia. This field population was resistant to the diagnostic doses of malathion (average mortality of 71.9%) and permethrin (77.4%), with pools of survivors and unexposed individuals analyzed using Illumina RNA-sequencing, alongside insecticide susceptible reference strains. This population also demonstrated deltamethrin resistance but complete susceptibility to alpha-cypermethrin, bendiocarb and propoxur, providing a phenotypic basis for detecting insecticide-specific resistance mechanisms. Transcriptomic data revealed overexpression of genes including cytochrome P450s, glutathione-s-transferases and carboxylesterases (including CYP4C36, CYP6AA1, CYP6M2, CYP6M3, CYP6P4, CYP9K1, CYP9L1, GSTD3, GSTE2, GSTE3, GSTE4, GSTE5, GSTE7 and two carboxylesterases) that were shared between malathion and permethrin survivors. We also identified nineteen highly overexpressed cuticular-associated proteins (including CYP4G16, CYP4G17 and chitinase) and eighteen salivary gland proteins (including D7r4 short form salivary protein), which may be contributing to a non-specific resistance phenotype by either enhancing the cuticular barrier or promoting binding and sequestration of insecticides, respectively. These findings provide novel insights into the molecular basis of insecticide resistance in this lesser well-characterized major malaria vector species.

First report of the L1014F kdr mutation in wild populations of Anopheles arabiensis in Cabo Verde, West Africa

BACKGROUND

Due to the lack of vaccines, malaria control mainly involves the control of anopheline vectors (Anopheles spp.) using chemical insecticides. However, the prolonged and indiscriminate use of these compounds has led to the emergence of resistance in Anopheles populations in Africa. Insecticide resistance surveillance programs are less frequent in Cabo Verde than in other African countries. This study aimed to investigate the circulation of the L1014F and L1014S alleles in natural populations of Anopheles arabiensis collected from two sampling sites in the city of Praia, Cabo Verde.

METHODS

Anopheles larvae were collected from the two sampling sites and reared in the laboratory until the adult stage. Mosquitoes were first morphologically identified by classical taxonomy and then by molecular species identification using molecular markers. All Anopheles arabiensis were subjected to PCR analysis to screen for mutations associated to resistance in the Na_v gene.

RESULTS

A total of 105 mosquitoes, all belonging to the Anopheles gambiae complex, were identified by classical taxonomy as well as by molecular taxonomy. Molecular identification showed that 100% of the An. gambiae senso lato specimens analyzed corresponded to An. arabiensis. Analysis of the Na_v gene revealed the presence of L1014S and L1014F alleles with frequencies of 0.10 and 0.19, respectively.

CONCLUSIONS

Our data demonstrated, for the first time, the presence of the L1014F allele in the An. arabiensis population from Cabo Verde, as well as an increase in the frequency of the kdr L1014S allele reported in a previous study. The results of this study demonstrate the need to establish new approaches in vector control programs in Cabo Verde.

Highly sensitive droplet digital PCR-based diagnostics for the surveillance of malaria vector populations in low transmission and incipient resistance settings

BACKGROUND

Sensitive monitoring of Plasmodium infective mosquitoes in low malaria transmission settings is of high priority for disease control. Early detection of insecticide resistance at low frequencies is also key for vector monitoring nowadays, when new insecticides are launched to control vector populations.

RESEARCH DESIGN AND METHODS

An. gambiae mosquitoes with predetermined infection and resistance status were used to produce populations with various malaria infection rates and mutant allelic frequencies (MAFs) of target site insecticide resistance traits. Total RNA and gDNA were isolated and used in droplet Digital PCR (ddPCR) and Reverse Transcription (RT) ddPCR performed in the QX200 ddPCR System.

RESULTS

We developed a novel ddPCR for detecting P. falciparum DNA in pooled mosquito head-thoraces with infective rate as low as 1.0%. A dissection-free RT-ddPCR assay for specific infective-stage detection was additionally developed and validated (accuracy = 100%) in mosquito pools with infective rates down to 1.0%. A novel ddPCR assay for insecticide resistant alleles, which was able to reliably quantify MAFs as low as 0.050% in pooled mosquito specimens, is also reported.

CONCLUSIONS

We developed highly sensitive and efficient (RT-) ddPCR assays for contemporary operational needs that require monitoring of low malaria transmission and emerging insecticide resistance.

CRISPR/Cas9 modified An. gambiae carrying kdr mutation L1014F functionally validate its contribution in insecticide resistance and combined effect with metabolic enzymes

Insecticide resistance in Anopheles mosquitoes is a major obstacle in maintaining the momentum in reducing the malaria burden; mitigating strategies require improved understanding of the underlying mechanisms. Mutations in the target site of insecticides (the voltage gated sodium channel for the most widely used pyrethroid class) and over-expression of detoxification enzymes are commonly reported, but their relative contribution to phenotypic resistance remain poorly understood. Here we present a genome editing pipeline to introduce single nucleotide polymorphisms in An. gambiae which we have used to study the effect of the classical kdr mutation L1014F (L995F based on An. gambiae numbering), one of the most widely distributed resistance alleles. Introduction of 1014F in an otherwise fully susceptible genetic background increased levels of resistance to all tested pyrethroids and DDT ranging from 9.9-fold for permethrin to >24-fold for DDT. The introduction of the 1014F allele was sufficient to reduce mortality of mosquitoes after exposure to deltamethrin treated bednets, even as the only resistance mechanism present. When 1014F was combined with over-expression of glutathione transferase Gste2, resistance to permethrin increased further demonstrating the critical combined effect between target site resistance and detoxification enzymes in vivo. We also show that mosquitoes carrying the 1014F allele in homozygosity showed fitness disadvantages including increased mortality at the larval stage and a reduction in fecundity and adult longevity, which can have consequences for the strength of selection that will apply to this allele in the field.

Escalation of pyrethroid resistance in Anopheles mosquitoes threatens to reduce the effectiveness of our most important tools in malaria control. Studying the mechanisms underlying insecticide resistance is critical to design mitigation strategies. Here, using genome modified mosquitoes, we functionally characterize the most prevalent mutation in resistant mosquitoes, showing that it confers substantial levels of resistance to all tested pyrethroids and undermines the performance of pyrethroid-treated nets. Furthermore, we show that combining this mutation with elevated levels of a detoxification enzyme further increases resistance. The pipeline we have developed provides a robust approach to quantifying the contribution of different combinations of resistance mechanisms to the overall phenotype, providing the missing link between resistance monitoring and predictions of resistance impact.

Mutations associated with pyrethroid resistance in Varroa mite, a parasite of honey bees, are widespread across the United States

BACKGROUND

Managed honey bees are key pollinators of many crops and play an essential role in the United States food production. For more than ten years, beekeepers in the United States have been reporting high rates of colony losses. One of the drivers of these losses is the parasitic mite Varroa destructor. Maintaining healthy honey bee colonies in the United States is dependent on a successful control of this mite. The pyrethroid tau-fluvalinate (Apistan®) was among the first synthetic varroacides registered in the United States. With over 20 years of use, mites resistant to Apistan® have emerged, and so it is unsurprising that treatment failures have been reported. Resistance to tau-fluvalinate in US mite populations is associated with point mutations at position 925 of the voltage-gated sodium channel.

RESULTS

Here, we have generated a distribution map of pyrethroid resistance alleles in Varroa samples collected from US apiaries in 2016 and 2017, using a high throughput allelic discrimination assay based on TaqMan®. Our results evidence that knockdown resistance (kdr)-type mutations are widely distributed in Varroa populations across the country showing high variability among apiaries. We used these data to predict the phenotype of the mites in the case of treatments with pyrethroids.

CONCLUSION

We highlight the relevance of monitoring the resistance in mite populations to achieve an efficient control of this pest. We also put forward the benefits of implementing this methodology to provide data for designing pest management programs aiming to control Varroa. © 2021 Society of Chemical Industry.

First report of N1575Y mutation in Anopheles gambiae in Sierra Leone

The resistance of mosquito vectors to insecticides is one of the biological obstacles in the fight against malaria. Understanding of the status and mechanisms underlying the insecticide resistance in Anopheles gambiae species is necessary for success of vector control efforts. The study aimed to determine the molecular forms of An. gambiae from four districts in Sierra Leone during May and June 2018, and the level of N1575Y mutation. The molecular form identification of adult female An. gambiae mosquitoes reared from larvae were carried out using polymerase chain reaction and sequencing. And the N1575Y mutations were detected using SNaPshot and sequencing. As a result, significant differences were found in the distribution of An. gambiae molecular forms among regions (P < 0.001). And a total of 638 An. gambiae sensu stricto, 106 An. coluzzi, and 4 hybrid individuals were identified. Moreover, the overall N1575Y mutation frequency was 10.2% with no statistical difference among regions (χ² = 3.009, P = 0.390). In addition, no significant differences in N1575Y mutation frequency were found among different An. gambiae molecular forms (P = 0.383). In conclusion, the N1575Y mutation in An. gambiae populations in Sierra Leone was reported for the first time in the present study. It provides key evidence for the necessity of monitoring vector susceptibility levels to insecticides used in this country.

Identification of main malaria vectors and their insecticide resistance profile in internally displaced and indigenous communities in Eastern Democratic Republic of the Congo (DRC)

BACKGROUND

Malaria remains a major public health concern in the Democratic Republic of the Congo (DRC) and its control is affected by recurrent conflicts. Médecins Sans Frontières (MSF) initiated several studies to better understand the unprecedented incidence of malaria to effectively target and implement interventions in emergency settings. The current study evaluated the main vector species involved in malaria transmission and their resistance to insecticides, with the aim to propose the most effective tools and strategies for control of local malaria vectors.

METHODS

This study was performed in 52 households in Shamwana (Katanga, 2014), 168 households in Baraka (South Kivu, 2015) and 269 households in Kashuga (North Kivu, 2017). Anopheles vectors were collected and subjected to standardized Word Health Organization (WHO) and Center for Disease Control (CDC) insecticide susceptibility bioassays. Mosquito species determination was done using PCR and Plasmodium falciparum infection in mosquitoes was assessed by ELISA targeting circumsporozoite protein.

RESULTS

Of 3517 Anopheles spp. mosquitoes collected, Anopheles gambiae sensu lato (s.l.) (29.6%) and Anopheles funestus (69.1%) were the main malaria vectors. Plasmodium falciparum infection rates for An. gambiae s.l. were 1.0, 2.1 and 13.9% for Shamwana, Baraka and Kashuga, respectively. Anopheles funestus showed positivity rates of 1.6% in Shamwana and 4.4% in Baraka. No An. funestus were collected in Kashuga. Insecticide susceptibility tests showed resistance development towards pyrethroids in all locations. Exposure to bendiocarb, malathion and pirimiphos-methyl still resulted in high mosquito mortality.

CONCLUSIONS

This is one of only few studies from these conflict areas in DRC to report insecticide resistance in local malaria vectors. The data suggest that current malaria prevention methods in these populations are only partially effective, and require additional tools and strategies. Importantly, the results triggered MSF to consider the selection of a new insecticide for indoor residual spraying (IRS) and a new long-lasting insecticide-treated net (LLIN). The reinforcement of correct usage of LLINs and the introduction of targeted larviciding were also included as additional vector control tools as a result of the studies.

Protective effects of Olyset® Net on Plasmodium falciparum infection after three years of distribution in western Kenya

Background

Several types of insecticides, treating technologies and materials are available for long-lasting insecticide-treated nets (LLINs). The variations may result in different efficacies against mosquitoes and correspondingly infection risks for the Plasmodium falciparum malaria parasite. This cross-sectional study investigated whether infection risk varied among children who slept under different LLIN brands in rural villages of western Kenya.

Methods

Children sleeping under various types of LLINs were tested for P. falciparum infection using a diagnostic polymerase chain reaction (PCR) assay. Data were collected for other potential factors associated with infection risk: sleeping location (with bed/without bed), number of persons sharing the same net, dwelling wall material, gap of eaves (open/close), proportional hole index, socio-economic status, and density of indoor resting anophelines. Bed-net efficacy against the Anopheles gambiae susceptible strain was estimated using the WHO cone test and the tunnel test. The residual insecticide content on nets was measured.

Results

Seven LLIN brands were identified, and deltamethrin-based DawaPlus® 2.0 was the most popular (48%) followed by permethrin-based Olyset® Net (28%). The former LLIN was distributed in the area about six months before the present study was conducted, and the latter net was distributed at least three years before. Of 254 children analysed, P. falciparum PCR-positive prevalence was 58% for DawaPlus® 2.0 users and 38% for Olyset® users. The multiple regression analysis revealed that the difference was statistically significant (adjusted OR: 0.67, 95% credible interval: 0.45–0.97), whereas the confounders were not statistically important. Among randomly selected net samples, all DawaPlus® 2.0 (n = 20) and 95% of Olyset® (n = 19) passed either the cone test or the tunnel test.

Conclusions

Olyset® was more effective in reducing infection risk compared with DawaPlus® 2.0. Although the data from the present study were too limited to explain the mechanism clearly, the results suggest that the characteristics of the former brand are more suitable for the conditions, such as vector species composition, of the study area.

Fine scale spatial investigation of multiple insecticide resistance and underlying target-site and metabolic mechanisms in Anopheles gambiae in central Côte d'Ivoire

Routine monitoring of occurrence, levels and mechanisms of insecticide resistance informs effective management strategies, and should be used to assess the effect of new tools on resistance. As part of a cluster randomised controlled trial evaluating a novel insecticide-based intervention in central Côte d’Ivoire, we assessed resistance and its underlying mechanisms in Anopheles gambiae populations from a subset of trial villages. Resistance to multiple insecticides in An. gambiae s.s. and An. coluzzii was detected across villages, with dose–response assays demonstrating extremely high resistance intensity to the pyrethroid deltamethrin (> 1,500-fold), and mortality following exposure to pyrethroid-treated bednets was low (< 30% mortality in cone bioassays). The 1014F kdr mutation was almost fixed (≥ 90%) in all villages but the 1575Y kdr-amplifying mutation was relatively rare (< 15%). The carbamate and organophosphate resistance-associated Ace-1 G119S mutation was also detected at moderate frequencies (22–43%). Transcriptome analysis identified overexpression of P450 genes known to confer pyrethroid resistance (Cyp9K1, Cyp6P3, and Cyp6M2), and also a carboxylesterase (COEAE1F) as major candidates. Cyp6P3 expression was high but variable (up to 33-fold) and correlated positively with deltamethrin resistance intensity across villages (r² = 0.78, P = 0.02). Tools and strategies to mitigate the extreme and multiple resistance provided by these mechanisms are required in this area to avoid future control failures.

Genomic Variant Analyses in Pyrethroid Resistant and Susceptible Malaria Vector, Anopheles sinensis

Anopheles sinensis is a major malaria vector in Southeast Asia. Resistance to pyrethroid insecticides in this species has impeded malaria control in the region. Previous studies found that An. sinensis populations from Yunnan Province, China were highly resistant to deltamethrin and did not carry mutations in the voltage-gated sodium channel gene that cause knockdown resistance. In this study, we tested the hypothesis that other genomic variants are associated with the resistance phenotype. Using paired-end whole genome sequencing (DNA-seq), we generated 108 Gb of DNA sequence from deltamethrin -resistant and -susceptible mosquito pools with an average coverage of 83.3× depth. Using a stringent filtering method, we identified a total of 916,926 single nucleotide variants (SNVs), including 32,240 non-synonymous mutations. A total of 958 SNVs differed significantly in allele frequency between deltamethrin -resistant and -susceptible mosquitoes. Of these, 43 SNVs were present within 37 genes that code for immunity, detoxification, cuticular, and odorant proteins. A subset of 12 SNVs were randomly selected for genotyping of individual mosquitoes by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and showed consistent allele frequencies with the pooled DNA-seq derived allele frequencies. In addition, copy number variations (CNVs) were detected in 56 genes, including 33 that contained amplification alleles and 23 that contained deletion alleles in resistant mosquitoes compared to susceptible mosquitoes. The genomic variants described here provide a useful resource for future studies on the genetic mechanism of insecticide resistance in this important malaria vector species.

Asymptomatic Submicroscopic Plasmodium Infection Is Highly Prevalent and Is Associated with Anemia in Children Younger than 5 Years in South Kivu/Democratic Republic of Congo

One of the most important problems in controlling malaria is the limited access to effective and accurate diagnosis of malaria parasitemia. In the Democratic Republic of Congo (DRC), malaria is one of the leading causes of morbidity and mortality. The purpose of this study was to assess the prevalence of anemia and the relationship with asymptomatic submicroscopic Plasmodium infection. A cross-sectional study was carried out among 1,088 apparently healthy children aged between 6 and 59 months selected at random in the health zone of Miti Murhesa in South Kivu/DRC. Capillary blood was obtained for hemoglobin (Hb) concentration measurement by Hemocue^® Hb 301. Malaria detection was performed by microscopy and the loop-mediated isothermal amplification (LAMP) assay. Anemia was defined as Hb < 11g/dL. We applied the chi-square test for comparisons, and multiple logistic regression was used to identify the risk factors for anemia and submicroscopic Plasmodium infection. The prevalence of anemia was 39.6%, and the prevalence of parasitemia was 15.9% and 34.0% using microscopy and LAMP test, respectively. Submicroscopic Plasmodium infection was found in 22.3% of the children. The independent risk factors for anemia are Plasmodium infection, children younger than 24 months, low middle-upper arm circumference, and history of illness two weeks before. Otherwise, children with submicroscopic malaria infection have a significantly increased risk for anemia, with a need of transfusion. The prevalence of malaria infection was underestimated, when microscopy was used to diagnose malaria. Children with low parasitemia detected by LAMP but not by microscopy showed a significantly increased prevalence of anemia.

Mapping trends in insecticide resistance phenotypes in African malaria vectors

Mitigating the threat of insecticide resistance in African malaria vector populations requires comprehensive information about where resistance occurs, to what degree, and how this has changed over time. Estimating these trends is complicated by the sparse, heterogeneous distribution of observations of resistance phenotypes in field populations. We use 6,423 observations of the prevalence of resistance to the most important vector control insecticides to inform a Bayesian geostatistical ensemble modelling approach, generating fine-scale predictive maps of resistance phenotypes in mosquitoes from the Anopheles gambiae complex across Africa. Our models are informed by a suite of 111 predictor variables describing potential drivers of selection for resistance. Our maps show alarming increases in the prevalence of resistance to pyrethroids and DDT across sub-Saharan Africa from 2005 to 2017, with mean mortality following insecticide exposure declining from almost 100% to less than 30% in some areas, as well as substantial spatial variation in resistance trends.

Intensity of pyrethroid resistance in Anopheles gambiae before and after a mass distribution of insecticide-treated nets in Kinshasa and in 11 provinces of the Democratic Republic of Congo

BACKGROUND

Between 2011 and 2018, an estimated 134.8 million pyrethroid-treated long-lasting insecticidal nets (LLINs) were distributed nationwide in the Democratic Republic of Congo (DRC) for malaria control. Pyrethroid resistance has developed in DRC in recent years, but the intensity of resistance and impact on LLIN efficacy was not known. Therefore, the intensity of resistance of Anopheles gambiae sensu lato (s.l.) to permethrin and deltamethrin was monitored before and after a mass distribution of LLINs in Kinshasa in December 2016, and in 6 other sites across the country in 2017 and 11 sites in 2018.

METHODS

In Kinshasa, CDC bottle bioassays using 1, 2, 5, and 10 times the diagnostic dose of permethrin and deltamethrin were conducted using An. gambiae s.l. collected as larvae and reared to adults. Bioassays were conducted in four sites in Kinshasa province 6 months before a mass distribution of deltamethrin-treated LLINs and then two, six, and 10 months after the distribution. One site in neighbouring Kongo Central province was used as a control (no mass campaign of LLIN distribution during the study). Nationwide intensity assays were conducted in six sites in 2017 using CDC bottle bioassays and in 11 sites in 2018 using WHO intensity assays. A sub-sample of An. gambiae s.l. was tested by PCR to determine species composition and frequency of kdr-1014F and 1014S alleles.

RESULTS

In June 2016, before LLIN distribution, permethrin resistance intensity was high in Kinshasa; the mean mortality rate was 43% at the 5× concentration and 73% at the 10× concentration. Bioassays at 3 time points after LLIN distribution showed considerable variation by site and time and there was no consistent evidence for an increase in pyrethroid resistance intensity compared to the neighbouring control site. Tests of An. gambiae s.l. in 6 sites across the country in 2017 and 11 sites in 2018 showed all populations were resistant to the diagnostic doses of 3 pyrethroids. In 2018, the intensity of resistance varied by site, but was generally moderate for all three pyrethroids, with survivors at ×5 the diagnostic dose. Anopheles gambiae sensu stricto (s.s.) was the most common species identified across 11 sites in DRC, but in Kinshasa, An. gambiae s.s. (91%) and Anopheles coluzzii (8%) were sympatric.

CONCLUSIONS

Moderate or high intensity pyrethroid resistance was detected nationwide in DRC and is a serious threat to sustained malaria control with pyrethroid LLINs. Next generation nets (PBO nets or bi-treated nets) should be considered for mass distribution.

A high throughput multi-locus insecticide resistance marker panel for tracking resistance emergence and spread in Anopheles gambiae

The spread of resistance to insecticides in disease-carrying mosquitoes poses a threat to the effectiveness of control programmes, which rely largely on insecticide-based interventions. Monitoring mosquito populations is essential, but obtaining phenotypic measurements of resistance is laborious and error-prone. High-throughput genotyping offers the prospect of quick and repeatable estimates of resistance, while also allowing resistance markers to be tracked and studied. To demonstrate the potential of highly-mulitplexed genotypic screening for measuring resistance-association of mutations and tracking their spread, we developed a panel of 28 known or putative resistance markers in the major malaria vector Anopheles gambiae, which we used to screen mosquitoes from a wide swathe of Sub-Saharan Africa (Burkina Faso, Ghana, Democratic Republic of Congo (DRC) and Kenya). We found resistance association in four markers, including a novel mutation in the detoxification gene Gste2 (Gste2-119V). We also identified a duplication in Gste2 combining a resistance-associated mutation with its wild-type counterpart, potentially alleviating the costs of resistance. Finally, we describe the distribution of the multiple origins of kdr resistance, finding unprecedented diversity in the DRC. This panel represents the first step towards a quantitative genotypic model of insecticide resistance that can be used to predict resistance status in An. gambiae.

Development of molecular assays to detect target-site mechanisms associated with insecticide resistance in malaria vectors from Latin America

BACKGROUND

Malaria remains an important public health problem in Latin America, and the development of insecticide resistance in malaria vectors poses a major threat to malaria elimination efforts. Monitoring of insecticide susceptibility and the determination of the mechanisms involved in insecticide resistance are needed to effectively guide the deployment of appropriate vector control measures. Here, molecular assays have been developed to screen for mutations associated with insecticide resistance on the voltage-gated sodium channel (VGSC) and acetylcholinesterase-1 (Ace-1) genes in four malaria vectors from Latin America.

METHODS

Degenerate primers were designed to amplify a partial fragment on the VGSC and Ace-1 genes. Wild-caught individuals for Anopheles albimanus (also historical samples and individuals from a laboratory strain), Anopheles darlingi, Anopheles vestitipennis and Anopheles pseudopunctipennis were used to optimize the PCR assays. All samples were sequenced to validate the PCR results and DNA alignments were constructed for each gene using the unique haplotypes observed.

RESULTS

Primers designed successfully amplified the VGSC gene in An. albimanus, An. darlingi, An. vestitipennis and An. pseudopunctipennis, and the Ace-1 gene in both An. albimanus and An. darlingi. DNA sequencing revealed that compared with Anopheles gambiae, there were a total of 29, 28, 21 and 24 single nucleotide polymorphisms (SNPs) on the VGSC gene for An. albimanus (308 bp), An. darlingi (311 bp), An. pseudopunctipennis (263 bp) and An. vestitipennis (254 bp), respectively. On the 459 bp fragment of the Ace-1 gene, a total of 70 SNPs were detected in An. darlingi and 59 SNPs were detected in An. albimanus compared with An. gambiae. The SNPs detected on the VGSC gene were all synonymous. On the Ace-1 gene, non-synonymous substitutions were identified on three different codons. All species showed the homozygous wild-type kdr allele (coding for leucine) at codon 995 (formerly reported as codon 1014) on the VGSC gene, but one sample was heterozygous at codon 280 (formerly reported as codon 119) on the Ace-1 gene, coding for both the resistant (serine) and susceptible (glycine) amino acids.

CONCLUSIONS

New molecular assays to amplify and screen the regions of the VGSC and Ace-1 genes associated with insecticide resistance are reported for An. albimanus, An. darlingi, An. vestitipennis, and An. pseudopunctipennis. The development of these PCR assays presents an important advance in the analysis of target-site resistance in malaria vectors in the Americas, and will further facilitate the characterization of insecticide resistance mechanisms in these species.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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