The goal of malaria elimination necessitates an improved understanding of any fine-scale geographic variations in transmission risk so that complementary vector control tools can be integrated into current vector control programmes as supplementary measures that are spatially targeted to maximize impact upon residual transmission. This study examines the distribution of host-seeking malaria vectors at households within two villages in rural Tanzania.
Host-seeking mosquitoes were sampled from 72 randomly selected households in two villages on a monthly basis throughout 2008 using CDC light-traps placed beside occupied nets. Spatial autocorrelation in the dataset was examined using the Moran’s I statistic and the location of any clusters was identified using the Getis-Ord Gi* statistic. Statistical associations between the household characteristics and clusters of mosquitoes were assessed using a generalized linear model for each species.
For both Anopheles gambiae sensu lato and Anopheles funestus, the density of host-seeking females was spatially autocorrelated, or clustered. For both species, houses with low densities were clustered in the semi-urban village centre while houses with high densities were clustered in the periphery of the villages. Clusters of houses with low or high densities of An. gambiae s.l. were influenced by the number of residents in nearby houses. The occurrence of high-density clusters of An. gambiae s.l. was associated with lower elevations while An. funestus was also associated with higher elevations. Distance from the village centre was also positively correlated with the number of household occupants and having houses constructed with open eaves.
The results of the current study highlight that complementary vector control tools could be most effectively targeted to the periphery of villages where the households potentially have a higher hazard (mosquito densities) and vulnerability (open eaves and larger households) to malaria infection.
The frontline vector tools deployed in the battle against malaria transmission are long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS)
[1, 2]. These tools are highly effective and their use has led to a significant reduction of transmission in many parts of Africa, including places that were historically holoendemic
[3–6]. In response to such success, the international community has now prioritized regional and national malaria elimination, with a long-term goal of malaria eradication
. However, vector control that solely targets insecticides to the inside of houses is unlikely to be sufficient to achieve elimination
. Thus there is a need for complementary vector control tools to target a range of alternative stages in the mosquito cycle, such as the larval stage, mating or sugar feeding. Such complementary tools will target specific ecosystems and understanding the fine-scale geographic variations in Anopheles mosquitoes and transmission risk will enable tools to be developed and effectively integrated into current vector control programmes. This need for detailed geographic research has been highlighted by recent literature demonstrating that malaria transmission is highly heterogeneous across the landscape
Changes in malaria transmission risk can be measured by the entomological inoculation rate (EIR)
[13, 14], which is the product of the anopheline biting rate and the proportion of infectious females (sporozoite rate). With regard to the anopheline biting rate, household-level characteristics have been demonstrated to influence biting rates, such as the number of occupants, screened windows, closed eaves or ceilings
[15–17]. Further, some level of geographic clustering has been observed where houses closest to breeding sites tend to experience higher adult biting rates
[10, 12, 18–21] which is supported with mathematical modelling
. However, in peri-urban or rural situations where houses are spread over large distances (many kms) and often inter-dispersed with larval habitats, the heterogeneity of anopheline biting rates may be influenced by household characteristics in addition to distance from the nearest breeding site. Understanding the spatial clustering of anopheline biting rates at the household level is important because households are the focal point of many predictive malaria models and, most importantly, are perhaps the easiest of targets for delivery of vector control measures through either vertical or horizontal delivery strategies.
The current study therefore examines the geographic relationships of anopheline host-seeking patterns at a household and village level in East Africa. The null hypothesis tested was that the adult biting rate would be randomly distributed across the households within villages. For this analysis, household-level characteristics (elevation, the number of occupants, closed eaves, the presence of bed nets, and distance from the village centre) were taken into account.
The study was conducted in the neighbouring villages of Namawala and Idete, located in the Kilombero Valley (8.1°S and 36.6°E), south-eastern Tanzania (Figure
1). These communities experience hyperendemic malaria transmission; primarily vectored by large populations of Anopheles gambiae sensu lato. In this area, this species complex has been historically represented by the two morphologically identical sibling species: Anopheles gambiae sensu stricto and Anopheles arabiensis[10, 24] but it should be noted that the proportional contribution of the former has been dramatically reduced following scale up of LLINs in the area
[24, 25]. A third, locally important vector species is Anopheles funestus sensu stricto. The ecosystem is dominated by a low-lying river valley, 150 km long and up to 40 km wide, which is interdispersed with villages and rice farms. Annual rains (December to May) create large quantities of ephemeral aquatic habitat suitable for An. gambiae s.l. oviposition and larval development. Both villages have semi-urban town centres with many people also residing in the rural farming regions on the outskirts (Idete = 1,229 and Namawala = 767 households). Idete village has been constructed in a slightly more elevated area compared with Namawala (Figure
1). There are subtle cultural differences between the villages, with more pastoral farmers residing in Namawala; for example, during 2008, there were 306 head of cattle in Idete and 6,667 in Namawala
Mosquito population sampling
Seventy-two households were randomly selected for mosquito sampling in both villages using census information from the Ifakara Health Institute (IHI) Demographic Surveillance System (Figure
1). For each household, the location and elevation was recorded using a handheld GPS unit (eTrex, Vista, Garmin Inc, USA). The distance of each household from the village centre was calculated using ArcGIS 10 software (ESRI, Redlands, CA, USA). Additionally, the physical structure of the eave openings was recorded directly by observation. The use of bed nets (either LLINs or insecticide-treated nets [ITNs]) and the number of household occupants was recorded during a one-time survey of all household heads.
Each house was visited once a month (six houses/day, four days/week and three weeks/month) over a period of 12 months (January to December 2008). Mosquitoes were collected inside each house using one CDC light trap for 12 hours (7 pm to 7 am)
. The light trap, fitted with an incandescent bulb, was placed 1–1.5 m above the floor and close to the feet of an LLIN occupant
. The LLIN used was provided to each participating household by the study team (Olyset, A to Z Textile Mills Ltd, Tanzania). Although permethrin-treated bed nets exhibit modest excito-repellency, they have surprising little effect on the relative efficiency of light traps when compared with untreated bed nets
[28, 29]. Traps were inspected each morning and all mosquitoes were morphologically identified to species complex or group and classified by sex
. The sibling species identity of the An. gambiae s.l. complex was identified using PCR
. For both An. gambiae s.l. and An. funestus, any specimens that contained sporozoites in the salivary glands were identified using ELISA
. Owing to the large number of female mosquitoes caught per trapping effort (up to approx 1,500), separate random subsamples, each averaging approximately 10% of the total in each trap, were used for molecular analysis. In cases where the catch was less than 10 females, molecular analysis was conducted for all individuals. Prior to molecular analysis, mosquitoes were stored at −20°C in micro centrifuge tubes containing a small amount of silica drying agent. The age structure of the An. gambiae s.l. population was estimated using parity dissections on a subset of the samples caught
Spatial and statistical analysis
The spatial and statistical analyses were conducted for An. gambiae s.l. and An. funestus. Regarding An. gambiae s.l., the ratio of An. arabiensis to An. gambiae s.s. was examined with a binomial GLMM with a categorical explanatory variable for week and adjusted for multiple comparisons with Dunnett’s test.
The spatial patterns in the dataset were analysed using the geographical information systems software ArcGIS 10 with Spatial Analysis and Statistical Tools (ESRI, Redlands, CA, USA). To assess any spatial patterns over time, the 12-month sampling period was broken into four time periods: January to February, March to April, May to June and July to December
. The time periods were selected to reflect the temporality of the system where the mosquito densities undergo large fluctuations during the wet season (January to June; thus three divisions) with less variation during the dry season (July to December; thus one division). For each month, one light trap sample was collected from each household; to enable spatial analysis with each household being a point of interest, the total number of mosquitoes collected from all trapping efforts (within the specified time frame) were summed.
Initially, the spatial patterns of An. gambiae s.l. and An. funestus densities within the two villages were mapped. Next any spatial autocorrelation patterns, i e, clustered, dispersed, random, were analysed using the Moran’s I statistic
. Sequentially, the localities of clustered households with high or low anopheline densities were identified using the Getis-Ord Gi* statistic
. Statistically significant (at a level of 0.05) clusters of households with high densities of anophelines were identified with Z scores >1.96, or vice versa, clustered households with low densities of anophelines were identified with Z scores < −1.96. For all spatial analyses, the spatial relationship among houses was conceptualized using the inverse distance, which is most appropriate for continuous point datasets because closer houses have larger influences on the computations for each target house than houses that are further away. Also, the distance between neighbouring features was calculated using the Euclidean distance and were run separately for each village.
Next, any association between the location of anopheline clustering and household characteristics was investigated. In Namawala and Idete, there was a clear socio-economic gradient from the semi-urban centres towards the rural village outskirts, and thus there were potential correlations between the household characteristics with distance from the village centre. The correlation of household characteristics, which were ordinal data (the number of occupants and the number of bed nets per person), was investigated using Pearson’s correlation coefficient. The presence of eaves was a binary factor (open or closed), and the correlation of this parameter with distance was investigated with a generalized linear model (GLM) with a binomial distribution and a logit link function. Sequentially, statistical associations between the household characteristics and clusters of mosquitoes were assessed using a GLM for each species. To account for spatial autocorrelations, the GLM was run with the Getis-Ord Gi* Z Score as the dependent factor and with a normal distribution. The independent factors incorporated in the model were: elevation (m), the presence of eaves, the number of occupants, the number of bed nets per person and distance from the village centre (m). Interaction terms were included for any correlated independent variables. This analysis was conducted using R statistical software (ver.2.14.2).
Ethical approval for the study was obtained from the IHI Institutional Review Board (IHRDC/IRB/No. A-32) and the Medical Research Coordination Committee of the National Institute for Medical Research (NIMR/HQ/R.8a/Vol. IX/764) in Tanzania. When the study commenced, permission was obtained from each household owner who, after consenting, signed an informed consent form stating their willingness to participate in the study.
During the 12-month, mosquito-sampling period, 1,648 light-trap nights of sampling were conducted. A total of 97,437 female mosquitoes were caught. Of these mosquitoes, 30.9% (n = 30,111) were An. gambiae s.l., comprising 85.8% Anopheles arabiensis and 14.2% An. gambiae s.s. (n = 2,924 PCR amplifications). The remaining mosquitoes were 2.0% An. funestus (n = 1,950), 62.0% Culex spp (n = 60,442), 2.4% Mansonia spp (n = 2,302) and 2.7% other species including Aedes and Coquillettidia spp (n = 2,605). The ratio of An. arabiensis to An. gambiae s.s. remained constant throughout the study (binomial GLMM, p >0.05, Figure
2). Therefore the detailed spatial analysis was conducted for the species complex overall. The average number of An. gambiae s.l. per light-trap night was 18.3 ± 2.3 and for An. funestus was 1.2 ± 0.1. The parity rate was only calculated for An. gambiae s.l. and the portion of the population that were parous was 43% (397/916).
Extremely high densities of An. gambiae s.l. occurred during the wet season, especially during March to April; whereas the density of An. funestus was relatively consistent throughout the year (Figure
). Mapping the spatial patterns of An. gambiae s.l. and An. funestus densities indicated that households with high or low densities of anophelines tended to be closer together (Figure
3) and this was confirmed with the Moran’s I spatial analysis. For both An. gambiae s.l. and An. funestus, the existence of spatial autocorrelation, or clustering, of host-seeking densities was evidenced with positive z scores (Table
1). Specifically for An. gambiae s.l., clustering was evident in Namawala during all time periods, except May to June, and was significant overall. For An. gambiae s.l. in Idete, clustering was only significant during July to December. With regard to An. funestus, clustering was significant in both villages during all time periods, except March to April, and was significant overall. The age structure (parity) of An. gambiae s.l. and houses with sporozoite positive An. gambiae s.l. and An. funestus specimens were randomly distributed across the landscape (Table
2). Thus, the biting rate was the only component of the EIR which demonstrated any spatial autocorrelation, this is supported by Magbity and Lines
. The localities of high and low clusters of anopheline densities were identified with the Z scores computed by Getis-Ord Gi* (Figure
4). For An. gambiae s.l., there were nine households with high densities and 46 with low densities. For An. funestus, there were seven households with high densities and 96 with low densities.
Moran’sIspatial autocorrelation indices forAnopheles gambiaes.l. andAnopheles funestusdensities over time in Namawala and Idete villages
Village and period
An. gambiae s.l.
Moran’sIspatial autocorrelation indices for the sporozite positivity and parity ofAnopheles gambiaes.l. andAnopheles funestusspecimens collected in Namawala and Idete villages
Sporozoite positive –An. gambiaes.l.
Sporozoite positive –An. funestus
Parity –An. gambiaes.l.
a No sporozoite positive An. funestus were detected in Namawala village.
The influence of household characteristics on the localities of anopheline clustering was sequentially investigated (Table
3). The independent household characteristics that were correlated with distance from the village centre were: number of occupants (t = 2.662; p = 0.0087), elevation (t = −4.535; p = <0.0001) and eaves (z = 2.883; p = 0.0039) (see Figure
5). The mean distance of houses with open eaves to the centre of the village was 1,881 m, while the mean distance of houses with closed eaves was 629 m. The total number of bed nets owned by each household was not correlated with distance from the village centre (t = 0.440; p = 0.6604).
The mean characteristics of households that were in clusters of high or low densities ofAnopheles gambiaes.l. andAnopheles funestus
Eaves – Closed
Eaves – Closed
aClusters of households with high densities of anophelines were identified with Z scores >1.96, or vice versa, clustered households with low densities of anophelines were identified with Z scores < −1.96.
bThe distance of each household from the village centre was calculated using ArcGIS software (ESRI, Redlands, CA, USA).
The multivariate GLM examining the association between household characteristics and clustering of mosquitoes found that distance from the village centre significantly influenced the occurrence of anopheline clusters (Table
6). For both An. gambiae s.l. and An. funestus, houses with low densities were clustered in the semi-urban centre of Idete and Namawala, and there was almost no variability in the location of these low-density clusters over time. Conversely, houses with high densities of An. gambiae s.l. and An. funestus were clustered in the rural outskirts of both villages. There was some seasonal variability in the locality of the high-density clusters. Broadly, households located towards the periphery of each village had a higher chance of being located with a cluster of households that had higher densities of anophelines.
Multivariate statistical comparison of household factors with clustering ofAnopheles gambiaes.l. andAnopheles funestusdensities
Eaves – Closed
Distance from centre (m)
Distance from centre (m) × Elevation (m)
Distance from centre (m) × Eaves – Closed
Distance from centre (m) × Occupants
Eaves – Closed
Distance from centre (m)
Distance from centre (m) × Elevation (m)
Distance from centre (m) × Eaves – Closed
Distance from centre (m) × Occupants
The association of anopheline clusters (Getis-Ord Gi* Z Score) with household factors was compared with a generalized linear model (GLM) with a normal distribution.
With regard to the remaining household characteristics, the occurrence of An. gambiae s.l. clusters was also associated with elevation and the number of occupants; both of these factors interacted with distance from the village centre (Table
4). Elevation was negatively associated with clusters. Notably households with high densities of An. gambiae s.l. occurred in the south of Namawala, being on the flood-plain and in close association with larval habitats (Figure
5). The number of occupants in a household was positively associated with high An. gambiae s.l. densities. It is important to note that households with higher numbers of occupants were generally clustered outside of the semi-urban centres (Figure
5). For An. funestus, elevation also significantly influenced the location of clusters and the influence of this factor interacted with distance from the village centre. The influence of this interaction can be seen in Figure
6: the high-density houses have diverged from the general pattern for distance from the village centre.
Interestingly, the number of bed nets/person was not associated with the spatial clustering of either An. gambiae s.l. or An. funestus. This occurred because the bed nets were fairly evenly distributed across the landscape with minimal or no evidence of clustering (Figure
5). This represents the equity of the national bed net distribution system as it operated in the Kilombero Valley
 and nationally
. The nets represented various distribution schemes and 46.8% of nets in use were either LLINs or had been treated within 1 year (for more detail see
Previously, models have demonstrated that the proportion of mosquitoes that are infectious is influenced by the age structure of the population. Specifically, the highest proportion of positive mosquitoes was modelled to be in localities with older mosquitoes, usually being the middle of villages and away from breeding sites
. This study did not find any statistical evidence for spatial clustering of the age structure of mosquitoes or the localities where sporozoite-infected, host-seeking female mosquitoes were caught. Consistent with other studies
, of the two factors which comprise EIR (biting rate and sporozoite rate), only the biting rate expressed strong spatial heterogeneity. This study, along with the published literature
[9–12] indicates that host-seeking mosquito densities are clustered and thus the risk of malaria infection at a household level is relatively more similar for close-by neighbours and not necessarily similar to the risk experienced by households further away but still situated in the same village. Such heterogeneous clustering of anopheline densities within villages
[23, 39, 40] has usually been associated with the proximity to larval sites at this fine-scale level
The number of household occupants positively influenced the densities of An. gambiae s.l. This is supported by previous research demonstrating that the number of household occupants does influence mosquito densities at an individual household level
. Thus the combined body odour of the many occupants may have attracted more mosquitoes to the area
[41, 42] and caused the host-seeking mosquitoes to aggregate
. In the current study, houses with more occupants, and also open eaves, tended to occur towards the periphery of the villages where higher densities of anophelines were also observed. Both the number of occupants and house construction are proxy indicators of socioeconomic status
[43, 44] and both have previously been associated with increased malaria risk
[15, 18, 21, 40, 45, 46]. Thus, this study demonstrates that as households are located towards the periphery of each village there is a gradient of increasing hazard (mosquito densities) and vulnerability (open eaves and larger households) which coincide in parallel from the village centre moving outwards.
For An. gambiae s.l., there was some mild variation in the locality of the high density clusters over time. This is most likely related to the ephemeral nature of larval habitats that changed drastically throughout the year. In the Kilombero Valley, the occurrence of An. gambiae s.l. larval sites is closely related to elevation: the flood plain occurs at lower elevations where the gradient of the land is gentler and water is able to pool and form larval sites, being one of the most influential factors associated with the location of high-density clusters
. For An. funestus, the localities of high-density clusters were localized with little spatial-temporal variability, as such it is plausible that the locality of the high clusters was related to the availability of larval sites. Especially since the larvae of An. funestus utilize large permanent or semi-permanent vegetated aquatic habitats such as stream pools and swamps
. For both An. gambiae s.l. and An. funestus, the locality of low-density clusters was consistent over time and may have been influenced by the closing of household eaves in the low-density clusters. The elevation of the study houses varied by less than 100 m (lowest house: 244 m; highest house: 334 m), thus it is unlikely that the temperature or relative humidity experienced by individual houses varied greatly; even though these factors are known to influence the distribution of vectors at broader scales
Apart from the LLINs that were provided to the study houses to comply with ethical protocols, LLINs and other bed nets were fairly evenly distributed across the landscape and therefore did not contribute to within village spatial clustering of anopheline densities. Nonetheless, the distribution of LLINs has impacted on the community-wide transmission in the villages by reducing the host-seeking density, sporozoite rate and the EIR of the anopheline populations overall
. Particularly in this study area, the scale-up in LLINs has seen a dramatic reduction in the density of vectors which feed indoors during the middle of the night, such as An. gambiae s.s.
[24, 25]. Consequently, any residual transmission is maintained by mosquitoes which express exophagic feeding behaviours, such as An. arabiensis, which can feed outdoors at dusk or dawn
[25, 50, 51]. Thus exists a need for complementary vector control interventions to help further suppress malaria transmission by targeting stages of the mosquito cycle, other than indoor biting, such as the larval stages, mating or sugar feeding
. The results of the current study highlight that such complementary vector control tools could be most effectively targeted to the periphery of villages where the households potentially have a higher hazard (mosquito densities) and vulnerability (open eaves and larger households) to malaria infection.
We thank the staff at the Ifakara Health Institute especially Japhet Kihonda, Peter Pazia, Daniel Lugiko, Nuru Nchimbi, Matrida Kideule and Cesilia Mchopa for technical assistance in conducting entomological surveillance and Deodatus Maliti, Bernard Liseki, Fransis Allen for assistance with molecular analysis. We thank the residents of Idete and Namawala village for their patience and cooperation in the entomology surveys. These entomological surveys were funded by the Addessium Foundation (Reeuwijk, The Netherlands). The manuscript has been approved by the National Institute for Medical Research in Tanzania (NIMR/HQ/P.12 VOL.XII/34).
Ifakara Health Institute, Environmental Sciences Thematic Group
Liverpool School of Tropical Medicine, Vector Group
Faculty of Medicine, Health and Molecular Sciences, James Cook University
Laboratory of Entomology, Wageningen University and Research Centre
Liverpool School of Tropical Medicine, Centre for Neglected Tropical Diseases
WHO: Global malaria control and elimination: report of a technical review. Geneva: World Health Organization; 2008.
WHO: World malaria report 2011. Geneva: World Health Organization; 2011.
Bhattarai A, Ali AS, Kachur SP, Mårtensson A, Abbas AK, Khatib R, Al-mafazy A-w, Ramsan M, Rotllant G, Gerstenmaier JF, Molteni F, Abdulla S, Montgomery SM, Kaneko A, Björkman A: Impact of artemisinin-based combination therapy and insecticide-treated nets on malaria burden in Zanzibar.PLoS Med 2007, 4:e309.PubMedView Article
Ceesay SJ, Casals-Pascual C, Nwakanma DC, Walther M, Gomez-Escobar N, Fulford AJC, Takem EN, Nogaro S, Bojang KA, Corrah T, Jaye MC, Taal MA, Sonko AAJ, Conway DJ: Continued decline of malaria in The Gambia with implications for elimination.PLoS One 2010, 5:e12242.PubMedView Article
D’Acremont V, Lengeler C, Genton B: Reduction in the proportion of fevers associated with Plasmodium falciparum parasitaemia in Africa: a systematic review.Malar J 2010, 9:240.PubMedView Article
O’Meara WP, Mangeni JN, Steketee R, Greenwood B: Changes in the burden of malaria in sub-Saharan Africa.Lancet Infect Dis 2010, 10:545–555.PubMedView Article
Feachem R, Sabot O: A new global malaria eradication strategy.Lancet 2008, 371:1633–1635.PubMedView Article
The malERA Consultative Group on Vector Control: A research agenda for malaria eradication: vector control.PLoS Med 2011, 8:e1000401.View Article
Smith DL, Dushoff J, McKenzie FE: The risk of a mosquito-borne infection in a heterogeneous environment.PLoS Biol 2004, 2:1957–1964.View Article
Smith T, Charlwood JD, Takken W, Tanner M, Spiegelhalter DJ: Mapping the densities of malaria vectors within a single village.Acta Trop 1995, 59:1–18.PubMedView Article
Trape JF, Lefebvre-Zante E, Legros F, Ndiaye G, Bouganali H, Druilhe P, Salem G: Vector density gradients and the epidemiology of urban malaria in Dakar, Senegal.Am J Trop Med Hyg 1992, 47:181–189.PubMed
Lindsay SW, Armstrong Schellenberg JRM, Zeiler HA, Daly RJ, Salum FM, Wilkins HA: Exposure of Gambian children to Anopheles gambiae malaria vectors in an irrigated rice production area.Med Vet Entomol 1995, 9:50–58.PubMedView Article
Shaukat A, Breman J, McKenzie FE: Using the entomological inoculation rate to assess the impact of vector control on malaria parasite transmission and elimination.Malar J 2010, 9:122.PubMedView Article
Kelly-Hope L, McKenzie FE: The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa.Malar J 2009, 8:19.PubMedView Article
Kirby MJ, Green C, Milligan PM, Sismanidis C, Jasseh M, Conway DJ, Lindsay SW: Risk factors for house-entry by malaria vectors in a rural town and satellite villages in The Gambia.Malar J 2008, 7:2.PubMedView Article
Ogoma SB, Lweitoijera DW, Ngonyani H, Furer B, Russell TL, Mukabana WR, Killeen GF, Moore SJ: Screening mosquito house entry points as a potential method for integrated control of endophagic filariasis, arbovirus and malaria vectors.PLoS Negl Trop Dis 2010, 4:e773.PubMedView Article
Kirby MJ, Ameh D, Bottomley C, Green C, Jawara M, Milligan PM, Snell PC, Conway DJ, Lindsay SW: Effect of two different house screening interventions on exposure to malaria vectors and on anaemia in children in The Gambia: a randomised control trial.Lancet 2009, 6736:61078.
Gunawardena DM, Wickremasinghe AR, Muthuwatta L, Weerasingha S, Rajakaruna J, Senanayaka T, Kotta PK, Attanayake N, Carter R, Mendis KN: Malaria risk factors in an endemic region of Sri Lanka, and the impact and cost implications of risk factor-based interventions.Am J Trop Med Hyg 1998, 58:533–542.PubMed
Clarke SE, Bøgh C, Brown RC, Walraven GEL, Thomas CJ, Lindsay SW: Risk of malaria attacks in Gambian children is greater away from malaria vector breeding sites.Trans R Soc Trop Med Hyg 2002, 96:499–506.PubMedView Article
Cohen J, Ernst K, Lindblade K, Vulule J, John C, Wilson M: Local topographic wetness indices predict household malaria risk better than land-use and land-cover in the western Kenya highlands.Malar J 2010, 9:328.PubMedView Article
Zhou G, Munga S, Minakawa N, Githeko AK, Yan G: Spatial relationship between adult malaria vector abundance and environmental factors in western Kenya highlands.Am J Trop Med Hyg 2007, 77:29–35.PubMed
Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, Bousema T, Drakeley CJ, Ferguson NM, Basáñez M-G, Ghani AC: Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies.PLoS Med 2010, 7:e1000324.PubMedView Article
Bousema T, Drakeley C, Gesase S, Hashim R, Magesa S, Mosha F, Otieno S, Carneiro I, Cox J, Msuya E, Kleinschmidt I, Maxwell C, Greenwood B, Riley E, Sauerwein R, Chandramohan D, Gosling R: Identification of hot spots of malaria transmission for targeted malaria control.J Infect Dis 2010, 201:1764–1774.PubMedView Article
Russell TL, Lwetiojera DW, Maliti D, Chipwaza B, Kihonda J, Charlwood JD, Smith TA, Lengeler C, Mwanyangala MA, Nathan R, Knols BGJ, Takken W, Killeen GF: Impact of promoting longer-lasting insecticide treatment of bednets upon malaria transmission in a rural Tanzanian setting with pre-existing high coverage of untreated nets.Malar J 2010, 9:187.PubMedView Article
Russell TL, Govella NJ, Azizi S, Drakeley CJ, Kachur SP, Killeen GF: Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania.Malar J 2011, 10:80.PubMedView Article
Lines JD, Curtis CF, Wilkes TJ, Njunwa KJ: Monitoring human-biting mosquitoes (Diptera: Culicidae) in Tanzania with light-traps hung beside mosquito nets.Bull Entomol Res 1991, 81:77–84.View Article
Mboera LEG, Kihonda J, Braks MA, Knols BGJ: Short report: Influence of centers for disease control light trap position, relative to a human-baited bed net, on catches of Anopheles gambiae and Culex quinquefasciatus in Tanzania.Am J Trop Med Hyg 1998, 59:595–596.PubMed
Magbity EB, Lines JD, Marbiah MT, David K, Peterson E: How reliable are light traps in estimating biting rates of adult Anopheles gambiae s.l. (Diptera: Culicidae) in the presence of treated bed nets?Bull Entomol Res 2002, 92:71–76.PubMedView Article
Killeen GF, Tami A, Kihonda J, Okumu FO, Kotas ME, Grundmann H, Kasigudi N, Ngonyani H, Mayagaya V, Nathan R, Abdulla S, Charlwood JD, Smith TA, Lengeler C: Cost-sharing strategies combining targeted public subsidies with private-sector delivery achieve high bednet coverage and reduced malaria transmission in Kilombero Valley, southern Tanzania.BMC Infect Dis 2007, 7:121.PubMedView Article
Gillies MT, Coetzee M: A supplement to the Anophelinae of Africa south of the Sahara (Afrotropical region). Johannesburg: South African Institute for Medical Research; 1987.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction.Am J Trop Med Hyg 1993, 49:520–529.PubMed
Burkot TR, Williams JL, Schneider I: Identification of Plasmodium falciparum infected mosquitoes by a double antibody enzyme-linked immunosorbent assay.Am J Trop Dis Hyg 1984, 33:783–788.
Gillies MT: A modified technique for the age-grading of populations of Anopheles gambiae.Ann Trop Med Parasitol 1958, 52:261–273.PubMed
Moran PAP: The interpretation of statistical maps.J R Stat Soc Series B Stat Methodol 1948, 10:243–251.
Getis A, Ord JK: The anlaysis of spatial association by use of distance statistics.Geogr Anal 1992, 24:189–206.View Article
Armstrong Schellenberg JRM, Abdulla S, Minja H, Nathan R, Mukasa O, Marchant T, Mponda H, Kikumbih N, Lyimo E, Manchester T, Tanner M, Lengeler C: KINET: a social marketing programme of treated nets and net treatment for malaria control in Tanzania, with evaluation of child health and long-term survival.Trans R Soc Trop Med Hyg 1999, 93:225–231.View Article
Bonner K, Mwita A, McElroy P, Omari S, Mzava A, Lengeler C, Kaspar N, Nathan R, Ngegba J, Mtung’e R, Brown N: Design, implementation and evaluation of a national campaign to distribute nine million free LLINs to children under five years of age in Tanzania.Malar J 2011, 10:73.PubMedView Article
Magbity EB, Lines JD: Spatial and temporal distribution of Anopheles gambiae s.l. (Diptera: Culicidae) in two Tanzanian villages: implications for designing mosquito sampling routines.Bull Entomol Res 2002, 92:483–488.PubMedView Article
Brooker S, Clarke S, Njagi JK, Polack S, Mugo B, Estambale B, Muchiri E, Magnussen P, Cox J: Spatial clustering of malaria and associated risk factors during an epidemic in a highland area of western Kenya.Trop Med Int Health 2004, 9:757–766.PubMedView Article
Gamage-Mendis AC, Carter R, Mendis C, De Zoysa AP, Herath PR, Mendis KN: Clustering of malaria infections within an endemic population: risk of malaria associated with the type of housing construction.Am J Trop Med Hyg 1991, 45:77–85.PubMed
Kelly DW: Why are some people bitten more than others?Trends Parasitol 2001, 17:578–581.PubMedView Article
Somi MF, Butler JR, Vahid F, Njau JD, Kachur SP, Abdulla S: Use of proxy measures in estimating socioeconomic inequalities in malaria prevalence.Trop Med Int Health 2008, 13:354–364.PubMedView Article
Somi MF, Butler JRG, Vahid F, Njau J, Kachur SP, Abdulla S: Is There evidence for dual causation between malaria and socioeconomic status? Findings from rural Tanzania.Am J Trop Med Hyg 2007, 77:1020–1027.PubMed
Ghebreyesus TA, Haile M, Witten KH, Getachew A, Yohannes M, Lindsay SW, Byass P: Household risk factors for malaria among children in the Ethiopian highlands.Trans R Soc Trop Med Hyg 2000, 94:17–21.PubMedView Article
Haddow AJ: The mosquito fauna and climate of native huts at Kisumu, Kenya.Bull Entomol Res 1942, 33:91–142.View Article
Kelly-Hope L, Hemingway J, McKenzie FE: Environmental factors associated with the malaria vectors Anopheles gambiae and Anopheles funestus in Kenya.Malar J 2009, 8:268.PubMedView Article
Sinka M, Bangs M, Manguin S, Coetzee M, Mbogo C, Hemingway J, Patil A, Temperley W, Gething P, Kabaria C, Okara R, Van Boeckel T, Godfray HC, Harbach R, Hay S: The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic precis.Parasit Vectors 2010, 3:117.PubMedView Article
De Souza D, Kelly-Hope L, Lawson B, Wilson M, Boakye D: Environmental factors associated with the distribution of Anopheles gambiae s.s in Ghana; an important vector of lymphatic filariasis and malaria.PLoS One 2010, 5:e9927.PubMedView Article
Bayoh MN, Mathias D, Odiere M, Mutuku F, Kamau L, Gimnig J, Vulule J, Hawley W, Hamel M, Walker E: Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya.Malar J 2010, 9:62.PubMedView Article
Mutuku F, King C, Mungai P, Mbogo C, Mwangangi J, Muchiri E, Walker E, Kitron U: Impact of insecticide-treated bed nets on malaria transmission indices on the south coast of Kenya.Malar J 2011, 10:356.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.