- Open Access
Dynamics of transmission of Plasmodium falciparum by Anopheles arabiensis and the molecular forms M and S of Anopheles gambiae in Dielmo, Senegal
Malaria Journalvolume 7, Article number: 136 (2008)
The adaptation of Anopheles gambiae to humans and its environment involves an ongoing speciation process that can be best demonstrated by the existence of various chromosomal forms adapted to different environments and of two molecular forms known as incipient taxonomic units.
The aim of this study was to compare the epidemiologic role of Anopheles arabiens is and the molecular forms M and S of Anopheles gambiae in the transmission of Plasmodium in a rural areas of southern Senegal, Dielmo. The sampling of mosquitoes was carried out monthly between July and December 2004, during the rainy season, by human volunteers and pyrethrum spray catches.
Anopheles arabiensis, An. gambiae M and S forms coexisted during the rainy season with a predominance of the M form in September and the peak of density being observed in August for the S form. Similar parity rates were observed in An. arabiensis [70.9%] (n = 86), An. gambiae M form [68.7%] (n = 64) and An. gambiae S form [81.1%] (n = 156). The circumsporozoite protein (CSP) rates were 2.82% (n = 177), 3.17% (n = 315) and 3.45% (n = 405), with the mean anthropophilic rates being 71.4% (n = 14), 86.3% (n = 22) and 91.6% (n = 24) respectively for An. arabiensis and An. gambiae M and S forms. No significant difference was observed either in host preference or in Plasmodium falciparum infection rates between sympatric M and S populations.
No difference was observed either in host preference or in Plasmodium falciparum infection rates between sympatric M and S populations, but they present different dynamics of population. These variations are probably attributable to different breeding conditions.
The Anopheles gambiae complex consists of at least seven species among which Anopheles gambiae s.s. is one of the most anthropophilic malaria vectors in Africa . The adaptation of An. gambiae to humans and its environment involves an ongoing speciation process that can be best demonstrated by the existence of a number of incipient taxonomic units, characterized by the presence of paracentric inversions leading to different chromosomal arrangements . This speciation process is primarily observed in West Africa, where five chromosomal forms of An. gambiae s.s. have been described and designated with a non-Linnean nomenclature: bamako, bissau, forest, mopti and savanna [3, 4].
During the last few years, several research teams have settled on a molecular approach to address speciation in An. gambiae s.s.. Various degrees of gene flow restriction were demonstrated between chromosomal forms, with strong hybrid heterokaryotype deficits in the areas of sympatry. Analysis of the rDNA intergenic spacers, located on the X-chromosome, revealed fixed sequence differences between sympatric and synchronous savanna/bamako and mopti populations in Mali, Burkina Faso and Cameroon [2, 5, 6]. To provide more insight into their taxonomic status, recent efforts have focused on the pattern of variation observed with molecular markers. This revealed the existence of two genetic variants referred to as the molecular M and S forms [7, 8].
However, whatever the geographical region, it has been clearly demonstrated that the gene flow between M and S forms is very limited, revealing a current speciation phenomenon. The genetic characteristics of these forms and their known geographical distribution have recently been reviewed . Studies carried out so far have shown that the M and S forms may have different habitat even in sympatric areas [5, 10].
The aim of this study was to compare the dynamics of transmission of Anopheles arabiensis and An. gambiae M and S molecular forms in a Senegalese village, where the two forms coexist.
The village of Dielmo (13°45N, 16°25W) is located in an area of Sudan-type savanna, 280 km Southeast of Dakar and about 15 km north of the Gambian border. About 300 inhabitants are living in the village. Rainfall occurs during a four-month period, from June to October. The average annual rainfall during our study period in 2004 was 642.4 mm. Dielmo is situated on the marshy bank of a small permanent stream, the Nema, where anopheline larval sites are present all year round. Only few cattle are living in this area .
Adult mosquitoes were collected monthly from July to December 2004. Two sampling methods were used: night landing catches (NLC) and pyrethrum spray catches (PSC). Hourly NLC were made on adult volunteers from 19.00 to 07.00 hours at the same two sites for three consecutive nights. Two collectors, one indoors and one outdoors, were positioned at each site. A total of twelve person-night of capture was done every month during the monthly collections. Six pyrethrum spray catches (PSC) was done during July to September early in the morning inside a total of twelve bedrooms, belonging to both types of houses, in different houses from those used for NLC.
Anophelines were identified to species level using morphological characteristics according to the identification keys of Gillies & De Meillon  and Gillies & Coetzee . Ovaries from a portion of females of An. gambiae s.l. captured on NLC were dissected to determine parity, by observing the degree of coiling of ovarian tracheoles . A total of 50–60 were randomly selected monthly. All mosquitoes belonging to the An. gambiae complex, dissected or not, were stored in individual tubes with silicagel and preserved at -20°C in the laboratory.
Bloodmeal sources of a sample of the females captured by PSC were tested by enzyme-linked immunosorbent assay (ELISA) to identify bovine, ovine, caprine (sheep and goat), equine (horse and donkey), or chicken hosts . Malaria infections were determined on the crushed head and thorax of all anopheline specimens by ELISA using monoclonal antibodies against Plasmodium falciparum circumsporozoite protein (CSP), as described by Wirtz et al : CSP rates and 95% confidence intervals were calculated. Species identification of the An. gambiae complex was done by PCR . Specimens identified as An. gambiae s.s. were tested for M and S molecular forms using the diagnostic PCR assay of Favia et al .
Parity, circumsporozoite protein and anthropophilic rates were estimated. Comparisons of these proportions between M and S molecular forms and the Hardy Weinberg equilibrium were performed by X2 test (P-values < 0.05 were considered to be statistically significant). The Human Biting Rate (HBR) was expressed as the average number of mosquito to bites per person per night during each month. The Entomological Inoculation Rate (EIR) was calculated by multiplying the (HBR) of each An. arabiensis and An. gambiae M and S molecular forms by the respective CSP rate, for each month.
A total of 1,109 An. gambiae s.l. were caught during 72 collection nights between July and December 2004 (889 by NLC and 220 by PSC). PCR for species identification was performed on 920 (83%) mosquitoes belonging to the gambiae complex, resulting in 177 An. arabiensis, 316 An. gambiae molecular form M, 405 An. gambiae molecular form S and 22 hybrids M/S (Table 1).
Parity and anthropophily rate
The mean parity rate was: 70.9% (n = 86), [CI: 60.1–80.2] for An. arabiensis, 68.7% (n = 64), [CI: 55.9–79.7] and 80.1% (n = 156), [CI: 73.0–86.1] for An. gambiae M and S forms, respectively. There was no significant difference between parity rates (X2 = 4.3; p = 0.11). The proportion of mosquitoes fed on human was 71.4% (n = 14) [CI: 41.9–91.6], 86.3% (n = 22) [CI: 65.1–97.1] and 91.6% (n = 24) [CI: 73.0–98.9], respectively, for An. arabiensis, An. gambiae M form and S form.
Nocturnal biting cycle
Details on hourly nocturnal biting cycle of the different species of An. gambiae s.l. are plotted in figure 1. HBR showed a consistent increased during night time reaching peaks between 11 p.m. and 5 a.m. depending on the species (Figure 1).
Human biting rate
An. gambiae S molecular form biting density was maximum in August (16.5 BPN), while the M molecular form was most abundant in September (17.6 BPN). An arabiensis, was less abundant and the maximum was recorded in August with 4.6 BPN (Figure 2).
Circumsporozoite protein and entomological inoculation rate
The results of CSP-ELISA assays are shown in Table 2. The mean CSP rate varied from 2.82% to 3.45%, according to the species. The differences between infection rates were not signicant (X2 = 0.16; p = 0.92). Among the hybrid species M/S specimens, no CSP antigen was found. During this 6 months study period, the mean EIR was 70 infective bites/person. An. arabiensis and An. gambiae forms M and S were responsible for 17%, 35% and 48% of P. falciparum transmission respectively.
Hardy Weinberg equilibrium
In total, 22 hybrids M/S were obtained. The frequencies observed were 0.4 and 0.6, respectively for the molecular forms M and S. At the Hardy Weinberg equilibrium, the number of M/S hybrids expected was 366 M/S heterozygous. Difference highly significant (p < 0.0001). The number of M/S hybrids obtained was signicantly inferior in the case of panmictic crossbreeding.
Anopheline densities varied during this study period. This strong variation is classical in savannah areas and is related to the presence of temporary ponds during the rainy season . For An. arabiensis this phenomenon was well known in Dielmo [1, 18]. The maximum biting rate by An. arabiensis, An. gambiae M and S forms occurred after 10 p.m., when most people are in bed, as generally observed throughout Africa . The late biting behaviour of the more dangerous mosquitoes is a significant finding with respect to the increase use of insecticide-treated mosquitoes nets (ITN) for the reduction of malaria. The overall parity rate is comparable to that reported by Lemasson et al  in Senegal. This indicates that the population is long lived and has the capacity to be an efficient vector of malaria during a relatively long period. The high parity rate can partially explain the CSP rates of An. arabiensis and An. gambiae in the present study which are greater than those previously recorded in other parts of Senegal [1, 18, 20, 21]. The high CSP rate in An. arabiensis is probably explained by the high antropophilic behaviour in this particular setting. Highton et al  and Joshi et al  in the same region Kisumu, found high and similar CSP rate for An. gambiae s.s. (5.3% and 7.5% respectively), while the CSP rate was very different for An. arabiensis (0.3% and 7.5% respectively). The difference observed in An.arabiensis was due to cattle abundance and movements.
A recent analysis of published and unpublished data on the molecular forms of An. gambiae has demonstrated that the M form shows a more latitudinal range in West Africa than the S form, being the only form recorded in the Sahelien region of northern Senegal [9, 24]. In Dielmo, the M form was mainly observed in rainy season from August to October, with a maximum in September. The S form was observed in the rainy season as well, with a maximum in August, but was also found during the dry season. An. arabiensis had an approximately constant density during the rainy season and the beginning of the dry season, without an obvious peak. Thus rainfall alone can not explain fluctuations observed in M and S forms densities.
Cytogenetic studies were conducted in 1991 in Dielmo and revealed an important chromosomic polymorphism. Thirty-two half-gravid females collected in January and February 1991 showed chromosomal inversions characteristic of bissau chromosomal form which belongs to the M molecular form . Twenty-four females collected in September 1991 mainly revealed savanna cytotype which belongs either to M form or S form . In this study, the replacement of S molecular form by M molecular form in the late rainy season might reflect the replacement of savanna chromosomal form by bissau chromosomal form. Climate (rainfall, temperature, relative humidity) as well as nature of anopheline larval development sites are dramatically change throughout the year during the rainy season. As chromosomal inversions are known to be involved in adaptation to climate and environment , the correspondence of molecular and chromosomal forms may explain the distribution pattern of molecular forms.
In West Africa, there is evidence of varying levels of hybridation between M and S forms, a mechanism by which adaptive genes may flow from one to the other, including those conferring insecticide resistance [25, 26]. No, or very few, M/S hybrids were observed throughout Cameroon [6, 27], Ghana , Mali and Burkina Faso . In this study, 22 M/S hybrids were observed (3.05%). Despite this number of hybrids, the An. gambiae population is far from Hardy Weinberg equilibrium suggesting restricted gene flow between M and S form. The occurrence of M/S hybrids in the field, however, has been reported several times. Della Torre et al  found three hybrids patterns out of 1,161 An. gambiae adult mosquitoes tested from throughout Africa (hybrid frequency = 0.26%). The authors, therefore, mentioned contamination as a possible cause of the hybrid patterns they reported. Further evidence for the viability of M/S hybrids in the wild was provided by Taylor et al , who reported occurrence of M/S hybrid larvae at a frequency of 0–1.29% in Banambani (Mali). Tripet et al  identified an inseminated female showing the M/S pattern, which demonstrated that M/S hybrids could be produced in the field, survive up to the adult stage and are reproductively active.
Anopheles arabiensis and the M and S molecular forms of An. gambiae coexist in Dielmo village. No difference was observed either in host preference or in Plasmodium falciparum infection rates between sympatric M and S populations, but they present different dynamics of transmission: the S form is the major vector during the first part of the rainy season and is replaced by M form later in the season. These variations are probably attributable to different breeding conditions.
Coluzzi M: Spatial distribution of chromosomal inversions and speciation in Anopheline mosquitoes. Prog Clin Biol Res. 1982, 96: 143-153.
Favia G, Lanfrancotti A, Spanos L, Siden-Kiamos I, Louis C: Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 19-23. 10.1046/j.1365-2583.2001.00236.x.
Coluzzi M, Petrarca V, Di Deco MA: Chromosomal inversion intergradation an incipient speciation in Anopheles gambiae. Boll Zool. 1985, 52: 45-63.
Toure YT, Petrarca V, Traore SF, Coulibaly A, Maiga HM, Sankare O, Sow M, Di Deco 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.
Diabate A, Dabire RK, Kim EH, Dalton R, Millogo N, Baldet T, Simard F, Gimnig JE, Hawley WA, Lehmann T: Larval development of the molecular forms of Anopheles gambiae (Diptera: Culicidae) in different habitats: a transplantation experiment. J Med Entomol. 2005, 42: 548-553. 10.1603/0022-2585(2005)042[0548:LDOTMF]2.0.CO;2.
Wondji C, Frederic S, Petrarca V, Etang J, Santolamazza F, Della Torre A, Fontenille D: Species and populations of the Anopheles gambiae complex in Cameroon with special emphasis on chromosomal and molecular forms of Anopheles gambiae s.s. J Med Entomol. 2005, 42: 998-1005. 10.1603/0022-2585(2005)042[0998:SAPOTA]2.0.CO;2.
Wondji C, Simard F, Fontenille D: Evidence for genetic differentiation between the molecular forms M and S within the Forest chromosomal form of Anopheles gambiae in an area of sympatry. Insect Mol Biol. 2002, 11: 11-19. 10.1046/j.0962-1075.2001.00306.x.
Della Torre A, Fanello C, Akogbeto M, Dossou-yovo J, Favia G, Petrarca V, Coluzzi M: Molecular evidence of incipient speciation within Anopheles gambiae s.s. in West Africa. Insect Mol Biol. 2001, 10: 9-18. 10.1046/j.1365-2583.2001.00235.x.
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.
Yawson AE, Weetman D, Wilson MD, Donnelly J: Ecological zones rather than molecular forms predict genetic differentiation in the malaria vector Anopheles gambiae s.s. in Ghana. Genetics. 2007, 175: 751-761. 10.1534/genetics.106.065888.
Trape JF, Rogier C, Konate L, Diagne N, Bouganali H, Canque B, Legros F, Badji A, Ndiaye G, Ndiaye P, Brahimi K, Faye O, Druilhe P, da Silva LP: The Dielmo project: a longitudinal study of natural malaria infection and the mechanisms of protective immunity in a community living in a holoendemic area of Senegal. Am J Trop Med Hyg. 1994, 51: 123-137.
Gillies MT, De Meillon D: The Anophelinae of Africa South of the Sahara. Publ South Afri Inst Med Res. 1968, 54: 343-
Gillies MT, Coetzee M: A supplement to the Anophelinae of Africa South of the Sahara, 2nd edn. Publ South Afri Inst Med Res. 1987, 55: 143-
Detinova TS, Gillies MT: Observations on the Determination of the Age Composition and Epidemiological Importance of Populations of Anopheles gambiae Giles and Anopheles funestus Giles in Tanganyika. Bull World Health Organ. 1964, 30: 23-28.
Fontenille D, Meunier JY, Nkondjio CA, Tchuinkam T: Use of circumsporozoite protein enzyme-linked immunosorbent assay compared with microscopic examination of salivary glands for calculation of malaria infectivity rates in mosquitoes (Diptera: Culicidae) from Cameroon. J Med Entomol. 2001, 38: 451-454.
Wirtz RA, Ballou WR, Schneider I, Chedid L, Gross MJ, Young JF, Hollingdale M, Diggs CL, Hockmeyer WT: Plasmodium falciparum: immunogenicity of circumsporozoite protein constructs produced in Escherichia coli. Exp Parasitol. 1987, 63: 166-172. 10.1016/0014-4894(87)90158-5.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Fontenille D, Lochouarn L, Diatta M, Sokhna C, Dia I, Diagne N, Lemasson JJ, Ba K, Tall A, Rogier C, Trape JF: Four years' entomological study of the transmission of seasonal malaria in Senegal and the bionomics of Anopheles gambiae and Anopheles arabiensis. Trans R Soc Trop Med Hyg. 1997, 91: 647-652. 10.1016/S0035-9203(97)90506-X.
Taye A, Hadis M, Adugna N, Tilahun D, Wirtz RA: Biting behavior and Plasmodium infection rates of Anopheles arabiensis from Sille, Ethiopia. Acta Trop. 2006, 97: 50-54. 10.1016/j.actatropica.2005.08.002.
Lemasson JJ, Fontenille D, Lochouarn L, Dia I, Simard F, Ba K, Diop A, Diatta M, Molez JF: Comparison of behavior and vector efficiency of Anopheles gambiae and Anopheles arabiensis (Diptera:Culicidae) in Barkedji, a Sahelian area of Senegal. J Med Entomol. 1997, 34: 396-403.
Dia I, Diallo D, Duchemin JB, Ba Y, Konate L, Costantini C, Diallo M: Comparisons of human-landing catches and odor-baited entry traps for sampling malaria vectors in Senegal. J Med Entomol. 2005, 42: 104-109. 10.1603/0022-2585(2005)042[0104:COHCAO]2.0.CO;2.
Highton RB, Bryan JH, Boreham PFL, Chandler JA: Studies on the sibling species Anopheles gambiae Giles and Anopheles arabiensis Patton (Diptera: Culicidae) in the Kisimu area, Kenya. Bull Ent Res. 1979, 69: 43-53.
Joshi GP, Service MW, Pradhan GD: A survey of species A and B of the Anopheles gambiae Giles complex in the Kisumu area of Kenya prior to insecticidal spraying with OMS-43 (Fenitrothion). Ann Trop Med Parasitol. 1975, 69: 91-104.
Carrara GC, Calzetta M, Santolamazza F, Di Deco MA, della Torre A, Petrarca V: Preliminary data on the distribution of species and forms of the Anopheles gambiae complex (Diptera: Culicidae) at sites of Angola. Parassitologia. 2004, 46: 85-
Fanello C, Santolamazza F, della Torre A: Simultaneous identification of species and molecular forms of the Anopheles gambiae complex by PCR-RFLP. Med Vet Entomol. 2002, 16: 461-464. 10.1046/j.1365-2915.2002.00393.x.
Yawson AE, McCall PJ, Wilson MD, Donnelly MJ: Species abundance and insecticide resistance of Anopheles gambiae in selected areas of Ghana and Burkina Faso. Med Vet Entomol. 2004, 18: 372-377. 10.1111/j.0269-283X.2004.00519.x.
Favia G, della Torre A, Bagayoko M, Lanfrancotti A, Sagnon N, Toure YT, Coluzzi M: Molecular identification of sympatric chromosomal forms of Anopheles gambiae and further evidence of their reproductive isolation. Insect Mol Biol. 1997, 6: 377-383. 10.1046/j.1365-2583.1997.00189.x.
Taylor C, Toure YT, Carnahan J, Norris DE, Dolo G, Traore SF, Edillo FE, Lanzaro GC: Gene flow among populations of the malaria vector, Anopheles gambiae, in Mali, West Africa. Genetics. 2001, 157: 743-750.
Tripet F, Toure YT, Taylor CE, Norris DE, Dolo G, Lanzaro GC: DNA analysis of transferred sperm reveals significant levels of gene flow between molecular forms of Anopheles gambiae. Mol Ecol. 2001, 10: 1725-1732. 10.1046/j.0962-1083.2001.01301.x.
We thank Charles Bouganali, Abdoulaye Gaye and Moise Gnafouna for their technical assistance, the villagers in Dielmo, Anna Cohuet for very helpful suggestions, and the DSF (Department support and formation of the south communities of IRD).
This work was supported by the programme Pal + programme of the French Ministry of Research and the Institute for Research and Development (IRD). This programme inclues three units of research: IRD-UR077 (Laboratoire de Paludologie Afrotropicale) – Dakar; UCAD (Université Cheikh Anta Diop) – Dakar and IRD-UR 016 (Laboratoire de Lutte contre les Insectes Nuisibles).
MON and DF have equally contributed to the design, acquisition, analysis, interpretation of data and manuscript drafting. LK contributed to conception of study and contributed markedly to the analysis of entomological data. CB for field activities and Molecular biology. JFT and CB participated in the conception and coordination of the study and helped to draft the manuscript. CS provided the scientific supervision in Dielmo.