Inversion polymorphism and incipient speciation in Anopheles gambiae s.str. in The Gambia, West Africa

In Silico Karyotyping of Chromosomally Polymorphic Malaria Mosquitoes in the Anopheles gambiae Complex

Chromosomal inversion polymorphisms play an important role in adaptation to environmental heterogeneities. For mosquito species in the Anopheles gambiae complex that are significant vectors of human malaria, paracentric inversion polymorphisms are abundant and are associated with ecologically and epidemiologically important phenotypes. Improved understanding of these traits relies on determining mosquito karyotype, which currently depends upon laborious cytogenetic methods whose application is limited both by the requirement for specialized expertise and for properly preserved adult females at specific gonotrophic stages. To overcome this limitation, we developed sets of tag single nucleotide polymorphisms (SNPs) inside inversions whose biallelic genotype is strongly correlated with inversion genotype. We leveraged 1,347 fully sequenced An. gambiae and Anopheles coluzzii genomes in the Ag1000G database of natural variation. Beginning with principal components analysis (PCA) of population samples, applied to windows of the genome containing individual chromosomal rearrangements, we classified samples into three inversion genotypes, distinguishing homozygous inverted and homozygous uninverted groups by inclusion of the small subset of specimens in Ag1000G that are associated with cytogenetic metadata. We then assessed the correlation between candidate tag SNP genotypes and PCA-based inversion genotypes in our training sets, selecting those candidates with >80% agreement. Our initial tests both in held-back validation samples from Ag1000G and in data independent of Ag1000G suggest that when used for in silico inversion genotyping of sequenced mosquitoes, these tags perform better than traditional cytogenetics, even for specimens where only a small subset of the tag SNPs can be successfully ascertained.

Entomological impact of indoor residual spraying with pirimiphos-methyl: a pilot study in an area of low malaria transmission in Senegal

Background

Scaling-up of effective anti-malarial control strategies in Central-West region of Senegal has resulted in the sharp decline in malaria prevalence in this area. However, despite these strategies, residual malaria transmission has been observed in some villages (hot spots). The objective of this study was to assess the impact of indoor residual spraying (IRS) with pirimiphos-methyl on malaria transmission in hot spot areas.

Methods

The malaria vector population dynamics were monitored in each of the six selected villages (4 of which used IRS, 2 were unsprayed control areas) using overnight human landing catches (HLC) and pyrethrum spray catches (PSC). The host source of blood meals from freshly fed females collected using PSC was identified using the direct ELISA method. Females caught through HLC were tested by ELISA for the detection of Plasmodium falciparum circumsporozoite protein and Anopheles gambiae complex was identified using PCR.

Results

Preliminary data shown that the densities of Anopheles populations were significantly lower in the sprayed areas (179/702) compared to the control. Overall, malaria transmission risk was 14 times lower in the intervention zone (0.94) compared to the control zone (12.7). In the control areas, three Anopheles species belonging to the Gambiae complex (Anopheles arabiensis, Anopheles coluzzii and Anopheles melas) maintained the transmission, while only An. coluzzii was infective in the sprayed areas.

Conclusion

The preliminary data from this pilot study showed that IRS with the CS formulation of pirimiphos-methyl is likely very effective in reducing malaria transmission risk. However, additional studies including further longitudinal entomological surveys as well as ecological and ethological and genetical characterization of vectors species and their populations are needed to better characterize the entomological impact of indoor residual spraying with pirimiphos-methyl in the residual transmission areas of Senegal.

Massive introgression drives species radiation at the range limit of Anopheles gambiae

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae.

The contribution of dietary restriction to extended longevity in the malaria vector Anopheles coluzzii

BACKGROUND

Variation in longevity has long been of interest in vector biology because of its implication in disease transmission through vectorial capacity. Recent studies suggest that Anopheles coluzzii adults persist during the ~7 month dry season via aestivation. Recently there has been a growing body of evidence linking dietary restriction and low ratio of dietary protein to carbohydrate with extended longevity of animals. Here, we evaluated the effects of dietary restriction and the protein : carbohydrate ratio on longevity of An. coluzzii.

RESULTS

In our experiment, we combined dietary regimes with temperature and relative humidity to assess their effects on An. coluzzii longevity, in an attempt to simulate aestivation under laboratory conditions. Our results showed significant effects of both the physical and the dietary variables on longevity, but that diet regimen had a considerably greater effect than those of the physical conditions. Higher temperature and lower humidity reduced longevity. At 22 °C dietary protein (blood) shortened longevity when sugar was not restricted (RH = 85%), but extended longevity when sugar was restricted (RH = 50%).

CONCLUSIONS

Dietary restriction extended longevity in accord with predictions, but protein : carbohydrate ratio had a negligible effect. We identified conditions that significantly extend longevity in malaria vectors, however, the extent of increase in longevity was insufficient to simulate aestivation.

Characterisation of Species and Diversity of Anopheles gambiae Keele Colony

Anopheles gambiae sensu stricto was recently reclassified as two species, An. coluzzii and An. gambiae s.s., in wild-caught mosquitoes, on the basis of the molecular form, denoted M or S, of a marker on the X chromosome. The An. gambiae Keele line is an outbred laboratory colony strain that was developed around 12 years ago by crosses between mosquitoes from 4 existing An. gambiae colonies. Laboratory colonies of mosquitoes often have limited genetic diversity because of small starting populations (founder effect) and subsequent fluctuations in colony size. Here we describe the characterisation of the chromosomal form(s) present in the Keele line, and investigate the diversity present in the colony using microsatellite markers on chromosome 3. We also characterise the large 2La inversion on chromosome 2. The results indicate that only the M-form of the chromosome X marker is present in the Keele colony, which was unexpected given that 3 of the 4 parent colonies were probably S-form. Levels of diversity were relatively high, as indicated by a mean number of microsatellite alleles of 6.25 across 4 microsatellites, in at least 25 mosquitoes. Both karyotypes of the inversion on chromosome 2 (2La/2L+^a) were found to be present at approximately equal proportions. The Keele colony has a mixed M- and S-form origin, and in common with the PEST strain, we propose continuing to denote it as an An. gambiae s.s. line.

Decline in frequency of the 2La chromosomal inversion in Anopheles gambiae (s.s.) in Western Kenya: correlation with increase in ownership of insecticide-treated bed nets

Background

The 2La chromosomal inversion, a genetic polymorphism in An. gambiae (sensu stricto) (s.s.), is associated with adaptation to microclimatic differences in humidity and desiccation resistance and mosquito behaviors. Ownership of insecticide-treated bed nets (ITNs) for malaria control has increased markedly in western Kenya in the last 20 years. An increase in the frequency of ITNs indoors could select against house entering or indoor resting of Anopheles mosquitoes. Thus, the frequency of the 2La inversion is postulated to change in An. gambiae (s.s.) with the increase of ITN ownership over time.

Methods

Anopheles gambiae mosquitoes were sampled between 1994 and 2011 using pyrethrum knockdown, bednet traps and human landing catches (HLC) from Asembo and Seme, western Kenya. The 2La inversion was detected by a PCR assay with primers designed for proximal breakpoints of the 2La/a and 2L+^a/+^a chromosomal conformation. Mosquitoes were tested for malaria parasite infection by sporozoite ELISA.

Results

The frequency of the 2La chromosomal inversion declined from 100 % of all chromosomes in 1994 to 17 % in 2005 and remained low through 2011 (21 %). ITN ownership increased from 0 to > 90 % of houses in the study area during this interval. The decline in the frequency of the 2La chromosomal inversion was significantly, negatively correlated with year (r = -0.93) and with increase in ITN ownership (r = -0.96). The frequency of the homo- and heterokaryotypes departed significantly from Hardy-Weinberg equilibrium, suggesting that 2La/a karyotype was under selection, earlier in its favor and later, against it. Precipitation and maximum monthly temperature did not vary over time, therefore there was no trend in climate that could account for the decline. There was no significant difference in frequency of the 2La inversion in An. gambiae (s.s.) females sampled indoors or outdoors in HCL in 2011, nor was there an association between the 2La inversion and infection with Plasmodium falciparum sporozoites.

Conclusions

The increase in ITN ownership in the study area was negatively correlated with the frequency of 2La inversion. The decline in 2La frequency in western Kenya is postulated to be due to differential impacts of ITNs on mosquitoes with different 2La karyotypes, possibly mediated by differences in behavior associated with the 2La karyotypes. Further research is required to determine if this is a widespread phenomenon, to further determine the association of the 2La karyotypes with mosquito behavior, and to assess whether ITNs are exerting selection mediated by differences in behavior on the different karyotypes.

Does insecticide resistance contribute to heterogeneities in malaria transmission in The Gambia?

BACKGROUND

Malaria hotspots, areas with consistently higher than average transmission, may become increasingly common as malaria declines. This phenomenon, currently observed in The Gambia, may be caused by several factors, including some related to the local vectors, whose contribution is poorly understood.

METHODS

Using WHO susceptibility bioassays, insecticide resistance status was determined in vector populations sampled from six pairs of villages across The Gambia, each pair contained a low and high prevalence village.

RESULTS

Three vector species were observed (23.5% Anopheles arabiensis, 31.2% Anopheles gambiae, 43.3% Anopheles coluzzii and 2.0% An. coluzzii × An. gambiae hybrids). Even at a fine scale, significant differences in species composition were detected within village pairs. Resistance to both DDT and deltamethrin was more common in An. gambiae, most markedly in the eastern part of The Gambia and partly attributable to differing frequencies of resistance mutations. The Vgsc-1014F target site mutation was strongly associated with both DDT (OR = 256.7, (95% CI 48.6-6374.3, p < 0.001) and deltamethrin survival (OR = 9.14, (95% CI 4.24-21.4, p < 0.001). A second target site mutation, Vgsc-1575Y, which co-occurs with Vgsc-1014F, and a metabolic marker of resistance, Gste2-114T, conferred additional survival benefits to both insecticides. DDT resistance occurred significantly more frequently in villages with high malaria prevalence (p = 0.025) though this did not apply to deltamethrin resistance.

CONCLUSION

Whilst causality of relationships requires further investigation, variation in vector species and insecticide resistance in The Gambia is associated with malaria endemicity; with a notably higher prevalence of infection and insecticide resistance in the east of the country. In areas with heterogeneous malaria transmission, the role of the vector should be investigated to guide malaria control interventions.

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.

Characterization of swarming and mating behaviour between Anopheles coluzzii and Anopheles melas in a sympatry area of Benin

The swarm structure of two sibling species, Anopheles gambiae coluzzii and Anopheles melas, was characterize to explore the ecological and environmental parameters associated with the formation of swarms and their spatial distribution. Swarms and breeding sites were searched and sampled between January and December 2010, and larval and adult samples were identified by PCR. During the dry season, 456 swarms of An. gambiae s.l. were sampled from 38 swarm sites yielding a total of 23,274 males and 76 females. Of these 38 swarming sites, 18 were composed exclusively of An. gambiae coluzzii and 20 exclusively of An. melas, presenting clear evidence of reproductive swarm segregation. The species makeup of couples sampled from swarms also demonstrated assortative mating. The swarm site localization was close to human dwellings in the case of the An. gambiae coluzzii and on salt production sites for An. melas. At the peak of the rainy season, swarms of An. melas were absent. These findings offer evidence that the ecological speciation of these two sibling species of mosquitoes is associated with spatial swarm segregation and assortative mating, providing strong support for the hypothesis that mate recognition is currently maintaining adaptive differentiation and promoting ecological speciation. Further studies on the swarming and mating systems of An. gambiae, with the prospect of producing a predictive model of swarm distribution, are needed to inform any future efforts to implement strategies based on the use of GMM or SIT.

Diversity and transmission competence in lymphatic filariasis vectors in West Africa, and the implications for accelerated elimination of Anopheles-transmitted filariasis

Lymphatic Filariasis (LF) is targeted for elimination by the Global Programme for the Elimination of Lymphatic Filariasis (GPELF). The strategy adopted is based on the density dependent phenomenon of Facilitation, which hypothesizes that in an area where the vector species transmitting Wuchereria bancrofti are Anopheles mosquitoes, it is feasible to eliminate LF using Mass Drug Administration (MDA) because of the inability of Anopheles species to transmit low-density microfilaraemia. Even though earlier studies have shown Anopheles species can exhibit the process of Facilitation in West Africa, observations point towards the process of Limitation in certain areas, in which case vector control is recommended. Studies on Anopheles species in West Africa have also shown genetic differentiation, cryptic taxa and speciation, insecticide resistance and the existence of molecular and chromosomal forms, all of which could influence the vectorial capacity of the mosquitoes and ultimately the elimination goal. This paper outlines the uniqueness of LF vectors in West Africa and the challenges it poses to the 2020 elimination goal, based on the current MDA 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.

Cytogenetic map for Anopheles nili: application for population genetics and comparative physical mapping

Anopheles nili is one of the major malaria vectors in Africa with a wide geographic distribution. However, the taxonomic and population genetic studies on this species are scarce. New research tools are urgently needed to genetically characterize this important malaria vector. In this study, a high-resolution cytogenetic map was developed for An. nili polytene chromosomes. Chromosomes were straightened and subdivided into 46 numbered divisions according to the banding pattern. Population analysis of An. nili females collected in Burkina Faso revealed the presence of two highly polymorphic inversions on the 2R chromosomal arm. A statistically significant departure from Hardy-Weinberg equilibrium due to a deficit in heterozygotes was detected for inversion 2Rb. To determine chromosome homologies and gene order conservation between An. nili and other major malaria vectors, PCR probes based on the An. gambiae coding sequences were mapped to An. nili chromosomes. Comparative mapping demonstrated that An. nili chromosomes have an An. stephensi-like arm association and that whole-arm translocations and paracentric inversions were the major types of rearrangement in evolution of these mosquitoes. The minimum number of fixed inversions among An. nili, An. gambiae, and An. stephensi was calculated using the Multiple Genome Rearrangements (MGR), Genome Rearrangements In Man and Mouse (GRIMM), and Sorting Permutation by Reversals and block-INterchanGes (SPRING) programs. The data suggest that the An. nili is, at least, as diverged from An. gambiae as An. stephensi. We provide evidence that 2La/a arrangement of An. gambiae is present in outgroup species An. nili and An. stephensi confirming the ancestral status of the 2La inversion in the An. gambiae complex. Availability of the new polytene chromosome map, polymorphic inversions, and physically mapped DNA markers for An. nili will further stimulate population genetic, taxonomic, and genomic studies of this neglected malaria vector.

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.

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.

High-resolution cytogenetic map for the African malaria vector Anopheles gambiae

Cytogenetic and physical maps are indispensible for precise assembly of genome sequences, functional characterization of chromosomal regions, and population genetic and taxonomic studies. We have created a new cytogenetic map for Anopheles gambiae by using a high-pressure squash technique that increases overall band clarity. To link chromosomal regions to the genome sequence, we attached genome coordinates, based on 302 markers of bacterial artificial chromosome, cDNA clones, and PCR-amplified gene fragments, to the chromosomal bands and interbands at approximately a 0.5-1 Mb interval. In addition, we placed the breakpoints of seven common polymorphic inversions on the map and described the chromosomal landmarks for the arm and inversion identification. The map's increased resolution can be used to further enhance physical mapping, improve genome assembly, and stimulate epigenomic studies of malaria vectors.

Spatial swarm segregation and reproductive isolation between the molecular forms of Anopheles gambiae

Anopheles gambiae, the major malaria vector in Africa, can be divided into two subgroups based on genetic and ecological criteria. These two subgroups, termed the M and S molecular forms, are believed to be incipient species. Although they display differences in the ecological niches they occupy in the field, they are often sympatric and readily hybridize in the laboratory to produce viable and fertile offspring. Evidence for assortative mating in the field was recently reported, but the underlying mechanisms awaited discovery. We studied swarming behaviour of the molecular forms and investigated the role of swarm segregation in mediating assortative mating. Molecular identification of 1145 males collected from 68 swarms in Donéguébougou, Mali, over 2 years revealed a strict pattern of spatial segregation, resulting in almost exclusively monotypic swarms with respect to molecular form. We found evidence of clustering of swarms composed of individuals of a single molecular form within the village. Tethered M and S females were introduced into natural swarms of the M form to verify the existence of possible mate recognition operating within-swarm. Both M and S females were inseminated regardless of their form under these conditions, suggesting no within-mate recognition. We argue that our results provide evidence that swarm spatial segregation strongly contributes to reproductive isolation between the molecular forms in Mali. However this does not exclude the possibility of additional mate recognition operating across the range distribution of the forms. We discuss the importance of spatial segregation in the context of possible geographic variation in mechanisms of reproductive isolation.

Exploring the origin and degree of genetic isolation of Anopheles gambiae from the islands of São Tomé and Príncipe, potential sites for testing transgenic-based vector control

The evolutionary processes at play between island and mainland populations of the malaria mosquito vector Anopheles gambiae sensu stricto are of great interest as islands may be suitable sites for preliminary application of transgenic-based vector control strategies. São Tomé and Príncipe, located off the West African coast, have received such attention in recent years. This study investigates the degree of isolation of An. gambiae s.s. populations between these islands and the mainland based on mitochondrial and ribosomal DNA molecular data. We identify possible continental localities from which these island populations derived. For these purposes, we used F_ST values, haplotype networks, and nested clade analysis to estimate migration rates and patterns. Haplotypes from both markers are geographically widespread across the African continent. Results indicate that the populations from São Tomé and Príncipe are relatively isolated from continental African populations, suggesting they are promising sites for test releases of transgenic individuals. These island populations are possibly derived from two separate continental migrations. This result is discussed in the context of the history of the African slave trade with respect to São Tomé and Príncipe.

Spatial distribution of the chromosomal forms of anopheles gambiae in Mali

Background

Maps of the distribution of malaria vectors are useful tools for stratification of malaria risk and for selective vector control strategies. Although the distribution of members of the Anopheles gambiae complex is well documented in Africa, a continuous map of the spatial distribution of the chromosomal forms of An. gambiae s.s. is not yet available at country level to support control efforts.

Methods

Bayesian geostatistical methods were used to produce continuous maps of the spatial distribution of the chromosomal forms of An. gambiae s.s. (Mopti, Bamako, Savanna and their hybrids/recombinants) based on their relative frequencies in relation to climatic and environmental factors in Mali.

Results

The maps clearly show that each chromosomal form favours a particular defined eco-climatic zone. The Mopti form prefers the dryer northern Savanna and Sahel and the flooded/irrigated areas of the inner delta of the Niger River. The Savanna form favours the Sudan savanna areas, particularly the South and South-Eastern parts of the country (Kayes and Sikasso regions). The Bamako form has a strong preference for specific environmental conditions and it is confined to the Sudan savanna areas around urban Bamako and the Western part of Sikasso region. The hybrids/recombinants favour the Western part of the country (Kayes region) bordering the Republic of Guinea Conakry.

Conclusion

The maps provide valuable information for selective vector control in Mali (insecticide resistance management) and may serve as a decision support tool for the basis for future malaria control strategies including genetically manipulated mosquitoes.

Anopheles gambiae complex along The Gambia river, with particular reference to the molecular forms of An. gambiae s.s

BACKGROUND

The geographic and temporal distribution of M and S molecular forms of the major Afrotropical malaria vector species Anopheles gambiae s.s. at the western extreme of their range of distribution has never been investigated in detail.

MATERIALS AND METHODS

Collections of indoor-resting An. gambiae s.l. females were carried out along a ca. 400 km west to east transect following the River Gambia from the western coastal region of The Gambia to south-eastern Senegal during 2005 end of rainy season/early dry season and the 2006 rainy season. Specimens were identified to species and molecular forms by PCR-RFLP and the origin of blood-meal of fed females was determined by ELISA test.

RESULTS

Over 4,000 An. gambiae s.l. adult females were collected and identified, 1,041 and 3,038 in 2005 and 2006, respectively. M-form was mainly found in sympatry with Anopheles melas and S-form in the western part of the transect, and with Anopheles arabiensis in the central part. S-form was found to prevail in rural Sudan-Guinean savannah areas of Eastern Senegal, in sympatry with An. arabiensis. Anopheles melas and An. arabiensis relative frequencies were generally lower in the rainy season samples, when An. gambiae s.s. was prevailing. No large seasonal fluctuations were observed for M and S-forms. In areas where both M and S were recorded, the frequency of hybrids between them ranged from to 0.6% to 7%.

DISCUSSION

The observed pattern of taxa distribution supports the hypothesis of a better adaptation of M-form to areas characterized by water-retaining alluvial deposits along the Gambia River, characterized by marshy vegetation, mangrove woods and rice cultivations. In contrast, the S-form seems to be better adapted to free-draining soil, covered with open woodland savannah or farmland, rich in temporary larval breeding sites characterizing mainly the eastern part of the transect, where the environmental impact of the Gambia River is much less profound and agricultural activities are mainly rain-dependent. Very interestingly, the observed frequency of hybridization between the molecular forms along the whole transect was much higher than has been reported so far for other areas.

CONCLUSION

The results support a bionomic divergence between the M and S-forms, and suggest that the western extreme of An. gambiae s.s. geographical distribution may represent an area of higher-than-expected hybridization between the two molecular forms.

The molecular forms of Anopheles gambiae: a phenotypic perspective

The African malaria mosquito Anopheles gambiae is undergoing speciation, being split into the M and S molecular forms. Speciation is the main process promoting biological diversity, thus, new vector species might complicate disease transmission. Genetic differentiation between the molecular forms has been extensively studied, but phenotypic differences between them, the evolutionary forces that generated divergence, and the mechanisms that maintain their genetic isolation have only recently been addressed. Here, we review recent studies suggesting that selection mediated by larval predation and competition promoted divergence between temporary and permanent freshwater habitats. These differences explain the sharp discontinuity in distribution of the molecular forms between rice fields and surrounding savanna, but they can also explain the concurrent cline between humid and arid environments due to the dependence on permanent habitats in the latter. Although less pronounced, differences in adult body size, reproductive output, and longevity also suggest that the molecular forms have adapted to distinct niches. Reproductive isolation between the molecular forms is achieved by spatial swarm segregation, although within-swarm mate recognition appears to play a role in certain locations. The implications of these results to disease transmission and control are discussed and many of the gaps in our understanding are highlighted.

Localization of candidate regions maintaining a common polymorphic inversion (2La) in Anopheles gambiae

Chromosomal inversion polymorphisms are thought to play a role in adaptive divergence, but the genes conferring adaptive benefits remain elusive. Here we study 2La, a common polymorphic inversion in the African malaria vector Anopheles gambiae. The frequency of 2La varies clinally and seasonally in a pattern suggesting response to selection for aridity tolerance. By hybridizing genomic DNA from individual mosquitoes to oligonucleotide microarrays, we obtained a complete map of differentiation across the A. gambiae genome. Comparing mosquitoes homozygous for the 2La gene arrangement or its alternative (2L+^a), divergence was highest at loci within the rearranged region. In the 22 Mb included within alternative arrangements, two ∼1.5 Mb regions near but not adjacent to the breakpoints were identified as being significantly diverged, a conclusion validated by targeted sequencing. The persistent association of both regions with the 2La arrangement is highly unlikely given known recombination rates across the inversion in 2La heterozygotes, thus implicating selection on genes underlying these regions as factors responsible for the maintenance of 2La. Polymorphism and divergence data are consistent with a model in which the inversion is maintained by migration-selection balance between multiple alleles inside these regions, but further experiments will be needed to fully distinguish between the epistasis (coadaptation) and local adaptation models for the maintenance of 2La.

A chromosomal inversion occurs when part of the chromosome breaks, rotates 180 degrees, and rejoins the broken chromosome. The result is a chromosome carrying a segment whose gene order is reversed. Whereas the physical rearrangement itself may have no direct consequences on gene function, recombination between alleles in the rearranged and wild type segments is suppressed. If multiple alleles inside the inverted or original orientations are well adapted to contrasting environmental conditions, suppressed recombination provides a mechanism to keep beneficial allelic combinations from being shuffled between different genetic backgrounds. Working with wild populations of flies, Dobzhansky provided the first evidence that selection was key to maintaining inversion polymorphism. Subsequently, examples of inversion polymorphisms under selection in other organisms have been found, notably in the mosquito that transmits most cases of human malaria, Anopheles gambiae. However, the genes or gene regions conferring fitness advantages have yet to be discovered. In this study, the authors used modern genomics tools to map such regions in an inversion at an unprecedented level of detail, and show that these regions are likely to be responsible for the maintenance of the inversion polymorphism in natural populations. This study lays the groundwork for future efforts to identify the genes themselves and their role in adaptation.

Evidence for late Pleistocene population expansion of the malarial mosquitoes, Anopheles arabiensis and Anopheles gambiae in Nigeria

Anopheles gambiae Giles s.s. and Anopheles arabiensis Patton (Diptera: Culicidae) are major vectors of malaria in Nigeria. We used 1115 bp of the mitochondrial COI gene to assess their population genetic structures based on samples from across Nigeria (n = 199). The mtDNA neighbour-joining tree, based on F(ST) estimates, separated An. gambiae M and S forms, except that samples of An. gambiae M from Calabar clustered with all the An. gambiae S form. Anopheles arabiensis and An. gambiae could be combined into a single star-shaped, parsimonious haplotype network, and shared three haplotypes. Haplotype diversity values were high in An. arabiensis and An. gambiae S, and intermediate in An. gambiae M; all nucleotide diversities were relatively low. Taken together, patterns of haplotype diversity, the star-like genealogy of haplotypes, five of seven significant neutrality tests, and the violation of the isolation-by-distance model indicate population expansion in An. arabiensis and An. gambiae S, but the signal was weak in An. gambiae M. Selection is supported as an important factor shaping genetic structure in An. gambiae in Nigeria. There were two geographical subdivisions in An. arabiensis: one included all southern localities and all but two central localities; the other included all northern and two central localities. Re-analysing an earlier microsatellite dataset of An. arabiensis using a Bayesian method determined that there were two distinctive clusters, northern and southern, that were fairly congruent with the mtDNA subdivisions. There was a trend towards decreasing genetic diversity in An. arabiensis from the northern savannah to the southern rainforest that corroborated previous data from microsatellites and polytene chromosomes.

Malaria in Africa: vector species' niche models and relative risk maps

A central theoretical goal of epidemiology is the construction of spatial models of disease prevalence and risk, including maps for the potential spread of infectious disease. We provide three continent-wide maps representing the relative risk of malaria in Africa based on ecological niche models of vector species and risk analysis at a spatial resolution of 1 arc-minute (9 185 275 cells of approximately 4 sq km). Using a maximum entropy method we construct niche models for 10 malaria vector species based on species occurrence records since 1980, 19 climatic variables, altitude, and land cover data (in 14 classes). For seven vectors (Anopheles coustani, A. funestus, A. melas, A. merus, A. moucheti, A. nili, and A. paludis) these are the first published niche models. We predict that Central Africa has poor habitat for both A. arabiensis and A. gambiae, and that A. quadriannulatus and A. arabiensis have restricted habitats in Southern Africa as claimed by field experts in criticism of previous models. The results of the niche models are incorporated into three relative risk models which assume different ecological interactions between vector species. The “additive” model assumes no interaction; the “minimax” model assumes maximum relative risk due to any vector in a cell; and the “competitive exclusion” model assumes the relative risk that arises from the most suitable vector for a cell. All models include variable anthrophilicity of vectors and spatial variation in human population density. Relative risk maps are produced from these models. All models predict that human population density is the critical factor determining malaria risk. Our method of constructing relative risk maps is equally general. We discuss the limits of the relative risk maps reported here, and the additional data that are required for their improvement. The protocol developed here can be used for any other vector-borne disease.

Reproductive output of female Anopheles gambiae (Diptera: Culicidae): comparison of molecular forms

Knowledge of ecological differences between the molecular forms of Anopheles gambiae Giles (Diptera: Culicidae) might lead to understanding of their unique contribution to disease transmission, to better vector control, and to identification of the forces that have separated them. We compared female fecundity measured as egg batch size in relation to body size between the molecular forms in Mali and contrasted them with their sibling species, Anopheles arabiensis Patton. To determine whether eggs of different egg batches are of similar "quality," we compared the total protein content of first-stage larvae (L1s), collected < 2 h after hatching in deionized water. Egg batch size significantly varied between An. gambiae and An. arabiensis and between the molecular forms of An. gambiae (mean batch size was 186.3, 182.5, and 162.0 eggs in An. arabiensis and the M and the S molecular form of An. gambiae, respectively). After accommodating female body size, however, the difference in batch size was not significant. In the S molecular form, egg protein content was not correlated with egg batch size (r = -0.08, P > 0.7) nor with female body size (r = -0.18, P > 0.4), suggesting that females with more resources invest in more eggs rather than in higher quality eggs. The mean total protein in eggs of the M form (0.407 microg per L1) was 6% higher than that of the S form (0.384 microg per L1), indicating that the M form invests a greater portion of her resources into current (rather than future) reproduction. A greater investment per offspring coupled with larger egg batch size may reflect an adaptation of the M form to low productivity larval sites as independent evidence suggests.

Malaria vector analysis and control

Antivector measures in malaria control should aim for a cost-effective reduction of the transmission potential ideally to below the critical level for sustained transmission. The available measures include those that decrease vector abundance, vector-human contact and vector survival rate or that increase the length of the sporogonic cycle. These have widely different impact on malaria transmission, as shown by epidemiological modelling. Direct modification of vector receptivity to Plasmodium is also hypothetically attainable by the use of transmission-blocking vaccines or by genetic manipulation and replacement of the vector population. Vector analysis constitutes the essential prerequisite for basic malaria epidemiology as well as for the development, planning and evaluation of antivector measures. The rationale, the problems and the perspectives of vector analysis are reviewed here by Mario Coluzzi, on the basis of his experience with Afrotropical and Mediterranean malaria vectors.

Multilevel analyses of genetic differentiation in Anopheles gambiae s.s. reveal patterns of gene flow important for malaria-fighting mosquito projects

Malaria control projects based on the introduction and spread of transgenes into mosquito populations depend on the extent of isolation between those populations. On the basis of the distribution of paracentric inversions, Anopheles gambiae has been subdivided into five subspecific chromosomal forms. Estimating gene flow between and within these forms of An. gambiae presents a number of challenges. We compared patterns of genetic divergence (F(ST)) between sympatric populations of the Bamako and Mopti forms at five sites. We used microsatellite loci within the j inversion on chromosome 2, which is fixed in the Bamako form but absent in the Mopti form, and microsatellites on chromosome 3, a region void of inversions. Estimates of genetic diversity and F(ST)'s suggest genetic exchanges between forms for the third chromosome but little for the j inversion. These results suggest a role for the inversion in speciation. Extensive gene flow within forms among sites resulted in populations clustering according to form despite substantial gene flow between forms. These patterns underscore the low levels of current gene flow between chromosomal forms in this area of sympatry. Introducing refractoriness genes in areas of the genome void of inversions may facilitate their spread within forms but their passage between forms may prove more difficult than previously thought.

Species abundance and insecticide resistance of Anopheles gambiae in selected areas of Ghana and Burkina Faso

The Ghanaian National Malaria Control Programme has prioritized insecticide-treated materials as a key strategy for malaria control. We report on a survey of the distribution of the molecular forms of Anopheles gambiae Giles (Diptera: Culicidae) and insecticide resistance (the kdr mutation), carried out by sampling mosquitoes from 11 locations in Ghana and one additional site in Burkina Faso. The molecular M and S forms of An. gambiae were found to occur in sympatry in southern Ghana. The S form predominated throughout its distribution in the coastal savannah, except at one location in the strand and mangrove zone where rice was cultivated. The M form was the only form collected in northern Ghana and was the predominant form (97.5%) in Burkina Faso. No M/S hybrids were detected. The kdr mutation was observed at very high frequencies (98-100%) within the S form but reached a maximum of only 3.38% in the M form in one population at an irrigation scheme in the Ghanaian coastal savannah zone.

The distribution of M and S molecular forms of Anopheles gambiae in Nigeria

The distribution of M and S molecular forms of Anopheles gambiae sensu stricto across Nigeria was determined. The molecular form of 40 to 45 specimens per locality from 9 localities was determined using mostly the same specimens from our recent study of genetic differentiation of A. gambiae across Nigeria (Onyabe & Conn, 2001). These samples were previously genotyped at 10 microsatellite loci, 5 located within chromosome inversions and 5 outside inversions. Both molecular forms occurred throughout the country, with no apparent relationship to the ecological transition from dry savannah in the north to humid forest in southern Nigeria. In all localities, however, 1 form or the other occurred virtually exclusively. No hybrids between forms were found. Across all loci, F(ST) values were as high within molecular forms as between forms. Regardless of molecular form, F(ST) values calculated across loci within inversions were much higher (range 0.0016 to 0.1988) than those calculated across loci outside inversions (range -0.0035 to 0.0260). Genetic distance was not significantly correlated with geographical distance within either form (P> 0.05). These observations suggest that, in addition to partial reproductive barriers between molecular forms, selection is a major factor shaping genetic differentiation of A. gambiae across Nigeria.

The sympatric occurrence of two molecular forms of the malaria vector Anopheles gambiae Giles sensu stricto in Kanyemba, in the Zambezi Valley, Zimbabwe

The complexity of Anopheles gambiae Giles sensu stricto, a sibling species of the A. gambiae complex, has been shown by the extent of chromosomal inversion polymorphisms and more recently, by divergence at the molecular level. Analysis by PCR of a field population of A. gambiae complex from the Zambezi Valley area of Kanyemba, northern Zimbabwe, showed the co-existence of three sibling species: A. arabiensis (90%), A. gambiae (4%), and A. quadriannulatus (6%). Molecular characterization of the A. gambiae s.s. material revealed the sympatric presence of the two molecular forms S and M within a sample of eight. This is the first record of the co-existence of these two molecular forms in southern Africa.

Localized breeding of the Anopheles gambiae complex (Diptera: Culicidae) along the River Gambia, West Africa

A study was undertaken to identify the major larval habitats of the Anopheles gambiae (Giles) complex in rural Gambia. Mosquito larvae and pupae were sampled along transects and in specific habitats in the central region of the country during the rainy seasons of 1996 and 1997. The sampling showed that the major breeding sites were located on the flooded alluvial soils bordering the river. The largest numbers of larvae were found during September, one month after the peak rains. Polymerase chain reaction analysis of specimens showed that Anopheles melas (Theobald) was the dominant species in the flooded areas (81.5%), followed by A. gambiae sensu stricto (Giles) (18.0%) and A. arabiensis (Patton) (0.5%). By sampling in specific habitats it was evident that A. arabiensis was mainly breeding in rain-fed rice fields along the edge of the alluvial soils. Anopheles melas and A. gambiae s.s. often coexisted but whereas A. melas were found in water with a salinity of up to 72% sea water (25.2 g NaCl l(-1)), A. gambiae s.s. only occurred in water with up to 30% sea water (10.5 g NaCl l(-1)). Anopheles melas larvae were found in association with plant communities dominated by sedges and grasses (Eleocharis sp., Paspalum sp., Sporobolus sp.) and sea-purslane Sesuvium portulacastrum (L.) and the presence of cattle hoof prints, whereas A. gambiae s.s. larvae mainly occurred in association with Paspalum sp. and Eleocharis sp. The study showed that even during the peak rainy season, breeding of the A. gambiae complex is almost entirely restricted to the extensive alluvial areas along the river.

Genetic differentiation of the malaria vector Anopheles gambiae across Nigeria suggests that selection limits gene flow

Gene flow was investigated in Anopheles gambiae from eight localities that span the ecological zones of Nigeria (arid savanna zones in the north gradually turn into humid forest zones in the south). Genetic differentiation was measured over 10 microsatellite loci and, to determine any effects of selection, five loci were located within chromosome inversions and the other five were outside inversions. Over all loci, the largest estimates of differentiation were in comparisons between localities in the savanna vs. forest zones (range FST 0.024-0.087, Nm 2.6-10.1; RST 0.014-0.100, Nm 2.2-16.4). However, three loci located within inversions on chromosome II, whose frequencies varied clinically from north to south, were responsible for virtually all of the differentiation. When the three loci were removed, genetic distances across the remaining seven loci were markedly reduced even between localities in the forest and savanna zones (range FST 0.001-0.019, Nm 12.7-226.1) or no longer significant (P > 0.05) in the case of RST. Although tests of isolation by distance gave seemingly equivocal results, geographical distance does not appear to limit gene flow. These observations suggest that gene flow is extensive across the country but that selection on genes located within some inversions on chromosome II counters the homogenizing effects of gene flow.

Mapping distributions of chromosomal forms of Anopheles gambiae in West Africa using climate data

The mosquito Anopheles gambiae Giles sensu stricto (Diptera: Culicidae), the principal vector of malaria in West Africa, comprises several chromosomal forms (e.g. Bissau, Forest, Mopti, Savanna) associated with climatic zones. Here we show how climate data can be used to map the geographical distribution of these chromosomal forms. The climate at 144 sites surveyed for mosquitoes in West Africa between 1971 and 92 was determined using computerized climate surfaces. Forest and Bissau forms occurred at relatively wet sites: median annual precipitation 1325 mm and 1438 mm, respectively, interquartile ranges (IQR) 1144-1858 mm and 1052-1825 mm), whilst the Mopti form was found at dry sites (annual 938 mm, IQR 713-1047 mm) and the Savanna form at sites intermediate between the wet and dry forms (annual 1067 mm, IQR 916-1279). Logistic regression analyses of the climate variables were carried out on a stratified random sample of half the sites. The resulting models correctly classified over 80% of the sites for presence or absence of each chromosomal form. When these models were tested against excluded sites they were also correct at over 80% of sites. The combined data produced models that were correct at over 86% of sites. Mean annual precipitation, evapotranspiration, minimum temperature and maximum temperature were the most important climate variables correlated with the distribution of these forms of An. gambiae. We used the logistic models to map the distribution of each chromosomal form within the reported range for An. gambiae s.s. in West Africa employing a geographical information system. Our maps indicate that each chromosomal form favours particular climate envelopes in well-defined ecoclimatic zones, although these forms are sympatric at the edges of their ranges. This study demonstrates that climate can be used to map the distribution of chromosomal forms of insects across large areas.

DNA analysis of transferred sperm reveals significant levels of gene flow between molecular forms of Anopheles gambiae

Anopheles gambiae populations in west Africa are complex, being composed of multiple, sympatric subpopulations. Recent studies have failed to reveal significant genetic differences among subpopulations, stimulating a debate regarding the levels of gene flow among them. The observed homogeneity may be the consequence of substantial contemporary gene flow or it may be that reproductive isolation is complete, but too recent for the accumulation of significant levels of genic divergence. Here, we report the results of a study estimating contemporary levels of gene flow between An. gambiae subpopulations by analysing females and transferred sperm removed from their reproductive systems. A total of 251 female and associated sperm extracts was analysed from a single site in Mali. Two molecular forms of An. gambiae, the M- and S-forms, occurred in sympatry at this site. Overall, we found very strong positive assortative mating within forms, however, we did observe significant hybridization between forms. In the M subpopulation 2/195 females (1.03%) contained sperm from S-form males and in 55 S-form females we found one female containing M-form sperm (1.82%). We also identified a mated M xS hybrid adult female. From mating frequencies, we estimate the Nem between the M- and S-form at 16.8, and from the adult hybrid frequency at 5.6. These values are consistent with our earlier estimate, based on FST for 21 microsatellite loci in which Nem = 5.8. We conclude that the general lack of genetic divergence between the M and S subpopulations of An. gambiae can be explained entirely by contemporary gene flow.

Microsatellite DNA polymorphism and heterozygosity among field and laboratory populations of Anopheles gambiae ss (Diptera: Culicidae)

We compared microsatellite polymorphism at nine loci located on chromosome 3 among two colonies and a field population of Anopheles gambiae sensu stricto Giles mosquitoes. Numbers of microsatellite alleles observed at each locus and mean heterozygosities were drastically reduced among laboratory colonies. Genetic analysis of the field population used in this study revealed an unprecedented frequency of rare alleles (<0.05). In contrast, colony samples revealed large numbers of alleles with frequencies >0.50. Partitioning of field data to assess the impact of rare alleles, null alleles, and sample size on estimates of mean heterozygosity revealed the plasticity of this measurement and suggests that heterozygosity may be reliably estimated from relatively small collections using microsatellites.

Gene flow among populations of the malaria vector, Anopheles gambiae, in Mali, West Africa

The population structure of the Anopheles gambiae complex is unusual, with several sibling species often occupying a single area and, in one of these species, An. gambiae sensu stricto, as many as three "chromosomal forms" occurring together. The chromosomal forms are thought to be intermediate between populations and species, distinguishable by patterns of chromosome gene arrangements. The extent of reproductive isolation among these forms has been debated. To better characterize this structure we measured effective population size, N(e), and migration rates, m, or their product by both direct and indirect means. Gene flow among villages within each chromosomal form was found to be large (N(e)m > 40), was intermediate between chromosomal forms (N(e)m approximately 3-30), and was low between species (N(e)m approximately 0.17-1.3). A recently developed means for distinguishing among certain of the forms using PCR indicated rates of gene flow consistent with those observed using the other genetic markers.

Attempts to molecularly distinguish cryptic taxa in Anopheles gambiae s.s

Analyses of inversions in polytene chromosomes indicate that, in West Africa, Anopheles gambiae (sensu stricto) may be a complex of more than a single taxonomic unit, and these units have been called chromosomal forms. In order to determine whether this genetic discontinuity extends to the rest of the genome, as would be expected if reproductive isolation exists, we have sequenced several regions of both the nuclear and mitochondrial genomes. With one exception, we were unable to identify any nucleotide sites that differentiate the chromosomal forms. The exception was the internal transcribed spacer (ITS) of the ribosomal DNA (rDNA). Three sites in this region distinguish Mopti chromosomal form from Savanna and Bamako in Mali and Burkina Faso. However, outside these two countries, the association between chromosomal form and rDNA type does not always hold. Together with the variants in the rDNA intergenic spacer (IGS) described in the accompanying papers (della Torre et al., 2001; Favia et al., 2001), we can recognize two major types of rDNA, Type I and Type II (corresponding to molecular forms S and M in della Torre et al., 2001). Type I is widespread in West Africa and is the only type found outside of West Africa (i.e. Tanzania and Madagascar). Type II is confined to West Africa. We were unable to detect any heterozygosity for the ITS types even in five collections containing both types. A sample from the island of São Tomé could not be classified into either Type I or Type II as the rDNA had characteristics of both. In general, our results confirm that An. gambiae is not a single pan-mictic unit, but exactly how to define any new taxa remains problematic. Finally, we have found minor variants of the major rDNA types fixed in local populations; contrary to most previous studies, this suggests restricted gene flow among populations of this species.

Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa

We karyotyped and identified by polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) analysis Anopheles gambiae s.s. samples collected in several African countries. The data show the existence of two non-panmictic molecular forms, named S and M, whose distribution extended from forest to savannahs. Mosquitoes of the S and M forms are homosequential standard for chromosome-2 inversions in forest areas. In dry savannahs, S is characterized mainly by inversion polymorphisms typical of Savanna and Bamako chromosomal forms, while M shows chromosome-2 arrangements typical of Mopti and/or Savanna and/or Bissau, depending on its geographical origin. Chromosome-2 inversions therefore seem to be involved in ecotypic adaptation rather than in mate-recognition systems. Strong support for the reproductive isolation of S and M in Ivory Coast comes from the observation that the kdr allele is found at high frequencies in S specimens and not at all in chromosomal identical M specimens. However, the kdr allele does not segregate with molecular forms in Benin.

Distribution of genetic variation among chromosomal forms of Anopheles gambiae s.s: introgressive hybridization, adaptive inversions, or recent reproductive isolation?

A series of four papers in this issue explores the reproductive status of the five chromosomal forms of An. gambiae s.s. using molecular techniques to examine the variation among twelve different genes located throughout the An. gambiae s.s. genome. Results of these and previous studies are consistent with a hypothesis of at least partial barriers to gene flow between some chromosomal forms in the Ivory Coast and other West African countries to the north and west, but introgression between S and M types in Benin and countries to the east. Collectively, these studies indicate the need for a broader geographical sampling of An. gambiae s.s., increased research on mechanisms of prezygotic reproductive isolation and field-based studies of survival and fecundity in hybrids to test for postzygotic reproductive isolation.

The kdr mutation occurs in the Mopti form of Anopheles gambiae s.s. through introgression

Anopheles gambiae s.s. is a complex of sibling taxa characterized by various paracentric inversions. In west and central Africa, where several taxa are sympatric, a kdr mutation responsible for pyrethroid resistance has been described in only one (the S taxon), suggesting an absence of gene flow between them. Following a thorough sampling, we have found a kdr mutation in another taxon (M). To establish whether this mutation is the same event or not, the large intron upstream of the kdr mutation was sequenced to find polymorphic sites in susceptible/resistant and M/S mosquitoes. The low genetic diversity found in this DNA region indicates that a local genetic sweep has recently occurred. However, some polymorphic sites were found, and it is therefore concluded that the kdr mutation in the M taxon is not an independent mutation event, and is best explained by an introgression from the S taxon. These results are discussed within the context of possible gene flow between members of An. gambiae s.s. taxa, and with the possible spread of the kdr mutation in other closely related malaria vectors of the An. gambiae complex.

Complexities in the genetic structure of Anopheles gambiae populations in west Africa as revealed by microsatellite DNA analysis

Chromosomal forms of Anopheles gambiae, given the informal designations Bamako, Mopti, and Savannah, have been recognized by the presence or absence of four paracentric inversions on chromosome 2. Studies of karyotype frequencies at sites where the forms occur in sympatry have led to the suggestion that these forms represent species. We conducted a study of the genetic structure of populations of An. gambiae from two villages in Mali, west Africa. Populations at each site were composed of the Bamako and Mopti forms and the sibling species, Anopheles arabiensis. Karyotypes were determined for each individual mosquito and genotypes at 21 microsatellite loci determined. A number of the microsatellites have been physically mapped to polytene chromosomes, making it possible to select loci based on their position relative to the inversions used to define forms. We found that the chromosomal forms differ at all loci on chromosome 2, but there were few differences for loci on other chromosomes. Geographic variation was small. Gene flow appears to vary among different regions within the genome, being lowest on chromosome 2, probably due to hitchhiking with the inversions. We conclude that the majority of observed genetic divergence between chromosomal forms can be explained by forces that need not involve reproductive isolation, although reproductive isolation is not ruled out. We found low levels of gene flow between the sibling species Anopheles gambiae and Anopheles arabiensis, similar to estimates based on observed frequencies of hybrid karyotypes in natural populations.

[Malaria transmission in the rural zone of Niakhar, Senegal]

The anopheline bioecology and the malaria transmission were studied from January to December 1995 in three villages of the sahelian rural area of Niakhar, Senegal. This area of 29000 inhabitants, has been for several decades, a regional observatory for population and health. The three methods used for collecting mosquitoes were the collection at larval stages, the all night human biting collection, and the pyrethrum spray catch in houses during afternoons. The anophelines collected were, by numerical importance: Anopheles arabiensis, An. rufipes, An. gambiae, An. pharoensis, An. funestus and An. coustani. In the An. gambiae complex, An. arabiensis represented 97% of man biting females and 98% of half gravid resting females (difference not significant); the other reminding species of this complex was always An. gambiae. These two species belonging to the An. gambiae complex were responsible for the totality of the transmission. The anthropophilic index, obtained from half gravid indoor resting An. gambiae s.l., was 83%. The annual biting rate of An. gambiae s.l. varied from 512 to 1558 bites per man per night, depending on the villages. Vectors were observed all year long but their densities were low during the dry season. Vector population presented a notable increase due to the rains, with a maximum of about 10 bites per man per night in September or at the beginning of October; during September the biting rate represented 48% of the annual biting rate. The sporozoitic index of An. gambiae s.l., obtained by ELISA revealing the circumsporozoite protein, was 1.6% for human biting females and 1.8% for half-gravid resting females (difference not significant). Plasmodium falciparum was the only plasmodial species observed among infected anophelines. The annual transmission in the vo villages representative of the Niakhar area were 9 and 12 bites of infected anophelines per man, occurring mainly from August to October. In the third village, not representative of the area regarding permanent breeding places, the transmission was 26 bites of infected anopheline per man per year. These results were discussed in the Senegambian and sahelian contexts.

Four years' entomological study of the transmission of seasonal malaria in Senegal and the bionomics of Anopheles gambiae and A. arabiensis

From 1993 to 1996, an entomological survey was conducted in the village of Ndiop, Senegal, as part of a research programme on malaria epidemiology and the mechanisms of protective immunity. Mosquitoes were captured on human bait and by indoor spraying. Species from the Anopheles gambiae complex were identified using the polymerase chain reaction, and Plasmodium falciparum infections were detected by enzyme-linked immunosorbent assay for circumsporozoite protein. The vector species identified were A. gambiae (33.9%), A. arabiensis (63.2%), A. melas (0.3%) and A. funestus (2.5%). Similar proportions of A. gambiae (74.2%) and A. arabiensis (73.8%) contained human blood; 27.0% of A. gambiae and 28.3% of A. arabiensis had fed on cattle. The sporozoite rates were similar for A. gambiae (3.2%) and A. arabiensis (3.7%). The annual entomological inoculation rates varied greatly depending on the year. There were 63, 17, 37 and 7 infected bites per person per year in 1993, 1994, 1995 and 1996 respectively. Transmission was highly seasonal, from July to October. A. arabiensis was responsible for 66% of malaria transmission, A. gambiae for 31%, and A. funestus for 3%.

Ecological genetic studies in the chromosomal form Mopti of Anopheles gambiae s.str. in Mali, west Africa

Among the sibling species of the Afrotropical Anopheles gambiae complex, the nominal taxon (An. gambiae s.str.) is the major malaria vector. Its bionomics suggest a man-dependent speciation process which involves, in West Africa, various incipient species chromosomally recognized by different combinations of 2R paracentric inversions. One of the most recent evolutionary steps of such a speciation process appears to be the chromosomal form Mopti, which is associated with dry season irrigation in arid zones, and is characterized by a remarkable ecological flexibility related to three 2R alternative arrangements, namely bc, u and +, whose expected karyotypes are found in Hardy-Weinberg equilibrium. The study of this chromosomal polymorphism in samples from a 16-locality transect in Mali shows wide variations and highly significant correlation with both temporal and spatial climatic differences. Mosquitoes homokaryotypic for 2Rbc are the actual dry season and arid areas breeders. The regular rise of 2Rbc frequency, up to fixation, during each dry season, corresponds to the South-North clinal increase of the same arrangement along the transect, from about 30% in the humid savanna to near fixation in the South-Saharan zone. This coherent ecological genetics case provides full support to the hypothesis of the adaptive nature of paracentric inversions. Moreover, the very peculiar system of combinations of contiguous 2R inversions, utilized by Mopti as well as by other chromosomal forms of An. gambiae, suggests a process of polygenic reorganization based on linkage disequilibria and involving the inversions as driving selection units.

A malaria control trial using insecticide-treated bed nets and targeted chemoprophylaxis in a rural area of The Gambia, west Africa. 1. A review of the epidemiology and control of malaria in The Gambia, west Africa

Malaria was recognized as an important cause of death among early European visitors to The Gambia, but the infection was first studied systematically in the local population only in the 1950s. Studies undertaken in the village of Keneba at that time showed that nearly all children under the age of 5 years had parasitaemia throughout the year. More recent surveys in rural areas of The Gambia have shown much lower levels of parasitaemia, probably as a result of a decline in rainfall in The Gambia during the past 30 years and because of an increase in the availability of anti-malarial drugs. Nevertheless, community surveys and reviews of hospital statistics show that malaria is still one of the most important causes of death among Gambian children; about 1 in 25 rural Gambian children die from malaria before reaching the age of 5 years. Until recently, malaria control in The Gambia relied upon prompt treatment of clinical attacks, first with quinine and more recently with chloroquine, and upon some limited vector control in the capital, Banjul. However, during the past few years, it has been shown that mortality in rural children can be reduced substantially by means of chemoprophylaxis given by village health workers. Bed nets (mosquito nets) are used widely in The Gambia and epidemiological surveys have shown an association between the use of bed nets and protection against malaria. This observation led to a series of small scale intervention trials.(ABSTRACT TRUNCATED AT 250 WORDS)

Malaria in a peri-urban area of The Gambia

A clinical and entomological survey of malaria was carried out in Bakau, a peri-urban coastal settlement in The Gambia, from June 1988-May 1989. Only 41 of a cohort of 560 children, aged from three months to nine-years-old, experienced a clinical episode of malaria during the observation period. The majority of cases were identified at clinics and not by regular community surveillance. In Bakau Old Town episodes of malaria were more common on the periphery of the settlement, adjacent to typical anopheline breeding sites, than in the centre. Overall malaria cases were not significantly clustered in space and time, although three pairs of cases among children sleeping in the same room at the same time were identified. A cross-sectional survey in November, at the end of the rainy season, revealed a point prevalence parasitaemia of 2.0% and a spleen rate of 0.3%. All malariometric parameters measured were much lower than any found in comparable studies undertaken in rural areas of the country, reflecting the low number of malaria vectors, Anopheles gambiae complex mosquitoes, found in Bakau. Chloroquine consumption, sleeping under bednets, houses with ceilings, the use of insecticide aerosols and burning traditional mosquito repellents may all have contributed to the low prevalence of malaria in the study area.