A global map of dominant malaria vectors
Parasites & Vectorsvolume 5, Article number: 69 (2012)
Global maps, in particular those based on vector distributions, have long been used to help visualise the global extent of malaria. Few, however, have been created with the support of a comprehensive and extensive evidence-based approach.
Here we describe the generation of a global map of the dominant vector species (DVS) of malaria that makes use of predicted distribution maps for individual species or species complexes.
Our global map highlights the spatial variability in the complexity of the vector situation. In Africa, An. gambiae, An. arabiensis and An. funestus are co-dominant across much of the continent, whereas in the Asian-Pacific region there is a highly complex situation with multi-species coexistence and variable species dominance.
The competence of the mapping methodology to accurately portray DVS distributions is discussed. The comprehensive and contemporary database of species-specific spatial occurrence (currently available on request) will be made directly available via the Malaria Atlas Project (MAP) website from early 2012.
Global malaria vector maps, by necessity, must simplify a complex diversity of numerous interacting and sympatric anopheline species. Such simplification refines the information down to a minimum, indicating only the primary vector(s) at each location and provides users, such as public health officials, modellers and opinion formers, with a global and regional picture that is easy to digest and utilise for scientific, operational and advocacy purposes.
Global maps have long been used to aid in visualising the malaria problem. These include the vector species map of May  and the 12 zones of malaria epidemiology described by Macdonald , determined using broad climatic ranges and physical land features, as well as consideration of the known distribution of the major anopheline vectors at the time. More recently Mouchet et al.  updated Macdonald's map, reassigning the 12 zones into more conventional biogeographical regions. This history of malaria vector (or vector-associated) visualisation indicates a past appetite for such maps, continuing more recently with Kiszewski et al.  publishing a global distribution map for the major malaria vectors in 2004. Their map was created to aid the authors in the development of a malaria transmission 'stability' map, but has since been adopted widely within the malaria research community and reproduced in many publications (their paper is listed as being cited 81 times in Web of Science and 37 times in PubMED). There is, therefore, a substantial and continuing demand for global maps of the major vectors of malaria.
Human malarial protozoa are transmitted by mosquitoes of the genus Anopheles, which includes 465 formally recognised species and more than 50 unnamed members of species complexes . Approximately 70 of these species have the capacity to transmit human malaria parasites  and 41 are considered here to be dominant vector species/species complexes (DVS), capable of transmitting malaria at a level of major concern to public health .
A comprehensive database of contemporary occurrence data for the 41 DVS was compiled over two years, beginning in January 2008 [7–10]. Using these and other data (see methods), distribution maps were produced for each species or species complex, which have been made available for download via the Malaria Atlas Project (MAP) website . This paper describes the production of a global map of dominant malaria vectors using these individual species maps. No other published vector map has had the benefit of the extensive and comprehensive evidence base that underlies the work presented here.
A list of 41 DVS were identified by consulting a number of authoritative reviews [3, 4, 12–14] ( now updated and translated ), including that of Kiszewski et al. , and included all those species or species complexes that were identified as 'principal', 'main', or 'dominant' vectors of malaria. Occurrence data for these 41 DVS were assembled into the MAP  vector database, incorporating published records of contemporary (post 1985) species-specific location information (Table 1). The database includes over 4800 sources relating to 15837 occurrence data points and also holds the other elements required to produce predictive distribution maps for all species/species complexes including a suite of open access environmental or climatic variables and expert opinion (EO) range maps.
Using the boosted regression tree (BRT) modelling methodology , predicted distribution maps were produced for each DVS. Nine species/species complex maps were produced for the Americas, seven for Africa, six for Europe and the Middle East and 19 for the Asian-Pacific region.
Building the global map
In ArcMap , each predictive species map was buffered at a 50 km limit beyond the EO boundary and only those pixels with a presence probability greater than 0.5 were included.
On a country-by-country basis, and by region, a list of all the DVS known and predicted to occur in each malaria endemic country was created. The lists were circulated to the project Technical Advisory Group (TAG) who identified the three most important DVS per country (where there were three or more species) and ranked these species by their relative importance. Due to the known complexity of DVS found in the Asian-Pacific region, all DVS found in each country within the region were ranked by importance. Additionally, for some countries (e.g. Indonesia), where the importance of a vector species can vary greatly across a country's geographical extent, more detailed species-specific spatial information was gathered. These rankings were used to guide the creation of the multi-species maps where every attempt was made to ensure that the top ranked species in specific countries and regions were uppermost in the map species layers (i.e. were displayed preferentially over the less highly ranked species). In those areas where the dominance of one species is not clear (e.g. An. arabiensis, An. funestus, An. gambiae across sub-Saharan Africa), the predictive maps for each were merged to indicate the presence of a combination of equally dominant species.
In the Asian-Pacific region, An. subpictus s.l., An. aconitus and An. annularis, which were all included in the original list of 41 DVS (Table 1), are not shown on the multi-species maps. Anopheles aconitus and An. annularis both tend to play only a focal role in malaria transmission within their respective ranges and are often considered secondary or incidental [18–20], however, in 'ideal conditions' An. aconitus can be a major DVS and similarly, An. annularis is important only in selected areas in India, Sri Lanka and Nepal. Both species are essentially zoophilic [19, 21, 22], as is An. subpictus s.l. and with this latter DVS, there is question about the reported identification of specimens based only on morphological characteristics. Thus in the assessment of dominance in this region, these three species were ranked lowest and therefore when represented on the multi-species map, were overlaid completely by other more important species and were removed.
Due to the fine detail available for each of the predictive maps, in addition to the global distribution, maps are also presented for those regions that contain countries with a high burden of malaria, i.e. the Americas, Africa and the Asian-Pacific region. Moreover, additional maps are also provided which highlight areas where there is a particularly high diversity of vector species (e.g. Central America, South-East Asia and Pacific).
Finally, two maps are presented for the African region. The first illustrates the distribution of those species considered to be of major importance, even within the confines of the DVS ranking (i.e. An. gambiae, An. arabiensis, An. funestus). Each of these species has a large range that would obscure the other DVS present in Africa. Moreover, as these species tend to be the focus of most vector control efforts, it seemed prudent to provide a map that only indicates their distributions. However, this does not mean that the other DVS in Africa should be overlooked. For example, in the forested areas of western/central Africa An. moucheti is known to be a highly anthropophilic and efficient vector. Therefore a second map of Africa dedicated to showing the distributions of these 'secondary' DVS (in comparison to the 'top three') is also presented.
The number of presence data per species or species complex within the database included in the original model is given in Table 1. The global map indicating the distribution of 34 DVS is presented in Figure 1 (Figure 1 is also provided as a downloadable poster in Additional file 1). More detail is provided in the regional maps (Figures 2, 3, 4), which present a larger scale image of the information given in the global map. Figures 5 & 6 provide maps focused specifically on areas of high DVS diversity including Central America (Figure 5) and South-East Asia and Pacific Islands (Figure 6). Figure 7 indicates the secondary DVS found in Africa that were not included on the global map (Figure 1).
Figure 1 highlights the variability in the complexity of the malaria vector species communities and their distribution on a global scale. For example, comparing the Asian-Pacific region with Africa clearly demonstrates the highly complex and diverse nature of the DVS in Asia, whereas Africa shows a relatively simple picture, with the three main species co-dominating western (except the forested west of the Democratic Republic of Congo (DRC), Republic of Congo, Angola etc. where An. arabiensis is not found) and southeastern areas. The ability of An. arabiensis to utilise drier environments than An. gambiae or An. funestus is also clearly indicated, with the distribution of An. arabiensis extending farther north into the Sahel, east into Ethiopia, the southwestern corner of the Arabian Peninsula, Kenya and Somalia and south into the desert and steppe environments of Namibia and Botswana in southern Africa.
Across Europe and the Middle East, the most striking feature is the extent of the range of An. messeae that extends from the United Kingdom in the west across western and eastern Europe into Asia. This species is also the most northerly distributed of the 34 DVS.
North America (Figure 2) shows a very simple vector profile, with only An. freeborni found in the northwest and An. quadrimaculatus s.l. in the southeast with some minimal overlap with An. pseudopunctipennis in the very south of the continent. In South America, An. darlingi is shown as mainly dominant, but in Central America (where An. darlingi is also present but more focally distributed) its dominance is superseded by An. albimanus and An. pseudopunctipennis. Anopheles aquasalis is not a particularly efficient vector, but its ability to oviposit and utilise saline larval habitats means it remains the 'dominant' species in coastal areas of Central and South America.
An additional and important contribution of the current DVS map in South America is that it highlights the presence of An. marajoara, an important emerging vector in both the Guyana shield and the Amazon basin, as well as the distribution of the often overlooked An. albitarsis complex in the savanna ecoregion .
In Asia (Figure 4 & 6), the co-dominance of species is even more pronounced than found on the African continent, with the An. culicifacies complex, An. stephensi and the An. fluviatilis complex sympatric in India. The An. dirus and An. minimus complexes dominate together across much of Southeast Asia. However, along the Thai/Malaysian peninsula where these species both occur, they appear to have diverged and no longer occupy the same locations or ecological niche.
The Pacific islands (Figure 6) show a highly complex vector situation. On the island of New Guinea, members of the Punctulatus Group dominate, including the An. farauti complex, An. koliensis and the An. punctulatus complex, but only the An. farauti complex extends eastward to the Solomon Islands. Members of this complex are also found on the northern coast of Australia in Queensland and the Northern Territory.
In Indonesia, there appears to be high diversity and sympatry of vector species on the major islands. For example, in Sumatra, An. sinensis is found inland along with the An. barbirostris complex, An. leucosphyrus/An. latens and the An. minimus complex. A number of other species also exist on Sumatra, for example, the An. maculatus group and An. flavirostris, but none are considered dominant on the island; hence they are overlaid by the other, more dominant species. Alongside the An. sundaicus complex distributed along the coast, An. flavirostris does increase in relative 'dominance', by virtue of a reduced presence of other species, extending southward through Java until it is the only DVS found in the Lesser Sunda islands. In Sumatra, there is very little overlap amongst the dominant species found, suggesting that each occupies a separate niche on the island. Anopheles balabacensis dominates across most of Borneo, with some impact by the An. barbirostris complex and An. leucosphyrus/latens inland and the An. sundaicus complex on the coast.
In more northern areas of Asia (Figure 4), including China and Mongolia, An. sinensis and An. lesteri (syn. An. anthropophagus) are the only DVS. They appear to be sympatric in much of China and Korea, but this may be an artefact due to mis-identification of the two species in some areas.
The maps presented here show overlaid areas of species-specific predicted occurrence based on the climatic and environmental variables provided to the BRT model. Each species map included in the composite maps only included those pixels where the model predicted a probability of presence greater than 0.5. As with all species mapping, the quality of the output depends, for the most part, on the amount and quality of the data input into the model. Species occurrence data are often poorly distributed spatially [7–10] or are limited numerically (e.g. An. leucosphyrus/An. latens n = 12 (Table 1)). The modelling methodology allowed these data to be supplemented with randomly assigned (and therefore more spatially dispersed) pseudo-presence points taken from within the EO area of the species' range . These pseudo-data were weighted at half that of the 'true' occurrence data. However, where the occurrence data were limited, the pseudo-data may have exerted a greater influence on the final model, and therefore on the area of predicted presence. This can be seen in the predicted species occurrence on New Guinea Island. The EO ranges for An. farauti s.l., An. punctulatus s.l. and An. koliensis, indicate a blanket coverage across the whole island without considering the highland areas that run across almost the entire central length of the island. Members of the Punctulatus Group, which include An. farauti s.l., An. punctulatus s.l. and An. koliensis, are not known to occur at altitudes higher than 2300 m (Bangs, unpub obs) and the highlands on this island peak with Puncak Jaya (Mt. Carstensz) at 4884 m . The range of these three DVS centres on New Guinea island with limited spread to some of the other smaller neighbouring islands, and in the case of An. farauti s.l., to Northern Australia . This small range may have focussed the pseudo-presence points which may have fallen within both the lower and higher altitude locations, and thus the model was unable to establish altitude as a limiting factor for these species.
The quality of the occurrence data also relies on accurate species identifications reported in the source literature. The data were faithfully abstracted from each source and no assumptions were made, however this will have introduced some varying level of error. For example 'An. funestus' was rarely reported as a species complex, but also rarely subjected to the additional molecular methods of identification (e.g. Polymerase Chain Reaction (PCR)) [25, 26] necessary to identify accurately the members of the complex. Moreover, it is possible that some studies were actually reporting more than one member of the Funestus Group or Subgroup rather than An. funestus s.s. or even the An. funestus complex. The same may also be said for the An. maculatus group in Asia.
For some species there is also current debate about their taxonomy; for example, the identity and vectorial capacity of An. messeae is currently in question, with some suggestion that An. daciae may be responsible for malaria transmission previously attributed to An. messeae, and may be sympatric with An. messeae across much of its range, which might also explain the apparent high polymorphism associated with An. messeae[27, 28].
Despite some uncertainty in species classifications that cannot be corrected, the presence points for each species were carefully examined by the TAG, and those points that were clearly unreliable or related to dubious species identifications were removed at an early stage in the mapping process.
The maps presented here show the predicted occurrence of the DVS. They do not, however, indicate the probability of presence, although this information does underlie the distribution of positive and negative pixels (and is indicated on the original species maps [8–10]). A pixel is marked as 'present' where the BRT model indicated a probability of presence greater than 0.5. Therefore within these 'positives' the probability will range from > 0.5 to ≤ 1. Similarly, a pixel is marked as 'absent' where the BRT model indicated a probability of presence less than 0.5, but will include probabilities from 0 up to 0.5. These probability values are defined by the interaction of the environmental and climatic variables that are identified as predictors by the BRT model indicating where the environment is suitable for the species to exist. Hence such probabilities provide no direct information about potential species abundance but are simply the full output of the analysis. However, as these probabilities may indicate increasing or decreasing environmental suitability, it is feasible that these measures could be used to estimate species abundance at a specified location [29–31]. Further work is needed to try and establish a quantifiable link between these probabilities and DVS abundance.
Figure 1 provides the best currently available evidence-based global picture of the distributions of the main DVS. However, there will always be locations where the process has resulted in an oversimplification and the models do not pick up areas where a species may or may not be present. For example, in Africa, there is some question regarding the extensive predicted presence of An. funestus (species or members of the complex) within the highland areas of Ethiopia (Kiszewski, pers com). Indeed, elsewhere in the country, even where it is found, members of the Funestus Subgroup are rarely considered dominant, with An. arabiensis regarded as the major vector species . Only one known study has conducted PCR identification of Funestus Group specimens from Ethiopia and only reported An. parensis, a non-vector, as present .
A lack of data across a large swath of central Africa should also be noted, for example only 3 sites reporting DVS occurrence were found for the Central African Republic, 2 sites in Congo and 23 in DRC . Such areas therefore may not be being accurately represented by the model, especially where variable or unique environments and ecologies exist.
The large number of islands in the Asian-Pacific region, and those elsewhere of small size, can be problematic to accurately predict species occurrence. Overall, the models appear to have done well (based on TAG expert opinion), however there are a few cases where the model is not picking up areas of known presence. For example, on Grenada Island in the Americas the occurrence of An. pseudopunctipennis has been reported (see ) yet the model is not indicating a presence. However, An. aquasalis is correctly predicted to occur on this island. Similarly, An. barbirostris s.l. on the Lesser Sunda Island chain (including Flores, Sumbawa, Sumba, Timor and others) is not fully represented despite the existence of a published data point on Flores, and the islands being clearly within the EO range of this DVS. Anopheles barbirostris s.l. demonstrates dramatic varying behavioural attributes and vector importance over its geographical range in Indonesia, being of little or no epidemiological significance in Java and Sumatra in contrast to its role as a primary malaria vector in the eastern regions of the archipelago ([34, 35] Bangs, unpub obs), thereby illustrating some of the difficulties with certain species and the finer details for interpreting distribution maps.
The scale of these regional and global maps can also limit the visibility of some areas of presence on the smaller islands. For example, the Maluku Island chain in eastern Indonesia, where An. farauti s.l. is an important vector along the coasts and An. punctulatus s.l. a vector inland, does indicate the presence of these vectors, but mostly as sporadic individual pixels, and thus their presence is easy to overlook.
The maps presented here show the predicted distributions of a number of species complexes without reference to the sibling species they represent. Moreover, the molecular forms (M and S) of An. gambiae are not distinguished despite reported behavioural differences between them. This is due to a lack of spatially dispersed data providing accurate and defendable sibling species or form identification. It is hoped that such data will become increasingly available as the importance of correctly and fully identifying these species becomes more widely accepted, thus allowing for updated and detailed species-specific maps to be produced in the future.
Despite the known limitations and caveats given here, the global map represents the best currently available indication of the distributions of the dominant vectors of malaria. In line with the open access principles of the MAP, and via the ROADMAP initiative, the comprehensive dataset that was compiled to create the original species maps (currently available on request) will be directly available via the updated MAP website  from the beginning of 2012. The same site also hosts the individual species maps and bionomics information published in the earlier papers [8–10] and will also hold the maps published here in a format available for download.
The global and regional maps presented here have been created with a range of users in mind including researchers, modellers, public health officials and vector control managers. The maps provide a good level of basic location information and highlight the most important vectors present in a particular area. However, for many of these stakeholders to be able to make use of this information fully, for example, in deciding the relevant control efforts needed at a particular location, they also need information on how these species and species complexes behave at that particular location.
An ongoing project conducted by MAP as part of the VECNet consortium  aims to address this need. Quantified and geo-located data on biology, propensity to infection, behaviour and larval site characteristics are being compiled, initially for An. gambiae but expanding to include other important DVS. These data will be made available via an online tool that will use both data from the maps presented here and data on species behaviour at different locations to facilitate control decision-making.
Boosted Regression Trees
Dominant Vector Species
Malaria Atlas Project
Polymerase Chain Reaction
Repository for Open Access Data collected by the Malaria Atlas Project
Technical Advisory Group.
May JM: Map of the world distribution of malaria vectors. Geograph Review. 1951, 41: 638-639. 10.2307/210709.
Macdonald G: Local features of malaria. The epidemiology and control of malaria. 1957, London: Oxford University Press, 63-99.
Mouchet J, Carnevale P, Coosemans M, Julvez J, Manguin S, Richard-Lenoble D, Sircoulon J: Biodiversité du paludisme dans le monde. 2004, Montrouge, France: John Libbey Eurotext
Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J: A global index representing the stability of malaria transmission. Am J Trop Med Hyg. 2004, 70: 486-498.
Harbach RE: Genus ANOPHELES Meigen, 1818. Mosquito Taxonomic Inventory. Book Genus ANOPHELES Meigen, 1818 Mosquito Taxonomic Inventory City. 2011
Service MW, Townson H: The Anopheles vector. Essential Malariology. Edited by: Gilles HM, Warrell DA. 2002, London: Arnold, 59-84. Fourth
Hay SI, Sinka ME, Okara RM, Kabaria CW, Mbithi PM, Tago CC, Benz D, Gething PW, Howes RE, Patil AP, Temperley WH, Bangs MJ, Chareonviriyaphap T, Elyazar IR, Harbach RE, Hemingway J, Manguin S, Mbogo CM, Rubio-Palis Y, Godfray HC: Developing global maps of the dominant Anopheles vectors of human malaria. PLoS Med. 2010, 7: e1000209-10.1371/journal.pmed.1000209.
Sinka ME, Bangs MJ, Manguin S, Chareonviriyaphap T, Patil AP, Temperley WH, Gething PW, Elyazar IR, Kabaria CW, Harbach RE, Hay SI: The dominant Anopheles vectors of human malaria in the Asia-Pacific region: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2011, 4: 89-10.1186/1756-3305-4-89.
Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, Patil AP, Temperley WH, Gething PW, Kabaria CW, Okara RM, Van Boeckel T, Godfray HCJ, Harbach RE, Hay SI: The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010, 3: 117-10.1186/1756-3305-3-117.
Sinka ME, Rubio-Palis Y, Manguin S, Patil AP, Temperley WH, Gething PW, Van Boeckel T, Kabaria CW, Harbach RE, Hay SI: The dominant Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010, 3: 72-10.1186/1756-3305-3-72.
Malaria Atlas Project (MAP). [http://www.map.ox.ac.uk/]
White GB: Malaria. Geographical distribution of arthropod-borne diseases and their principal vectors Geneva: World Health Organization, Division of Vector Biology and Control. 1989, 7-22. WHO/VBC/89967
Service MW: The Anopheles vector. Bruce-Chwatt's Essential Malariology. Edited by: Gilles HM, Warrell DA. 1993, London: Edward Arnold, 96-123. Third
Service MW: Appendix II. Characteristics of some major Anopheles vectors of human malaria. Bruce-Chwatt's Essential Malariology. Edited by: Gilles HM, Warrell DA. 1993, London: Edward Arnold, 305-310. Third
Manguin S, Carnevale P, Mouchet J, Coosemans M, Julvez J, Richard-Lenoble D, Sircoulon J: Biodiversity of malaria in the world. 2008, Montrouge, France: John Libbey Eurotext
Elith J, Leathwick JR, Hastie T: A working guide to boosted regression trees. J Anim Ecol. 2008, 77: 802-813. 10.1111/j.1365-2656.2008.01390.x.
ArcMAP (ArcGIS) ESRI. [http://www.esri.com/products/index.html]
Reid JA: Anopheline mosquitoes of Malaya and Borneo. 1968, Malaysia
Harrison BA: Medical entomology studies - XIII. The Myzomyia Series of Anopheles (Cellia) in Thailand, with emphasis on intra-interspecific variations (Diptera: Culicidae). Contrib Am Entomol Inst (Ann Arbor). 1980, 17: 1-195.
Rao TR: The anophelines of India. 1984, New Delhi: Malaria Research Centre, Indian Council of Medical Research, xvi + 518, 2
Trung HD, Bortel WV, Sochantha T, Keokenchanh K, Briet OJ, Coosemans M: Behavioural heterogeneity of Anopheles species in ecologically different localities in Southeast Asia: a challenge for vector control. Trop Med Int Health. 2005, 10: 251-262. 10.1111/j.1365-3156.2004.01378.x.
Parida SK, Hazra RK, Marai N, Tripathy HK, Mahapatra N: Host feeding patterns of malaria vectors of Orissa, India. J Am Mosq Control Assoc. 2006, 22: 629-634. 10.2987/8756-971X(2006)22[629:HFPOMV]2.0.CO;2.
Rubio-Palis Y, Zimmerman RH: Ecoregional classification of malaria vectors in the neotropics. J Med Entomol. 1997, 34: 499-510.
Koekemoer LL, Kamau L, Hunt RH, Coetzee M: Cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002, 66: 804-811.
Cohuet A, Simard F, Toto JC, Kengne P, Coetzee M, Fontenille D: Species identification within the Anopheles funestus group of malaria vectors in Cameroon and evidence for a new species. Am J Trop Med Hyg. 2003, 69: 200-205.
Di Luca M, Boccolini D, Marinuccil M, Romi R: Intrapopulation polymorphism in Anopheles messeae (An. maculipennis complex) inferred by molecular analysis. J Med Entomol. 2004, 41: 582-586. 10.1603/0022-2585-41.4.582.
Linton YM, Lee AS, Curtis C: Discovery of a third member of the Maculipennis Group in SW England. Euro Mosq Bull. 2005, 19: 5-9.
Hutchinson GE: Population Studies - Animal Ecology and Demography - Concluding Remarks. Cold Spring Harbor Symposia on Quantitative Biology. 1957, 22: 415-427. 10.1101/SQB.1957.022.01.039.
Kulhanek SA, Leung B, Ricciardi A: Using ecological niche models to predict the abundance and impact of invasive species: application to the common carp. Ecol Appl. 2011, 21: 203-213. 10.1890/09-1639.1.
VanDerWal J, Shoo LP, Johnson CN, Williams SE: Abundance and the Environmental Niche: Environmental Suitability Estimated from Niche Models Predicts the Upper Limit of Local Abundance. Am Nat. 2009, 174: 282-291. 10.1086/600087.
Nyanjom SR, Chen H, Gebre-Michael T, Bekele E, Shililu J, Githure J, Beier JC, Yan G: Population genetic structure of Anopheles arabiensis mosquitoes in Ethiopia and Eritrea. J Hered. 2003, 94: 457-463. 10.1093/jhered/esg100.
Weeto MM, Koekemoer LL, Kamau L, Hunt RH, Coetzee M: Evaluation of a species-specific PCR assay for the Anopheles funestus group from eleven African countries and Madagascar. Trans R Soc Trop Med Hyg. 2004, 98: 142-147. 10.1016/S0035-9203(03)00019-1.
Reid JA: The Anopheles barbirostris group (Diptera: Culicidae). Bull Entomol Res. 1962, 53: 1-57. 10.1017/S0007485300047945.
Reid JA, Harrison BA, Atmosoedjono S: Variation and vector status in Anopheles barbirostris. Mosq Syst. 1979, 11: 235-251.
We wish to thank Robi Okara, Rosalind Howes, Edward Haynes, Philip Mbithi, Owen Yang, Carolynn Tago and Elisabeth Thiveyrat for primary data abstraction and to Katherine Battle and David Pigott for proof reading the manuscript. We are very grateful for insightful comments and encouragement from Anthony Kiszewski on the global and regional composite maps. We also thank unreservedly our Technical Advisory Group for their extended support over the duration of the project. The initial vector database and species distribution maps were created while MES was funded by a project grant from the Wellcome Trust (#083534) to SIH. The current work was conducted with support and funding from the Bill and Melinda Gates Foundation via the VECNet consortium http://www.vecnet.org/, which also provides support for TRB. SIH is funded by a Senior Research Fellowship from the Wellcome Trust (#079091), which also supported CWK and supports PWG. APP and WHT are funded by a Wellcome Trust Biomedical Resources Grant (#091835). SIH and PWG also acknowledge support from the RAPIDD program of the Science & Technology Directorate, Department of Homeland Security, and the Fogarty International Center, National Institutes of Health http://www.fic.nih.gov. This work forms part of the output of the Malaria Atlas Project (MAP, http://www.map.ox.ac.uk), principally funded by the Wellcome Trust, U.K. MAP also acknowledges the support of the Global Fund to fight AIDS, Tuberculosis, and Malaria http://www.theglobalfund.org and grants from the University of Oxford--Li Ka Shing Foundation Global Health Program and the Oxford Tropical Network. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that they have no competing interests.
SIH conceived the study. MES compiled the global and regional multi-species maps, wrote the first draft of the manuscript and assembled the occurrence data, with assistance from CWK and others (see acknowledgements), and ran the models to create the individual DVS species maps. CWK also digitised and edited all the expert opinion maps. WHT designed and maintained the databases. APP implemented the BRT scripts for predictive mapping. PWG processed the environmental and climatic data grids. TAG members (MJB, SM, YR-P, REH, TC, MC, CMM and JH) provided data and advice in updating the EO range maps and in ranking the importance of the DVS per country. TRB provided additional comments on the maps. All authors participated in the writing and editing of the manuscript and read and approved the final submitted version. REH also served as the taxonomy advisor/consultant to the TAG and throughout this vector mapping project.