The interplay between malaria vectors and human activity accounts for high residual malaria transmission in a Burkina Faso village with universal ITN coverage

BACKGROUND

Mosquito and human behaviour interaction is a key determinant of the maximum level of protection against malaria that can be provided by insecticide-treated nets (ITNs). Nevertheless, scant literature focuses on this interaction, overlooking a fundamental factor for efficient malaria control. This study aims to estimate malaria transmission risk in a Burkina Faso village by integrating vector biting rhythms with some key information about human habits.

METHODS

Indoor/outdoor human landing catches were conducted for 16 h (16:00-08:00) during 8 nights (September 2020) in Goden village. A survey about net usage and sleeping patterns was submitted to half the households (October-December 2020). A subsample of collected specimens of Anopheles gambiae sensu lato was molecularly processed for species identification, Plasmodium detection from heads-thoraxes and L1014F pyrethroid-resistance allele genotyping. Hourly mosquito abundance was statistically assessed by GLM/GAM, and the entomological inoculation rate (EIR) was corrected for the actual ITN usage retrieved from the questionnaire.

RESULTS

Malaria transmission was mainly driven by Anopheles coluzzii (68.7%) followed by A. arabiensis (26.2%). The overall sporozoite rate was 2% with L1014F estimated frequency of 0.68 (N = 1070 out of 15,201 A. gambiae s.l. collected). No major shift in mosquito biting rhythms in response to ITN or differences between indoor and outdoor catches were detected. Impressive high biting pressure (mean 30.3 mosquitoes/person/hour) was exerted from 20:00 to 06:00 with a peak at 4:00. Human survey revealed that nearly all inhabitants were awake before 20:00 and after 7:00 and at least 8.7% had no access to bednets. Adjusting for anthropological data, the EIR dropped from 6.7 to 1.2 infective bites/person/16 h. In a scenario of full net coverage and accounting only for the human sleeping patterns, the daily malaria transmission risk not targetable by ITNs was 0.69 infective bites.

CONCLUSIONS

The high mosquito densities and interplay between human/vector activities means that an estimated 10% of residual malaria transmission cannot be prevented by ITNs in the village. Locally tailored studies, like the current one, are essential to explore the heterogeneity of human exposure to infective bites and, consequently, to instruct the adoption of new vector control tools strengthening individual and community protection.

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.

Pyriproxyfen-treated bed nets reduce reproductive fitness and longevity of pyrethroid-resistant Anopheles gambiae under laboratory and field conditions

Background

The efficacy of insecticide-treated nets (ITNs) containing the insect growth regulator pyriproxyfen (PPF) and pyrethroid insecticides (PPF-ITNs) is being assessed in clinical trials to determine whether they provide greater protection from malaria than standard pyrethroid-treated ITNs in areas where mosquitoes are resistant to pyrethroids. Understanding the entomological mode of action of this new ITN class will aide interpretation of the results from these trials.

Methods

Anopheles gambiae sensu lato (s.l.) mosquitoes from a susceptible laboratory strain were exposed to PPF-treated netting 24 h, 6 h, and immediately prior to, or 24 h post blood feeding, and the impact on fecundity, fertility and longevity recorded. Pyrethroid-resistant populations were exposed to nets containing permethrin and PPF (PPF-ITNs) in cone bioassays and daily mortality recorded. Mosquitoes were also collected from inside houses pre- and post-distribution of PPF-ITNs in a clinical trial conduced in Burkina Faso; female An. gambiae s.l. were then assessed for fecundity and fertility.

Results

PPF exposure reduced the median adult lifespan of insecticide-susceptible mosquitoes by 4 to 5 days in all exposure times (p < 0.05) other than 6 h pre-blood meal and resulted in almost complete lifelong sterilization. The longevity of pyrethroid-resistant mosquitoes was also reduced by at least 5 days after exposure to PPF-ITNs compared to untreated nets, but was unaffected by exposure to standard pyrethroid only ITNs. A total of 386 blood-fed or gravid An. gambiae s.l. females were collected from five villages between 1 and 12 months before distribution of PPF-ITNs. Of these mosquitoes, 75% laid eggs and the remaining 25% appeared to have normal ovaries upon dissection. In contrast, only 8.6% of the 631 blood-fed or gravid An. gambiae s.l. collected post PPF-ITN distribution successfully oviposited; 276 (43.7%) did not oviposit but had apparently normal ovaries upon dissection, and 301 (47.7%) did not oviposit and had abnormal eggs upon dissection. Egg numbers were also significantly lower (average of 138/female prior distribution vs 85 post distribution, p < 0.05).

Conclusion

Exposure to a mixture of PPF and pyrethroids on netting shortens the lifespan of mosquitoes and reduces reproductive output. Sterilization of vectors lasted at least one year under operational conditions. These findings suggest a longer effective lifespan of PPF-pyrethroid nets than reported previously.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03794-z.

Comprehensive Ecological and Geographic Characterization of Eukaryotic and Prokaryotic Microbiomes in African Anopheles

Exposure of mosquitoes to numerous eukaryotic and prokaryotic microbes in their associated microbiomes has probably helped drive the evolution of the innate immune system. To our knowledge, a metagenomic catalog of the eukaryotic microbiome has not been reported from any insect. Here we employ a novel approach to preferentially deplete host 18S ribosomal RNA gene amplicons to reveal the composition of the eukaryotic microbial communities of Anopheles larvae sampled in Kenya, Burkina Faso and Republic of Guinea (Conakry). We identified 453 eukaryotic operational taxonomic units (OTUs) associated with Anopheles larvae in nature, but an average of 45% of the 18S rRNA sequences clustered into OTUs that lacked a taxonomic assignment in the Silva database. Thus, the Anopheles microbiome contains a striking proportion of novel eukaryotic taxa. Using sequence similarity matching and de novo phylogenetic placement, the fraction of unassigned sequences was reduced to an average of 4%, and many unclassified OTUs were assigned as relatives of known taxa. A novel taxon of the genus Ophryocystis in the phylum Apicomplexa (which also includes Plasmodium) is widespread in Anopheles larvae from East and West Africa. Notably, Ophryocystis is present at fluctuating abundance among larval breeding sites, consistent with the expected pattern of an epidemic pathogen. Species richness of the eukaryotic microbiome was not significantly different across sites from East to West Africa, while species richness of the prokaryotic microbiome was significantly lower in West Africa. Laboratory colonies of Anopheles coluzzii harbor 26 eukaryotic OTUs, of which 38% (n = 10) are shared with wild populations, while 16 OTUs are unique to the laboratory colonies. Genetically distinct An. coluzzii colonies co-housed in the same facility maintain different prokaryotic microbiome profiles, suggesting a persistent host genetic influence on microbiome composition. These results provide a foundation to understand the role of the Anopheles eukaryotic microbiome in vector immunity and pathogen transmission. We hypothesize that prevalent apicomplexans such as Ophryocystis associated with Anopheles could induce interference or competition against Plasmodium within the vector. This and other members of the eukaryotic microbiome may offer candidates for new vector control tools.

Leucine-Rich Immune Factor APL1 Is Associated With Specific Modulation of Enteric Microbiome Taxa in the Asian Malaria Mosquito Anopheles stephensi

The commensal gut microbiome is contained by the enteric epithelial barrier, but little is known about the degree of specificity of host immune barrier interactions for particular bacterial taxa. Here, we show that depletion of leucine-rich repeat immune factor APL1 in the Asian malaria mosquito Anopheles stephensi is associated with higher midgut abundance of just the family Enterobacteraceae, and not generalized dysbiosis of the microbiome. The effect is explained by the response of a narrow clade containing two main taxa related to Klebsiella and Cedecea. Analysis of field samples indicate that these two taxa are recurrent members of the wild Anopheles microbiome. Triangulation using sequence and functional data incriminated relatives of C. neteri and Cedecea NFIX57 as candidates for the Cedecea component, and K. michiganensis, K. oxytoca, and K.sp. LTGPAF-6F as candidates for the Klebsiella component. APL1 presence is associated with host ability to specifically constrain the abundance of a narrow microbiome clade of the Enterobacteraceae, and the immune factor may promote homeostasis of this clade in the enteric microbiome for host benefit.

Influence of genetic polymorphism on transcriptional enhancer activity in the malaria vector Anopheles coluzzii

Enhancers are cis-regulatory elements that control most of the developmental and spatial gene expression in eukaryotes. Genetic variation of enhancer sequences is known to influence phenotypes, but the effect of enhancer variation upon enhancer functional activity and downstream phenotypes has barely been examined in any species. In the African malaria vector, Anopheles coluzzii, we identified candidate enhancers in the proximity of genes relevant for immunity, insecticide resistance, and development. The candidate enhancers were functionally validated using luciferase reporter assays, and their activity was found to be essentially independent of their physical orientation, a typical property of enhancers. All of the enhancers segregated genetically polymorphic alleles, which displayed significantly different levels of functional activity. Deletion mutagenesis and functional testing revealed a fine structure of positive and negative regulatory elements that modulate activity of the enhancer core. Enhancer polymorphisms segregate in wild A. coluzzii populations in West Africa. Thus, enhancer variants that modify target gene expression leading to likely phenotypic consequences are frequent in nature. These results demonstrate the existence of naturally polymorphic A. coluzzii enhancers, which may help explain important differences between individuals or populations for malaria transmission efficiency and vector adaptation to the environment.

A novel nested polymerase chain reaction assay targeting Plasmodium mitochondrial DNA in field-collected Anopheles mosquitoes

Sensitive techniques for the detection of Plasmodium (Aconoidasida: Plasmodiidae) sporozoites in field-collected malaria vectors are essential for the correct assessment of risk for malaria transmission. A real-time polymerase chain reaction (RT-PCR) protocol targeting Plasmodium mtDNA proved to be much more sensitive in detecting sporozoites in mosquitoes than the widely used enzyme-linked immunosorbent assay targeting Plasmodium circumsporozoite protein (CSP-ELISA). However, because of the relatively high costs associated with equipment and reagents, RT-PCRs are mostly used to assess the outcomes of experimental infections in the frame of research experiments, rather than in routine monitoring of mosquito infection in the field. The present authors developed a novel mtDNA-based nested PCR protocol, modified from a loop-mediated isothermal amplification (LAMP) assay for Plasmodium recognition in human blood samples, and compared its performance with that of routinely used CSP-ELISAs in field-collected Anopheles coluzzii (Diptera: Culicidae) samples. The nested PCR showed 1.4-fold higher sensitivity than the CSP-ELISA. However, nested PCR results obtained in two laboratories and in different replicates within the same laboratory were not 100% consistent, probably because the copy number of amplifiable Plasmodium mtDNA was close in some specimens to the threshold of nested PCR sensitivity. This implies that Plasmodium-positive specimens should be confirmed by a second nested PCR to avoid false positives. Overall, the results emphasize the need to use molecular approaches to obtain accurate estimates of the actual level of Plasmodium circulation within malaria vector populations.

Unexpectedly high Plasmodium sporozoite rate associated with low human blood index in Anopheles coluzzii from a LLIN-protected village in Burkina Faso

Despite the effectiveness of mass distribution of long-lasting insecticidal nets (LLINs) in reducing malaria transmission in Africa, in hyperendemic areas such as Burkina Faso the burden of malaria remains high. We here report the results of a 4-month survey on the feeding habits and Plasmodium infection in malaria vectors from a village in Burkina Faso one year following a national LLIN distribution programme. Low values of human blood index (HBI) observed in the major malaria vectors in the area (Anopheles coluzzii: N = 263, 20.1%; An. arabiensis: 5.8%, N = 103) are consistent with the hypothesis that LLINs reduced the availability of human hosts to mosquitoes. A regression meta-analysis of data from a systematic review of published studies reporting HBI and sporozoite rates (SR) for An. gambiae complex revealed that the observed SR values (An. coluzzii: 7.6%, N = 503; An. arabiensis: 5.3%, N = 225) are out of the ranges expected based on the low HBI observed. We hypothesize that a small fraction of inhabitants unprotected by bednets acts as a "core group" repeatedly exposed to mosquito bites, representing the major Plasmodium reservoir for the vectors, able to maintain a high risk of transmission even in a village protected by LLINs.

Predicting the likelihood and intensity of mosquito infection from sex specific Plasmodium falciparum gametocyte density

Understanding the importance of gametocyte density on human-to-mosquito transmission is of immediate relevance to malaria control. Previous work (Churcher et al., 2013) indicated a complex relationship between gametocyte density and mosquito infection. Here we use data from 148 feeding experiments on naturally infected gametocyte carriers to show that the relationship is much simpler and depends on both female and male parasite density. The proportion of mosquitoes infected is primarily determined by the density of female gametocytes though transmission from low gametocyte densities may be impeded by a lack of male parasites. Improved precision of gametocyte quantification simplifies the shape of the relationship with infection increasing rapidly before plateauing at higher densities. The mean number of oocysts per mosquito rises quickly with gametocyte density but continues to increase across densities examined. The work highlights the importance of measuring both female and male gametocyte density when estimating the human reservoir of infection.

The Anopheles gambiae 2La chromosome inversion is associated with susceptibility to Plasmodium falciparum in Africa

Chromosome inversions suppress genetic recombination and establish co-adapted gene complexes, or supergenes. The 2La inversion is a widespread polymorphism in the Anopheles gambiae species complex, the major African mosquito vectors of human malaria. Here we show that alleles of the 2La inversion are associated with natural malaria infection levels in wild-captured vectors from West and East Africa. Mosquitoes carrying the more-susceptible allele (2L+^a) are also behaviorally less likely to be found inside houses. Vector control tools that target indoor-resting mosquitoes, such as bednets and insecticides, are currently the cornerstone of malaria control in Africa. Populations with high levels of the 2L+^a allele may form reservoirs of persistent outdoor malaria transmission requiring novel measures for surveillance and control. The 2La inversion is a major and previously unappreciated component of the natural malaria transmission system in Africa, influencing both malaria susceptibility and vector behavior.

Evolution of GOUNDRY, a cryptic subgroup of Anopheles gambiae s.l., and its impact on susceptibility to Plasmodium infection

The recent discovery of a previously unknown genetic subgroup of Anopheles gambiae sensu lato underscores our incomplete understanding of complexities of vector population demographics in Anopheles. This subgroup, named GOUNDRY, does not rest indoors as adults and is highly susceptible to Plasmodium infection in the laboratory. Initial description of GOUNDRY suggested it differed from other known Anopheles taxa in surprising and sometimes contradictory ways, raising a number of questions about its age, population size and relationship to known subgroups. To address these questions, we sequenced the complete genomes of 12 wild-caught GOUNDRY specimens and compared these genomes to a panel of Anopheles genomes. We show that GOUNDRY is most closely related to Anopheles coluzzii, and the timing of cladogenesis is not recent, substantially predating the advent of agriculture. We find a large region of the X chromosome that has swept to fixation in GOUNDRY within the last 100 years, which may be an inversion that serves as a partial barrier to contemporary gene flow. Interestingly, we show that GOUNDRY has a history of inbreeding that is significantly associated with susceptibility to Plasmodium infection in the laboratory. Our results illuminate the genomic evolution of one of probably several cryptic, ecologically specialized subgroups of Anopheles and provide a potent example of how vector population dynamics may complicate efforts to control or eradicate malaria.

An Evolution-Based Screen for Genetic Differentiation between Anopheles Sister Taxa Enriches for Detection of Functional Immune Factors

Nucleotide variation patterns across species are shaped by the processes of natural selection, including exposure to environmental pathogens. We examined patterns of genetic variation in two sister species, Anopheles gambiae and Anopheles coluzzii, both efficient natural vectors of human malaria in West Africa. We used the differentiation signature displayed by a known coordinate selective sweep of immune genes APL1 and TEP1 in A. coluzzii to design a population genetic screen trained on the sweep, classified a panel of 26 potential immune genes for concordance with the signature, and functionally tested their immune phenotypes. The screen results were strongly predictive for genes with protective immune phenotypes: genes meeting the screen criteria were significantly more likely to display a functional phenotype against malaria infection than genes not meeting the criteria (p = 0.0005). Thus, an evolution-based screen can efficiently prioritize candidate genes for labor-intensive downstream functional testing, and safely allow the elimination of genes not meeting the screen criteria. The suite of immune genes with characteristics similar to the APL1-TEP1 selective sweep appears to be more widespread in the A. coluzzii genome than previously recognized. The immune gene differentiation may be a consequence of adaptation of A. coluzzii to new pathogens encountered in its niche expansion during the separation from A. gambiae, although the role, if any of natural selection by Plasmodium is unknown. Application of the screen allowed identification of new functional immune factors, and assignment of new functions to known factors. We describe biochemical binding interactions between immune proteins that underlie functional activity for malaria infection, which highlights the interplay between pathogen specificity and the structure of immune complexes. We also find that most malaria-protective immune factors display phenotypes for either human or rodent malaria, with broad specificity a rarity.

Anopheles gambiae and Anopheles coluzzii are the primary mosquito vectors of human malaria in West Africa. Both of these closely related species efficiently transmit the disease, although they display ecological differences. Previous work showed that A. coluzzii displays distinct genetic patterns in genes important for mosquito immunity. Here, we use this genetic pattern as a filter to examine a panel of potential immune genes, and show that the genetic pattern is strongly predictive for genes that play a functional role in immunity when tested with malaria parasites.

Reticulate Speciation and Barriers to Introgression in the Anopheles gambiae Species Complex

Speciation as a process remains a central focus of evolutionary biology, but our understanding of the genomic architecture and prevalence of speciation in the face of gene flow remains incomplete. The Anopheles gambiae species complex of malaria mosquitoes is a radiation of ecologically diverse taxa. This complex is well-suited for testing for evidence of a speciation continuum and genomic barriers to introgression because its members exhibit partially overlapping geographic distributions as well as varying levels of divergence and reproductive isolation. We sequenced 20 genomes from wild A. gambiae s.s., Anopheles coluzzii, Anopheles arabiensis, and compared these with 12 genomes from the “GOUNDRY” subgroup of A. gambiae s.l. Amidst a backdrop of strong reproductive isolation, we find strong evidence for a speciation continuum with introgression of autosomal chromosomal regions among species and subgroups. The X chromosome, however, is strongly differentiated among all taxa, pointing to a disproportionately large effect of X chromosome genes in driving speciation among anophelines. Strikingly, we find that autosomal introgression has occurred from contemporary hybridization between A. gambiae and A. arabiensis despite strong divergence (∼5× higher than autosomal divergence) and isolation on the X chromosome. In addition to the X, we find strong evidence that lowly recombining autosomal regions, especially pericentromeric regions, serve as barriers to introgression secondarily to the X. We show that speciation with gene flow results in genomic mosaicism of divergence and introgression. Such a reticulate gene pool connecting vector taxa across the speciation continuum has important implications for malaria control efforts.

Association mapping by pooled sequencing identifies TOLL 11 as a protective factor against Plasmodium falciparum in Anopheles gambiae

Background

The genome-wide association study (GWAS) techniques that have been used for genetic mapping in other organisms have not been successfully applied to mosquitoes, which have genetic characteristics of high nucleotide diversity, low linkage disequilibrium, and complex population stratification that render population-based GWAS essentially unfeasible at realistic sample size and marker density.

Methods

We designed a novel mapping strategy for the mosquito system that combines the power of linkage mapping with the resolution afforded by genetic association. We established founder colonies from West Africa, controlled for diversity, linkage disequilibrium and population stratification. Colonies were challenged by feeding on the infectious stage of the human malaria parasite, Plasmodium falciparum, mosquitoes were phenotyped for parasite load, and DNA pools for phenotypically similar mosquitoes were Illumina sequenced. Phenotype-genotype mapping was carried out in two stages, coarse and fine.

Results

In the first mapping stage, pooled sequences were analysed genome-wide for intervals displaying relativereduction in diversity between phenotype pools, and candidate genomic loci were identified for influence upon parasite infection levels. In the second mapping stage, focused genotyping of SNPs from the first mapping stage was carried out in unpooled individual mosquitoes and replicates. The second stage confirmed significant SNPs in a locus encoding two Toll-family proteins. RNAi-mediated gene silencing and infection challenge revealed that TOLL 11 protects mosquitoes against P. falciparum infection.

Conclusions

We present an efficient and cost-effective method for genetic mapping using natural variation segregating in defined recent Anopheles founder colonies, and demonstrate its applicability for mapping in a complex non-model genome. This approach is a practical and preferred alternative to population-based GWAS for first-pass mapping of phenotypes in Anopheles. This design should facilitate mapping of other traits involved in physiology, epidemiology, and behaviour.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2009-z) contains supplementary material, which is available to authorized users.

The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing

Background

Members of the Anophelesgambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anophelescoluzzii by introgression carries point mutations of the voltage-gated sodium channel gene para, including the L1014F kdr mutation associated with insensitivity to pyrethroid insecticides. Carriage of L1014F kdr is also correlated with higher susceptibility to infection with Plasmodium falciparum. However, the genetic mechanism and causative gene(s) underlying the parasite susceptibility phenotype are not known.

Methods

Mosquitoes from the wild Burkina Faso population were challenged by feeding on natural P. falciparum gametocytes. Oocyst infection phenotypes were determined and were tested for association with SNP genotypes. Candidate genes in the detected locus were prioritized and RNAi-mediated gene silencing was used to functionally test for gene effects on P. falciparum susceptibility.

Results

A genetic locus, Pfin6, was identified that influences infection levels of P. falciparum in mosquitoes. The locus segregates as a ~3 Mb haplotype carrying 65 predicted genes including the para gene. The haplotype carrying the kdr allele of para is linked to increased parasite infection prevalence, but many single nucleotide polymorphisms on the haplotype are also equally linked to the infection phenotype. Candidate genes in the haplotype were prioritized and functionally tested. Silencing of para did not influence P. falciparum infection, while silencing of a predicted immune gene, serine protease ClipC9, allowed development of significantly increased parasite numbers.

Conclusions

Genetic variation influencing Plasmodium infection in wild Anopheles is linked to a natural ~3 megabase haplotype on chromosome 2L that carries the kdr allele of the para gene. Evidence suggests that para gene function does not directly influence parasite susceptibility, and the association of kdr with infection may be due to tight linkage of kdr with other gene(s) on the haplotype. Further work will be required to determine if ClipC9 influences the outcome of P. falciparum infection in nature, as well as to confirm the absence of a direct influence by para.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0924-8) contains supplementary material, which is available to authorized users.

Evaluation of a protocol for remote identification of mosquito vector species reveals BG-Sentinel trap as an efficient tool for Anopheles gambiae outdoor collection in Burkina Faso

Background

Feasibility and costs of monitoring efforts aimed to monitor mosquito species are strictly dependent on the presence of skilled entomologists directly in the field. However, in several contexts this is not possible or easy to organize, thus limiting the possibility to obtain crucial information on presence/abundance of potential disease vectors and of new invasive species. Digital imaging approaches could be extremely useful in the frame of medical entomology to overcome this limit. This work describes a surveillance approach to collect and morphologically identify host-seeking malaria vectors based on remote transmission of digital images of specimens collected with ad hoc modified traps.

Methods

CDC light trap (CDC) and the BG-Sentinel trap (BG), both baited with BG-lure and CO₂, were modified in order to have collected mosquitoes immobilized on a bi-dimensional surface. The performance of the two traps in the field was comparatively tested by Latin-square experiments in two villages of Burkina Faso under low and high mosquito densities. The efficiency of identifications based the inspection of digital images versus microscopic identifications of collected specimens was compared.

Results

A total of 1,519 mosquitoes belonging to 16 species were collected, of which 88.5% were microscopically identified as Anopheles gambiae s.l. (mainly Anopheles coluzzii, 85.7%). During dry season BG collected 15 times more females than CDC outdoors, whereas indoors the BG collected 0.4 times less than CDC. During rainy season the ratio BG/CDC was 6.4 and 0.7 outdoors and indoors, respectively. The efficiency of digital images versus microscopic identifications of collected specimens was 97.9%, 95.6% and 81.5% for Culicidae, Anophelinae and An. gambiae s.l., respectively.

Conclusions

Results strongly encourage the use of BG-trap for collecting host-seeking An. gambiae particularly in the outdoor environment, providing new perspectives to the challenge of collecting this fraction of the biting population, whose epidemiological relevance is increasing due to the success of large-scale implementation of indoor malaria vector control strategies. Moreover, results show that the transmission of digital images of specimens collected by the ad hoc modified host-seeking traps efficiently allows identification of malaria vectors, thus opening the perspective to easily carry out mosquito monitoring also in the absence of entomologists directly in the field.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0674-7) contains supplementary material, which is available to authorized users.

Mosquito genomics. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

The Sticky Resting Box, a new tool for studying resting behaviour of Afrotropical malaria vectors

BACKGROUND

Monitoring densities of adult mosquito populations is a major challenge in efforts to evaluate the epidemiology of mosquito-borne diseases, and their response to vector control interventions. In the case of malaria, collection of outdoor-resting Anophelines is rarely incorporated into surveillance and control, partially due to the lack of standardized collection tools. Such an approach, however, is increasingly important to investigate possible changes in mosquito behaviour in response to the scale up of Insecticide Treated Nets and Indoor Residual Spraying. In this study we evaluated the Sticky Resting Box (SRB) - i.e. a sticky variant of previously investigated mosquito Resting Box, which allows passive collection of mosquitoes entering the box - and compared its performance against traditional methods for indoor and outdoor resting mosquito sampling.

METHODS

Daily collections were carried out in two neighbouring villages of Burkina Faso during rainy season 2011 and dry season 2012 by SRB located indoors and outdoors, and by Back-Pack aspiration inside houses (BP) and in ad hoc built outdoor pit-shelters (PIT).

RESULTS

Overall, almost 20,000 Culicidae specimens belonging to 16 species were collected and morphologically identified. Malaria vectors included Anopheles coluzzii (53%), An. arabiensis (12%), An. gambiae s.s. (2.0%) and An. funestus (4.5%). The diversity of species collected in the two villages was similar for SRB and PIT collections outdoors, and significantly higher for SRB than for BP indoors. The population dynamics of An. gambiae s.l. mosquitoes, as obtained by SRB-collections was significantly correlated with those obtained by the traditional methods. The predicted mean estimates of An. gambiae s.l. specimens/sampling-unit/night-of-collections was 6- and 5-times lower for SRB than for BP indoors and PIT outdoors, respectively.

CONCLUSIONS

Overall, the daily performance of SRB in terms of number of malaria vectors/trap was lower than that of traditionally used approaches for in- and outdoor collections. However, unlike these methods, SRB could be set up to collect mosquitoes passively over at least a week. This makes SRB a promising tool for passively monitoring anopheline resting populations, with data presented here providing guidance for how to set up SRB-based collections to acquire information comparable to those obtained with other methods.

Behavioural divergence of sympatric Anopheles funestus populations in Burkina Faso

Background

In Burkina Faso, two chromosomal forms of the malaria vector Anopheles funestus, Folonzo and Kiribina, are distinguished by contrasting frequencies of shared polymorphic chromosomal inversions. Sympatric and synchronous populations of Folonzo and Kiribina mate assortatively, as indicated by a significant deficit of heterokaryotypes, and genetic associations among inversions on independently segregating chromosome arms. The present study aimed to assess, by intensive longitudinal sampling, whether sympatric Folonzo and Kiribina populations are characterized by behavioural differences in key malaria vectorial parameters.

Methods

The study was conducted in two adjacent villages near Ouagadougou, in the dry savanna of central Burkina Faso. Mosquito adult resting behaviour of both forms was compared based on parallel indoor/outdoor collections across six breeding seasons; 8,235 fully karyotyped samples of half-gravid females were analysed in total. Additionally, indoor/outdoor human biting behaviour, host selection, and Plasmodium falciparum sporozoite rate was assessed and compared between chromosomal forms.

Results

The Kiribina form was numerically predominant in the area. However, the Folonzo form was significantly over-represented in indoor resting collections and showed stronger post-prandial endophily, while Kiribina predominated outdoors. Neither form was statistically distinguishable in human biting behaviour, and both were more likely to seek human blood meals indoors than outside. The human blood index and sporozoite rate were comparably high in both chromosomal forms in indoor collections (>89% and >8%, respectively).

Conclusions

Both Kiribina and Folonzo chromosomal forms are formidable malaria vectors in Burkina Faso. However, the significantly greater tendency for the Kiribina form to rest outdoors despite its pronounced anthropophily suggests that uniform exposure of the overall An. funestus population to indoor-based vector control tools cannot be expected; Kiribina is more likely to evade indoor interventions and escape unharmed outdoors, reducing the efficacy of malaria control. Accordingly, more efficient methods to detect Kiribina and Folonzo, and a more complete understanding of their distribution and behaviour in Africa are advocated.

Experimental hut evaluation of Fendona 6SC®-treated bednets and Interceptor® long-lasting nets against Anopheles gambiae s.l. in Burkina Faso

BACKGROUND & OBJECTIVES

Malaria prevention relies heavily on insecticide-treated bednets. Even though the benefits of bednets have been proven that in most of the studies carried out in Africa, their efficacy remains dependent on local conditions. In this study, under field conditions, we evaluated the efficacy of two LLINs (PermaNet® vs Interceptor® ) and two bednet treatment kits (K-O TAB® vs Fendona 6SC® ) against Anopheles gambiae s.l.

METHODS

Bednets were evaluated using experimental huts in the village of Pissy located in the Saponé health district of Burkina Faso. Treatments and sleepers were randomly rotated between huts. Results are expressed in terms of induced exophily, mortality after 24 h and blood-feeding inhibition.

RESULTS

A total of 1392 An. gambiae s.l. mosquitoes were collected during 120 nights in the experimental huts. The overall mortality rates were 85.4% (CL: 79.7-91.4) and 77.5% (CL: 56.9-97.3) for PermaNet® and Interceptor® ,respectively. For the conventionally treated bednets, the mortality was 78.2% (CL: 63.13-96.7) with the Fendona 6SC®-treated nets and 75.5% (CL: 61.2-93) with the K-O TAB®-treated nets. The proportion of blood-fed mosquitoes was significantly higher in the untreated bednet arm than in the treated one, as well as for long-lasting nets than for conventionally treated nets. The entry rate did not vary significantly according to the bednet type, but the treated bednets increased the level of exophily by at least 43%.

CONCLUSION

In the field, the Fendona 6SC® kit and the Interceptor bednets showed comparable efficacy to the already used K-O TAB® kit and PermaNet® bednets. These results could help National Malaria Program managers to formulate appropriate policy for effective vector control.

Three years of insecticide resistance monitoring in Anopheles gambiae in Burkina Faso: resistance on the rise?

Background and methods

A longitudinal Anopheles gambiae s.l. insecticide-resistance monitoring programme was established in four sentinel sites in Burkina Faso. For three years, between 2008 and 2010, WHO diagnostic dose assays were used to measure the prevalence of resistance to all the major classes of insecticides at the beginning and end of the malaria transmission season. Species identification and genotyping for target site mutations was also performed and the sporozoite rate in adults determined.

Results

At the onset of the study, resistance to DDT and pyrethroids was already prevalent in An. gambiae s.l. from the south-west of the country but mosquitoes from the two sites in central Burkina Faso were largely susceptible. Within three years, DDT and permethrin resistance was established in all four sites. Carbamate and organophosphate resistance remains relatively rare and largely confined to the south-western areas although a small number of bendiocarb survivors were found in all sites by the final round of monitoring. The ace-1^R target site resistance allele was present in all localities and its frequency exceeded 20% in 2010 in two of the sites. The frequency of the 1014F kdr mutation increased throughout the three years and by 2010, the frequency of 1014F in all sites combined was 0.02 in Anopheles arabiensis, 0.56 in An. gambiae M form and 0.96 in An. gambiae S form. This frequency did not differ significantly between the sites. The 1014S kdr allele was only found in An. arabiensis but its frequency increased significantly throughout the study (P = 0.0003) and in 2010 the 1014S allele frequency was 0.08 in An. arabiensis. Maximum sporozoite rates (12%) were observed in Soumousso in 2009 and the difference between sites is significant for each year.

Conclusion

Pyrethroid and DDT resistance is now established in An. gambiae s.l. throughout Burkina Faso. Results from diagnostic dose assays are highly variable within and between rounds of testing, and hence it is important that resistance monitoring is carried out on more than one occasion before decisions on insecticide procurement for vector control are made. The presence of 1014S in An. gambiae s.l., in addition to 1014F, is not unexpected given the recent report of 1014S in Benin but highlights the importance of monitoring for both mutations throughout the continent. Future research must now focus on the impact that this resistance is having on malaria control in Burkina Faso.

Development of an allele-specific, loop-mediated, isothermal amplification method (AS-LAMP) to detect the L1014F kdr-w mutation in Anopheles gambiae s. l

Background

Malaria control relies heavily on treated bed nets and indoor residual spraying with pyrethroid insecticides. Unfortunately, the resistance to pyrethroid insecticides, mainly due to the kdr mutation, is spreading in the main malaria vector Anopheles gambiae s.l., decreasing the insecticides’ efficacy. To manage the insecticide resistance rapidly and flexibly, simple and effective tools for the early detection of resistant mosquitoes are needed. This study aimed to develop an allele-specific, loop-mediated, isothermal amplification (AS-LAMP) method to detect the West African-type kdr mutation (kdr-w; L1014F) in field-collected mosquitoes.

Methods

DNA fragments of the wild-type and the mutated kdr gene were used to select the primers and develop the method. The primers were designed with the mutation at the 5’ end of the backward inner primer (BIP). The AS-LAMP method was compared to the AS-PCR method using the genomic DNA of 120 field-collected mosquitoes.

Results

The AS-LAMP method could discriminate between the wild-type homozygote, the heterozygote, and the kdr-w homozygote within 75 min. The AS-LAMP method has the advantage of being faster and at least as sensitive and specific as the AS-PCR method.

Conclusions

The AS-LAMP method can be used to detect the kdr mutation for quick decision-making, even in less well-equipped laboratories.

A cryptic subgroup of Anopheles gambiae is highly susceptible to human malaria parasites

Population subgroups of the African malaria vector Anopheles gambiae have not been comprehensively characterized owing to the lack of unbiased sampling methods. In the arid savanna zone of West Africa, where potential oviposition sites are scarce, widespread collection from larval pools in the peridomestic human habitat yielded a comprehensive genetic survey of local A. gambiae population subgroups, independent of adult resting behavior and ecological preference. A previously unknown subgroup of exophilic A. gambiae is sympatric with the known endophilic A. gambiae in this region. The exophilic subgroup is abundant, lacks differentiation into M and S molecular forms, and is highly susceptible to infection with wild Plasmodium falciparum. These findings might have implications for the epidemiology of malaria transmission and control.

De novo transcriptome sequencing in Anopheles funestus using Illumina RNA-seq technology

BACKGROUND

Anopheles funestus is one of the primary vectors of human malaria, which causes a million deaths each year in sub-Saharan Africa. Few scientific resources are available to facilitate studies of this mosquito species and relatively little is known about its basic biology and evolution, making development and implementation of novel disease control efforts more difficult. The An. funestus genome has not been sequenced, so in order to facilitate genome-scale experimental biology, we have sequenced the adult female transcriptome of An. funestus from a newly founded colony in Burkina Faso, West Africa, using the Illumina GAIIx next generation sequencing platform.

METHODOLOGY/PRINCIPAL FINDINGS

We assembled short Illumina reads de novo using a novel approach involving iterative de novo assemblies and "target-based" contig clustering. We then selected a conservative set of 15,527 contigs through comparisons to four Dipteran transcriptomes as well as multiple functional and conserved protein domain databases. Comparison to the Anopheles gambiae immune system identified 339 contigs as putative immune genes, thus identifying a large portion of the immune system that can form the basis for subsequent studies of this important malaria vector. We identified 5,434 1:1 orthologues between An. funestus and An. gambiae and found that among these 1:1 orthologues, the protein sequence of those with putative immune function were significantly more diverged than the transcriptome as a whole. Short read alignments to the contig set revealed almost 367,000 genetic polymorphisms segregating in the An. funestus colony and demonstrated the utility of the assembled transcriptome for use in RNA-seq based measurements of gene expression.

CONCLUSIONS/SIGNIFICANCE

We developed a pipeline that makes de novo transcriptome sequencing possible in virtually any organism at a very reasonable cost ($6,300 in sequencing costs in our case). We anticipate that our approach could be used to develop genomic resources in a diversity of systems for which full genome sequence is currently unavailable. Our An. funestus contig set and analytical results provide a valuable resource for future studies in this non-model, but epidemiologically critical, vector insect.

Seasonal distribution of Anopheles funestus chromosomal forms from Burkina Faso

Background

Previous studies of Anopheles funestus chromosomal inversion polymorphisms in Burkina Faso showed large departures from Hardy-Weinberg equilibrium and linkage disequilibrium among inversions located on different chromosomes, implying the existence of two taxonomic units ("chromosomal forms") with limited genetic flow. One chromosomal form, named Folonzo, is highly polymorphic for alternative rearrangements of 3Ra, 3Rb, 2Ra, and 3La; the other, Kiribina, is predominantly characterized by the standard arrangement of these inversions. To investigate the temporal distribution of these chromosomal forms, further collections were carried out in two villages near Ouagadougou where they are found in sympatry.

Methods

Chromosomal karyotypes were determined from indoor-resting, half-gravid females sampled within and across six breeding seasons, from December 1998 to April 2007.

Results

As expected, the pattern of chromosomal polymorphism in An. funestus was consistent with assortatively mating Folonzo and Kiribina forms. When samples were assigned to each chromosomal form, their relative abundance varied within successive breeding seasons in a repeating pattern of temporal variability. Relative abundance of the Folonzo form was correlated with climatic variables related to temperature and rainfall.

Conclusion

The relative abundance of Folonzo and Kiribina forms of An. funestus likely reflects different larval ecologies that are linked to varying climatic conditions. Further analysis of the bionomics of these vectors is recommended in light of its relevance to vector control.

Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae

BACKGROUND

Ongoing lineage splitting within the African malaria mosquito Anopheles gambiae is compatible with ecological speciation, the evolution of reproductive isolation by divergent natural selection acting on two populations exploiting alternative resources. Divergence between two molecular forms (M and S) identified by fixed differences in rDNA, and characterized by marked, although incomplete, reproductive isolation is occurring in West and Central Africa. To elucidate the role that ecology and geography play in speciation, we carried out a countrywide analysis of An. gambiae M and S habitat requirements, and that of their chromosomal variants, across Burkina Faso.

RESULTS

Maps of relative abundance by geostatistical interpolators produced a distinct pattern of distribution: the M-form dominated in the northernmost arid zones, the S-form in the more humid southern regions. Maps of habitat suitability, quantified by Ecological Niche Factor Analysis based on 15 eco-geographical variables revealed less contrast among forms. M was peculiar as it occurred proportionally more in habitat of marginal quality. Measures of ecological niche breadth and overlap confirmed the mismatch between the fundamental and realized patterns of habitat occupation: forms segregated more than expected from the extent of divergence of their environmental envelope--a signature of niche expansion. Classification of chromosomal arm 2R karyotypes by multilocus genetic clustering identified two clusters loosely corresponding to molecular forms, with 'mismatches' representing admixed individuals due to shared ancestral polymorphism and/or residual hybridization. In multivariate ordination space, these karyotypes plotted in habitat of more marginal quality compared to non-admixed, 'typical', karyotypes. The distribution of 'typical' karyotypes along the main eco-climatic gradient followed a consistent pattern within and between forms, indicating an adaptive role of inversions at this geographical scale.

CONCLUSION

Ecological segregation between M and S is consistent with niche expansion into marginal habitats by chromosomal inversion variants during early lineage divergence; presumably, this process is promoted by inter-karyotype competition in the higher-quality core habitat. We propose that the appearance of favourable allelic combinations in other regions of suppressed recombination (e.g. pericentromeric portions defining speciation islands in An. gambiae) fosters development of reproductive isolation to protect linkage between separate chromosomal regions.

Divergence with gene flow in Anopheles funestus from the Sudan Savanna of Burkina Faso, West Africa

Anopheles funestus is a major vector of malaria across Africa. Understanding its complex and nonequilibrium population genetic structure is an important challenge that must be overcome before vector populations can be successfully perturbed for malaria control. Here we examine the role of chromosomal inversions in structuring genetic variation and facilitating divergence in Burkina Faso, West Africa, where two incipient species (chromosomal forms) of A. funestus, defined principally by rearrangements of chromosome 3R, have been hypothesized. Sampling across an approximately 300-km east-west transect largely contained within the Sudan-Savanna ecoclimatic zone, we analyzed chromosomal inversions, 16 microsatellite loci distributed genomewide, and 834 bp of the mtDNA ND5 gene. Both molecular markers revealed high genetic diversity, nearly all of which was accounted for by within-population differences among individuals, owing to recent population expansion. Across the study area there was no correlation between genetic and geographic distance. Significant genetic differentiation found between chromosomal forms on the basis of microsatellites was not genomewide but could be explained by chromosome 3R alone on the basis of loci inside and near inversions. These data are not compatible with complete reproductive isolation but are consistent with differential introgression and sympatric divergence between the chromosomal forms, facilitated by chromosome 3R inversions.

Chromosomal evidence of incipient speciation in the Afrotropical malaria mosquito Anopheles funestus

The analysis of chromosomal polymorphism of paracentric inversions in anopheline mosquitoes has often been instrumental to the discovery of sibling species complexes and intraspecific genetic heterogeneities associated with incipient speciation processes. To investigate the population structure of Anopheles funestus Giles (Diptera: Culicidae), one of the three most important vectors of human malaria in sub-Saharan Africa, a three-year survey of chromosomal polymorphism was carried out on 4,638 karyotyped females collected indoors and outdoors from two villages of central Burkina Faso. Large and temporally stable departures from Hardy-Weinberg equilibrium due to significant deficits of heterokaryotypes were found irrespective of the place of capture, and of the spatial and temporal units chosen for the analysis. Significant linkage disequilibrium was observed among inversion systems on independently assorting chromosomal arms, indicating the existence of assortative mating phenomena. Results were consistent with the existence of two chromosomal forms characterized by contrasting degrees of inversion polymorphism maintained by limitations to gene flow. This hypothesis was supported by the reestablishment of Hardy-Weinberg and linkage equilibria when individual specimens were assigned to each chromosomal form according to two different algorithms. This pattern of chromosomal variability is suggestive of an incipient speciation process in An. funestus populations from Burkina Faso.

Molecular differentiation between chromosomally defined incipient species of Anopheles funestus

Anopheles funestus Giles is one of the most important vectors of malaria in sub-Saharan Africa. The population structure of this mosquito in Burkina Faso, West Africa based on chromosomal inversion data led to the description of two chromosomal forms, Kiribina and Folonzo. Because both forms co-occur in the same locales yet differ significantly, both in the frequency of inverted arrangements on chromosome arms 3R and 2R and in vectorial capacity, they were hypothesized to be emerging species with at least partial barriers to gene flow. This hypothesis would be strengthened by molecular evidence of differentiation between Kiribina and Folonzo at loci outside chromosomal inversions. We surveyed molecular variation in sympatric populations of the two forms using sequences from the mitochondrial ND5 gene and genotypes at sixteen microsatellite loci distributed across the genome. Both classes of marker revealed slight but significant differentiation between the two forms (mtDNA F(ST) = 0.023, P < 0.001; microsatellite F(ST) = 0.004, P < 0.001; R(st) = 0.009, P = 0.002). Locus-by-locus analysis of the microsatellite data showed that significant differentiation was not genome-wide, but could be attributed to five loci on chromosome 3R (F(ST) = 0.010, P < 0.001; R(st) = 0.016, P = 0.002). Importantly, three of these loci are outside of, and in linkage equilibrium with, chromosomal inversions, suggesting that differentiation between chromosomal forms extends beyond the inversions themselves. The slight overall degree of differentiation indicated by both marker classes is likely an underestimate because of recent population expansion inferred for both Folonzo and Kiribina. The molecular evidence from this study is consistent with the hypothesis of incipient speciation between Kiribina and Folonzo.

A microsatellite map of the African human malaria vector Anopheles funestus

Microsatellite markers and chromosomal inversion polymorphisms are useful genetic markers for determining population structure in Anopheline mosquitoes. In Anopheles funestus (2N = 6), only chromosome arms 2R, 3R, and 3L are known to carry polymorphic inversions. The physical location of microsatellite markers with respect to polymorphic inversions is potentially important information for interpreting population genetic structure, yet none of the available marker sets have been physically mapped in this species. Accordingly, we mapped 32 polymorphic A. funestus microsatellite markers to the polytene chromosomes using fluorescent in situ hybridization (FISH) and identified 16 markers outside of known polymorphic inversions. Here we provide an integrated polytene chromosome map for A. funestus that includes the breakpoints of all known polymorphic inversions as well as the physical locations of microsatellite loci developed to date. Based on this map, we suggest a standard set of 16 polymorphic microsatellite markers that are distributed evenly across the chromosome complement, occur predominantly outside of inversions, and amplify reliably. Adoption of this set by researchers working in different regions of Africa will facilitate metapopulation analyses of this primary malaria vector.

Ovicidal and larvicidal activity against Aedes aegypti and Anopheles gambiae complex mosquitoes of essential oils extracted from three spontaneous plants of Burkina Faso

Essential oils extracted from dried leaves of three spontaneous plants naturally growing in Burkina Faso, i.e. Cymbopogon proximus, Lippia multiflora and Ocimum canum, exhibited larvicidal activity by the WHO standard protocol against 3rd and 4th instar F1-larvae of field-collected mosquitoes vectors of human disease, namely Aedes aegypti and members of the Anopheles gambiae complex, An. arabiensis and An. gambiae. The median lethal concentration (LC50) for Ae. aegypti and An. gambiae s.l. larvae ranged between 53.5-258.5 ppm and 61.9-301.6 ppm, respectively. The LC90 estimates ranged 74.8-334.8 ppm for Ae. aegypti, and 121.6-582.9 ppm for An. gambiae s.l. Ovicidal activity against eggs of An. gambiae s.l. was also demonstrated. The LC50 values for An. gambiae s.l. eggs ranged between 17.1-188.7 ppm, while LC90 values ranged between 33.5-488 ppm. Lippia multiflora showed the highest activity against An. gambiae s.l. eggs and Ae. aegypti larvae, whereas no difference was found among C. proximus and L. multiflora in their activity against An. gambiae s.l. larvae. Of the three plants, essential oils from O. canum had the lowest activity against both eggs and larvae. Eggs were more susceptible than larvae. Ae. aegypti larvae were more susceptible than larvae of An. gambiae s.l.