Spatial variation of Anopheles-transmitted Wuchereria bancrofti and Plasmodium falciparuminfection densities in Papua New Guinea
© Alexander et al; licensee BioMed Central Ltd. 2003
Received: 06 June 2003
Accepted: 14 September 2003
Published: 14 September 2003
The spatial variation of Wuchereria bancrofti and Plasmodium falciparum infection densities was measured in a rural area of Papua New Guinea where they share anopheline vectors. The spatial correlation of W. bancrofti was found to reduce by half over an estimated distance of 1.7 km, much smaller than the 50 km grid used by the World Health Organization rapid mapping method. For P. falciparum, negligible spatial correlation was found. After mass treatment with anti-filarial drugs, there was negligible correlation between the changes in the densities of the two parasites.
Geolocation and remote sensing technologies are increasingly being applied to the mapping and spatial analysis of infectious diseases, including lymphatic filariasis  and malaria . For such maps to reflect the real pattern of infection or disease, the sampling scale must be fine enough to register its variation. Lymphatic filariasis and malaria are both currently subject to renewed control programmes, and share anopheline vectors in some parts of the world. Filariasis is being targeted for elimination, but reduction in its infection intensity could conceivably, via the removal of infection-induced vector mortality, lead to more efficient transmission of malaria [3, 4]. Accurate mapping of the two infections can help monitor their control, including the detection of any unwanted interactions. Here, we present spatial analysis of both parasites in the same area of Papua New Guinea, using a recently developed technique which estimates the scale of variation, taking into account the typically highly skewed distribution of parasite counts .
In a rural area of the East Sepik Province, cross-sectional parasitological surveys were done of the population aged over five years, as part of a trial of diethylcarbamazine (DEC) plus ivermectin versus DEC alone against lymphatic filariasis . Although, with regard to malaria, the lower age threshold is a limitation, in the nearby Wosera area, peak prevalence was not reached till after five years . Wuchereria bancrofti microfilariae were counted microscopically on Nuclepore filters, through which 1 ml of night-collected blood had been passed. Asexual forms of Plasmodium falciparum were counted per 200 white cells from thick films. Anopheles mosquitoes, predominantly Anopheles punctulatus s.s., are the vectors of filariasis in the area .
Villages are divided into subunits, which we call hamlets, distinguished by names in the local languages. These were mapped with a hand-held Trimble Ensign Global Positioning System (GPS) machine. In the 1994 survey, immediately before the first round of treatment, blood samples were obtained from 2,219 people in 149 hamlets. Hamlets whose distance to their nearest neighbour was, due to GPS inaccuracy, measured as less than 10 m were combined with their nearest neighbour. There were two such instances, so reducing the total number of hamlets to 147.
As described elsewhere , spatial structure in the infection densities was fitted by a negative binomial model. The mean density is fitted as a log-linear function of age, sex, and the hamlet mean. In turn, the hamlet means are given a spatial structure, with closer hamlets being more highly correlated. Specifically, hamlet i adds a 'hamlet effect' u i to the logarithm of the mean density. When exponentiated, these are similar to standardized mortality ratios (SMRs) although relate to parasite densities rather than death rates, so can be called standardized parasite density ratios. These age- and sex-adjusted hamlet effect u il have a covariance matrix (1/φ)exp(-d il /α), where d il is the distance between hamlet i and l. The parameter α measures the scale of spatial correlation. More specifically, αloge2 is the distance over which the correlation reduces by half, which we call the 'half-distance'. The process can also be thought of a smoothing process of the raw hamlet means, adjusting for age and sex, and with a large spatial scale (α) corresponding to a greater degree of smoothing. The model was fitted by programs written in the C and FORTRAN languages, and checked in terms of ability to represent the spatial variability, and robustness to distributional assumptions.
One year after the first round of mass treatment with anti-filarial drugs (coverage 88%), all 14 randomized clusters showed a reduction in Williams mean W. bancrofti density, ranging from 47 to 96%, mean 80% . All 14 clusters showed an increase in Williams mean P. falciparum, the increases ranging from 2 to 255%, mean 91%. However, although DEC+ivermectin was more effective against W. bancrofti, no tendency was seen for those clusters to have larger increases in P. falciparum: in 3 of the 7 pairs the DEC+ivermectin cluster had the larger increase, with 4 having the opposite tendency . To further investigate the possibility that this increase was due to the treatment, the results of the spatial models were put into a hamlet-level analysis: an 'ecological' analysis in epidemiological terms. The difference in the post- and pre-treatment hamlet effects was calculated for both filariasis and malaria, and a linear relation tested using the Pearson correlation coefficient. There was negligible correlation between the changes in the densities of the two infections (r = 0.026, p = 0.75).
The World Health Organization methodology of RAGFIL (Rapid Assessment of the Geographical Distribution of Bancroftian Filariasis) is based on a 50 × 50 km grid . Gyapong et al. have compared this to a finer 25 km grid in Ghana, with both leading to operationally similar conclusions . When extended to three other countries of West Africa, the method showed filariasis to have an unexpectedly wide geographical distribution . The current results suggest that, at the opposite end of the spatial scale, it may be possible for foci to persist within the interstices of a 50 × 50 km grid, or even a 25 × 25 km one; a conclusion similar to that reached by Srividya et al. in India . However, the overall importance of any such effects will also depend on the relative magnitudes of variation at different scales, which we are unable to measure beyond the extent of our current study area. The need for accurate geographical monitoring may increase if the current elimination campaign reduces filariasis over a wide scale and proceeds from an 'attack' to a 'consolidation' phase (to borrow malaria eradication terminology).
Previous studies have found a very small scale of spatial variation of malaria [13, 14]. For example, Thompson et al., found the risk of malaria varying by a factor of 6 over 500 m . However, the scale of variation found in the current study was still smaller than expected, with even the closest hamlets showing little correlation. This is despite malaria and filariasis sharing the same anopheline vectors. The much smaller scale of variation of malaria is probably related to its rapid variation over time. By comparison, the risk of acquiring W. bancrofti infection per mosquito bite seems to be much smaller , but infections can last much longer. A single cross-sectional survey therefore reflects a longer-term aggregate of filariasis exposure history.
Our estimation of the hamlet effects can be seen as a smoothing method, in which the degree of smoothing is determined by the data. In this example, the west-east trend in filariasis density was not easily discernible in the raw data, while the low spatial correlation of malaria meant that the smoothed map is similar to the crude one. Any pattern presumably reflects spatial variation in factors which affect the parasite, vector or host, but which are not included in the model. In the current study, such factors are likely to include those related to vector density, such as distance to breeding sites, but this information is not available for the whole study area, and the lack of pattern in malaria infection suggests that other factors must also be involved. We view the technique as potentially useful in identifying such factors, as well as identifying 'hot spots' of infection. Our method takes account of the extreme skewness shown by most parasite distributions, which can cause problems even at the exploratory stage of analysis. For example, we found that using the common log(x+1) transformation to plot the hamlet 'geometric means' induced a spurious spatial correlation in malaria infection.
Excess mosquito mortality caused by W. bancrofti would make it feasible for filariasis control to enhance malaria transmission. An 'ecological' analysis did not show such an effect in this study, but monitoring should continue in those areas where Anopheles are vectors of both infections. For the elimination of filariasis, mapping at small as well as large scales will be necessary, including urban areas, especially as efforts proceed beyond the 'attack' phase.
List of abbreviations
Global Positioning System
Rapid Assessment of the Geographical Distribution of Bancroftian Filariasis
standardized mortality ratio
The drug trial was supported financially by the World Health Organization (grant number TDR 910466) and the government of Papua New Guinea. Julian Stander was partially supported by the European Union TMR network ERB-FMRX-CT96-0095. Bryan Grenfell was supported by the Wellcome Trust and the Biotechnology and Biological Sciences Research Council.
- Gyapong JO, Kyelem D, Kleinschmidt I, Agbo K, Ahouandogbo F, Gaba J, Owusu-Banahene G, Sanou S, Sodahlon YK, Biswas G, Kale OO, Molyneux DH, Roungou JB, Thomson MC, Remme J: The use of spatial analysis in mapping the distribution of bancroftian filariasis in four West African countries. Annals of Tropical Medicine and Parasitology. 2002, 96: 695-705. 10.1179/000349802125001735.View ArticlePubMedGoogle Scholar
- Snow RW, Craig MH, Deichmann U, le Sueur D: A preliminary continental risk map for malaria mortality among African children. Parasitology Today. 1999, 15: 99-104. 10.1016/S0169-4758(99)01395-2.View ArticlePubMedGoogle Scholar
- Alexander NDE: Heterogeneity and the Epidemiology of Lymphatic Filariasis. Cambridge: University of Cambridge; PhD thesis. 1998Google Scholar
- Pichon G, Leonard J, Gaillard FO, Pion S: Assessment of ELF program consequences on malaria transmission in West Africa. In: Vector-borne Diseases Control in the New Era: Science, Policy and Action; Vector Control Research Centre, Pondicherry, India. 2000, [http://www.bondy.ird.fr/~pichon/biblio/Pondich1.html]Google Scholar
- Alexander N, Moyeed R, Stander J: Spatial modelling of individual-level parasite counts using the negative binomial distribution. Biostatistics. 2000, 1: 453-463. 10.1093/biostatistics/1.4.453.View ArticlePubMedGoogle Scholar
- Bockarie MJ, Alexander NDE, Hyun P, Dimber Z, Bockarie F, Ibam E, Alpers MP, Kazura JW: Randomised community-based trial of annual single-dose diethylcarbamazine with or without ivermectin against Wuchereria bancrofti infection in human beings and mosquitoes. Lancet. 1998, 351: 162-168. 10.1016/S0140-6736(97)07081-5.View ArticlePubMedGoogle Scholar
- Genton B, Al-Yaman F, Beck H-P, Hii J, Mellor S, Narara A, Gibson N, Smith T, Alpers MP: The epidemiology of malaria in the Wosera area, East Sepik Province, Papua New Guinea, in preparation for vaccine trials. I. Malariometrc indices and immunity. Annals of Tropical Medicine and Parasitology. 1995, 89: 359-376.PubMedGoogle Scholar
- Bockarie M, Kazura J, Alexander N, Dagoro H, Bockarie F, Perry R, Alpers M: Transmission dynamics of Wuchereria bancrofti in East Sepik Province, Papua New Guinea. American Journal of Tropical Medicine and Hygiene. 1996, 54: 577-581.PubMedGoogle Scholar
- Alexander NDE, Grenfell BT: The effect of pregnancy on Wuchereria bancrofti microfilarial load in humans. Parasitology. 1999, 119: 151-156. 10.1017/S0031182099004527.View ArticlePubMedGoogle Scholar
- World Health Organization: Update on rapid assessment of Bancroftian filariasis. TDR News. 1999, 10-Google Scholar
- Gyapong JO, Remme JH: The use of grid sampling methodology for rapid assessment of the distribution of bancroftian filariasis. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2001, 95: 681-686.View ArticlePubMedGoogle Scholar
- Srividya A, Michael E, Palaniyandi M, Pani SP, Das PK: A geostatistical analysis of the geographic distribution of lymphatic filariasis prevalence in southern India. American Journal of Tropical Medicine and Hygiene. 2002, 67: 480-489.PubMedGoogle Scholar
- Cattani JA, Moir JS, Gibson FD, Ginny M, Paino J, Davidson W, Alpers MP: Small-area variations in the epidemiology of malaria in Madang Province. Papua New Guinea Medical Journal. 1986, 29: 11-17.PubMedGoogle Scholar
- Greenwood BM: The microepidemiology of malaria and its importance to malaria control. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1989, 83 (supplement): 25-29.View ArticlePubMedGoogle Scholar
- 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. American Journal of Tropical Medicine and Hygiene. 1997, 57: 550-559.PubMedGoogle Scholar
- Southgate BA: Intensity and efficiency of transmission and the development of microfilaraemia and disease: their relationship in lymphatic filariasis. Journal of Tropical Medicine and Hygiene. 1992, 95: 1-12.PubMedGoogle 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.