The relationship between the Plasmodium falciparum parasite ratio in childhood and climate estimates of malaria transmission in Kenya
© Omumbo et al; licensee BioMed Central Ltd. 2004
Received: 17 May 2004
Accepted: 17 June 2004
Published: 17 June 2004
Plasmodium falciparum morbid and fatal risks are considerably higher in areas supporting parasite prevalence ≥25%, when compared with low transmission areas supporting parasite prevalence below 25%. Recent descriptions of the health impacts of malaria in Africa are based upon categorical descriptions of a climate-driven fuzzy model of suitability (FCS) for stable transmission developed by the Mapping Malaria Risk in Africa collaboration (MARA).
An electronic and national search was undertaken to identify community-based parasite prevalence surveys in Kenya. Data from these surveys were matched using ArcView 3.2 to extract spatially congruent estimates of the FCS values generated by the MARA model. Levels of agreement between three classes used during recent continental burden estimations of parasite prevalence (0%, >0 – <25% and ≥25%) and three classes of FCS (0, >0 – <0.75 and ≥0.75) were tested using the kappa (k) statistic and examined as continuous variables to define better levels of agreement.
Two hundred and seventeen independent parasite prevalence surveys undertaken since 1980 were identified during the search. Overall agreement between the three classes of parasite prevalence and FCS was weak although significant (k = 0.367, p < 0.0001). The overall correlation between the FCS and the parasite ratio when considered as continuous variables was also positive (0.364, p < 0.001). The margins of error were in the stable, endemic (parasite ratio ≥25%) class with 42% of surveys represented by an FCS <0.75. Reducing the FCS value criterion to ≥0.6 improved the classification of stable, endemic parasite ratio surveys. Zero values of FCS were not adequate discriminators of zero parasite prevalence.
Using the MARA model to categorically distinguish populations at differing intensities of malaria transmission in Kenya may under-represent those who are exposed to stable, endemic transmission and over-represent those at no risk. The MARA approach to defining FCS values of suitability for stable transmission represents our only contemporary continental level map of malaria in Africa but there is a need to redefine Africa's population at risk in accordance with both climatic and non-climatic determinants of P. falciparum transmission intensity to provide a more informed approach to estimating the morbid and fatal consequences of infection across the continent.
In recent years there has been a renaissance in mapping malaria distribution at the national [1–4], continental [5, 6] and global scales . These maps have adopted a variety of approaches and data sources in their construction. The most widely cited, contemporary continental resolution map of Plasmodium falciparum transmission distribution for Africa was developed by the Mapping Malaria Risk in Africa (MARA) collaboration ; http://www.mara.org.za. It has formed the basis for several reports by the Roll Back Malaria partnership [8, 9] and was used during several recent estimates of the pan-African public health malaria burden to identify population's at-risk [10–12]. The MARA model used a fuzzy membership approach, assigning 5 × 5 km areas to a suitability estimate for stable P. falciparum transmission based upon simple rainfall and temperature determinants of the parasite's sporogonic development and mosquito survival. The model did not attempt to define intensities of transmission; rather it determined the likelihood that stable transmission could occur. Using the MARA model, areas of low stable and high stable endemicity across the continent have been distinguished categorically by assuming that the greater the climatic likelihood of stable transmission, the more likely areas would support higher intensities of malaria transmission [11, 12]. In this paper, the validity of these assumptions and the likely margins of error are examined by comparing MARA climate suitability values with empirical P. falciparum parasite prevalence survey data in Kenya.
Plasmodium falciparum prevalence surveys data among children aged 0–15 years in Kenya
A search of published and unpublished literature related to malaria infection prevalence surveys in Kenya was conducted as part of the MARA project [13, 14]. In brief, electronic database searches were performed using Medline® (SilverPlatter International, Boston, MA, USA 2000), Popline® (Johns Hopkins School of Hygiene & Public Health, Baltimore, MD, USA, 2000) and EmBase® (Elsevier Science Little Rock, Arkansas, USA 1999–2000). The following keywords were used in the search: Kenya, malaria, parasite and malaria and Kenya, Plasmodium falciparum, parasite rate, parasite prevalence and malaria transmission. For each publication, bibliographies were cross-referenced to identify additional sources of information from other studies. Where additional details could not be identified through the published sources, authors were contacted to provide more information on geographical location, survey dates and age-specific characteristics of the parasitological data. Postgraduate theses held in the libraries of four departments of the University of Nairobi were also searched (Medical school campuses at the Kenyatta National Hospital and Chiromo, Community Health Department and the Faculty of Science). Annual reports, journals and conference proceedings of national medical research institutes, and non-governmental organisations were reviewed at respective institute's libraries. Results of routine parasite prevalence surveys undertaken by the Ministry of Health's (MoH) Vector Control Department were manually searched in archives at national headquarters and at seven Provincial offices in Kenya.
Using a variety of sources: 1:50 000 scale topographic maps , digital maps of administrative units in Kenya  or public domain digital gazetteers [17, 18], a longitude and latitude was ascertained for each parasitological survey in decimal degrees. These geo-references were imported into a geographical information system platform ArcView 3.2 (ESRI, Redlands, CA, USA), mapped and overlaid on administrative boundary maps for Kenya . The national administrative boundary maps were first warped within ArcView 3.2 to United Nations approved national boundaries  and then used to check for inconsistencies in spatial positioning and to define the spatial coverage of each survey (see below).
All surveys undertaken in Kenya (n = 923) were subjected to a number of selection criteria for inclusion in the analysis. First, to allow for a contemporary assessment of infection risk, surveys were only selected if they were undertaken between 1980 and 2003 (n = 657). The historical data (1927–1979) will be described elsewhere (Omumbo & Snow, in preparation). Second, surveys were excluded if the survey formed part of clinic visits, drug sensitivity testing or included intervention arms of controlled trials. Only total population, randomly sampled or longitudinal community-based surveys were included. Several community-based surveys, particularly those undertaken by the Ministry of Health, provided little detail on the sampling strategy used to select the childhood populations. These surveys were retained but were coded separately in the database to distinguish them from truly random selections should there be differences in estimates of parasite ratio dependent upon sampling methods (n = 74). Third, surveys were excluded if they covered infants only or an age range that extended into adulthood (>=15 years, (n = 12)). Fourthly, a minimum survey sample size of 50 was imposed on the selection to allow for adequate precision in the estimates of infection prevalence , or surveys were excluded if there where no details of the denominator or numerator (n = 21). Finally, repeat surveys by the same investigators within a twenty-four month period were combined into a single estimate. Surveys undertaken by several investigators at different times in the same location were reduced to one estimate by selecting the most recent survey, or the one with the largest sample size.
MARA Fuzzy Climate Suitability Index
The MARA model describes climatic conditions that range from unsuitable (0) to completely suitable (1) for stable P. falciparum transmission . The Fuzzy Climate Suitability (FCS) index is defined by a series of curves
where x is a climate parameter, U is the value of x when conditions are unsuitable, and S is the value of x when conditions are suitable. When S is greater than U the suitability (1-y), increases with x; when S is less than U the suitability y, decreases as x increases. The model defines a monthly increasing curve (S = 22 C, U = 18 C) and decreasing curve (S = 22 C, U = 40 C) for mean diurnal air temperature, a monthly increasing curve (S = 80 mm, U = 0 mm) for rainfall, and a single increasing curve (S = 6 C, U = 4 C) for annual minimum temperature.
The FCS values for each parasite ratio data point were extracted from the FCS model developed by Craig et al. . To provide a spatially congruent FCS value for each parasite prevalence survey, surveys were classified as representing one of five spatial dimensions: First, for surveys representing a single village, the central longitude and latitude was used to create a 2.5 km buffer surrounding the village using ArcView 3.2 (ESRI, Redlands, CA, USA). The size of the buffer was defined by the average flight range of Anopheles gambiae s.l. . Using the MILA Utilities 3.2 update extension within ArcView 3.2 http://www.esri.com/arcscripts, the average of all the 5 × 5 km pixel FCS values within this buffer was extracted to represent the average FCS value for the respective village parasite prevalence survey. The second spatial classification reflected surveys that sampled from several villages but presented the data as a single parasite ratio estimate (n = 5). In this case, a polygon was created to connect the villages and a 5 km buffer created around the polygon to represent the wider spatial sampling. Average FCS values within this buffered polygon were summarised. The third, fourth and fifth spatial criteria corresponded to surveys undertaken at the 5th (sub-location), 4th (location) or 3rd (division) administrative unit levels in Kenya. These spatially distinct polygons were created within ArcView 3.2 using a digital administrative map of Kenya  and a 2.5 km buffer was created around each polygon extent before extraction of the average FCS values.
Data entry and statistical methods
All parasite prevalence survey data were entered twice in Microsoft® Access version 7 (Microsoft Corporation 1989–1996; Seattle, Washington, USA). Data entry errors and range inconsistencies were checked and verified against the original material. Data were summarized for matching to other data sources using Microsoft® Excel 2000 version 9.0 and then analysed using SPSS (SPSS, v10.01, 1999, Chicago, Illinois, USA).
Agreement between parasite prevalence categories among 217 surveys and spatially congruent categorical values of the Fuzzy Climate Suitability (FCS) using the 0.75 threshold for stable endemic malaria.
FCS = 0.00
FCS >0 – <0.75
FCS >= 0.75
Parasite prevalence 0%
Parasite prevalence >0 – <25%
Parasite prevalence >= 25%
The significance of agreement between the three classes of parasite ratio and FCS was tested using the kappa (k) statistic , which is a measure of the agreement between two classifications discounting for the probability that the agreement could be due to chance. Values of kappa range from 0 (no agreement) to 1 (full agreement) and Landis & Koch  suggest the following interpretation of agreement defined by the kappa statistic: poor where k ≤ 0.4, good where 0.4 < k < 0.75 and excellent where kappa ≥0.75.
The simple correlation between the contiguous measures of the FCS values and their respective estimates of the parasite prevalence was also tested. The continuous association was further tested using weighted least-squares regression with the model weighting the dependent variable (parasite ratio) for its precision using the Log10+1 of the standard error of the parasite ratio. Covariates in this model included the end-year of the study (pre-1990, 1990–94, and 1995–2002), the survey sampling methodology (random, longitudinal or total population surveys versus surveys without precise details) and the maximum age range covered in the sample (up to 5 years, up to 10 years and up to 15 years). The proportion of variation in the parasite ratio explained by each variable in the model was calculated by comparing overall proportions of explained variation (R2) between subsequent models with additional variables entered in a stepwise fashion allowing for all variables already in the model.
Agreement between parasite prevalence categories among 217 surveys and spatially congruent categorical values of the Fuzzy Climate Suitability (PCS) using the 0.6 threshold for stable endemic malaria.
FCS = 0.00
FCS >0 – <0.6
FCS >= 0.6
Parasite prevalence 0%
Parasite prevalence >0 – <25%
Parasite prevalence >= 25%
The positive, albeit weak, correlation with measures of the intensity of transmission lends some support for the MARA model's ability to define populations at-risk of differing intensities of malaria transmission. There were too few surveys reporting zero infection prevalence (n = 8) to argue whether the FCS model can correctly distinguish areas of no transmission, however, 5 areas reporting zero prevalence did have a FCS value greater than zero. More striking was the ability of the FCS categories >0 and <0.75 to correctly identify populations at low risk of malaria infection with parasite prevalence's between >0 – <25% (81%: Table 1). During estimations of malaria burden it has been assumed that populations residing in these areas experience much lower risks of malaria-specific morbidity and mortality compared to populations located in areas described by an FCS ≥0.75 [11, 12]. Conversely areas described as supporting stable, endemic transmission with parasite prevalence ≥25% were less well described using the criteria FCS ≥0.75. Only 58% of the parasite prevalence surveys reporting a parasite ratio ≥25% were classified as having a FCS value ≥0.75, the remaining 42% of surveys were classified as having FCS values >0 but less than 0.75. Altering the FCS criteria to ≥0.6 for stable, endemic transmission improved the classifications of parasite prevalence surveys in Kenya with 68% of surveys correctly identified (Table 2).
The results of this study suggest that applying climate suitability modelled estimates of transmission potential to distinguishing populations at differing levels of malaria infection intensity has several limitations. First it may not adequately distinguish populations at no risk of infection. Second, populations exposed to stable, endemic transmission may be poorly defined using criteria of ≥0.75 FCS. These results relate only to Kenya and similar validations are required in other settings in Africa. Nevertheless, for Kenya the disease and mortality burden of malaria will have been under-estimated using the criteria shown in Figure 1. The extent of higher intensity transmission resulting in higher malaria burdens could be wider and the extent of no risk might be more restricted. The model developed by Craig et al.  represents the only available continental scale map upon which to define populations at risk of P. falciparum infection. Our analysis suggests that, at a crude level, it does distinguish communities exposed to differing levels of malaria transmission intensity. Given the large number of assumptions made regarding the paucity of attribute morbidity and mortality data to define continental scale DALY's for malaria, defining the spatial extents of population denominators is only one part of the problem that is at least soluble. New malaria risk models are being developed which go beyond the climatic determinants of P. falciparum transmission and new iterations of malaria risk models will hopefully become more robust with the inclusion of new global scale data on population settlement, land use and ecology.
This study received financial support from The Wellcome Trust, UK; International Development Research Centre, Canada; the South African Medical Research Council, the MARA/ARMA collaboration and the Kenya Medical Research Institute. The Wellcome Trust supports JAO, SIH and RWS as part of their Prize Studentship (#060063), Research Career Development (#069045) and Senior Research Fellow (#058992) programmes respectively. The authors acknowledge the contribution of malaria control personnel and scientists and in Kenya who shared unpublished data. Particular thanks are due to the Division of Malaria Control (Kenya Ministry of Health), Vicki Marsh and Tabitha Mwangi (KEMRI), Dejan Zurovac (MSF-France) and Sporton Kwanya Agak, (KEDHAR). The authors also wish to thank Greg Fegan and Hesbon Ooko for statistical support and Abdisalan M Noor for help with the spatial mapping of the FCS values. This paper is published with the permission of the director of KEMRI.
- Snow RW, Gouws E, Omumbo JA, Rapuoda BA, 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 Roy Soc Trop Med Hyg. 1998, 92: 601-606. 10.1016/S0035-9203(98)90781-7.View ArticlePubMedGoogle Scholar
- Kleinschmidt I, Omumbo JA, Briët O, van de Giesen N, Sogoba N, Mensah N, Windmeijer P, Moussa M, Teuscher T: An empirical malaria distribution map for West Africa. Trop Med Int Health. 2001, 6: 779-786. 10.1046/j.1365-3156.2001.00790.x.View ArticlePubMedGoogle Scholar
- Briët OJT, Gunawardena DM, van der Hoek W, Amerasinghe FP: Sri Lanka malaria maps. Malar J. 2003, 2: 22-32. 10.1186/1475-2875-2-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Hay SI: An overview of remote sensing and geodesy for epidemiology and public health application. Adv Parasitol. 2000, 47: 1-35. 10.1016/S0065-308X(00)47005-3.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- PAHO: Status of malaria programmes in the Americas (based on 2002 data). Report of the 44th directing council, 55th session of the regional committee. Pan American Health Organization (PAHO)/World Health Organization (WHO): Washington D.C. 2003Google Scholar
- Hay SI, Guerra CA, Tatem AJ, Noor AM, Snow RW: The global distribution and population at risk of malaria: past, present and future. Lancet Infect Dis. 2004, 4: 327-336. 10.1016/S1473-3099(04)01043-6.PubMed CentralView ArticlePubMedGoogle Scholar
- World Health Organization: The Africa Malaria Report 2003. World Health Organization/United Nations Children's Fund: Geneva/New York. 2003, 120-Google Scholar
- World Health Organisation: Rolling back malaria. The world health report 1999: making a difference. World Health Organization: Geneva. 1999, 49-63.Google Scholar
- Snow RW, Craig MH, Deichmann U, Marsh K: Estimating mortality an disability due to malaria among Africa's non-pregnant population. Bull World Health Organ. 1999, 77: 624-640.PubMed CentralPubMedGoogle Scholar
- Snow RW, Craig MH, Newton CRJC, Steketee RW: The public health burden of Plasmodium falciparum malaria in Africa: deriving the numbers. Working Paper 11, Disease Control Priorities Project. Fogarty International Center, National Institutes of Health: Bethseda, Maryland. 2003, [http://www.fic.nih.gov/dcpp]Google Scholar
- Snow RW, Omumbo JA: Malaria mortality in sub-Saharan Africa (in press). In Disease and mortality in sub-Saharan Africa. Edited by: Jamison D, Feacham R, Makgoba W, Hofman K, Rogo K. 2004, Oxford University Press & World Bank,Google Scholar
- Omumbo JA, 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.View ArticlePubMedGoogle Scholar
- MARA/ARMA: Towards an atlas of malaria risk in Africa. Fist technical report of the MARA/ARMA collaboration. MARA/ARMA. Durban. 1998Google Scholar
- Directorate of Overseas Surveys: 1:50,000 topographical maps. Series Y731. Survey of Kenya: Nairobi, Kenya. 1971Google Scholar
- Central Bureau of Statistics: 1999 population and housing census: counting our people for development. Volume 1.1 Population distribution by administrative and urban centres. Ministry of Finance & Planning, Nairobi, Kenya. 2001Google Scholar
- GDE Systems: Populated places: GeoName Digital Gazetteer v.1 (CD-ROM). Geographic Information Department (MZ 1211-M), GDE Systems Inc., P.O. Box 509009, San Diego, California. 1995Google Scholar
- World Resources Institute: Topographic data: Africa Data Sampler (CD-ROM). World Resources Institute, USA. 1995Google Scholar
- SALB: Second Administrative Level Boundaries (SALB) Project overviews: concepts, progress and future. 2004, [http://www3.who.int/whosis/gis/salb/salbpo.htm]Google Scholar
- Snow RW, Gilles HM: The epidemiology of malaria. In Bruce Chwatt's Essential Malariology. Edited by: Warrell D, Gilles H. 2002, Arnold Publishers: London, 85-106.Google Scholar
- Service M: Mosquito (Diptera: Culicidae) dispersal – the long and short of it. J Med Entomol. 1997, 34: 579-588.View ArticlePubMedGoogle Scholar
- Cohen J: A co-efficient of agreement for nominal scale. Educ Psychol Measurements. 1960, 20: 37-46.View ArticleGoogle Scholar
- Landis J, Koch G: The measurement of observer agreement for categorical data. Biometrics. 1977, 33: 159-174.View ArticlePubMedGoogle Scholar
- Metselaar D, Van Theil P: Classification of malaria. Trop Geog Med. 1959, 11: 157-161.Google Scholar
- Small J, Goetz SJ, Hay SI: Climatic suitability for malaria transmission in Africa, 1911–1995. Proc Natl Acad USA. 2003, 100: 15341-15345. 10.1073/pnas.2236969100.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.