Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Malaria in Africa: Vector Species' Niche Models and Relative Risk Maps

  • Alexander Moffett,

    Affiliation Section of Integrative Biology, University of Texas at Austin, Austin, Texas, United States of America

  • Nancy Shackelford,

    Affiliation Section of Integrative Biology, University of Texas at Austin, Austin, Texas, United States of America

  • Sahotra Sarkar

    To whom correspondence should be addressed. E-mail: sarkar@mail.utexas.edu

    Affiliation Section of Integrative Biology, University of Texas at Austin, Austin, Texas, United States of America

Abstract

A central theoretical goal of epidemiology is the construction of spatial models of disease prevalence and risk, including maps for the potential spread of infectious disease. We provide three continent-wide maps representing the relative risk of malaria in Africa based on ecological niche models of vector species and risk analysis at a spatial resolution of 1 arc-minute (9 185 275 cells of approximately 4 sq km). Using a maximum entropy method we construct niche models for 10 malaria vector species based on species occurrence records since 1980, 19 climatic variables, altitude, and land cover data (in 14 classes). For seven vectors (Anopheles coustani, A. funestus, A. melas, A. merus, A. moucheti, A. nili, and A. paludis) these are the first published niche models. We predict that Central Africa has poor habitat for both A. arabiensis and A. gambiae, and that A. quadriannulatus and A. arabiensis have restricted habitats in Southern Africa as claimed by field experts in criticism of previous models. The results of the niche models are incorporated into three relative risk models which assume different ecological interactions between vector species. The “additive” model assumes no interaction; the “minimax” model assumes maximum relative risk due to any vector in a cell; and the “competitive exclusion” model assumes the relative risk that arises from the most suitable vector for a cell. All models include variable anthrophilicity of vectors and spatial variation in human population density. Relative risk maps are produced from these models. All models predict that human population density is the critical factor determining malaria risk. Our method of constructing relative risk maps is equally general. We discuss the limits of the relative risk maps reported here, and the additional data that are required for their improvement. The protocol developed here can be used for any other vector-borne disease.

Introduction

A central theoretical goal of epidemiology is the construction of spatial models of disease prevalence and risk, including maps for the potential spread of infectious disease [1], [2]. In particular, Snow et al. [3] have emphasized the need for risk maps for malaria in Africa which accounts for an estimated 85% of the 1 million annual deaths due to this disease [4]. Risk maps can be used to identify appropriate strategies of response to disease outbreaks including targeted vaccination [5] and vector, reservoir, or agent control [6]. Risk maps have been constructed using a variety of techniques including reports of disease cases [7] and distributions of disease agents, reservoirs, or vectors, based on surveys and expert opinion [8]. In recent years, these methods have been extended to use ecological models to predict the potential spatial spread of disease [9]. The use of ecological niche models for quantitative prediction of geographical distributions of agent, reservoir, and vector species has been advocated to augment traditional mapping methods such as splining and kriging [9]. Here we report a systematic attempt to construct niche models for all vectors of malaria in Africa for which data are available.

Niche models predict the “fundamental niche” of a species which identifies the region in ecological space that the species would occupy were its movement unrestricted [9], [10]. Within biodiversity studies, in which the concern is typically with the conservation of the extant habitat of a species, this fundamental niche (representing a potential distribution) is typically restricted to a realized geographical niche by using additional information such as confirmed occurrence records for a species within each accepted contiguous piece of modelled suitable habitat [10]. Such a choice is conservative in the sense that a species may occur outside its predicted realized niche. In typical contexts of biodiversity conservation planning such conservatism is appropriate; it ensures that areas selected for conservation management are maximally likely to contain the species predicted for them. However, in epidemiological contexts, especially if the interest is in identifying risk, the geographical extent of the fundamental niche is more relevant as this range defines the areas to which agent, reservoir, and vector species may potentially spread because of ecological suitability. (This means that, in epidemiological contexts, the fundamental niche or potential distribution should not be clipped to a smaller realized niche in most cases.) In this analysis, fundamental niches as predicted by niche models are used as the basis for the construction of relative risk maps.

We construct niche models for 10 Anopheles species recognized as vectors of malaria in Africa using a maximum entropy method based on known species' occurrences and environmental layers. These niche models predict geographic distributions of the species. For seven of these species the results presented here appear to be the first niche models reported in the literature. We used the Maxent software package [11] for the maximum entropy modelling. Within biodiversity studies, in which niche modelling is a standard technique [12][16] this maximum entropy method has emerged as one of the three most reliable techniques for predicting species' distributions [11], [17]. The other two most reliable methods are genetic algorithms (GARP [18]) and regression trees [19]. The advantage of Maxent over GARP is that it is much faster and allows for the simultaneous modelling of an indefinite number of species. The advantage of maximum entropy methods over both genetic algorithms and regression methods is that Maxent predicts relative probabilities of occurrence rather than the simple presence or absence of a species. This permits a finer (more nuanced) risk assessment than what can be achieved from presence-absence predictions alone. Our analysis appears to be the first use of Maxent in an epidemiological context.

Using the niche models, spatial information on human population densities, and the human blood index (HBI) values of Anopheles species for which these values were available, we construct preliminary relative risk maps for malaria in Africa. These maps report the relative risk of malaria occurrence at different geographical locations. We discuss in detail what data are necessary to make such maps more accurate and how our methods can be generalized to other vector-borne diseases.

There have been five recent vector-based attempts to construct risk maps for malaria. Kiszewski et al. [8] produced a global risk map based on the stability of malaria parasite transmission. They divided the world into 260 malarious regions and then identified the dominant vector in each region. The dominant vector was defined as the most abundant Anopheles species in the region that was a malaria vector, contained sporozoites frequently, and fed predominantly on humans. The stability index of a region was defined using the human feeding rate, daily survival rate, and length of the extrinsic incubation period of the dominant vector. The effects of temperature and precipitation on these parameters were quantified and environmental data were used to produce a worldwide projection of malaria transmission.

Lindsay et al. [20] used nonlinear regression to relate known occurrences of A. arabiensis and A. gambiae to environmental parameters. The regression was used to predict the relative proportion of the two species throughout Africa. Kuhn et al. [21] based their map on a database of occurrence records of six Anopheles species in Europe. Statistical analyses correlated these occurrence points with environmental parameters. The correlations were used to predict the distribution of the Anopheles species across Europe. Rogers et al. [22] used a maximum likelihood analysis to identify the optimal environmental parameters to predict the occurrence of five Anopheles species in Africa. Levine et al. [23], [24] used niche models to map malaria risk in Africa, Central America, and the United States. Using species' occurrence and environmental data in a genetic algorithm (in the GARP software package), they predicted the distribution of three Anopheles species in Africa and five species in the United States.

In addition to these vector-based approaches, risk maps have also been constructed by mapping the distribution of the malaria parasite. Kleinschmidt et al. [25], [26] used regression analysis to determine the relationship between the malaria parasite prevalence of an area and environmental variables. The regression model was refined using kriging for spatial interpolation thus producing a map of malaria prevalence. Diggle et al. [27] used a generalized linear mixed model to determine the relationship between malaria presence in children and their age and bed net use along with available medical services and land cover. Gemperli et al. [7] developed a Bayesian model to calculate parasite transmission intensity on the basis of malaria survey information and land cover, temperature, and rainfall data.

We map malaria risk using distributions of malaria vectors rather than malaria parasites. The underlying assumption is that malaria vector abundance is an adequate surrogate for malaria risk without explicit inclusion of parasite abundance. However, even if this assumption is invalid, risk maps based on vector distributions are still useful to augment risk analyses based only on parasite distributions. (In future work, we plan to construct risk maps using more sophisticated models of transmission that include interactions between parasites, vectors, and variable human susceptibility [28][30].)

This analysis differs from and extends previous efforts in five ways. First, our niche models are constructed using a maximum entropy method. Second, we model a larger number of species than previously attempted. Third, we include human population distribution data in the risk analysis and associated relative risk models. Fourth, we perform a sensitivity analysis to evaluate the robustness of our results. Fifth, we show how additional data can lead to a more sophisticated risk analysis for malaria and other vector-borne diseases.

Methods

Data

A 1 arc-minute (0.01666°×0.01666° longitude×latitude, approximately 4 km2 at the equator) grid was used to divide Africa into 9 185 275 cells. Twenty environmental layers were obtained from the WorldClim database [31]. Each layer was available at a resolution of 30 arc-seconds (0.008333°×0.008333°) and was resampled at a resolution of 1 arc-minute. These layers are listed in Table 1. (All spatial data manipulation used ArcMap GIS [32].)

thumbnail
Table 1. Environmental Parameters Used in Niche Modelling

https://doi.org/10.1371/journal.pone.0000824.t001

Land cover data were obtained from the Global Land Cover Facility [33]. AVHRR satellite data acquired between 1981 and 1994 were used to derive 14 classes of land cover: water; evergreen needleleaf forest; evergreen broadleaf forest; deciduous needleleaf forest; deciduous broadleaf forest; mixed forest; woodland; wooded grassland; closed shrubland; open shrubland; grassland; cropland; bare ground; urban and built. The data were initially available at a resolution of 1 km2 and resampled at a resolution of 1 arc-minute.

A review of the global distribution of malaria vectors [34] was used to identify those African Anopheles species that are capable of transmitting malaria. The 29 species used for this analysis are listed in Table 2. An extensive literature search was performed to obtain records of vector occurrences. Besides using the Mapping Malaria Risk in Africa (MARA) database [3] both PubMed and Google Scholar were searched using “Africa” and “distribution” in conjunction with the names of each of the Anopheles species listed in Table 2. References from those papers so identified were also searched. This resulted in a data set of 3 342 records of 22 malaria vectors, with 2634 of the records drawn from the MARA database.

Each record was georeferenced to the nearest arc minute and assigned to a corresponding cell. For those cells containing more than one record of a given species only the most recent record was kept. So as to increase the likelihood that current values of the selected environmental parameters represent the environment as it existed when the records were obtained, only those records reporting observations made after 1980 were included. The resulting data set consisted of 977 records of 22 Anopheles species with 367 of these records drawn from the MARA database. Table 2 provides a summary of these records.

Human population density data for the year 2000 were obtained from the Gridded Population of the World database [35]. These data were provided at a resolution of 2.5 arc-minutes (0.041666°×0.041666°) and resampled at a resolution of 1 arc-minute. Figure 1 shows a map of normalized population densities.

thumbnail
Figure 1. Population density in Africa.

The population densities have been normalized so as to range over the unit interval.

https://doi.org/10.1371/journal.pone.0000824.g001

HBI values were obtained from the literature. Both PubMed and Google Scholar were searched using the terms “Anopheles” and “human blood index” with each of the species names. References from those papers so identified were also searched. Of the 10 species for which we were able to construct reliable niche models (see below), HBI values were available for nine species, with no value available for A. paludis.

Niche Models

Previous niche models of Anopheles species have used genetic algorithms, as implemented in GARP [23], [24]. In this analysis we use instead a maximum entropy technique implemented in the Maxent software package [11]. Maxent was used instead of GARP because Maxent provides relative probabilities of presence rather than only presence-absence output and because it has, in some recent studies, been shown to outperform GARP with respect to its predictive success [11], [17]. However, there are some preliminary data suggesting that GARP may perform better than Maxent at extrapolating from occurrence data (A. Townsend Peterson, personal communication). This means that models constructed using Maxent may be excessively conservative thus predicting false absences more often than GARP. If this is true then while the areas of high relative risk identified in our maps are probably reliable, those with low relative risk should not be entirely discounted. Regression tree methods (for instance those implemented in the RandomForest software package [19]) have so far not been used in an epidemiological context.

Given a set of records of species occurrences and values of selected environmental variables defined over a chosen geographical region, Maxent predicts the distribution of species in that space by finding the distribution of maximum entropy subject to the constraint that the expectation of the distribution of each species matches its observed average over the sample locations [11]. The distribution thus produced is a relative probability distribution for each species over all cells. If each relative probability in a cell is divided by the maximum such probability across the landscape, then the result is a normalized relative probability distribution which assumes that the species is certainly present in the cell in which it had its highest predicted relative probability. This analysis used this normalization.

Computer memory limitations prevented Maxent runs using as input all 9 185 175 cells simultaneously. Instead, Maxent was run using 100 sets of 10 000 cells drawn randomly from the complete set. Sets of 10 000 cells were used because Maxent performance does not significantly improve when more than 10 000 cells are used (S. Phillips, personal communication), while 100 sets were selected to sample widely from the complete set of cells. The 100 different niche models produced for each species were averaged and subsequently normalized as described above. As we were concerned with the geographical range of the fundamental niche of the vector species, these output maps were not further refined.

The accuracy of the niche models thus produced was evaluated by constructing the models using only 75% of the records with the remaining 25% set aside for testing. The accuracy of each model was then determined by performing both a threshold-dependent binomial test of omission and a threshold-independent receiver operating characteristic analysis with those cells set aside during model development [11].

In the threshold-dependent binomial test of omission a threshold of 0.10 was used to classify each vector as either present or absent in each cell, with a vector present in a given cell if the niche model assigned it a value greater than 0.10 in that cell and absent otherwise. This threshold transformed the continuous data produced by the niche models into binary data to allow a binomial test to be performed. For each vector the number of cells in which the vector was predicted to be present was compared to the number of cells known to contain the vector. A one tailed binomial test was used to determine whether the niche model outperformed a random model predicting the vector to be present in the same number of cells [11]. Maxent provides test statistics for binomial tests for 10 different threshold values. The value of 0.10 used in this analysis was arbitrarily selected from this set of possible threshold values.

In the threshold-independent receiver operating characteristic analysis, the sensitivity and specificity of the niche models were calculated at all possible thresholds. The sensitivity of a model at a threshold was defined as the percentage of species occurrences that were correctly predicted by the model at the threshold, while the specificity of the model at a threshold was defined as the percentage of correctly predicted species absences at the threshold [11]. By calculating the sensitivity and specificity of a model at all possible thresholds a receiver operating characteristic curve was produced with sensitivity plotted on the y-axis and (1–specificity) plotted on the x-axis. The area under the curve (AUC) of the resulting plot provides a measure of model performance independent of the choice of any particular threshold. An optimal model, one that predicted each occurrence of a species and for which each prediction was accurate, would have an AUC of 1.0 while a model that predicted species occurrences at random would have an AUC of 0.5.

These two tests were used to restrict attention to those models that performed significantly better than random. Only those niche models possessing both a p value less than 0.05 for the binomial test of omission and an AUC greater than 0.75 were used in this analysis. The same protocol for model retention has previously been used by Pawar et al. [16].

Risk Models

Three different relative risk models were constructed in which a value between 0 and 1 was assigned to each cell representing the relative risk of malaria posed to the human population residing within it. These models only used the nine species for which an HBI value was available. The models only incorporate the risk from ecological and demographic factors and ignore the modulation of risk through human intervention such as measures to control the spread of parasites or vectors. They also assume that parasites are present at sufficient densities to be capable of spreading whenever vectors are present. These assumptions are generally appropriate for Africa given the continued prevalence of malaria within it. However they would not be appropriate for regions such as northern Australia from which malaria parasites have been eliminated though malaria vectors remain

Let alik be the relative abundance of the l-th vector for the i-th parasite in cell k (relative to other cells in the landscape). Let pk be the human population of cell k. Let hlik be the HBI of the l-th vector for the i-th parasite in cell k. Let elik be the transmission efficiency of the l-th vector for the i-th parasite in cell k, measured by the relative likelihood of parasites being transferred to the human agent with each bite. We use a simple multiplicative model for the relative risk, φlik, due to the l-th vector for the i-th parasite in cell k:(1)Constraints on available data force further simplification of the model: (i) differences between the two primary parasites for malaria in Africa (Plasmodium vivax and P. falciparum) were ignored; (ii) the HBI was interpreted as an intrinsic property of the vector species that does not vary over geographical space; and (iii) differences in transmission efficiency between vectors were ignored. This results in the simplified relative risk model:(2)This equation expresses the relative number of individuals in cell k to whom the malaria parasite is expected to be transmitted by the l-th vector. Finally, a linear correlation is assumed between the expected relative abundance of a vector and its relative probability of presence as predicted by its niche model. This assumption must be tested in the field. Because of these four assumptions, the relative risk model presented here must be regarded as very preliminary and treated with caution.

To construct a relative risk map from this model, the relative risk from the different vectors must be compounded for each cell in the landscape. Three different models for compounding relative risk were used for this purpose:

(i) Additive model.

The relative risk due to the different vectors were added together and then normalized on a scale of 0 to 1. Let Φu be the relative risk of cell k. Then:(3)This model thus assumes that there is no interaction between vector species, that each vector is able to inhabit the full extent of its fundamental niche and that there is no competition between vectors for human blood meals. If there is no such interaction then it should be possible to add the relative risks posed by the vectors in arriving at the overall relative risk posed to a given cell. However, the presumed lack of interaction between vectors has been questioned [36] and probably does not hold for all vector species. Moreover, given the methodology of this paper, there is an additional problem with the additive model. The relative probabilities provided by Maxent sum to 1 over the landscape. Thus, even if the relative abundances of a vector in the cells are linearly correlated with these probabilities, there is no way to estimate the absolute abundance. Adding the relative risk values makes the assumption that the highest absolute abundance of a vector among the cells of the landscape is the same for each vector. This assumption is ecologically suspect. Consequently, the next two models for compounding relative risk are more plausible than the additive model.

(ii) Minimax model.

The relative risk of each cell was defined as the maximum relative risk from any one of the vectors. Thus:(4)This is being called a minimax model since, ultimately, in epidemiology, the goal is to minimize the risk of disease while what is being used as a relative risk measure is the maximum risk associated with all the vectors in a given cell.

(iii) Competitive exclusion model.

The relative risk associated with the vector of highest relative abundance is identified with the relative risk of the cell. Let . Then:(5)Since the relative abundances were identified with relative probabilities of occurrence, this model assumes that the vector that has the highest predicted probability of occurrence will displace all others. This assumption is similar to that about a dominant vector incorporated in the risk maps of Kiszewski et al. [8]. However, as studies of malaria transmission have identified the presence of more than one significant malaria vector within some regions [37], [38] there is again reason to doubt the validity of this assumption in many contexts.

Sensitivity Analysis

A sensitivity analysis was performed to determine the robustness of our results to variation in the values of the parameters included in the models. For each of the 100 niche models produced for each of the Anopheles species, 10 sets of HBI values were obtained by drawing randomly from a uniform distribution over the interval defined by the minimum and maximum measured HBI values for each species. For each of the three models, 1 000 relative risk values were produced for each cell using the 100 niche models, the random HBI values, and the actual human population densities. The standard deviation of the relative risk values of each cell was calculated and used to identify areas of high sensitivity.

Results

The need to set aside 25% of the records of each Anopheles species for testing purposes resulted in the production of niche models for only those 17 (out of 29) species for which four or more records were available. Data representing the accuracy of the niche models are provided in Table 3. Of these 17 niche models all but one possessed an AUC greater than 0.75. However, only 10 of the niche models possessed a p value less than 0.05 at a threshold of 0.10. As would be expected, the p values of the niche models were closely correlated with the number of records upon which the models were based. Niche models with an AUC greater than 0.75 and a p value less than 0.05 were produced for A. arabiensis, A. coustani, A. funestus, A. gambiae, A. melas, A. merus, A. moucheti, A. nili, and A. quadriannulatus. Figure 2 presents maps of the predicted geographical distributions of these species. Darker regions are those of greater relative probability of occurrence while lighter areas are those in which the relative probability of occurrence is small.

thumbnail
Figure 2. The distributions of 10 malaria vectors in Africa.

Distributions are provided for: (a) A. arabiensis; (b) A. coustani; (c) A. funestus; (d) A. gambiae; (e) A. melas; (f) A. merus; (g) A. moucheti; (h) A. nili; (i) A. paludis; (j) A. quadriannulatus.

https://doi.org/10.1371/journal.pone.0000824.g002

Two different tests were performed to determine the contributions of each of the environmental parameters to the niche models. In the first test, the AUC of each niche model was calculated using each of the environmental parameters individually. Those parameters that resulted in the highest AUC were interpreted as those which possessed the most information regarding the niche of a species. In the second test, the AUC of each niche model was calculated after omitting each of the environmental parameters one at a time. The effects of omitting each parameter were determined by comparing the resulting AUC to the actual AUC of the model. Those parameters for which the difference between these values was highest were interpreted as possessing the most information not present in the other environmental parameters.

The results of these two tests are provided in Table 4. As can be seen, there was no clear pattern in the contributions of the environmental parameters. The niches of different species appeared to be determined by different parameters. Over the entire data set of ten species, none of the parameters appeared to be significantly more important than any other parameter.

Data representing the HBI values of the Anopheles species are provided in Table 5. HBI values were available for 9 of the 10 species for which niche models were produced. An HBI value for A. plaudis could not be found in the literature. The HBI values of the remaining species were observed to vary significantly, with mean values ranging between 0.011 for A. quadriannulatus and 1.00 for A. merus. The values averaged to obtain the mean HBI value were likewise found to vary for different species, with standard deviations as high as 0.269 for A. melas and 0.241 for A. arabiensis.

Figure 3 depicts the relative risk maps constructed using the additive, minimax, and competitive exclusion models. Cells in each map were assigned a value between 0 and 1, with 0 representing no risk and 1 representing maximal relative risk. The resulting relative risk values were found to congregate closely to either 0 or 1. To ease the visual discernment of the relative risk faced in different regions, the maps in Figure 3 plot the natural logarithm of the relative risk values.

thumbnail
Figure 3. The distribution of malaria relative risk in Africa.

Three different types of risk were calculated as follows: (a) the probability of occurrence of each vector in each cell was multiplied by both the human population density of the cell and the HBI of the vector. The relative risk of malaria in the cell was calculated as the sum of these values; (b) the vector possessing the maximum probability of occurrence was identified for each cell. Its probability of occurrence was multiplied by its HBI and the human population density of the cell. The relative risk of malaria in the cell was calculated as the product of these three values; and (c) the probability of occurrence of each vector in each cell was multiplied by the human population density of the cell and the HBI of the vector. The relative risk of malaria in the cell was calculated as the maximum of these values. The maps plot the natural logarithm of the relative risk.

https://doi.org/10.1371/journal.pone.0000824.g003

In general the relative risk maps produced by the three relative risk models were quite similar. As can be seen in Figure 3, the maps produced from the minimax and competitive exclusion models were nearly identical. These maps differed slightly from the relative risk map produced using the additive model as the relative risk in the additive model was more closely restricted to the areas of high population density. As can be seen in comparing Figure 1 with Figure 3, population density appears to have been the primary determinant of the relative risk of malaria in each of the three models. Areas of high human population density were found to be those in which the relative risk of malaria was greatest. The sole exception to this was in North Africa in which high human population densities were accompanied by relatively low relative risk values.

The sensitivity of the additive model was observed to be less than that of the competitive and minimax models. The standard deviations in the relative risk values ranged from 0 to 0.144 for the additive model to 0 to 0.280 for the competitive model and 0 to 0.279 for minimax model. In each model the majority of cells possessed a standard deviation less than 0.001. In the additive model 98.9% of cells had a standard deviation less than 0.001, while the percentages for the competitive and minimax models were 97.7% and 98.5%, respectively. Most cells were thus robust to variation in the parameter values. There was no clear pattern observed in the distribution of sensitivity across the landscape. While relative risk was concentrated in areas of high human population density, sensitivity was not.

Discussion

As explained below, our relative risk maps are preliminary and must be treated with due caution. However, the methodology developed here can be used to construct relative risk maps for other infectious diseases, with the risk models modified to reflect the presence of multiple agents and reservoirs, more complicated modes of transmission, and other relevant ecological factors. Because, as we and others have shown [9], [11], [16], [17], niche models of acceptable accuracy based on climatic and topographic parameters (as determined by the internal tests of software packages as well as from comparisons of maps produced by experts—see below for more discussion) can be constructed with sparse data, this methodology shows promise for the construction of relative risk maps for infectious diseases even when data are limited and not easily collected. Some suggestive recent work argues that, besides climatic and topographic factors, land use and land cover change influence the spread of malaria [39], [40]. The models presented here used land cover as one of the explanatory variables in the niche models; in future work we plan to include land use and to explore the effects of land cover in more detail, using a finer classification of land cover types.

A variety of different climatic and topographical factors were important for predicting the distributions of different species. Levine et al. [23] report similar variable results for the three Anopheles species that they model. These results have the implication that, to predict malaria vector species distributions, simple ecological heuristic rules (for instance, those incorporating optimal precipitation and/or temperature) will likely be unreliable and should be replaced with predictions of niche models based on as complete a set of environmental parameters as possible. Rogers and Randolph [41] have also noted the lack of relatively simple heuristic rules.

A perhaps not unexpected result is that the most important determinant of relative risk for malaria was human population density, assuming as we do here for Africa that (i) ecological and demographic factors rather than control measures are the determinants of risk and (ii) parasite densities are sufficiently high for disease spread. The extent of the effect of population density in this analysis is likely a result of the multiplicative nature of our risk models. However, this effect cannot be regarded solely as an artifact of the admitted simplicity of these models as any plausible model of infectious disease transmission should include the human population as a multiplicative factor. This dependence on population density has the implication that the expected local increase of population density due to increased urbanization [42], especially in Africa [43], will increasingly exacerbate the risk of malaria unless control measures are implemented [39].

Turning to the details of our results, with a few notable exceptions, the distributions for the niche models in Figure 2 closely follow the previously proposed distributions of most modeled species. This can be seen by comparing by region the vector distributions provided in Figure 2 with those previously presented in the literature.

The niche models of vectors within West Africa appear to coincide quite closely with the expert-based distributions of the modeled vectors within this region. Haworth [34] identified A. arabiensis, A. coustani, A. funestus, A. gambiae, and A. nili as primary and secondary malaria vectors in West Africa. This is consistent with the niche models of Figure 2 in which these three species are shown to possess extremely high relative probabilities of occurrence within this region. Haworth [34] also identified A. melas, A. moucheti, and A. paludis as vectors with limited presence in West Africa. This is consistent with the models in Figure 2 in which these species have limited distributions in parts of that region. In addition, Gillies and Coetzee [44] presented this region as containing A. arabiensis, A. gambiae, and A. melas, while White [45] reported the region as containing A. arabiensis, A. gambiae, A. melas, and A. funestus. The distributions of these species within this region provided by these two sources were thus consistent with the niche models of Figure 2.

In Southeast Africa and Madagascar the niche models again correspond closely to expert-based distributions. Haworth [34] identified A. arabiensis, A. coustani, A. funestus, A. gambiae, A. merus, and A. paludis as the primary and secondary malaria vectors within the region. In Figure 2 these species have high relative probabilities of occurrence in the region. Moreover, the distributions within this region of A. arabiensis, A. gambiae, and A. merus, as presented by Gillies and Coetzee [44], and the distributions of A. arabiensis, A. gambiae, A. funestus, and A. merus, as presented by White [45], were consistent with the niche models in Figure 2.

While the correspondence between the niche models produced in this analysis and the expert-based distributions of the Anopheles species was quite strong in most regions, discrepancies between these distributions are found in Central Africa. Gillies and Coetzee [44] presented the distributions of A. arabiensis and A. gambiae as stretching across Central Africa. However, these distributions conflict with the niche models of Figure 2, in which both vectors have low probabilities of occurrence throughout this region. In slight contrast to Gillies and Coetzee, White [45] depicted both A. gambiae and A. funestus as present in Central Africa, while depicting A. arabiensis as largely absent from the region. While this distribution of A. arabiensis agrees with the niche models of Figure 2, the proposed distributions of A. gambiae and A. funestus do not. Levine et al. [23] found both A. arabiensis and A. gambiae to be distributed throughout Central Africa, unlike what is seen in Figure 2. Similar results were obtained by Lindsay et al. [20] who found the climate throughout Central Africa to be suitable for both of these species. Rogers et al. [22] obtained results similar to those of both Levine et al. and Lindsay et al. with respect to the distribution of A. gambiae, yet found A. arabiensis to be largely absent from Central Africa (which thus agrees with Figure 2).

These discrepancies may simply be a result of the general lack of distributional data for Anopheles species in Central Africa. Fewer than 10 occurrence records within this region were available for use in the present analysis, with the other analyses likewise lacking much data for the region. If the niche models produced by Maxent are excessively conservative, then the low probability of occurrence associated with A. gambiae in Central Africa may be a consequence of the lack of any occurrence data for the species in this region. However, given the current absence of records of Anopheles species within Central Africa, the attribution of a low relative probability of occurrence to both A. arabiensis and A. gambiae within the region may be a correct prediction. The high AUC values associated with the niche models produced in this analysis support this claim. Field data from Central Africa are required to resolve the discrepancies between our predictions and other distributional maps..

Some other discrepancies appear to indicate that our models are an improvement over previous proposed distributions. For example, Coetzee [46] criticized the niche models of Levine et al. [23] for predicting the occurrence of A. quadriannulatus in areas of central Botswana that are supposed to be too arid for the species. She also questioned proposed distributions of Lindsay et al. [20] for including A. arabiensis in arid parts of South Africa. In contrast with these distributions, the niche models in Figure 2 predict small relative probabilities of occurrence for A. quadriannulatus and A. arabiensis in central Botswana and South Africa, respectively. While our results thus avoid these criticisms, since Levine et al.'s criticisms were meant to illustrate the general shortcomings of computer modeling, our results also argue against excessive skepticism about computational approaches in epidemiology.

Finally, the possibly questionable assumptions underlying our preliminary relative risk maps deserve explicit emphasis. The most important of these is that the relative abundances of species can be estimated using the relative probabilities of occurrence predicted by niche models. The relative risk maps based on the minimax model and, especially, the competitive exclusion model are less affected by this assumption than that based on the additive model as, in the former models, the relative probabilities of occurrence of the various vectors species were not aggregated. However, before the predictions of any of these models are fully accepted, the relationship between the predicted relative probabilities of occurrence of these species and their actual distributions must be empirically tested. The importance of such a test extends beyond epidemiology and will be relevant to all disciplines that use niche models (including conservation biology). It is essential to test this assumption for a wide variety of vector species before recommending its adoption. If the assumption holds, then absolute abundances, drawn from a portion of the area under investigation, can be used to calibrate the relative abundance predictions of the niche models, thus allowing for the determination of absolute abundances across the landscape for use in the risk models.

If the assumption does not hold, the situation is more difficult. In almost all epidemiological contexts, at the landscape or at larger scales, it is unlikely that the absolute abundances of vector species can be empirically measured for the entirety of the region of concern. It may be possible to obtain some measurements, and use traditional interpolation techniques such as splining and kriging to acquire estimates for other regions within the geographical boundaries of the measured cells. These can then be used along with the niche models to extrapolate across the landscape using statistical techniques such as regression to other areas that are predicted to be suitable (by the niche models) for the vector species. However, the accuracy of these traditional interpolation techniques in epidemiological contexts must likewise be empirically tested.

In addition to this primary assumption a number of secondary assumptions were made regarding the relationships between the vectors considered in this analysis. One such assumption was that the variable risk posed by P. vivax and P. falciparum could be ignored. As these two parasites differ both in the risk that they pose to humans [47] and in their geographical distribution [48], this assumption should be questioned. In this analysis the assumption was mandated by a lack of data on the varying abilities of most vectors to transmit these two parasites. A more complete analysis will require the explicit incorporation of the differential risk posed by these two parasites.

It was also assumed that the efficiency with which malaria parasites were transmitted, both to and from humans, is the same both for different vectors and in different regions. Since transmission efficiency varies between both species [49], [50] and regions [51], it should be explicitly included for a credible assessment of malaria risk. It was not included in this analysis due to a lack of data.

A similar caveat is needed with respect to the HBI values. Though each HBI value was treated as a species-specific parameter, these values are environmentally influenced. Measured HBI values for a species obtained at one region vary substantially from those measured at other regions, with the variation often depending primarily on the availability of human blood meals [52]. In addition, the data used to derive the HBI values used here were not obtained by random sampling; thus, there is reason to question whether these HBI values adequately reflect the varying anthrophilicity of the vectors. The spatial heterogeneity of HBI values also questions the use of a single HBI value for a vector. However, our sensitivity analysis shows that including such spatial heterogeneity would likely have had little impact upon the results. Relative risk would still have been shown to be concentrated in areas of high population density. Nevertheless, a more precise representation of malaria risk should consider such heterogeneity.

The questionable simplifying assumptions made in the construction of the preliminary relative risk maps presented in this paper show that to construct large-scale risk maps that go beyond the predicted relative probabilities of vector species presence will require much more data, and data of different types (including, for instance, abundance and transmission efficiency data) than what are now available in the literature. Furthermore, the use of more sophisticated transmission models than those on which Equations (1) and (2) are based will only require even more data. For disease risk analysis, research geared towards the acquisition of such data remain a high priority.

Acknowledgments

For discussion and help in developing this work we thank Trevon Fuller, Nelson Guda, Samraat Pawar, Steven Phillips, and Victor Sánchez-Cordero.

Author Contributions

Conceived and designed the experiments: SS AM NS. Performed the experiments: AM. Analyzed the data: SS AM. Wrote the paper: SS AM. Other: Wrote the first draft of the paper: SS.

References

  1. 1. Waller LA, Carlin BP, Xia H, Gelfand AE (1997) Hierarchical spatio-temporal mapping of disease rates. J Am Stat Assoc 92: 607–617.
  2. 2. MacNab YC, Dean CB (2001) Autoregressive spatial smoothing and temporal spline smoothing for mapping rates. Biometrics 57: 949–956.
  3. 3. Snow RW, Marsh K, Le Sueur D (1996) The need for maps of transmission intensity to guide malaria control in Africa. Parasitol Today 12: 455–457.
  4. 4. World Health Organization (2003) International travel and health: situations as on 1 January 2003. Geneva: World Health Organization.
  5. 5. Haydon DT, Randall DA, Matthews L, Knobel DL, Tallents LA, et al. (2006) Low-coverage vaccination strategies for the conservation of endangered species. Nature 443: 692–695.
  6. 6. Caldas de Castro M, Yamagat Y, Mtasiwa D, Tanner M, Utzinger J, et al. (2004) Integrated urban malaria control: a case study in Dar es Salaam, Tanzania. Am J Trop Med Hyg 71: 103–117.
  7. 7. Gemperli A, Sogoba N, Fondjo E, Mabaso M, Bagayoko M, et al. (2006) Mapping malaria transmission in West and Central Africa. Trop Med Int Health 11: 1032–1046.
  8. 8. Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs S, et al. (2004) A global index representing the stability of malaria transmission. Am J Trop Med Hyg 70: 486–498.
  9. 9. Peterson AT (2006) Ecological niche modeling and spatial patterns of disease transmission. Emerg Infect Dis 12: 1822–1826.
  10. 10. Soberón J, Peterson AT (2005) Interpretation of models of fundamental ecological niches and species' distributional areas. Biodiversity Informatics 2: 1–10.
  11. 11. Phillips SJ, Anderson RP, Schapire RE (2006) Maximum entropy modeling of species geographic distributions. Ecol Model 190: 231–259.
  12. 12. Illoldi-Rangel P, Sanchez-Cordero V, Peterson AT (2004) Predicting distributions of Mexican mammals using ecological niche modeling. J Mammal 85: 658–662.
  13. 13. Ortega-Huerta MA, Peterson AT (2004) Modelling spatial patterns of biodiversity for conservation prioritization in North-Eastern Mexico. Divers Distrib 10: 39–54.
  14. 14. Fuller T, Munguía M, Mayfield M, Sánchez-Cordero V, Sarkar S (2006) Incorporating connectivity into conservation planning: A multi-criteria case study from Central Mexico. Biol Conserv 133: 131–142.
  15. 15. Fuller T, Sánchez-Cordero V, Illoldi-Rangel P, Linaje M, Sarkar S (2007) The cost of postponing biodiversity conservation in Mexico. Biol Conserv 134: 593–600.
  16. 16. Pawar SS, Koo MS, Kelley C, Ahmed MF, Choudhury S, et al. (2007) Conservation assessment and prioritization of areas in Northeast India: priorities for amphibians and reptiles. Biol Conserv.. In press.
  17. 17. Elith J, Graham CH, Anderson RP, Dudik M, Ferrier S, et al. (2006) Novel methods improve predictions of species' distributions from occurrence data. Ecography 29: 129–151.
  18. 18. Stockwell DRB, Peters D (1999) The GARP modeling system: problems and solutions to automated spatial prediction. Int J Geogr Inf Sci 13: 143–158.
  19. 19. Breiman L (2001) Random Forests. Mach Learn 45: 5–32.
  20. 20. Lindsay SW, Parson L, Thomas CJ (1998) Mapping the ranges and relative abundance of the two principle African malaria vectors, Anopheles gambiae sensu stricto and An. arabiensis, using climate data. Proc Biol Sci 265: 847–854.
  21. 21. Kuhn KG, Campbell-Lendrum DH, Davies CR (2002) A continental risk map for malaria mosquito (Diptera: Culicidae) vectors in Europe.” J Med Entomol 39: 621–630.
  22. 22. Rogers DJ, Randolph SE, Snow RW, Hay SI (2002) Satellite imagery in the study and forecast of malaria. Nature 415: 710–715.
  23. 23. Levine RS, Peterson AT, Benedict MQ (2004) Geographic and ecologic distributions of the Anopheles gambiae complex predicted using a genetic algorithm. Am J Trop Med Hyg 70: 105–109.
  24. 24. Levine RS, Peterson AT, Benedit MQ (2004) Distribution of members of Anopheles quadrimaculatus Say s.l. (Diptera: Culicidae) and implications for their roles in malaria transmission in the United States. J Med Entomol 41: 607–613.
  25. 25. Kleinschmidt I, Omumbo J, Briet O, van de Giesen N, Sogoba N, et al. (2001) An empirical malaria distribution map for West Africa. Trop Med Int Health 6: 779–786.
  26. 26. Kleinschmidt I, Bagayoko M, Clarke GPY, Craig M, Le Sueur D (2000) A spatial statistical approach to malaria mapping. Int J Epidemiol 29: 355–361.
  27. 27. Diggle P, Moyeed R, Rowlingston B, Thomson M (2002) Childhood malaria in the Gambia: a case study in model-based geostatistics.” J Roy Stat Soc B 51: 493–506.
  28. 28. Smith DL, Dushoff J, McKenzie FE (2004) The risk of a mosquito-borne infection in a heterogeneous environment. PLoS Bio 2: e368.
  29. 29. Smith DL, McKenzie FE (2004) Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J 3: 13.
  30. 30. Black WC, Moore CG (2005) Population biology as a tool to study vector-borne diseases. In: Marquardt WH, Black WC, Freier JE, Hagedorn HH, Hemingway J, editors. Biology of disease vectors second edition. Burlington: Elsevier Academic Press.
  31. 31. Robert HJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 25: 1965–1978.
  32. 32. ESRI (2004) ArcGIS 9.1. Redlands: ESRI.
  33. 33. Hansen MC, Defries RS, Townshend JRG, Sohlberg R (2000) Global land cover classification at 1 km spatial resolution using a classification tree approach. Int J Remote Sens 21: 1331–1364.
  34. 34. Haworth J (1988) The global distribution of malaria and the present control effort. In: Wernsdorfer WH, McGregor I, editors. Malaria: principles and practice of malariology. Edinburgh: Churchill Livingstone. pp. 1379–1419.
  35. 35. Center for International Earth Science Information Network, Columbia University and Centro Internacional de Agricultura Tropical (2004) Gridded Population of the World (GPW), Version 3. Palisades: Columbia University. Available: http://beta.sedac.ciesin.columbia.edu/gpw.
  36. 36. Schneider P, Takken W, McCall PJ (2000) Interspecific competition between sibling species larvae of Anopheles arabiensis and An. gambiae. Med Vet Entomol 14: 165–170.
  37. 37. Antonio-Nkondjio C, Kerah CH, Simard F, Awono-Ambene P, Chouaibou , et al. (2006) Complexity of the malaria vectorial system in Cameroon: contribution of secondary vectors to malaria transmission. J Med Entomol 43: 1215–1221.
  38. 38. Cano J, Descalzo MA, Moreno M, Chen Z, Nzambo S, et al. (2006) Spatial variability in the density, distribution and vectorial capacity of anopheline species in a high transmission village (Equatorial Guinea). Malar j 5: 21.
  39. 39. Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW (2005) Urbanization, malaria transmission, and disease burden in Africa. Nat Rev Microbio 3: 81–90.
  40. 40. Guerra CA, Snow RW, Hay SI (2006) A global assessment of closed forest, deforestation, and malaria risk. Ann Trop Med Parasitol 100: 189–204.
  41. 41. Rogers DJ, Randolph SE (2000) The global spread of malaria in a future, warmer world. Science 289: 1763–1766.
  42. 42. United Nations (2002) World urbanization prospects: the 2001 revision. Data, tables, and highlights. New York: United Nations.
  43. 43. United Nations (2004) World urbanization prospects: the 2003 revision. Data, tables and highlights. New York: United Nations.
  44. 44. Gillies MT, Coetzee M (1987) A supplement to the Anophelinae of Africa South of the Sahara. Johannesburg: The South African Institute for Medical Research.
  45. 45. White GB (1989) Malaria. Geographical Distribution of Arthropod Borne Diseases and their Principal Vectors. Geneva: World Health Organization, WHO/VBC. pp. 7–22.
  46. 46. Coetzee M (2004) Distribution of the African malaria vectors of the Anopheles. Am J Trop Med Hyg 70: 103–104.
  47. 47. Collins FH, Paskewitz SM (1995) Malaria: current and future prospects for control. Annu Rev Entomol 40: 195–219.
  48. 48. Mendis K, Sina BJ, Marchesini P, Carter R (2001) The neglected burden of Plasmodium vivax malaria. Am J Trop Med Hyg 64: 97–106.
  49. 49. Gillies MT, Wilkes TJ (1965) A study of the age-composition of populations of Anopheles gambiae Giles and A. funestus Giles in North-Eastern Tanzania. Bull Entomol Res 56: 237–262.
  50. 50. Haji H, Smith T, Meuwissen JT, Sauerwein R, Charlwood JD (1996) Estimation of the infectious reservoir of Plasmodium falciparum in natural vector populations based on oocyst size. Trans R Soc Trop Med Hyg 90: 494–497.
  51. 51. Bonnet S, Gouagna LC, Paul RE, Safeukui I, Meunier JY, et al. (2003) Estimation of malaria transmission from humans to mosquitoes in two neighbouring villages in South Cameroon: evaluation and comparison of several indices. Trans R Soc Trop Med Hyg 97: 53–59.
  52. 52. Garrett-Jones C (1963) The human blood index of some anopheline mosquitoes, with reference to epidemiological assessment during malaria eradication. Geneva: World Health Organization, WHO/Mal/391.
  53. 53. Ameneshewa BB, Service AW (1996) Resting habits of Anopheles arabiensis in the Awash river valley of Ethiopia. Ann Trop Med Parasitol 90: 515–521.
  54. 54. Bryan JH, Petrarca V, Di Deco MA, Coluzzi M (1987) Adult behavior of members of the Anopheles gambiae complex in the Gambia with special reference to An. melas and its chromsomal variants. Parassitologia 29: 221–249.
  55. 55. Bryan JH, Di Deco MA, Petrarca V, Coluzzi M (1982) Inversion polymorphism and incipient speciation in Anopheles gambiae s.str. in The Gambia West Africa. Genetica 59: 167–176.
  56. 56. Cambournac FJ, Petrarca V, Coluzzi M (1982) Anopheles arabiensis in the Cape Verde archipelago. Parassitologia 24: 265–267.
  57. 57. Charlwood JD, Edoh D (1996) Polymerase chain reaction used to describe larval habitat use by Anopheles gambiae complex (Diptera: Culicidae) in the environs of Ifakara, Tanzania. J Med Entomol 33: 202–204.
  58. 58. Clarke D, Odialla H, Ouma J, Kenny V, MacCabe R, et al. (1996) A malariometric survey in Turkana District Kenya: chemosensitivity in vivo of Plasmodium falciparum infections and identity of the vector. Trans R Soc Trop Med Hyg 90: 302–304.
  59. 59. Coetzee M (1989) Comparative morphology and multivariate analysis for the discrimination of four members of the Anopheles gambiae group in Southern Africa. Mosq Sys 21: 100–116.
  60. 60. Coetzee M, Hunt RH, Baack L, Davidson G (1993) Distribution of mosquitoes belonging to the Anopheles gambiae complex, including malaria vectors, south of latitude 15°S. S Afr J Sci 89: 227–231.
  61. 61. Constantini C, Li SG, Della-Torre A, Sagnon N, Coluzzi M, et al. (1996) Density, survival and dispersal of Anopheles gambiae complex mosquitoes in a west African Savana village. Med Vet Entomol 10: 203–219.
  62. 62. Coosemans M, Petrarca V, Barutwanayo M, Coluzzi M (1989) Species of the Anopheles gambiae complex and chromosomal polymorphism in a rice-growing area of the Rusizi Valley (Republic of Burundi). Parassitologia 31: 113–122.
  63. 63. Faye O, Konate L, Mouchet J, Fontenille D, Sy N, et al. (1997) Indoor resting by outdoor biting females of Anopheles gambiae complex (Diptera: Culicidae) in the Sahel of northern Senegal. J Med Entomol 34: 285–289.
  64. 64. Faye O, Gaye O, Faye O, Diallo S (1994) Malaria transmission in villages, near or far of a mangrove area in Senegal. B Soc Pathol Exot 87: 157–163.
  65. 65. Fontenille D, Faye O, Konate L, Sy N, Collins FH (1993) Comparison of PCR and cytogenetic methods for the identification of mosquito species of the Anopheles gambiae complex in Senegal. Ann Parasitol Hum Comp 68: 239–240.
  66. 66. Fontenille D, Lepers JP, Coluzzi M, Campbell GH, Rakotoarivony I, et al. (1992) Malaria transmission and vector biology on Sainte Marie Island, Madagascar. J Med Entomol 29: 197–202.
  67. 67. Githeko AK, Adungo NI, Karanja DM, Hawley WA, Vulule JM, et al. (1996) Some observations on the biting behavior of Anopheles gambiae s.s Anopheles arabiensis and Anopheles funestus and their implication for malaria control. Exp Parasitol 82: 306–315.
  68. 68. Githeko AK, Brandling-Bennett AD, Beier M, Mbogo CM, Atieli FK, et al. (1993) Confirmation that Plasmodium falciparum has aperiodic infectivity to Anopheles gambiae. Med Vet Entomol 7: 373–376.
  69. 69. Green CA (1982) Cladistic analysis of mosquito chromosome data (Anopheles (Cellia) Myzomyia). J Hered 73: 3–11.
  70. 70. Hogg JC, Thomson MC, Hurd H (1996) Comparative fecundity and associated factors for two sibling species of the Anopheles gambiae complex occurring sympatrically in The Gambia. Med Vet Entomol 10: 385–391.
  71. 71. Ijumba JN, Mwangi RW, Beier JC (1990) Malaria transmission potential of Anopheles mosquitoes in the Mwea-Tebere irrigation scheme, Kenya. Med Vet Entomol 4: 425–432.
  72. 72. La Grange JJP (1995) Survey of anopheline mosquitoes (Diptera: Culicidae) in a malarious area of Swaziland. Afr Entomol 3: 217–219.
  73. 73. Mbogo CN, Baya NM, Ofulla AV, Githure JI, Snow RW (1996) The impact of permethrin-impregnated bednets on malaria vectors of the Kenyan coast. Med Vet Entomol 10: 251–259.
  74. 74. Mekuria Y, Petrarca V, Tesfamariam T (1982) Cytogenetic studies on the malaria vector mosquito Anopheles arabiensis Patton in the Awash Valley, Ethiopia. Parassitologia 24: 237–243.
  75. 75. Mnzava AE, Rwegoshora RT, Wilkes TJ, Curtis CF (1995) Anopheles arabiensis and An. gambiae chromosomal inversion polymorphism feeding and resting behaviour in relation to insecticide house-spraying in Tanzania. Med Vet Entomol 9: 316–324.
  76. 76. Mnzava AE, Mutinga MJ, Staak C (1994) Host blood meals and chromosomal inversion polymorphism in Anopheles arabiensis in the Baringo District of Kenya. J Am Mosq Control Assoc 10: 507–510.
  77. 77. Mnzava AE, Kilama WL, Kasigwa PF (1989) Application of a biochemical key to study transmission of malaria and Bancroftian filariasis in sibling species of the Anopheles gambiae complex in north-eastern Tanzania. Acta Trop 46: 323–333.
  78. 78. Mnzava AE, Kilama WL (1986) Observations on the distribution of the Anopheles gambiae complex in Tanzania. Acta Trop 43: 277–282.
  79. 79. Mosha FW, Petrarca V (1983) Ecological studies on Anopheles gambiae complex sibling species on the Kenya coast. Trans R Soc Trop Med Hyg 77: 344–345.
  80. 80. Mosha FW, Subra R (1982) Ecological studies on Anopheles gambiae complex sibling species in Kenya 1 preliminary observations on their geographical distribution and chromosomal polymorphic inversions. Document WHO/VBC/82867. Geneva: World Health Organization.
  81. 81. Mpofu SM (1985) Seasonal vector density and disease incidence patterns of malaria in an area of Zimbabwe. Trans R Soc Trop Med Hyg 79: 169–175.
  82. 82. Nigatu W, Curtis CF, Lulu M (1995) Test for association of DDT resistance with inversion polymorphism in Anopheles arabiensis from Ethiopia. J Am Mosq Control Assoc 11: 238–240.
  83. 83. Paskewitz SM, Ng K, Coetzee M, Hunt RH (1993) Evaluation of the polymerase chain reaction method for identifying members of the Anopheles gambiae (Diptera: Culicidae) complex in southern Africa. J Med Entomol 30: 953–957.
  84. 84. Petrarca V, Sabatinelli G, Youre YT, Di Deco MA (1998) Morphometric multivariate analysis of field samples of adult Anopheles arabiensis and An. gambiae s.s. (Diptera: Culicidae). J Med Entomol 35: 16–25.
  85. 85. Petrarca V, Beier JC, Onyango F, Koros J, Asiago C, et al. (1991) Species composition of the Anopheles gambiae complex (Diptera: Culicidae) at two sites in western Kenya. J Med Entomol 28: 307–313.
  86. 86. Petrarca V, Vercruysse J, Coluzzi M (1987) Observations on the Anopheles gambiae complex in the Senegal River Basin West Africa. Med Vet Entomol 1: 303–312.
  87. 87. Petrarca V, Petrangeli G, Rossi P, Sabatinelli G (1986) Chromosomal study of Anopheles gambiae and Anopheles arabiensis in Ouagadougou (Burkina Faso) and various neighboring villages. Parassitologia 28: 41–61.
  88. 88. Petrarca V, Carrara GC, Di Deco MA, Petrangeli G (1984) Observazioni citogenetiche e biometriche sui membri del complesso Anopheles gambiae in Mozambico. Parassitologia 26: 247–259.
  89. 89. Petrarca V, Carrara GC, Di Deco MA, Petrangeli G (1983) Il complesso Anopheles gambiae in Guinea Bissau. Parassitologia 25: 29–39.
  90. 90. Robert V, Carnevale P (1991) Influence of deltamethrin treatment of bed nets on malaria transmission in the Kou valley Burkina Faso. Bull World Health Organ 69: 735–740.
  91. 91. Robert V, Verhave JP, Carnevale P (1990) Plasmodium falciparum infection does not increase the precocious mortality rate of Anopheles gambiae. Trans R Soc Trop Med Hyg 84: 346–347.
  92. 92. Robert V, Petrarca V, Carnevale P, Ovazza L, Coluzzi M (1989) Analyse cytogenetique du complexe Anopheles gambiae dans la region de Bobo-Dioulasso (Burkina Faso). Ann Parasitol Hum Comp 64: 290–311.
  93. 93. Sharp BL, Le Sueur D, Wilken GB, Bredenkamp BL, Nqxongo S, et al. (1993) Assessment of the residual efficiency of lambda-cyhalothrin 2 A Comparison with DDT for the intradomiciliary control of Anopheles arabiensis in South Africa. J Am Mosq Control Assoc 9: 414–420.
  94. 94. Sharp BL, Le Sueur D (1991) Behavioral variation of Anopheles arabiensis (Diptera: Culicidae) populations in Natal, South Africa. Bull Entomol Res 81: 107–110.
  95. 95. Smits A, Roelants P, Van Bortel W, Coosemans M (1996) Enzyme polymorphisms in the Anopheles gambiae (Diptera: Culicidae) complex related to feeding and resting behavior in the Imbo Valley, Burundi. J Med Entomol 33: 545–553.
  96. 96. Taylor CE, Toure YT, Coluzzi M, Petrarca V (1993) Effective population size and persistence of Anopheles arabiensis during the dry season in West Africa. Med Vet Entomol 7: 351–357.
  97. 97. Thomson MC, Connor SJ, Quinones ML, Jawara M, Todd J, et al. (1995) Movement of Anopheles gambiae sl malaria vectors between villages in The Gambia. Med Vet Entomol 9: 413–419.
  98. 98. Toure YT, Traore SF, Sankare O, Sow MY, Coulibaly A, et al. (1996) Perennial transmission of malaria by the Anopheles gambiae complex in a north Sudan savanna area of Mali. Med Vet Entomol 10: 197–199.
  99. 99. Van Rensburg AJ, Hunt RH, Koekemoer LL, Coetzee M, Shiff CJ, et al. (1996) The polymerase chain reaction method as a tool for identifying members of the Anopheles gambiae complex (Diptera: Culicidae) in northeastern Tanzania. J Am Mosq Control Assoc 12: 271–274.
  100. 100. Walker TW, Robert LL, Copeland RA, Githeko AZ, Wirtz RA, et al. (1996) Field evaluation of arthropod repellents deet and a piperidine compound AI3-37220 against Anopheles funestus and Anopheles arabiensis in western Kenya. J Am Mosq Control Assoc 12: 172–176.
  101. 101. Adugna N, Petros B (1996) Determination of the human blood index of some anopheline mosquitos by using ELISA. Ethiop Med J 34: 1–10.
  102. 102. Aranda C, Aponte JJ, Saute F, Casimiro S, Pinto J, et al. (2005) Entomological characteristics of malaria transmission in Manhica a rural area in Southern Mozambique. J Med Entomol 42: 180–186.
  103. 103. Awolola TS, Okwa O, Hunt RH, Ogunrinade AF, Coetzee M (2002) Dynamics of the malaria-vector populations in coastal Lagos, south-western Nigeria. Ann Trop Med Parasitol 96: 75–82.
  104. 104. Carnahan J, Zheng L, Taylor CE, Toure YT, Norris DE, et al. (2002) Genetic differentiation of Anopheles gambiae s.s. populations in Mali West Africia using microsatellite loci. J Hered 93: 249–253.
  105. 105. Chen H, Githeko AK, Zhou G, Githure JI, Yan G (2006) New records of Anopheles arabiensis breeding on the Mount Kenya Highlands indicate indigenous malaria transmission. Malar J 5: 17.
  106. 106. Constantini C, Sagnon N, Della Torre A, Diallo M, Brady J, et al. (1998) Odor-mediated host preferences of West African mosquitoes with particular reference to malaria vectors. Am J Trop Med Hyg 58: 56–63.
  107. 107. Dekker T, Takken W (1998) Differential responses of mosquito sibling species Anopheles arabiensis and An. quadriannulatus to carbon dioxide a man or a calf. Med Vet Entomol 12: 136–140.
  108. 108. Fontenille D, Lochouarn L, Diatta M, Sokhna C, Dia I, et al. (1997) Four years' entomological study of the transmission of seasonal malaria in Senegal and the bionomics of Anopheles gambiae and A. arabiensis. Trans R Soc Trop Med Hyg 91: 647–652.
  109. 109. Habtewold T, Prior A, Torr SJ, Gibson G (2004) Could insecticide-treated cattle reduce Afrotropical malaria transmission? Effects of deltamethrin-treated Zebu on Anopheles arabiensis behavior and survival in Ethiopia. Med Vet Entomol 18: 408–417.
  110. 110. Hunt RH, Coetzee M, Fettene M (1998) The Anopheles gambiae complex: a new species from Ethiopia. Trans R Soc Trop Med Hyg 92: 231–235.
  111. 111. Ijumba JN, Mosha FW, Lindsay SW (2002) Malaria transmission risk variations derived from different agricultural practices in an irrigated area of northern Tanzania. Med Vet Entomol 16: 28–38.
  112. 112. Klinkenberg E, Takken W, Huibers F, Toure YT (2003) The phenology of malaria mosquitoes in irrigated rice fields in Mali. Acta Trop 85: 71–82.
  113. 113. Leong Pock Tsy J, Duchemin J, Marrama L, Rabarison P, Le Goff G, et al. (2003) Distribution of the species of the Anopheles gambiae complex and first evidence of Anopheles merus as a malaria vector in Madagascar. Malar J 2: 33.
  114. 114. Mendis C, Jacobsen JL, Gamage-Mendis A, Bule E, Dgedge M, et al. (2000) Anopheles arabiensis and An. funestus are equally important vectors of malaria in Matola coastal suburb of Maputo southern Mozambique. Med Vet Entomol 14: 171–180.
  115. 115. Muturi E, Shililu J, Jacob B, Gu W, Githure J, et al. (2006) Mosquito species diversity and abundance in relation to land use in a riceland agroecosystem in Mwea, Kenya. J Vector Ecol 31: 129–137.
  116. 116. Onyabe DJ, Conn JE (2001) The distribution of two major malaria vectors Anopheles gambiae and Anopheles arabiensis in Nigeria. Mem Inst Oswaldo Cruz 96: 1081–1084.
  117. 117. Petrarca V, Nugud AD, Elkarim Ahmed MA, Haridi AM, Di Deco MA, et al. (2000) Cytogenetics of the Anopheles gambiae complex in Sudan with special reference to An. arabiensis: relationships with East and West African populations. Med Vet Entomol 14: 149–164.
  118. 118. Robert V, Awono-Ambene HP, Thioulouse J (1998) Ecology of larval mosquitoes with special reference to Anopheles arabiensis (Diptera: Culicidae) in market-garden wells in urban Dakar, Senegal. J Med Entomol 35: 948–955.
  119. 119. Simard F, Lehmann T, Lemasson J, Diatta M, Fontenille D (2000) Persistence of Anopheles arabiensis during the severe dry season conditions in Senegal: an indirect approach using microsatellite loci. Insect Mol Biol 9: 467–479.
  120. 120. Taye A, Hadis M, Adugna N, Tilahun D, Wirtz RA (2006) Biting behavior and Plasmodium infection rates of Anopheles arabiensis from Sille, Ethiopia. Acta Trop 97: 50–54.
  121. 121. Toure YT, Petrarca V, Traore SF, Coulibaly A, Maiga HM, et al. (1998) The distribution and inversion polymporphism of chromosomally recognized taxa of the Anopheles gambiae complex in Mali, West Africa. Parasitologia 40: 477–511.
  122. 122. Tripet F, Dolo G, Lanzaro GC (2005) Multilevel analyses of genetic differentiation in Anopheles gambiae s.s. reveal patterns of gene flow important for malaria-fighting mosquito projects. Genetics 169: 313–324.
  123. 123. Tripet F, Toure YT, Dolo G, Lanzaro GC (2003) Frequency of multiple inseminations in field-collected Anopheles gambiae females revealed by DNA analysis of transferred sperm. Am J Trop Med Hyg 68: 1–5.
  124. 124. Wondji C, Simard F, Petrarca V, Etang J, Santolamazza F, et al. (2005) 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 42: 998–1005.
  125. 125. Woyessa A, Gebre-Michael T, Ali A (2004) An indigenous malaria transmission in the outskirts of Addis Ababa, Akaki Town and its environs. Ethiop J Health Dev 18: 2–7.
  126. 126. Yawson AE, McCall PJ, Wilson MD, Donnelly MJ (2004) Species abundance and insecticide resistance of Anopheles gambiae in selected areas of Ghana and Burkina Faso. Med Vet Entomol 18: 372–377.
  127. 127. Coene J (1993) Malaria in urban and rural Kinshasa: the entomological input. Med Vet Entomol 7: 127–137.
  128. 128. Diallo M, Ba Y, Sall AA, Diop OM, Ndione JA, et al. (2003) Amplification of the sylvatic cycle of dengue virus type 2 Senegal 1999–2000: entomological findings and epidemiological considerations. Emerg Infect Dis 9: 362–367.
  129. 129. Antonio-Nkondjio C, Awono-Ambene P, Toto JC, Meunier JY, Zebaze-Kemleu S, et al. (2002) High malaria transmission intensity in a village close to Yaounde the capital city of Cameroon. J Med Entomol 39: 350–355.
  130. 130. Coetzee M (1994) Anopheles crypticus, new species from South Africa is distinguished from Anopheles coustani (Diptera: Culicidae). Mosq Syst 26: 125–131.
  131. 131. Dia I, Diallo D, Duchemin J, Ba Y, Konate L, et al. (2005) Comparisons of human-landing catches and odor-baited entry traps for sampling malaria vectors in Senegal. J Med Entomol 42: 104–109.
  132. 132. Drakeley C, Schellenberg D, Kihonda J, Sousa CA, Arez AP, et al. (2003) An estimation of the entomological inoculation rate for Ifakara: a semi-urban area in a region of intense malaria transmission in Tanzania. Trop Med Int Health 8: 767–774.
  133. 133. Elissa N, Karch S, Bureau P, Ollomo B, Lawoko M, et al. (1999) Malaria transmission in a region of savanna-forest mosaic Haut-Ogooue Gabon. J Am Mosq Control Assoc 15: 15–23.
  134. 134. Gimnig JE, Ombok M, Kamau L, Hawley WA (2001) Characteristics of larval anopheline (Diptera: Culicidae) habitats in western Kenya. J Med Entomol 38: 282–288.
  135. 135. Logan TM, Linthicum KJ, Thande PC, Wagateh JN, Roberts CR (1991) Mosquito species collected from a marsh in western Kenya during the long rains. J Am Mosq Control Assoc 7: 395–399.
  136. 136. Mboera LEG, Takken W (1999) Odour-mediated host preference of Culex quinquefasciatus in Tanzania. Entomol Exp Appl 92: 83–88.
  137. 137. Sylla EHK, Kun JFJ, Kremsner PG (2000) Mosquito distribution and entomological inoculation rates in three malaria-endemic areas in Gabon. Trans R Soc Trop Med Hyg 94: 652–656.
  138. 138. Antonio-Nkondjio C, Simard F, Awono-Ambene P, Ngassam P, Toto J, et al. (2005) Malaria vectors and urbanization in the equatorial forest region of south Cameroon. Trans R Soc Trop Med Hyg 99: 347–354.
  139. 139. Traore-Lamizana M, Fontenille D, Diallo M, Ba Y, Zeller HG, et al. (2001) Arbovirus surveillance from 1990 to 1995 in the Barkedji area (Ferlo) of Senegal, a possible natural focus of rift valley fever virus. J Med Entomol 38: 480–492.
  140. 140. Ayala D, Le Goff G, Robert V, de Jong P, Takken W (2006) Population structure of the malaria vector Anopheles funestus (Dipetera: Culicidae) in Madagascar and Comoros. Acta Trop 97: 292–300.
  141. 141. Cohuet A, Dia I, Simard F, Raymond M, Fontenille D (2004) Population structure of the malaria vector Anopheles funestus in Senegal based on microsatellite and cytogenetic data. Insect Mol Biol 13: 251–258.
  142. 142. Cohuet A, Simard F, Wondji CS, Antonio-Nkondjio C, Awono-Ambene P, et al. (2004) High malaria transmission intensity due to Anopheles funestus (Diptera: Culicidae) in a village of savannah-forest transition area in Cameroon. J Med Entomol 41: 901–905.
  143. 143. Cohuet A, Simard F, Toto J, Kengne P, Coetzee M, et al. (2003) Species identification within the Anopheles funestus group of malaria vectors in Cameroon and evidence for a new species. Am J Trop Med Hyg 69: 200–205.
  144. 144. Cuamba N, Choi KS, Townson H (2006) Malaria vectors in Angola: distribution of species and molecular forms of the Anopheles gambiae complex their pyrethroid insecticide knockdown resistance (kdr) status and Plasmodium falciparum sporozoite rates. Malar J 5: 2.
  145. 145. Dossou-yovo J, Doannio JMC, Riviere F, Chauvancy G (1995) Malaria in Côte d'Ivoire wet savannah region: the entomological input. Trop Med Parasitol 46: 263–269.
  146. 146. Duchemin JB, Leong Pock Tsy JM, Rabarison P, Roux J, Coluzzi M, et al. (2001) Zoophily of Anopheles arabiensis and An. gambiae in Madagascar demonstrated by odourbaited entry traps. Med Vet Entomol 15: 50–57.
  147. 147. Fillinger U, Sonye G, Killeen GF, Knols BGJ, Becker N (2004) The practical importance of permanent and semipermanent habitats for controlling aquatic stages of Anopheles gambiae sensu lato mosquitoes: operational observations from a rural town in western Kenya. Trop Med Int Health 9: 1274–1289.
  148. 148. Kamau L, Munyekenye GO, Koekemoer LL, Hunt RH, Coetzee M (2003) A survey of Anopheles funestus (Diptera: Culicidae) group of mosquitoes from 10 sites in Kenya with special emphasis on population genetic structure based on chromosomal inversion karyotypes. J Med Entomol 40: 664–671.
  149. 149. Kamau L, Hunt R, Coetzee M (2002) Analysis of the population structure of Anopheles funestus (Diptera: Culicidae) from western and coastal Kenya using paracentric chromosomal inversion frequencies. J Med Entomol 39: 78–83.
  150. 150. Mukabayire O, Boccolini D, Lochouarn L, Fontenille D, Besansky NJ (1999) Mitochondrial and ribosomal internal transcribed spacer (ITS2) diversity of the African malaria vector Anopheles funestus. Mol Ecol 8: 289–297.
  151. 151. Quakyi IA, Leke RGF, Befidi-Mengue R, Tsafack M, Bomba-Nkolo D, et al. (2000) The epidemiology of Plasmodium falciparum malaria in two Cameroon villages: Simbok and Etoa. Am J Trop Med Hyg 63: 222–230.
  152. 152. Zhou G, Minakawa N, Githeko A, Yan G (2004) Spatial distribution patterns of malaria vectors and sample size determination in spatially heterogeneous environments: a case study in the West Kenyan highland. J Med Entomol 41: 1001–1009.
  153. 153. Akogbeto M, Di Deco MA, Romano R, Coluzzi M (1988) Analisi comparativa del ruolo vettore per la malaria de Anopheles gambiae e Anopheles melas nella zona lagunare costiera del Benin. Parassitologia 30: 7–8.
  154. 154. Appawu MA, Baffoe-Wilmot A, Afari EA, Nkrumah FK, Petrarca V (1994) Species composition and inversion polymorphism of the Anopheles gambiae complex in some sites of Ghana west Africa. Acta Trop 56: 15–23.
  155. 155. Bockarie MJ, Service MW, Barnish G, Maude GH, Greenwood BM (1994) Malaria in a rural area of Sierra Leone III Vector ecology and disease transmission. Ann Trop Med Parasitol 88: 251–262.
  156. 156. Coluzzi M, Petrarca V, Di Deco MA (1985) Chromosomal inversion intergradation and incipient speciation in Anopheles gambiae. Boll Zool 52: 45–63.
  157. 157. Jaenson TGT, Gomes MJ, Barreto dos Santos RC, Petrarca V, Fortini , et al. (1994) Control of endophagic Anopheles mosquitoes and human malaria in Guinea Bissau West Africa by permethrin-treated bednets. Trans R Soc Trop Med Hyg 88: 620–624.
  158. 158. Matola YG, White GB, Magayuka SA (1987) The changed pattern of malaria endemicity and transmission at Amani in the eastern Usambara mountains noth-eastern Tanzania. J Trop Med Hyg 90: 127–134.
  159. 159. Petrarca V, Sabatinelli G, Di Deco MA, Papakay M (1990) The Anopheles gambiae complex in the Federal Islamic Republic of Comoros (Indian Ocean): some cytogenetic and biometric data. Parassitologia 32: 371–380.
  160. 160. Robert V, Le Goff G, Toto JC, Mulder L, Fondjo E, et al. (1993) Anthrophilic mosquitoes and malaria transmission at Edea Cameroon. Trop Med Parasitol 44: 14–18.
  161. 161. Smits A, Coosemans M, Van Bortel W, Barutwanayo M, Delacollette C (1995) Readjustment of the malaria vector control strategy in the Rusizi Valley, Burundi. Bull Entomol Res 85: 541–548.
  162. 162. Toure YT, Petrarca V, Coluzzi M (1983) Nuove entita del complesso Anopheles gambiae in Mali. Parassitologia 25: 367–370.
  163. 163. Vulule JM, Beach RF, Atieli FK, Roberts JM, Mount DL, et al. (1994) Reduced susceptibility of Anopheles gambiae to permethrin associated with the use of permethrin-impregnated bednets and curtains in Kenya. Med Vet Entomol 8: 71–75.
  164. 164. Afrane YA, Lawson BW, Githeko AK, Yan G (2005) Effects of microclimatic changes caused by land use and land cover on duration of gonotrophic cycles of Anopheles gambiae (Diptera: Culicidae) in Western Kenya Highlands. J Med Entomol 42: 974–980.
  165. 165. Awolola TS, Oyewole IO, Amajoh CN, Idowu ET, Ajayi MB, et al. (2005) Distribution of the molecular forms of Anopheles gambiae and pyrethroid knock down resistance gene in Nigeria. Acta Trop 95: 204–209.
  166. 166. Cano J, Angel Desalzo M, Moreno M, Chen Z, Nzambo S, et al. (2006) Spatial variability in the density distribution and vectorial capacity of anopheline species in a high transmission village (Equatorial Guinea). Malar J 5: 21.
  167. 167. Charalambous M, Townson H, Harbach RE, Mukwaya LG (1999) Electrophoretic and DNA identification of Anopheles bwambae and A. gambiae (Diptera: Culicidae) in western Uganda. Bull Entomol Res 89: 111–117.
  168. 168. Charlwood JD, Pinto J, Sousa CA, Ferreira C, Petrarca V, et al. (2003) ‘A mate or a meal’–Pre-gravid behaviour of female Anopheles gambiae from the islands of Sao Tome and Principe West Africa. Malar J 2: 9.
  169. 169. della Toree A, Fanello C, Akogbeto M, Dossou-yovo J, Favia G, et al. (2001) Molecular evidence of incipient speciation within Anopheles gambiae s.s.. in West Africa. Insect Mol Biol 10: 9–18.
  170. 170. Lehman T, Licht M, Elissa N, Maega BTA, Chimumbwa JM, et al. (2003) Population structure of Anopheles gambiae in Africa. J Hered 94: 133–147.
  171. 171. Munga S, Minakawa N, Zhou G, Barrack OJ, Githeko AK, et al. (2005) Oviposition site preference and egg hatchability of Anopheles gambiae: effects of land cover types. J Med Entomol 42: 993–997.
  172. 172. Pinto J, Lynd A, Elissa N, Donnelly MJ, Costa C, et al. (2006) Co-occurrence of East and West African kdr mutations suggests high levels of resistance to pyrethroid insecticides in Anopheles gambiae from Libreville, Gabon. Med Vet Entomol 20: 27–32.
  173. 173. Lambert D (1983) A population genetical study of the African mosquito Anopheles marshallii (Theobald). Evolution 37: 484–495.
  174. 174. Sharp BL (1983) Anopheles merus (Donitz) its biting cycle in relation to environmental parameters. J Entomol Soc South Afr 46: 367–374.
  175. 175. Sharp BL, Quicke FC, Jansen EJ (1984) Aspects of the behavior of 5 anopheline species in the endemic malaria area of Natal South Africa. J Entomol Soc South Afr 47: 251–258.
  176. 176. Gad AM, el Said S, Hassan AN, Shoukry A (1987) The distribution and ecology of the mosquitoes in the Red Sea governorate, Egypt. J Egypt Soc Parasitol 17: 207–221.
  177. 177. Awone-Ambene HP, Kengne P, Simard F, Antonio-Nkondjio C, Fontenille D (2004) Description and bionomics of Anopheles (Cellia) ovengensis (Diptera: Culicidae) a new malaria vector species of the Anopheles nili group from South Cameroon. J Med Entomology 41: 561–568.
  178. 178. Kenge P, Awono-Ambene P, Antonio-Nkondjio C, Simard F, Fontenille D (2003) Molecular identification of the Anopheles nili group of African malaria vectors. Med Vet Entomol 17: 67–74.
  179. 179. el Said S, Beier JC, Kenawy MA, Morsy ZS, Merdan AI (1986) Anopheles population dynamics in two malaria endemic villages in Faiyum Governorate, Egypt. J Am Mosq Control Assoc 2: 158–163.
  180. 180. Fettene M, Koekemoer LL, Hunt RH, Coetzee M (2002) PCR assay for identification of Anopheles quadriannulatus species B from Ethiopia and other sibling species of the Anopheles gambiae complex. Med Vet Entomol 16: 214–217.
  181. 181. Bogh C, Clarke SE, Pinder M, Sanyang F, Lindsay SW (2001) Effect of passive zooprophylaxis on malaria transmission in the Gambia. J Med Entomol 38: 822–828.
  182. 182. Charlwood JD, Qassim M, Elnsur EI, Donnelly M, Petrarca V, et al. (2001) The impact of indoor spraying with malathion on malaria in refugee camps in eastern Sudan. Acta Trop 80: 1–8.
  183. 183. Fettene M, Hunt RH, Coetzee M, Tessema F (2004) Behaviour of Anopheles arabiensis and An. quadriannulatus sp. B mosquitoes and malaria transmission in southwestern Ethiopia. Aft Entomol 12: 83–87.
  184. 184. Garrett-Jones C, Boreham PFL, Pant CP (1980) Feeding habits of Anophelines (Diptera: Culicidae) in 1971-1978 with reference to the human blood index: a review. Bull Entomol Res 70: 165–185.
  185. 185. Githeko AK, Service MW, Mbogo CM, Atieli FK, Juma FO (1994) Origin of blood meals in indoor and outdoor resting malaria vectors in western Kenya. Acta Trop 58: 307–316.
  186. 186. Habtewold T, Walker AR, Curtis CF, Osir EO, Thapa N (2001) The feeding behavior and Plasmodium infection of Anopheles mosquitoes in southern Ethiopia in relation to use of insecticide-treated livestock for malaria control. Trans R Soc Trop Med Hyg 95: 584–586.
  187. 187. Hadis M, Lulu M, Makonnen Y, Asfaw T (1997) Host choice by indoor-resting Anopheles arabiensis in Ethiopia. Trans R Soc Trop Med Hyg 91: 376–378.
  188. 188. Highton RB, Bryan JH, Boreham PFL, Chandler JA (1979) Studies on the sibling species Anopheles gambiae Giles and Anopheles arabiensis Patton (Diptera: Culicidae) in the Ksumu area. Kenya Bull Entomol Res 69: 43–53.
  189. 189. Joshi GP, Service MW, Pradhan GD (1975) 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 69: 91–104.
  190. 190. Lemasson JJ, Fontenille D, Lochouarn L, Dia I, Simard F, et al. (1997) Comparison of behavior and vector efficiency of Anopheles gambiae and An. arabiensis (Diptera: Culicidae) in Barkedji a Sahelian area of Senegal. J Med Entomol 34: 396–403.
  191. 191. Mwangangi JM, Mbogo CM, Nzovu JG, Githure JI, Yan G, et al. (2003) Blood-meal analysis for anopheline mosquitoes sampled along the Kenyan coast. J Am Mosq Control 19: 371–375.
  192. 192. Service MW Joshi GP, Pradhan GD (1978) A Survey of Anopheles gambiae (Species A) and Anopheles arabiensis (species B) of the Anopheles gambiae Giles complex in the Kisumu area of Kenya following insecticidal spraying with OMS-43 (Fenitrothion). Ann Trop Med Parasitol 72: 377–386.
  193. 193. White GB, Rosen P (1973) Comparative studies on sibling species of the Anopheles gambiae Giles complex (Diptera: Culicidae). II. ecology of species A and B in Savanna around Kaduna, Nigeria, during transition from wet to dry season. Bull Entomol Res 62: 613–625.
  194. 194. White GB, Magayuka SA, Boreham PFL (1972) Comparative studies on sibling species of the Anopheles gambiae Giles complex (Diptera: Culicidae): bionomics and vectorial activity of species A and species B at Segara, Tanzania. Bull Entomol Res 62: 295–317.
  195. 195. Beier JC, Perkins PV, Wirtz RA, Koros J, Diggs D, et al. (1988) Bloodmeal identification by direct enzyme linked immunosorbent assay (ELISA) tested on Anopheles (Diptera: Culicidae) in Kenya. J Med Entomol 25: 9–16.
  196. 196. Bogh C, Pedersen EM, Mukoko DA, Ouma JH (1998) Permethrin-impregnated bednet effects on resting and feeding behavior of lymphatic filariasis vector mosquitoes in Kenya. Med Vet Entomol 12: 52–59.
  197. 197. Chandler JA, Boreham PFL, Highton RB, Hill MN (1975) A study of the host selection patterns of the mosquitoes of the Kisumu area of Kenya. Trans R Soc Trop Med Hyg 69: 415–425.
  198. 198. Mbogo CN, Kabiru EW, Muiruri SK, Nzovu JM, Ouma JH, et al. (1993) Bloodfeeding behavior of Anopheles gambiae s.l. and Anopheles funestus in Kilifi District, Kenya. J Am Mosq Control Assoc 9: 225–227.
  199. 199. White GB (1971) Blood feeding habits of mosquitoes in the South Pare District of Tanzania ten years after cessation of a dieldrin residual spraying campaign. East Afr Med J 48: 120–134.
  200. 200. Sousa CA, Pinto J, Almeida AP, do Rosario VE, Charlwood JD (2001) Dogs as a favored host choice of Anopheles gambiae sensu stricto (Diptera: Culicidae) of São Tomé West Africa. J Med Entomol 38: 122–125.