Ecological zones rather than molecular forms predict genetic differentiation in the malaria vector Anopheles gambiae s.s. in Ghana

Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS

Resistance to insecticides and adaptation to a diverse range of environments present challenges to Anopheles gambiae s.l. mosquito control efforts in sub-Saharan Africa. Whole-genome-sequencing is often employed for identifying the genomic basis underlying adaptation in Anopheles, but remains expensive for large-scale surveys. Reduced coverage whole-genome-sequencing can identify regions of the genome involved in adaptation at a lower cost, but is currently untested in Anopheles mosquitoes. Here, we use reduced coverage WGS to investigate population genetic structure and identify signatures of local adaptation in Anopheles mosquitoes across southern Ghana. In contrast to previous analyses, we find no structuring by ecoregion, with Anopheles coluzzii and Anopheles gambiae populations largely displaying the hallmarks of large, unstructured populations. However, we find signatures of selection at insecticide resistance loci that appear ubiquitous across ecoregions in An. coluzzii, and strongest in forest ecoregions in An. gambiae. Our study highlights resistance candidate genes in this region, and validates reduced coverage WGS, potentially to very low coverage levels, for population genomics and exploratory surveys for adaptation in Anopheles taxa.

Genetic Diversity, Haplotype Relationships, and kdr Mutation of Malaria Anopheles Vectors in the Most Plasmodium knowlesi-Endemic Area of Thailand

Plasmodium knowlesi, a malaria parasite that occurs naturally in long-tailed macaques, pig-tailed macaques, and banded leaf monkeys, is currently regarded as the fifth of the human malaria parasites. We aimed to investigate genetic diversity based on the cytochrome c oxidase subunit I (COI) gene, detect Plasmodium parasites, and screen for the voltage-gated sodium channel (VGSC)-mutation-mediated knockdown resistance (kdr) of Anopheles mosquitoes in Ranong province, which is the most P. knowlesi-endemic area in Thailand. One hundred and fourteen Anopheles females belonging to eight species, including An. baimaii (21.05%), An. minimus s.s. (20.17%), An. epiroticus (19.30%), An. jamesii (19.30%), An. maculatus s.s. (13.16%), An. barbirostris A3 (5.26%), An. sawadwongporni (0.88%), and An. aconitus (0.88%), were caught in three geographical regions of Ranong province. None of the Anopheles mosquitoes sampled in this study were infected with Plasmodium parasites. Based on the sequence analysis of COI sequences, An. epiroticus had the highest level of nucleotide diversity (0.012), followed by An. minimus (0.011). In contrast, An. maculatus (0.002) had the lowest level of nucleotide diversity. The Fu's Fs and Tajima's D values of the Anopheles species in Ranong were all negative, except the Tajima's D values of An. minimus (0.077). Screening of VGSC sequences showed no presence of the kdr mutation of Anopheles mosquitoes. Our results could be used to further select effective techniques for controlling Anopheles populations in Thailand's most P. knowlesi-endemic area.

Population genetic structure of the malaria vector Anopheles minimus in Thailand based on mitochondrial DNA markers

BACKGROUND

The malaria vector Anopheles minimus has been influenced by external stresses affecting the survival rate and vectorial capacity of the population. Since An. minimus habitats have continuously undergone ecological changes, this study aimed to determine the population genetic structure and the potential gene flow among the An. minimus populations in Thailand.

METHODS

Anopheles minimus was collected from five malaria transmission areas in Thailand using Centers for Disease Control and Prevention (CDC) light traps. Seventy-nine females from those populations were used as representative samples. The partial mitochondrial cytochrome c oxidase subunit I (COI), cytochrome c oxidase subunit II (COII) and cytochrome b (Cytb) gene sequences were amplified and analyzed to identify species and determine the current population genetic structure. For the past population, we determined the population genetic structure from the 60 deposited COII sequences in GenBank of An. minimus collected from Thailand 20 years ago.

RESULTS

The current populations of An. minimus were genetically divided into two lineages, A and B. Lineage A has high haplotype diversity under gene flow similar to the population in the past. Neutrality tests suggested population expansion of An. minimus, with the detection of abundant rare mutations in all populations, which tend to arise from negative selection.

CONCLUSIONS

This study revealed that the population genetic structure of An. minimus lineage A was similar between the past and present populations, indicating high adaptability of the species. There was substantial gene flow between the eastern and western An. minimus populations without detection of significant gene flow barriers.

Population structure and diversity of Plasmodium falciparum in children with asymptomatic malaria living in different ecological zones of Ghana

BACKGROUND

Genetic diversity in Plasmodium falciparum populations can be used to describe the resilience and spatial distribution of the parasite in the midst of intensified intervention efforts. This study used microsatellite analysis to evaluate the genetic diversity and population dynamics of P. falciparum parasites circulating in three ecological zones of Ghana.

METHODS

A total of 1168 afebrile children aged between 3 to 13 years were recruited from five (5) Primary schools in 3 different ecological zones (Sahel (Tamale and Kumbungu), Forest (Konongo) and Coastal (Ada and Dodowa)) of Ghana. Asymptomatic malaria parasite carriage was determined using microscopy and PCR, whilst fragment analysis of 6 microsatellite loci was used to determine the diversity and population structure of P. falciparum parasites.

RESULTS

Out of the 1168 samples examined, 16.1 and 39.5% tested positive for P. falciparum by microscopy and nested PCR respectively. The genetic diversity of parasites in the 3 ecological zones was generally high, with an average heterozygosity (He) of 0.804, 0.787 and 0.608 the rainy (peak) season for the Sahel, Forest and Coastal zones respectively. The mean He for the dry (off-peak) season were 0.562, 0.693 and 0.610 for the Sahel, Forest and Coastal zones respectively. Parasites from the Forest zone were more closely related to those from the Sahel than from the Coastal zone, despite the Coastal zone being closer in physical distance to the Forest zone. The fixation indexes among study sites ranged from 0.049 to 0.112 during the rainy season and 0.112 to 0.348 during the dry season.

CONCLUSION

A large asymptomatic parasite reservoir was found in the school children during both rainy and dry seasons, especially those in the Forest and Sahel savannah zones where parasites were also found to be related compared to those from the Coastal zone. Further studies are recommended to understand why despite the roll out of several malaria interventions in Ghana, high transmission still persist.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Genetic analysis and population structure of the Anopheles gambiae complex from different ecological zones of Burkina Faso

The Anopheles gambiae complex (Diptera: Culicidae) is the most important vector for malaria in Sub-Saharan Africa, besides other vectors such as Anopheles funestus. Malaria vector control should encompass specific identification, genetic diversity and population structure of An. gambiae to design vector control strategies. The aim of this study was to determine the distribution of sibling species of the An. gambiae complex according to climatic regions related to cotton-growing or cotton-free areas by using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Then, variation in mitochondrial cytochrome c oxidase 1 (COI) was used to assess the genetic structure within and between populations from our selected ecological zones. At the sibling species level, the following proportions were found across all samples (n = 180): An. coluzzii 65.56%, An. gambiae stricto sensu (s.s). 21.11%, and An. arabiensis 3.33%. Hybrids between An. gambiae s.s. and An. coluzzii (7.78%) and hybrids between An. coluzzii and An. arabiensis (2.22%) were found. The phylogenetic tree and Integer Neighbour-Joining (IntNJ) haplotype network did not reveal any distinct genetic structure pattern related to climatic or agricultural conditions in Burkina Faso. The Fst (Wright's F-statistic) values close to zero showed a free gene flow and no differentiation in An. gambiae complex populations. Furthermore, neutrality indices calculated by Tajima's D, Fu and Li's D^⁎, Fu and Li's F^⁎, Fu's Fs tests suggested an excess of rare mutations in the populations. Overall, variation in the proportions of An. gambiae s.s., An. coluzzii and An. arabiensis was found according to climatic regions, but COI analysis did not evidence any population structuring of the An. gambiae complex. These scientific contributions can be used as a basis for further in-depth study of the genetic diversity of the An. gambiae complex for epidemiological risk assessment of malaria in Burkina Faso.

Larval ecology of Anopheles coluzzii in Cape Coast, Ghana: water quality, nature of habitat and implication for larval control

Background

There is a growing interest in larval control intervention to supplement existing malaria control strategies, particularly in urban areas. However, effective implementation requires a good understanding of habitat ecology of Anopheles mosquitoes. Clean water bodies have long been reported by several studies as a preferred breeding habitat for Anopheles gambiae. Other studies have also reported the breeding of An. gambiae in polluted water bodies. However, the term clean or polluted is mostly based on visual examination and is not well defined. This study was conducted with the aim of assessing water quality in Anopheles breeding habitats and the practicability of larval control in Cape Coast, Ghana.

Methods

A larval survey was conducted for 15 months in Cape Coast. In individual breeding habitats, habitat characteristics, physicochemical parameters and bacterial fauna were measured in both Anopheles positive breeding (APL) habitats and habitats colonized by only Culex species. The sibling species of An. gambiae were identified using PCR assay.

Results

Anopheles coluzzii dominated in almost all the APL habitats found in this study. The habitats had high levels of salinity and ammonium ions. However, ammonium ions were significantly higher (p = 0.001) in habitats colonized by only Culex larvae compared to APL habitats. About 47 % of the habitats that were colonized by only Culex larvae had no measurable dissolved oxygen while An. coluzzii was absent in such habitats. High concentration of faecal bacteria confirmed faecal contamination in both groups of breeding habitats.

Conclusions

From the results, it was evident that larval stages of An. coluzzii have tolerance to high levels of salinity and organic pollution in breeding habitats. However, its level of tolerance to organic pollution is probably lower than Culex larvae. The nature of breeding habitats found in the city demonstrates the opportunistic behaviour of An. coluzzii and how its breeding requirements are so intimately intertwined with the haphazard and uncontrolled human activities in the urban area. Considering the nature of APL habitats, larval control intervention could greatly reduce Anopheles population. However, improving basic hygiene and sanitation in the city could even make larval control intervention more practical and cost effective.

Habitat segregation and ecological character displacement in cryptic African malaria mosquitoes

Understanding how divergent selection generates adaptive phenotypic and population diversification provides a mechanistic explanation of speciation in recently separated species pairs. Towards this goal, we sought ecological gradients of divergence between the cryptic malaria vectors Anopheles coluzzii and An. gambiae and then looked for a physiological trait that may underlie such divergence. Using a large set of occurrence records and eco-geographic information, we built a distribution model to predict the predominance of the two species across their range of sympatry. Our model predicts two novel gradients along which the species segregate: distance from the coastline and altitude. Anopheles coluzzii showed a ‘bimodal’ distribution, predominating in xeric West African savannas and along the western coastal fringe of Africa. To test whether differences in salinity tolerance underlie this habitat segregation, we assessed the acute dose–mortality response to salinity of thirty-two larval populations from Central Africa. In agreement with its coastal predominance, Anopheles coluzzii was overall more tolerant than An. gambiae. Salinity tolerance of both species, however, converged in urban localities, presumably reflecting an adaptive response to osmotic stress from anthropogenic pollutants. When comparing degree of tolerance in conjunction with levels of syntopy, we found evidence of character displacement in this trait.

Insecticide susceptibility of Anopheles coluzzii and Anopheles gambiae mosquitoes in Ibadan, Southwest Nigeria

The emergence of insecticide resistance in Anopheles (Diptera: Culicidae) mosquitoes has great implications for malaria control in Nigeria. This study aimed to determine the dynamics of insecticide susceptibility levels and the frequency of knock-down resistance (kdr) mutations (L1014F) in wild Anopheles coluzzii Coetzee & Wilkerson sp. n. and Anopheles gambiae Giles from the Ojoo and Bodija areas of Ibadan, in southwest Nigeria. Insecticide susceptibility to pyrethroids, organophosphates, carbamates and organochlorines was assessed using World Health Organization (WHO) bioassays. A subset of the mosquitoes exposed to pyrethroids and DDT was used for species and molecular form identification; kdr genotyping was determined using the TaqMan real-time polymerase chain reaction assay. The mosquitoes were resistant to pyrethroids and DDT but completely susceptible to organophosphates and carbamates. Bodija samples (n = 186) consisted of An. gambiae (91.4%) and An. coluzzii (8.1%) and included one An. coluzzii/An. gambiae hybrid specimen. All mosquitoes screened in Ojoo (n = 26) were An. gambiae. The 1014F kdr mutation was detected at frequencies of 24.5 and 5.8% in Bodija and Ojoo, respectively. No correlation was observed between kdr genotypes and resistance phenotypes. The results indicate that metabolic resistance probably plays an important role in the development of resistance and highlight the need to implement insecticide resistance management strategies.

Islands and stepping-stones: comparative population structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and implications for the spread of insecticide resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise FST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average FST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.

Effects of the kdr resistance mutation on the susceptibility of wild Anopheles gambiae populations to Plasmodium falciparum: a hindrance for vector control

Background

In the context of generalization of insecticide resistance, the hypothesis that insecticide resistance has a positive impact on the capacity of mosquitoes to transmit malaria constitutes a hindrance for malaria elimination. The aim of this study was to investigated populations of Anopheles coluzzii and Anopheles gambiae S molecular form to assess whether different genotypes at the kdr locus are responsible for different susceptibility to Plasmodium falciparum infection.

Methods

F3 progeny of An. gambiae s.l. collected in Dielmo were infected by direct membrane feeding with P. falciparum gametocyte-containing blood sampled from volunteer patients. The presence of oocysts was determined by light microscopy after seven days, and the presence of sporozoites by ELISA after 14 days. Mosquito species and molecular forms were identified by PCR. Generalized linear models were performed using the R software to test the effect of explanatory variables including the genotype at the kdr locus on infection rate and density.

Results

The odds of being infected with oocysts and sporozoites were greater in RS and RR groups than in SS groups (χ² = 42.8, df = 1, P(>χ²) = 6.1e-11). The density of infection was also dependent on genotype, with RR and RS genotypes showing denser infection than SS genotypes. Pairwise comparisons of oocyst number and absorbance indicated sometime a small betwen species (i.e. between An. gambiae S form, and An. coluzzii), but the effect of genotype was much more important.

Conclusion

The presence of the resistance allele at the kdr locus increases susceptibility to Plasmodium not only at the oocyst stage but also at the sporozoite stage in non-genetically modified wild mosquitoes. These results have significant implications and should be taken into account in the development of strategies for malaria control.

Prominent intraspecific genetic divergence within Anopheles gambiae sibling species triggered by habitat discontinuities across a riverine landscape

The Anopheles gambiae complex of mosquitoes includes malaria vectors at different stages of speciation, whose study enables a better understanding of how adaptation to divergent environmental conditions leads to evolution of reproductive isolation. We investigated the population genetic structure of closely related sympatric taxa that have recently been proposed as separate species (An. coluzzii and An. gambiae), sampled from diverse habitats along the Gambia river in West Africa. We characterized putatively neutral microsatellite loci as well as chromosomal inversion polymorphisms known to be associated with ecological adaptation. The results revealed strong ecologically associated population subdivisions within both species. Microsatellite loci on chromosome-3L revealed clear differentiation between coastal and inland populations, which in An. coluzzii is reinforced by a unusual inversion polymorphism pattern, supporting the hypothesis of genetic divergence driven by adaptation to the coastal habitat. A strong reduction of gene flow was observed between An. gambiae populations west and east of an extensively rice-cultivated region apparently colonized exclusively by An. coluzzii. Notably, this 'intraspecific' differentiation is higher than that observed between the two species and involves also the centromeric region of chromosome-X which has previously been considered a marker of speciation within this complex, possibly suggesting that the two populations may be at an advanced stage of differentiation triggered by human-made habitat fragmentation. These results confirm ongoing ecological speciation within these most important Afro-tropical malaria vectors and raise new questions on the possible effect of this process in malaria transmission.

Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation

Adaptive introgression can provide novel genetic variation to fuel rapid evolutionary responses, though it may be counterbalanced by potential for detrimental disruption of the recipient genomic background. We examine the extent and impact of recent introgression of a strongly selected insecticide-resistance mutation (Vgsc-1014F) located within one of two exceptionally large genomic islands of divergence separating the Anopheles gambiae species pair. Here we show that transfer of the Vgsc mutation results in homogenization of the entire genomic island region (~1.5% of the genome) between species. Despite this massive disruption, introgression is clearly adaptive with a dramatic rise in frequency of Vgsc-1014F and no discernable impact on subsequent reproductive isolation between species. Our results show (1) how resilience of genomes to massive introgression can permit rapid adaptive response to anthropogenic selection and (2) that even extreme prominence of genomic islands of divergence can be an unreliable indicator of importance in speciation.

Highly divergent genomic islands segregate between a species pair of the mosquito, Anopheles gambiae. Here Clarkson et al. show that loss of one of the largest islands, driven by adaptive introgression of an insecticide-resistance mutation, had no impact on reproductive isolation.

Acetylcholinesterase (Ace-1) target site mutation 119S is strongly diagnostic of carbamate and organophosphate resistance in Anopheles gambiae s.s. and Anopheles coluzzii across southern Ghana

BACKGROUND

With high DDT resistance present throughout much of West Africa, carbamates and organophosphates are increasingly important alternatives to pyrethroids for indoor residual spraying (IRS). Though less widespread, resistance to both of these alternative insecticide classes has also been documented within the Anopheles gambiae species pair (formerly the M and S molecular forms) in West Africa. To manage insecticide efficacy, it is important to predict how and where resistance is likely to occur and spread, which could be aided by using molecular diagnostics with high predictive value.

METHODS

Anopheles coluzzii and An. gambiae s.s. were collected from 18 sites throughout southern Ghana and bioassayed with bendiocarb, the most commonly applied carbamate, and an organophosphate, fenitrothion. The Ace-1 target site substitution G119S was genotyped by qPCR.

RESULTS

Fenitrothion induced higher mortality than bendiocarb, though phenotypes correlated strongly across populations. Ace-1 119S was found at much higher frequency in An. gambiae s.s than An. coluzzii, exceeding 90% in a population from Greater Accra, the highest frequency reported to date. Ace-1 G119S was very strongly associated with resistance to both insecticides, providing high predictive power for diagnosis, though with some evidence for a differential effect between molecular forms for bendiocarb. Sequencing of the gene revealed a lack of variation in resistant alleles precluding determination of origin, but Ace-1 copy number variation was detected for the first time in Ghana.

CONCLUSIONS

The results validate G119S as a useful diagnostic of organophosphate and carbamate resistance within and among populations, whilst highlighting the potential for an aggregate nature of Ace-1 genotypes, which may comprise both single-copy and duplicated genes. Further work is now required to determine the distribution and resistance-association of Ace-1 duplication.

Insect speciation rules: unifying concepts in speciation research

The study of speciation is concerned with understanding the connection between causes of divergent evolution and the origin and maintenance of barriers to gene exchange between incipient species. Although the field has historically focused either on examples of recent divergence and its causes or on the genetic basis of reproductive isolation between already divergent species, current efforts seek to unify these two approaches. Here we integrate these perspectives through a discussion of recent progress in several insect speciation model systems. We focus on the evolution of speciation phenotypes in each system (i.e., those phenotypes causally involved in reducing gene flow between incipient species), drawing an explicit connection between cause and effect (process and pattern). We emphasize emerging insights into the genomic architecture of speciation as well as timely areas for future research.

Geographic population structure of the African malaria vector Anopheles gambiae suggests a role for the forest-savannah biome transition as a barrier to gene flow

The primary Afrotropical malaria mosquito vector Anopheles gambiae sensu stricto has a complex population structure. In west Africa, this species is split into two molecular forms and displays local and regional variation in chromosomal arrangements and behaviors. To investigate patterns of macrogeographic population substructure, 25 An. gambiae samples from 12 African countries were genotyped at 13 microsatellite loci. This analysis detected the presence of additional population structuring, with the M-form being subdivided into distinct west, central, and southern African genetic clusters. These clusters are coincident with the central African rainforest belt and northern and southern savannah biomes, which suggests restrictions to gene flow associated with the transition between these biomes. By contrast, geographically patterned population substructure appears much weaker within the S-form.

Physiological correlates of ecological divergence along an urbanization gradient: differential tolerance to ammonia among molecular forms of the malaria mosquito Anopheles gambiae

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Feeding frequency and survival of Anopheles gambiae in a rice-growing area in Ghana

Mortality rates, determined by dissection, of predominantly M form female Anopheles gambiae (Diptera: Culicidae) were estimated. Mosquitoes were collected in tent traps and light traps in an irrigation project village in Ghana in June and July 2010, when much of the area was flooded. Both M and S form larvae were collected from rice fields (74 of 80 specimens were M form). Adults were collected in equal proportions from the two traps (90 of 107 specimens from the light trap and 106 of 116 specimens from the tent trap were M form). During the study, collection numbers rose from 105 to 972 per night. A total of 1787 of the 15 431 An. gambiae collected were dissected. Of these, 953 (53%) were found to have taken their first bloodmeal, either as virgins or following mating. The age profiles of mosquitoes collected alive and dead, respectively, were similar. Eighteen of 2933 (0.61 ± 0.49%) specimens were found to be positive for sporozoites in an enzyme-linked immunosorbent assay (ELISA). Lagged cross correlations among the different age groups implied that the mosquitoes fed on days 2 and 4 following emergence prior to oviposition and every 2.65 ± 0.17 days thereafter. The best model to describe the observed population patterns implied a daily mortality of 84%. The results are discussed in relation to possible mosquito control measures for the village.

Anthropogenic habitat disturbance and ecological divergence between incipient species of the malaria mosquito Anopheles gambiae

BACKGROUND

Anthropogenic habitat disturbance is a prime cause in the current trend of the Earth's reduction in biodiversity. Here we show that the human footprint on the Central African rainforest, which is resulting in deforestation and growth of densely populated urban agglomerates, is associated to ecological divergence and cryptic speciation leading to adaptive radiation within the major malaria mosquito Anopheles gambiae.

METHODOLOGY/PRINCIPAL FINDINGS

In southern Cameroon, the frequency of two molecular forms--M and S--among which reproductive isolation is strong but still incomplete, was correlated to an index of urbanisation extracted from remotely sensed data, expressed as the proportion of built-up surface in each sampling unit. The two forms markedly segregated along an urbanisation gradient forming a bimodal cline of ∼6-km width: the S form was exclusive to the rural habitat, whereas only the M form was present in the core of densely urbanised settings, co-occurring at times in the same polluted larval habitats of the southern house mosquito Culex quinquefasciatus--a species association that was not historically recorded before.

CONCLUSIONS/SIGNIFICANCE

Our results indicate that when humans create novel habitats and ecological heterogeneities, they can provide evolutionary opportunities for rapid adaptive niche shifts associated with lineage divergence, whose consequences upon malaria transmission might be significant.

Studies on the behaviour of peridomestic and endophagic M form Anopheles gambiae from a rice growing area of Ghana

The 'paddy paradox', the occurrence of large populations of vectors but low amounts of malaria transmission where irrigated rice is grown, was investigated in a village in Ghana where M form Anopheles gambiae are common. Peridomestic and indoor host-seeking mosquitoes were collected in tent traps and light traps over 21 consecutive nights at the start of the rainy season in June 2009 when the population increased exponentially from less than 100 per night to over 1000. Infection rates in the overall mosquito population were 0.3% and in the estimated parous population were 1.9%. Numbers of An. gambiae in the tent trap peaked between midnight and 02:40 am. The majority of insects were taking their first blood meal, as virgins or shortly after mating. More than expected were collected in the light trap during a rainstorm at the start of the rains but overall numbers were not affected. Fewer than expected were collected after a subsequent storm. Recruitment to the adult population decreased over the following days. It is hypothesised that the 'paddy paradox' is due to young pre-gravid insects dispersing more widely than gravid ones, not necessarily to low survival in the mosquito.

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.

Studies of Anopheles gambiae s.l (Diptera: Culicidae) exhibiting different vectorial capacities in lymphatic filariasis transmission in the Gomoa district, Ghana

Background

Two lymphatic filariasis endemic communities Mampong and Hwida in Ghana have been regularly monitored for impact on transmission after annual mass drug administration (MDA) with albendazole and ivermectin. After six MDAs even though the ABR for Mampong was 55883/person/year and that of Hwida was 2494/person/year, they both had ATPs of 15.21 infective larvae/person/year. Interestingly the human microfilaraemia levels had reduced significantly from 14% to 0% at Mampong and 12% to 3% at Hwida. In an attempt to understand this anomaly, we collected mosquitoes over a 5-month period using human landing catches to determine the species composition, the number of cibarial teeth, the lengths and widths of the cibarium and the cibarial dome of the vector populations.

Results

Out of 2553 mosquitoes caught at Mampong, 42.6% were An. gambiae s.l. All 280 identified further by PCR were An. gambiae s.s (275 M and 5 S molecular forms). At Hwida, 112 mosquitoes were obtained; 67 (59.8%) were An. gambiae s.l, comprised of 40 (59.7%) An. melas, 24 (35.8%) An. gambiae s.s (17 and 5 M and S molecular forms respectively) and 3 (4.5%) unidentified. The mean number of teeth for An. melas was 14.1 (median = 14, range = 12-15), An. gambiae s.s., 15.7 (median = 15, range = 13-19) M form 15.5 (median = 15 range = 13-19) and S form 16 (median = 16, range 15-17). The observed differences in teeth numbers were significantly different between An. melas and An. gambiae s.s (p = 0.004), and the M form (p = 0.032) and the S form (p = 0.002).

Conclusions

In this study, An. gambiae s.s was the main vector at Mampong and was found to possess significantly more cibarial teeth than An. melas, the principal vector at Hwida. We postulate that the different impact observed after 6 MDAs may be due to An. gambiae s.s exhibiting 'facilitation' at Mampong and at Hwida An. melas the main vector exhibits 'limitation'. Thus it may be necessary to compliment MDA with vector control to achieve interruption of transmission in areas where An. melas may exhibit limitation.

Population genetic structure of Anopheles gambiae and Anopheles arabiensis in Niger

The increasing usage of long-lasting insecticide-treated nets allows protection of millions of people from malaria infection. Monitoring studies should be planned during any wide-scale malaria control program integrating insecticide-treated materials, to evaluate their effects and effectiveness on epidemiologically relevant parameters. Such operational control interventions may be challenged by insecticide resistance spread within vector populations, as a result of wide insecticide pressure. A nationwide distribution of long-lasting insecticidal nets was implemented throughout Niger in 2005. We studied the population genetic structure of major malaria vectors across Nigerien Sahel, and investigated potential effects of this large malaria control intervention. Wild-caught Anopheles gambiae sensu lato females from seven villages and two wet seasons were genotyped at 12 microsatellite loci. The genetic diversity within both species appeared homogenous between villages and years. The estimated genetic differentiation among samples was very low within both species, indicating high gene flow across the area. An absence of differentiation was also found between 2005 and 2006 wet seasons, for all samples but one, showing that the net distribution did not impact significantly the genetic diversity and structure of vector populations in a single year. We provide valuable results participating to document effects of large malaria control programs, to maximize the efficiency of available tools in future interventions.

Environmental factors associated with the distribution of Anopheles gambiae s.s in Ghana; an important vector of lymphatic filariasis and malaria

Anopheles gambiae s.s mosquitoes are important vectors of lymphatic filariasis (LF) and malaria in Ghana. To better understand their ecological aspects and influence on disease transmission, we examined the spatial distribution of the An. gambiae (M and S) molecular forms and associated environmental factors, and determined their relationship with disease prevalence. Published and current data available on the An. gambiae species in Ghana were collected in a database for analysis, and the study sites were georeferenced and mapped. Using the An. gambiae s.s sites, environmental data were derived from climate, vegetation and remote-sensed satellite sources, and disease prevalence data from existing LF and malaria maps in the literature. The data showed that An. gambiae M and S forms were sympatric in most locations. However, the S form predominated in the central region, while the M form predominated in the northern and coastal savanna regions. Bivariate and multiple regression analyses identified temperature as a key factor distinguishing their distributions. An. gambiae M was significantly correlated with LF, and 2.5 to 3 times more prevalent in the high LF zone than low to medium zones. There were no significant associations between high prevalence An. gambiae s.s locations and malaria. The distribution of the An. gambiae M and S forms and the diseases they transmit in Ghana appear to be distinct, driven by different environmental factors. This study provides useful baseline information for disease control, and future work on the An. gambiae s.s in Ghana.

Field, genetic, and modeling approaches show strong positive selection acting upon an insecticide resistance mutation in Anopheles gambiae s.s

Alleles subject to strong, recent positive selection will be swept toward fixation together with contiguous sections of the genome. Whether the genomic signatures of such selection will be readily detectable in outbred wild populations is unclear. In this study, we employ haplotype diversity analysis to examine evidence for selective sweeps around knockdown resistance (kdr) mutations associated with resistance to dichlorodiphenyltrichloroethane and pyrethroid insecticides in the mosquito Anopheles gambiae. Both kdr mutations have significantly lower haplotype diversity than the wild-type (nonresistant) allele, with kdr L1014F showing the most pronounced footprint of selection. We complement these data with a time series of collections showing that the L1014F allele has increased in frequency from 0.05 to 0.54 in 5 years, consistent with a maximum likelihood-fitted selection coefficient of 0.16 and a dominance coefficient of 0.25. Our data show that strong, recent positive selective events, such as those caused by insecticide resistance, can be identified in wild insect populations.

Inferring selection in the Anopheles gambiae species complex: an example from immune-related serine protease inhibitors

Background

Mosquitoes of the Anopheles gambiae species complex are the primary vectors of human malaria in sub-Saharan Africa. Many host genes have been shown to affect Plasmodium development in the mosquito, and so are expected to engage in an evolutionary arms race with the pathogen. However, there is little conclusive evidence that any of these mosquito genes evolve rapidly, or show other signatures of adaptive evolution.

Methods

Three serine protease inhibitors have previously been identified as candidate immune system genes mediating mosquito-Plasmodium interaction, and serine protease inhibitors have been identified as hot-spots of adaptive evolution in other taxa. Population-genetic tests for selection, including a recent multi-gene extension of the McDonald-Kreitman test, were applied to 16 serine protease inhibitors and 16 other genes sampled from the An. gambiae species complex in both East and West Africa.

Results

Serine protease inhibitors were found to show a marginally significant trend towards higher levels of amino acid diversity than other genes, and display extensive genetic structuring associated with the 2La chromosomal inversion. However, although serpins are candidate targets for strong parasite-mediated selection, no evidence was found for rapid adaptive evolution in these genes.

Conclusion

It is well known that phylogenetic and population history in the An. gambiae complex can present special problems for the application of standard population-genetic tests for selection, and this may explain the failure of this study to detect selection acting on serine protease inhibitors. The pitfalls of uncritically applying these tests in this species complex are highlighted, and the future prospects for detecting selection acting on the An. gambiae genome are discussed.

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.

Ecological and genetic relationships of the Forest-M form among chromosomal and molecular forms of the malaria vector Anopheles gambiae sensu stricto

BACKGROUND

Anopheles gambiae sensu stricto, one of the principal vectors of malaria, has been divided into two subspecific groups, known as the M and S molecular forms. Recent studies suggest that the M form found in Cameroon is genetically distinct from the M form found in Mali and elsewhere in West Africa, suggesting further subdivision within that form.

METHODS

Chromosomal, microsatellite and geographic/ecological evidence are synthesized to identify sources of genetic polymorphism among chromosomal and molecular forms of the malaria vector Anopheles gambiae s.s.

RESULTS

Cytogenetically the Forest M form is characterized as carrying the standard chromosome arrangement for six major chromosomal inversions, namely 2La, 2Rj, 2Rb, 2Rc, 2Rd, and 2Ru. Bayesian clustering analysis based on molecular form and chromosome inversion polymorphisms as well as microsatellites describe the Forest M form as a distinct population relative to the West African M form (Mopti-M form) and the S form. The Forest-M form was the most highly diverged of the An. gambiae s.s. groups based on microsatellite markers. The prevalence of the Forest M form was highly correlated with precipitation, suggesting that this form prefers much wetter environments than the Mopti-M form.

CONCLUSION

Chromosome inversions, microsatellite allele frequencies and habitat preference all indicate that the Forest M form of An. gambiae is genetically distinct from the other recognized forms within the taxon Anopheles gambiae sensu stricto. Since this study covers limited regions of Cameroon, the possibility of gene flow between the Forest-M form and Mopti-M form cannot be rejected. However, association studies of important phenotypes, such as insecticide resistance and refractoriness against malaria parasites, should take into consideration this complex population structure.

High levels of hybridization between molecular forms of Anopheles gambiae from Guinea Bissau

In the malaria vector Anopheles gambiae Giles sensu stricto, two molecular forms denoted M and S are considered units of incipient speciation within this species. Very low hybrid frequencies and significant genetic differentiation have been found in sympatric M- and S-form populations. We studied the molecular form composition and the degree of genetic differentiation at 15 microsatellites in two samples of An. gambiae collected in two consecutive years from Bissau, Guinea Bissau. High frequencies of M/S hybrids (19-24%) were found in this area. Coincidently, very low levels of genetic differentiation were detected between forms when analysis involved microsatellites mapped at chromosome-3 (mean Fst, 0.000-0.002). The single exception was the X-linked AGXH678, for which high differentiation was measured (Fst, 0.158-0.301). This locus maps near the centromere of chromosome X, a low recombination region in which selection is likely to promote divergence between M and S forms. These results strongly suggest that the degree of isolation between M and S forms, considered the units of incipient speciation within An. gambiae, is not homogenous throughout the species distribution range.

Dynamics of transmission of Plasmodium falciparum by Anopheles arabiensis and the molecular forms M and S of Anopheles gambiae in Dielmo, Senegal

Background

The adaptation of Anopheles gambiae to humans and its environment involves an ongoing speciation process that can be best demonstrated by the existence of various chromosomal forms adapted to different environments and of two molecular forms known as incipient taxonomic units.

Methods

The aim of this study was to compare the epidemiologic role of Anopheles arabiens is and the molecular forms M and S of Anopheles gambiae in the transmission of Plasmodium in a rural areas of southern Senegal, Dielmo. The sampling of mosquitoes was carried out monthly between July and December 2004, during the rainy season, by human volunteers and pyrethrum spray catches.

Results

Anopheles arabiensis, An. gambiae M and S forms coexisted during the rainy season with a predominance of the M form in September and the peak of density being observed in August for the S form. Similar parity rates were observed in An. arabiensis [70.9%] (n = 86), An. gambiae M form [68.7%] (n = 64) and An. gambiae S form [81.1%] (n = 156). The circumsporozoite protein (CSP) rates were 2.82% (n = 177), 3.17% (n = 315) and 3.45% (n = 405), with the mean anthropophilic rates being 71.4% (n = 14), 86.3% (n = 22) and 91.6% (n = 24) respectively for An. arabiensis and An. gambiae M and S forms. No significant difference was observed either in host preference or in Plasmodium falciparum infection rates between sympatric M and S populations.

Conclusion

No difference was observed either in host preference or in Plasmodium falciparum infection rates between sympatric M and S populations, but they present different dynamics of population. These variations are probably attributable to different breeding conditions.

Assessment of alternative mating strategies in Anopheles gambiae: Does mating occur indoors?

Mating in Anopheles gambiae has been observed only in outdoor swarms. Here we evaluate whether mating also occurs indoors. Mark-release-recapture of virgin males and females in natural houses showed that mating occurred over a single day even when mosquitoes can leave the house through exit traps and without adaptation to laboratory conditions. In these experiments, insemination rate in the M molecular form of An. gambiae (and An. arabiensis) was higher than that of the S form (15 versus 6%). Under these conditions, smaller females of the M form mated more frequently than larger females of that form. Sampling mosquitoes throughout the day showed that both sexes enter houses around sunrise and leave around sunset, staying indoors together from dawn to dusk. In an area dominated by the M form, the daily rate of insemination in samples from exit traps was approximately 5% higher than in those from entry traps, implying that mating occurred indoors. Importantly, frequency of cross mating between the molecular forms was as high as that between members of the same form, indicating that, indoors, assortative mating breaks down. Altogether, these results suggest that indoor mating is an alternative mating strategy of the M molecular form of An. gambiae. Because naturally occurring mating couples have not yet been observed indoors, this conclusion awaits validation.

High genetic differentiation between the M and S molecular forms of Anopheles gambiae in Africa

Background

Anopheles gambiae, a major vector of malaria, is widely distributed throughout sub-Saharan Africa. In an attempt to eliminate infective mosquitoes, researchers are trying to develop transgenic strains that are refractory to the Plasmodium parasite. Before any release of transgenic mosquitoes can be envisaged, we need an accurate picture of the differentiation between the two molecular forms of An. gambiae, termed M and S, which are of uncertain taxonomic status.

Methodology/Principal Findings

Insertion patterns of three transposable elements (TEs) were determined in populations from Benin, Burkina Faso, Cameroon, Ghana, Ivory Coast, Madagascar, Mali, Mozambique, Niger, and Tanzania, using Transposon Display, a TE-anchored strategy based on Amplified Fragment Length Polymorphism. The results reveal a clear differentiation between the M and S forms, whatever their geographical origin, suggesting an incipient speciation process.

Conclusions/Significance

Any attempt to control the transmission of malaria by An. gambiae using either conventional or novel technologies must take the M/S genetic differentiation into account. In addition, we localized three TE insertion sites that were present either in every individual or at a high frequency in the M molecular form. These sites were found to be located outside the chromosomal regions that are suspected of involvement in the speciation event between the two forms. This suggests that these chromosomal regions are either larger than previously thought, or there are additional differentiated genomic regions interspersed with undifferentiated regions.

X and Y chromosome inheritance and mixtures of rDNA intergenic spacer regions in Anopheles gambiae

We observed Anopheles gambiae sensu stricto stocks that contained both Mopti (M) and Savanna (S) rDNA intergenic spacers (IGS). ASEMBO1 male IGS sequences consistently had a mixture. A diagnostic M and S Hha I restriction enzyme site in these fragments was concordant with two SNPs associated with M and S. Standard M and S stocks demonstrated X-chromosome-only inheritance of the rDNA form, but the ASEMBO1 males showed X and Y chromosome linkage of mixed rDNA. The metaphase Y chromosomes of ASEMBO1 contained a significantly larger amount of DNA relative to the X than a standard S stock. Analysis of wild A. gambiae males from the putative location of origin of the ASEMBO1 stock did not detect the same pattern of polymorphism observed in the laboratory stock but did detect heterogeneous arrays including some missing a diagnostic Hha I restriction site. These results demonstrate that M and S IGS types can occur within the rDNA arrays of a chromatid in laboratory A. gambiae stocks, and some A. gambiae s.s. have rDNA on the Y chromosome.

Multiple origins of knockdown resistance mutations in the Afrotropical mosquito vector Anopheles gambiae

How often insecticide resistance mutations arise in natural insect populations is a fundamental question for understanding the evolution of resistance and also for modeling its spread. Moreover, the development of resistance is regarded as a favored model to study the molecular evolution of adaptive traits. In the malaria vector Anopheles gambiae two point mutations (L1014F and L1014S) in the voltage-gated sodium channel gene, that confer knockdown resistance (kdr) to DDT and pyrethroid insecticides, have been described. In order to determine whether resistance alleles result from single or multiple mutation events, genotyping of the kdr locus and partial sequencing of the upstream intron-1 was performed on a total of 288 A. gambiae S-form collected from 28 localities in 15 countries. Knockdown resistance alleles were found to be widespread in West Africa with co-occurrence of both 1014S and 1014F in West-Central localities. Differences in intron-1 haplotype composition suggest that kdr alleles may have arisen from at least four independent mutation events. Neutrality tests provided evidence for a selective sweep acting on this genomic region, particularly in West Africa. The frequency and distribution of these kdr haplotypes varied geographically, being influenced by an interplay between different mutational occurrences, gene flow and local selection. This has important practical implications for the management and sustainability of malaria vector control programs.

Genetic population structure of Anopheles gambiae in Equatorial Guinea

BACKGROUND

Patterns of genetic structure among mosquito vector populations in islands have received particular attention as these are considered potentially suitable sites for experimental trials on transgenic-based malaria control strategies. In this study, levels of genetic differentiation have been estimated between populations of Anopheles gambiae s.s. from the islands of Bioko and Annobón, and from continental Equatorial Guinea (EG) and Gabon.

METHODS

Genotyping of 11 microsatellite loci located in chromosome 3 was performed in three island samples (two in Bioko and one in Annobón) and three mainland samples (two in EG and one in Gabon). Four samples belonged to the M molecular form and two to the S-form. Microsatellite data was used to estimate genetic diversity parameters, perform demographic equilibrium tests and analyse population differentiation.

RESULTS

High levels of genetic differentiation were found between the more geographically remote island of Annobón and the continent, contrasting with the shallow differentiation between Bioko island, closest to mainland, and continental localities. In Bioko, differentiation between M and S forms was higher than that observed between island and mainland samples of the same molecular form.

CONCLUSION

The observed patterns of population structure seem to be governed by the presence of both physical (the ocean) and biological (the M-S form discontinuity) barriers to gene flow. The significant degree of genetic isolation between M and S forms detected by microsatellite loci located outside the "genomic islands" of speciation identified in A. gambiae s.s. further supports the hypothesis of on-going incipient speciation within this species. The implications of these findings regarding vector control strategies are discussed.

Population genetics of Plasmodium resistance genes in Anopheles gambiae: no evidence for strong selection

Anopheles mosquitoes are the primary vectors for malaria in Africa, transmitting the disease to more than 100 million people annually. Recent functional studies have revealed mosquito genes that are crucial for Plasmodium development, but there is presently little understanding of which genes mediate vector competence in the wild, or evolve in response to parasite-mediated selection. Here, we use population genetic approaches to study the strength and mode of natural selection on a suite of mosquito immune system genes, CTL4, CTLMA2, LRIM1, and APL2 (LRRD7), which have been shown to affect Plasmodium development in functional studies. We sampled these genes from two African populations of An. gambiae s.s., along with several closely related species, and conclude that there is no evidence for either strong directional or balancing selection on these genes. We highlight a number of challenges that need to be met in order to apply population genetic tests for selection in Anopheles mosquitoes; in particular the dearth of suitable outgroup species and the potential difficulties that arise when working within a closely-related species complex.