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

First report of the L1014S kdr mutation in wild populations of Anopheles gambiae M and S molecular forms in Burkina Faso (West Africa)

We investigated the occurrence of the L1014F and L1014S kdr mutations in malaria vector populations in Burkina Faso (West Africa). A cross-sectional survey was conducted at 10 sites all located in cotton cultivation areas which are assumed to be the major insecticide resistance selection foci in Burkina Faso. The hot ligation method was used to detect the two kdr mutations in field collected Anopheles gambiae s.l. samples. For the first time in Burkina Faso the L1014S mutation was identified in both M and S forms of An. gambiae s.s. populations collected from the site of Koupela in the central-eastern region at low frequency. Furthermore, the L1014S mutation was also found in one specimen of An. arabiensis collected from the Dano site. The data generated in this study provides additional evidence of the spread of the L1014S mutation into An. gambiae s.l. populations in West Africa. It is now important to evaluate the role of the L1014S mutation in the pyrethroid resistance phenotype and assess its potential impact on the efficacy of pyrethroid-based control measures in West Africa where several resistance mutations now coexist.

Application of a qPCR assay in the investigation of susceptibility to malaria infection of the M and S molecular forms of An. gambiae s.s. in Cameroon

Plasmodium falciparum is the causative agent of malaria, a disease that kills almost one million persons each year, mainly in sub-Saharan Africa. P. falciparum is transmitted to the human host by the bite of an Anopheles female mosquito, and Anopheles gambiae sensus stricto is the most tremendous malaria vector in Africa, widespread throughout the afro-tropical belt. An. gambiae s.s. is subdivided into two distinct molecular forms, namely M and S forms. The two molecular forms are morphologically identical but they are distinct genetically, and differ by their distribution and their ecological preferences. The epidemiological importance of the two molecular forms in malaria transmission has been poorly investigated so far and gave distinct results in different areas. We have developed a real-time quantitative PCR (qPCR) assay, and used it to detect P. falciparum at the oocyst stage in wild An. gambiae s.s. mosquitoes experimentally infected with natural isolates of parasites. Mosquitoes were collected at immature stages in sympatric and allopatric breeding sites and further infected at the adult stage. We next measured the infection prevalence and intensity in female mosquitoes using the qPCR assay and correlated the infection success with the mosquito molecular forms. Our results revealed different prevalence of infection between the M and S molecular forms of An. gambiae s.s. in Cameroon, for both sympatric and allopatric populations of mosquitoes. However, no difference in the infection intensity was observed. Thus, the distribution of the molecular forms of An. gambiae s.s. may impact on the malaria epidemiology, and it will be important to monitor the efficiency of malaria control interventions on the two M and S forms.

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.

Longitudinal follow-up of malaria transmission dynamics in two villages in a Sahelian area of Niger during a nationwide insecticide-treated bednet distribution programme

Malaria transmission was monitored in two villages in the Sahel zone of Niger over 4 years. During this period, a nationwide vector control programme was carried out in which insecticide-treated bednets were distributed free to mothers of children aged <5 years. Anopheles gambiae and Anopheles arabiensis (Diptera: Culicidae) were found to be the major malaria vectors. The dynamics of An. gambiae s.l. did not vary dramatically over the study period although the proportion of female mosquitoes found resting indoors decreased in both villages and, in one village, the parity rate and sporozoite index were significantly reduced after bednet distribution. By contrast with An. gambiae, the dynamics of Anopheles funestus altered greatly after the bednet distribution period, when adult density, endophagous rate and sporozoite rates decreased dramatically. Our observations highlight the importance of quantifying and monitoring the dynamics and infections of malaria vectors during large-scale vector control interventions.

Physiology and development of the M and S molecular forms of Anopheles gambiae in Burkina Faso (West Africa)

In West Africa, M and S molecular forms of Anopheles gambiae sensu stricto (Diptera: Culicidae) Giles, frequently occur together, although with different population bionomics. The S form typically breeds in rain-dependant water collections and is present during the rainy season only whereas the M form can thrive all year long in areas with permanent breeding opportunities. In the present study, we explored physiological and developmental trade-offs at play in laboratory colonies and field populations of the M and S forms that originated from an area of sympatry in Burkina Faso, where M and S larvae exhibit such habitat segregation. In the laboratory, larvae of the M form developed slower than the S form (mean values 9.51 and 8.85 days, respectively, Wilcoxon's test, P < 0.001). Although wing length and dry weight at emergence showed large variations, M females were on average 8% heavier than S females of similar wing length. Higher nutritional reserves (proteins and lipids) in teneral adults explained part of this weight difference, reflecting a better ability of the M form to garner resources at the larval stage. Furthermore, a higher rate of ovarian maturation was observed in the M form after a single bloodmeal. The relevance of these findings for parasite transmission is discussed.

Evidence for population-specific positive selection on immune genes of Anopheles gambiae

Host-pathogen interactions can be powerful drivers of adaptive evolution, shaping the patterns of molecular variation at the genes involved. In this study, we sequenced alleles from 28 immune-related loci in wild samples of multiple genetic subpopulations of the African malaria mosquito Anopheles gambiae, obtaining unprecedented sample sizes and providing the first opportunity to contrast patterns of molecular evolution at immune-related loci in the recently discovered GOUNDRY population to those of the indoor-resting M and S molecular forms. In contrast to previous studies that focused on immune genes identified in laboratory studies, we centered our analysis on genes that fall within a quantitative trait locus associated with resistance to Plasmodium falciparum in natural populations of A. gambiae. Analyses of haplotypic and genetic diversity at these 28 loci revealed striking differences among populations in levels of genetic diversity and allele frequencies in coding sequence. Putative signals of positive selection were identified at 11 loci, but only one was shared by two subgroups of A. gambiae. Striking patterns of linkage disequilibrium were observed at several loci. We discuss these results with respect to ecological differences among these strata as well as potential implications for disease transmission.

Multiple insecticide resistance in Anopheles gambiae s.l. populations from Burkina Faso, West Africa

Malaria control programs are being jeopardized by the spread of insecticide resistance in mosquito vector populations. The situation in Burkina Faso is emblematic with Anopheles gambiae populations showing high levels of resistance to most available compounds. Although the frequency of insecticide target-site mutations including knockdown resistance (kdr) and insensitive acetylcholinesterase (Ace-1^R) alleles has been regularly monitored in the area, it is not known whether detoxifying enzymes contribute to the diversity of resistance phenotypes observed in the field. Here, we propose an update on the phenotypic diversity of insecticide resistance in An. gambiae populations sampled from 10 sites in Burkina Faso in 2010. Susceptibility to deltamethrin, permethrin, DDT, bendiocarb and fenithrotion was assessed. Test specimens (N = 30 per locality) were identified to species and molecular form and their genotype at the kdr and Ace-1 loci was determined. Detoxifying enzymes activities including non-specific esterases (NSEs), oxydases (cytochrome P450) and Glutathione S-Transferases (GSTs) were measured on single mosquitoes (N = 50) from each test locality and compared with the An. gambiae Kisumu susceptible reference strain. In all sites, mosquitoes demonstrated multiple resistance phenotypes, showing reduced mortality to several insecticidal compounds at the same time, although with considerable site-to-site variation. Both the kdr 1014L and Ace-1^R 119S resistant alleles were detected in the M and the S forms of An. gambiae, and were found together in specimens of the S form. Variation in detoxifying enzyme activities was observed within and between vector populations. Elevated levels of NSEs and GSTs were widespread, suggesting multiple resistance mechanisms segregate within An. gambiae populations from this country. By documenting the extent and diversity of insecticide resistance phenotypes and the putative combination of their underlying mechanisms in An. gambiae mosquitoes, our work prompts for new alternative strategies to be urgently developed for the control of major malaria vectors in Burkina Faso.

Use of a mixture statistical model in studying malaria vectors density

Vector control is a major step in the process of malaria control and elimination. This requires vector counts and appropriate statistical analyses of these counts. However, vector counts are often overdispersed. A non-parametric mixture of Poisson model (NPMP) is proposed to allow for overdispersion and better describe vector distribution. Mosquito collections using the Human Landing Catches as well as collection of environmental and climatic data were carried out from January to December 2009 in 28 villages in Southern Benin. A NPMP regression model with “village” as random effect is used to test statistical correlations between malaria vectors density and environmental and climatic factors. Furthermore, the villages were ranked using the latent classes derived from the NPMP model. Based on this classification of the villages, the impacts of four vector control strategies implemented in the villages were compared. Vector counts were highly variable and overdispersed with important proportion of zeros (75%). The NPMP model had a good aptitude to predict the observed values and showed that: i) proximity to freshwater body, market gardening, and high levels of rain were associated with high vector density; ii) water conveyance, cattle breeding, vegetation index were associated with low vector density. The 28 villages could then be ranked according to the mean vector number as estimated by the random part of the model after adjustment on all covariates. The NPMP model made it possible to describe the distribution of the vector across the study area. The villages were ranked according to the mean vector density after taking into account the most important covariates. This study demonstrates the necessity and possibility of adapting methods of vector counting and sampling to each setting.

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.

Population dynamics of Anopheles gambiae s.l. in Bobo-Dioulasso city: bionomics, infection rate and susceptibility to insecticides

BACKGROUND

Historical studies have indicated that An. gambiae s.s. is the predominant malaria vector species in Bobo-Dioulasso the second biggest city of Burkina Faso (West Africa). However, over the last decade, An. arabiensis appears to be replacing An. gambiae s.s. as the most prevalent malaria vector in this urban setting. To investigate this species transition in more detail the present study aims to provide an update on the malaria vector composition in Bobo-Dioulasso, and also the Plasmodium infection rates and susceptibility to insecticides of the local An. gambiae s.l. population.

METHODS

An entomological survey was carried out from May to December 2008 in Dioulassoba and Kodeni (central and peripheral districts respectively), which are representative of the main ecological features of the city. Sampling consisted of the collection of larval stages from water bodies, and adults by monthly indoor residual spraying (IRS) using aerosol insecticides. Insecticide susceptibility tests were performed using the WHO filter paper protocol on adults emerged from larvae. PCR was used to determine vector species and to identify resistance mechanisms (kdr and ace-1(R)). The Plasmodium infection rate was estimated by ELISA performed on female mosquitoes collected indoors by IRS.

RESULTS

An. arabiensis was found to be the major malaria vector in Bobo-Dioulasso, comprising 50 to 100% of the vector population. The sporozoite infection rate for An. arabiensis was higher than An. gambiae s.s. at both Dioulassoba and Kodeni. An. gambiae s.l. was resistant to DDT and cross-resistant to pyrethroids at the two sites with higher levels of resistance observed in An. gambiae s.s. than An. arabiensis. Resistance to 0.1% bendiocarb was observed in the An. gambiae s.s. S form but not the M form or in An. arabiensis. The L1014F kdr mutation was detected in the two molecular forms of An. gambiae s.s. at varying frequencies (0.45 to 0.92), but was not detected in An. arabiensis, suggesting that other mechanisms are involved in DDT resistance in this species. The ace-1(R) mutation was only detected in the S molecular form and was observed at the two sites at similar frequency (0.3).

CONCLUSIONS

Over the last ten years, An. arabiensis has become the major malaria vector in Bobo-Dioulasso city where it was formerly present only at low frequency. However, the ecological determinant that enhances the settlement of this species into urban and peri-urban areas of Bobo-Dioulasso remains to be clarified. The impact of the changing An. gambiae s.l. population in this region for vector control including resistance management strategies is discussed.

Evidence for X-linked introgression between molecular forms of Anopheles gambiae from Angola

The M and S molecular forms of the African malaria vector Anopheles gambiae (Diptera: Culicidae) are morphologically identical incipient species in which reproductive isolation is incomplete, enabling low-level gene flow between forms. In an attempt to find differences between the M and S forms, sequence variation was studied at loci along the X chromosome in adult female An. gambiae from Angola. A high proportion of M form specimens from Angola (79% of the 456 X chromosomes sampled) were found to contain a 16-bp insertion in intron 4 of the X-linked GPRCCK1 locus, relative to the AgamP3 release of the An. gambiae PEST genome sequence. The insertion was in Hardy-Weinberg equilibrium in Angolan M form populations. The same insertion was found in all S form specimens examined, regardless of where in Africa they were sampled, but was absent from a sample of M form specimens collected in Ghana, Bioko and Mali. In M form specimens from Angola, there was an association between alleles at the GPRCCK1 locus and those at a microsatellite locus, AGXH678, close to the centromere of the X chromosome, with significant linkage disequilibrium between loci separated by 0.472 Mbp (P < 0.033). We show that the insertion results from introgression from the S form into the M form, rather than from the retention of an ancestral character. Gene flow from the S to M form could allow genes of adaptive value to be transferred, including those conferring insecticide resistance and others influencing ecology and behaviour, and thus malaria transmission and control. We discuss factors that may have led to this introgression event.

Midgut microbiota of the malaria mosquito vector Anopheles gambiae and interactions with Plasmodium falciparum infection

The susceptibility of Anopheles mosquitoes to Plasmodium infections relies on complex interactions between the insect vector and the malaria parasite. A number of studies have shown that the mosquito innate immune responses play an important role in controlling the malaria infection and that the strength of parasite clearance is under genetic control, but little is known about the influence of environmental factors on the transmission success. We present here evidence that the composition of the vector gut microbiota is one of the major components that determine the outcome of mosquito infections. A. gambiae mosquitoes collected in natural breeding sites from Cameroon were experimentally challenged with a wild P. falciparum isolate, and their gut bacterial content was submitted for pyrosequencing analysis. The meta-taxogenomic approach revealed a broader richness of the midgut bacterial flora than previously described. Unexpectedly, the majority of bacterial species were found in only a small proportion of mosquitoes, and only 20 genera were shared by 80% of individuals. We show that observed differences in gut bacterial flora of adult mosquitoes is a result of breeding in distinct sites, suggesting that the native aquatic source where larvae were grown determines the composition of the midgut microbiota. Importantly, the abundance of Enterobacteriaceae in the mosquito midgut correlates significantly with the Plasmodium infection status. This striking relationship highlights the role of natural gut environment in parasite transmission. Deciphering microbe-pathogen interactions offers new perspectives to control disease transmission.

During their development in the mosquito vector, Plasmodium parasites undergo complex developmental steps and incur severe bottlenecks. The largest parasite losses occur in the mosquito midgut where robust immune responses are activated. Variability in P. falciparum infection levels indicates that parasite transmission is the result of complex interactions between vectors and parasites, which rely on both genetic and environmental factors. However, in contrast to genetically encoded factors, the role of environmental factors in parasite transmission has received little attention. In this study, we characterized the midgut microbiota of mosquitoes derived from diverse breeding sites using pyrosequencing. We show that the composition of the midgut microbiota in adult mosquitoes exhibits great variability, which is likely determined by bacterial richness of the larval habitats. When field mosquitoes were collected at late immature stages in natural breeding sites and the emerging females challenged with Plasmodium falciparum in the laboratory, significant correlation was observed between P. falciparum infection and the presence of Enterobacteriaceae in the mosquito midgut. Greater understanding of these malaria-bacteria interactions may lead to novel malaria control strategies.

No evidence for biased co-transmission of speciation islands in Anopheles gambiae

Genome-scale scans have revealed highly heterogeneous levels of divergence between closely related taxa in many systems. Generally, a small number of regions show high differentiation, with the rest of the genome showing no or only low levels of divergence. These patterns have been interpreted as evidence for ongoing speciation-with-gene-flow, with introgression homogenizing the whole genome except loci involved in reproductive isolation. However, as the number of selected loci increases, the probability of introgression at unselected loci decreases unless there is a transmission ratio distortion causing an over-representation of specific combinations of alleles. Here we examine the transmission of three 'speciation islands' that contain fixed differences between the M and S forms of the mosquito, Anopheles gambiae. We made reciprocal crosses between M and S parents and genotyped over 2000 F(2) individuals, developing a hierarchical likelihood model to identify specific genotypes that are under- or over-represented among the recombinant offspring. Though our overall results did not match the expected number of F(2) genotypes, we found no biased co-transmission among M or S alleles in the three islands. Our likelihood model did identify transmission ratio distortion at two of the three islands, but this distortion was small (approx. 3%) and in opposite directions for the two islands. We discuss how our results impinge on hypotheses of current gene flow between M and S and ongoing speciation-with-gene-flow in this system.

Differential Plasmodium falciparum infection of Anopheles gambiae s.s. molecular and chromosomal forms in Mali

Background

Anopheles gambiae sensu stricto (s.s.) is a primary vector of Plasmodium falciparum in sub-Saharan Africa. Although some physiological differences among molecular and chromosomal forms of this species have been demonstrated, the relative susceptibility to malaria parasite infection among them has not been unequivocally shown. The objective of this study was to investigate P. falciparum circumsporozoite protein infection (CSP) positivity among An. gambiae s.s. chromosomal and molecular forms.

Methods

Wild An. gambiae from two sites Kela (n = 464) and Sidarebougou (n = 266) in Mali were screened for the presence of P. falciparum CSP using an enzyme-linked immunosorbent assay (ELISA). Samples were then identified to molecular form using multiple PCR diagnostics (n = 713) and chromosomal form using chromosomal karyotyping (n = 419).

Results

Of 730 An. gambiae sensu lato (s.l.) mosquitoes, 89 (12.2%) were CSP ELISA positive. The percentage of positive mosquitoes varied by site: 52 (11.2%) in Kela and 37 (13.9%) in Sidarebougou. Eighty-seven of the positive mosquitoes were identified to molecular form and they consisted of nine Anopheles arabiensis (21.4%), 46 S (10.9%), 31 M (12.8%), and one MS hybrid (14.3%). Sixty of the positive mosquitoes were identified to chromosomal form and they consisted of five An. arabiensis (20.0%), 21 Savanna (15.1%), 21 Mopti (30.4%), 11 Bamako (9.2%), and two hybrids (20.0%).

Discussion

In this collection, the prevalence of P. falciparum infection in the M form was equivalent to infection in the S form (no molecular form differential infection). There was a significant differential infection by chromosomal form such that, P. falciparum infection was more prevalent in the Mopti chromosomal forms than in the Bamako or Savanna forms; the Mopti form was also the most underrepresented in the collection. Continued research on the differential P. falciparum infection of An. gambiae s.s. chromosomal and molecular forms may suggest that Plasmodium – An. gambiae interactions play a role in malaria transmission.

Adaptation to aridity in the malaria mosquito Anopheles gambiae: chromosomal inversion polymorphism and body size influence resistance to desiccation

Chromosomal inversions are thought to confer a selective advantage in alternative habitats by protecting co-adapted alleles from recombination. The frequencies of two inversions (2La and 2Rb) of the afro-tropical malaria mosquito Anopheles gambiae change gradually along geographical clines, increasing in frequency with degree of aridity. Such clines can result from gene flow and local selection acting upon alternative karyotypes along the cline, suggesting that these inversions may be associated with tolerance to xeric conditions. Since water loss represents a major challenge in xeric habitats, it can be supposed that genes inside these inversions are involved in water homeostasis. To test this hypothesis, we compared the desiccation resistance of alternative karyotypes from a colonised 2Rb/2La polymorphic population of A. gambiae from Cameroon. The strain included only the molecular form S, one of the genetic units marking incipient speciation in this taxon. Day-old mosquitoes of both sexes were assayed individually for time to death in a dry environment and the karyotype of each was determined post-mortem using molecular diagnostic assays for each inversion. In agreement with expectations based on their eco-geographical distribution, we found that 2La homokaryotypes survived significantly longer (1.3 hours) than the other karyotypes. However, there was weak support for the effect of 2Rb on desiccation resistance. Larger mosquitoes survived longer than smaller ones. Median survival of females was greater than males, but the effect of sex on desiccation resistance was weakly supported, indicating that differential survival was correlated to differences between sexes in average size. We found weak evidence for a heterotic effect of 2La karyotype on size in females. These results support the notion that genes located inside the 2La inversion are involved in water balance, contributing towards local adaptation of A. gambiae to xeric habitats, beyond the adaptive value conferred by a larger body size.

Behavioural responses of Anopheles gambiae sensu stricto M and S molecular form larvae to an aquatic predator in Burkina Faso

Background

Predation of aquatic immature stages has been identified as a major evolutionary force driving habitat segregation and niche partitioning in the malaria mosquito Anopheles gambiae sensu stricto in the humid savannahs of Burkina Faso, West Africa. Here, we explored behavioural responses to the presence of a predator in wild populations of the M and S molecular forms of An. gambiae that typically breed in permanent (e.g., rice field paddies) and temporary (e.g., road ruts) water collections.

Methods

Larvae used in these experiments were obtained from eggs laid by wild female An. gambiae collected from two localities in south-western Burkina Faso during the 2008 rainy season. Single larvae were observed in an experimental arena, and behavioural traits were recorded and quantified a) in the absence of a predator and b) in the presence of a widespread mosquito predator, the backswimmer Anisops jaczewskii. Differences in the proportion of time allocated to each behaviour were assessed using Principal Component Analysis and Multivariate Analysis of Variance.

Results

The behaviour of M and S form larvae was found to differ significantly; although both forms mainly foraged at the water surface, spending 60-90% of their time filtering water at the surface or along the wall of the container, M form larvae spent on average significantly more time browsing at the bottom of the container than S form larvae (4.5 vs. 1.3% of their overall time, respectively; P < 0.05). In the presence of a predator, larvae of both forms modified their behaviour, spending significantly more time resting along the container wall (P < 0.001). This change in behaviour was at least twice as great in the M form (from 38.6 to 66.6% of the time at the wall in the absence and presence of the predator, respectively) than in the S form (from 48.3 to 64.1%). Thrashing at the water surface exposed larvae to a significantly greater risk of predation by the notonectid (P < 0.01), whereas predation occurred significantly less often when larvae were at the container wall (P < 0.05) and might reflect predator vigilance.

Conclusions

Behavioural differences between larvae of the M and S form of An. gambiae in response to an acute predation risk is likely to be a reflection of different trade-offs between foraging and predator vigilance that might be of adaptive value in contrasting aquatic ecosystems. Future studies should explore the relevance of these findings under the wide range of natural settings where both forms co-exist in Africa.

Larval habitat segregation between the molecular forms of the mosquito Anopheles gambiae in a rice field area of Burkina Faso, West Africa

In West Africa, lineage splitting between the M and S molecular forms of the major Afro-tropical malaria mosquito, Anopheles gambiae (Diptera: Culicidae), is thought to be driven by ecological divergence, occurring mainly at the larval stage. Here, we present evidence for habitat segregation between the two molecular forms in and around irrigated rice fields located within the humid savannahs of western Burkina Faso. Longitudinal sampling of adult mosquitoes emerging from a range of breeding sites distributed along a transect extending from the heart of the rice field area into the surrounding savannah was conducted from June to November 2009. Analysis revealed that the two molecular forms and their sibling species Anopheles arabiensis are not randomly distributed in the area. A major ecological gradient was extracted in relation to the perimeter of the rice fields. The M form was associated with larger breeding sites mostly consisting of rice paddies, whereas the S form and An. arabiensis were found to depend upon temporary, rain-filled breeding sites. These results support hypotheses about larval habitat segregation and confirm the suggestion that the forms have different larval habitat requirements. Segregation appears to be clearly linked to anthropogenic permanent habitats and the community structure they support.

Spatially explicit analyses of anopheline mosquitoes indoor resting density: implications for malaria control

Background

The question of sampling and spatial aggregation of malaria vectors is central to vector control efforts and estimates of transmission. Spatial patterns of anopheline populations are complex because mosquitoes' habitats and behaviors are strongly heterogeneous. Analyses of spatially referenced counts provide a powerful approach to delineate complex distribution patterns, and contributions of these methods in the study and control of malaria vectors must be carefully evaluated.

Methodology/Principal Findings

We used correlograms, directional variograms, Local Indicators of Spatial Association (LISA) and the Spatial Analysis by Distance IndicEs (SADIE) to examine spatial patterns of Indoor Resting Densities (IRD) in two dominant malaria vectors sampled with a 5×5 km grid over a 2500 km² area in the forest domain of Cameroon. SADIE analyses revealed that the distribution of Anopheles gambiae was different from regular or random, whereas there was no evidence of spatial pattern in Anopheles funestus (Ia = 1.644, Pa<0.05 and Ia = 1.464, Pa>0.05, respectively). Correlograms and variograms showed significant spatial autocorrelations at small distance lags, and indicated the presence of large clusters of similar values of abundance in An. gambiae while An. funestus was characterized by smaller clusters. The examination of spatial patterns at a finer spatial scale with SADIE and LISA identified several patches of higher than average IRD (hot spots) and clusters of lower than average IRD (cold spots) for the two species. Significant changes occurred in the overall spatial pattern, spatial trends and clusters when IRDs were aggregated at the house level rather than the locality level. All spatial analyses unveiled scale-dependent patterns that could not be identified by traditional aggregation indices.

Conclusions/Significance

Our study illustrates the importance of spatial analyses in unraveling the complex spatial patterns of malaria vectors, and highlights the potential contributions of these methods in malaria control.

Plasmodium falciparum produce lower infection intensities in local versus foreign Anopheles gambiae populations

Both Plasmodium falciparum and Anopheles gambiae show great diversity in Africa, in their own genetic makeup and population dynamics. The genetics of the individual mosquito and parasite are known to play a role in determining the outcome of infection in the vector, but whether differences in infection phenotype vary between populations remains to be investigated. Here we established two A. gambiae s.s. M molecular form colonies from Cameroon and Burkina Faso, representing a local and a foreign population for each of the geographical sites. Experimental infections of both colonies were conducted in Cameroon and Burkina Faso using local wild P. falciparum, giving a sympatric and allopatric vector-parasite combination in each site. Infection phenotype was determined in terms of oocyst prevalence and intensity for at least nine infections for each vector-parasite combination. Sympatric infections were found to produce 25% fewer oocysts per midgut than allopatric infections, while prevalence was not affected by local/foreign interactions. The reduction in oocyst numbers in sympatric couples may be the result of evolutionary processes where the mosquito populations have locally adapted to their parasite populations. Future research on vector-parasite interactions must take into account the geographic scale of adaptation revealed here by conducting experiments in natural sympatric populations to give epidemiologically meaningful results.

Patterns of genomic differentiation between ecologically differentiated M and S forms of Anopheles gambiae in West and Central Africa

Anopheles gambiae M and S are thought to be undergoing ecological speciation by adapting to different larval habitats. Toward an improved understanding of the genetic determinants and evolutionary processes shaping their divergence, we used a 400,000 single-nucleotide polymorphism (SNP) genotyping array to characterize patterns of genomic differentiation between four geographically paired M and S population samples from West and Central Africa. In keeping with recent studies based on more limited genomic or geographic sampling, divergence was not confined to a few isolated "speciation islands." Divergence was both widespread across the genome and heterogeneous. Moreover, we find consistent patterns of genomic divergence across sampling sites and mutually exclusive clustering of M and S populations using genetic distances based on all 400,000 SNPs, implying that M and S are evolving collectively across the study area. Nevertheless, the clustering of local M and S populations using genetic distances based on SNPs from genomic regions of low differentiation is consistent with recent gene flow and introgression. To account for these data and reconcile apparent paradoxes in reported patterns of M-S genomic divergence and hybridization, we propose that extrinsic ecologically based postmating barriers vary in strength as environmental conditions fluctuate or change.

Ecological genomics of Anopheles gambiae along a latitudinal cline: a population-resequencing approach

The association between fitness-related phenotypic traits and an environmental gradient offers one of the best opportunities to study the interplay between natural selection and migration. In cases in which specific genetic variants also show such clinal patterns, it may be possible to uncover the mutations responsible for local adaptation. The malaria vector, Anopheles gambiae, is associated with a latitudinal cline in aridity in Cameroon; a large inversion on chromosome 2L of this mosquito shows large differences in frequency along this cline, with high frequencies of the inverted karyotype present in northern, more arid populations and an almost complete absence of the inverted arrangement in southern populations. Here we use a genome resequencing approach to investigate patterns of population divergence along the cline. By sequencing pools of individuals from both ends of the cline as well as in the center of the cline- where the inversion is present in intermediate frequency- we demonstrate almost complete panmixia across collinear parts of the genome and high levels of differentiation in inverted parts of the genome. Sequencing of separate pools of each inversion arrangement in the center of the cline reveals large amounts of gene flux (i.e., gene conversion and double crossovers) even within inverted regions, especially away from the inversion breakpoints. The interplay between natural selection, migration, and gene flux allows us to identify several candidate genes responsible for the match between inversion frequency and environmental variables. These results, coupled with similar conclusions from studies of clinal variation in Drosophila, point to a number of important biological functions associated with local environmental adaptation.

Oviposition habitat selection by Anopheles gambiae in response to chemical cues by Notonecta maculata

A number of mosquito species avoid predator-inhabited oviposition sites by detecting predator-released kairomones. In the laboratory, we found that when offered de-ionized water and de-ionized water conditioned with Notonecta maculata, gravid Anopheles gambiae females preferentially oviposited into the former. We then conducted further experiments using two chemical components found in Notonecta-conditioned water, chemically pure n-tricosane and/or n-heneicosane, that was previously shown to repel oviposition by Culiseta longiareolata. These hydrocarbons failed to deter oviposition by An. gambiae females. Thus, different mosquito species may rely on distinct chemical cues to avoid predators. Identification and chemical characterization of such kairomones could facilitate innovative, environmentally sound mosquito control.

Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization

The suggestion that genetic divergence can arise and/or be maintained in the face of gene flow has been contentious since first proposed. This controversy and a rarity of good examples have limited our understanding of this process. Partially reproductively isolated taxa have been highlighted as offering unique opportunities for identifying the mechanisms underlying divergence with gene flow. The African malaria vector, Anopheles gambiae s.s., is widely regarded as consisting of two sympatric forms, thought by many to represent incipient species, the M and S molecular forms. However, there has been much debate about the extent of reproductive isolation between M and S, with one view positing that divergence may have arisen and is being maintained in the presence of gene flow, and the other proposing a more advanced speciation process with little realized gene flow because of low hybrid fitness. These hypotheses have been difficult to address because hybrids are typically rare (<1%). Here, we assess samples from an area of high hybridization and demonstrate that hybrids are fit and responsible for extensive introgression. Nonetheless, we show that strong divergent selection at a subset of loci combined with highly asymmetric introgression has enabled M and S to remain genetically differentiated despite extensive gene flow. We propose that the extent of reproductive isolation between M and S varies across West Africa resulting in a 'geographic mosaic of reproductive isolation'; a finding which adds further complexity to our understanding of divergence in this taxon and which has considerable implications for transgenic control strategies.

Evolution of Anopheles gambiae in Relation to Humans and Malaria

The closely related and morphologically indistinguishable mosquito species in the Afrotropical Anopheles gambiae complex differ dramatically in their contribution to malaria transmission, ranging from major vectors through minor or locally important vectors and nonvectors. Radiation of the A. gambiae complex and ongoing diversification within its nominal species appears to be a product of recent and rapid adaptation to environmental heterogeneities, notably those of anthropogenic origin. Polytene chromosome and genomic analyses suggest that paracentric chromosomal inversions and possibly other low-recombination regions have played instrumental roles in this process by facilitating ecotypic differentiation both within and across semipermeable species boundaries. Forthcoming complete genome sequences from several members of the A. gambiae complex will provide powerful tools to accelerate ongoing investigation of how genetic diversification of populations and species has shaped behavioral and physiological traits, such as vector competence, that bear on vectorial importance.

Morphological differentiation may mediate mate-choice between incipient species of Anopheles gambiae s.s

The M and S molecular forms of Anopheles gambiae s.s. have been considered incipient species for more than ten years, yet the mechanism underlying assortative mating of these incipient species has remained elusive. The discovery of the importance of harmonic convergence of wing beat frequency in mosquito mating and its relation to wing size have laid the foundation for exploring phenotypic divergence in wing size of wild populations of the two forms. In this study, wings from field collected mosquitoes were measured for wing length and wing width from two parts of the sympatric distribution, which differ with respect to the strength of assortative mating. In Mali, where assortative mating is strong, as evidenced by low rates of hybridization, mean wing lengths and wing widths were significantly larger than those from Guinea-Bissau. In addition, mean wing widths in Mali were significantly different between molecular forms. In Guinea-Bissau, assortative mating appears comparatively reduced and wing lengths and widths did not differ significantly between molecular forms. The data presented in this study support the hypothesis that wing beat frequency may mediate assortative mating in the incipient species of A. gambiae and represent the first documentation of a morphological difference between the M and S molecular forms.

Improving the population genetics toolbox for the study of the African malaria vector Anopheles nili: microsatellite mapping to chromosomes

BACKGROUND

Anopheles nili is a major vector of malaria in the humid savannas and forested areas of sub-Saharan Africa. Understanding the population genetic structure and evolutionary dynamics of this species is important for the development of an adequate and targeted malaria control strategy in Africa. Chromosomal inversions and microsatellite markers are commonly used for studying the population structure of malaria mosquitoes. Physical mapping of these markers onto the chromosomes further improves the toolbox, and allows inference on the demographic and evolutionary history of the target species.

RESULTS

Availability of polytene chromosomes allowed us to develop a map of microsatellite markers and to study polymorphism of chromosomal inversions. Nine microsatellite markers were mapped to unique locations on all five chromosomal arms of An. nili using fluorescent in situ hybridization (FISH). Probes were obtained from 300-483 bp-long inserts of plasmid clones and from 506-559 bp-long fragments amplified with primers designed using the An. nili genome assembly generated on an Illumina platform. Two additional loci were assigned to specific chromosome arms of An. nili based on in silico sequence similarity and chromosome synteny with Anopheles gambiae. Three microsatellites were mapped inside or in the vicinity of the polymorphic chromosomal inversions 2Rb and 2Rc. A statistically significant departure from Hardy-Weinberg equilibrium, due to a deficit in heterozygotes at the 2Rb inversion, and highly significant linkage disequilibrium between the two inversions, were detected in natural An. nili populations collected from Burkina Faso.

CONCLUSIONS

Our study demonstrated that next-generation sequencing can be used to improve FISH for microsatellite mapping in species with no reference genome sequence. Physical mapping of microsatellite markers in An. nili showed that their cytological locations spanned the entire five-arm complement, allowing genome-wide inferences. The knowledge about polymorphic inversions and chromosomal locations of microsatellite markers has been useful for explaining differences in genetic variability across loci and significant differentiation observed among natural populations of An. nili.

Variation in metabolic rate of Anopheles gambiae and A. arabiensis in a Sahelian village

In the Sahel, the Anopheles gambiae complex consists of Anopheles arabiensis and the M and S molecular forms of A. gambiae sensu stricto. However, the composition of these malaria vectors varies spatially and temporally throughout the region and is thought to be linked to environmental factors such as rainfall, larval site characteristics and duration of the dry season. To examine possible physiological divergence between these taxa, we measured metabolic rates of mosquitoes during the wet season in a Sahelian village in Mali. To our knowledge, this study provides the first measurements of metabolic rates of A. gambiae and A. arabiensis in the field. The mean metabolic rate of A. arabiensis was higher than that of M-form A. gambiae when accounting for the effects of female gonotrophic status, temperature and flight activity. However, after accounting for their difference in body size, no significant difference in metabolic rate was found between these two species (whilst all other factors were found to be significant). Thus, body size may be a key character that has diverged in response to ecological differences between these two species. Alternatively, these species may display additional differences in metabolic rate only during the dry season. Overall, our results indicate that changes in behavior and feeding activity provide an effective mechanism for mosquitoes to reduce their metabolic rate, and provide insight into the possible strategies employed by aestivating individuals during the dry season. We hypothesize that female mosquitoes switch to sugar feeding while in dormancy because of elevated metabolism associated with blood digestion.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

Comparative susceptibility to Plasmodium falciparum of the molecular forms M and S of Anopheles gambiae and Anopheles arabiensis

Background

The different taxa belonging to Anopheles gambiae complex display phenotypic differences that may impact their contribution to malaria transmission. More specifically, their susceptibility to infection, resulting from a co-evolution between parasite and vector, might be different. The aim of this study was to compare the susceptibility of M and S molecular forms of Anopheles gambiae and Anopheles arabiensis to infection by Plasmodium falciparum.

Methods

F3 progenies of Anopheles gambiae s.l. collected in Senegal were infected, using direct membrane feeding, with P. falciparum gametocyte-containing blood sampled on volunteer patients. The presence of oocysts was determined by light microscopy after 7 days, and the presence of sporozoite by ELISA after 14 days. Mosquito species and molecular forms were identified by PCR.

Results

The oocyst rate was significantly higher in the molecular S form (79.07%) than in the M form (57.81%, Fisher's exact test p < 0.001) and in Anopheles arabiensis (55.38%, Fisher's exact test vs. S group p < 0.001). Mean ± s.e.m. number of oocyst was greater in the An. gambiae S form (1.72 ± 0.26) than in the An. gambiae M form (0.64 ± 0.04, p < 0.0001) and in the An. arabiensis group (0.58 ± 0.04, vs. S group, p < 0.0001). Sporozoite rate was also higher in the molecular form S (83.52%) than in form M (50.98%, Fisher's exact test p < 0.001) and Anopheles arabiensis 50.85%, Fisher's exact test vs. S group p < 0.001).

Conclusion

Infected in the same experimental conditions, the molecular form S of An. gambiae is more susceptible to infection by P. falciparum than the molecular form M of An. gambiae and An. arabiensis.

Dynamics of insecticide resistance in malaria vectors in Benin: first evidence of the presence of L1014S kdr mutation in Anopheles gambiae from West Africa

Background

Insecticide resistance monitoring is essential to help national programmers to implement more effective and sustainable malaria control strategies in endemic countries. This study reported the spatial and seasonal variations of insecticide resistance in malaria vectors in Benin, West Africa.

Methods

Anopheles gambiae s.l populations were collected from October 2008 to June 2010 in four sites selected on the basis of different use of insecticides and environment. WHO susceptibility tests were carried out to detect resistance to DDT, fenitrothion, bendiocarb, permethrin and deltamethrin. The synergist piperonyl butoxide was used to assess the role of non-target site mechanisms in pyrethroid resistance. Anopheles gambiae mosquitoes were identified to species and to molecular M and S forms using PCR techniques. Molecular and biochemical assays were carried out to determine kdr and Ace.1^R allelic frequencies and activity of the detoxification enzymes.

Results

Throughout the surveys very high levels of mortality to bendiocarb and fenitrothion were observed in An. gambiae s.l. populations. However, high frequencies of resistance to DDT and pyrethroids were seen in both M and S form of An. gambiae s.s. and Anopheles arabiensis. PBO increased the toxicity of permethrin and restored almost full susceptibility to deltamethrin. Anopheles gambiae s.l. mosquitoes from Cotonou and Malanville showed higher oxidase activity compared to the Kisumu susceptible strain in 2009, whereas the esterase activity was higher in the mosquitoes from Bohicon in both 2008 and 2009. A high frequency of 1014F kdr allele was initially showed in An. gambiae from Cotonou and Tori-Bossito whereas it increased in mosquitoes from Bohicon and Malanville during the second year. For the first time the L1014S kdr mutation was found in An. arabiensis in Benin. The ace.1^R mutation was almost absent in An. gambiae s.l.

Conclusion

Pyrethroid and DDT resistance is widespread in malaria vector in Benin and both metabolic and target site resistance are implicated. Resistance was not correlated with a change of malaria species and/or molecular forms. The 1014S kdr allele was first identified in wild population of An. arabiensis hence confirming the expansion of pyrethroid resistance alleles in Africa.

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.

Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiaemolecular form identification

BACKGROUND

Anopheles gambiae M and S molecular forms, the major malaria vectors in the Afro-tropical region, are ongoing a process of ecological diversification and adaptive lineage splitting, which is affecting malaria transmission and vector control strategies in West Africa. These two incipient species are defined on the basis of single nucleotide differences in the IGS and ITS regions of multicopy rDNA located on the X-chromosome. A number of PCR and PCR-RFLP approaches based on form-specific SNPs in the IGS region are used for M and S identification. Moreover, a PCR-method to detect the M-specific insertion of a short interspersed transposable element (SINE200) has recently been introduced as an alternative identification approach. However, a large-scale comparative analysis of four widely used PCR or PCR-RFLP genotyping methods for M and S identification was never carried out to evaluate whether they could be used interchangeably, as commonly assumed.

RESULTS

The genotyping of more than 400 A. gambiae specimens from nine African countries, and the sequencing of the IGS-amplicon of 115 of them, highlighted discrepancies among results obtained by the different approaches due to different kinds of biases, which may result in an overestimation of MS putative hybrids, as follows: i) incorrect match of M and S specific primers used in the allele specific-PCR approach; ii) presence of polymorphisms in the recognition sequence of restriction enzymes used in the PCR-RFLP approaches; iii) incomplete cleavage during the restriction reactions; iv) presence of different copy numbers of M and S-specific IGS-arrays in single individuals in areas of secondary contact between the two forms.

CONCLUSIONS

The results reveal that the PCR and PCR-RFLP approaches most commonly utilized to identify A. gambiae M and S forms are not fully interchangeable as usually assumed, and highlight limits of the actual definition of the two molecular forms, which might not fully correspond to the two A. gambiae incipient species in their entire geographical range. These limits are discussed and operational suggestions on the choice of the most convenient method for large-scale M- and S-form identification are provided, also taking into consideration technical aspects related to the epidemiological characteristics of different study areas.

Divergent transcriptional response to thermal stress by Anopheles gambiae larvae carrying alternative arrangements of inversion 2La

The African malaria mosquito Anopheles gambiae is polymorphic for chromosomal inversion 2La, whose frequency strongly correlates with degree of aridity across environmental gradients. Recent physiological studies have associated 2La with resistance to desiccation in adults and thermal stress in larvae, consistent with its proposed role in aridity tolerance. However, the genetic basis of these traits remains unknown. To identify genes that could be involved in the differential response to thermal stress, we compared global gene expression profiles of heat-hardened 2La or 2L+(a) larvae at three time points, for up to eight hours following exposure to the heat stress. Treatment and control time series, replicated four times, revealed a common and massive induction of a core set of heat-shock genes regardless of 2La orientation. However, clear differences between the 2La and 2L+(a) arrangements emerged at the earliest (0.25 h) time point, in the intensity and nature of the stress response. Overall, 2La was associated with the more aggressive response: larger numbers of genes were heat responsive and up-regulated. Transcriptionally induced genes were enriched for functions related to ubiquitin-proteasomal degradation, chaperoning and energy metabolism. The more muted transcriptional response of 2L+(a) was largely repressive, including genes involved in proteolysis and energy metabolism. These results may help explain the maintenance of the 2La inversion polymorphism in An. gambiae, as the survival benefits offered by high thermal sensitivity in harsh climates could be offset by the metabolic costs of such a drastic response in more equable climates.

The "far-west" of Anopheles gambiae molecular forms

The main Afrotropical malaria vector, Anopheles gambiae sensu stricto, is undergoing a process of sympatric ecological diversification leading to at least two incipient species (the M and S molecular forms) showing heterogeneous levels of divergence across the genome. The physically unlinked centromeric regions on all three chromosomes of these closely related taxa contain fixed nucleotide differences which have been found in nearly complete linkage disequilibrium in geographic areas of no or low M-S hybridization. Assays diagnostic for SNP and structural differences between M and S forms in the three centromeric regions were applied in samples from the western extreme of their range of sympatry, the only area where high frequencies of putative M/S hybrids have been reported. The results reveal a level of admixture not observed in the rest of the range. In particular, we found: i) heterozygous genotypes at each marker, although at frequencies lower than expected under panmixia; ii) virtually all possible genotypic combinations between markers on different chromosomes, although genetic association was nevertheless detected; iii) discordant M and S genotypes at two X-linked markers near the centromere, suggestive of introgression and inter-locus recombination. These results could be indicative either of a secondary contact zone between M and S, or of the maintenance of ancestral polymorphisms. This issue and the perspectives opened by these results in the study of the M and S incipient speciation process are discussed.

A cryptic subgroup of Anopheles gambiae is highly susceptible to human malaria parasites

Population subgroups of the African malaria vector Anopheles gambiae have not been comprehensively characterized owing to the lack of unbiased sampling methods. In the arid savanna zone of West Africa, where potential oviposition sites are scarce, widespread collection from larval pools in the peridomestic human habitat yielded a comprehensive genetic survey of local A. gambiae population subgroups, independent of adult resting behavior and ecological preference. A previously unknown subgroup of exophilic A. gambiae is sympatric with the known endophilic A. gambiae in this region. The exophilic subgroup is abundant, lacks differentiation into M and S molecular forms, and is highly susceptible to infection with wild Plasmodium falciparum. These findings might have implications for the epidemiology of malaria transmission and control.

Chromosomal inversions, natural selection and adaptation in the malaria vector Anopheles funestus

Chromosomal polymorphisms, such as inversions, are presumably involved in the rapid adaptation of populations to local environmental conditions. Reduced recombination between alternative arrangements in heterozygotes may protect sets of locally adapted genes, promoting ecological divergence and potentially leading to reproductive isolation and speciation. Through a comparative analysis of chromosomal inversions and microsatellite marker polymorphisms, we hereby present biological evidence that strengthens this view in the mosquito Anopheles funestus s.s, one of the most important and widespread malaria vectors in Africa. Specimens were collected across a wide range of geographical, ecological, and climatic conditions in Cameroon. We observed a sharp contrast between population structure measured at neutral microsatellite markers and at chromosomal inversions. Microsatellite data detected only a weak signal for population structuring among geographical zones (F(ST) < 0.013, P < 0.01). By contrast, strong differentiation among ecological zones was revealed by chromosomal inversions (F(ST) > 0.190, P < 0.01). Using standardized estimates of F(ST), we show that inversions behave at odds with neutral expectations strongly suggesting a role of environmental selection in shaping their distribution. We further demonstrate through canonical correspondence analysis that heterogeneity in eco-geographical variables measured at specimen sampling sites explained 89% of chromosomal variance in A. funestus. These results are in agreement with a role of chromosomal inversions in ecotypic adaptation in this species. We argue that this widespread mosquito represents an interesting model system for the study of chromosomal speciation mechanisms and should provide ample opportunity for comparative studies on the evolution of reproductive isolation and speciation in major human malaria vectors.

Adaptive divergence between incipient species of Anopheles gambiae increases resistance to Plasmodium

The African malaria mosquito Anopheles gambiae is diversifying into ecotypes known as M and S forms. This process is thought to be promoted by adaptation to different larval habitats, but its genetic underpinnings remain elusive. To identify candidate targets of divergent natural selection in M and S, we performed genomewide scanning in paired population samples from Mali, followed by resequencing and genotyping from five locations in West, Central, and East Africa. Genome scans revealed a significant peak of M-S divergence on chromosome 3L, overlapping five known or suspected immune response genes. Resequencing implicated a selective target at or near the TEP1 gene, whose complement C3-like product has antiparasitic and antibacterial activity. Sequencing and allele-specific genotyping showed that an allelic variant of TEP1 has been swept to fixation in M samples from Mali and Burkina Faso and is spreading into neighboring Ghana, but is absent from M sampled in Cameroon, and from all sampled S populations. Sequence comparison demonstrates that this allele is related to, but distinct from, TEP1 alleles of known resistance phenotype. Experimental parasite infections of advanced mosquito intercrosses demonstrated a strong association between this TEP1 variant and resistance to both rodent malaria and the native human malaria parasite Plasmodium falciparum. Although malaria parasites may not be direct agents of pathogen-mediated selection at TEP1 in nature--where larvae may be the more vulnerable life stage--the process of adaptive divergence between M and S has potential consequences for malaria transmission.

The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis

Background

This is the second in a series of three articles documenting the geographical distribution of 41 dominant vector species (DVS) of human malaria. The first paper addressed the DVS of the Americas and the third will consider those of the Asian Pacific Region. Here, the DVS of Africa, Europe and the Middle East are discussed. The continent of Africa experiences the bulk of the global malaria burden due in part to the presence of the An. gambiae complex. Anopheles gambiae is one of four DVS within the An. gambiae complex, the others being An. arabiensis and the coastal An. merus and An. melas. There are a further three, highly anthropophilic DVS in Africa, An. funestus, An. moucheti and An. nili. Conversely, across Europe and the Middle East, malaria transmission is low and frequently absent, despite the presence of six DVS. To help control malaria in Africa and the Middle East, or to identify the risk of its re-emergence in Europe, the contemporary distribution and bionomics of the relevant DVS are needed.

Results

A contemporary database of occurrence data, compiled from the formal literature and other relevant resources, resulted in the collation of information for seven DVS from 44 countries in Africa containing 4234 geo-referenced, independent sites. In Europe and the Middle East, six DVS were identified from 2784 geo-referenced sites across 49 countries. These occurrence data were combined with expert opinion ranges and a suite of environmental and climatic variables of relevance to anopheline ecology to produce predictive distribution maps using the Boosted Regression Tree (BRT) method.

Conclusions

The predicted geographic extent for the following DVS (or species/suspected species complex*) is provided for Africa: Anopheles (Cellia) arabiensis, An. (Cel.) funestus*, An. (Cel.) gambiae, An. (Cel.) melas, An. (Cel.) merus, An. (Cel.) moucheti and An. (Cel.) nili*, and in the European and Middle Eastern Region: An. (Anopheles) atroparvus, An. (Ano.) labranchiae, An. (Ano.) messeae, An. (Ano.) sacharovi, An. (Cel.) sergentii and An. (Cel.) superpictus*. These maps are presented alongside a bionomics summary for each species relevant to its control.

Humming in tune: sex and species recognition by mosquitoes on the wing

Mosquitoes are more sensitive to sound than any other insect due to the remarkable properties of their antennae and Johnston's organ at the base of each antenna. Male mosquitoes detect and locate female mosquitoes by hearing the female's flight tone, but until recently we had no idea that females also respond to male flight tones. Our investigation of a novel mechanism of sex recognition in Toxorhynchites brevipalpis revealed that male and female mosquitoes actively respond to the flight tones of other flying mosquitoes by altering their own wing-beat frequencies. Male-female pairs converge on a shared harmonic of their respective fundamental flight tones, whereas same sex pairs diverge. Most frequency matching occurs at frequencies beyond the detection range of the Johnston's organ but within the range of mechanical responsiveness of the antennae. We have shown that this is possible because the Johnston's organ is tuned to, and able to detect difference tones in, the harmonics of antennal vibrations which are generated by the combined input of flight tones from both mosquitoes. Acoustic distortion in hearing organs exists usually as an interesting epiphenomenon. Mosquitoes, however, appear to use it as a sensory cue that enables male-female pairs to communicate through a signal that depends on auditory interactions between them. Frequency matching may also provide a means of species recognition. Morphologically identical but reproductively isolated molecular forms of Anopheles gambiae fly in the same mating swarms, but rarely hybridize. Extended frequency matching occurs almost exclusively between males and females of the same molecular form, suggesting that this behavior is associated with observed assortative mating.

Breakpoint structure of the Anopheles gambiae 2Rb chromosomal inversion

Background

Alternative arrangements of chromosome 2 inversions in Anopheles gambiae are important sources of population structure, and are associated with adaptation to environmental heterogeneity. The forces responsible for their origin and maintenance are incompletely understood. Molecular characterization of inversion breakpoints provides insight into how they arose, and provides the basis for development of molecular karyotyping methods useful in future studies.

Methods

Sequence comparison of regions near the cytological breakpoints of 2Rb allowed the molecular delineation of breakpoint boundaries. Comparisons were made between the standard 2R+^b arrangement in the An. gambiae PEST reference genome and the inverted 2Rb arrangements in the An. gambiae M and S genome assemblies. Sequence differences between alternative 2Rb arrangements were exploited in the design of a PCR diagnostic assay, which was evaluated against the known chromosomal banding pattern of laboratory colonies and field-collected samples from Mali and Cameroon.

Results

The breakpoints of the 7.55 Mb 2Rb inversion are flanked by extensive runs of the same short (72 bp) tandemly organized sequence, which was likely responsible for chromosomal breakage and rearrangement. Application of the molecular diagnostic assay suggested that 2Rb has a single common origin in An. gambiae and its sibling species, Anopheles arabiensis, and also that the standard arrangement (2R+^b) may have arisen twice through breakpoint reuse. The molecular diagnostic was reliable when applied to laboratory colonies, but its accuracy was lower in natural populations.

Conclusions

The complex repetitive sequence flanking the 2Rb breakpoint region may be prone to structural and sequence-level instability. The 2Rb molecular diagnostic has immediate application in studies based on laboratory colonies, but its usefulness in natural populations awaits development of complementary molecular tools.

Widespread divergence between incipient Anopheles gambiae species revealed by whole genome sequences

The Afrotropical mosquito Anopheles gambiae sensu stricto, a major vector of malaria, is currently undergoing speciation into the M and S molecular forms. These forms have diverged in larval ecology and reproductive behavior through unknown genetic mechanisms, despite considerable levels of hybridization. Previous genome-wide scans using gene-based microarrays uncovered divergence between M and S that was largely confined to gene-poor pericentromeric regions, prompting a speciation-with-ongoing-gene-flow model that implicated only about 3% of the genome near centromeres in the speciation process. Here, based on the complete M and S genome sequences, we report widespread and heterogeneous genomic divergence inconsistent with appreciable levels of interform gene flow, suggesting a more advanced speciation process and greater challenges to identify genes critical to initiating that process.

Association mapping of insecticide resistance in wild Anopheles gambiae populations: major variants identified in a low-linkage disequilbrium genome

Background

Association studies are a promising way to uncover the genetic basis of complex traits in wild populations. Data on population stratification, linkage disequilibrium and distribution of variant effect-sizes for different trait-types are required to predict study success but are lacking for most taxa. We quantified and investigated the impacts of these key variables in a large-scale association study of a strongly selected trait of medical importance: pyrethroid resistance in the African malaria vector Anopheles gambiae.

Methodology/Principal Findings

We genotyped ≈1500 resistance-phenotyped wild mosquitoes from Ghana and Cameroon using a 1536-SNP array enriched for candidate insecticide resistance gene SNPs. Three factors greatly impacted study power. (1) Population stratification, which was attributable to co-occurrence of molecular forms (M and S), and cryptic within-form stratification necessitating both a partitioned analysis and genomic control. (2) All SNPs of substantial effect (odds ratio, OR>2) were rare (minor allele frequency, MAF<0.05). (3) Linkage disequilibrium (LD) was very low throughout most of the genome. Nevertheless, locally high LD, consistent with a recent selective sweep, and uniformly high ORs in each subsample facilitated significant direct and indirect detection of the known insecticide target site mutation kdr L1014F (OR≈6; P<10⁻⁶), but with resistance level modified by local haplotypic background.

Conclusion

Primarily as a result of very low LD in wild A. Gambiae, LD-based association mapping is challenging, but is feasible at least for major effect variants, especially where LD is enhanced by selective sweeps. Such variants will be of greatest importance for predictive diagnostic screening.

Genome-wide analysis of transcriptomic divergence between laboratory colony and field Anopheles gambiae mosquitoes of the M and S molecular forms

Our knowledge of Anopheles gambiae molecular biology has mainly been based on studies using inbred laboratory strains. Differences in the environmental exposure of these and natural field mosquitoes have inevitably led to physiological divergences. We have used global transcript abundance analyses to probe into this divergence, and identified transcript abundance patterns of genes that provide insight on specific adaptations of caged and field mosquitoes. We also compared the gene transcript abundance profiles of field mosquitoes belonging to the two morphologically indistinguishable but reproductively isolated sympatric molecular forms, M and S, from two different locations in the Yaoundé area of Cameroon. This analysis suggested that environmental exposure has a greater influence on the transcriptome than does the mosquito's molecular form-specific genetic background.

Mitotic and polytene chromosomes analysis of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

The Oriental fruit fly, Batrocera dorsalis s.s. (Hendel) is one of the most destructive agricultural pests, belonging to a large group of difficult to distinguish morphologically species, referred as the B. dorsalis complex. We report here a cytogenetic analysis of two laboratory strains of the species and provide a photographic polytene chromosome map from larval salivary glands. The mitotic complement consists of six chromosome pairs including a heteromorphic sex (XX/XY) chromosome pair. Analysis of the polytene complement has shown a total of five polytene chromosomes (10 polytene arms) that correspond to the five autosomes. The most important landmarks of each polytene chromosome and characteristic asynapsis at a specific chromosomal region are presented and discussed. Chromosomal homology between B. dorsalis and Ceratitis capitata has been determined by comparing chromosome banding patterns. The detection of chromosome inversions in both B. dorsalis strains is shown and discussed. Our results show that the polytene maps presented here are suitable for cytogenetic analysis of this species and can be used for comparative studies among species of the Tephritidae family. They also provide a diagnostic tool that could accelerate species identification within the B. dorsalis complex and could shed light on the ongoing speciation in this complex. Polytene chromosome maps can facilitate the development of biological control methods and support the genome mapping project of the species that is currently in progress.

Major effect genes or loose confederations? The development of insecticide resistance in the malaria vector Anopheles gambiae

Insecticide use in public health and agriculture presents a dramatic adaptive challenge to target and non-target insect populations. The rapid development of genetically modulated resistance to insecticides is postulated to develop in two distinct ways: By selection for single major effect genes or by selection for loose confederations in which several factors, not normally associated with each other, inadvertently combine their effects to produce resistance phenotypes. Insecticide resistance is a common occurrence and has been intensively studied in the major malaria vector Anopheles gambiae, providing a useful model for examining how insecticide resistance develops and what pleiotropic effects are likely to emerge as a consequence of resistance. As malaria vector control becomes increasingly reliant on successfully managing insecticide resistance, the characterisation of resistance mechanisms and their pleiotropic effects becomes increasingly important.

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.