Major African Malaria Vector Research Articles

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are associated with insecticide resistance in the major malaria vectors Anopheles gambiae s.l. and Anopheles funestus

Malaria remains one of the highest causes of morbidity and mortality, with 249 million cases and over 608,000 deaths in 2022. Insecticides, which target the Anopheles mosquito vector, are the primary method to control malaria. The widespread nature of resistance to the most important insecticide class, the pyrethroids, threatens the control of this disease. To reverse the stall in malaria control there is urgent need for new vector control tools, which necessitates understanding the molecular basis of pyrethroid resistance. In this study we utilised multi-omics data to identify uridine-diphosphate (UDP)-glycosyltransferases (UGTs) potentially involved in resistance across multiple Anopheles species. Phylogenetic analysis identifies sequence similarities between Anopheline UGTs and those involved in agricultural pesticide resistance to pyrethroids, pyrroles and spinosyns. Expression of five UGTs was characterised in An. gambiae and An. coluzzii to determine constitutive over-expression, induction, and tissue specificity. Furthermore, a UGT inhibitor, sulfinpyrazone, restored susceptibility to pyrethroids and DDT in An. gambiae, An. coluzzii, An. arabiensis and An. funestus, the major African malaria vectors. Taken together, this study provides clear association of UGTs with pyrethroid resistance as well as highlighting the potential use of sulfinpyrazone as a novel synergist for vector control.

Longitudinal survey of insecticide resistance in a village of central region of Burkina Faso reveals co-occurrence of 1014F, 1014S and 402L mutations in Anopheles coluzzii and Anopheles arabiensis

BackgroundPyrethroid resistance is one of the major threats for effectiveness of insecticide-treated bed nets (ITNs) in malaria vector control. Genotyping of mutations in the voltage gated sodium channel (VGSC) gene is widely used to easily assess the evolution and spread of pyrethroid target-site resistance among malaria vectors. L1014F and L1014S substitutions are the most common and best characterized VGSC mutations in major African malaria vector species of the Anopheles gambiae complex. Recently, an additional substitution involved in pyrethroid resistance, i.e. V402L, has been detected in Anopheles coluzzii from West Africa lacking any other resistance alleles at locus 1014. The evolution of target-site resistance mutations L1014F/S and V402L was monitored in An. coluzzii and Anopheles arabiensis specimens from a Burkina Faso village over a 10-year range after the massive ITN scale-up started in 2010.MethodsAnopheles coluzzii (N = 300) and An. arabiensis (N = 362) specimens collected both indoors and outdoors by different methods (pyrethrum spray catch, sticky resting box and human landing collections) in 2011, 2015 and 2020 at Goden village were genotyped by TaqMan assays and sequencing for the three target site resistance mutations; allele frequencies were statistically investigated over the years.ResultsA divergent trend in resistant allele frequencies was observed in the two species: 1014F decreased in An. coluzzii (from 0.76 to 0.52) but increased in An. arabiensis (from 0.18 to 0.70); 1014S occurred only in An. arabiensis and slightly decreased over time (from 0.33 to 0.23); 402L increased in An. coluzzii (from 0.15 to 0.48) and was found for the first time in one An. arabiensis specimen. In 2020 the co-occurrence of different resistance alleles reached 43% in An. coluzzii (alleles 410L and 1014F) and 32% in An. arabiensis (alleles 1014F and 1014S).ConclusionsOverall, an increasing level of target-site resistance was observed among the populations with only 1% of the two malaria vector species being wild type at both loci, 1014 and 402, in 2020. This, together with the co-occurrence of different mutations in the same specimens, calls for future investigations on the possible synergism between resistance alleles and their phenotype to implement local tailored intervention strategies.

Discovery of knock-down resistance in the major African malaria vector Anopheles funestus.

A major mechanism of insecticide resistance in insect pests is knock-down resistance (kdr) caused by mutations in the voltage-gated sodium channel (Vgsc) gene. Despite being common in most malaria Anopheles vector species, kdr mutations have never been observed in Anopheles funestus, the principal malaria vector in Eastern and Southern Africa. While monitoring 10 populations of An. funestus in Tanzania, we unexpectedly found resistance to DDT, a banned insecticide, in one location. Through whole-genome sequencing of 333 An. funestus samples from these populations, we found 8 novel amino acid substitutions in the Vgsc gene, including the kdr variant, L976F (L1014F in An. gambiae), in tight linkage disequilibrium with another (P1842S). The mutants were found only at high frequency in one region, with a significant decline between 2017 and 2023. Notably, kdr L976F was strongly associated with survivorship to the exposure to DDT insecticide, while no clear association was noted with a pyrethroid insecticide (deltamethrin). Further study is necessary to identify the origin and spread of kdr in An. funestus, and the potential threat to current insecticide-based vector control in Africa.

Characterization of the Tissue and Strain-Specific Microbiota of Anopheles funestus Giles (Diptera: Culicidae).

The mosquito microbiota is a critical determinant of mosquito life history. It is therefore a target for novel vector control strategies like paratransgenesis. However, the microbiota in Anopheles funestus, a major African malaria vector, is poorly characterized. Thus, the study aimed to investigate the overall bacterial landscape in the salivary glands, ovaries and midguts of three laboratory strains of An. funestus differing in insecticide-resistant phenotype by sequencing the V3-V4 hypervariable region of bacterial 16S rRNA genes. When examining alpha diversity, the salivary glands harbored significantly more bacteria in terms of species richness and evenness compared to ovaries and midguts. On the strain level, the insecticide-susceptible FANG strain had significantly lower bacterial diversity than the insecticide-resistant FUMOZ and FUMOZ-R strains. When looking at beta diversity, the compositions of microbiota between the three tissues as well as between the strains were statistically different. While there were common bacteria across all three tissues and strains of interest, each tissue and strain did exhibit differentially abundant bacterial genera. However, overall, the top five most abundant genera across all tissues and strains were Elizabethkingia, Acinetobacter, Aeromonas, Cedecea and Yersinia. The presence of shared microbiota suggests a core microbiota that could be exploited for paratransgenesis efforts.

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

BackgroundThe 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.MethodsA 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.ResultsIn 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.ConclusionsDespite 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.

The insecticidal activity of essential oil constituents against pyrethroid-resistant Anopheles funestus (Diptera: Culicidae).

Malaria vector control relies on the use of insecticides for indoor residual spraying and long-lasting bed nets. However, insecticide resistance to pyrethroids among others, has escalated. Anopheles funestus, one of the major African malaria vectors, has attained significant levels of resistance to pyrethroids. Overexpressed P450 monooxygenases have been previously identified in pyrethroid resistant An. funestus. The escalating resistance against conventional insecticides signals an urgent need for identification of novel insecticides. Essential oils have gained recognition as promising sources of alternative natural insecticides. This study investigated six essential oil constituents, farnesol, (-)-α-bisabolol, cis-nerolidol, trans-nerolidol, methyleugenol, santalol (α and β isomers) and essential oil of sandalwood, for the adulticidal effects against pyrethroid-resistant An. funestus strain. The susceptibility against these terpenoids were evaluated on both pyrethroid-susceptible and resistant An. funestus. Furthermore, the presence of overexpressed monooxygenases in resistant An. funestus was confirmed. Results showed that both the pyrethroid-susceptible and resistant An. funestus were susceptible to three EOCs; cis-nerolidol, trans-nerolidol and methyleugenol. On the other hand, the pyrethroid-resistant An. funestus survived exposure to both farnesol and (-)-α-bisabolol. This study however does not show any direct association of the overexpressed Anopheles monooxygenases and the efficacy of farnesol and (-)-α-bisabolol. The enhanced activity of these terpenoids against resistant An. funestus that has been pre-exposed to a synergist, piperonyl butoxide, suggests their potential effectiveness in combination with monooxygenase inhibitors. This study proposes that cis-nerolidol, trans-nerolidol and methyleugenol are potential agents for further investigation as novel bioinsecticides against pyrethroid-resistant An. funestus strain.

Influence of GST- and P450-based metabolic resistance to pyrethroids on blood feeding in the major African malaria vector Anopheles funestus (lay summary)

Influence of GST- and P450-based metabolic resistance to pyrethroids on blood feeding in the major African malaria vector Anopheles funestus (lay summary)

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.

Manure and mosquitoes: life history traits of two malaria vector species enhanced by larval exposure to cow dung, whilst chicken dung has a strong negative effect

BackgroundMalaria vectors have a strong ecological association with rice agroecosystems, which can provide abundant aquatic habitats for larval development. Climate-adapted rice cultivation practices, such as the System of Rice Intensification (SRI), are gaining popularity in malaria-endemic countries seeking to expand rice production; however, the potential impact of these practices on vector populations has not been well characterised. In particular, SRI encourages the use of organic fertilisers (OFs), such as animal manures, as low-cost and environmentally friendly alternatives to industrially produced inorganic fertilisers. We therefore set out to understand the effects of two common manure-based OFs on the life history traits of two major African malaria vectors, Anopheles arabiensis and Anopheles gambiae sensu stricto (s.s.).MethodsLarvae of An. arabiensis and An. gambiae s.s. were reared from first instar to emergence in water containing either cow or chicken dung at one of four concentrations (0.25, 0.5, 0.75, and 1.0 g/100 ml), or in a clean water control. Their life history traits were recorded, including survival, development rate, adult production, and adult wing length.ResultsExposure to cow dung significantly increased the development rate of An. gambiae s.s. independent of concentration, but did not affect the overall survival and adult production of either species. Chicken dung, however, significantly reduced survival and adult production in both species, with a greater effect as concentration increased. Interestingly, An. arabiensis exhibited a relative tolerance to the lowest chicken dung concentration, in that survival was unaffected and adult production was not reduced to the same extent as in An. gambiae s.s. The effects of chicken dung on development rate were less clear in both species owing to high larval mortality overall, though there was some indication that it may reduce development rate. Adult wing lengths in males and females increased with higher concentrations of both cow and chicken dung.ConclusionsOur findings suggest that manure-based OFs significantly alter the life history traits of An. gambiae s.s. and An. arabiensis. In both species, exposure to cow dung may improve fitness, whereas exposure to chicken dung may reduce it. These findings have implications for understanding vector population dynamics in rice agroecosystems and may inform the use of OFs in SRI, and rice agriculture more widely, to avoid their adverse effects in enhancing vector fitness.Graphical

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii.

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences. Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Harmonic convergence coordinates swarm mating by enhancing mate detection in the malaria mosquito Anopheles gambiae

The mosquito Anopheles gambiae is a major African malaria vector, transmitting parasites responsible for significant mortality and disease burden. Although flight acoustics are essential to mosquito mating and present promising alternatives to insecticide-based vector control strategies, there is limited data on mosquito flight tones during swarming. Here, for the first time, we present detailed analyses of free-flying male and female An. gambiae flight tones and their harmonization (harmonic convergence) over a complete swarm sequence. Audio analysis of single-sex swarms showed synchronized elevation of male and female flight tones during swarming. Analysis of mixed-sex swarms revealed additional 50 Hz increases in male and female flight tones due to mating activity. Furthermore, harmonic differences between male and female swarm tones in mixed-sex swarms and in single-sex male swarms with artificial female swarm audio playback indicate that frequency differences of approximately 50 Hz or less at the male second and female third harmonics (M2:F3) are maintained both before and during mating interactions. This harmonization likely coordinates male scramble competition by maintaining ideal acoustic recognition within mating pairs while acoustically masking phonotactic responses of nearby swarming males to mating females. These findings advance our knowledge of mosquito swarm acoustics and provide vital information for reproductive control strategies.

Combined over-expression of two cytochrome P450 genes exacerbates the fitness cost of pyrethroid resistance in the major African malaria vector Anopheles funestus

Metabolic resistance driven by multiple P450 genes is worsening insecticide resistance in malaria vectors. However, it remains unclear whether such multiple over-expression imposes an additive fitness cost in the vectors. Here, we showed that two highly over-expressed P450 genes (CYP6P9a and CYP6P9b) combine to impose additive fitness costs in pyrethroid-resistant Anopheles funestus. Genotyping of the CYP6P9b resistance allele in hybrid mosquitoes from a pyrethroid-resistant FUMOZ-R and the susceptible FANG strains revealed that this gene imposes a fitness cost in resistant mosquitoes similar to CYP6P9a. Homozygote susceptible CYP6P9b_S (SS) significantly lay more eggs than the resistant (OR = 2.2, P = 0.04) and with greater hatching rate (p < 0.04). Homozygote resistant larvae CYP6P9b_R (RR) developed significantly slower than homozygote susceptible from L1-L4 (χ2 = 7.2; P = 0.007) with a late pupation observed for RR compared to both heterozygotes and homozygotes susceptible (χ2 = 11.17; P = 0.0008). No difference was observed between genotypes for adult longevity with no change in allele frequency and gene expression across the lifespan. Furthermore, we established that CYP6P9b combines with CYP6P9a to additively exacerbate the fitness cost of pyrethroid resistance with a greater reduction in fecundity/fertility and increased developmental time of double homozygote resistant mosquitoes. Moreover, an increased proportion of double homozygote susceptible individuals was noted over 10 generations in the insecticide-free environment (χ2 = 6.3; P = 0.01) suggesting a reversal to susceptibility in the absence of selection. Such greater fitness cost imposed by multiple P450 genes shows that resistance management strategy based on rotation could help slow the spread of resistance.

Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii.

Mosquito control remains a central pillar of efforts to reduce malaria burden in sub-Saharan Africa. However, insecticide resistance is entrenched in malaria vector populations, and countries with a high malaria burden face a daunting challenge to sustain malaria control with a limited set of surveillance and intervention tools. Here we report on the second phase of a project to build an open resource of high-quality data on genome variation among natural populations of the major African malaria vector species Anopheles gambiae and Anopheles coluzzii. We analyzed whole genomes of 1142 individual mosquitoes sampled from the wild in 13 African countries, as well as a further 234 individuals comprising parents and progeny of 11 laboratory crosses. The data resource includes high-confidence single-nucleotide polymorphism (SNP) calls at 57 million variable sites, genome-wide copy number variation (CNV) calls, and haplotypes phased at biallelic SNPs. We use these data to analyze genetic population structure and characterize genetic diversity within and between populations. We illustrate the utility of these data by investigating species differences in isolation by distance, genetic variation within proposed gene drive target sequences, and patterns of resistance to pyrethroid insecticides. This data resource provides a foundation for developing new operational systems for molecular surveillance and for accelerating research and development of new vector control tools. It also provides a unique resource for the study of population genomics and evolutionary biology in eukaryotic species with high levels of genetic diversity under strong anthropogenic evolutionary pressures.

Influence of GST- and P450-based metabolic resistance to pyrethroids on blood feeding in the major African malaria vector Anopheles funestus.

Insecticide resistance genes are often associated with pleiotropic effects on various mosquito life-history traits. However, very little information is available on the impact of insecticide resistance on blood feeding process in mosquitoes. Here, using two recently detected DNA-based metabolic markers in the major malaria vector, An. funestus, we investigated how metabolic resistance genes could affect the blood meal intake. After allowing both the field F1 and lab F8 Anopheles funestus strains to feed on the human arm for 30 minutes, we assessed the association between key parameters of blood meal process including, probing time, feeding duration, blood feeding success, blood meal size, and markers of glutathione S-transferase (L119F-GSTe2) and cytochrome P450 (CYP6P9a_R)-mediated metabolic resistance. None of the parameters of blood meal process was associated with L119F-GSTe2 genotypes. By contrast, for CYP6P9a_R, homozygous resistant mosquitoes were significantly more able to blood-feed than homozygous susceptible (OR = 3.3; CI 95%: 1.4-7.7; P = 0.01) mosquitoes. Moreover, the volume of blood meal ingested by CYP6P9a-SS mosquitoes was lower than that of CYP6P9a-RS (P<0.004) and of CYP6P9a-RR (P<0.006). This suggests that CYP6P9a gene is inked with the feeding success and blood meal size of An. funestus. However, no correlation was found in the expression of CYP6P9a and that of genes encoding for salivary proteins involved in blood meal process. This study suggests that P450-based metabolic resistance may influence the blood feeding process of Anopheles funestus mosquito and consequently its ability to transmit malaria parasites.

Cytochrome P450 metabolic resistance (CYP6P9a) to pyrethroids imposes a fitness cost in the major African malaria vector Anopheles funestus

Metabolic resistance threatens the sustainability of pyrethroid-based malaria control interventions. Elucidating the fitness cost and potential reversal of metabolic resistance is crucial to design suitable resistance management strategies. Here, we deciphered the fitness cost associated with the CYP6P9a (P450-mediated metabolic resistance) in the major African malaria vector Anopheles funestus. Reciprocal crosses were performed between a pyrethroid susceptible (FANG) and resistant (FUMOZ-R) laboratory strains and the hybrid strains showed intermediate resistance. Genotyping the CYP6P9a-R resistance allele in oviposited females revealed that CYP6P9a negatively impacts the fecundity as homozygote susceptible mosquitoes (CYP6P9a-SS) lay more eggs than heterozygote (OR = 2.04: P = 0.01) and homozygote resistant mosquitoes. CYP6P9a also imposes a significant fitness cost on the larval development as homozygote resistant larvae (CYP6P9a-RR) developed significantly slower than heterozygote and homozygote susceptible mosquitoes (χ2 = 11.2; P = 0.0008). This fitness cost was further supported by the late pupation of homozygote resistant than susceptible mosquitoes (OR = 2.50; P < 0.01). However, CYP6P9a does not impact the longevity as no difference was observed in the life span of mosquitoes with different genotypes (χ2 = 1.6; P = 0.9). In this hybrid strain, a significant decrease of the resistant CYP6P9a-RR genotype was observed after ten generations (χ2 = 6.6; P = 0.01) suggesting a reversal of P450-based resistance in the absence of selection. This study shows that the P450-mediated metabolic resistance imposes a high fitness cost in malaria vectors supporting that a resistance management strategy based on rotation could help mitigate the impact of such resistance.

Evidence of extrinsic factors dominating intrinsic blood host preferences of major African malaria vectors

One of the key determinants of a haematophagous vector’s capacity to transmit pathogens is its selection of which host to secure a blood meal from. This choice is influenced by both intrinsic (genetic) and extrinsic (environmental) factors, but little is known of their relative contributions. Blood fed Anopheles mosquitoes were collected from a malaria endemic village in Ghana. Collections were conducted across a range of different host availabilities and from both indoor and outdoor locations. These environmental factors were shown to impact dramatically the host choice of caught malaria vectors: mosquitoes caught indoors were ten-fold more likely to have sourced their blood meal from humans; and a halving in odds of being human-fed was found for mosquitoes caught only 25 m from the centre of the village. For the first time, we demonstrate that anthropophagy was better explained by extrinsic factors (namely, local host availability and indoor/outdoor trapping location) than intrinsic factors (namely, the (sibling) species of the mosquito caught) (respective Akaike information criterion estimates: 243.0 versus 359.8). Instead of characterizing biting behaviour on a taxonomic level, we illustrate the importance of assessing local entomology. Accounting for this behavioural plasticity is important, both in terms of measuring effectiveness of control programmes and in informing optimal disease control strategies.

Functional genetic validation of key genes conferring insecticide resistance in the major African malaria vector, Anopheles gambiae

Resistance in Anopheles gambiae to members of all 4 major classes (pyrethroids, carbamates, organochlorines, and organophosphates) of public health insecticides limits effective control of malaria transmission in Africa. Increase in expression of detoxifying enzymes has been associated with insecticide resistance, but their direct functional validation in An. gambiae is still lacking. Here, we perform transgenic analysis using the GAL4/UAS system to examine insecticide resistance phenotypes conferred by increased expression of the 3 genes-Cyp6m2, Cyp6p3, and Gste2-most often found up-regulated in resistant An. gambiae We report evidence in An. gambiae that organophosphate and organochlorine resistance is conferred by overexpression of GSTE2 in a broad tissue profile. Pyrethroid and carbamate resistance is bestowed by similar Cyp6p3 overexpression, and Cyp6m2 confers only pyrethroid resistance when overexpressed in the same tissues. Conversely, such Cyp6m2 overexpression increases susceptibility to the organophosphate malathion, presumably due to conversion to the more toxic metabolite, malaoxon. No resistant phenotypes are conferred when either Cyp6 gene overexpression is restricted to the midgut or oenocytes, indicating that neither tissue is involved in insecticide resistance mediated by the candidate P450s examined. Validation of genes conferring resistance provides markers to guide control strategies, and the observed negative cross-resistance due to Cyp6m2 gives credence to proposed dual-insecticide strategies to overcome pyrethroid resistance. These transgenic An. gambiae-resistant lines are being used to test the "resistance-breaking" efficacy of active compounds early in their development.

A chromosome-scale assembly of the major African malaria vector Anopheles funestus.

Background Anopheles funestus is one of the 3 most consequential and widespread vectors of human malaria in tropical Africa. However, the lack of a high-quality reference genome has hindered the association of phenotypic traits with their genetic basis in this important mosquito.FindingsHere we present a new high-quality A. funestus reference genome (AfunF3) assembled using 240× coverage of long-read single-molecule sequencing for contigging, combined with 100× coverage of short-read Hi-C data for chromosome scaffolding. The assembled contigs total 446 Mbp of sequence and contain substantial duplication due to alternative alleles present in the sequenced pool of mosquitos from the FUMOZ colony. Using alignment and depth-of-coverage information, these contigs were deduplicated to a 211 Mbp primary assembly, which is closer to the expected haploid genome size of 250 Mbp. This primary assembly consists of 1,053 contigs organized into 3 chromosome-scale scaffolds with an N50 contig size of 632 kbp and an N50 scaffold size of 93.811 Mbp, representing a 100-fold improvement in continuity versus the current reference assembly, AfunF1.ConclusionThis highly contiguous and complete A. funestus reference genome assembly will serve as an improved basis for future studies of genomic variation and organization in this important disease vector.

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

Metabolic resistance to insecticides such as pyrethroids in mosquito vectors threatens control of malaria in Africa. Unless it is managed, recent gains in reducing malaria transmission could be lost. To improve monitoring and assess the impact of insecticide resistance on malaria control interventions, we elucidated the molecular basis of pyrethroid resistance in the major African malaria vector, Anopheles funestus We showed that a single cytochrome P450 allele (CYP6P9a_R) in A. funestus reduced the efficacy of insecticide-treated bednets for preventing transmission of malaria in southern Africa. Expression of key insecticide resistance genes was detected in populations of this mosquito vector throughout Africa but varied according to the region. Signatures of selection and adaptive evolutionary traits including structural polymorphisms and cis-regulatory transcription factor binding sites were detected with evidence of selection due to the scale-up of insecticide-treated bednet use. A cis-regulatory polymorphism driving the overexpression of the major resistance gene CYP6P9a allowed us to design a DNA-based assay for cytochrome P450-mediated resistance to pyrethroid insecticides. Using this assay, we tracked the spread of pyrethroid resistance and found that it was almost fixed in mosquitoes from southern Africa but was absent from mosquitoes collected elsewhere in Africa. Furthermore, a field study in experimental huts in Cameroon demonstrated that mosquitoes carrying the resistance CYP6P9a_R allele survived and succeeded in blood feeding more often than did mosquitoes that lacked this allele. Our findings highlight the need to introduce a new generation of insecticide-treated bednets for malaria control that do not rely on pyrethroid insecticides.

Fitness Costs of the Glutathione S-Transferase Epsilon 2 (L119F-GSTe2) Mediated Metabolic Resistance to Insecticides in the Major African Malaria Vector Anopheles Funestus.

Metabolic resistance to insecticides threatens malaria control. However, little is known about its fitness cost in field populations of malaria vectors, thus limiting the design of suitable resistance management strategies. Here, we assessed the association between the glutathione S-transferase GSTe2-mediated metabolic resistance and life-traits of natural populations of Anopheles funestus. A total of 1200 indoor resting blood-fed female An. funestus (F0) were collected in Mibellon, Cameroon (2016/2017), and allowed to lay eggs individually. Genotyping of F1 mosquitoes for the L119F-GSTE2 mutation revealed that L/L119-homozygote susceptible (SS) mosquitoes significantly laid more eggs than heterozygotes L119F-RS (odds ratio (OR) = 2.06; p < 0.0001) and homozygote resistant 119F/F-RR (OR = 2.93; p < 0.0001). L/L119-SS susceptible mosquitoes also showed the higher ability for oviposition than 119F/F-RR resistant (OR = 2.68; p = 0.0002) indicating a reduced fecundity in resistant mosquitoes. Furthermore, L119F-RS larvae developed faster (nine days) than L119F-RR and L119F-SS (11 days) (X2 = 11.052; degree of freedom (df) = 4; p = 0.02) suggesting a heterozygote advantage effect for larval development. Interestingly, L/L119-SS developed faster than 119F/F-RR (OR = 5.3; p < 0.0001) revealing an increased developmental time in resistant mosquitoes. However, genotyping and sequencing revealed that L119F-RR mosquitoes exhibited a higher adult longevity compared to RS (OR > 2.2; p < 0.05) and SS (OR > 2.1; p < 0.05) with an increased frequency of GSTe2-resistant haplotypes in mosquitoes of D30 after adult emergence. Additionally, comparison of the expression of GSTe2 revealed a significantly increased expression from D1-D30 after emergence of adults (Anova test (F) = 8; df= 3; p = 0.008). The negative association between GSTe2 and some life traits of An. funestus could facilitate new resistance management strategies. However, the increased longevity of GSTe2-resistant mosquitoes suggests that an increase in resistance could exacerbate malaria transmission.