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.

Speciation requires the acquisition of reproductive isolation, and the circumstances under which this could evolve are of great interest. Are new species formed after the acquisition of generalized incompatibility arising between physically separated populations, or may they arise as a result of the action of disruptive selection beginning with the divergence of a rather restricted set of gene loci? Here we apply the technique of amplified fragment length polymorphism (AFLP) analysis to an intertidal snail whose populations display a cline in shell shape across vertical gradients on rocky shores. We compare the FST values for 306 AFLP loci with the distribution of FST estimated from a simulation model using values of mutation and migration derived from the data. We find that about 5% of these loci show greater differentiation than expected, providing evidence of the effects of selection across the cline, either direct or indirect through linkage. This is consistent with expectations from nonallopatric speciation models that propose an initial divergence of a small part of the genome driven by strong disruptive selection while divergence at other loci is prevented by gene flow. However, the pattern could also be the result of differential introgression after secondary contact.

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.

In the last decade there have been marked reductions in malaria incidence in sub-Saharan Africa. Sustaining these reductions will rely upon insecticides to control the mosquito malaria vectors. We report that in the primary African malaria vector, Anopheles gambiae sensu stricto, a single enzyme, CYP6M2, confers resistance to two classes of insecticide. This is unique evidence in a disease vector of cross-resistance associated with a single metabolic gene that simultaneously reduces the efficacy of two of the four classes of insecticide routinely used for malaria control. The gene-expression profile of a highly DDT-resistant population of A. gambiae s.s. from Ghana was characterized using a unique whole-genome microarray. A number of genes were significantly overexpressed compared with two susceptible West African colonies, including genes from metabolic families previously linked to insecticide resistance. One of the most significantly overexpressed probe groups (false-discovery rate-adjusted P < 0.0001) belonged to the cytochrome P450 gene CYP6M2. This gene is associated with pyrethroid resistance in wild A. gambiae s.s. populations) and can metabolize both type I and type II pyrethroids in recombinant protein assays. Using in vitro assays we show that recombinant CYP6M2 is also capable of metabolizing the organochlorine insecticide DDT in the presence of solubilizing factor sodium cholate.

Insecticide resistance is an ideal model to study the emergence and spread of adaptative variants. In the African malaria mosquito, Anopheles gambiae, this is complemented by a strong public health rationale. In this insect, resistance to pyrethroid and DDT insecticides is strongly associated with the mutations L1014F and L1014S within the para voltage-gated sodium channel (VGSC). Across much of West Africa, 1014F frequency approaches fixation. Here, we document the emergence of a mutation, N1575Y, within the linker between domains III-IV of the VGSC. In data extending over 40 kbp of the VGSC 1575Y occurs on only a single long-range haplotype, also bearing 1014F. The 1014F-1575Y haplotype was found in both M and S molecular forms of An. gambiae in West/ Central African sample sites separated by up to 2,000 km. In Burkina Faso M form, 1575Y allele frequency rose significantly from 0.053 to 0.172 between 2008 and 2010. Extended haplotype homozygosity analysis of the wild-type 1575N allele showed rapid decay of linkage disequilibrium (LD), in sharp contrast to the extended LD exhibited by 1575Y. A haplotype with long-range LD and high/increasing frequency is a classical sign of strong positive selection acting on a recent mutant. 1575Y occurs ubiquitously on a 1014F haplotypic background, suggesting that the N1575Y mutation compensates for deleterious fitness effects of 1014F and/or confers additional resistance to insecticides. Haplotypic tests of association suggest the latter: The 1014F-1575Y haplotype confers a significant additive benefit above 1014F-1575N for survival to DDT (M form P = 0.03) and permethrin (S form P = 0.003).kdr | selective sweep | inactivation particle