Cytogenetics of the Anopheles gambiae complex in Sudan, with special reference to An. arabiensis: relationships with East and West African populations

The species composition of malaria vector mosquitoes belonging to the Anopheles gambiae complex (Diptera: Culicidae) from >40 localities in Sudan, representing most ecological situations, was determined by analysis of ovarian polytene chromosomes. Of 2162 females, 93% were identified as An. arabiensis Patton and 7% were An. gambiae Giles sensu stricto. No hybrids were found between the two species. Anopheles arabiensis occurred in all but two sites, whereas An. gambiae s.s. was effectively limited to the southernmost, more humid localities. For chromosomal paracentric inversions, the degree of polymorphism was low in An. gambiae s.s. (inversions 2La, 2Rb and 2Rd), higher in An. arabiensis (inversions Xe, 2Ra, b, bc, d1, s; 3Ra, d). Anopheles gambiae samples from Sudan were all apparently panmictic, i.e. they did not show restricted gene flow such as observed among West African populations (interpreted as incipient speciation). Chromosomal inversion patterns of An. gambiae in southern Sudan showed characteristics of intergrading Savanna/Forest populations similar to those observed in comparable eco-climatic situations of West Africa. Anopheles arabiensis was polymorphic for inversion systems recorded in West Africa (2Ra, 2Rb, 2Rdl, 3Ra) and for a novel 2Rs polymorphism, overlapping with inversion systems 2Rb and 2Rd1. Samples carrying the 2Rs inversion were mostly from Khashm-el-Girba area in central-eastern Sudan. In the great majority of the samples all polymorphic inversions were found to be in Hardy-Weinberg equilibrium. Sudan populations of An. arabiensis should therefore be considered as generally panmictic. Anopheles arabiensis shows more inversion polymorphism in west than in east African populations. Sudan populations have more evident similarities with those from westwards than those from eastwards of the Great Rift Valley. The possible influence of the Rift on evolution of An. arabiensis is discussed.     

The Anopheles gambiae complex in Guinea Bissau

Three species of the Anopheles gambiae complex were identified in Guinea Bissau (West Africa) by chromosomal analysis. They were An. melas, An. arabiensis and An. gambiae s.s. An melas was observed in coastal and insular zones of the study area as well as in areas where the rivers are tidal and brackish and bordered by mangroves. For this reason, the species occurs also in inland riverine localities such as Farim and Bissorà. An. arabiensis apparently occurs only in low numbers in a very limited inland area during the dry season. An gambiae s.s. was observed nearly everywhere in the study area. In the samples of An. melas three inversion polymorphisms occurred: one on the chromosomal arm 2R (2Rn) and two on the arm 3R (3Rc and 3Re). It was observed that the frequencies of the inverted arrangements 2Rn and 3Re were significantly higher in the islands as compared to the continental sampling localities. The An. arabiensis sample was characterized by the presence of three inversion polymorphisms: 2Ra, 2Rb and 3Ra. A very high degree of polymorphism was shown by the An. gambiae s.s. samples. Four inversion polymorphisms were observed (three on chromosomal arm 2R and one on arm 2L) with very different frequencies of the alternative arrangements in different zones of the study area. From these data it seemed possible to split the species into three populations, each of them apparently linked with a peculiar ecological situation. The first population, characterized by high frequencies of 2Rd arrangement, is present on the coastal zones and in the islands; the second one is present in the northern inland areas particularly during the dry season and it is characterized by high frequencies of 2Rb and 2La arrangements. The third population is present only in the inland zones and it is characterized by high frequencies of 2Rjb, 2Rjd and 2Rjbd arrangements.     

Selective Introgression of Paracentric Inversions between Two Sibling Species of the Anopheles Gambiae Complex

The Anopheles gambiae complex includes the major vectors of malaria in sub-Saharan Africa where >80% of all world-wide cases occur. These mosquitoes are characterized by chromosomal inversions associated to the speciation process and to intraspecific ecological and behavioral flexibility. It has been postulated that introgressive hybridization has selectively transferred inversions on the second chromosome between A. gambiae and A. arabiensis, the two most important vectors of malaria. Here we directly test this hypothesis with laboratory experiments in which hybrid populations were established and the fate of chromosomal inversions were followed. Consistent with the hypothesis, ``foreign' X chromosomes were eliminated within two generations, while some ``foreign' second chromosomes persisted for the duration of the experiments and, judging from the excess of heterozygotes, established stable heterotic polymorphisms. Only those second chromosome inversions found naturally in the species could be introgressed.     

Observations on the Anopheles gambiae complex in the Senegal River Basin, West Africa

1. Three sibling species of mosquitoes of the Anopheles gambiae complex are found in the Senegal River Basin: An. melas Theobald, An. gambiae Giles and An. arabiensis Patton. 2. An. melas is restricted to the river delta and environs where saltwater breeding places are present. 3. An. gambiae and An. arabiensis are sympatric in the study area; An. arabiensis predominates in coastal zones where it breeds also during the dry season; An. gambiae predominates in inland areas where breeding is mostly restricted to the rainy season (July-September). 4. An. arabiensis is chromosomally polymorphic all over the study area, with much variation of inversion frequencies, particularly for the 2Ra arrangement. 5. An. gambiae is characterized by a very high frequency of the 2Rb-2La inversion arrangement which is typical of the Savanna chromosomal form.     

Exceptional diversity, maintenance of polymorphism, and recent directional selection on the APL1 malaria resistance genes of Anopheles gambiae

The three-gene APL1 locus encodes essential components of the mosquito immune defense against malaria parasites. APL1 was originally identified because it lies within a mapped QTL conferring the vector mosquito Anopheles gambiae natural resistance to the human malaria parasite, Plasmodium falciparum, and APL1 genes have subsequently been shown to be involved in defense against several species of Plasmodium. Here, we examine molecular population genetic variation at the APL1 gene cluster in spatially and temporally diverse West African collections of A. gambiae. The locus is extremely polymorphic, showing evidence of adaptive evolutionary maintenance of genetic variation. We hypothesize that this variability aids in defense against genetically diverse pathogens, including Plasmodium. Variation at APL1 is highly structured across geographic and temporal subpopulations. In particular, diversity is exceptionally high during the rainy season, when malaria transmission rates are at their peak. Much less allelic diversity is observed during the dry season when mosquito population sizes and malaria transmission rates are low. APL1 diversity is weakly stratified by the polymorphic 2La chromosomal inversion but is very strongly subdivided between the M and S "molecular forms." We find evidence that a recent selective sweep has occurred at the APL1 locus in M form mosquitoes only. The independently reported observation of a similar M-form restricted sweep at the Tep1 locus, whose product physically interacts with APL1C, suggests that epistatic selection may act on these two loci causing them to sweep coordinately.     

Dieldrin resistance in the malaria vector Anopheles gambiae in Ghana

Anopheles gambiae Giles s.s. (Diptera: Culicidae) is one of the principal vectors of malaria in the Ashanti region of central Ghana. High levels of resistance to dieldrin were recorded in a wild-caught sample from Obuasi (south of Kumasi) as well as a laboratory colony established using material from the wild population. Cytogenetic analysis of wild-caught and laboratory samples revealed chromosomal polymorphism for inversions 2La and 2Rb. Although inversion 2La has previously been shown to be associated with dieldrin resistance in certain other laboratory strains originating from West Africa, there was no obvious association between inversion karyotype assortment and the resistance phenotype in the Obuasi population. In addition, polymerase chain reaction analysis indicated the presence of the alanine296 to glycine mutation in the GABA (gamma amino-butyric acid) receptor (which has been mapped to a chromosomal position within inversion 2La). This mutation has previously been shown to be associated with dieldrin resistance in the same An. gambiae laboratory strains of West African origin. Our data show only a weak association between the dieldrin resistance phenotype and the presence of this mutation, suggesting that another dieldrin resistance mechanism is operational in the Obuasi population. Biochemical and synergist exposure assays suggest a metabolic component, probably mediated by monooxygenase P450 enzymes. We conclude that dieldrin resistance in the An. gambiae population of the Obuasi region occurs at a high level - most likely in the absence of selection - and that control of the resistance phenotype is polyfactorial and must include components other than mutations in the GABA receptor locus.     

Gene flow between chromosomal forms of the malaria vector Anopheles funestus in Cameroon, Central Africa, and its relevance in malaria fighting

Knowledge of population structure in a major vector species is fundamental to an understanding of malaria epidemiology and becomes crucial in the context of genetic control strategies that are being developed. Despite its epidemiological importance, the major African malaria vector Anopheles funestus has received far less attention than members of the Anopheles gambiae complex. Previous chromosomal data have shown a high degree of structuring within populations from West Africa and have led to the characterization of two chromosomal forms, "Kiribina" and "Folonzo." In Central Africa, few data were available. We thus undertook assessment of genetic structure of An. funestus populations from Cameroon using chromosomal inversions and microsatellite markers. Microsatellite markers revealed no particular departure from panmixia within each local population and a genetic structure consistent with isolation by distance. However, cytogenetic studies demonstrated high levels of chromosomal heterogeneity, both within and between populations. Distribution of chromosomal inversions was not random and a cline of frequency was observed, according to ecotypic conditions. Strong deficiency of heterokaryotypes was found in certain localities in the transition area, indicating a subdivision of An. funestus in chromosomal forms. An. funestus microsatellite genetic markers located within the breakpoints of inversions are not differentiated in populations, whereas in An. gambiae inversions can affect gene flow at marker loci. These results are relevant to strategies for control of malaria by introduction of transgenes into populations of vectors.     

Anopheles bwambae sp.n., a malaria vector in the Semliki Valley, Uganda, and its relationships with other sibling species of the An.gambiae complex (Diptera: Culicidae)

Abstract. The name Anopheles bwambae is proposed for the taxon previously called species D of the An.gambiae complex. This sibling species is known only from the vicinity of Buranga hot springs in Bwamba County, Toro District, Uganda, where it breeds in brackish water from geothermal springs together with other halophilic mosquitoes (Ae.albocephalus, Ae.natwnius, Cx.tenagius). An.bwambae adults inhabit the Semliki Forest, where the natural hosts are unknown, but they also enter houses in nearby villages and bite human beings avidly. Plasmodium sporozoites and developing larvae of Wuchereria bancrofti were found in An.bwambae females, so it is assumed that this species of vector contributes locally to transmission of human malaria and filariasis (together with An.funestus, An.gambiae and An.arabiensis). The specific diagnosis of An.bwambae depends upon (i) the presence of fixed inversion 3La and polymorphic inversions 2R1 and 3Rb of the karyotype, as interpreted from ovarian polytene chromosomes, (ii) a fast fixed allozyme of super-oxide dismutase (Sod-105), (iii) female palpi with broader pale apical band and narrower dark sub-apical band than for other sibling species. Cytotaxonomically An.bwambae has the same X-chromosome banding pattern as An.quadriannulatus and An.melas; the latter also has autosomal inversion 3La and therefore appears to be the sister-species of An.bwambae. An identification key is given to the six named sibling species of the An.gambiae complex.     

An integrated chromosome map of microsatellite markers and inversion breakpoints for an Asian malaria mosquito, Anopheles stephensi

Anopheles stephensi is one of the major vectors of malaria in the Middle East and Indo-Pakistan subcontinent. Understanding the population genetic structure of malaria mosquitoes is important for developing adequate and successful vector control strategies. Commonly used markers for inferring anopheline taxonomic and population status include microsatellites and chromosomal inversions. Knowledge about chromosomal locations of microsatellite markers with respect to polymorphic inversions could be useful for better understanding a genetic structure of natural populations. However, fragments with microsatellites used in population genetic studies are usually too short for successful labeling and hybridization with chromosomes. We designed new primers for amplification of microsatellite loci identified in the A. stephensi genome sequenced with next-generation technologies. Twelve microsatellites were mapped to polytene chromosomes from ovarian nurse cells of A. stephensi using fluorescent in situ hybridization. All microsatellites hybridized to unique locations on autosomes, and 7 of them localized to the largest arm 2R. Ten microsatellites were mapped inside the previously described polymorphic chromosomal inversions, including 4 loci located inside the widespread inversion 2Rb. We analyzed microsatellite-based population genetic data available for A. stephensi in light of our mapping results. This study demonstrates that the chromosomal position of microsatellites may affect estimates of population genetic parameters and highlights the importance of developing physical maps for nonmodel organisms.     

Identifying targets of selection in mosaic genomes with machine learning: applications in Anopheles gambiae for detecting sites within locally adapted chromosomal inversions

Chromosomal inversions are important structural changes that may facilitate divergent selection when they capture co‐adaptive loci in the face of gene flow. However, identifying selection targets within inversions can be challenging. The high degrees of differentiation between heterokaryotypes, as well as the differences in demographic histories of collinear regions compared with inverted ones, reduce the power of traditional outlier analyses for detecting selected loci. Here, we develop a new approach that uses discriminant functions informed from inversion‐specific expectations to classify loci that are under selection (or drift). Analysis of RAD sequencing data we collected in a classic dipteran species with polymorphic inversion clines—Anopheles gambiae, a malaria vector species from sub‐Saharan Africa—demonstrates the benefits of the approach compared with traditional outlier analyses. We focus specifically on two polymorphic inversions, the 2La and 2Rb arrangements that predominate in dry habitats and the 2L+^a and 2R+^b arrangements in wet habitats, which contrast with the minimal geographic structure of SNPs from collinear regions. With our approach, we identify two strongly selected regions within 2La associated with dry habitat. Moreover, we also show that the prevalence of selection is greater in the arrangement 2L+^a that is associated with wet habitat (unlike presumed importance of selective divergence associated with the shift of the mosquitoes into dry habitats). We discuss the implications of these results with respect to studies of rapid adaptation in these malaria vectors, and in particular, the insights our newly developed approach offers for identifying not only potential targets of selection, but also the population that has undergone adaptive change.     

Association mapping desiccation resistance within chromosomal inversions in the African malaria vector Anopheles gambiae

Inversion polymorphisms are responsible for many ecologically important phenotypes and are often found under balancing selection. However, the same features that ensure their large role in local adaptation—especially reduced recombination between alternate arrangements—mean that uncovering the precise loci within inversions that control these phenotypes is unachievable using standard mapping approaches. Here, we take advantage of long‐term balancing selection on a pair of inversions in the mosquito Anopheles gambiae to map desiccation tolerance via pool‐GWAS. Two polymorphic inversions on chromosome 2 of this species (denoted 2La and 2Rb) are associated with arid and hot conditions in Africa and are maintained in spatially and temporally heterogeneous environments. After measuring thousands of wild‐caught individuals for survival under desiccation stress, we used phenotypically extreme individuals homozygous for alternative arrangements at the 2La inversion to construct pools for whole‐genome sequencing. Genomewide association mapping using these pools revealed dozens of significant SNPs within both 2La and 2Rb, many of which neighboured genes controlling ion channels or related functions. Our results point to the promise of similar approaches in systems with inversions maintained by balancing selection and provide a list of candidate genes underlying the specific phenotypes controlled by the two inversions studied here.     

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.     

Seasonal climate effects anemotaxis in newly emerged adult Anopheles gambiae Giles in Mali, West Africa

The direction and magnitude of movement by the malaria vector Anopheles gambiae Giles has been of great interest to medical entomologists for over 70 years. This direction of movement is likely to be affected by many factors, from environmental conditions and stage of life history of the mosquito to the existence of attractants in the vicinity. We report here the direction of movement of newly emerged An. gambiae in nature, around the village of Donéguébougou, Mali. We assessed the direction of movement for individual mosquitoes by placing them in a novel enclosure with exit traps oriented in the direction of the cardinal and intermediate points of the compass. We consistently found predominantly Southward directions of movement during 2009 and 2010, with an additional Eastward component during the dry season and a Westward one during the wet season. Our data indicate that wind has an important effect on the direction of movement, but that this effect varied by season: Average directions of movement were downwind during the dry season and upwind during the wet season. A switch in anemotactic response suggests that the direction of movement of An. gambiae relative to the wind immediately after emergence under varying conditions of humidity should be further investigated under controlled conditions.     

Population structure, speciation, and introgression in the Anopheles gambiae complex

We review here what is known about the population structure and evolutionary dynamics of members of the Anopheles gambiae complex with emphasis on the situation in West Africa. First, the importance of the 2nd chromosome inversion polymorphism is demonstrated especially in adaptation to levels of aridity, a major environmental variable in Africa. This affects the distribution of karyotypes on both a macro- and micro-geographic scale as well as temporally. Such differentiation leads to karyotypes being differentially effective transmitters of malaria and differentially susceptible to indoor residual spraying of insecticides. Second, we review the evidence that cryptic taxa, especially in An. gambiae s.s., exist. This observation stems from both karyotype studies and molecular studies. It is abundantly clear that West African populations of An. gambiae s.s. are often not panmictic units, with premating factors evidently acting to maintain distinct genetic forms. Third, we review phylogenetic studies that have revealed the presence of introgression between the two most important vectors, An. gambiae and An. arabiensis. This is most evident for the 2nd chromosome inversions. This interpretation of phylogenetic data is consistent with a direct laboratory study indicating inversions in this chromosome are stably maintained in back-crossed populations. All of this information has led to the view that members of the An. gambiae complex are highly variable with an abundance of adaptive genetic variation. This presents a significant challenge to vector control programs designed to reduce malaria in sub-Saharan Africa.     

Are chromosomal inversions induced by transposable elements? A paradigm from the malaria mosquito Anopheles gambiae

Chromosomal rearrangements abound in nature and can be studied in detail in organisms with polytene chromosomes. In Drosophila and in Anopheline mosquitoes most speciation processes seem to be associated with the establishment of chromosomal rearrangements, particularly of paracentric inversions. It is not known what triggers inversions in natural populations. In the laboratory inversions are commonly generated by X-rays, mutagens or after the activity of certain transposable elements (TEs). The Anopheles gambiae complex is comprised of six sibling species, each one characterized by the presence of fixed paracentric inversions on their chromosomes. Two of these, An. gambiae s.s. and An. arabiensis, are the most important vectors of human malaria and are structured into sub-populations, each carrying a characteristic set of polymorphic chromosomal inversions. We have cloned the breakpoints of the naturally occurring polymorphic inversion In(2R)d' of An. arabiensis. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a transposable element that we called Odysseus. Characteristics of Odysseus' terminal region and its cytological distribution in different strains as well as within the same strain indicate that Odysseus is an actively transposing element. The presence of Odysseus at the junction of the naturally occurring inversion In(2R)d' suggests that the inversion may be the result of the TEs activity. Cytological evidence from Drosophila melanogaster has also implicated the hobo transposable element in the generation of certain Hawaiian endemic inversions. This picture supports the hypothesis of the important role of TEs in generating natural inversions.     

Molecular identification of chromosomal forms of Anopheles gambiae sensu stricto

The Afrotropical malaria vectors Anopheles gambiae sensu stricto and An. arabiensis, responsible for more than 3/4 of the world Plasmodium falciparum inoculations, are members of the Anopheles gambiae complex, a group consisting of six sibling species. The nominal species (An. gambiae s.s.) is by far the most anthropophilic vector and its adaptation to man and his environment involves further ongoing speciation processes. This fact is shown by the existence of a number of incipient taxonomic units characterised by different chromosomal arrangements derived from the presence of polymorphic paracentric inversions. This speciation process is centered in West Africa, where five so-called 'chromosomal forms' have been described, designated with a non-Linnean nomenclature: Forest, Bissau, Savanna, Bamako, and Mopti. Studies on the distribution and the ecology of these incipient species have highlighted specific adaptations to eco-ethological parameters, which might reflect on their relative efficiency as malaria vectors. Cytogenetic analysis, in spite of some inherent difficulties, has proved to be a powerful tool for the identification of An. gambiae sibling species and the individual chromosomal forms. Yet, modern molecular tools are now available, providing alternative faster low-cost technologies, and we discuss here their relative merits.     

Linear and spatial organization of polytene chromosomes of the African malaria mosquito Anopheles funestus

Anopheles funestus Giles is one of the major malaria vectors in Africa, but little is known about its genetics. Lack of a cytogenetic map characterized by regions has hindered the progress of genetic research with this important species. This study developed a cytogenetic map of An. funestus using ovarian nurse cell polytene chromosomes. We demonstrate an important application with the cytogenetic map for characterizing various chromosomal inversions for specimens collected from coastal Kenya. The linear and spatial organization of An. funestus polytene chromosomes was compared with the best-studied malaria mosquito, An. gambiae Giles. Comparisons of chromosome morphology between the two species have revealed that the most extensive chromosomal rearrangement occurs in pericentromeric heterochromatin of autosomes. Differences in pericentromeric heterochromatin types correlate with nuclear organization differences between An. funestus and An. gambiae. Attachments of chromosomes to the nuclear envelope strongly depend on the presence of diffusive beta-heterochromatin. Thus, An. funestus and An. gambiae exhibit species-specific characteristics in chromosome-linear and -spatial organizations.