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

The demographic histories of the M and S molecular forms of Anopheles gambiae s.s

Anopheles gambiae is a primary vector of Plasmodium falciparum, a human malaria parasite that causes over a million deaths each year in sub-Saharan Africa. Population genetic tests have been employed to detect natural selection at suspected A. gambiae antimalaria genes, but these tests have generally been compromised by the lack of demographically correct null models. Here, we used a coalescent simulation approach within a maximum likelihood framework to fit population growth, bottleneck, and migration models to polymorphism data from Cameroonian A. gambiae. The best-fit models for both the "M" and the "S" molecular forms of A. gambiae included ancient population growth and a high rate of migration from an unsampled subpopulation. After correcting for differences in effective population size, our models suggest that the molecular forms expanded at different times and both expansions significantly predate the advent of agriculture. We show that correcting null models for demography increases the power to detect natural selection in A. gambiae.

Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae

Previous efforts to uncover the genetic underpinnings of ongoing ecological speciation of the M and S forms of the African malaria vector Anopheles gambiae revealed two centromere-proximal islands of genetic divergence on X and chromosome 2. Under the assumption of considerable ongoing gene flow between M and S, these persistently divergent genomic islands were widely considered to be 'speciation islands'. In the course of microarray-based divergence mapping, we discovered a third centromere-associated island of divergence on chromosome 3, which was validated by targeted re-sequencing. To test for genetic association between the divergence islands on all three chromosomes, SNP-based assays were applied in four natural populations of M and S spanning West, Central and East Africa. Genotyping of 517 female M and S mosquitoes revealed nearly complete linkage disequilibrium between the centromeres of the three independently assorting chromosomes. These results suggest that despite the potential for inter-form gene flow through hybridization, actual (realized) gene flow between M and S may be substantially less than commonly assumed and may not explain most shared variation. Moreover, the possibility of very low gene flow calls into question whether diverged pericentromeric regions-characterized by reduced levels of variation and recombination-are in fact instrumental rather than merely incidental to the speciation process.

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

Disruptive selection mediated by predation on aquatic immature stages has been proposed as a major force driving ecological divergence and fostering speciation between the M and S molecular forms of the African malaria mosquito, Anopheles gambiae. In the dry savannahs of West Africa where both molecular forms co-occur, the S form thrives in temporary pools filled with rainwater, whereas the M form preferentially breeds in permanent freshwater habitats where predator pressure is higher. Here, we explored the proximal mechanisms by which predation may contribute to habitat segregation between molecular forms using progeny of female mosquitoes captured in Burkina Faso. We show that the S form suffers higher predation rates than the M form when simultaneously exposed to the widespread predator, Anisops jaczewskii in an experimental arena. Furthermore, behavioral plasticity induced by exposure to the predator was observed in the M form, but not in the S form, and may partially explain its habitat use and ecological divergence from the S form. We discuss the role of adaptive phenotypic plasticity in allowing successful colonization of a new ecological niche by the M form and highlight further research areas that need to be addressed for a better understanding of the ultimate mechanisms underlying ecological speciation in this pest of major medical importance.

"Singing on the wing" as a mechanism for species recognition in the malarial mosquito Anopheles gambiae

Anopheles gambiae, responsible for the majority of malaria deaths annually, is a complex of seven species and several chromosomal/molecular forms. The complexity of malaria epidemiology and control is due in part to An. gambiae's remarkable genetic plasticity, enabling its adaptation to a range of human-influenced habitats. This leads to rapid ecological speciation when reproductive isolation mechanisms develop [1-6]. Although reproductive isolation is essential for speciation, little is known about how it occurs in sympatric populations of incipient species [2]. We show that in such a population of "M" and "S" molecular forms, a novel mechanism of sexual recognition (male-female flight-tone matching [7-9]) also confers the capability of mate recognition, an essential precursor to assortative mating; frequency matching occurs more consistently in same-form pairs than in mixed-form pairs (p = 0.001). [corrected] Furthermore, the key to frequency matching is "difference tones" produced in the nonlinear vibrations of the antenna by the combined flight tones of a pair of mosquitoes and detected by the Johnston's organ. By altering their wing-beat frequencies to minimize these difference tones, mosquitoes can match flight-tone harmonic frequencies above their auditory range. This is the first description of close-range mating interactions in incipient An. gambiae species.

Habitat suitability and ecological niche profile of major malaria vectors in Cameroon

BACKGROUND

Suitability of environmental conditions determines a species distribution in space and time. Understanding and modelling the ecological niche of mosquito disease vectors can, therefore, be a powerful predictor of the risk of exposure to the pathogens they transmit. In Africa, five anophelines are responsible for over 95% of total malaria transmission. However, detailed knowledge of the geographic distribution and ecological requirements of these species is to date still inadequate.

METHODS

Indoor-resting mosquitoes were sampled from 386 villages covering the full range of ecological settings available in Cameroon, Central Africa. Using a predictive species distribution modeling approach based only on presence records, habitat suitability maps were constructed for the five major malaria vectors Anopheles gambiae, Anopheles funestus, Anopheles arabiensis, Anopheles nili and Anopheles moucheti. The influence of 17 climatic, topographic, and land use variables on mosquito geographic distribution was assessed by multivariate regression and ordination techniques.

RESULTS

Twenty-four anopheline species were collected, of which 17 are known to transmit malaria in Africa. Ecological Niche Factor Analysis, Habitat Suitability modeling and Canonical Correspondence Analysis revealed marked differences among the five major malaria vector species, both in terms of ecological requirements and niche breadth. Eco-geographical variables (EGVs) related to human activity had the highest impact on habitat suitability for the five major malaria vectors, with areas of low population density being of marginal or unsuitable habitat quality. Sunlight exposure, rainfall, evapo-transpiration, relative humidity, and wind speed were among the most discriminative EGVs separating "forest" from "savanna" species.

CONCLUSIONS

The distribution of major malaria vectors in Cameroon is strongly affected by the impact of humans on the environment, with variables related to proximity to human settings being among the best predictors of habitat suitability. The ecologically more tolerant species An. gambiae and An. funestus were recorded in a wide range of eco-climatic settings. The other three major vectors, An. arabiensis, An. moucheti, and An. nili, were more specialized. Ecological niche and species distribution modelling should help improve malaria vector control interventions by targeting places and times where the impact on vector populations and disease transmission can be optimized.

Seasonal distribution of Anopheles funestus chromosomal forms from Burkina Faso

Background

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

Methods

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

Results

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

Conclusion

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

The population genomics of trans-specific inversion polymorphisms in Anopheles gambiae

In the malaria mosquito Anopheles gambiae polymorphic chromosomal inversions may play an important role in adaptation to environmental variation. Recently, we used microarray-based divergence mapping combined with targeted resequencing to map nucleotide differentiation between alternative arrangements of the 2La inversion. Here, we applied the same technique to four different polymorphic inversions on the 2R chromosome of An. gambiae. Surprisingly, divergence was much lower between alternative arrangements for all 2R inversions when compared to the 2La inversion. For one of the rearrangements, 2Ru, we successfully mapped a very small region (approximately 100 kb) of elevated divergence. For the other three rearrangements, we did not identify any regions of significantly high divergence, despite ample independent evidence from natural populations of geographic clines and seasonal cycling, and stable heterotic polymorphisms in laboratory populations. If these inversions are the targets of selection as hypothesized, we suggest that divergence between rearrangements may have escaped detection due to retained ancestral polymorphism in the case of the youngest 2R rearrangements and to extensive gene flux in the older 2R inversion systems that segregate in both An. gambiae and its sibling species An. arabiensis.

Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation

BACKGROUND

Speciation 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.

RESULTS

Species 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.

CONCLUSION

Our 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.