Skip to main content

Preliminary survey on Anopheles species distribution in Botswana shows the presence of Anopheles gambiae and Anopheles funestus complexes



Botswana is one of the four front line malaria elimination countries in Southern Africa, with malaria control activities that include routine vector control. Past and recent studies have shown that Anopheles arabiensis is the only known vector of Plasmodium parasites in the country. This report presents a preliminary evaluation on Anopheles species composition in seven districts of Botswana with some inferences on their vectorial role.


Overall, 404 Anopheles mosquito females were collected, of which 196 were larvae collected from several breeding sites, and 208 were adults obtained from indoor pyrethrum spray catches (PSC). Anopheles arabiensis (58.9%) accounted for the highest relative frequency in 5 out of 7 districts sampled. The other species collected, among those identified, were barely represented: Anopheles longipalpis type C (16.3%), Anopheles parensis (8.9%), Anopheles quadriannulatus (5.4%), and Anopheles leesoni (0.2%). PCR test for human β-globin on mosquitoes collected by PSC showed that An. arabiensis and An. parensis had bitten human hosts. Moreover, An. arabiensis showed a non-negligible Plasmodium falciparum infection rate in two sites (3.0% and 2.5% in Chobe and Kweneng West districts, respectively).


This work provides first time evidence of Anopheles diversity in several areas of Botswana. Anopheles arabiensis is confirmed to be widespread in all the sampled districts and to be vector of P. falciparum. Moreover, the presence of Anopheles funestus group in Botswana has been assessed. Further research, entomological surveillance activities and possibly vector control programmes need to be better developed and implemented as well as targeting outdoors resting vectors.


Botswana is one of the four Southern African countries on the nearing malaria elimination (together with South Africa, Namibia and Swaziland). Therefore, knowledge of the transmission dynamics is critical in moving forward. Malaria mainly occurs in five Northern and Eastern districts (Okavango, Ngami, Chobe, Boteti and Tutume) with other districts being affected occasionally due to local outbreaks/epidemics, because the country’s ecosystem is receptive to malaria [1]. These regions experience active malaria circulation, especially during the peak of malaria vector breeding season that spans the summer months (November–April). Each of these districts have developed malaria control programmes, including routine vector control, which is primarily based on the application of indoor residual insecticide spraying (IRS) and use of long-lasting insecticide-treated nets (LLINs) [2]. In Botswana, the only known malaria vector is Anopheles arabiensis, a mosquito belonging to the Anopheles gambiae complex known to be a major vector of Plasmodium falciparum and Plasmodium vivax parasites. The presence of An. arabiensis as the main malaria vector in Botswana has been recently assessed on the Okavango region in North-West Botswana [35]. No further information is available from other areas of the country. Moreover, nothing is known about the presence of other species, such as those of the Anopheles funestus complex, which includes major malaria vector species of Southern Africa [68].

In general, several of these species are often found to occur in sympatry and their importance in malaria transmission varies depending on behaviour, seasonal prevalence and vectorial capacity. These peculiarities contribute to the varied malaria epidemiological patterns observed in a particular geographical region. Therefore, accurate identification of malaria vector system in a defined area is important in the understanding of transmission dynamics scenario. There is a clear need for Botswana to carry out vector surveillance studies in order to obtain baseline knowledge on seasonal prevalence and vectorial capacity in areas where malaria is of unstable endemicity. This can help in understanding the malaria vectorial system and, therefore, identify the risk areas and the local foci of potential transmission where malaria could be re-introduced. It is important to point out that Botswana has high vulnerability to malaria due to the fact that rainfall anomalies are widely considered to be a major driver of inter-annual variability of malaria incidence in particular in semi-arid areas of Southern Africa, following the El Nino Southern Oscillation (ENSO) pattern [2, 9, 10].

The aim of this work is to present the spectrum of Anopheles species composition of specimen collected in 2015 from several districts of Botswana, where malaria is of unstable endemicity, and from one Southern district at the time of a malaria outbreak in 2012. For some areas presented here, this is the first report of malaria vectors presence in recent years. Moreover, the information on the detection of human blood and the Plasmodium positivity rates in the specimens collected is provided.


Anopheles specimens were collected in 7 districts of Botswana (Kweneng West in March 2012; Okavango, Ngami, Chobe, Boteti, Tutume, Bobirwa in February–March 2015) (see Fig. 1). In 2011, seven sentinel sites were established to develop a surveillance system for malaria vector population and monitor vector susceptibility to insecticides. The sites represent three epidemiological zones: A-endemic malaria transmission (Okavango, Ngami, Chobe), B-moderate malaria transmission (Boteti, Tutume, Bobirwa), C-malaria free zones but prone to outbreaks (Kweneng West) [11]. Entomological surveillance was done in those 7 sites to provide insight information on vector density, monitor vector susceptibility to insecticides and assess quality of vector control interventions (IRS/LLINs). The locations were chosen after assessing the availability of potential breeding sites for mosquito sampling for each location (with radius of 2–10 km). Mosquitoes were collected both as larvae from several breeding sites and as resting adults by indoor Pyrethrum Spray Catches (PSC). The PSC were performed for 2–3 days in each site and 40 houses per site were sampled. Collected larvae were brought to the insectary and allowed to emerge as adult before morphological and genetic analysis was performed. Adult female specimens were first identified morphologically [12], before molecular identification in case of specimens belonging to sibling species complexes.

Fig. 1

Anopheles species composition and geographical distribution in 7 districts of Botswana

The protocol of Scott et al. [13] was applied for the species identification of An. gambiae sensu lato (s.l.) specimens and the protocol of Koekemoer et al. [6] for species identification of those belonging to An. funestus s.l. complex. Finally, to further identify the presence of Anopheles longipalpis type C in the An. funestus s.l. samples, the protocol of Choi et al. [14] was adopted.

As no information was taken on the gonotrophic stage of the specimens collected as adults, all specimens were tested for the presence of human blood through the detection of human β-globin DNA specific sequence, according to the protocol of Quinones et al. [15]. Moreover, the same specimens were also tested for the presence of P. falciparum using the molecular detection for pfmdr1 parasite gene through a nested-PCR approach [16]. For both tests we used appropriate positive controls (human DNA, and HB3 and DD2 P. falciparum parasite strains, respectively).


Four hundred and four (404) female Anopheles mosquitoes were collected in seven districts of the country, as shown in Table 1. Of these, 196 were larvae collected from several breeding sites, and 208 were adults obtained by indoor PSC. Forty one (41) out of 404 mosquitoes were not identified both morphologically and by molecular analysis because of the bad preservation status of the specimens (identification success rate 89.9%). In 3 out of 7 districts (Boteti, Tutume, Bobirwa) the proportion of unidentified specimens varied from 31 to 50% (Table 1). There was a large variability of species collected in the different sites (Fig. 1), although An. arabiensis (58.9%), was identified with higher relative frequency in 5 out of 7 districts. The other species collected, among those identified, were far less represented: Anopheles quadriannulatus (5.4%) which is a non-vector species belonging the An. gambiae complex and An. longipalpis type C (16.3%), Anopheles parensis (8.9%) and Anopheles leesoni (0.2%) all belonging to the An. funestus complex. The adult mosquitoes collected by PSC and tested for human blood showed that An. arabiensis and An. parensis had bitten human hosts (Table 2). In addition, An. arabiensis showed a noticeable P. falciparum infection rate in two sites (3.0% in Chobe and 2.5% in Kweneng West, Table 3), with one single mosquito found to be positive to parasite DNA in both sites.

Table 1 Scheme of Anopheles mosquitoes collected in the different districts sampled
Table 2 Human blood positivity in PSC Anopheles species
Table 3 Plasmodium falciparum positivity rate in PSC Anopheles species


In this report, we provide information about Anopheles diversity in several areas of Botswana, covering the Northern to the Southern parts of the country. The results, although with relatively small sample size, confirm the widespread presence of An. arabiensis in all the sampled districts. Diversity higher than previously reported [17] was observed, and in particular a noticeable presence of some species of the An. funestus group, which are unknown or potential vectors of Plasmodium parasites [18]. Furthermore, the presence of human DNA in indoor resting An. arabiensis and An. parensis in several districts was detected, and it was confirmed the possible vector role of An. arabiensis for P. falciparum in Chobe and Kweneng West districts. In fact, the P. falciparum infected An. arabiensis from Chobe district was not positive for human β-globin DNA, supporting the hypothesis that the sporogonic cycle of the malaria parasite was not at an early stage (blood meal digestion), but progressed towards oocyst formation or salivary glands invasion [19, 20]. Instead, for the An. arabiensis positive for P. falciparum and human β-globin DNA from Khudumelapye (Kweneng West) it is possible to infer either a recent P. falciparum-positive blood meal (with or without gametocytes) or a possible second blood meal of an already infected mosquito.

In general, the data does not provide any information about anthropophily of the species collected, but confirm that humans are among the hosts of indoor resting An. arabiensis and An. parensis in Botswana. While the role of An. arabiensis as malaria vector is generally known, this paper presents some evidence of anthropophagic behaviour of An. parensis in Botswana. Moreover, the latter data is supported by the record of P. falciparum infected specimens of An. parensis species in neighbouring South Africa [21].

Anopheles arabiensis populations exhibit both exophilic and endophilic behaviours and various degrees of anthropophily depending on the prevailing ecological conditions [22]. This has implications for the planning of vector control measures. Although in Botswana it has been broadly proclaimed that malaria transmission is mostly due to An. arabiensis, it is not known how the transmission process plays out (indoor or outdoor transmission). This is important in the selection of the most effective tools for combating the vector which indirectly affects the effectiveness of control measures for blocking transmission. Moreover, after more than 70 years of indoor residual spraying in the Botswana [23], it could be that indoor biting behaviour of the An. arabiensis population may have shifted from indoor to outdoor biting preference with possible impact on malaria transmission, as already observed elsewhere [24, 25]. Recently, a genetic component linked to the biting behaviour has been postulated [26], congruently with the hypothesis that a possible shift toward a higher degree of zoophily in An. arabiensis could be due to selection by indoor vector control activities. Conversely, insecticide resistance could have played an opposite role in selecting for insecticide resistant populations, thus negatively impacting on indoor residual spraying activity and ultimately in malaria control.

Unfortunately, Botswana lacks almost all data on local vector species and their susceptibility to insecticides, as well as on vector and human behaviours that may allow vectors to avoid contact with interventions and then maintain residual transmission. There is also a critical need for punctual monitoring of the coverage, usage, quality and durability of vector-control interventions such as IRS and LLINs. Evaluation of the impact of interventions on malaria outcomes should also be undertaken. Moreover, Botswana urges to invest better in public health entomology capacity, to support the control and elimination of malaria [27].

Despite the malaria elimination campaigns, Botswana still experiences a significant number of reported cases with deaths [28]. In the malaria elimination setting it is crucial to know where to target the interventions. In Botswana interventions are primarily focused on vector control but clearly additional measures, especially those targeting outdoor resting vectors, need to be evaluated. Entomological surveillance activities, including routine insecticide resistance monitoring, need to be scaled up in order to prevent new outbreaks driven by resistant and new vector populations.


This study provides a first report on the spread and abundance of An. gambiae and An. funestus complex in Botswana. Additionally, there is evidence of An. arabiensis as malaria vector in Botswana and a possible potential role of An. parensis. Future strategies on vector control must take into consideration more tools that target both indoor and outdoor transmission of Plasmodium species.


  1. 1. Accessed 14 Feb 2017.

  2. 2.

    Thomson MC, Mason SJ, Phindela T, Connor SJ. Use of rainfall and sea surface temperature monitoring for malaria early warning in Botswana. Am J Trop Med Hyg. 2005;73:214–21.

    PubMed  Google Scholar 

  3. 3.

    Chirebvu E, Chimbari MJ, Ngwenya BN. Assessment of risk factors associated with malaria transmission in tubu village, northern botswana. Malar Res Treat. 2014;2014:403069.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    Chirebvu E, Chimbari MJ. Characteristics of Anopheles arabiensis larval habitats in Tubu village, Botswana. J Vector Ecol. 2015;40:129–38.

    Article  PubMed  Google Scholar 

  5. 5.

    Chirebvu E, Chimbari MJ. Characterization of an indoor-resting population of Anopheles arabiensis (Diptera: Culicidae) and the implications on malaria transmission in Tubu village in Okavango subdistrict, Botswana. J Med Entomol. 2016;53:569–76.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66:804–11.

    CAS  PubMed  Google Scholar 

  7. 7.

    Coetzee M, Fontenille D. Advances in the study of Anopheles funestus, a major vector of malaria in Africa. Insect Biochem Mol Biol. 2004;34:599–605.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Choi KS, Koekemoer LL, Coetzee M. Population genetic structure of the major malaria vector Anopheles funestus s.s. and allied species in southern Africa. Parasit Vectors. 2012;5:283.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Mabaso ML, Kleinschmidt I, Sharp B, Smith T. El Niño Southern Oscillation (ENSO) and annual malaria incidence in Southern Africa. Trans R Soc Trop Med Hyg. 2007;101:326–30.

    Article  PubMed  Google Scholar 

  10. 10.

    MacLeod DA, Jones A, Di Giuseppe F, Caminade C, Morse AP. Demonstration of successful malaria forecasts for Botswana using an operational seasonal climate model. Environ Res Lett. 2015;10:044005.

    Article  Google Scholar 

  11. 11.

    Motshoge T, Ababio GK, Aleksenko L, Read J, Peloewetse E, Loeto M, et al. Molecular evidence of high rates of asymptomatic P. vivax infection and very low P. falciparum malaria in Botswana. BMC Infect Dis. 2016;16:520.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Gillies MT, Coetzee M. A supplement to the Anophelinae of Africa south of the Sahara (Afrotropical Region). Publ S Afr Inst Med Res. 1987;55:1–143.

    Google Scholar 

  13. 13.

    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–9.

    CAS  PubMed  Google Scholar 

  14. 14.

    Choi KS, Coetzee M, Koekemoer LL. Simultaneous identification of the Anopheles funestus group and Anopheles longipalpis type C by PCR-RFLP. Malar J. 2010;9:316.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Quiñones ML, Drakeley CJ, Müller O, Lines JD, Haywood M, Greenwood BM. Diversion of Anopheles gambiae from children to other hosts following exposure to permethrin-treated bednets. Med Vet Entomol. 2000;14:369–75.

    Article  PubMed  Google Scholar 

  16. 16.

    Djimdé A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourté Y, et al. A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001;344:257–63.

    Article  PubMed  Google Scholar 

  17. 17.

    Ministry of Health. Guidelines for malaria vector control in Botswana. Gaborone: Department of Public Health, National Malaria Control Programme, Ministry of Health; 2007.

    Google Scholar 

  18. 18.

    Dia I, Guelbeogo MW, Ayala D. Advances and perspectives in the study of the malaria mosquito Anopheles funestus. In: Manguin S, editor. Anopheles mosquitoes—new insights into malaria vectors. InTech Publ. 2013.

  19. 19.

    Mukabana WR, Takken W, Seda P, Killeen GF, Hawley WA, Knols BG. Extent of digestion affects the success of amplifying human DNA from blood meals of Anopheles gambiae (Diptera: Culicidae). Bull Entomol Res. 2002;92:233–9.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Oshaghi MA, Chavshin AR, Vatandoost H, Yaaghoobi F, Mohtarami F, Noorjah N. Effects of post-ingestion and physical conditions on PCR amplification of host blood meal DNA in mosquitoes. Exp Parasitol. 2006;112:232–6.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Mouatcho JC, Hargreaves K, Koekemoer LL, Brooke BD, Oliver SV, Hunt RH, et al. Indoor collections of the Anopheles funestus group (Diptera: Culicidae) in sprayed houses in northern KwaZulu-Natal, South Africa. Malar J. 2007;6:30.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Costantini C, Sagnon N, Della Torre A, Coluzzi M. Mosquito behavioural aspects of vector-human interactions in the Anopheles gambiae complex. Parassitologia. 1999;41:209–17.

    CAS  PubMed  Google Scholar 

  23. 23.

    Mabaso ML, Sharp B, Lengeler C. Historical review of malarial control in southern African with emphasis on the use of indoor residual house-spraying. Trop Med Int Health. 2004;9:846–56.

    Article  PubMed  Google Scholar 

  24. 24.

    Takken W. Do insecticide-treated bednets have an effect on malaria vectors? Trop Med Int Health. 2002;7:1022–130.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Russell TL, Govella NJ, Azizi S, Drakeley CJ, Kachur SP, Killeen GF. Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania. Malar J. 2011;10:80.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Main BJ, Lee Y, Ferguson HM, Kreppel KS, Kihonda A, Govella NJ, et al. The genetic basis of host preference and resting behavior in the major African malaria vector, Anopheles arabiensis. PLoS Genet. 2016;12:e1006303.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Mnzava AP, Macdonald MB, Knox TB, Temu EA, Shiff CJ. Malaria vector control at a crossroads: public health entomology and the drive to elimination. Trans R Soc Trop Med Hyg. 2014;108:550–4.

    Article  PubMed  Google Scholar 

  28. 28.

    WHO. World Malaria Report 2016. Geneva: World Health Organization; 2016.

Download references

Authors’ contributions

TL performed the experiments, data analysis, drafted and revised the manuscript. RP performed the experiments, data analysis, drafted and revised the manuscript. WK performed the sample collection, data analysis and revised the manuscript. MCW coordinated the lab work and revised the manuscript. MN coordinates the work and revised the manuscript. NDS supervised the field survey and revised the manuscript. QIK critically revised the manuscript. PM participated in the design, coordination, data analysis and manuscript revision. PGM set and design the study, supervised all the laboratory activities, analysed the data and draft and revised the manuscript. All authors read and approved the final manuscript.


This work was based at the University of Botswana under the supervision and financial support from the Ministry of Health and in collaboration with the University of Pennsylvania (USA) and “Sapienza” University of Rome (Italy). Positive controls for all the members of the An. gambiae and An. funestus groups were kindly provided from M. Pombi and L. Koekemoer, at “Sapienza” University of Rome (Italy) and at University of the Witwatersrand in Johannesburg (South Africa), respectively.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.


This publication was made possible through core services and support from Botswana Ministry of Health and the Penn Center for AIDS Research (CFAR), an NIH-funded program (P30 AI 045008).

Author information



Corresponding author

Correspondence to Giacomo Maria Paganotti.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tawe, L., Ramatlho, P., Waniwa, K. et al. Preliminary survey on Anopheles species distribution in Botswana shows the presence of Anopheles gambiae and Anopheles funestus complexes. Malar J 16, 106 (2017).

Download citation


  • Anopheles gambiae
  • Anopheles funestus
  • Botswana
  • Malaria
  • Plasmodium falciparum