key: cord-0022387-63ggfz9m
authors: Buxton, Mmabaledi; Nyamukondiwa, Casper; Wasserman, Ryan J.; Othenin-Girard, Victor; Pigeault, Romain; Christe, Philippe; Glaizot, Olivier
title: Surveillance Studies Reveal Diverse and Potentially Pathogenic-Incriminated Vector Mosquito Species across Major Botswana Touristic Hotspots
date: 2021-10-06
journal: Insects
DOI: 10.3390/insects12100913
sha: 44bf5bf8a399e2b1a5c34feab580c2b13be3f08e
doc_id: 22387
cord_uid: 63ggfz9m

SIMPLE SUMMARY: Mosquitoes vector pathogens that cause burdening diseases in humans, livestock and wildlife worldwide. Spatially and temporally, mosquito diversity varies considerably in response to bio-physical environments. As such, there is a need for mosquito diversity and distribution studies, as well as monitoring programmes, to inform on the risk of associated diseases. This survey assessed mosquito species in three major touristic areas of Botswana that are likely to harbour pathogens across prevailing hosts. The results revealed that all regions surveyed had important mosquito groups (Anopheles, Aedes and Culex) that are threats to public, wildlife and livestock health globally, including the arid Central Kalahari Game Reserve. The findings represent useful species inventories for future surveys and monitoring programmes. ABSTRACT: Vector mosquitoes contribute significantly to the global burden of diseases in humans, livestock and wildlife. As such, the spatial distribution and abundance of mosquito species and their surveillance cannot be ignored. Here, we surveyed mosquito species across major tourism hotspots in semi-arid Botswana, including, for the first time, the Central Kalahari Game Reserve. Our results reported several mosquito species across seven genera, belonging to Aedes, Anopheles, Culex, Mansonia, Mimomyia, Coquillettidia and Uranotaenia. These results document a significant species inventory that may inform early warning vector-borne disease control systems and likely help manage the risk of emerging and re-emerging mosquito-borne infections.

Arthropods are economically important organisms given their role in the global transmission of disease to humans, livestock and wildlife [1, 2] . Amongst arthropods, mosquitoes are by far the most important vector species contributing to global human and animal health burdens [3, 4] . Indeed, Aedes, Anopheles, Culex and Mansonia have been implicated as the main vector mosquito groups of medical and veterinary importance [5] . Vector mosquitoes transmit pathogens (e.g., viruses, protozoans and helminths) that cause diseases such as Zika, dengue, West Nile fever, yellow fever, Rift Valley fever, types of encephalitis (such as Japanese and St. Louis), chikungunya, filariasis and human and avian malaria [6] [7] [8] [9] . Owing to the socio-economic burdens brought about by vector mosquitoes, adequate and updated knowledge on their species diversity, spatio-temporal distribution and abundance are thus key in monitoring and evaluating prevalence and the risk of their associated diseases.

The proliferation of mosquito vectors and associated pathogens poses a challenge in health and epidemiological systems, mainly in Africa [10] . Given the increase in invasive vector species with climate change [11] , increased human populations, trade and globalisation [12, 13] , invasive vector mosquito species could exacerbate this problem. Indeed, mosquitoes contribute immensely to the global socio-economic burden of a global crisis in rising invasive insect species, and their economic costs are increasing significantly [14, 15] . Mosquito species have the potential to successfully establish in diverse landscapes influenced by foraging (e.g., host availability, nectar dietary resources [16] or climatic niches [17] ). Thus, in areas where there is a broad feeding preference (e.g., humans, livestock and wildlife), and favourable abiotic factors, mosquito species are likely to proliferate [18] . As such, novel environmental monuments such as game reserves and national parks with potentially diverse animal hosts, human influx, vegetation and perennial water supplies become a crucial ecological platform for vector-pathogen and host interaction fostering emerging and re-emerging mosquito-borne disease transmission dynamics [19] . Therefore, there is need for robust and continuous monitoring across major mosquito-vector-species risk pathways to update the inventory status [18] .

In Africa, the most prominent mosquito-borne disease is malaria [20] (but also see Willcox et al. [21] ; Mwanyika et al. [22] ), accounting for~400,000 deaths annually [23] . As such, most of the published works have mainly focused on Anopheles vector species and their pathogen incrimination surveillance [24] (but see Braack et al. [9] ). Similarly, in Botswana much work is focused on human malaria vectors, neglecting other species [25] . In addition, mosquito research is also spatially skewed, mostly based in the malaria endemic zone (northern part) of the country and around permanent human habitation [26] . As such, there is a need for more holistic vector and pathogen surveillance that collectively covers all potential vector/parasite species and across different spatial scales.

Here, we assessed the diversity of mosquito vector species in major tourism hotspot regions of Botswana, including the Okavango Delta, the Chobe enclave and, for the first time, the Central Kalahari Game Reserve (CKGR). Currently, tourism is the second largest foreign exchange earner in Botswana, after diamonds [27] . Ecotourism comprises about 5% of the gross domestic product of the country's economy [28] , with the Okavango Delta, the Chobe enclave and the CKGR being the main centres of tourist activities [29] . Whilst previous works on mosquito diversity and vector competence have been conducted in parts of the Okavango and Chobe areas and other parts of the country [26, 30, 31] (see Table 1 ), there is still much exploration required given the potential anthropogenic and climate change likely to impact vector distribution and abundance [32] . This study thus provides an important species diversity database that may help strengthen vector surveillance, disease risk analysis and early warning systems in high-risk regions that are frequented by travelling humans. Table 2 for details of stations. Thumbnail picture shows the location of Botswana in Africa.

In February 2020, 2 stations west of Okavango (North-West District) were prospected during 3 and 2 nights, respectively. Sampling was done again at the VTR Research Centre for 5 nights (Figure 1 ; Table 2 ). All samples were collected during the summer season of the southern hemisphere, when vector mosquitoes are most abundant [38] . The CKGR is a hot desert region with a distinct dry (April to October) and wet (November to March) season. Semi-permanent but scarce artificial water bodies are found across the region. During the wet season, the landscape becomes green and dry with no rains [39] . The Okavango Delta is largely characterised by permanent streams and swampy basins to seasonal floodplains. Although the wet and dry seasons follow the same pattern as the CKGR, the Okavango Delta is rich in fauna and flora; however, both regions remain unique in attracting tourists from diverse local and international destinations [40] . Whilst the Chobe enclave consists of permanent river systems in the northern part of the district, which sustain wildlife, the southern region is relatively dry [41] . Sampling was carried out on the dry area (southern part) away from the Chobe River, i.e., around the VTR Centre (9 km south of Parakarungu village). Whilst the Chobe enclave and the Okavango Delta are largely characterised by a livestock and wildlife interface, the CKGR site is predominantly characterised by wildlife [42] .

Six BG-Sentinel carbon-dioxide-baited traps (Biogents AG, Regensberg, Germany) were set in the afternoon (1600 h) and the mosquitoes were collected in the morning (0600 h). The specimens were preserved individually in 2 mL Eppendorf tubes with 80% ethanol. However, because of overwhelming adult mosquito catches in the Okavango (up to 1500/night), the mosquitoes were first separated into their genera in the field using gross morphology [43] and then stored by identified taxon in 25 or 50 mL vials.

All specimens were identified morphologically following the protocols of Jupp [43] under a stereo binocular microscope (Leica M205, Leica microsystems, Wetzlar, Germany). Further, 1 to 40 samples per taxon, depending on the number of individuals collected, were subjected to molecular identification following Buxton et al. [36] . Here, the cytochrome oxidase subunit I (COI) gene was amplified using LCO1490 (5 -GGTCAACAAATCATAAAGATATTGG-3 ) and HCO2198 (5 -TAAACTTCAGGGTGACC AAAAAATCA-3 ) primers [44] . Sequences were examined with the BLAST algorithm in GenBank (accessed in June 2021). Molecular identification was used whenever blast identity was above 98%. When sequences were not available in GenBank (or blast identity was ≤ 98%), we used morphological identities. Unidentified specimens (morphological and molecular) were given to the genera level by groups of morphologically identical specimens (e.g., sp. 1, sp. 2).

Mosquito species richness (i.e., total species represented) per site was determined and abundances were analysed as proportions using percentages across sampled areas to compare dominance. The Shannon-Wiener diversity index was determined (equation 1), which takes into account the species abundance and evenness (i.e., the proportion of individuals amongst the different species sampled) following protocols of Khoobdel et al. [45] . Thus, the Shannon diversity index (H ) was calculated as the summation of the product of species abundance probability (Pi) and its natural logarithm (lnPi).

A total of 5486 mosquitoes from 32 taxa of seven genera (Aedes, Anopheles, Culex, Mansonia, Mimomyia, Coquillettidia and Uranotaenia) were collected (Table 3 ). In the Chobe enclave (CHO), 173 mosquitoes of 9 species in 3 genera were recorded (Aedes: Ae. aegypti; Anopheles: An. gambiae s.s, An. squamosus gp. 1; and Culex: Cx. naevei/simpsoni, Cx. perexiguus, Cx. pipiens, Cx. simpliciforceps, Cx. wigglesworthi and Cx. univittatus). Two unidentified Aedes specimens (Ae. sp. 1 1 and Ae. sp. 3 1 ) and eight unidentified Culex (Cx. sp. 1 1 ) were collected in Chobe as well (Table 3) . Culex pipiens was the most dominant species (75.72%), followed by An. gambiae s.s (11.56%). Other mosquito taxa reported were few, ranging from one (0.58%) to eight (4.62%) specimens per taxa (Figure 2a) . Table 3 . Mosquito identification information of species collected from major tourism destinations (Central Kalahari Game Reserve (CKGR), Okavango (OKA) and Chobe enclave (CHO)) in Botswana. Specimens were identified morphologically (morph.), molecularly (mol.) or both (morph/mol.) as described in the Materials and Methods section. The GenBank reference is given with pairwise identities (%). The GenBank corresponding reference, if available, is given whenever morphological identity is undefined at the species level. 1 Specimens unidentified morphologically with no or weak match in GenBank. 2 Morphologically identified but species absent from GenBank. 3 Specimens from CKGR and CHO with bad sequences, morphological identification.

The Okavango Delta (OKA) was by far the richest region with 5266 mosquitoes collected (96% of all catches), representing 17 mosquito species and a further 5 unidentified taxa (Table 3, Figure 2b ). The area was largely dominated by two Mansonia species (Ma. uniformis and Ma. africana), which account for 60.9% and 24.5% of the mosquitoes caught in the region, respectively. Other frequent species were Cx. naevei/simpsoni (4.35%) and Cx. poicilipes (3.11%). Seven other species of Culex were caught in the station with lower abundances, ranging from 1 to 111 individuals ( Table 3 ). The proportion of Cx. pipiens (less than 1%) was relatively low compared to their frequency in other regions. Other species were Ae. mcintoshi, seven species of Anopheles (An. gambiae being the most prevalent with 1.41% of the catches), three species of Coquillettidia (Cq. fuscopennata, Cq. metallica, Cq. microannulata), two species each of Mimomyia (Mi. mimomyiaformis, Mi. splendens) and Uranotaenia (Uranotaenia sp. 1, Uranotaenia sp. 2) ( Table 3) .

The Central Kalahari Game Reserve (CKGR) reported one species of Aedes (Ae. aegypti, 2.13%) and three species of Culex (Cx. univittatus, Cx. perexiguus and Cx. pipiens) representing 2.13%, 6.38% and 55.32% of the specimens collected, respectively (Table 3, Figure 2c) . It was the least diverse region sampled, and the majority and most abundant species (Cx. pipiens) was collected at the gate (CKGR_1), where the rangers' camp was located. Other individuals were of three unidentified species of Aedes and one species of Culex, all being relatively rare (three to five specimens).

Okavango reported the highest species richness (22) , the Kalahari had the lowest (8) and Chobe was intermediate (12) . The species diversity index was highest in the Kalahari (1.491), whereas Chobe had the lowest diversity (0.973) with Okavango being intermediate (1.207) . Overall, all the sampled areas had globally important mosquito genera of medical and veterinary importance (Aedes, Anopheles and Culex) with Cx. pipiens species reported in all study sites (Table 4 ). 

The current study documents 32 mosquito taxa (22 identified at the species level and 10 at the genera level) belonging to 7 genera in Botswana's tourism hotspots. To the best of our knowledge, this is the first study to date documenting mosquito species from the Kalahari Desert and one of the few in other major touristic hotspots in Botswana (but see Pachka et al. [30] ; Cornel et al. [31] ). Congruent to Cornel et al. [31] , we found mosquito species of important medical and veterinary concern (see Table 4 ). In particular, the results showed that Cx. pipiens dominated other species in the Chobe enclave and CKGR. Most of the Cx. pipiens in CKGR were collected in a station next to the park entrance (CKGR_1, Figure 1 ) where people from the park are living, in accordance with its anthropogenic habits such as the colonisation of sewage systems. Mansonia uniformis and, to a lesser extent, Ma. africana by far dominated the Okavango area. Mansonia uniformis is known to feed mostly on mammals (e.g., cattle, humans) but readily bites other vertebrate hosts, including avian species [64, 65] . The alarming high numbers of the duo species should serve as a potential early warning sign of the risk of emerging-re-emerging diseases and pathogen circulation across an interface of humans, livestock and wildlife in the region. This study is also the first to report on the invasive species Ae. aegypti in the Chobe enclave and in the CKGR. It has been observed in Botswana for the first time in the south-east of the country, close to the South African border [36] . Apart from transmitting several pathogens of medical importance (e.g., Zika, dengue, yellow fever and chikungunya viruses), Ae. aegypti has also been reported to transmit lumpy skin virus in cattle [47] and potentially to wildlife. The unidentified species of Aedes (and other unidentified taxa) have implications for the risk of associated diseases in touristic areas; moreover, their pathogenicity status has not yet been assessed. Furthermore, we also report mosquito species belonging to Coquillettidia, Mansonia, Mimomyia and Uranotaenia only in the swampy Okavango, suggesting spatial environmental context effects in determining mosquito species diversity, given the highest richness reported in the area [66] .

Differences between the three study sites are likely driven by differences in environmental conditions. The CKGR is more arid, with no permanent water sources supporting mosquito breeding, than to the Chobe enclave and Okavango. Thus, mosquito breeding in the CKGR during dry seasons may be restricted to artificial water holes designed to water wildlife. This may partly account for the fewer mosquito species recorded in the CKGR. Conversely, Okavango consists of semi-permanent basins that are likely to sustain mosquito-breeding success, leading to elevated species richness in the area. For instance, Coquillettidia, Mimomyia and Mansonia species were reported only in the Okavango likely supported by the swampy wetland systems with large clear water bodies, whereby their larvae attach to the submerged aquatic plants for sustained feeding and oxygen intake [67] . The Chobe enclave has permanent water bodies likely to support mosquito species, congruent with Buxton et al. [33] . However, sampling was performed in a relatively arid part of the Chobe enclave, where breeding opportunities are mostly related to the presence of human infrastructure, as indicated by the majority of anthropogenic species found in this region (Ae. aegypti, An. gambiae s.s and Cx. pipiens). Large water bodies, through their semi-permanent swamps, and water that collects in animal hoof prints are a major source of facilitating mosquito-breeding refuge [33, 68] . Although habitat characterisation was not included in this study, resource portioning and species distribution remains critical for mosquito vector success. This has an implication on the risk of pathogens they may potentially transmit and the distribution of 'previously unidentified' mosquitoes to species level across touristic regions. It is well known that increasing human populations, international travel and globalisation have implications on importing and/or exporting vectors, pathogens and parasites from native ranges or areas of most prevalence [69] . Despite this, however, mosquito diversity from tourism hotspots in Botswana remains relatively unexplored [30, 31, 70] .

The presence of numerous vector species in tourist areas in Botswana, including the invasive Ae. aegypti, a major contributor of dengue infections in the tropics and subtropics globally [71] , highlights the need for rigorous mosquito vector surveillance in these areas. Therefore, the use of different trapping techniques in time and space may be essential in exploring a full spectrum of mosquito species diversity across varying seasonal and environmental niches [72] . More robust and spatially even sampling protocols with standardised effort may also be needed to get more accurate diversity indices across the studied sites. Furthermore, pathogen screening in vectors and prevailing vertebrate hosts (humans, livestock and wildlife) is pivotal in assessing transmission dynamics and the potential risk of emerging and re-emerging diseases in the country. This study also reported Cx. pipiens mosquitoes, which exist as a species complex. Thus, future studies may explore genetic variation among sibling species and consider surveillance of the diversity of other arthropod vectors capable of incriminating humans, livestock and wildlife with their associated pathogens. This study nevertheless provides significant baseline vector diversity data that may play a significant role in identifying infection risk pathways in pest and disease risk analysis and under emerging and re-emerging vector-borne infections. 

Data Availability Statement: The datasets during and/or analysed during the current study are available from the corresponding author on reasonable request.

Arthropods as vectors of transmissible diseases in Spain

Arthropods: Why It Is So Crucial to Know Their Biodiversity? In Measuring Arthropod Biodiversity

Mosquitoes (Diptera: Culicidae) and their relevance as disease vectors in the city of Vienna

Parasites and wildlife in a changing world: The vector-host-pathogen interaction as a learning case

A detailed review of the mosquitoes (Diptera: Culicidae) of Iran and their medical and veterinary importance

A global map of dominant malaria vectors

Identifying avian malaria vectors: Sampling methods influence outcomes. Parasites Vectors

Distribution and abundance of key vectors of Rift Valley fever and other arboviruses in two ecologically distinct counties in Kenya

Mosquito-borne arboviruses of African origin: Review of key viruses and vectors. Parasites Vectors

Changes in malaria epidemiology in Africa and new challenges for elimination

No saturation in the accumulation of alien species worldwide

Global consumption and international trade in deforestation-associated commodities could influence malaria risk

Human practices promote presence and abundance of disease-transmitting mosquito species

High and rising economic costs of biological invasions worldwide

The magnitude, diversity, and distribution of the economic costs of invasive terrestrial invertebrates worldwide

Implications of diet on mosquito life history traits and pathogen transmission

Niche conservatism of Aedes albopictus and Aedes aegypti-two mosquito species with different invasion histories

Surveillance of Culicine Mosquitoes in Six Villages of Taita-Taveta County, Kenya, With Host Determinations From Blood-Fed Females

Species Composition and Risk of Transmission of Aedes-Borne Arboviruses around the Mole Game Reserve in Northern Ghana

Current Situation of Malaria in Africa

Seroepidemiology of dengue, Zika, and yellow fever viruses among children in the Democratic Republic of the Congo

Dengue Virus Infection and Associated Risk Factors in Africa: A Systematic Review and Meta-Analysis. Viruses

African malaria deaths set to dwarf covid-19 fatalities as pandemic hits control efforts, WHO warns

Implicating cryptic and novel anophelines as malaria vectors in Africa

Distribution of Anopheles mosquito species, their vectorial role and profiling of knock-down resistance mutations in Botswana

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

Developing Tourism in Botswana: Progress and Challenges; Botswana Documents; Botswana institute for Development Policy Analysis: Gaborone

Tourism in Botswana in the last 50 years: A review

The dynamics of tourist visitation to national parks and game reserves in Botswana

Rift Valley fever vector diversity and impact of meteorological and environmental factors on Culex pipiens dynamics in the Okavango Delta

Mosquito community composition in South Africa and some neighboring countries. Parasites Vectors

Projected shifts in the distribution of malaria vectors due to climate change

Disease Vector Relative Spatio-Temporal Abundances to Water Bodies and Thermal Fitness Across Malaria Endemic Semi-Arid Areas

Past and current biological factors affecting malaria in the low transmission setting of Botswana: A review

Bionomics and cytogenetics of Anopheles seretsei Abdulla-Khan, Coetzee, and Hunt, a new species from northern Botswana

First record of Aedes (Stegomyia) aegypti (Linnaeus, 1762) (Diptera: Culicidae) in Botswana. Bioinvasions Rec

Cattle-induced eutrophication favours disease-vector mosquitoes

Mosquito community composition and abundance at contrasting sites in northern South Africa

Mapping vegetation morphology types in southern Africa savanna using MODIS time-series metrics: A case study of central Kalahari

Vander Post, C. Species diversity of the Okavango delta

Resource partitioning within a browsing guild in a key habitat, the Chobe Riverfront

Preliminary Assessment of Bovine Tuberculosis at the Livestock/Wildlife Interface in two Protected Areas of Northern Botswana

Mosquitoes of Southern Africa: Culicinae and toxorhynchitinae

DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates

Species diversity of mosquitoes of the Genus Culex (Diptera, Culicidae) in the coastal areas of the Persian Gulf

Aedes mosquitoes and Aedes-borne arboviruses in Africa: Current and future threats

Transmission of lumpy skin disease virus: A short review

Predicting Abundances of Aedes mcintoshi, a primary Rift Valley fever virus mosquito vector

Differential contribution of Anopheles coustani and Anopheles arabiensis to the transmission of Plasmodium falciparum and Plasmodium vivax in two neighbouring villages of Madagascar. Parasites Vectors 2020

A review of mosquitoes associated with Rift Valley fever virus in Madagascar

The relationship between size and longevity of the malaria vector Anopheles gambiae (ss) depends on the larval environment. Parasites Vectors

Rickettsia species in African Anopheles mosquitoes

Unexpected anthropophily in the potential secondary malaria vectors Anopheles coustani sl and Anopheles squamosus in Macha, Zambia. Vector Borne Zoonotic Dis

Mosquitoes of western Uganda

Nonspecific patterns of vector, host and avian malaria parasite associations in a central African rainforest

Coquillettidia (Culicidae, Diptera) mosquitoes are natural vectors of avian malaria in Africa

Feeding behaviour of potential vectors of West Nile virus in Senegal. Parasites Vectors

The role of different Culex mosquito species in the transmission of West Nile virus and avian malaria parasites in Mediterranean areas

Lumpy Skin Disease Virus

Elsevier Science publishers: Amsterdam, The Netherlands

Assessment of the risk of mosquitoes nuisance and appearance of vector-borne diseases in rubber production areas in Southern Côte d'Ivoire

Mansonia africana and Mansonia uniformis are vectors in the transmission of Wuchereria bancrofti lymphatic filariasis in Ghana

Identification of Wolbachia strains in mosquito disease vectors

Isolations of West Nile and Bagaza viruses from mosquitoes (Diptera: Culicidae) in central Senegal (Ferlo)

Seasonal abundance and biting behaviour of Mansonia annulifera. M. uniformis and M. indiana and their relative role in the transmission of malayan filariasis in shertallai (Kerala state)

Biodiversity of mosquitoes and Mansonia uniformis as a potential vector of Wuchereria bancrofti in Hulu Sungai Utara District

Landscape structure affects distribution of potential disease vectors (Diptera: Culicidae). Parasites Vectors

Colonization of Coquillettidia linealis (Skuse) with reference to other Coquillettidia and Mansonia species

Spatial distribution and habitat characterization of mosquito species during the dry season along the Mara River and its tributaries, in Kenya and Tanzania

Mapping population and pathogen movements

Characteristics of Anopheles arabiensis larval habitats in Tubu village

Adult Aedes abundance and risk of dengue transmission

Traps and trapping techniques for adult mosquito control

We thank the Botswana International University of Science and Technology (BIUST) for providing funding and infrastructure support to M.B., C.N. and R.J.W., the University of Lausanne and the Swiss National Science Foundation, grant number 31003A_179378, for funding R.P. and P.C. and the Museum of Zoology in Lausanne for funding O.G. The VanThuyne-Ridge (VTR) Research Centre is acknowledged for providing research facilities and infrastructure. In addition, we thank the Ministry of Environment, Natural Resources Conservation and Tourism (Botswana) for the issuing of a research permit (ENT 8/36/4 XXX II [10] ). Further, we sincerely acknowledge the contributions made by Chiswaniso Picture Kakambi, Ali Muchuka Mainga, Julie Glaizot, Kealeboga Matshameko and Mogomotsi Mokgadi in data collection and field work proceedings. We also thank Molly Baur and Jérôme Wassef for the molecular identification of mosquito species.

All authors declare no conflict of interests.