- Open Access
Seasonal distribution of Anopheles funestus chromosomal forms from Burkina Faso
Malaria Journalvolume 8, Article number: 239 (2009)
Previous studies of Anopheles funestus chromosomal inversion polymorphisms in Burkina Faso showed large departures from Hardy-Weinberg equilibrium and linkage disequilibrium among inversions located on different chromosomes, implying the existence of two taxonomic units ("chromosomal forms") with limited genetic flow. One chromosomal form, named Folonzo, is highly polymorphic for alternative rearrangements of 3Ra, 3Rb, 2Ra, and 3La; the other, Kiribina, is predominantly characterized by the standard arrangement of these inversions. To investigate the temporal distribution of these chromosomal forms, further collections were carried out in two villages near Ouagadougou where they are found in sympatry.
Chromosomal karyotypes were determined from indoor-resting, half-gravid females sampled within and across six breeding seasons, from December 1998 to April 2007.
As expected, the pattern of chromosomal polymorphism in An. funestus was consistent with assortatively mating Folonzo and Kiribina forms. When samples were assigned to each chromosomal form, their relative abundance varied within successive breeding seasons in a repeating pattern of temporal variability. Relative abundance of the Folonzo form was correlated with climatic variables related to temperature and rainfall.
The relative abundance of Folonzo and Kiribina forms of An. funestus likely reflects different larval ecologies that are linked to varying climatic conditions. Further analysis of the bionomics of these vectors is recommended in light of its relevance to vector control.
Malaria remains a major health concern in Africa today. This situation is due to the presence of three efficient vectors in subgenus Cellia--Anopheles gambiae, its sibling species Anopheles arabiensis and Anopheles funestus. These species co-occur geographically across most of sub-Saharan Africa and can inhabit the same villages, shelter in the same houses, and blood-feed on the same individuals. Vector control efforts are complicated not only by the spread of resistance to insecticides, but also by behavioural and genetic variation within and between species. Among the challenges of efficient vector control are cryptic barriers to gene flow among populations of the same species that can arise as a consequence of differential adaptations to heterogeneities in the environment. This challenge applies not only to the introduction and spread of "refractoriness genes" in wild mosquito populations, but also to the use of bed nets whose efficacy could be lessened by exophilic and exophagic vector sub-populations . Sustained vector control efforts across Africa will benefit from knowledge of the extent of genetic diversity in natural populations, how it is distributed in time and space, and how it is generated and maintained among vector populations [2–4].
An. funestus abounds during the dry season and is less rain-dependant than An. gambiae s.l., owing to the tendency to breed in permanent or semi-permanent swamps or pools . As a result, it is considered as a vector that bridges malaria transmission across the dry season. However, it has long been noted that An. funestus is a highly polymorphic species whose populations could be structured. Evidence of genetic heterogeneities within this mosquito were revealed during early cytogenetic investigations . Polymorphic inversions are found in several populations in East , South [6, 8], Central [3, 9], and West Africa [10–12]. Some of the inversion systems are specific to each locality. Echoing the discovery of complexities in other isomorphic anopheline species, cytogenetic analysis of An. funestus in Burkina Faso revealed large departures from Hardy-Weinberg equilibrium and linkage disequilibrium among inversions located on different chromosomes. These results led to the recognition of two chromosomal forms, namely "Folonzo" and "Kiribina," characterized by contrasting degrees of chromosomal polymorphism and limited amounts of gene flow, considered indicative of an incipient speciation process . This hypothesis was reinforced in a detailed cytogenetic analysis during three consecutive seasons in a pair of villages less than 2 km apart near Ouagadougou . More recently, further evidence for genetic differentiation between the two chromosomal forms in Burkina Faso was found using molecular markers: 16 microsatellite loci dispersed over the chromosome complement, and the ND5 region of mtDNA . The geographic extent of this diversification in An. funestus is uncertain, as studies of the two chromosomal forms in Cameroon and Senegal detected no evidence of genetic differentiation beyond that accounted for by geographical isolation [3, 15].
The spatio-temporal pattern of these diverging entities and the underlying factors governing the pattern remain unclear. Ultimately, assessment of both spatial and temporal components of the genetic structure of a species is necessary to identify factors that influence the genetic variability and relationships among its populations. Monitoring populations over time, the objective of the present study, allows determination of the degree of stability of spatial structuring within a population or local geographical area. Toward this end, a temporal analysis of An. funestus chromosomal forms was performed during seven breeding seasons in villages where the two forms are sympatric.
Mosquito collections were carried out in two rural villages near Ouagadougou, Burkina Faso: Koubri (12°11'54"N; 1°23'43"W) and Kuiti (12°11'36"N; 1°23'11"W). The two villages lie in the West African sudan savannah vegetation belt, and are located ~2 km apart on opposite margins of a swamp [for detailed map, see ]. The area receives 750 mm of rain per year, mainly falling during the rainy season from July until September. The breeding season of An. funestus starts at the end of the rainy season (September), extends throughout the cool dry season (October-February), and ends in April, mid-way through the hot dry season (March-May).
For detailed analysis of chromosomal forms within breeding seasons, indoor resting collections were made in three consecutive years from August 1999 to March 2002 (see Figure 1). Mosquitoes were collected three times per week manually or by spray-catch. Further details of these collections are presented in references [13, 14]. Based on preliminary karyotype analysis of these collections, further collections were performed during three additional years, at the anticipated peak period of chromosomal polymorphism beginning in October and continuing through January or beyond.
Adult An. funestus group mosquitoes were sorted morphologically in the field [5, 16]. Half-gravid females were immediately dissected; ovaries were preserved individually in 1.5 ml microcentrifuge tubes using Carnoy's fixative and the carcass was placed in a correspondingly labelled tube with desiccant. After fixation, polytene chromosomes were squashed . The slides were examined under a phase-contrast microscope, and chromosomal arrangements were scored with reference to the cytogenetic map . Assignment to chromosomal form was based on the algorithm of Costantini et al  modified by Guelbeogo et al .
Climatic variables (monthly temperature, relative humidity, and rainfall) were provided by the National Meteorological Agency through their meteorological station located in the village of Koubri.
As previous studies showed genetic homogeneity within chromosomal forms collected from Koubri and Kuiti [13, 14], analysis was performed on pooled samples from these two adjacent villages. The Spearman Rank correlation coefficient was used to test the association between the relative abundance (i.e. relative frequency--hereafter used synonymously) of Folonzo and several climatic variables: total annual and monthly rainfall, monthly relative humidity (maximum and minimum), and monthly temperature (minimum, maximum and mean). Correlation was tested between the pooled relative abundance of the An. funestus Folonzo form in Oct-Dec of each year with respect to the corresponding climatic variables during the preceding (Jan-Sep) or concurrent (Oct-Dec) months. Where a significant correlation was found, the relationship between those variables and Folonzo relative abundance was parameterized by logistic regression. Analyses were run in R v. 2.9.0 .
The study protocol was reviewed and approved by the institutional health ethical review board of Burkina Faso (code N° 2007-034). Participants were explained the study procedures, benefits and risks and their consent was obtained before collecting mosquito in their settlements.
A total of 6,582 indoor-resting An. funestus adults were successfully scored chromosomally. Of these, 1,959 and 4,623 were assigned to Folonzo and Kiribina forms, respectively.
In the three consecutive breeding seasons for which monthly samples were taken, relative frequencies of the chromosomal forms cycled in a repeatable pattern (Figure 1). The Folonzo form achieved its greatest relative frequencies in October-December. Although other seasonal collections were less comprehensive, the chromosomal data were consistent with this same pattern of Folonzo peaking in relative abundance following the rains.
Beyond this general trend, the average relative frequency of the two chromosomal forms varied between breeding seasons. These temporal patterns of abundance suggest the influence of climatic variables. The possible relationship between Folonzo relative frequency during the peak months of abundance (Oct-Dec) and rainfall, relative humidity and temperature was explored with reference to months preceding the breeding season in question. A significant correlation was found with four monthly variables (Table 1). During January, maximum relative humidity and temperature were negatively correlated with Folonzo relative abundance in the following breeding season. Similarly in March, mean and minimum temperatures were positively correlated with Folonzo relative frequency afterward. Overall, annual rainfall amount was positively correlated with Folonzo relative abundance. Figure 2 shows the relationships between these five variables and the Folonzo form frequency as described by logistic regression lines. The corresponding equations are provided in Table 2. It is acknowledged that some of these significant associations might be spurious because the climatic variables used are not independent, and a large number of tests were performed. However, in adopting the approach recommended by Moran  for ecological research, a correction for multiple statistical tests was not applied at this exploratory stage. The fact that five tests gave fairly low probability values suggests that a general pattern of association between climate and Folonzo relative abundance is likely to be valid, but further studies will be required to evaluate the relative importance of these and other variables to the observed pattern.
The three main components of the malaria vectorial system in tropical Africa have been observed to cycle in relative abundance according to seasonal ecological conditions , in a pattern often described as a relay due to the temporal succession of An. gambiae followed by An. arabiensis and finally An. funestus. For the first two mosquito vectors, the positive relationship between rainfall and abundance is unsurprising, because these species typically breed in temporary rain-dependent pools and puddles. Because An. funestus characteristically breeds in large permanent or semi-permanent pools (e.g., swamps and irrigated rice fields), rain-dependent fluctuations in its abundance may seem less obvious. However, emergent vegetation--particularly around the margins of breeding sites--is considered one of the more important ecological features for this species . Especially in the dry savannah where rainfall is restricted to one season, both the physical extent of the site as well as the growth of vegetation depends upon rainfall and the level of the water table. Accordingly, potential An. funestus larval breeding sites are not uniformly distributed during the year.
Available breeding sites for An. funestus in the immediate study area of Kuiti and Koubri include both a permanent swamp and rice fields. There is considerable modification of the swamp at the end of rainy season because decreasing water levels lead to the appearance and disappearance of different types of breeding sites. In addition, human activities surrounding rice cultivation--the flooding of the fields in June/July and their harvest in October/November--destroy and create new breeding sites that may favour one or the other of the alternative chromosomal forms. Because the mosquitoes analyzed during this study were collected as adults, there is no direct evidence linking the increased relative abundance of the Folonzo form to features of these breeding sites. However, by analogy to habitat segregation in the An. gambiae complex [21–24], it is possible to hypothesize that niche partitioning exists between the chromosomal variants of An. funestus. Consistent with this idea, Costantini et al  reported high frequencies of Folonzo in villages associated with water reservoirs containing natural emerging and floating vegetation, whereas high frequencies of Kiribina were documented in villages with large scale rice crop areas. This suggests that Kiribina prefers breeding in rice fields or out-competes Folonzo in these artificial, anthropogenic habitats; the converse may be true of Folonzo in sites containing natural vegetation, such as swamps. Testing of this hypothesis crucially depends upon the discovery and application of molecular markers that can distinguish Folonzo and Kiribina forms, for two reasons. From a practical standpoint, only the polytene chromosomes extracted from half-gravid adult females can be used in the recognition of the An. funestus chromosomal forms presently , posing a serious impediment if not roadblock to research on their larval ecology. Of potential theoretical importance is the correspondence between chromosomally recognized variants and true assortatively mating and biologically relevant reproductive units. As is apparent in the case of An. gambiae, such correspondence is far from perfect, reflecting the fact that different reproductive units share chromosomal variants to some extent [23–26]. If a similar situation exists within An. funestus as hypothesized, the discovery of molecular markers better reflecting reproductive discontinuities will be an important prerequisite to further research defining the ecological mechanisms driving diversification in this lineage. It is hoped that the recently approved whole genome sequencing of several Anopheles species including An. funestus  will offer the requisite tools to help overcome this problem.
Coluzzi [2, 28] has emphasized the significance of such a genetically heterogeneous malaria vectorial system to malaria epidemiology and control in Africa, because it results in very high and stable disease transmission. The operational value of distinguishing between isomorphic vector and non-vector species can be taken for granted, but Coluzzi [2, 28] also stressed the less obvious relevance of distinctions between reproductively isolated siblings or incipient species presumed to share high vectorial potential. The core of his argument is that the stable co-existence of two taxa implies some form of biological divergence, the absence of which would lead to costly and ultimately unsustainable competition. Biological divergence could be reflected in such epidemiologically relevant factors as vectorial capacity and response to vector control strategies. In fact, it is conceivable that the Folonzo and Kiribina forms do differ in their behavioural responses to environmental heterogeneities. In a remarkable parallel with An. gambiae [25, 29, 30], inversion frequencies in An. funestus populations from Cameroon vary clinally from mesic forest regions in the south where inverted arrangements are fixed to xeric regions in the north where the standard karyotype is fixed . These findings are in agreement with the positive association between Folonzo relative frequency and total rainfall that was found in the present study. This suggests that Folonzo would be better adapted to environments with a higher degree of water vapour saturation and Kiribina would be better adapted to periods with a higher saturation deficit, which could impact epidemiologically important behaviours such as indoor resting and biting. This pattern of chromosomal variation is also consistent with the present findings in an area of sympatry, where the Folonzo form breeds preferentially during the months following the end of the rainy season. Answering Coluzzi's  call for vector analysis in the case of An. funestus thus will require a wide range of genotypic and phenotypic information preferably based on parallel adult and larvae sampling.
Molineaux L, Grammicia G: The Garki Project. Research on the Epidemiology and Control of Malaria in the Sudan Savanna of West Africa. 1980, Geneva, Switzerland: World Health Organization
Coluzzi M: Malaria vector analysis and control. Parasitol Today. 1992, 8: 113-118. 10.1016/0169-4758(92)90277-9.
Cohuet A, Dia I, Simard F, Raymond M, Rousset F, Antonio-Nkondjio C, Awono-Ambene PH, Wondji CS, Fontenille D: Gene flow between chromosomal forms of the malaria vector Anopheles funestus in Cameroon, Central Africa, and its relevance in malaria fighting. Genetics. 2005, 169: 301-311. 10.1534/genetics.103.025031.
Tripet F, Dolo G, Lanzaro GC: Multilevel analyses of genetic differentiation in Anopheles gambiae s.s. reveal patterns of gene flow important for malaria-fighting mosquito projects. Genetics. 2005, 169: 313-324. 10.1534/genetics.104.026534.
Gillies MT, De Meillon B: The Anophelinae of Africa South of the Sahara. 1968, Johannesburg: South African Institute for Medical Research, 2
Green CA, Hunt RH: Interpretation of variation in ovarian polytene chromosomes of Anopheles funestus Giles, A. parensis Gillies, and A. aruni?. Genetica. 1980, 51: 187-195. 10.1007/BF00121610.
Kamau L, Hunt R, Coetzee M: Analysis of the population structure of Anopheles funestus (Diptera: Culicidae) from western and coastal Kenya using paracentric chromosomal inversion frequencies. J Med Entomol. 2002, 39: 78-83. 10.1603/0022-2585-39.1.78.
Boccolini D, Carrara GC, Dia I, Fortes F, Cani PJ, Costantini C: Chromosomal differentiation of Anopheles funestus from Luanda and Huambo Provinces, western and central Angola. Am J Trop Med Hyg. 2005, 73: 1071-1076.
Dia I, Boccolini D, Antonio-Nkondjio C, Costantini C, Fontenille D: Chromosomal inversion polymorphism of Anopheles funestus from forest villages of South Cameroon. Parassitologia. 2000, 42: 227-229.
Costantini C, Sagnon NF, Ilboudo-Sanogo E, Coluzzi M, Boccolini D: Chromosomal and bionomic heterogeneities suggest incipient speciation in Anopheles funestus from Burkina Faso. Parassitologia. 1999, 41: 595-611.
Dia I, Lochouarn L, Boccolini D, Costantini C, Fontenille D: Spatial and temporal variations of the chromosomal inversion polymorphism of Anopheles funestus in Senegal. Parasite. 2000, 7: 179-184.
Lochouarn L, Dia I, Boccolini D, Coluzzi M, Fontenille D: Bionomical and cytogenetic heterogeneities of Anopheles funestus in Senegal. Trans R Soc Trop Med Hyg. 1998, 92: 607-612. 10.1016/S0035-9203(98)90782-9.
Guelbeogo WM, Grushko O, Boccolini D, Ouedraogo PA, Besansky NJ, Sagnon NF, Costantini C: Chromosomal evidence of incipient speciation in the Afrotropical malaria mosquito Anopheles funestus. Med Vet Entomol. 2005, 19: 458-469. 10.1111/j.1365-2915.2005.00595.x.
Michel AP, Guelbeogo WM, Grushko O, Schemerhorn BJ, Kern M, Willard MB, Sagnon N, Costantini C, Besansky NJ: Molecular differentiation between chromosomally defined incipient species of Anopheles funestus. Insect Mol Biol. 2005, 14: 375-387. 10.1111/j.1365-2583.2005.00568.x.
Cohuet A, Dia I, Simard F, Raymond M, Fontenille D: Population structure of the malaria vector Anopheles funestus in Senegal based on microsatellite and cytogenetic data. Insect Mol Biol. 2004, 13: 251-258. 10.1111/j.0962-1075.2004.00482.x.
Gillies MT, Coetzee M: A Supplement to the Anophelinae of Africa South of the Sahara. 1987, Johannesburg: The South African Institute for Medical Research
della Torre A: Polytene chromosome preparation from anopheline mosquitoes. Molecular Biology of Disease Vectors: A Methods Manual. Edited by: Crampton JM, Beard CB, Louis C. 1997, London: Chapman & Hall, 329-336.
Sharakhov IV, Sharakhova MV, Mbogo CM, Koekemoer LL, Yan G: Linear and spatial organization of polytene chromosomes of the African malaria mosquito Anopheles funestus. Genetics. 2001, 159: 211-218.
R_Development_Core_Team: R: A language and environment for statistical computing. 2009, R Foundation for Statistical Computing, Vienna, Austria, [http://www.R-project.org]
Moran MD: Arguments for rejecting the sequential Bonferroni in ecological studies. Oikos. 2003, 100: 403-405. 10.1034/j.1600-0706.2003.12010.x.
Coluzzi M, Petrarca V, DiDeco MA: Chromosomal inversion intergradation and incipient speciation in Anopheles gambiae. Boll Zool. 1985, 52: 45-63.
Coluzzi M, Sabatini A, Della Torre A, Di Deco MA, Petrarca V: A polytene chromosome analysis of the Anopheles gambiae species complex. Science. 2002, 298: 1415-1418. 10.1126/science.1077769.
Simard F, Ayala D, Kamdem GC, Etouna J, Ose K, Fotsing J-M, Fontenille D, Besansky NJ, Costantini C: Ecological niche partitioning between the M and S molecular forms of Anopheles gambiae in Cameroon: the ecological side of speciation. BMC Ecol. 2009, 9: 17-10.1186/1472-6785-9-17.
Costantini C, Ayala D, Guelbeogo WM, Pombi M, Some CY, Bassole IHN, Ose K, Fotsing J-M, Sagnon NF, Fontenille D: Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae. BMC Ecol. 2009, 9: 16-10.1186/1472-6785-9-16.
Toure YT, Petrarca V, Traore SF, Coulibaly A, Maiga HM, Sankare O, Sow M, DiDeco MA, Coluzzi M: The distribution and inversion polymorphism of chromosomally recognized taxa of the Anopheles gambiae complex in Mali, West Africa. Parassitologia. 1998, 40: 477-511.
Della Torre A, Tu Z, Petrarca V: On the distribution and genetic differentiation of Anopheles gambiae s.s. molecular forms. Insect Biochem Mol Biol. 2005, 35: 755-769. 10.1016/j.ibmb.2005.02.006.
Besansky NJ: Genome Analysis Of Vectorial Capacity In Major Anopheles Vectors Of Malaria Parasites. 2008, Eukaryotic Pathogen and Disease Vector Sequencing Project, NIAID and NHGRI
Coluzzi M: Heterogeneities of the malaria vectorial system in tropical Africa and their significance in malaria epidemiology and control. Bull World Health Organ. 1984, 62: 107-113.
Toure YT, Petrarca V, Traore SF, Coulibaly A, Maiga HM, Sankare O, Sow M, Di Deco MA, Coluzzi M: Ecological genetic studies in the chromosomal form Mopti of Anopheles gambiae s.str. in Mali, West Africa. Genetica. 1994, 94: 213-223. 10.1007/BF01443435.
Coluzzi M, Sabatini A, Petrarca V, Di Deco MA: Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979, 73: 483-497. 10.1016/0035-9203(79)90036-1.
We thank the inhabitants of Koubri and Kuiti for their collaboration, and the entomological team of CNRFP for their efforts during the course of this study. This research was financially supported by grant R01 AI48842 from the US National Institutes of Health to NJB, re-entry grant A 41413 to GWM from WHO/TDR; and 50090 from the MIM/TDR.
The authors declare that they have no competing interests.
CC, NS and NB conceived the study. WG performed the field collection. NS provided logistical support throughout the study. WG, OG, DB and MAY performed the karyotyping. WG and CC analysed results. WG, NB and CC wrote the manuscript. All authors read and approved the final manuscript.