A behavioral mechanism underlying ecological divergence in the malaria mosquito Anopheles gambiae

BackgroundLimitations in the ability of organisms to tolerate environmental stressors affect their fundamental ecological niche and constrain their distribution to specific habitats. Evolution of tolerance, therefore, can engender ecological niche dynamics. Forest populations of the afro-tropical malaria mosquito Anopheles gambiae have been shown to adapt to historically unsuitable larval habitats polluted with decaying organic matter that are found in densely populated urban agglomerates of Cameroon. This process has resulted in niche expansion from rural to urban environments that is associated with cryptic speciation and ecological divergence of two evolutionarily significant units within this taxon, the molecular forms M and S, among which reproductive isolation is significant but still incomplete. Habitat segregation between the two forms results in a mosaic distribution of clinally parapatric patches, with the M form predominating in the centre of urban agglomerates and the S form in the surrounding rural localities. We hypothesized that development of tolerance to nitrogenous pollutants derived from the decomposition of organic matter, among which ammonia is the most toxic to aquatic organisms, may affect this pattern of distribution and process of niche expansion by the M form.ResultsAcute toxicity bioassays indicated that populations of the two molecular forms occurring at the extremes of an urbanization gradient in Yaounde, the capital of Cameroon, differed in their response to ammonia. The regression lines best describing the dose-mortality profile differed in the scale of the explanatory variable (ammonia concentration log-transformed for the S form and linear for the M form), and in slope (steeper for the S form and shallower for the M form). These features reflected differences in the frequency distribution of individual tolerance thresholds in the two populations as assessed by probit analysis, with the M form exhibiting a greater mean and variance compared to the S form.ConclusionsIn agreement with expectations based on the pattern of habitat partitioning and exposure to ammonia in larval habitats in Yaounde, the M form showed greater tolerance to ammonia compared to the S form. This trait may be part of the physiological machinery allowing forest populations of the M form to colonize polluted larval habitats, which is at the heart of its niche expansion in densely populated human settlements in Cameroon.

Summary Predation is a major evolutionary force driving speciation. Identifying the stimuli prompting anti‐predator responses is essential for unravelling the proximate mechanisms of anti‐predator adaptations and for understanding how predation impacts species diversification. Here, we explore for the first time the divergence in the use of cues of predation risk by different incipient species of mosquitoes within the Anopheles gambiae complex. We examined the anti‐predator responses to various predation cues by the aquatic larvae of Anopheles arabiensis and the two molecular forms, M and S, of Anopheles gambiae s.s. sampled from wild populations in Burkina Faso naturally exposed to different predator densities. The larvae altered their behaviour in response to various predation cues, shifting from the surface of the water to a less exposed location on the walls of the experimental arena. There were important differences in the use of predation cues between populations sharing the same larval development sites and, within the M form, related to predation pressure at their site of origin. Anopheles arabiensis larvae, which develop in small temporary waterbodies, relied only on physical cues signalling the presence of a predator to mount a vigilance response, whereas the sympatric larvae of the M and S forms of An. gambiae used both physical and chemical cues to fine‐tune their response. M‐form larvae, developing in permanent aquatic habitats with high densities of predators, responded chiefly to chemical stimuli ensuing from acute predation. Our results suggest that predation might play a role as a disruptive selective force promoting ecological divergence within the An. gambiae complex and, more recently, between the M and S molecular forms. Notably, beyond interspecific differences in the use of predation cues, our study also documented phenotypic differences within the M form based on predation pressure at their site of origin, suggesting that a process of ecological divergence is generating species ‘forerunners’ within the An. gambiae complex in the wild.

The molecular forms of Anopheles gambiae: A phenotypic perspective

BackgroundSpeciation among members of the Anopheles gambiae complex is thought to be promoted by disruptive selection and ecological divergence acting on sets of adaptation genes protected from recombination by polymorphic paracentric chromosomal inversions. However, shared chromosomal polymorphisms between the M and S molecular forms of An. gambiae and insufficient information about their relationship with ecological divergence challenge this view. We used Geographic Information Systems, Ecological Niche Factor Analysis, and Bayesian multilocus genetic clustering to explore the nature and extent of ecological and chromosomal differentiation of M and S across all the biogeographic domains of Cameroon in Central Africa, in order to understand the role of chromosomal arrangements in ecological specialisation within and among molecular forms.ResultsSpecies distribution modelling with presence-only data revealed differences in the ecological niche of both molecular forms and the sibling species, An. arabiensis. The fundamental environmental envelope of the two molecular forms, however, overlapped to a large extent in the rainforest, where they occurred in sympatry. The S form had the greatest niche breadth of all three taxa, whereas An. arabiensis and the M form had the smallest niche overlap. Correspondence analysis of M and S karyotypes confirmed that molecular forms shared similar combinations of chromosomal inversion arrangements in response to the eco-climatic gradient defining the main biogeographic domains occurring across Cameroon. Savanna karyotypes of M and S, however, segregated along the smaller-scale environmental gradient defined by the second ordination axis. Population structure analysis identified three chromosomal clusters, each containing a mixture of M and S specimens. In both M and S, alternative karyotypes were segregating in contrasted environments, in agreement with a strong ecological adaptive value of chromosomal inversions.ConclusionOur data suggest that inversions on the second chromosome of An. gambiae are not causal to the evolution of reproductive isolation between the M and S forms. Rather, they are involved in ecological specialization to a similar extent in both genetic backgrounds, and most probably predated lineage splitting between molecular forms. However, because chromosome-2 inversions promote ecological divergence, resulting in spatial and/or temporal isolation between ecotypes, they might favour mutations in other ecologically significant genes to accumulate in unlinked chromosomal regions. When such mutations occur in portions of the genome where recombination is suppressed, such as the pericentromeric regions known as speciation islands in An. gambiae, they would contribute further to the development of reproductive isolation.

Vector-borne infectious diseases are responsible for the deaths of over 700,000 people annually, than 400,000 of them resulting from malaria. The mosquito Anopheles gambiae is one of the dominant vector species of human malaria transmission. A significant issue of the conventional insecticides which target the arthropod borne infectious diseases is their induced resistance. To overcome this inconvenience, insecticides with new modes of action are required. One of the most promising targets for the development of new potential insecticides as evidenced by current studies is the D1-like dopamine receptor (DAR). To get a deeper understanding of the structural information of this receptor, the 3D homology model was built. The possible sites within the protein were identified and the most probable binding site was highlighted. The homology model along with a series of DAR antagonists with known activity against Anopheles gambiae larvae were used in docking experiments to gain insight into their intermolecular interactions. Furthermore, virtual screening of the natural compounds from the SPECS database led to the prediction of toxicity and environmental hazards for one potential new insecticide against the Anopheles gambiae mosquito.

BackgroundGenes involved in post-mating processes of multiple mating organisms are known to evolve rapidly due to coevolution driven by sexual conflict among male-female interacting proteins. In the malaria mosquito Anopheles gambiae - a monandrous species in which sexual conflict is expected to be absent or minimal - recent data strongly suggest that proteolytic enzymes specifically expressed in the female lower reproductive tissues are involved in the processing of male products transferred to females during mating. In order to better understand the role of selective forces underlying the evolution of proteins involved in post-mating responses, we analysed a cluster of genes encoding for three serine proteases that are down-regulated after mating, two of which specifically expressed in the atrium and one in the spermatheca of A. gambiae females.ResultsThe analysis of polymorphisms and divergence of these female-expressed proteases in closely related species of the A. gambiae complex revealed a high level of replacement polymorphisms consistent with relaxed evolutionary constraints of duplicated genes, allowing to rapidly fix novel replacements to perform new or more specific functions. Adaptive evolution was detected in several codons of the 3 genes and hints of episodic selection were also found. In addition, the structural modelling of these proteases highlighted some important differences in their substrate specificity, and provided evidence that a number of sites evolving under selective pressures lie relatively close to the catalytic triad and/or on the edge of the specificity pocket, known to be involved in substrate recognition or binding. The observed patterns suggest that these proteases may interact with factors transferred by males during mating (e.g. substrates, inhibitors or pathogens) and that they may have differently evolved in independent A. gambiae lineages.ConclusionsOur results - also examined in light of constraints in the application of selection-inference methods to the closely related species of the A. gambiae complex - reveal an unexpectedly intricate evolutionary scenario. Further experimental analyses are needed to investigate the biological functions of these genes in order to better interpret their molecular evolution and to assess whether they represent possible targets for limiting the fertility of Anopheles mosquitoes in malaria vector control strategies.

BackgroundPredation of aquatic immature stages has been identified as a major evolutionary force driving habitat segregation and niche partitioning in the malaria mosquito Anopheles gambiae sensu stricto in the humid savannahs of Burkina Faso, West Africa. Here, we explored behavioural responses to the presence of a predator in wild populations of the M and S molecular forms of An. gambiae that typically breed in permanent (e.g., rice field paddies) and temporary (e.g., road ruts) water collections.MethodsLarvae used in these experiments were obtained from eggs laid by wild female An. gambiae collected from two localities in south-western Burkina Faso during the 2008 rainy season. Single larvae were observed in an experimental arena, and behavioural traits were recorded and quantified a) in the absence of a predator and b) in the presence of a widespread mosquito predator, the backswimmer Anisops jaczewskii. Differences in the proportion of time allocated to each behaviour were assessed using Principal Component Analysis and Multivariate Analysis of Variance.ResultsThe behaviour of M and S form larvae was found to differ significantly; although both forms mainly foraged at the water surface, spending 60-90% of their time filtering water at the surface or along the wall of the container, M form larvae spent on average significantly more time browsing at the bottom of the container than S form larvae (4.5 vs. 1.3% of their overall time, respectively; P < 0.05). In the presence of a predator, larvae of both forms modified their behaviour, spending significantly more time resting along the container wall (P < 0.001). This change in behaviour was at least twice as great in the M form (from 38.6 to 66.6% of the time at the wall in the absence and presence of the predator, respectively) than in the S form (from 48.3 to 64.1%). Thrashing at the water surface exposed larvae to a significantly greater risk of predation by the notonectid (P < 0.01), whereas predation occurred significantly less often when larvae were at the container wall (P < 0.05) and might reflect predator vigilance.ConclusionsBehavioural differences between larvae of the M and S form of An. gambiae in response to an acute predation risk is likely to be a reflection of different trade-offs between foraging and predator vigilance that might be of adaptive value in contrasting aquatic ecosystems. Future studies should explore the relevance of these findings under the wide range of natural settings where both forms co-exist in Africa.

BackgroundGermline specific promoters are an essential component of potential vector control strategies which function by genetic drive, however suitable promoters are not currently available for the main human malaria vector Anopheles gambiae.ResultsWe have identified the Anopheles gambiae vasa-like gene and found its expression to be specifically localized to both the male and female gonads in adult mosquitoes. We have functionally characterised using transgenic reporter lines the regulatory regions required for driving transgene expression in a pattern mirroring that of the endogenous vasa locus. Two reporter constructs indicate the existence of distinct vasa regulatory elements within the 5' untranslated regions responsible not only for the spatial and temporal but also for the sex specific germline expression. vasa driven eGFP expression in the ovary of heterozygous mosquitoes resulted in the progressive accumulation of maternal protein and transcript in developing oocytes that were then detectable in all embryos and neonatal larvae.ConclusionWe have characterized the vasa regulatory regions that are not only suited to drive transgenes in the early germline of both sexes but could also be utilized to manipulate the zygotic genome of developing embryos via maternal deposition of active molecules. We have used computational models to show that a homing endonuclease-based gene drive system can function in the presence of maternal deposition and describe a novel non-invasive control strategy based on early vasa driven homing endonuclease expression.

Mosquitoes genetically engineered to be resistant to Plasmodium parasites represent a promising novel approach in the fight against malaria. The insect immune system itself is a source of anti-parasitic genes potentially exploitable for transgenic designs. The Anopheles gambiae thioester containing protein 1 (TEP1) is a potent anti-parasitic protein. TEP1 is secreted and circulates in the mosquito hemolymph, where its activated cleaved form binds and eliminates malaria parasites. Here we investigated whether TEP1 can be used to create malaria resistant mosquitoes. Using a GFP reporter transgene, we determined that the fat body is the main site of TEP1 expression. We generated transgenic mosquitoes that express TEP1r, a potent refractory allele of TEP1, in the fat body and examined the activity of the transgenic protein in wild-type or TEP1 mutant genetic backgrounds. Transgenic TEP1r rescued loss-of-function mutations, but did not increase parasite resistance in the presence of a wild-type susceptible allele. Consistent with previous reports, TEP1 protein expressed from the transgene in the fat body was taken up by hemocytes upon a challenge with injected bacteria. Furthermore, although maturation of transgenic TEP1 into the cleaved form was impaired in one of the TEP1 mutant lines, it was still sufficient to reduce parasite numbers and induce parasite melanization. We also report here the first use of Transcription Activator Like Effectors (TALEs) in Anopheles gambiae to stimulate expression of endogenous TEP1. We found that artificial elevation of TEP1 expression remains moderate in vivo and that enhancement of endogenous TEP1 expression did not result in increased resistance to Plasmodium. Taken together, our results reveal the difficulty of artificially influencing TEP1-mediated Plasmodium resistance, and contribute to further our understanding of the molecular mechanisms underlying mosquito resistance to Plasmodium parasites.

BackgroundOngoing lineage splitting within the African malaria mosquito Anopheles gambiae is compatible with ecological speciation, the evolution of reproductive isolation by divergent natural selection acting on two populations exploiting alternative resources. Divergence between two molecular forms (M and S) identified by fixed differences in rDNA, and characterized by marked, although incomplete, reproductive isolation is occurring in West and Central Africa. To elucidate the role that ecology and geography play in speciation, we carried out a countrywide analysis of An. gambiae M and S habitat requirements, and that of their chromosomal variants, across Burkina Faso.ResultsMaps of relative abundance by geostatistical interpolators produced a distinct pattern of distribution: the M-form dominated in the northernmost arid zones, the S-form in the more humid southern regions. Maps of habitat suitability, quantified by Ecological Niche Factor Analysis based on 15 eco-geographical variables revealed less contrast among forms. M was peculiar as it occurred proportionally more in habitat of marginal quality. Measures of ecological niche breadth and overlap confirmed the mismatch between the fundamental and realized patterns of habitat occupation: forms segregated more than expected from the extent of divergence of their environmental envelope – a signature of niche expansion. Classification of chromosomal arm 2R karyotypes by multilocus genetic clustering identified two clusters loosely corresponding to molecular forms, with 'mismatches' representing admixed individuals due to shared ancestral polymorphism and/or residual hybridization. In multivariate ordination space, these karyotypes plotted in habitat of more marginal quality compared to non-admixed, 'typical', karyotypes. The distribution of 'typical' karyotypes along the main eco-climatic gradient followed a consistent pattern within and between forms, indicating an adaptive role of inversions at this geographical scale.ConclusionEcological segregation between M and S is consistent with niche expansion into marginal habitats by chromosomal inversion variants during early lineage divergence; presumably, this process is promoted by inter-karyotype competition in the higher-quality core habitat. We propose that the appearance of favourable allelic combinations in other regions of suppressed recombination (e.g. pericentromeric portions defining speciation islands in An. gambiae) fosters development of reproductive isolation to protect linkage between separate chromosomal regions.

In West Africa, lineage splitting between the M and S molecular forms of the major Afro-tropical malaria mosquito, Anopheles gambiae (Diptera: Culicidae), is thought to be driven by ecological divergence, occurring mainly at the larval stage. Here, we present evidence for habitat segregation between the two molecular forms in and around irrigated rice fields located within the humid savannahs of western Burkina Faso. Longitudinal sampling of adult mosquitoes emerging from a range of breeding sites distributed along a transect extending from the heart of the rice field area into the surrounding savannah was conducted from June to November 2009. Analysis revealed that the two molecular forms and their sibling species Anopheles arabiensis are not randomly distributed in the area. A major ecological gradient was extracted in relation to the perimeter of the rice fields. The M form was associated with larger breeding sites mostly consisting of rice paddies, whereas the S form and An. arabiensis were found to depend upon temporary, rain-filled breeding sites. These results support hypotheses about larval habitat segregation and confirm the suggestion that the forms have different larval habitat requirements. Segregation appears to be clearly linked to anthropogenic permanent habitats and the community structure they support.

Molecular characterization of insulin-like peptides in the yellow fever mosquito, Aedes aegypti: Expression, cellular localization, and phylogeny

Anopheles gambiae, the major malaria vector in Africa, can be divided into two subgroups based on genetic and ecological criteria. These two subgroups, termed the M and S molecular forms, are believed to be incipient species. Although they display differences in the ecological niches they occupy in the field, they are often sympatric and readily hybridize in the laboratory to produce viable and fertile offspring. Evidence for assortative mating in the field was recently reported, but the underlying mechanisms awaited discovery. We studied swarming behaviour of the molecular forms and investigated the role of swarm segregation in mediating assortative mating. Molecular identification of 1145 males collected from 68 swarms in Donéguébougou, Mali, over 2 years revealed a strict pattern of spatial segregation, resulting in almost exclusively monotypic swarms with respect to molecular form. We found evidence of clustering of swarms composed of individuals of a single molecular form within the village. Tethered M and S females were introduced into natural swarms of the M form to verify the existence of possible mate recognition operating within-swarm. Both M and S females were inseminated regardless of their form under these conditions, suggesting no within-mate recognition. We argue that our results provide evidence that swarm spatial segregation strongly contributes to reproductive isolation between the molecular forms in Mali. However this does not exclude the possibility of additional mate recognition operating across the range distribution of the forms. We discuss the importance of spatial segregation in the context of possible geographic variation in mechanisms of reproductive isolation.

We report the first finding of the knockdown Leu-Phe and Leu-Ser mutations associated with resistance to pyrethroids and DDT insecticides in the malaria mosquito Anopheles gambiae from Cameroon. The Leu-Phe mutation was found in both the M and S molecular forms of An. gambiae. Importantly, two specimens of the S molecular form were found to carry both mutations in a heterozygous state.

Chromosomal inversions are thought to confer a selective advantage in alternative habitats by protecting co-adapted alleles from recombination. The frequencies of two inversions (2La and 2Rb) of the afro-tropical malaria mosquito Anopheles gambiae change gradually along geographical clines, increasing in frequency with degree of aridity. Such clines can result from gene flow and local selection acting upon alternative karyotypes along the cline, suggesting that these inversions may be associated with tolerance to xeric conditions. Since water loss represents a major challenge in xeric habitats, it can be supposed that genes inside these inversions are involved in water homeostasis. To test this hypothesis, we compared the desiccation resistance of alternative karyotypes from a colonised 2Rb/2La polymorphic population of A. gambiae from Cameroon. The strain included only the molecular form S, one of the genetic units marking incipient speciation in this taxon. Day-old mosquitoes of both sexes were assayed individually for time to death in a dry environment and the karyotype of each was determined post-mortem using molecular diagnostic assays for each inversion. In agreement with expectations based on their eco-geographical distribution, we found that 2La homokaryotypes survived significantly longer (1.3 hours) than the other karyotypes. However, there was weak support for the effect of 2Rb on desiccation resistance. Larger mosquitoes survived longer than smaller ones. Median survival of females was greater than males, but the effect of sex on desiccation resistance was weakly supported, indicating that differential survival was correlated to differences between sexes in average size. We found weak evidence for a heterotic effect of 2La karyotype on size in females. These results support the notion that genes located inside the 2La inversion are involved in water balance, contributing towards local adaptation of A. gambiae to xeric habitats, beyond the adaptive value conferred by a larger body size.