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

The sustainability of malaria control in Africa is threatened by the rise of insecticide resistance in Anopheles mosquitoes that transmit the disease1. To gain a deeper understanding of how mosquito populations are evolving, we sequenced the genomes of 765 specimens of Anopheles gambiae and Anopheles coluzzii sampled from 15 locations across Africa, identifying over 50 million single nucleotide polymorphisms within the accessible genome. These data revealed complex population structure and patterns of gene flow, with evidence of ancient expansions, recent bottlenecks, and local variation in effective population size. Strong signals of recent selection were observed in insecticide resistance genes, with multiple sweeps spreading over large geographical distances and between species. The design of novel tools for mosquito control using gene drive will need to take account of high levels of genetic diversity in natural mosquito populations.

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.

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.