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

Distinct physiological, biochemical and morphometric adjustments in the malaria vectors Anopheles gambiae and A. coluzzii as means to survive dry season conditions in Burkina Faso

Aestivation and dispersive migration are the two strategies evoked in the literature to explain the way in which malaria vectors Anopheles coluzzii and A. gambiae survive the harsh climatic conditions of the dry season in sub-Saharan Africa. However, the physiological mechanisms regulating these two strategies are unknown. In the present study, mosquito species were exposed to controlled environmental conditions mimicking the rainy and dry seasons of south western Burkina Faso. Survival strategies were studied through morphometric (wing length), ecophysiological (respiratory gas exchanges), biochemical (cuticular hydrocarbons composition) and molecular (AKH mRNA expression levels) parameters, variations of which are usually considered to be hallmarks of aestivation and dispersion mechanisms in various insects. Our results showed that ecophysiological and morphometric adjustments are made in both species to prevent water losses during the dry season. However, the usual metabolic rate modifications expected as signatures of aestivation and migration were not observed, highlighting specific and original physiological mechanisms sustaining survival in malaria mosquitoes during the dry season. Differences in epicuticular hydrocarbon composition and AKH levels of expression were found between the permanent and temporary A. coluzzii populations, illustrating the great phenotypic plasticity of this mosquito species. Altogether, our work underlines the diverse and complex pattern of changes occurring in the two mosquito species and at the population level to cope with the dry season and highlights potential targets of future control tools.

Collision sensitive niche profile of the worst affected bird-groups at wind turbine structures in the Federal State of Brandenburg, Germany

Biodiversity-related impacts at wind energy facilities have increasingly become a cause of conservation concern, central issue being the collision of birds. Utilizing spatial information of their carcass detections at wind turbines (WTs), we quantified the detections in relation to the metric distances of the respective turbines to different land-use types. We used ecological niche factor analysis (ENFA) to identify combinations of land-use distances with respect to the spatial allocation of WTs that led to higher proportions of collisions among the worst affected bird-groups: Buntings, Crows, Larks, Pigeons and Raptors. We also assessed their respective similarities to the collision phenomenon by checking for overlaps amongst their distance combinations. Crows and Larks showed the narrowest “collision sensitive niche”; a part of ecological niche under higher risk of collisions with turbines, followed by that of Buntings and Pigeons. Raptors had the broadest niche showing significant overlaps with the collision sensitive niches of the other groups. This can probably be attributed to their larger home range combined with their hunting affinities to open landscapes. Identification of collision sensitive niches could be a powerful tool for landscape planning; helping avoid regions with higher risks of collisions for turbine allocations and thus protecting sensitive bird populations.

Dissecting functional components of reproductive isolation among closely related sympatric species of the Anopheles gambiae complex

Explaining how and why reproductive isolation evolves and determining which forms of reproductive isolation have the largest impact on the process of population divergence are major goals in the study of speciation. By studying recent adaptive radiations in incompletely isolated taxa, it is possible to identify barriers involved at early divergence before other confounding barriers emerge after speciation is complete. Sibling species of the Anopheles gambiae complex offer opportunities to provide insights into speciation mechanisms. Here, we studied patterns of reproductive isolation among three taxa, Anopheles coluzzii, An. gambiae s.s. and Anopheles arabiensis, to compare its strength at different spatial scales, to dissect the relative contribution of pre‐ versus postmating isolation, and to infer the involvement of ecological divergence on hybridization. Because F1 hybrids are viable, fertile and not uncommon, understanding the dynamics of hybridization in this trio of major malaria vectors has important implications for how adaptations arise and spread across the group, and in planning studies of the safety and efficacy of gene drive as a means of malaria control. We first performed a systematic review and meta‐analysis of published surveys reporting on hybrid prevalence, showing strong reproductive isolation at a continental scale despite geographically restricted exceptions. Second, we exploited our own extensive field data sets collected at a regional scale in two contrasting environmental settings, to assess: (i) levels of premating isolation; (ii) spatio/temporal and frequency‐dependent dynamics of hybridization, (iii) relationship between reproductive isolation and ecological divergence and (iv) hybrid viability penalty. Results are in accordance with ecological speciation theory predicting a positive association between the strength of reproductive isolation and degree of ecological divergence, and indicate that postmating isolation does contribute to reproductive isolation among these species. Specifically, only postmating isolation was positively associated with ecological divergence, whereas premating isolation was correlated with phylogenetic distance.

Pre-copula acoustic behaviour of males in the malarial mosquitoes Anopheles coluzzii and Anopheles gambiae s.s. does not contribute to reproductive isolation

We reveal that males of two members of the Anopheles gambiae s.l. species complex, Anopheles coluzzii and Anopheles gambiae s.s. (hereafter A. gambiae), which are both malaria vectors, perform a stereotypical acoustic behaviour in response to pure tones at frequencies that encompass the frequency range of the female's flight-tones. This behaviour resembles that described for Culex quinquefasciatus and consists of phonotactic flight initiated by a steep increase in wing-beat frequency (WBF) followed by rapid frequency modulation (RFM) of WBF when in close proximity to the sound source. RFM was elicited without acoustic feedback or the presence of a live female, but it appears to be a stereotypic behaviour in the immediate lead up to copula formation. RFM is an independent and different behavioural process from harmonic convergence interactions used by male-female pairs for mate recognition at earlier stages of mating. Acoustic threshold for RFM was used to plot behavioural audiograms from free-flying A coluzzii and A gambiae males. These audiograms were almost identical (minima ∼400 Hz) and encompassed the WBF ranges of A coluzzii (378-601 Hz) and A gambiae (373-590 Hz) females, indicating that males of the two species share similar frequency tuning and range. Furthermore, no differences were found between the two species in their WBFs, RFM behaviour or harmonic convergence ratios. These results indicate that assortative mating between A coluzzii and A gambiae is unlikely to be based on male-specific acoustic behaviours during RFM. The significance of these findings in relation to possible mechanisms for assortative mating is discussed.

Chromosome inversions and ecological plasticity in the main African malaria mosquitoes

Chromosome inversions have fascinated the scientific community, mainly because of their role in the rapid adaption of different taxa to changing environments. However, the ecological traits linked to chromosome inversions have been poorly studied. Here, we investigated the roles played by 23 chromosome inversions in the adaptation of the four major African malaria mosquitoes to local environments in Africa. We studied their distribution patterns by using spatially explicit modeling and characterized the ecogeographical determinants of each inversion range. We then performed hierarchical clustering and constrained ordination analyses to assess the spatial and ecological similarities among inversions. Our results show that most inversions are environmentally structured, suggesting that they are actively involved in processes of local adaptation. Some inversions exhibited similar geographical patterns and ecological requirements among the four mosquito species, providing evidence for parallel evolution. Conversely, common inversion polymorphisms between sibling species displayed divergent ecological patterns, suggesting that they might have a different adaptive role in each species. These results are in agreement with the finding that chromosomal inversions play a role in Anopheles ecotypic adaptation. This study establishes a strong ecological basis for future genome-based analyses to elucidate the genetic mechanisms of local adaptation in these four mosquitoes.

Characterisation of Species and Diversity of Anopheles gambiae Keele Colony

Anopheles gambiae sensu stricto was recently reclassified as two species, An. coluzzii and An. gambiae s.s., in wild-caught mosquitoes, on the basis of the molecular form, denoted M or S, of a marker on the X chromosome. The An. gambiae Keele line is an outbred laboratory colony strain that was developed around 12 years ago by crosses between mosquitoes from 4 existing An. gambiae colonies. Laboratory colonies of mosquitoes often have limited genetic diversity because of small starting populations (founder effect) and subsequent fluctuations in colony size. Here we describe the characterisation of the chromosomal form(s) present in the Keele line, and investigate the diversity present in the colony using microsatellite markers on chromosome 3. We also characterise the large 2La inversion on chromosome 2. The results indicate that only the M-form of the chromosome X marker is present in the Keele colony, which was unexpected given that 3 of the 4 parent colonies were probably S-form. Levels of diversity were relatively high, as indicated by a mean number of microsatellite alleles of 6.25 across 4 microsatellites, in at least 25 mosquitoes. Both karyotypes of the inversion on chromosome 2 (2La/2L+^a) were found to be present at approximately equal proportions. The Keele colony has a mixed M- and S-form origin, and in common with the PEST strain, we propose continuing to denote it as an An. gambiae s.s. line.

Comparative physiological plasticity to desiccation in distinct populations of the malarial mosquito Anopheles coluzzii

Background

In West Africa, populations of the malaria vector mosquito, Anopheles coluzzii, are seasonally exposed to strong desiccating conditions during the dry season. Their dynamics strictly follows the pace of the availability of suitable larval development sites (water collections). Accordingly, mosquitoes can reproduce all year long where permanent breeding is possible, or stop reproduction and virtually disappear at the onset of the dry season when surface water dries up, like observed in temporary habitats of dry savannah areas. This highlights the strong adaptive abilities of this mosquito species, which relies at least in part, upon physiological and molecular mechanisms of specific signatures.

Methods

Here, we analysed a range of physiological and molecular responses expressed by geographically different populations of An. coluzzii inhabiting permanent and temporary breeding sites from the north and the south-west of Burkina Faso. Four mosquito colonies, namely (i) Oursi, built from females breeding in permanent habitats of the north; (ii) Déou, from temporary northern habitats; (iii) Soumousso from south-western temporary breeding sites; and (iv) Bama, from permanent habitats of the same south-western zone, were reared in climatic chambers under contrasted environmental conditions, mimicking temperature, relative humidity and light regimen occurring in northern Burkina Faso. Female mosquitoes were analysed for the seasonal variation in their amounts of proteins, triglycerides and free-circulating metabolites. The expression level of genes coding for the adipokinetic (AKH-I) and the AKH/corazonin-related peptides (ACP) were also assessed and compared among populations and environmental conditions.

Results

Our analysis did not reveal an apparent pattern of physiological and molecular variations strictly correlated with either the larval ecotype or the geographical origin of the mosquitoes. However, specific distinct responses were observed among populations, suggesting that dry season survival may rely on more complex ecological parameters at a micro-habitat scale. Interestingly, the physiological and molecular data support the hypothesis that different aestivation abilities exist among populations of An. coluzzii inhabiting contrasted ecological settings. In particular, the striking metabotypes differentiation and the AKH mRNA expression level observed in females from temporary northern populations may suggest the existence of a “strong” aestivation strategy in these specimens.

Conclusion

Our work provides insights into the physiological and molecular basis of dry and rainy season responses in An. coluzzii, and highlights the important diversity of the mechanisms involved. Such results represent key data for understanding the ecophysiological mechanisms underpinning the strong adaptive potential of this malaria vector species, which undoubtedly contributes to the spreading of mosquito distribution areas in space and time.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1854-1) contains supplementary material, which is available to authorized users.

[New insight into Anopheles coluzzii Coetzee & Wilkerson, 2013 larval ecology in the Wouri estuary, Littoral-Cameroon]

Malaria vectors control is essentially based on the use of insecticides against adult mosquitoes. However because of the development of resistance to insecticides, there is now a renewed interest in the management of larval sources. The aim of the present study was to map and characterize the breeding sites of Anopheles coluzzii in the Wouri river estuary in Cameroon. Larval surveys were carried out between December 2013 and August 2014 in rural areas on the island of Manoka and urban area in Youpwe at Douala. Culicidae breeding sites identified were georeferenced and mapped. Their larval productivity was evaluated by the method of "dipping" and their physicochemical parameters measured by spectrophotometry and oximetry. Culicidae collected larvae were reared in the insectarium to the adult stage. Adult mosquitoes were subjected to morphological identification and those belonging to the Anopheles gambiae complex have subsequently been subjected to molecular identification by the PCR-RFLP technique. A total of 240 breeding sites were geo-referenced in the two sites, including 10 types. Abandoned containers and pools were the most frequent breeding sites respectively in Manoka and in Youpwe. After morphological and molecular identification, eleven mosquito species have been identified. Anopheles coluzzii and Culex quinquefasciatus were the most frequent species respectively in Manoka and in Youpwe. Mosquito density was higher in managed gutters and canoes respectively in Manoka and in Youpwe. Culex and Aedes genus were more frequent in the hollow palm and water wells respectively in Manoka and Youpwe. The productivity of breeding sites varied according to the physicochemical parameters. Species richness varied according to the type of breeding site. Anopheles coluzzii was observed for the first time in Cameroon in water storage containers, tires, discarded containers and canoes. This study highlighted diversity in the type of breeding site of An. coluzzii in the Wouri estuary, suggesting the adaptation of this species in its environment. These results could be used to develop an antilarval control strategy in Manoka and in Youpwe.

Interactive cost of Plasmodium infection and insecticide resistance in the malaria vector Anopheles gambiae

Insecticide resistance raises concerns for the control of vector-borne diseases. However, its impact on parasite transmission could be diverse when considering the ecological interactions between vector and parasite. Thus we investigated the fitness cost associated with insecticide resistance and Plasmodium falciparum infection as well as their interactive cost on Anopheles gambiae survival and fecundity. In absence of infection, we observed a cost on fecundity associated with insecticide resistance. However, survival was higher for mosquito bearing the kdr mutation and equal for those with the ace-1(R) mutation compared to their insecticide susceptible counterparts. Interestingly, Plasmodium infection reduced survival only in the insecticide resistant strains but not in the susceptible one and infection was associated with an increase in fecundity independently of the strain considered. This study provides evidence for a survival cost associated with infection by Plasmodium parasite only in mosquito selected for insecticide resistance. This suggests that the selection of insecticide resistance mutation may have disturbed the interaction between parasites and vectors, resulting in increased cost of infection. Considering the fitness cost as well as other ecological aspects of this natural mosquito-parasite combination is important to predict the epidemiological impact of insecticide resistance.

Spatial and Temporal Trends in Insecticide Resistance among Malaria Vectors in Chad Highlight the Importance of Continual Monitoring

Background

A longitudinal Anopheles gambiae s.l. insecticide resistance monitoring programme was established in four sentinel sites in Chad 2008–2010. When this programme ended, only sporadic bioassays were performed in a small number of sites.

Methods

WHO diagnostic dose assays were used to measure the prevalence of insecticide resistance to 0.1% bendiocarb, 4% DDT, 0.05% deltamethrin, 1% fenitrothion, and 0.75% permethrin in the main malaria vectors at the beginning and end of the malaria transmission season for three years 2008–2010, with subsequent collections in 2011 and 2014. Species and molecular identification of An. gambiae M and S forms and kdr genotyping was performed using PCR-RLFP; circumsporozoite status was assessed using ELISA.

Results

Between 2008 and 2010, significant changes in insecticide resistance profiles to deltamethrin and permethrin were seen in 2 of the sites. No significant changes were seen in resistance to DDT in any site during the study period. Testing performed after the period of routine monitoring had ended showed dramatic increases to DDT and pyrethroid resistance in 3 sites. No resistance to organophosphate or carbamate insecticides was detected. An. arabiensis was the predominate member of the An. gambiae complex in all 4 sites; adult collections showed temporal variation in species composition in only 1 site. Kdr analysis identified both 1014F and 1014S alleles in An. gambiae S only. Circumsporozoite analysis showed the highest vector infection rates were present in Donia, a site with extensive use of agricultural insecticides.

Conclusions

During the monitoring gap of four years, significant changes occurred in resistance prevalence in 3 of the 4 sites (p = <0.001), endangering the efficacy of currently implemented malaria control interventions. Significant changes in insecticide resistance profiles and a lack of kdr resistance alleles in adult populations highlight the urgent need for comprehensive entomological monitoring to be implemented and sustained in country.

Anopheles coluzzii larval habitat and insecticide resistance in the island area of Manoka, Cameroon

BACKGROUND

The effectiveness of Long-Lasting Insecticidal Nets and Indoor Residual Spraying in malaria vector control is threatened by vector resistance to insecticides. Knowledge of mosquito habitats and patterns of insecticide resistance would facilitate the development of appropriate vector control strategies. Therefore, we investigated An. coluzzii larval habitats and resistance to insecticides in the Manoka rural island area compared with the Youpwe suburban inland area, in Douala VI and II districts respectively.

METHODS

Anopheline larvae and pupae were collected from open water bodies in December 2013 and April 2014 and reared until adult emergence. Two to four day old emerging females were morphologically identified as belonging to the An. gambiae complex and used for WHO susceptibility tests with 4 % DDT, 0.75 % permethrin, and 0.05 % deltamethrin, with or without piperonyl butoxide (PBO) synergist. Control and surviving specimens were identified down to the species using a PCR-RFLP method. Survivors were genotyped for kdr L1014 mutations using Hot Oligonucleotide Ligation Assay.

RESULTS

In both study sites, ponds, residual puddles, boats, and drains were identified as the major An. gambiae s.l. larval habitats. A total of 1397 females, including 784 specimens from Manoka and 613 from Youpwe, were used for resistance testing. The two mosquito populations displayed resistance to DDT, permethrin and deltamethrin, with variable mortality rates from 1 % to 90 %. The knock-down times were also significantly increased (at least 2.8 fold). Pre-exposure of mosquitoes to PBO did not impact on their mortality to DDT, conversely the mortality rates to permethrin and deltamethrin were significantly increased (7.56 ≤ X(2) ≤ 48.63, df = 1, p < 0.01), suggesting involvement of P450 oxidases in pyrethroid resistance. A subsample of 400 An. gambiae s.l. specimens including 280 control and 120 survivors from bioassays were all found to be An. coluzzii species. Only the kdr 1014 F mutation was found in survivors, with 88.5 % (N = 76) and 75 % (N = 44) frequencies in Youpwe and Manoka respectively.

CONCLUSION

This is the first report of An. coluzzii resistance to insecticides in an insular area in Cameroon. Since permanent larval habitats have been identified, larval source management strategies may be trialed in this area as complementary vector control interventions.

Insecticide susceptibility of natural populations of Anopheles coluzzii and Anopheles gambiae (sensu stricto) from Okyereko irrigation site, Ghana, West Africa

BACKGROUND

The increasing spread of insecticide resistance in malaria vectors has been well documented across sub-Saharan Africa countries. The influence of irrigation on increasing vector resistance is poorly understood, and is critical to successful and ethical implementation of food security policies. This study investigated the insecticide resistance status of An. gambiae (s.l.) mosquitoes collected from the irrigated rice area of Okyereko, a village containing about 42 hectares of irrigated field within an irrigation project plan in the Central Region of Ghana. Large amounts of insecticides, herbicides and fertilizers are commonly used in the area to boost the annual production of the rice.

METHODS

Mosquito larvae were collected and adults were assayed from the F1 progeny. The resistance status, allele and genotype were characterized using WHO susceptibility testing and PCR methods respectively.

RESULTS

The An. gambiae (s.l.) populations from Okyereko are highly resistant to DDT and pyrethroid insecticides, with possible involvement of metabolic mechanisms including the elevation of P450 and GST enzyme as well as P-gp activity. The population was mostly composed of An. coluzzii specimens (more than 96 %) with kdr and ace-1 frequencies of 0.9 and 0.2 %, respectively.

CONCLUSION

This study brings additional information on insecticide resistance and the characterization of An. gambiae (s.l.) mosquitoes from Okyereko, which can be helpful in decision making for vector control programmes in the region.

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.

Does insecticide resistance contribute to heterogeneities in malaria transmission in The Gambia?

BACKGROUND

Malaria hotspots, areas with consistently higher than average transmission, may become increasingly common as malaria declines. This phenomenon, currently observed in The Gambia, may be caused by several factors, including some related to the local vectors, whose contribution is poorly understood.

METHODS

Using WHO susceptibility bioassays, insecticide resistance status was determined in vector populations sampled from six pairs of villages across The Gambia, each pair contained a low and high prevalence village.

RESULTS

Three vector species were observed (23.5% Anopheles arabiensis, 31.2% Anopheles gambiae, 43.3% Anopheles coluzzii and 2.0% An. coluzzii × An. gambiae hybrids). Even at a fine scale, significant differences in species composition were detected within village pairs. Resistance to both DDT and deltamethrin was more common in An. gambiae, most markedly in the eastern part of The Gambia and partly attributable to differing frequencies of resistance mutations. The Vgsc-1014F target site mutation was strongly associated with both DDT (OR = 256.7, (95% CI 48.6-6374.3, p < 0.001) and deltamethrin survival (OR = 9.14, (95% CI 4.24-21.4, p < 0.001). A second target site mutation, Vgsc-1575Y, which co-occurs with Vgsc-1014F, and a metabolic marker of resistance, Gste2-114T, conferred additional survival benefits to both insecticides. DDT resistance occurred significantly more frequently in villages with high malaria prevalence (p = 0.025) though this did not apply to deltamethrin resistance.

CONCLUSION

Whilst causality of relationships requires further investigation, variation in vector species and insecticide resistance in The Gambia is associated with malaria endemicity; with a notably higher prevalence of infection and insecticide resistance in the east of the country. In areas with heterogeneous malaria transmission, the role of the vector should be investigated to guide malaria control interventions.

Genetic Structure of a Local Population of the Anopheles gambiae Complex in Burkina Faso

Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. Population heterogeneities for ecological and behavioral attributes expand and stabilize malaria transmission over space and time, and populations may change in response to vector control, urbanization and other factors. There is a need for approaches to comprehensively describe the structure and characteristics of a sympatric local mosquito population, because incomplete knowledge of vector population composition may hinder control efforts. To this end, we used a genome-wide custom SNP typing array to analyze a population collection from a single geographic region in West Africa. The combination of sample depth (n = 456) and marker density (n = 1536) unambiguously resolved population subgroups, which were also compared for their relative susceptibility to natural genotypes of Plasmodium falciparum malaria. The population subgroups display fluctuating patterns of differentiation or sharing across the genome. Analysis of linkage disequilibrium identified 19 new candidate genes for association with underlying population divergence between sister taxa, A. coluzzii (M-form) and A. gambiae (S-form).

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.

Larval ecology of Anopheles coluzzii in Cape Coast, Ghana: water quality, nature of habitat and implication for larval control

Background

There is a growing interest in larval control intervention to supplement existing malaria control strategies, particularly in urban areas. However, effective implementation requires a good understanding of habitat ecology of Anopheles mosquitoes. Clean water bodies have long been reported by several studies as a preferred breeding habitat for Anopheles gambiae. Other studies have also reported the breeding of An. gambiae in polluted water bodies. However, the term clean or polluted is mostly based on visual examination and is not well defined. This study was conducted with the aim of assessing water quality in Anopheles breeding habitats and the practicability of larval control in Cape Coast, Ghana.

Methods

A larval survey was conducted for 15 months in Cape Coast. In individual breeding habitats, habitat characteristics, physicochemical parameters and bacterial fauna were measured in both Anopheles positive breeding (APL) habitats and habitats colonized by only Culex species. The sibling species of An. gambiae were identified using PCR assay.

Results

Anopheles coluzzii dominated in almost all the APL habitats found in this study. The habitats had high levels of salinity and ammonium ions. However, ammonium ions were significantly higher (p = 0.001) in habitats colonized by only Culex larvae compared to APL habitats. About 47 % of the habitats that were colonized by only Culex larvae had no measurable dissolved oxygen while An. coluzzii was absent in such habitats. High concentration of faecal bacteria confirmed faecal contamination in both groups of breeding habitats.

Conclusions

From the results, it was evident that larval stages of An. coluzzii have tolerance to high levels of salinity and organic pollution in breeding habitats. However, its level of tolerance to organic pollution is probably lower than Culex larvae. The nature of breeding habitats found in the city demonstrates the opportunistic behaviour of An. coluzzii and how its breeding requirements are so intimately intertwined with the haphazard and uncontrolled human activities in the urban area. Considering the nature of APL habitats, larval control intervention could greatly reduce Anopheles population. However, improving basic hygiene and sanitation in the city could even make larval control intervention more practical and cost effective.

Population Genetic Structure of the Malaria Vector Anopheles sinensis (Diptera: Culicidae) Sensu Stricto and Evidence for Possible Introgression in the Republic of Korea

Anopheles sinensis Wiedemann sensu stricto (s.s.) is a dominant mosquito and considered a secondary malaria vector in the Republic of Korea (ROK). Despite the potential significance for malaria control, population genetics studies have been conducted using only mitochondrial DNA (mtDNA), and studies of the genetics of hybridization have never been attempted. In this study, 346 specimens from 23 localities were subject to experiments. Among them, 305 An. sinensis s.s. specimens from 20 localities were used for mtDNA analysis, and 346 specimens comprising 341 An. sinensis s.s. from 22 localities and five Anopheles kleini Rueda from one locality were examined in the microsatellite study. Neighbor-joining analysis of pairwise FST and RST based on microsatellite results showed that the populations are divided into two groups, as did the mtDNA results. However, the Bayesian analysis and factorial correspondence analysis plots showed three distinct clusters. Among the mtDNA and microsatellite results, only microsatellites represented small but positive and significant isolation-by-distance patterns. Both molecular markers show the Taebaek and Sobaek Mountain ranges as barriers between the northern and southern parts of the ROK. The newly recognized third group suggests possible introgressive hybridization of An. sinensis s.s. with closely related species. The slightly different composition of populations in each group based on different markers is probably because of different population dynamics in each group. These results imply that there is restricted gene flow of epidemiologically important malaria-related genes between the northern and southern parts of the ROK.

Habitat segregation and ecological character displacement in cryptic African malaria mosquitoes

Understanding how divergent selection generates adaptive phenotypic and population diversification provides a mechanistic explanation of speciation in recently separated species pairs. Towards this goal, we sought ecological gradients of divergence between the cryptic malaria vectors Anopheles coluzzii and An. gambiae and then looked for a physiological trait that may underlie such divergence. Using a large set of occurrence records and eco-geographic information, we built a distribution model to predict the predominance of the two species across their range of sympatry. Our model predicts two novel gradients along which the species segregate: distance from the coastline and altitude. Anopheles coluzzii showed a ‘bimodal’ distribution, predominating in xeric West African savannas and along the western coastal fringe of Africa. To test whether differences in salinity tolerance underlie this habitat segregation, we assessed the acute dose–mortality response to salinity of thirty-two larval populations from Central Africa. In agreement with its coastal predominance, Anopheles coluzzii was overall more tolerant than An. gambiae. Salinity tolerance of both species, however, converged in urban localities, presumably reflecting an adaptive response to osmotic stress from anthropogenic pollutants. When comparing degree of tolerance in conjunction with levels of syntopy, we found evidence of character displacement in this trait.

Seasonal variation in wing size and shape between geographic populations of the malaria vector, Anopheles coluzzii in Burkina Faso (West Africa)

The mosquito, Anopheles coluzzii is a major vector of human malaria in Africa with widespread distribution throughout the continent. The species hence populates a wide range of environments in contrasted ecological settings often exposed to strong seasonal fluctuations. In the dry savannahs of West Africa, this mosquito population dynamics closely follows the pace of surface water availability: the species pullulates during the rainy season and is able to reproduce throughout the dry season in areas where permanent water bodies are available for breeding. The impact of such environmental fluctuation on mosquito development and the phenotypic quality of emerging adults has however not been addressed in details. Here, we examined and compared phenotypic changes in the duration of pre-imaginal development, body dry mass at emergence and wing size, shape and surface area in young adult females An. coluzzii originated from five distinct geographic locations when they are reared in two contrasting conditions mimicking those experienced by mosquitoes during the rainy season (RS) and at the onset of the dry season (ODS) in Burkina Faso (West Africa). Our results demonstrated strong phenotypic plasticity in all traits, with differences in the magnitude and direction of changes between RS and ODS depending upon the geographic origin, hence the genetic background of the mosquito populations. Highest heterogeneity within population was observed in Bama, where large irrigation schemes allow year-round mosquito breeding. Further studies are needed to explore the adaptive value of such phenotypic plasticity and its relevance for local adaptation in An. coluzzii.

Mosquito genomics. Extensive introgression in a malaria vector species complex revealed by phylogenomics

Introgressive hybridization is now recognized as a widespread phenomenon, but its role in evolution remains contested. Here, we use newly available reference genome assemblies to investigate phylogenetic relationships and introgression in a medically important group of Afrotropical mosquito sibling species. We have identified the correct species branching order to resolve a contentious phylogeny and show that lineages leading to the principal vectors of human malaria were among the first to split. Pervasive autosomal introgression between these malaria vectors means that only a small fraction of the genome, mainly on the X chromosome, has not crossed species boundaries. Our results suggest that traits enhancing vectorial capacity may be gained through interspecific gene flow, including between nonsister species.

VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases

VectorBase is a National Institute of Allergy and Infectious Diseases supported Bioinformatics Resource Center (BRC) for invertebrate vectors of human pathogens. Now in its 11th year, VectorBase currently hosts the genomes of 35 organisms including a number of non-vectors for comparative analysis. Hosted data range from genome assemblies with annotated gene features, transcript and protein expression data to population genetics including variation and insecticide-resistance phenotypes. Here we describe improvements to our resource and the set of tools available for interrogating and accessing BRC data including the integration of Web Apollo to facilitate community annotation and providing Galaxy to support user-based workflows. VectorBase also actively supports our community through hands-on workshops and online tutorials. All information and data are freely available from our website at https://www.vectorbase.org/.

Distinct population structure for co-occurring Anopheles goeldii and Anopheles triannulatus in Amazonian Brazil

To evaluate whether environmental heterogeneity contributes to the genetic heterogeneity in Anopheles triannulatus, larval habitat characteristics across the Brazilian states of Roraima and Pará and genetic sequences were examined. A comparison with Anopheles goeldii was utilised to determine whether high genetic diversity was unique to An. triannulatus. Student t test and analysis of variance found no differences in habitat characteristics between the species. Analysis of population structure of An. triannulatus and An. goeldii revealed distinct demographic histories in a largely overlapping geographic range. Cytochrome oxidase I sequence parsimony networks found geographic clustering for both species; however nuclear marker networks depicted An. triannulatus with a more complex history of fragmentation, secondary contact and recent divergence. Evidence of Pleistocene expansions suggests both species are more likely to be genetically structured by geographic and ecological barriers than demography. We hypothesise that niche partitioning is a driving force for diversity, particularly in An. triannulatus.

Genome-wide patterns of polymorphism in an inbred line of the African malaria mosquito Anopheles gambiae

Anopheles gambiae is a major mosquito vector of malaria in Africa. Although increased use of insecticide-based vector control tools has decreased malaria transmission, elimination is likely to require novel genetic control strategies. It can be argued that the absence of an A. gambiae inbred line has slowed progress toward genetic vector control. In order to empower genetic studies and enable precise and reproducible experimentation, we set out to create an inbred line of this species. We found that amenability to inbreeding varied between populations of A. gambiae. After full-sib inbreeding for ten generations, we genotyped 112 individuals--56 saved prior to inbreeding and 56 collected after inbreeding--at a genome-wide panel of single nucleotide polymorphisms (SNPs). Although inbreeding dramatically reduced diversity across much of the genome, we discovered numerous, discrete genomic blocks that maintained high heterozygosity. For one large genomic region, we were able to definitively show that high diversity is due to the persistent polymorphism of a chromosomal inversion. Inbred lines in other eukaryotes often exhibit a qualitatively similar retention of polymorphism when typed at a small number of markers. Our whole-genome SNP data provide the first strong, empirical evidence supporting associative overdominance as the mechanism maintaining higher than expected diversity in inbred lines. Although creation of A. gambiae lines devoid of nearly all polymorphism may not be feasible, our results provide critical insights into how more fully isogenic lines can be created.

Islands and stepping-stones: comparative population structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and implications for the spread of insecticide resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise FST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average FST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.

Composition of Anopheles coluzzii and Anopheles gambiae microbiota from larval to adult stages

During their immature life stages, malaria mosquitoes are exposed to a wide array of microbes and contaminants from the aquatic habitats. Although prior studies have suggested that environmental exposure shapes the microbial community structure in the adult mosquito, most reports have focused on laboratory-based experiments and on a single mosquito epithelium, the gut. In this study, we investigated the influence of the breeding site on the development of the Anopheles coluzzii and Anopheles gambiae microbiota in natural conditions. We characterized bacterial communities from aquatic habitats, at surface microlayer and subsurface water levels, to freshly emerge adult mosquitoes using multiplexed 16S rRNA gene pyrosequencing and we separately analyzed the microbiota associated with the different epithelia of adult individual, midguts, ovaries and salivary glands. We found that the distribution of bacterial communities in the aquatic habitats differed according to the depth of water collections. Inter-individual variation of bacterial composition was large in larvae guts but adult mosquitoes from a same breeding site shared quite similar microbiota. Although some differences in bacterial abundances were highlighted between the different epithelia of freshly emerged An. coluzzii and An. gambiae, an intriguing feature from our study is the particular similarity of the overall bacterial communities. Our results call for further investigations on the bacterial population dynamics in the different tissues to determine the distinctive characteristics of each microbiota during the mosquito lifespan and to identify specific interactions between certain key phyla or species and the insect life history traits.

Molecular characterization of DDT resistance in Anopheles gambiae from Benin

BACKGROUND

Insecticide resistance in the mosquito vector is the one of the main obstacles against effective malaria control. In order to implement insecticide resistance management strategies, it is important to understand the genetic factors involved. In this context, we investigated the molecular basis of DDT resistance in the main malaria vector from Benin.

METHODS

Anopheles gambiae mosquitoes were collected from four sites across Benin and identified to species/molecular form. Mosquitoes from Cotonou (M-form), Tori-Bossito (S-form) and Bohicon (S-form) were exposed to DDT 4% at a range of exposure times (30 min to 300 min). Another batch of mosquitoes from Cotonou and Malanville were exposed to DDT for 1 hour and the survivors 48 hours post exposure were used to quantify metabolic gene expression. Quantitative PCR assays were used to quantify mRNA levels of metabolic enzymes: GSTE2, GSTD3, CYP6P3 and CYP6M2. Expression (fold-change) was calculated using the ∆∆Ct method and compared to susceptible strains. Detection of target-site mutations (L1014F, L1014S and N1575Y) was performed using allelic discrimination TaqMan assays.

RESULTS

DDT resistance was extremely high in all populations, regardless of molecular form, with no observed mortality after 300 min exposure. In both DDT-survivors and non-exposed mosquitoes, GSTE2 and GSTD3 were over-expressed in the M form at 4.4-fold and 3.5-fold in Cotonou and 1.5-fold and 2.5-fold in Malanville respectively, when compared to the susceptible strain. The CYP6M2 and CYP6P3 were over-expressed at 4.6-fold and 3.8-fold in Cotonou and 1.2-fold and 2.5-fold in Malanville respectively. In contrast, no differences in GSTE2 and CYP6M2 were observed between S form mosquitoes from Tori-Bossito and Bohicon compared to susceptible strain. The 1014 F allele was fixed in the S-form and at high frequency in the M-form (0.7-0.914). The frequency of 1575Y allele was 0.29-0.36 in the S-form and nil in the M-form. The 1014S allele was detected in the S form of An. gambiae in a 1014 F/1014S heterozygous specimen.

CONCLUSION

Our results show that the kdr 1014 F, 1014S and 1575Y alleles are widespread in Benin and the expression of two candidate metabolic markers (GSTE2 and CYP6M2) are over-expressed specifically in the M-form.

Salivary gland proteome analysis reveals modulation of anopheline unique proteins in insensitive acetylcholinesterase resistant Anopheles gambiae mosquitoes

Insensitive acetylcholinesterase resistance due to a mutation in the acetylcholinesterase (ace) encoding ace-1 gene confers cross-resistance to organophosphate and carbamate insecticides in Anopheles gambiae populations from Central and West Africa. This mutation is associated with a strong genetic cost revealed through alterations of some life history traits but little is known about the physiological and behavioural changes in insects bearing the ace-1^R allele. Comparative analysis of the salivary gland contents between An. gambiae susceptible and ace-1^R resistant strains was carried out to charaterize factors that could be involved in modifications of blood meal process, trophic behaviour or pathogen interaction in the insecticide-resistant mosquitoes. Differential analysis of the salivary gland protein profiles revealed differences in abundance for several proteins, two of them showing major differences between the two strains. These two proteins identified as saglin and TRIO are salivary gland-1 related proteins, a family unique to anopheline mosquitoes, one of them playing a crucial role in salivary gland invasion by Plasmodium falciparum sporozoites. Differential expression of two other proteins previously identified in the Anopheles sialome was also observed. The differentially regulated proteins are involved in pathogen invasion, blood feeding process, and protection against oxidation, relevant steps in the outcome of malaria infection. Further functional studies and insect behaviour experiments would confirm the impact of the modification of the sialome composition on blood feeding and pathogen transmission abilities of the resistant mosquitoes. The data supports the hypothesis of alterations linked to insecticide resistance in the biology of the primary vector of human malaria in Africa.

Reproductive isolation among sympatric molecular forms of An. gambiae from inland areas of south-eastern Senegal

The Anopheles gambiae species complex includes at least seven morphologically indistinguishable species, one of which, Anopheles gambiae sensu stricto, is the primary mosquito vector responsible for the transmission of malaria across sub-Saharan Africa. Sympatric ecological diversification of An. gambiae s.s. is in progress within this complex, leading to the emergence of at least two incipient species (the M and S molecular forms now recognized as good species and named An. coluzzii and An. gambiae respectively) that show heterogeneous levels of divergence in most parts of Africa. However, this process seems to have broken down in coastal areas of West Africa at the extreme edge of the distribution. We undertook a longitudinal study to describe An. gambiae s.s. populations collected from two inland transects with different ecological characteristics in south-eastern Senegal. Analysis of samples collected from 20 sites across these two transects showed the M and S molecular forms coexisted at almost all sampled sites. Overall, similar hybridization rates (2.16% and 1.86%) were recorded in the two transects; sites with relatively high frequencies of M/S hybrids (up to 7%) were clustered toward the north-western part of both transects, often near urban settings. Estimated inbreeding indices for this putative speciation event varied spatially (range: 0.52-1), with hybridization rates being generally lower than expected under panmictic conditions. Such observations suggest substantial reproductive isolation between the M and S molecular forms, and further support the ongoing process of speciation in these inland areas. According to a recent reclassification of the An. gambiae complex, the M and S molecular forms from this zone correspond to An. coluzzii and An. gambiae, respectively. There is considerable evidence that these molecular forms differ in their behavioural and ecological characteristics. Detailed study of these characteristics will allow the development and implementation of better insect control strategies for combating malaria.

Contemporary gene flow between wild An. gambiae s.s. and An. arabiensis

Background

In areas where the morphologically indistinguishable malaria mosquitoes Anopheles gambiae Giles and An. arabiensis Patton are sympatric, hybrids are detected occasionally via species-diagnostic molecular assays. An. gambiae and An. arabiensis exhibit both pre- and post-reproductive mating barriers, with swarms largely species-specific and male F1 (first-generation) hybrids sterile. Consequently advanced-stage hybrids (back-crosses to parental species), which would represent a route for potentially-adaptive introgression, are expected to be very rare in natural populations. Yet the use of one or two physically linked single-locus diagnostic assays renders them indistinguishable from F1 hybrids and levels of interspecific gene flow are unknown.

Methods

We used data from over 350 polymorphic autosomal SNPs to investigate post F1 gene flow via patterns of genomic admixture between An. gambiae and An. arabiensis from eastern Uganda. Simulations were used to investigate the statistical power to detect hybrids with different levels of crossing and to identify the hybrid category significantly admixed genotypes could represent.

Results

A range of admixture proportions were detected for 11 field-collected hybrids identified via single-locus species-diagnostic PCRs. Comparison of admixture data with simulations indicated that at least seven of these hybrids were advanced generation crosses, with backcrosses to each species identified. In addition, of 36 individuals typing as An. gambiae or An. arabiensis that exhibited outlying admixture proportions, ten were identified as significantly mixed backcrosses, and at least four of these were second or third generation crosses.

Conclusions

Our results show that hybrids detected using standard diagnostics will often be hybrid generations beyond F1, and that in our study area around 5% (95% confidence intervals 3%-9%) of apparently ‘pure’ species samples may also be backcrosses. This is likely an underestimate because of rapidly-declining detection power beyond the first two backcross generations. Post-F1 gene flow occurs at a far from inconsequential rate between An. gambiae and An. arabiensis, and, especially for traits under strong selection, could readily lead to adaptive introgression of genetic variants relevant for vector control.

Electronic supplementary material

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

Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation

Adaptive introgression can provide novel genetic variation to fuel rapid evolutionary responses, though it may be counterbalanced by potential for detrimental disruption of the recipient genomic background. We examine the extent and impact of recent introgression of a strongly selected insecticide-resistance mutation (Vgsc-1014F) located within one of two exceptionally large genomic islands of divergence separating the Anopheles gambiae species pair. Here we show that transfer of the Vgsc mutation results in homogenization of the entire genomic island region (~1.5% of the genome) between species. Despite this massive disruption, introgression is clearly adaptive with a dramatic rise in frequency of Vgsc-1014F and no discernable impact on subsequent reproductive isolation between species. Our results show (1) how resilience of genomes to massive introgression can permit rapid adaptive response to anthropogenic selection and (2) that even extreme prominence of genomic islands of divergence can be an unreliable indicator of importance in speciation.

Highly divergent genomic islands segregate between a species pair of the mosquito, Anopheles gambiae. Here Clarkson et al. show that loss of one of the largest islands, driven by adaptive introgression of an insecticide-resistance mutation, had no impact on reproductive isolation.

Gene expression divergence between malaria vector sibling species Anopheles gambiae and An. coluzzii from rural and urban Yaoundé Cameroon

Divergent selection based on aquatic larval ecology is a likely factor in the recent isolation of two broadly sympatric and morphologically identical African mosquito species, the malaria vectors Anopheles gambiae and An. coluzzii. Population-based genome scans have revealed numerous candidate regions of recent positive selection, but have provided few clues as to the genetic mechanisms underlying behavioural and physiological divergence between the two species, phenotypes which themselves remain obscure. To uncover possible genetic mechanisms, we compared global transcriptional profiles of natural and experimental populations using gene-based microarrays. Larvae were sampled as second and fourth instars from natural populations in and around the city of Yaoundé, capital of Cameroon, where the two species segregate along a gradient of urbanization. Functional enrichment analysis of differentially expressed genes revealed that An. coluzzii--the species that breeds in more stable, biotically complex and potentially polluted urban water bodies--overexpresses genes implicated in detoxification and immunity relative to An. gambiae, which breeds in more ephemeral and relatively depauperate pools and puddles in suburbs and rural areas. Moreover, our data suggest that such overexpression by An. coluzzii is not a transient result of induction by xenobiotics in the larval habitat, but an inherent and presumably adaptive response to repeatedly encountered environmental stressors. Finally, we find no significant overlap between the differentially expressed loci and previously identified genomic regions of recent positive selection, suggesting that transcriptome divergence is regulated by trans-acting factors rather than cis-acting elements.

Cuticular differences associated with aridity acclimation in African malaria vectors carrying alternative arrangements of inversion 2La

BACKGROUND

Principal malaria vectors in Africa, An. gambiae and An. coluzzii, share an inversion polymorphism on the left arm of chromosome 2 (2La/2L+a) that is distributed non-randomly in the environment. Genomic sequencing studies support the role of strong natural selection in maintaining steep clines in 2La inversion frequency along environmental gradients of aridity, and physiological studies have directly implicated 2La in heat and desiccation tolerance, but the precise genetic basis and the underlying behavioral and physiological mechanisms remain unknown. As the insect cuticle is the primary barrier to water loss, differences in cuticle thickness and/or epicuticular waterproofing associated with alternative 2La arrangements might help explain differences in desiccation resistance.

METHODS

To test that hypothesis, two subcolonies of both An. gambiae and An. coluzzii were established that were fixed for alternative 2La arrangements (2La or 2L+a) on an otherwise homosequential and shared genetic background. Adult mosquitoes reared under controlled environmental conditions (benign or arid) for eight days post-eclosion were collected and analyzed. Measurements of cuticle thickness were made based on scanning electron microscopy, and cuticular hydrocarbon (CHC) composition was evaluated by gas chromatography-mass spectrometry.

RESULTS

After removing the allometric effects of body weight, differences in mean cuticle thickness were found between alternative 2La karyotypes, but not between alternative environments. Moreover, the thicker cuticle of the An. coluzzii 2La karyotype was contrary to the known higher rate of water loss of this karyotype relative to 2L+a. On the other hand, quantitative differences in individual CHCs and overall CHC profiles between alternative karyotypes and environmental conditions were consistent with expectation based on previous physiological studies.

CONCLUSIONS

Our results suggest that alternative arrangements of the 2La inversion are associated with differences in cuticle thickness and CHC composition, but that only CHC composition appears to be relevant for desiccation resistance. Differences in the CHC composition were consistent with previous findings of a lower rate of water loss for the 2L+a karyotype at eight days post-eclosion, suggesting that CHC composition is an important strategy for maintaining water balance in this genetic background, but not for 2La. Despite a higher rate of water loss at eight days, higher body water content of the 2La karyotype confers a level of desiccation resistance equivalent to that of the 2L+a karyotype.

Metabolomic and ecdysteroid variations in Anopheles gambiae s.l. mosquitoes exposed to the stressful conditions of the dry season in Burkina Faso, West Africa

This study explored the metabolic adjustments prompted by a switch between the rainy and dry season conditions in the African malaria mosquitoes Anopheles gambiae (M and S molecular forms) and Anopheles arabiensis. Mosquitoes were reared in contrasted experimental conditions reflecting environmental variation in Burkina Faso. Thirty-five metabolites (including sugars, polyols, and amino acids) were monitored in newly emerged males and females, and their ecdysteroid titers were determined. Metabolomic signatures were remarkably similar across species, when specimens of same age and sex were reared under identical experimental conditions. In males and females, amino acids (including glycine, leucine, phenylanine, serine, threonine, and valine) were accumulated in 1-h-old mosquitoes, then decreased 24 h after emergence, probably reflecting adult maturation and the amino acid-consuming process of cuticle sclerotisation. In turn, elevated amounts of alanine and proline in 24-h-old mosquitoes may assist the development of flight ability. Lower concentration of tricarboxylic acid cycle intermediates and isoleucine characterized older females reared under dry season conditions, suggesting metabolic and reproduction depression. In all cases, ecdysteroid concentration was much higher in males than in females, with significant seasonal variation in males. This might reflect a unique role of these hormones in shaping reproductive strategies and population demography in the An. gambiae s.l. species complex, further contributing to local adaptation in a highly fluctuating environment.

Microdistribution of the resistance of malaria vectors to deltamethrin in the region of Plateau (southeastern Benin) in preparation for an assessment of the impact of resistance on the effectiveness of Long Lasting Insecticidal Nets (LLINs)

BACKGROUND

This study aims to research two areas, one with a resistant and the other with a susceptible profile of An. gambiae to deltamethrin in the region of Plateau (southern Benin). In each area, eight localities were sought. Both areas were needed for the assessment of the impact of malaria vector resistance to pyrethroids on the effectiveness of Long Lasting Insecticidal Nets (LLINs). The susceptible area of An. gambiae to deltamethrin was used as a control.

METHODS

In total, 119 localities in the region of Plateau were screened by sampling An. gambiae s.l larvae. Female mosquitoes resulting from these larvae were exposed to 0.05% deltamethrin following WHO standards. PCR was used to identify species and molecular forms of the dead and alive mosquitoes. Finally, we identified kdr mutations (1014 F and 1014S) using the HOLA technique.

RESULTS

Fifty-six out of 119 prospected localities tested positive for Anopheles gambae s.l breeding sites. The results showed that An. gambiae was resistant to deltamethrin in 39 localities and susceptible in only 2 localities; resistance to deltamethrin was suspected in 15 localities. The HOLA technique confirmed the presence of kdr 1014 F mutation and the absence of kdr 1014S mutation. The kdr 1014 F mutation was found in both M and S molecular forms at relatively high frequencies therefore confirming the susceptibility tests.

CONCLUSION

We were unable to identify the eight susceptible areas due to the overall resistance of An. gambiae to deltamethrin in the region of Plateau. To implement the study, we kept two areas, one with high resistance (R+++) and the other with low resistance (R+) of An. gambiae to deltamethrin.

Larval competition between An. coluzzii and An. gambiae in insectary and semi-field conditions in Burkina Faso

Competition in mosquito larvae is common and different ecological context could change competitive advantage between species. Here, larval competition between the widely sympatric African malaria mosquitoes, Anopheles coluzzii and Anopheles gambiae were investigated in controlled insectary conditions using individuals from laboratory colonies and under ambient conditions using wild mosquitoes in a semi-field enclosure in western Burkina Faso. Larvae of both species were reared in trays at the same larval density and under the same feeding regimen in either single-species or mixed-species populations at varying species ratios reflecting 0%, 25%, 50% and 75% of competitor species. In the insectaries, where environmental variations are controlled, larvae of the An. coluzzii colony developed faster and with lower mortality than larvae of the An. gambiae colony (8.8±0.1 days and 21±3% mortality vs. 9.5±0.1 days and 32±3% mortality, respectively). Although there was no significant effect of competition on these phenotypic traits in any species, there was a significant trend for higher fitness of the An. coluzzii colony when competing with An. gambiae under laboratory conditions (i.e. lower development time and increased wing length at emergence, Cuzik's tests, P<0.05). In semi-field experiments, competition affected the life history traits of both species in a different way. Larvae of An. gambiae tended to reduce development time when in competition with An. coluzzii (Cuzick's test, P=0.002) with no impact either on mortality or size at emergence. On the other hand, An. coluzzii showed a significant trend for reduced larval mortality with increasing competition pressure (Cuzick's test, P=0.037) and production of smaller females when grown together with An. gambiae (Cuzick's test, P=0.002). Our results hence revealed that competitive interactions between larvae of the two species are context dependent. They further call for caution when exploring ecological processes using inbred laboratory colonies in this system of utmost medical importance.

Divergence in threat sensitivity among aquatic larvae of cryptic mosquito species

Predation is a major evolutionary force driving speciation. The threat-sensitive response hypothesis predicts that prey adjust and balance the time spent on a costly antipredator response with other activities that enhance their fitness. Thus, prey able to develop an antipredator response proportional to risk intensity should have a selective advantage. Knowledge on how evolution has shaped threat sensitivity among closely related species exposed to different predation pressures is scarce, prompting investigations to better predict and explain its effect on communities. We explored and compared the antipredator response of aquatic mosquito larvae in three sibling species of the Anopheles gambiae complex, with contrasting larval biologies in Burkina Faso. Anopheles arabiensis and An. gambiae sensu stricto breed in temporary water collections where predator densities are low, whereas Anopheles coluzzii is able to thrive in permanent pools where the predation pressure is much higher. We hypothesized that the increase and decline of behavioural antipredator responses might differ between the three species over time. To test this hypothesis, progenies of field-collected mosquitoes were experimentally exposed to a range of soluble predation cues and their response was monitored for up to 48 h. The three species were all threat sensitive but their reaction norms differed. For the range of concentrations tested, An. coluzzii larvae gradually increased in antipredator response, whereas An. gambiae larvae readily displayed antipredator behaviour at low concentrations leading to a saturation of the response for high cue concentrations. An. arabiensis displayed a narrower reaction norm with low response intensity. Larval instars did not differ in their threat sensitivity. The antipredator behaviour of the three species waned after about 1 h of exposure. Early instars tended to express antipredation behaviour for longer than did older instars. This study provides information on how aquatic prey species with an aerial adult stage manage larval predation risk over time according to cue concentrations and suggests that different predation pressures might play a role as a disruptive selective force fostering habitat segregation and speciation within the An. gambiae complex. The evolution of phenotypic plasticity is further discussed in the light of divergent predation pressures.

Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae

The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An "islands of speciation" model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative "incidental island" model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A "divergence island SNP" assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.

Acetylcholinesterase (Ace-1) target site mutation 119S is strongly diagnostic of carbamate and organophosphate resistance in Anopheles gambiae s.s. and Anopheles coluzzii across southern Ghana

BACKGROUND

With high DDT resistance present throughout much of West Africa, carbamates and organophosphates are increasingly important alternatives to pyrethroids for indoor residual spraying (IRS). Though less widespread, resistance to both of these alternative insecticide classes has also been documented within the Anopheles gambiae species pair (formerly the M and S molecular forms) in West Africa. To manage insecticide efficacy, it is important to predict how and where resistance is likely to occur and spread, which could be aided by using molecular diagnostics with high predictive value.

METHODS

Anopheles coluzzii and An. gambiae s.s. were collected from 18 sites throughout southern Ghana and bioassayed with bendiocarb, the most commonly applied carbamate, and an organophosphate, fenitrothion. The Ace-1 target site substitution G119S was genotyped by qPCR.

RESULTS

Fenitrothion induced higher mortality than bendiocarb, though phenotypes correlated strongly across populations. Ace-1 119S was found at much higher frequency in An. gambiae s.s than An. coluzzii, exceeding 90% in a population from Greater Accra, the highest frequency reported to date. Ace-1 G119S was very strongly associated with resistance to both insecticides, providing high predictive power for diagnosis, though with some evidence for a differential effect between molecular forms for bendiocarb. Sequencing of the gene revealed a lack of variation in resistant alleles precluding determination of origin, but Ace-1 copy number variation was detected for the first time in Ghana.

CONCLUSIONS

The results validate G119S as a useful diagnostic of organophosphate and carbamate resistance within and among populations, whilst highlighting the potential for an aggregate nature of Ace-1 genotypes, which may comprise both single-copy and duplicated genes. Further work is now required to determine the distribution and resistance-association of Ace-1 duplication.

Insect speciation rules: unifying concepts in speciation research

The study of speciation is concerned with understanding the connection between causes of divergent evolution and the origin and maintenance of barriers to gene exchange between incipient species. Although the field has historically focused either on examples of recent divergence and its causes or on the genetic basis of reproductive isolation between already divergent species, current efforts seek to unify these two approaches. Here we integrate these perspectives through a discussion of recent progress in several insect speciation model systems. We focus on the evolution of speciation phenotypes in each system (i.e., those phenotypes causally involved in reducing gene flow between incipient species), drawing an explicit connection between cause and effect (process and pattern). We emphasize emerging insights into the genomic architecture of speciation as well as timely areas for future research.

Assortative mating in mixed swarms of the mosquito Anopheles gambiae s.s. M and S molecular forms, in Burkina Faso, West Africa

The molecular form composition of Anopheles gambiae Giles s.s. (Diptera: Culicidae) mating swarms and the associated mating pairs (copulae) were investigated during two rainy seasons (July to October, 2005 and July to November, 2006) in the villages of Soumousso and Vallée du Kou (VK7). Although the habitats of these villages differ markedly, sympatric populations of M and S molecular forms of An. gambiae s.s. occur in both places periodically. The main aim was to assess the degree to which these molecular forms mate assortatively. In Soumousso, a wooded savannah habitat, the majority of swarm samples consisted of only S-form males (21/28), although a few M-form males were found in mixed M- and S-form swarms. In VK7, a rice growing area, the majority of swarm samples consisted of only M-form males (38/62), until October and November 2006, when there were nearly as many mixed-form as single-form swarms. Overall, ∼60% of M- and S-form swarms were temporally or spatially segregated; the two forms were effectively prevented from encountering each other. Of the remaining 40% of swarms, however, only about half were single-form and the rest were mixed-form. Of the 33 copulae collected from mixed-form swarms, only four were mixed-form pairs, significantly fewer than expected by random pairing between forms (χ(2) = 10.34, d.f. = 2, P < 0.01). Finally, all specimens of inseminated females were of the same form as the sperm contained within their spermatheca (n = 91), even for the four mixed-form copulae. These findings indicate that assortative mating occurs within mixed-form swarms, mediated most probably by close-range mate recognition cues.

Molecular characterization and identification of members of the Anopheles subpictus complex in Sri Lanka

BACKGROUND

Anopheles subpictus sensu lato is a major malaria vector in South and Southeast Asia. Based initially on polytene chromosome inversion polymorphism, and subsequently on morphological characterization, four sibling species A-D were reported from India. The present study uses molecular methods to further characterize and identify sibling species in Sri Lanka.

METHODS

Mosquitoes from Sri Lanka were morphologically identified to species and sequenced for the ribosomal internal transcribed spacer-2 (ITS2) and the mitochondrial cytochrome c oxidase subunit-I (COI) genes. These sequences, together with others from GenBank, were used to construct phylogenetic trees and parsimony haplotype networks and to test for genetic population structure.

RESULTS

Both ITS2 and COI sequences revealed two divergent clades indicating that the Subpictus complex in Sri Lanka is composed of two genetically distinct species that correspond to species A and species B from India. Phylogenetic analysis showed that species A and species B do not form a monophyletic clade but instead share genetic similarity with Anopheles vagus and Anopheles sundaicus s.l., respectively. An allele specific identification method based on ITS2 variation was developed for the reliable identification of species A and B in Sri Lanka.

CONCLUSION

Further multidisciplinary studies are needed to establish the species status of all chromosomal forms in the Subpictus complex. This study emphasizes the difficulties in using morphological characters for species identification in An. subpictus s.l. in Sri Lanka and demonstrates the utility of an allele specific identification method that can be used to characterize the differential bio-ecological traits of species A and B in Sri Lanka.

A preliminary investigation of the relationship between water quality and Anopheles gambiae larval habitats in Western Cameroon

Background

Water quality and anopheline habitat have received increasing attention due to the possibility that challenges during larval life may translate into adult susceptibility to malaria parasite infection and/or insecticide resistance.

Methods

A preliminary study of Anopheles gambiae s.s. larval habitats in the north-west and south-west regions of Cameroon was conducted in order to detect associations between An. gambiae s.s. molecular form and 2La inversion distributions with basic water quality parameters. Water quality was measured by temperature, pH, conductivity, total dissolved solids (TDS) at seven sites in Cameroon and one site in Selinkenyi, Mali.

Results

Principal components and correlation analyses indicated a complex relationship between 2La polymorphism, temperature, conductivity and TDS. Cooler water sites at more inland locations yielded more S form larvae with higher 2La inversion polymorphism while warmer water sites yielded more M form larvae with rare observations of the 2La inversion.

Discussion

More detailed studies that take into account the population genetics but also multiple life stages, environmental data relative to these life stages and interactions with both humans and the malaria parasite may help us to understand more about how and why this successful mosquito is able to adapt and diverge, and how it can be successfully managed.

Equivalent susceptibility of Anopheles gambiae M and S molecular forms and Anopheles arabiensis to Plasmodium falciparum infection in Burkina Faso

Background

The Anopheles gambiae sensu lato (s.l.) species complex in Burkina Faso consists of Anopheles arabiensis, and molecular forms M and S of Anopheles gambiae sensu stricto (s.s.). Previous studies comparing the M and S forms for level of infection with Plasmodium falciparum have yielded conflicting results.

Methods

Mosquito larvae were sampled from natural pools, reared to adulthood under controlled conditions, and challenged with natural P. falciparum by experimental feeding with blood from gametocyte carriers. Oocyst infection prevalence and intensity was determined one week after infection. DNA from carcasses was genotyped to identify species and molecular form.

Results

In total, 7,400 adult mosquitoes grown from wild-caught larvae were challenged with gametocytes in 29 experimental infections spanning four transmission seasons. The overall infection prevalence averaged 40.7% for A. gambiae M form, 41.4% for A. gambiae S form, and 40.1% for A. arabiensis. There was no significant difference in infection prevalence or intensity between the three population groups. Notably, infection experiments in which the population groups were challenged in parallel on the same infective blood displayed less infection difference between population groups, while infections with less balanced composition of population groups had lower statistical power and displayed apparent differences that fluctuated more often from the null average.

Conclusion

The study clearly establishes that, at the study site in Burkina Faso, there is no difference in genetic susceptibility to P. falciparum infection between three sympatric population groups of the A. gambiae s.l. complex. Feeding the mosquito groups on the same infective blood meal greatly increases statistical power. Conversely, comparison of the different mosquito groups between, rather than within, infections yields larger apparent difference between mosquito groups, resulting from lower statistical power and greater noise, and could lead to false-positive results. In making infection comparisons between population groups, it is more accurate to compare the different groups after feeding simultaneously upon the same infective blood.

Geographic population structure of the African malaria vector Anopheles gambiae suggests a role for the forest-savannah biome transition as a barrier to gene flow

The primary Afrotropical malaria mosquito vector Anopheles gambiae sensu stricto has a complex population structure. In west Africa, this species is split into two molecular forms and displays local and regional variation in chromosomal arrangements and behaviors. To investigate patterns of macrogeographic population substructure, 25 An. gambiae samples from 12 African countries were genotyped at 13 microsatellite loci. This analysis detected the presence of additional population structuring, with the M-form being subdivided into distinct west, central, and southern African genetic clusters. These clusters are coincident with the central African rainforest belt and northern and southern savannah biomes, which suggests restrictions to gene flow associated with the transition between these biomes. By contrast, geographically patterned population substructure appears much weaker within the S-form.

Insecticide resistance alleles affect vector competence of Anopheles gambiae s.s. for Plasmodium falciparum field isolates

The widespread insecticide resistance raises concerns for vector control implementation and sustainability particularly for the control of the main vector of human malaria, Anopheles gambiae sensu stricto. However, the extent to which insecticide resistance mechanisms interfere with the development of the malignant malaria parasite in its vector and their impact on overall malaria transmission remains unknown. We explore the impact of insecticide resistance on the outcome of Plasmodium falciparum infection in its natural vector using three An. gambiae strains sharing a common genetic background, one susceptible to insecticides and two resistant, one homozygous for the ace-1(R) mutation and one for the kdr mutation. Experimental infections of the three strains were conducted in parallel with field isolates of P. falciparum from Burkina Faso (West Africa) by direct membrane feeding assays. Both insecticide resistant mutations influence the outcome of malaria infection by increasing the prevalence of infection. In contrast, the kdr resistant allele is associated with reduced parasite burden in infected individuals at the oocyst stage, when compared to the susceptible strain, while the ace-1 (R) resistant allele showing no such association. Thus insecticide resistance, which is particularly problematic for malaria control efforts, impacts vector competence towards P. falciparum and probably parasite transmission through increased sporozoite prevalence in kdr resistant mosquitoes. These results are of great concern for the epidemiology of malaria considering the widespread pyrethroid resistance currently observed in Sub-Saharan Africa and the efforts deployed to control the disease.

Chromosome inversions, genomic differentiation and speciation in the African malaria mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.

Modelling the risk of being bitten by malaria vectors in a vector control area in southern Benin, west Africa

Background

The diversity of malaria vector populations, expressing various resistance and/or behavioural patterns could explain the reduced effectiveness of vector control interventions reported in some African countries. A better understanding of the ecology and distribution of malaria vectors is essential to design more effective and sustainable strategies for malaria control and elimination. Here, we analyzed the spatio-temporal risk of the contact between humans and the sympatric An. funestus and both M and S molecular forms of An. gambiae s.s. in an area of Benin with high coverage of vector control measures with an unprecedented level of resolution.

Methods

Presence-absence data for the three vectors from 1-year human-landing collections in 19 villages were assessed using binomial mixed-effects models according to vector control measures and environmental covariates derived from field and remote sensing data. After 8-fold cross-validations of the models, predictive maps of the risk of the contact between humans and the sympatric An. funestus and both molecular M and S forms of An. gambiae s.s. were computed.

Results

Model validations showed that the An. funestus, An. gambiae M form, and S form models provided an excellent (Area Under Curve>0.9), a good (AUC>0.8), and an acceptable (AUC>0.7) level of prediction, respectively. The distribution area of the probability of contact between human and An. funestus largely overlaps that of An. gambiae M form but this latter showed important seasonal variation. An. gambiae S form also showed seasonal variation but with different ecological preferences. Landscape data were useful to discriminate between the species’ distributions.

Conclusions

These results showed that available remote sensing data could help in predicting the human-vector contact for several species of malaria vectors at a village level scale. The predictive maps showed seasonal and spatial variations in the risk of human-vector contact for all three vectors. Such maps could help Malaria Control Programmes to implement more effective vector control strategy by taking into account to the dynamics of malaria vector species.