- Open Access
Using the SaTScan method to detect local malaria clusters for guiding malaria control programmes
Malaria Journal volume 8, Article number: 68 (2009)
Mpumalanga Province, South Africa is a low malaria transmission area that is subject to malaria epidemics. SaTScan methodology was used by the malaria control programme to detect local malaria clusters to assist disease control planning. The third season for case cluster identification overlapped with the first season of implementing an outbreak identification and response system in the area.
SaTScan™ software using the Kulldorf method of retrospective space-time permutation and the Bernoulli purely spatial model was used to identify malaria clusters using definitively confirmed individual cases in seven towns over three malaria seasons. Following passive case reporting at health facilities during the 2002 to 2005 seasons, active case detection was carried out in the communities, this assisted with determining the probable source of infection. The distribution and statistical significance of the clusters were explored by means of Monte Carlo replication of data sets under the null hypothesis with replications greater than 999 to ensure adequate power for defining clusters.
Results and discussion
SaTScan detected five space-clusters and two space-time clusters during the study period. There was strong concordance between recognized local clustering of cases and outbreak declaration in specific towns. Both Albertsnek and Thambokulu reported malaria outbreaks in the same season as space-time clusters. This synergy may allow mutual validation of the two systems in confirming outbreaks demanding additional resources and cluster identification at local level to better target resources.
Exploring the clustering of cases assisted with the planning of public health activities, including mobilizing health workers and resources. Where appropriate additional indoor residual spraying, focal larviciding and health promotion activities, were all also carried out.
Malaria is the most important parasitic disease of humans. Over three billion people live in malarious areas and the disease causes over 500 million cases with one to three million deaths per year [1, 2]. An estimated one hundred million people in Africa are at risk of malaria epidemics . In common with most vector-borne infectious diseases, malaria is heterogeneous in its distribution in time and space [4–6], and incidence can vary greatly between districts, towns and villages. This heterogeneity is affected by patterns of malaria vector distribution, human-vector contact, human host behavioural factors, house construction, and malaria prevention methods used [5, 7–10].
Characterization of malaria heterogeneity may allow prioritization of risk areas and allow focused control interventions. The ability to identify localized malaria clusters in remote highland areas of east Africa facilitated early intervention in the absence of early warning systems, as these malaria "hotspots" remained constant in epidemic and non epidemic years .
A number of models have been developed for describing malaria spatial distribution and seasonality [12–15] transmission [8, 16–18], mosquito distribution [19, 20] and risk factors associated with malaria [21–23]. However many malaria endemic environments are resource limited and tools to assist decision support need to take account of these limitations.
Malaria remains a public health problem in north-eastern South Africa, in the low altitude provinces of KwaZulu-Natal, Limpopo and Mpumalanga [24, 25]. Malaria risk is low compared to other hyper- and holo-endemic areas of sub-Saharan Africa and naturally acquired immunity does not develop in the local population .
The burden of malaria is well described in this region as definitive diagnosis using rapid diagnostic tests (RDTs) and mandatory reporting of malaria cases is universally practised in the public health system, which manages the vast majority of malaria cases [27, 28]. All confirmed malaria cases are entered into a computerized malaria surveillance system.
This paper reports on the analysis of malaria notification data from 2002–2005 from Mpumalanga Province to determine local clustering of cases and how that was used to direct local control efforts and enhance the district outbreak identification and response system .
Mpumalanga Province, which borders Mozambique and Swaziland, is a predominantly rural area with a population of approximately four million people. Malaria occurrence in this area is seasonal during the wet and humid months of October to April. Ninety-five percent of malaria cases are due to Plasmodium falciparum infection. This area has historically experienced malaria outbreaks and epidemics with relatively high mortality [29, 30]. The main vector in this region is Anopheles arabiensis and vector control is predominantly by indoor residual spraying (IRS) with DDT [31, 32].
The seven rural towns that experience the highest malaria risk in Mpumalanga were included in this study (Figure 1), with malaria incidence ranging from 30 to 61 per 1,000 persons between 1997 and 2005.
Original town maps acquired in 2000 were upgraded from digitized aerial photographs produced in 2002 (Cadnet, Nelspruit, South Africa) and updates were converted into MapInfo version 6.4 files (MapInfo Corporation, New York, USA) with Geomedia software (Symmetry Systems Inc., New York, USA). Unique stand (household) numbers were allocated during the digitizing process. A stand may include one or more structures belonging to a family unit on a designated piece of land.
Active malaria case investigation
Malaria case investigators followed up all malaria cases reported at local health facilities for these towns as part of malaria control program activities for the 2002/2003, 2003/2004 and 2004/2005 seasons. A malaria case is defined as a person diagnosed either by positive rapid test or blood slide. An investigator would identify the probable source of infection by confirming the travel history of the patient recorded on the notification.
Spatial and temporal clusters
SaTScan™ software, version 5.1.3, using the Kulldorf method of retrospective space-time permutation and the Bernoulli purely spatial model [34, 35], was used by malaria control programme operational staff to detect malaria clusters in individual towns in the three seasons under investigation and over the combined time period. Individual malaria cases were used and recorded against source of infection households. Households who did not seek health care for malaria or tested negative during active case detection were used as controls for the Bernoulli method. More than one case can be reported per household. A household in these areas is typically defined as a family unit with a single land owner.
The circular scan statistic is isotopic with respect to the rotation of the geographical area http://www.satscan.org. This method has previously been validated for plotting and understanding local malaria time-space-clusters [21, 36–39].
Observed cases in a cluster were compared to the distribution of expected cases if spatial and temporal locations of all cases were independent. The model adjusts for entirely spatial or entirely temporal clusters. With spatial adjustment, time remained dormant and during temporal analysis seasons were considered. The distribution and statistical significance of the clusters were explored by means of Monte Carlo replication of data sets under the null hypothesis with replications greater than 999 to ensure adequate power for defining clusters http://www.satscan.org. Clusters were prioritized for public health action according to their statistical significance. Hard copy town maps showing cluster areas were distributed to malaria field staff investigating cases allowing the coordination of intervention efforts within communities.
Four hundred and twenty two malaria cases were notified during the three seasons from 341 households across the seven towns. Malaria incidence differed significantly between the towns during the three seasons (χ2 = 6.442, p = 0.040) (Figure 2).
SaTScan analysis detected a number of clusters during the study period (Table 1). Albertsnek produced two space-clusters over the three-season period using the Bernoulli model, one in the northern (log likelihood ratio = 12.308, p = 0.003) and one in the south eastern part of town (log likelihood ratio = 12.187, p = 0.003) (Figure 3A). Three additional space-clusters were observed; in Goba (log likelihood ratio = 12.226, p = 0.001), Mbuzini (log likelihood ratio = 22.372, P = 0.001) and Thambokulu (Figure 4A) (log likelihood ratio = 22.372, p = 0.001). Only two towns, Albertsnek (test statistic = 5.548, p = 0.007) and Thambokulu (test statistic = 3.668, p = 0.004), had space-time clusters and both were during the 2004/2005 season (Figures 3B and 4B).
Cluster detection to guide control activities
The third season for case cluster identification, 2004/2005, overlapped with the first season of implementation of an outbreak identification and response system in the Nkomazi area . Thambokulu and Albertsnek triggered outbreak declarations during this season (Table 1). Analysis demonstrated clusters in both towns during this season. Real time cluster detection, monitoring as cases occur, triggered at the same time as the binomial outbreak identification thresholds .
The physical plotting of cases by household during this season (2004/2005) as part of the cluster identification activities at district level and on hard copy maps provided a visual distribution of risk in the towns and assisted with the timely planning of public health activities especially logistics, including active case detection, early diagnosis and treatment of positive cases in the areas of the clusters, additional indoor residual spraying, focal larviciding and health promotion activities.
Application of the SaTScan method successfully identified malaria clusters and clearly demonstrated malaria risk heterogeneity at local level. Four towns in this study experienced spatial clusters and two produced space-time clusters in the 2004/2005 malaria season, the latter resulting in targeted local control efforts. The high rates of treatment-seeking behaviour at primary health care level in Mpumalanga  combined with the use of passive notification and active case detection, strengthens disease surveillance and provides good quality data for cluster identification. SaTScan cluster identification has also proven valuable for targeting control strategies in the Kenyan highlands . Over a ten week period during a 2002 epidemic, spatial targeting of vector control interventions reduced the abundance of Anopheles mosquitoes. The investigation of causes for clusters within the areas reported here was not explored. The hypo-endemic nature of malaria in these areas does predict the unstable patterns of disease occurrence and more specifically outbreaks that can be clustered. The proximity of these areas to higher malaria risk areas such as Mozambique allows for the import of parasites and eventually local transmission.
Spatial clustering of infectious disease is enjoying renewed interest, particularly in areas of limited resources. Statistical methods have been used to investigate spatial clustering of dengue  encephalitis  and sleeping sickness but the application to malaria has been limited .
Gaudart et al  compared the oblique decision tree model, a complex statistical technique with Kulldorff's SaTScan™ cluster technique, which was used in this study. Gaudart et al  produced similar results using both methods in a village in West Africa to identify malaria risk clusters. This investigation confirmed the usefulness of the Kulldorff's scan statistic [44–46]. Due to the circular isotopic technique of Kulldorff's SaTScan™ it is a useful tool to detect clusters but has limitations on detecting irregular shaped clusters due to its fixed scan window [47, 48].
The seven Mpumalanga towns under cluster surveillance were included in a malaria outbreak identification and response system based on formal case reporting . There was strong concordance between recognized local clustering and outbreak identification in specific towns, with Albertsnek and Thambokulu both reporting malaria outbreaks in the same season as the time-space clusters. This synergy allows mutual validation of the two systems in confirming outbreaks, which demand additional resources, and cluster identification that could better target these resources within the affected town.
Guerra CA, Hay SI, Lucioparedes LS, Gikandi PW, Tatem AJ, Noor AM, Snow RW: Assembling a global database of malaria parasite prevalence for the Malaria Atlas Project. Malar J. 2007, 6: 17-10.1186/1475-2875-6-17.
Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434: 214-217. 10.1038/nature03342.
World Health Organization: Malaria epidemics: forecasting, prevention, early detection and control. From policy to practice. 2003, Leysin, Switzerland, WHO
Carter R, Mendis KN, Roberts D: Spatial targeting of interventions against malaria. Bull World Health Organ. 2000, 78: 1401-1411.
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.
Greenwood BM: The microepidemiology of malaria and its importance to malaria control. Trans R Soc Trop Med Hyg. 1989, 83: 25-29. 10.1016/0035-9203(89)90599-3.
Snow RW, Gouws E, Omumbo J, Rapuoda B, Craig MH, Tanser FC, le Sueur D, Ouma J: Models to predict the intensity of Plasmodium falciparum transmission: applications to the burden of disease in Kenya. Trans R Soc Trop Med Hyg. 1998, 92: 601-606. 10.1016/S0035-9203(98)90781-7.
Thomas CJ, Lindsay SW: Local-scale variation in malaria infection amongst rural Gambian children estimated by satellite remote sensing. Trans R Soc Trop Med Hyg. 2000, 94: 159-163. 10.1016/S0035-9203(00)90257-8.
Thompson R, Begtrup K, Cuamba N, Dgedge M, Mendis C, Gamage-Mendis A, Enosse SM, Barreto J, Sinden RE, Hogh B: The Matola malaria project: a temporal and spatial study of malaria transmission and disease in a suburban area of Maputo, Mozambique. Am J Trop Med Hyg. 1997, 57: 550-559.
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.
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-10.1186/1475-2875-5-78.
Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.
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. 10.1016/S0035-9203(98)90936-1.
Omumbo J, Ouma J, Rapuoda B, Craig MH, le Sueur D, Snow RW: Mapping malaria transmission intensity using geographical information systems (GIS): an example from Kenya. Ann Trop Med Parasitol. 1998, 92: 7-21.
Tanser FC, Sharp B, le Sueur D: Potential effect of climate change on malaria transmission in Africa. Lancet. 2003, 362: 1792-1798. 10.1016/S0140-6736(03)14898-2.
Kleinschmidt I, Bagayoko M, Clarke GP, Craig M, le SD: A spatial statistical approach to malaria mapping. Int J Epidemiol. 2000, 29: 355-361. 10.1093/ije/29.2.355.
Snow RW, Peshu N, Forster D, Bomu G, Mitsanze E, Ngumbao E, Chisengwa R, Schellenberg JR, Hayes RJ, Newbold CI, Marsh K: Environmental and entomological risk factors for the development of clinical malaria among children on the Kenyan coast. Trans R Soc Trop Med Hyg. 1998, 92: 381-385. 10.1016/S0035-9203(98)91056-2.
Yang G, Zhou X, Malone JB, McCarroll JC, Wang T, Liu J, Gao Q, Zhang X, Hong Q, Sun L: GIS prediction model of malaria transmission in Jiangsu province. Zhonghua Yu Fang Yi Xue Za Zhi. 2002, 36: 103-105.
Kitron U, Pener H, Costin C, Orshan L, Greenberg Z, Shalom U: Geographic information system in malaria surveillance: mosquito breeding and imported cases in Israel, 1992. Am J Trop Med Hyg. 1994, 50: 550-556.
Sithiprasasna R, Linthicum KJ, Liu GJ, Jones JW, Singhasivanon P: Use of GIS-based spatial modeling approach to characterize the spatial patterns of malaria mosquito vector breeding habitats in northwestern Thailand. Southeast Asian J Trop Med Public Health. 2003, 34: 517-528.
Bostoen K, Chalabi Z: Optimization of household survey sampling without sample frames. Int J Epidemiol. 2006, 35: 751-755. 10.1093/ije/dyl019.
Gemperli A, Vounatsou P, Kleinschmidt I, Bagayoko M, Lengeler C, Smith T: Spatial patterns of infant mortality in Mali: the effect of malaria endemicity. Am J Epidemiol. 2004, 159: 64-72. 10.1093/aje/kwh001.
Hu H, Singhasivanon P, Salazar NP, Thimasarn K, Li X, Wu Y, Yang H, Zhu D, Supavej S, Looarecsuwan S: Factors influencing malaria endemicity in Yunnan Province, PR China (analysis of spatial pattern by GIS). Geographical Information System. Southeast Asian J Trop Med Public Health. 1998, 29: 191-200.
leSueur D, Sharp BL, Gouws E, Ngxongo S: Malaria in South Africa. S Afr Med J. 1996, 86: 936-939.
Sharp BL, le SD: Malaria in South Africa – the past, the present and selected implications for the future. S Afr Med J. 1996, 86: 83-89.
Richie TL, Saul A: Progress and challenges for malaria vaccines. Nature. 2002, 415: 694-701. 10.1038/415694a.
Durrheim DN, La Grange JJ, Govere J, Mngomezulu NM: Accuracy of a rapid immunochromatographic card test for Plasmodium falciparum in a malaria control programme in South Africa. Trans R Soc Trop Med Hyg. 1998, 92: 32-33. 10.1016/S0035-9203(98)90944-0.
Durrheim DN, Harris BN, Speare R, Billinghurst K: The use of hospital-based nurses for the surveillance of potential disease outbreaks. Bull World Health Organ. 2001, 79: 22-27.
Coleman M, Coleman M, Mabuza A, Kok G, Coetzee M, Durrheim D: Evaluation of an operational malaria outbreak identification and response system in Mpumalanga Province, South Africa. Malar J. 2008, 7: 69-10.1186/1475-2875-7-69.
Govere J, Durrheim DN, Baker L, Hunt R, Coetzee M: Efficacy of three insect repellents against the malaria vector Anopheles arabiensis. Med Vet Entomol. 2000, 14: 441-444. 10.1046/j.1365-2915.2000.00261.x.
Kager PA: Malaria control: constraints and opportunities. Trop Med Int Health. 2002, 7: 1042-1046. 10.1046/j.1365-3156.2002.00981.x.
Kapp C: Hazard or help?. Lancet. 2004, 364: 1113-1114. 10.1016/S0140-6736(04)17122-5.
Booman M, Durrheim DN, La Grange K, Martin C, Mabuza AM, Zitha A, Mbokazi FM, Fraser C, Sharp BL: Using a geographical information system to plan a malaria control programme in South Africa. Bull World Health Organ. 2000, 78: 1438-1444.
Kulldorff M, Athas WF, Feurer EJ, Miller BA, Key CR: Evaluating cluster alarms: a space-time scan statistic and brain cancer in Los Alamos, New Mexico. Am J Public Health. 1998, 88: 1377-1380. 10.2105/AJPH.88.9.1377.
Kulldorff M, Heffernan R, Hartman J, Assuncao R, Mostashari F: A space-time permutation scan statistic for disease outbreak detection. PLoS Med. 2005, 2: e59-10.1371/journal.pmed.0020059.
Boscoe FP, McLaughlin C, Schymura MJ, Kielb CL: Visualization of the spatial scan statistic using nested circles. Health Place. 2003, 9: 273-277. 10.1016/S1353-8292(02)00060-6.
kulldorf m, Feuer EJ, Miller BA, Freeman LS: Breast Cancer in northeastern United States: a geographical analysis. Am J Epidemiol. 1997, 146: 161-170.
Nkhoma ET, Hsu CE, Hunt VI, Harris AM: Detecting spatiotemporal clusters of accidental poisoning mortality among texas counties, United States, 1980–2001. Int J Health Geogr. 2004, 3: 25-10.1186/1476-072X-3-25.
Odoi A, Martin SW, Michel P, Middleton D, Holt J, Wilson J: Investigation of clusters of giardiasis using GIS and a spatial scan statistic. Int J Health Geogr. 2004, 3: 11-10.1186/1476-072X-3-11.
Govere J, Durrheim D, La GK, Mabuza A, Booman M: Community knowledge and perceptions about malaria and practices influencing malaria control in Mpumalanga Province, South Africa. S Afr Med J. 2000, 90: 611-616.
Morrison AC, Getis A, Santiago M, Rigau-Perez JG, Reiter P: Exploratory space-time analysis of reported dengue cases during an outbreak in Florida, Puerto Rico, 1991–1992. Am J Trop Med Hyg. 1998, 58: 287-298.
Chadee DD, Kitron U: Spatial and temporal patterns of imported malaria cases and local transmission in Trinidad. Am J Trop Med Hyg. 1999, 61: 513-517.
Gaudart J, Poudiougou B, Ranque S, Doumbo O: Oblique decision trees for spatial pattern detection: optimal algorithm and application to malaria risk. BMC Med Res Methodol. 2005, 5: 22-10.1186/1471-2288-5-22.
Hjalmars U, Kulldorff M, Gustafsson G, Nagarwalla N: Childhood leukaemia in Sweden: using GIS and a spatial scan statistic for cluster detection. Stat Med. 1996, 15: 707-715. 10.1002/(SICI)1097-0258(19960415)15:7/9<707::AID-SIM242>3.0.CO;2-4.
Kulldorf M: A spatial scan statistic. Communications in Statistics: Theory and Methods. 1997, 26: 1481-1496. 10.1080/03610929708831995.
Shaheen P, Clart JD, Williams D: Small area analysis: a review and analysis of North American literature. J Health Polit Policy Law. 1987, 12: 741-809. 10.1215/03616878-12-4-741.
Aamodt G, Samuelsen SO, Skrondal A: A simulation study of three methods for detecting disease clusters. Int J Health Geogr. 2006, 5: 15-10.1186/1476-072X-5-15.
Patil GP, Taille C: Geographic and network surveillance via scan statistics for critical area detection. Stat Sci. 2003, 18: 457-465. 10.1214/ss/1081443229.
Special thanks and acknowledgement to the field personnel from the Malaria Control Programme in Mpumalanga province for their dedication and commitment to the cluster and outbreak activities.
The authors declare that they have no competing interests.
MaCol carried out all the investigative work and analysis under the guidance of MCoe and DND Field work would not have been feasible without the assistance of AMM and GK MiCol drafted the manuscript and corrected subsequent drafts.
About this article
Cite this article
Coleman, M., Coleman, M., Mabuza, A.M. et al. Using the SaTScan method to detect local malaria clusters for guiding malaria control programmes. Malar J 8, 68 (2009) doi:10.1186/1475-2875-8-68
- Malaria Case
- Indoor Residual Spray
- Malaria Incidence
- Malaria Control Programme