- Open Access
Spatio-temporal patterns of malaria infection in Bhutan: a country embarking on malaria elimination
© Wangdi et al; licensee BioMed Central Ltd. 2011
- Received: 3 November 2010
- Accepted: 16 April 2011
- Published: 16 April 2011
At the verge of elimination of malaria in Bhutan, this study was carried out to analyse the trend of malaria in the endemic districts of Bhutan and to identify malaria clusters at the sub-districts. The findings would aid in implementing the control activities. Poisson regression was performed to study the trend of malaria incidences at district level from 1994 to 2008. Spatial Empirical Bayesian smoothing was deployed to identify clusters of malaria at the sub-district level from 2004 to 2008.
Trend of the overall districts and most of the endemic districts have decreased except Pemagatshel, which has an increase in the trend. Spatial cluster-outlier analysis showed that malaria clusters were mostly concentrated in the central and eastern Bhutan in three districts of Dagana, Samdrup Jongkhar and Sarpang. The disease clusters were reported throughout the year. Clusters extended to the non-transmission areas in the eastern Bhutan.
There is significant decrease in the trend of malaria with the elimination at the sight. The decrease in the trend can be attributed to the success of the control and preventive measures. In order to realize the target of elimination of malaria, the control measure needs to be prioritized in these high-risk clusters of malaria.
- Malaria Transmission
- Malaria Case
- Indoor Residual Spray
- Malaria Incidence
Since the first malaria survey of Bhutan in 1962, malaria cases decreased gradually over the years. The first case of malaria was reported in 1965, there were 518 cases. The highest number of cases was recorded in 1994 with 39,852 cases and the highest death was recorded in 1994 with 62 deaths . The Annual Parasite Incidence (API) per 1,000 population at risk of malaria in Bhutan was 3.98% in the year 2006 . Bhutan contributed second lowest number of cases after Sri Lanka with 1,868 cases in SEARO region in the same year .
The control measures have also changed over the years. Indoor residual spraying (IRS) using dichlorodiphenyltrichloroethane (DDT) was implemented since 1964 with three rounds and later two rounds until 1994. With the reports of resistance to DDT in some parts of world, DDT was replaced by deltamethrin (synthetic pyrethoid) from 1994 till 1997 [3, 4]. Insecticide treated-bed nets (ITN) became the main control strategy from 1998 with focal IRS during outbreaks and emergencies, and high Plasmodium falciparum transmission areas with API >10%. Long-lasting insecticide-treated nets (LLINs) have proven to reduce the mortality and morbidity of malaria by reducing the man-mosquito contact , so Bhutan distributed over 100,000 LLINs with the grants from Global Fund to Fight against AIDS, Tuberculosis and Malaria (GFATM), as of 2006 . The diagnosis and treatment of malaria is integrated with the general health care system with early diagnosis and prompt treatment remaining the cornerstone of the malaria control. Microscopic facilities for early diagnosis are extended to all the health centres in these endemic districts. Microscopy is further supplemented with rapid diagnostic kits in most of the health facilities. Uncomplicated P. falciparum is treated with the combination of artemether and lumefantrine (Coartem®) with different regimen for the pregnant woman. The present regimen was initiated since 2006. While the treatment for Plasmodium vivax has not changed and is treated with chloroquine for three days and primaquine administered over 14-days [6, 7].
Disease surveillance connotes the ongoing systematic collection, analysis, and interpretation of out-come specific data for use in planning, implementation, and evaluation of public health practice . Its purpose is to detect changes in trends or distribution, so that necessary control and preventive measures can be initiated. Malaria surveillance includes laboratory confirmation of presumptive diagnosis, finding out the source of infection and identification of all cases and susceptible contacts and still others who are at risk, in order to prevent further spread of the disease. The ultimate objective of malaria surveillance is prevention and control of malaria in the community .
In Bhutan, malaria surveillance is a vertical system and through the passive reporting of the fever cases and the microscopically confirmed malaria cases. Therefore, every health centre; basic health unit II (BHU II), basic health unit I (BHU I), district hospital, regional and national referral hospital reports the malaria cases to Vector-borne Disease Control Programme (VDCP). Every week, the health centres in the endemic districts report the number of fever cases to the programme. At the end of every month, microscopically confirmed malaria cases are reported from the hospitals and health centres to VDCP. The fever surveillance serves as the early warning system so that prompt and early investigation can be initiated. In case of outbreaks, necessary interventions and control measures can be implemented.
As per the WHO guidelines, 10% of negative slides and 25% of positive cases are sent to the VDCP for quality control. As a national health policy, every fever patient that reports to health centres in the endemic districts and in the other non-endemic districts with a history of travel to malaria endemic areas in the past two weeks are needed to be tested for malaria either with the microscopy or rapid diagnostic test (RDT). However, the positive cases with RDT are further confirmed microscopically.
At the programme level, the reported fever and malaria cases are analysed such as the recent trends and information is fed back to the health centres. In case of unusually large numbers of fever and malaria cases, the health centres which report the unusual episodes are facilitated to carry out the investigations for possible causes of unusual fever or malaria cases for early containment and control.
Health facilities under the study in endemic districts.
Types of health centres
Overall malaria occurrence
The data for the study was extracted from the database of VDCP programme, Department of Public Health, Ministry of Health. All microscopically-confirmed malaria cases observed in any BHU IIs, BHU Is and hospitals are regularly reported to the VDCP. In this study, monthly malaria occurrences during 1994 to 2008 at the district level and 2004 to 2008 at the sub-district level were obtained from the database for the analysis.
Trend analysis in this study was assessed using Poisson regression. Poisson regression is a type of a Generalised Linear Model (GLM) model that follows a Poisson distribution. This regression model is used to model discrete events that occur infrequently in time or space; hence, it is sometimes called the distribution of rare events . In this study, the distribution of microscopically confirmed malaria cases (events) occurred under baseline population across years and districts were assumed to follow the Poisson distribution. Then the Poisson regression model was performed to calculate the incidence rate ratio and 95% confidence interval was calculated using the Poisson regression model to determine the change in malaria incidence over years, stratified by districts.
The sub-districts analysis of the malaria in the endemic districts and mapping of the malaria clusters was done using the five-year malaria surveillance data from 2004 to 2008. Spatial Empirical Bayesian (SEB) was used to minimize the phenomenon of the Modifiable Areal Unit Problem (MAUP). The clusters were observed at the sub-district level at the monthly period over five years.
The Moran's I statistics was used to calculate the autocorrelation or global clusters of malaria incidence rates of different sub-districts.
Where W ij = is a measure of the closeness of areas i and j; W ij = 1 if areas i and j are adjacent (have a common boundary) and 0 otherwise
Z i = could be the residuals (O i - E i ) or standard mortality or morbidity ratio (SMR) of an area/attribute value for area i
Z j = attribute value for area j
= mean value in study region
n = the numbers of areas
I = Moran's I
The spatial weights indicate whether regions neighbour each other. Spatial weights based on contiguity matrix in this research are queen contiguity.
δ = standard deviation of the variable X
Z i = Z-score of X i
W ij = spatial weights matrix on row-standardized form
LISA detects clusters of either similar or dissimilar disease frequency values around observations. The sum of the LISA for all observations is proportional to the global Moran's I statistic [12, 14].
Malaria cases of the endemic districts of Bhutan from 2004 to 2008.
Different age break down of the malaria cases from 2004 to 2007.
Incidence rate ratio, standard error and 95% confidence interval (CI) of seven endemic districts and the overall endemic districts.
All endemic districts
Overall malaria occurrence
Malaria in Bhutan is mostly reported by seven districts namely: Chukha, Dagana, Pemagatshel, Samdrup Jongkhar, Samtse, Sarpang and Zhemgang (Table 2 and Figure 1). Most of the transmission usually takes place in the summer months when ambient temperature is very favourable. These districts have sub-tropical climate with hot summers, temperatures reaching as high as 30°Celsius; most of the districts have abundant rainfall with humidity above 60% throughout the year . The transmission within the districts vary greatly for a number of reasons, such as varying in altitude, climatic factors and environmental factors of sub-districts that make up districts.
There appears to have a seasonal pattern with two peaks, first smaller peak just before the onset of the summer and the bigger peak during the summer (June-August) when rainy season starts. The small peak could be due to the increased in temperature just before the summer, which enhance the vector multiplication and development. Majority of rural population of Bhutan depends on agriculture for their livelihood. Therefore, people would be exposed to the mosquito bite during this time while they are involved in agricultural activities. The second bigger peak would be associated with the summer months, where rainfall is abundant, this providing the aquatic medium for the vector proliferation. With the onset of summer, irrigation in the paddy fields is usually intensified. The paddy fields and the irrigation channels had been identified as the source of the mosquitoes. Similar seasonal transmission patterns with two peaks had been reported in neighbouring Indian states on Assam and West Bengal [16–18].
The malaria trend decreased gradually with upsurge again in 1999 till 2003. This was closely related to the increased malaria cases in Assam and West Bengal [16–18]. Bhutan shares eastern half of borders with Assam, three districts that shares the borders are Samdrup Jongkhar, Pemagatshel and Sarpang. While western districts bordering West Bengal are Samtse, Chukha and Dagana. The malaria cases from Samdrup Jongkhar and Sarpang increased during this period. However, no significant increase in malaria cases in the districts bordering West Bengal was noted, despite open border between these districts and West Bengal, India.
The only district to have an increase in trend was Pemagatshel district (Table 4). Three endemic sub-districts of Dechenling, Norbugang and Chorkroling were added to Pemagatshel from Samdrup Jongkhar district in 2006. This could have lead to the increased malaria cases reported from Pemagatshel district. Other plausible reason is the migration of the people for job opportunity since the construction of one of the largest cement plant in Bhutan is under way. This job opportunity could lead to the migration of people both internal as well as expatriates. The risk of increased malaria cases should be anticipated in the coming years.
The identification of clusters in the sub-districts would provide important information on the malaria transmission especially with very few cases as in Bhutan. It had been reported in other studies that malaria transmission usually does not seems to take place at an altitude above 2000 meters [20, 21]. However in this study, malaria clusters were persistently observed in sub-districts which are located at an altitude above 2000 meters, like Serthi and Shingkhar Lauri under Samdrup Jongkhar district for the month of May and June over the study period (Figure 2). Malaria clusters in these high altitude districts could be as a result of imported cases from other endemic areas, because most of the businesses and administrations are located in sub-districts, where malaria transmission takes place. People could possibly get the infection from the sub-districts with malaria transmission while they travel for work opportunities and businesses, then manifest clinical signs and symptoms on returning to their villages. This finding could be important to the programme in mitigating the malaria transmission among the internal migrants.
In this study, malaria clusters occurred throughout the year. A number of plausible reasons were hypothesized; firstly, these areas are located in the environmentally favourable region with sub-tropical climate, this facilitates in vector multiplication. Secondly, most of these sub-districts have large areas as the paddy field/wet lands. These paddy fields and irrigation channels could provide breeding sites for the mosquitoes. The main vectors in Bhutan are Anopheles pseudowillmori (33%), Anopheles vagus (31.5%) and Anopheles peditaeniatus (10.3%) . Similar finding of Anopheles maculates multiplying in the paddy field had been report by Kankaew . The introduction of community fish farming in these districts in recent years could also be one of the factors for identification of malaria clusters. Common fish that are reared in these fisheries are Cirrhina cirrhosus (Mirgal), Ctenopharyngodon idella (Grass carp), Hypopthalmicthys molitrix (Silver carp), Cyprinus carpio (Common carp), Aristichthys nobilis (Big head carp), Catla catla (Catla), and Labeo rohita (Rohu) (personal communications). These fish farming provided aquatic environment for the development of mosquitoes. Fourthly, the introduction of non-immune population through migrations for various purposes could be one of the reasons. Cross-border malaria could be other reasons since borders are quite open and people of either Bhutan or India can easily cross the border. Malaria transmission on either side of border would have a direct effect on malaria cases on either side of border, since malaria transmission is spatially related. As seen in the malaria cluster map, clusters mostly occurred in the sub-districts that share the border with India (Figure 2). Two states that share the border with Bhutan in the east and west are Assam and West Bengal, respectively. Dagana, Chukha and Samtse share border with West Bengal, while Sarpang, Zhemgang, Pemagatshel and Samdrup Jongkhar share border with Assam. These districts particularly Sarpang and Samdrup Jongkhar showed malaria clusters consistently throughout the year. Assam is malaria endemic district and contributes more than 5% of total malaria in India . West Bengal is also considered as malaria endemic state , but only some of the sub-districts of Dagana showed clusters while other sub-districts of Chukha and Samtse districts did not have malaria clusters for some unknown reasons.
There is gradual decrease in the malaria trend in Bhutan with very low malaria incidence. The decrease in the trend could be attributed to the success of control and preventive measures. Malaria clusters vary amongst the districts of the endemic areas. The finding of malaria clusters in the areas where malaria transmissions usually do not occur suggests that the control measures in these regions need to be addressed. The elimination of malaria in Bhutan could be possible with the control measures reinforced in the areas of malaria clusters at the sub-districts.
We would like to thank Ministry of Health, Royal Government of Bhutan for allowing us to use the surveillance data. We also would like to thank the staff at the VDCP for helping to extract the data and providing all the necessary information during the process of analysis and writing this papers.
- Vector-borne Disease Control Programme: Annual Report. 2007Google Scholar
- WHO: Malaria Disease Burden in SEA Region. [http://www.searo.who.int/EN/Section10/Section21/Section340_4018.htm]
- World Health Organization: Sixth report of expert committee on malaria. Technical Reports Series. 1957, 123: 184-Google Scholar
- Roberts DR, Andre RG: Insecticide resistance issues in vector-borne disease control. Am J Trop Med Hyg. 1994, 50 (6 Suppl): 21-34.PubMedGoogle Scholar
- Lengeler C: Insecticide treated bed nets and curtains for preventing malaria. Cochrane review. 2004Google Scholar
- Vector-borne Disease Control Programme: Guideline on the management of malaria and vector control in Bhutan. 2008Google Scholar
- Vector-borne Disease Control Programme: Bhutan Indicator Framework 2008-2013. 2007Google Scholar
- Lombardo JS, Buckeridge DL: Disease surveillance a public health informatics approach. 2007, A John Wiley & Sons Inc, PublicationGoogle Scholar
- Pull JH: Malaria surveillance methods, their uses and limitations. Am J Trop Med Hyg. 1972, 21: 651-657.PubMedGoogle Scholar
- Bhutan health management information system, Ministry of Health, Royal Goverment of Bhutan: Annual Health Bulletin. 2010Google Scholar
- Pagango M, Gauvreau K: Principles of Biostatistics. 2000, 2Google Scholar
- Pfeiffer D, Robinson T, Stevenson M, Stevens K, Rogers D, Clements A: Spatial analysis in epidemiology. 2008, Oxford University PressView ArticleGoogle Scholar
- Elliott P, Wakefield JC, Best NG, Briggs DJ: Spatial epidemiology methods and applications. 2000, Oxford University PressGoogle Scholar
- Wong DWS, Lee J: Statistical analysis of geographic information with ArcView GIS and Arc GIS. 2005, John Wiley and Sons, IncGoogle Scholar
- Wangdi K, Singhasivanon P, Silawan T, Lawpoolsri S, White NJ, Kaewkungwal J: Development of temporal modelling for forecasting and prediction of malaria infection using time-series and ARIMAX: A case study in endemic districts of Bhutan. Malar J. 2010, 9: 251-10.1186/1475-2875-9-251.PubMed CentralView ArticlePubMedGoogle Scholar
- Dev V, Phookan S, Sharma VP, Anand SP: Physiographic and entamologic risk factors of malaria in Assam, India. Am J Trop Med Hyg. 2004, 71: 451-456.PubMedGoogle Scholar
- Dev V, Phookan S, Sharma VP, Dash AP, Anand SP: Malaria parasite burden and treatment seeking behavior in ethnic communities of Assam, Northeastern India. J Infect. 2006, 52: 131-139. 10.1016/j.jinf.2005.02.033.View ArticlePubMedGoogle Scholar
- Sharma PK, Ramakrishnan R, Hutin YJF, Gupte MD: Increasing incidence of malaria in Kurseong, Darjeeling District, West Bengal, India, 2000-2004. Trans R - Soc Trop Med Hyg J. 2009, 103: 691-697. 10.1016/j.trstmh.2008.07.019.View ArticleGoogle Scholar
- Zhou G, Sirichaisinthop J, Sattabongkot J, Jones J, Bjornstad ON, Yan G, Cui L: Spatio-temporal distribution of P. falciparum and P. vivax malaria in Thailand. Am J Trop Med Hyg. 2005, 72: 256-262.PubMedGoogle Scholar
- Brooker S, Leslie T, Kolaczinski K, Mohsen E, Mehboob N, Saleheen S, Khudonazarov J, Freeman T, Clements A, Rowland M, Kolaczinski J: Spatial epidemiology of P. vivax, Afghanistan. Emerg Infect Dis. 2006, 12: 1600-1602.PubMed CentralView ArticlePubMedGoogle Scholar
- Ernest KC, Adoka SO, Kowuor DO, Wilson ML, John CC: Malaria hotspot area in a highland Kenya site are consistent in epidemic and non-epidemic years and are assoicated with ecological factors. Malar J. 2006, 5: 78-10.1186/1475-2875-5-78.View ArticleGoogle Scholar
- Kankaew P, Krasaisub S, Sithiprasana R: Soil analysis around anophiline breeding habitats in north-western Thailand. Southeast Asian J Trop Med Public Health. 2005, 36: 1153-1161.PubMedGoogle Scholar
- Kumar A, Valecha N, Jain T, Dash AP: Burden of malaria in India: retrospective and prospective view. Am J Trop Med Hyg. 2007, 77: 69-78.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.