Skip to main content
  • Research article
  • Open access
  • Published:

Do hotspots fuel malaria transmission: a village-scale spatio-temporal analysis of a 2-year cohort study in The Gambia

Abstract

Background

Despite the biological plausibility of hotspots fueling malaria transmission, the evidence to support this concept has been mixed. If transmission spreads from high burden to low burden households in a consistent manner, then this could have important implications for control and elimination program development.

Methods

Data from a longitudinal cohort in The Gambia was analyzed. All consenting individuals residing in 12 villages across the country were sampled monthly from June (dry season) to December 2013 (wet season), in April 2014 (mid dry season), and monthly from June to December 2014. A study nurse stationed within each village recorded passively detected malaria episodes between visits. Plasmodium falciparum infections were determined by polymerase chain reaction and analyzed using a geostatistical model.

Results

Household-level observed monthly incidence ranged from 0 to 0.50 infection per person (interquartile range = 0.02–0.10) across the sampling months, and high burden households exist across all study villages. There was limited evidence of a spatio-temporal pattern at the monthly timescale irrespective of transmission intensity. Within-household transmission was the most plausible hypothesis examined to explain the observed heterogeneity in infections.

Conclusions

Within-village malaria transmission patterns are concentrated in a small proportion of high burden households, but patterns are stochastic regardless of endemicity. Our findings support the notion of transmission occurring at the household and village scales but not the use of a targeted approach to interrupt spreading of infections from high to low burden areas within villages in this setting.

Peer Review reports

Background

Within populations, heterogeneity in exposure to malaria has been widely documented; it is generally estimated that 20% of the population experience 80% of the disease burden [1,2,3]. The skewed distribution of exposure has been observed at every spatial scale, in different transmission landscapes, and is expected to be more pronounced when transmission is low [4]. Several studies have documented both spatial and spatio-temporal high burden areas of malaria, typically referred to as hotspots but here defined as clusters, and have fueled the notion of spatially targeting interventions for control and elimination [5,6,7].

The consistent presence of spatial clusters of high malaria burden within populations contributed to the hypothesis that there may be hotspots, or certain households, or subsets of households within foci (spatially discrete areas with sustained transmission) that fuel transmission [8]. The number and size of clusters within foci and the delineation of a foci itself will likely depend on the specific setting. For example, on the coast of Kenya, multiple clusters were identified per foci [2], whereas a single cluster was observed in a highland setting [6]. If such clusters are in fact hotspots, meaning they are drivers of malaria transmission, and they could be easily identified and targeted with interventions, then resources could be used more effectively and their impact on transmission intensity may be greater than that of a uniform approach [8, 9]. For a hotspot-driven approach at the sub-village level to be viable, it is critical to determine whether the observed heterogeneity at the village scale is a feature of malaria transmission and supports the notion of “hotspots” fueling transmission or whether it follows a more stochastic pattern [10].

The notion of hotspots as intrinsic drivers of malaria transmission being an inherent part of the transmission landscape is plausible with risk being driven by macroscale and microscale characteristics [11, 12]. For example, the observed seasonality in transmission is associated with climate, specifically the rainfall patterns and temperature [13, 14]. Similarly, at the local scale, malaria risk is known to be associated with microepidemiological variation in risk factors, including greater odds of infection in those residing in proximity to mosquito breeding sites (e.g., ponds or forests) or living with other infected individuals [15,16,17]. The observed spatial heterogeneity in infected individuals also has implications for quantifying and understanding transmission intensity [18]. As described as part of the hotspot model, the high burden households within an endemic area may amplify transmission by acting as a constant parasite reservoir, or equally they could absorb infectious bites, attenuating observed transmission events [19, 20]. If these households or groups of households are driving transmission within foci, then hotspot-targeted interventions would be justified [8, 21].

Although biologically plausible, the evidence to support the concept of hotspots, here considered as a single household or group of high burden households within foci, fueling transmission has been mixed. For example, a recent trial targeting serologically defined hotspots of exposure failed to observe any sustained reduction in transmission outside of the targeted area [22]. Transmission in the study area may have been too high for well-defined hotspots, hotspot boundaries may not have been effectively defined, or hotspots may not have contributed to maintaining transmission in this setting [23]. Despite the limited evidence to support the use of hotspot-targeted approaches, several malaria elimination programs have engaged in hotspot-inspired strategies [3, 21, 24].

In this study, we conducted a spatio-temporal analysis on a full population cohort distributed in six pairs of villages across The Gambia. The aim of this research was to establish if predicted risk of malaria transmission intensity exhibits a consistent pattern, meaning the risk of malaria moving from a high burden household or a group of households to neighboring households, over time. If the expected pattern exists, we aimed to identify at what transmission intensity this dynamic becomes apparent. In case of limited evidence to support the hotspot pattern, some potential drivers of any observed heterogeneity were explored.

Methods

Malaria transmission in The Gambia is highly seasonal and occurs during and soon after the rainy season, typically between August and December. Epidemiological data from the study cohort has been recently described [25]. Briefly, monthly blood samples were collected during the 2013 and 2014 malaria transmission seasons (June–December) from all people residing in every household in the study villages (Fig. 1). An additional blood sample was collected during the dry season, in April 2014. Village pairs were approximately 1–3 km apart and were considered as discrete spatial units. Populations ranged between 100 and 700 individuals per village, and all residents were included in the study. All households were geo-located. The number of households per village ranged from 13 to 69, and the distance between households within a village ranged from 0.4 to a maximum of 986.8 m (Table 1). Furthermore, one round of mass drug administration (MDA) with dihydroartemisinin-piperaquine was carried out in June 2014.

Fig. 1
figure 1

Map of The Gambia showing the location of the 12 study villages. The study villages are represented as circles and labeled A–H and J–M. The circles are colored according to the overall observed malaria prevalence

Table 1 Key characteristics of study villages including demographics and the observed malaria burden

Finger prick blood samples were collected on filter paper for identification of Plasmodium falciparum infections using polymerase chain reaction (PCR). All febrile individuals (auxiliary temperature ≥ 37.5 °C or history of fever in the last 24 h) were screened for malaria by rapid diagnostic test (RDT), and if positive they were treated with artemether-lumefantrine according to national guidelines. A study nurse was stationed within each village and recorded all malaria episodes between monthly visits, including administering an RDT and collecting a blood sample on filter paper.

Malaria parasites are transmitted to humans via the bite of an infected Anopheles mosquito and can be directly measured using the entomological inoculation rate (EIR) [26, 27]. The P. falciparum parasite rate (PfPR) is a known correlate to EIR; it provides a measure of transmission intensity and is a more operationally feasible metric to generate [28]. Using PCR infection as the dependent variable as a proxy for transmission intensity, geostatistical analysis was conducted using the PrevMap package in R (v3.3.2) to determine the predicted malaria prevalence per household per month within each village accounting for spatial autocorrelation as well as temporal trends [29]. A Bayesian geostatistical probit model was used to predict the spatial variation in malaria parasite prevalence within each village. More details on the model specification are provided in Additional file 1. Because the cohort was a full population sample, no interpolation at unsampled locations was required. Predicted prevalence per household was estimated using the median of the posterior distribution, and maps of the combined and monthly predicted prevalence were generated.

Models were adjusted for sample date, distance to road, distance to river, and mean monthly rainfall. The distance to river and road variables were determined by extracting the relevant features from pan-sharpened Landsat 8 imagery and using the gDistance function in the rgeos package [30] to estimate the straight-line distance in kilometers. Monthly rainfall was obtained from weather stations located in each of the six study regions across the country.

The observed overdispersed distribution of infection counts has been used to support the notion of malaria hotspots [1]. However, it is possible that the skewed distribution is due to measurement bias in how infections are defined. For example, PCR-detected infections were not treated in this study (until becoming symptomatic and detectable by RDT) and could represent an infection from a single infectious bite or repeated inoculations within the same individual until treatment is sought. For example, by considering each time point where a PCR infection is detected as unique would lead to counting a single infection detected at 5 sequential time points as 5 unique infections instead of 1, thereby driving the observed overdispersion. To demonstrate the degree of potential measurement bias as an alternative explanation to the skewed patterns of case counts, two different approaches for counting were employed. The first was to consider each time point when a PCR positive result was recorded as a unique infection irrespective of whether there was a confirmed treatment in between sampling. The second approach considered any infections detected at sequential time points as the same infection unless the individual had been treated for malaria as part of the study. Any negative sample between two PCR positive samples in a non-treated individual was assumed to be a false negative and considered as a single infection. Any subsequent infection detected after a known treatment event (e.g., symptomatic and RDT positive, or participated in the MDA) was considered as a new infection.

Results

In total, 41,548 monthly observations were available from 360 households across 14 sampling time points. The size of households ranged from a single person to 78 individuals, and the residents had a similar age distribution between villages (Table 1). The aggregated infection prevalence across the study period ranged from 2.6 to 18.3% across the 12 villages (Fig. 1). During the 2-year study period, 2877 samples were positive for malaria infection, with substantial heterogeneity between villages. The lowest transmission village recorded 34 infections in 10 households, whereas the village with the highest transmission had 845 infections in 42 households (Table 1). Across all time points, 12.5% (45/360) households did not record a single infection, while the number of households without any infection varied from 0 in village L to 12 in village F. Household-level observed monthly incidence ranged from 0 to 0.50 infection per person (interquartile range (IQR) = 0.02–0.10) across the sampling months.

The overall number of observed infections per individual (Fig. 2a) and per household (Fig. 2b) exhibit the expected overdispersion pattern, illustrating the considerable heterogeneity in malaria exposure experienced by this population. Results of the geostatistical model exhibited 100 m as the range of spatial autocorrelation, suggesting that village pairs were discrete transmission units. However, the geostatistical model failed to provide evidence of a pronounced spatial pattern within villages at either low or high transmission intensities (Fig. 3; see Additional file 1 for model output). Across all villages, only a single village (Fig. 3, village F) showed a pattern of high burden households grouping together. When the predicted household-level prevalence is plotted over time, there is no evidence that infection dynamics around high burden households exhibit a regular pattern around neighboring households at the monthly time step; the patterns appear stochastic (Additional file 2). Furthermore, as a group of high burden households was only evident in a single village with moderate transmission levels, the presence of hotspots within villages does not appear to be associated with transmission intensity. The spatial patterns were similar irrespective of whether an infection was symptomatic or asymptomatic (Additional file 3) [25].

Fig. 2
figure 2

Frequency distributions of malaria infections in the study population. Frequency of number of observed PCR positive infections a per individual and b per household

Fig. 3
figure 3

Overall predicted PCR prevalence per household (circles), per village (panels a-h, j-m, corresponding to the village code) according to the spatio-temporal model. The size and color of the circles are scaled according to prevalence. The black dots identify those households with zero malaria infections recorded during the study

Monthly predicted PCR prevalence per household (circles) for all study villages (panels) according to the spatio-temporal model. The size and color of the circles are scaled according to prevalence. Each household is identified by a black dot. Households with a predicted PCR prevalence between 0 and 1% are identified in greyscale. (MP4 285 kb)

As a consistent spatio-temporal dynamic of malaria around high burden households was not observed, the next step was to explore alternative explanations for the overdispersion pattern of malaria burden in the study population. The first explanation examined was measurement bias in how infections were defined. If we consider the most conservative definition and assume only new infections as those after a recorded treatment event, the distribution becomes less skewed, with fewer households experiencing multiple malaria episodes (Fig. 4a). As expected, the differences between methods for counting infections are more pronounced in high transmission settings (Fig. 4b, village M) compared to low transmission settings (Fig. 4b, village A). Although neither method of counting infections is expected to fully capture the number of “true” infections experienced in the population, the heterogeneity in malaria burden was still present despite the most extreme definition of counting infections being applied.

Fig. 4
figure 4

Density plot for the number of infections per compound according to the two definitions tested. Distributions according to definitions are provided for a the combined data and b an example of a low (A) and high (M) transmission village. The red curves show the distribution if each time point with an infection is counted as new. The blue curves show the distribution of unique infections assuming an infection is only counted as new if there is evidence of treatment at a prior time point

The second explanation for the observed heterogeneity in malaria that we explored was to consider the household as the relevant spatial unit of transmission. Patterns of infections appearing within households suggested that three scenarios are evident: there are cases when several individuals are infected within the same month, there are cases of infections appearing the month after another individual within the household becomes infected, and there are cases of stochastic introductions (Fig. 5). All patterns were observed in households in both the low (Fig. 5; village A) and high (Fig. 5; village M) transmission settings. However, parasite genetic data is required to confirm this hypothesis. See Additional file 4 for heat maps showing transmission dynamics within all study households.

Fig. 5
figure 5

Heat maps showing within-household transmission dynamics. Heat maps showing within-household transmission dynamics in a low transmission village (village A) and a high transmission village (village M). Each grid represents a household with each individual residing within the household shown in the rows. Each column within each grid represents a sampling month starting in June 2013 through December 2014. The color of each grid cell represents their infection and/or treatment status at that time point. Infection status is defined by those who are PCR positive with treatment being administered when there was a symptomatic infection confirmed by RDT in the field or the mass drug administration (MDA) administered between transmission seasons (June 2014)

Discussion

Heterogeneity in malaria burden is an inherent aspect of transmission, rooted in complex interactions between environmental, vector, and individual characteristics [9, 21, 31]. However, evidence on the importance of the observed heterogeneity within a village in maintaining or fueling transmission, consistent with the concept of hotspots, is required to support the use of such a strategy as part of control or elimination programs. In this study, we explored spatio-temporal trends of malaria transmission intensity to see if it shifted from high burden households to the surrounding area. Though high burden households within villages exist, they were not consistently the same, and the risk of malaria was not observed to spread from high to low burden households at the monthly timescale. Together, these findings suggest that the relevant operational unit for targeting transmission in this setting is the household or the entire village, depending on the program goals and interventions being employed.

As heterogeneity in malaria infections was observed in the data, we next explored non-spatial factors that could be driving the pattern. In this study, participants were only treated if they had a symptomatic, RDT positive infection or participated in the MDA. Therefore, we hypothesized that each observed infection is unlikely to represent a unique infection event, and the overdispersion in burden may be partly driven by measurement bias. Assuming that new infections are only those identified after documented antimalarial treatment decreased but did not eliminate the observed heterogeneity. The “extreme” assumptions we used, namely that all detected infections are new ones or that new infections are only those occurring after treatment, are unlikely to represent the true number of infection events, as individuals may have cleared them spontaneously, received treatment outside of the study, or experienced superinfections [32,33,34]. Being able to account for superinfections and identify the role of these individuals in fueling onward transmission would help refine methods for counting new or incident infections and determining which infections matter for maintaining transmission intensity [35]. The true incidence likely falls somewhere in between the two estimates used, but measurement bias is unlikely to contribute substantially to the levels of heterogeneity detected.

We next explored the extent to which transmission occurs within the household as a possible explanation for the observed overdispersion. Household-level risk has been identified in other settings whereby individuals residing within an infected house are more likely to also be or become infected [17, 25, 36, 37]. However, it is not known whether the increased burden is due to the aggregation of factors that increase risk of infection or because the household itself is the unit of transmission. In this setting, we observed sequential infections within households where new household members became infected in the month after the initial introduced infection. This pattern suggests that within-household transmission is plausible and supports the use of reactive case detection strategies, where households of any confirmed infection are visited and screened and/or treated for malaria to capture additional cases expected within households of index cases [17, 38]. Based on the limited spreading pattern observed, including neighboring households or those within a specific radius around index households would not be recommended in this setting. Also, a reactive approach for targeting residual infections within households is not likely to be appropriate in all settings. This is particularly true for those settings where transmission occurs outside of the household, for example, in forests, as is common in Southeast Asia [24, 39]. Furthermore, given the stochastic nature of infections across all villages, a reactive approach may not contribute to a reduction of transmission but may contribute to infections averted in household members, particularly if a drug with a longer prophylactic period is used. Given that all villages in this study are capable of supporting transmission and would therefore be considered as “active” according to the World Health Organization (WHO) definition of foci, one could argue that targeting the whole village population with interventions may be more appropriate as a way to accelerate malaria elimination [40].

It is possible that hotspots do exist and fuel transmission within foci, but it was not observed in this setting. It is unlikely that infections were missed, as routine sampling occurred every month during the transmission season with a study nurse capturing episodes between regular visits. Although the monthly time step was selected as it would account for the intrinsic and extrinsic incubation periods, it is possible that this temporal scale was not optimal or the monthly aggregated datasets too small to detect the spreading of infections between households. The treatment of detectable infections as part of both the passive and active screenings may have altered or masked spatio-temporal patterns. However, the expected rate of treatments required to interrupt transmission is much higher than was administered as part of routine surveillance. Secondly, the spatio-temporal patterns observed pre- and post-MDA were similar, despite the magnitude of transmission intensity being lower in the second year. Therefore, the role of treatment likely had a minimal impact on the ability to observe any patterns. Alternative spatio-modeling approaches such as point pattern or dispersion models may have yielded different results. However, the number of points per village limited any point-based analysis, and understanding whether infections cluster would not directly address the question of interest. Incorporating the parasite genetic data into this analysis to track infections within and between households may help us understand the extent of within- and between-household transmission dynamics [41]. The detailed genetic data required for this analysis was not available. However, recent work supports the notion of microepidemiological clustering of parasite strains [33]. Next, the non-response bias experienced in this 2-year cohort may have masked any hotspot dynamics. It is possible that the individuals missed could have better illustrated any spreading between households. However, the participation rate across all villages was reasonably high and was consistent between villages, so although possible, we do not consider this as likely.

Conclusions

Approaches for more efficient targeting of malaria control and elimination activities have shifted to incorporating spatial dynamics of transmission and identifying lingering foci. Although hotspots fueling malaria transmission within a village or foci are biologically plausible, the limited evidence in field settings puts their role in sustaining transmission into question. The results presented here further support this shift in thinking [40, 42]. This population-level cohort in 12 villages across The Gambia showed that there is considerable heterogeneity in transmission both within and between study villages. Our results suggest that spatio-temporal patterns of malaria risk are stochastic at all endemicities and are inconsistent with the idea of hotspots fueling malaria transmission. Transmission was more likely to occur within households in this setting, supporting the use of reactive case detection strategies targeting the household only or to target the entire village as a focus, but not an approach targeting hotspots with the goal of interrupting transmission from high to low burden areas.

Abbreviations

IQR:

Interquartile range

MDA:

Mass drug administration

PCR:

Polymerase chain reaction

RDT:

Rapid diagnostic test

References

  1. Woolhouse MEJ, Dye C, Etard J-F, Smith T, Charlwood JD, Garnett GP, Hagan P, Hii JLK, Ndhlovu PD, Quinnell RJ, et al. Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci U S A. 1997;94:338–42.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Bejon P, Williams TN, Nyundo C, Hay SI, Benz D, Gething PW, Otiende M, Peshu J, Bashraheil M, Greenhouse B, et al. A micro-epidemiological analysis of febrile malaria in coastal Kenya showing hotspots within hotspots. eLife. 2014;3:e02130.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bjorkman A, Cook J, Sturrock H, Msellem M, Ali A, Xu W, Molteni F, Gosling R, Drakeley C, Martensson A. Spatial distribution of falciparum malaria infections in Zanzibar: implications for focal drug administration strategies targeting asymptomatic parasite carriers. Clin Infect Dis. 2017;64(9):1236–43.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Mogeni P, Omedo I, Nyundo C, Kamau A, Noor A, Bejon P, Hotspot Group A. Effect of transmission intensity on hotspots and micro-epidemiology of malaria in sub-Saharan Africa. BMC Med. 2017;15(1):121.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Gaudart J, Poudiougou B, Dicko A, Ranque S, Toure O, Sagara I, Diallo M, Diawara S, Ouattara A, Diakite M, et al. Space-time clustering of childhood malaria at the household level: a dynamic cohort in a Mali village. BMC Public Health. 2006;6:286.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Ernst KC, Adoka SO, Kowuor DO, Wilson ML, John CC. Malaria hotspot areas in a highland Kenya site are consistent in epidemic and non-epidemic years and are associated with ecological factors. Malar J. 2006;5:78.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Cook J, Kleinschmidt I, Schwabe C, Nseng G, Bousema T, Corran PH, Riley EM, Drakeley CJ. Serological markers suggest heterogeneity of effectiveness of malaria control interventions on Bioko Island, Equatorial Guinea. PLoS One. 2011;6(9):e25137.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, Ghani A, Drakeley C, Gosling R. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012;9(1):e1001165.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Kangoye DT, Noor A, Midega J, Mwongeli J, Mkabili D, Mogeni P, Kerubo C, Akoo P, Mwangangi J, Drakeley C, et al. Malaria hotspots defined by clinical malaria, asymptomatic carriage, PCR and vector numbers in a low transmission area on the Kenyan coast. Malar J. 2016;15(1):213.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Paull SH, Song S, McClure KM, Sackett LC, Kilpatrick AM, Johnson PT. From superspreaders to disease hotspots: linking transmission across hosts and space. Front Ecol Environ. 2012;10(2):75–82.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Stresman GH. Beyond temperature and precipitation: ecological risk factors that modify malaria transmission. Acta Trop. 2010;116(3):167–72.

    Article  PubMed  Google Scholar 

  12. Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, Tatem AJ, Hay SI. A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011;10:378.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Paajmans KP, Read AF, Thomas MB. Understanding the link between malaria risk and climate. Proc Natl Acad Sci U S A. 2009;106(33):13844–9.

    Article  Google Scholar 

  14. Hay SI, Snow RW, Rogers DJ. Predicting malaria seasons in Kenya using multitemporal meteorological satellite sensor data. Trans R Soc Trop Med Hyg. 1998;92:12–20.

    Article  PubMed  CAS  Google Scholar 

  15. Bannister-Tyrrell M, Verdonck K, Hausmann-Muela S, Gryseels C, Muela Ribera J, Peeters Grietens K. Defining micro-epidemiology for malaria elimination: systematic review and meta-analysis. Malar J. 2017;16(1):164.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Tucker Lima JM, Vittor A, Rifai S, Valle D. Does deforestation promote or inhibit malaria transmission in the Amazon? A systematic literature review and critical appraisal of current evidence. Phil Trans R Soc London Ser B, Biol Sci. 2017;372(1722):20160125.

    Article  Google Scholar 

  17. Sturrock HJ, Novotny JM, Kunene S, Dlamini S, Zulu Z, Cohen JM, Hsiang MS, Greenhouse B, Gosling RD. Reactive case detection for malaria elimination: real-life experience from an ongoing program in Swaziland. PLoS One. 2013;8(5):e63830.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Smith DL, Dushoff J, McKenzie FE. The risk of a mosquito-borne infection in a heterogeneous environment. PLoS Biol. 2004;2(11):e368.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Smith DL, McKenzie FE, Snow RW, Hay SI. Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007;5(3):e42.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Smith DL, Perkins TA, Reiner RC Jr, Barker CM, Niu T, Chaves LF, Ellis AM, George DB, Le Menach A, Pulliam JR, et al. Recasting the theory of mosquito-borne pathogen transmission dynamics and control. Trans R Soc Trop Med Hyg. 2014;108(4):185–97.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Sturrock HJ, Hsiang MS, Cohen JM, Smith DL, Greenhouse B, Bousema T, Gosling RD. Targeting asymptomatic malaria infections: active surveillance in control and elimination. PLoS Med. 2013;10(6):e1001467.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Bousema T, Stresman G, Baidjoe AY, Bradley J, Knight P, Stone W, Osoti V, Makori E, Owaga C, Odongo W, et al. The impact of hotspot targeted interventions on malaria transmission in Rachuonyo South District in the western Kenyan highlands: a cluster-randomized controlled trial. PLoS Med. 2016;13(4):e1001993.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Stresman GH, Giorgi E, Baidjoe A, Knight P, Odongo W, Owaga C, Shagari S, Makori E, Stevenson J, Drakeley C, et al. Impact of metric and sample size on determining malaria hotspot boundaries. Sci Rep. 2017;7:45849.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. van Eijk AM, Ramanathapuram L, Sutton PL, Kanagaraj D, Sri Lakshmi Priya G, Ravishankaran S, Asokan A, Tandel N, Patel A, Desai N, et al. What is the value of reactive case detection in malaria control? A case-study in India and a systematic review. Malar J. 2016;15(1):67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Mwesigwa J, Achan J, Di Tanna GL, Affara M, Jawara M, Worwui A, Hamid-Adiamoh M, Kanuteh F, Ceesay S, Bousema T, et al. Residual malaria transmission dynamics varies across The Gambia despite high coverage of control interventions. PLoS One. 2017;12(11):e0187059.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Tusting LS, Bousema T, Smith DL, Drakeley C. Measuring changes in Plasmodium falciparum transmission: precision, accuracy and costs of metrics. Adv Parasitol. 2014;84:151–208.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Smith DL, Dushoff J, Snow RW, Hay SI. The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005;438(7067):492–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Smith DL, Guerra CA, Snow RW, Hay SI. Standardizing estimates of the Plasmodium falciparum parasite rate. Malar J. 2007;6:131.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Giorgi E, Diggle P. PrevMap: an R package for prevalence mapping. J Stat Softw. 2017;78:1–29.

    Article  Google Scholar 

  30. Bivand R, Rundel C, Pebesma E, Stuetz R, Hufthammer KO: Package 'rgeos'. The Comprehensive R Archive Network (CRAN) 2017(0.3–23).

  31. Hardy A, Mageni Z, Dongus S, Killeen G, Macklin MG, Majambare S, Ali A, Msellem M, Al-Mafazy AW, Smith M, et al. Mapping hotspots of malaria transmission from pre-existing hydrology, geology and geomorphology data in the pre-elimination context of Zanzibar, United Republic of Tanzania. Parasites & vectors. 2015;8:41.

    Article  Google Scholar 

  32. Chang HH, Childs LM, Buckee CO. Variation in infection length and superinfection enhance selection efficiency in the human malaria parasite. Sci Rep. 2016;6:26370.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Omedo I, Mogeni P, Bousema T, Rockett K, Amambua-Ngwa A, Oyier I, Stevenson JC, Baidjoe AY, de Villiers EP, Fegan G, et al. Micro-epidemiological structuring of Plasmodium falciparum parasite populations in regions with varying transmission intensities in Africa. Wellcome Open Research. 2017;2:10.

    PubMed  PubMed Central  Google Scholar 

  34. Karl S, White MT, Milne GJ, Gurarie D, Hay SI, Barry AE, Felger I, Mueller I. Spatial effects on the multiplicity of Plasmodium falciparum infections. PLoS One. 2016;11(10):e0164054.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Stone W, Goncalves BP, Bousema T, Drakeley C. Assessing the infectious reservoir of falciparum malaria: past and future. Trends Parasitol. 2015;31(7):287–96.

    Article  PubMed  Google Scholar 

  36. Stresman GH, Baidjoe AY, Stevenson J, Grignard L, Odongo W, Owaga C, Osoti V, Makori E, Shagari S, Marube E, et al. Focal screening to identify the subpatent parasite reservoir in an area of low and heterogeneous transmission in the Kenya highlands. J Infect Dis. 2015;212(11):1768–77.

    Article  PubMed  Google Scholar 

  37. Clark TD, Greenhouse B, Njama-Meya D, Nzarubara B, Maiteki-Sebuguzi C, Staedke SG, Seto E, Kamya MR, Rosenthal PJ, Dorsey G. Factors determining the heterogeneity of malaria incidence in children in Kampala, Uganda. J Infect Dis. 2008;198(3):393–400.

    Article  PubMed  Google Scholar 

  38. Stresman GH, Kamanga A, Moono P, Hamapumbu H, Mharakurwa S, Kobayashi T, Moss WJ, Shiff C. A method of active case detection to target reservoirs of asymptomatic malaria and gametocyte carriers in a rural area in Southern Province, Zambia. Malaria J. 2010;9:265.

    Google Scholar 

  39. Gerardin J, Bever CA, Bridenbecker D, Hamainza B, Silumbe K, Miller JM, Eisele TP, Eckhoff PA, Wenger EA. Effectiveness of reactive case detection for malaria elimination in three archetypical transmission settings: a modelling study. Malar J. 2017;16(1):248.

    Article  PubMed  PubMed Central  Google Scholar 

  40. WHO. A framework for malaria elimination. Geneva: World Health Organization; 2017.

  41. Reiner RC, Le Manach A, Kunene S, Ntshalintshali N, Hsiang MS, Perkins TA, Greenhouse B, Tatem AJ, Cohen JM, Smith DL. Mapping residual transmission for malaria elimination. eLife. 2015;4:e9520

  42. Lessler J, McKay HS, Moore SM, Azman AS. What is a hotspot anyway? Am J Trop Med Hygiene. 2017;96(6):1270–3.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank the teams involved in data collection as well as those in the laboratory who contributed to analyzing the generated samples. The authors would also like to acknowledge Dr. Nunu Sepulveda, who provided valuable input on the data and statistical analysis.

Funding

This study was funded by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement (MC_EX_MR/J002364/1) awarded to UD. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

The datasets analyzed during the current study and the R code are available through the corresponding author upon reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

JM, JA, GLdT, JvG, and UD were responsible for the cohort study design. JM, AW, and MJ were involved in data collection and database management. Analysis and interpretation of the data was performed by GS, EG, JM, CD, TB, and UD. GS wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Gillian H. Stresman.

Ethics declarations

Ethics approval and consent to participate

Informed written consent was obtained from all study participants or their guardians upon recruitment into the study cohort. This study was approved by the Scientific Coordinating Committee of the Medical Research Council, Gambia Unit (SSC 1318) and The Gambia Government/MRC Joint Ethics Committee.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Description of the geostatistical model including outputs. (PDF 181 kb)

Additional file 3:

Results of analysis stratified by whether an infection was symptomatic or asymptomatic. Symptomatic infections were those individuals with fever and a confirmed infection by rapid diagnostic test (RDT), whereas asymptomatic infections were those who were febrile but RDT negative and positive by polymerase chain reaction (PCR) or were afebrile and positive for malaria by PCR. The distribution of prevalence of infection type per month per village as well as the predicted monthly prevalence of each infection type per village are shown. (DOCX 4478 kb)

Additional file 4:

Heat maps for all villages, removing households with no infections and individuals with 5 or fewer time points sampled. Each panel represents a household, and each row within the panel represents an individual residing in that household. Individuals are ordered by increasing age with the youngest on top within each panel. Each column within each grid represents a sampling month starting in June 2013 through December 2014. (DOCX 4496 kb)

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Stresman, G.H., Mwesigwa, J., Achan, J. et al. Do hotspots fuel malaria transmission: a village-scale spatio-temporal analysis of a 2-year cohort study in The Gambia. BMC Med 16, 160 (2018). https://doi.org/10.1186/s12916-018-1141-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-018-1141-4

Keywords