Michel A. Slotman

INTRODUCTION Control of mosquito vectors has historically proven to be an effective means of eliminating malaria. Human malaria is transmitted only by mosquitoes in the genus Anopheles , but not all species within the genus, or even all members of each vector species, are efficient malaria vectors. Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. RATIONALE This variation in vectorial capacity suggests an underlying genetic/genomic plasticity that results in variation of key traits determining vectorial capacity within the genus. Sequencing the genome of Anopheles gambiae , the most important malaria vector in sub-Saharan Africa, has offered numerous insights into how that species became highly specialized to live among and feed upon humans and how susceptibility to mosquito control strategies is determined. Until very recently, similar genomic resources have not existed for other anophelines, limiting comparisons to individual genes or sets of genomic markers with no genome-wide data to investigate attributes associated with vectorial capacity across the genus. RESULTS We sequenced and assembled the genomes and transcriptomes of 16 anophelines from Africa, Asia, Europe, and Latin America, spanning ~100 million years of evolution and chosen to represent a range of evolutionary distances from An. gambiae , a variety of geographic locations and ecological conditions, and varying degrees of vectorial capacity. Genome assembly quality reflected DNA template quality and homozygosity. Despite variation in contiguity, the assemblies were remarkably complete and searches for arthropod-wide single-copy orthologs generally revealed few missing genes. Genome annotation supported with RNA sequencing transcriptomes yielded between 10,738 and 16,149 protein-coding genes for each species. Relative to Drosophila, the closest dipteran genus for which equivalent genomic resources exist, Anopheles exhibits a dynamic genomic evolutionary profile. Comparative analyses show a fivefold faster rate of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses in Anopheles . Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. We also document evidence of variation in important reproductive phenotypes, genes controlling immunity to Plasmodium malaria parasites and other microbes, genes encoding cuticular and salivary proteins, and genes conferring metabolic insecticide resistance. 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. CONCLUSIONS Anopheline mosquitoes exhibit a molecular evolutionary profile very distinct from Drosophila , and their genomes harbor strong evidence of functional variation in traits that determine vectorial capacity. These 16 new reference genome assemblies provide a foundation for hypothesis generation and testing to further our understanding of the diverse biological traits that determine vectorial capacity. Geography, vector status, and molecular phylogeny of the 16 newly sequenced anopheline mosquitoes and selected other dipterans. The maximum likelihood molecular phylogeny of all sequenced anophelines and two mosquito outgroups was constructed from the aligned protein sequences of 1085 single-copy orthologs. Shapes between branch termini and species names indicate vector status and are colored according to geographic ranges depicted on the map.

Understanding the processes by which species colonize and adapt to human habitats is particularly important in the case of disease-vectoring arthropods. The mosquito species Aedes aegypti, a major vector of dengue and yellow fever viruses, probably originated as a wild, zoophilic species in sub-Saharan Africa, where some populations still breed in tree holes in forested habitats. Many populations of the species, however, have evolved to thrive in human habitats and to bite humans. This includes some populations within Africa as well as almost all those outside Africa. It is not clear whether all domestic populations are genetically related and represent a single 'domestication' event, or whether association with human habitats has developed multiple times independently within the species. To test the hypotheses above, we screened 24 worldwide population samples of Ae. aegypti at 12 polymorphic microsatellite loci. We identified two distinct genetic clusters: one included all domestic populations outside of Africa and the other included both domestic and forest populations within Africa. This suggests that human association in Africa occurred independently from that in domestic populations across the rest of the world. Additionally, measures of genetic diversity support Ae. aegypti in Africa as the ancestral form of the species. Individuals from domestic populations outside Africa can reliably be assigned back to their population of origin, which will help determine the origins of new introductions of Ae. aegypti.

BACKGROUND: Vector control through indoor residual spraying (IRS) has been employed on Bioko Island, Equatorial Guinea, under the Bioko Island Malaria Control Project (BIMCP) since 2004. This study analyses the change in mosquito abundance, species composition and outdoor host-seeking proportions from 2009 to 2014, after 11 years of vector control on Bioko Island. METHODS: All-night indoor and outdoor human landing catches were performed monthly in the Bioko Island villages of Mongola, Arena Blanca, Biabia and Balboa from 2009 to 2014. Collected mosquitoes were morphologically identified and a subset of Anopheles gambiae sensu lato (s.l.) were later identified molecularly to their sibling species. Mosquito collection rates, species composition and indoor/outdoor host-seeking sites were analysed using generalized linear mixed models to assess changes in mosquito abundance and behaviour. RESULTS: The overall mosquito collection rate declined in each of the four Bioko Island villages. Anopheles coluzzii and Anopheles melas comprised the An. gambiae s.l. mosquito vector population, with a range of species proportions across the four villages. The proportion of outdoor host-seeking An. gambiae s.l. mosquitoes increased significantly in all four villages with an average increase of 58.8 % [57.9, 59.64 %] in 2009 to 70.0 % [67.8, 72.0 %] in 2014. Outdoor host-seeking rates did not increase in the month after an IRS spray round compared to the month before, suggesting that insecticide repellency has little impact on host-seeking behaviour. CONCLUSION: While vector control on Bioko Island has succeeded in substantial reduction in overall vector biting rates, populations of An. coluzzii and An. melas persist. Host-seeking behaviour has changed in these An. gambiae s.l. populations, with a shift towards outdoor host-seeking. During this study period, the proportion of host-seeking An. gambiae s.l. caught outdoors observed on Bioko Island increased to high levels, exceeding 80 % in some locations. It is possible that there may be a genetic basis underlying this large shift in host-seeking behaviour, in which case outdoor feeding could pose a serious threat to current vector control programmes. Currently, the BIMCP is preparing for this potential challenge by testing source reduction as a complementary control effort that also targets outdoor transmission.

Evidence for introgression between Anopheles gambiae and An. arabiensis has accumulated for some time. We examined the fate of introgressed DNA directly, using microsatellite markers located throughout the genome. Introgressed X chromosomes were removed within two generations. Furthermore, substantial differences in introgressive capacity between the two autosomes were found. After introgression from An. arabiensis into An. gambiae, most introgressed alleles at third chromosome markers, particularly those on 3R, decreased steadily, indicating selection against them. No such pattern was observed for 2L markers and several 2R markers. The frequency of introgressed alleles on 2L were close to the original frequency even after 19 generations, whereas only two 2R markers showed a modest decrease. Even though limited information was available on the reciprocal cross, the pattern appears to be identical. Although the decrease in frequency of the introgressed X chromosome can be attributed to the presence of sterility and inviability effects, the variation in introgressive capacity of the autosomes does not appear to be explained by the presence of inversion polymorphisms, or regions causing hybrid sterility and inviability. These results can have some important implications for the spread of insecticide resistance and the control of these vector populations via the release of transgenic mosquitoes.

Malaria vectors in sub-Saharan Africa have proven themselves very difficult adversaries in the global struggle against malaria. Decades of anti-vector interventions have yielded mixed results--with successful reductions in transmission in some areas and limited impacts in others. These varying successes can be ascribed to a lack of universally effective vector control tools, as well as the development of insecticide resistance in mosquito populations. Understanding the impact of vector control on mosquito populations is crucial for planning new interventions and evaluating existing ones. However, estimates of population size changes in response to control efforts are often inaccurate because of limitations and biases in collection methods. Attempts to evaluate the impact of vector control on mosquito effective population size (N(e)) have produced inconclusive results thus far. Therefore, we obtained data for 13-15 microsatellite markers for more than 1,500 mosquitoes representing multiple time points for seven populations of three important vector species--Anopheles gambiae, An. melas, and An. moucheti--in Equatorial Guinea. These populations were exposed to indoor residual spraying or long-lasting insecticidal nets in recent years. For comparison, we also analyzed data from two populations that have no history of organized vector control. We used Approximate Bayesian Computation to reconstruct their demographic history, allowing us to evaluate the impact of these interventions on the effective population size. In six of the seven study populations, vector control had a dramatic impact on the effective population size, reducing N(e) between 55%-87%, the exception being a single An. melas population. In contrast, the two negative control populations did not experience a reduction in effective population size. This study is the first to conclusively link anti-vector intervention programs in Africa to sharply reduced effective population sizes of malaria vectors.