Space–Time Clustering Characteristics of Malaria in Bhutan at the End Stages of Elimination
by
, , , , and
Kinley Wangdi
1,*
,
Kinley Penjor
2
,
Tobgyal
2,
Saranath Lawpoolsri
3
,
Ric N. Price
4,5,6,
Peter W. Gething
7,8,
Darren J. Gray
1,
Elivelton Da Silva Fonseca
9
and
Archie C. A. Clements
7,8 
1
Research School of Population Health, Australian National University, Canberra, ACT 2601, Australia
2
Vector-Borne Diseases Control Program, Department of Public Health, Ministry of Health, Gelephu 31101, Bhutan
3
Department of Tropical Hygiene, Faculty of Tropical Medicine, Mahidol University, Bangkok 10400, Thailand
4
Global and Tropical Health Division, Menzies School of Health Research, Charles Darwin University, Darwin, NT 0810, Australia
5
Centre for Tropical Medicine and Global Health, Nuffield Department of Clinical Medicine, University of Oxford, Oxford OX1 2JD, UK
6
Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok 10400, Thailand
7
Faculty of Health Sciences, Curtin University, Bentley, WA 6102, Australia
8
Telethon Kids Institute, Nedlands, WA 6009, Australia
9
Institute of Geography, Federal University of Uberlândia, Uberlândia, MG 38408-100, Brazil
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(11), 5553; https://doi.org/10.3390/ijerph18115553
Submission received: 14 April 2021
/
Revised: 16 May 2021
/
Accepted: 17 May 2021
/
Published: 22 May 2021
(This article belongs to the Special Issue Spatial Epidemiology of Emerging Infectious Diseases)
Abstract
:Malaria in Bhutan has fallen significantly over the last decade. As Bhutan attempts to eliminate malaria in 2022, this study aimed to characterize the space–time clustering of malaria from 2010 to 2019. Malaria data were obtained from the Bhutan Vector-Borne Disease Control Program data repository. Spatial and space–time cluster analyses of Plasmodium falciparum and Plasmodium vivax cases were conducted at the sub-district level from 2010 to 2019 using Kulldorff’s space–time scan statistic. A total of 768 confirmed malaria cases, including 454 (59%) P. vivax cases, were reported in Bhutan during the study period. Significant temporal clusters of cases caused by both species were identified between April and September. The most likely spatial clusters were detected in the central part of Bhutan throughout the study period. The most likely space–time cluster was in Sarpang District and neighboring districts between January 2010 to June 2012 for cases of infection with both species. The most likely cluster for P. falciparum infection had a radius of 50.4 km and included 26 sub-districts with a relative risk (RR) of 32.7. The most likely cluster for P. vivax infection had a radius of 33.6 km with 11 sub-districts and RR of 27.7. Three secondary space–time clusters were detected in other parts of Bhutan. Spatial and space–time cluster analysis identified high-risk areas and periods for both P. vivax and P. falciparum malaria. Both malaria types showed significant spatial and spatiotemporal variations. Operational research to understand the drivers of residual transmission in hotspot sub-districts will help to overcome the final challenges of malaria elimination in Bhutan.
1. Introduction
Historically, malaria in Bhutan has been reported from seven out of twenty districts: Chukha, Dagana, Pemagatshel, Samdrup Jongkhar, Samtse, Sarpang and Zhemgang [1]. These districts lie in the southern foothills of the Himalayas bordering the Indian states of Assam and West Bengal. The main malaria parasite species are Plasmodium falciparum and Plasmodium vivax, and the main suspected vector species are Anopheles pseudowillmori and Anopheles culicifacies [2]. Malaria cases have been falling in recent years. As a result, Bhutan announced a national strategy to eliminate malaria by 2022 [3,4]. The success of malaria control in Bhutan can be attributed to the implementation of a three-pronged approach, namely: (1) universal access to early diagnosis and prompt treatment with artemisinin-based combination therapy (ACT); (2) the protection of at-risk populations with indoor residual spraying (IRS) and long-lasting insecticide nets (LLINs), and (3) integrated vector management (IVM).
As the malaria map shrinks, cases tend to be confined to hard-to-reach populations, often in border areas [5,6]. Identifying these last pockets of transmission is imperative for the program to eliminate the parasite and prevent its reintroduction. Spatial epidemiological tools (including geographical information systems (GIS) and spatial analytic methods) can be used to quantify spatial and spatiotemporal patterns of malaria, including clustering [7,8,9,10]. The identification of malaria clusters can help in the delineation of problem areas and the application of targeted program interventions suitable in the elimination phase [11]. Focused interventions in areas with a higher risk of malaria are likely to be more cost-effective than uniform resource allocation, particularly in resource-constrained settings for sustainable elimination programs [9,12,13]. In addition, knowledge of seasonal patterns will aid in the estimation of time for malaria transmission for initiating suitable and appropriate control measures [14,15].
The aims of this study were to quantify the temporal, spatial and spatial–temporal patterns of both species of malaria in sub-districts in Bhutan using retrospective surveillance data.
2. Methods
2.1. Study Site
Bhutan is a small Himalayan country in South Asia, neighboring India in the east, south and west, and China in the north. The country is divided into 20 districts and 205 sub-districts. The total population of Bhutan in 2017 was estimated at 681,720 [16]. There is a wide variation in the climate, with sub-tropical conditions in the southern foothills neighboring India to snow-capped mountains rising above 4500 m from sea level in the north. Malaria transmission occurs mainly in southern Bhutan in 82 sub-districts within the seven historically malaria-endemic districts (Figure 1).
2.2. Data Source
Malaria data for the study were obtained from the national malaria surveillance system, housed within the Bhutan Vector-borne Diseases Control Program (VDCP). Monthly malaria cases for 2010–2019 from 82 sub-districts were extracted, including those with zero cases. In these malaria-transmitting sub-districts, all febrile patients presenting to health facilities (hospitals and primary health centers) undergo a mandatory blood test with RDTs or microscopy. Malaria cases stratified by species as P. falciparum and P. vivax are reported through the national malaria surveillance system at the end of each month. All patients in public health facilities are offered treatment free of charge according to national guidelines. The yearly sub-district population was calculated using the district population growth rate (Chukha 1.7, Dagana 1.9, Pemagatshel 1.2, Samdrup Jongkhar 1.8, Samtse 1.6, Sarpang 2.0 and Zhemgang 1.4) as outlined by the National Statistical Bureau and the Office of the Census Commissioner of Bhutan [17,18]. An electronic map of sub-district boundaries in shapefile format was obtained from the Global Administrative Areas database (http://www.gadm.org/country, accessed on 21 December 2020).
2.3. Exploration of Seasonal Patterns and Inter-Annual Patterns
The malaria time series (January 2010–December 2019) was decomposed using seasonal trend decomposition based on locally (STL) weighted regression. The STL model was structured as follows:
where Yt represents numbers of local malaria cases with logarithmic transformation, St is the additive seasonal component, Tt is the trend, and Rt is the “remainder component”; t is time in months [19,20,21].
2.4. Spatial and Spatiotemporal Cluster Analysis
Kulldorff’s spatial and space–time scan statistics were applied to identify malaria incidence clusters by species (P. falciparum and P. vivax) at a sub-district level over the 10-year study period. The analysis was implemented using the SaTScanTM software v.9.4.2 (M Kulldorf and Information Management Services Inc., Calverton, MD, USA) [22]. Three data sets were inputted to the program for analysis: the shapefile centroids of the sub-districts, using the latitude and longitude in the World Geodetic Coordinate System of 1984, yearly populations and reported cases of 82 sub-districts.
Spatial cluster analysis for each year were identified using a Poisson variant of the spatial scan statistic. The following configurations were used in the analysis: clusters of maximum size equal to 25% of the exposed population, no geographical overlapping of clusters and 999 Monte Carlo replications for the testing of statistical significance to ensure adequate power for defining clusters [23,24]. In the space–time cluster analysis, the unit of time was a month, with same configurations as the spatial analysis. SaTScan uses moving scanning windows (circular or elliptical) of varying sizes to estimate the probability that the frequency of positive cases within a window is greater than what is expected by chance. It takes into account the observed number of cases inside and outside the windows, to detect clusters and estimate relative risk (RR). For each location, the window size with the highest log-likelihood ratio (LLR) was considered the most probable cluster, i.e., the cluster that is least likely to have occurred randomly. The SaTScan output for statistically significant clusters includes the location of the center of the scanning window, the radius of the scanning window, the number of observed and expected positives within the circle and the relative risk, LLR and p value of the cluster. Statistically significant clusters were defined as those with a p < 0.05. Centroids of the districts and significant clusters were mapped using the GIS ArcMap 10.5.1 (ESRI, Redlands, CA, USA).
3. Results
3.1. Characteristics of the Study Population
Between 2010 and 2019, a total of 768 confirmed malaria cases were reported, including 314 (41%) patients infected with P. falciparum and 454 (59%) with P. vivax. Forty-three percent (330) were in the age group of 21–45 years and only 1.3% (10) were children under one year of age. Males made up 65.1% (500) of patients and farmers and students constituted 47.3% (363) and 25.0% (192), respectively. The years with the highest and lowest numbers of cases were 2010 and 2013, with 52.2% (401) and 1.6% (12), respectively. Sarpang, Samdrup Jongkhar and Samtse reported 68.1% (523), 10.3% (79) and 6.5% (65) of all cases (Table 1). The yearly annual parasite incidence (API) of sub-districts (cases of sub-district divided by population per 1000) is shown in Supplementary Figures S1 and S2.
3.2. Time-Series Decompositions
The time-series decompositions show a clear seasonal pattern which is evident in the raw data for both P. falciparum and P. vivax. A large peak occurred for P. falciparum in May, followed by two smaller peaks later in the year. For P. vivax, a large peak occurred in May, followed by a plateauing for the subsequent three months and dropping off during the remainder of the rainy season. The inter-annual pattern showed a large peak in 2010, with a lower incidence in subsequent years for both species of malaria (Figure 2).
3.3. Purely Spatial Analysis
The most likely clusters of P. falciparum infections were in the central part of southern Bhutan for eight of the ten years of the study period (Table 2). In 2014 and 2019, the most likely clusters shifted to the eastern part of Bhutan, in the sub-district of Langchenphug. There were three secondary clusters in 2010, 2012 and 2014, respectively (Figure 3). No cluster was observed in 2016 or 2017.
Similar to P. falciparum, the most likely clusters of P. vivax infection were in the central part of Bhutan, except in 2017, when it was in the western part of Bhutan (Table 3). There was a secondary cluster in 2010 in the eastern part of Bhutan (Figure 4). The largest radius cluster was observed in 2010 for both species.
3.4. Spatiotemporal Analysis
The most likely spatiotemporal cluster was in the central part of the country between January 2010 to June 2012 for both species (Figure 5). The most likely cluster of P. falciparum infection had a radius of 50.4 km and included 26 sub-districts with an RR of 32.7. Two secondary clusters in April 2014 had a radius of 0 (indicating a single sub-district; RR = 25.4) and 56.8 km (RR = 16.7), respectively. Another secondary cluster was detected from January 2010 to December 2014 in the far east of Bhutan (one district; RR 25.4).
The location and period of the most likely cluster for P. vivax was similar to P. falciparum but it had a smaller radius of 33.6 km, contained 11 sub-districts and had an RR of 27.7. Two secondary clusters were observed from March–May 2010 with a radius of 9.5 km (two sub-districts; RR = 104.7) and 13.1 km (seven sub-districts; RR = 17.9). The third secondary cluster was observed in August–November 2017 with one sub-district and an RR of 21.3 (Table 4 and Table 5).
4. Discussion
This 10-year retrospective analysis of malaria surveillance data from Bhutan provides an in-depth characterization of significant spatial and spatiotemporal clustering of both P. falciparum and P. vivax malaria incidence. Both species showed a clear seasonal pattern with a large peak in the early part of the year. Malaria clusters of both species were observed in the central parts of Bhutan along the international border throughout the study period. The most likely spatiotemporal clusters were identified from January 2010–June 2012 for both species of malaria in Sarpang District. There were three secondary clusters for both species of malaria in other districts.
The seasonal pattern corresponds with the monsoon in Bhutan, which lasts for five months from May to September [25]. The monsoon season supports the breeding and growth of the main malaria vectors in Bhutan; Anopheles pseudowillmori and Anopheles Culicifacies [2,26]. Therefore, VDCP can utilize this knowledge to implement focused interventions of malaria including the distribution of LLINs and IRS. Even though LLIN use was reported to be high in Bhutan [27], monitoring the use of LLINs should be strengthened, especially during these months, to interrupt local transmission.
The spatial cluster analysis identified high-risk sub-districts along the international border with India. As Bhutan approaches malaria elimination, a key challenge will be residual malaria in the population living in the border area, highlighting the need for novel approaches to address the problem. Foreign nationals entering Bhutan through the land crossings in the southern districts include two groups: those staying in Bhutan on a long-term basis (such as Indians working on developmental projects) and those visiting Bhutan during day visits for business and employment [2]. Between 2006–2014, daily visitors crossing the southern border with India for employment accounted for 13% of all malaria cases [2]. In light of ongoing transmission in Sarpang District, the VDCP will need to monitor the incidence of malaria in this migrant population closely, using both passive and active case detection to better understand local transmission dynamics. Further, active blood surveys in the hotspot sub-districts, followed up by targeted interventions, should be started in earnest to disrupt local transmission.
The lack of spatial clusters of P. falciparum in 2016 and 2017 could be attributed to the effects of LLINs which are distributed every three years, and which might have disrupted the natural transmission dynamics of malaria in that year. Since the start of the LLIN distribution program in 2006, five rounds of LLINs were completed in 2006, 2010, 2014, 2017 and 2020. In 2014, LLIN coverage was 99.0%, and 94% of people slept regularly under LLINs [27].
The space–time analysis identified a primary cluster in Sarpang District from January 2010–June 2012. Sarpang continues to be the location of the remaining hotspot of malaria in Bhutan [1,25]. The elimination of malaria from Sarpang has proven more difficult as compared to other districts. A major reason is likely to be the proximity of villages in Sarpang to the very long and porous border with the Indian states of Assam [2,28]. The risk of transmission of malaria across the border within the flight range of infected mosquitoes in these villages is high. Control measures across the international border (in Assam) are sporadically implemented and transmission is maintained due to the local populations living in villages near to forest fringe areas [29]. Cross-border movement has been acknowledged as one of the key challenges in malaria elimination [5,13,28], leading to a series of cross-border meetings being held to strengthen collaboration between India and Bhutan to achieve malaria elimination [30]. However, the presence of a P. falciparum cluster in Langchenphu (in Samdrup Jongkhar District) should precipitate a call for the VDCP to invest time and resources to study the local drivers of malaria.
A secondary space–time cluster of P. vivax from August to November 2017 was observed in Samtse District. This was due to an outbreak of P. vivax in Samtse with 16 cases reported over three months. A WHO external review of malaria surveillance reported that the outbreak was due to the late diagnosis of indigenous cases with a consequent delay in response. As malaria cases dwindle in Bhutan, people are less likely to seek care for malaria as well as use protective measures regularly. Therefore, the prompt diagnosis of malaria and adherence to the regular use of LLINs should be reinforced through regular health education. Further, primaquine is used for the radical cure of P. vivax in Bhutan [31], although no studies on adherence to the recommended 14-day radical cure regimen have been undertaken. This will be an important step as the country attempts to eliminate malaria by 2022. Poor adherence has been reported to be an important factor that drives P. vivax transmission in other parts of the world [32]. It is timely for the VDCP to undertake operational research on the barriers and enablers of the current radical cure regimen, including directly observed treatment (DOT) recommended for P. vivax in Bhutan, to develop strategies to improve adherence.
Our study has a number of important limitations. The main limitations are inherent to the method used for the space–time analysis, especially the scale (sub-district) at which data were available. Even though SaTScan is a powerful statistical tool frequently used to analyze spatial and spatio-temporal patterns of disease, it is a challenge to determine an optimal set of scaling parameters for the SaTScan analysis, especially for the maximum window size. Indeed, large maximum window sizes (such as 50% of the total population) can hide small, homogeneous clusters within the larger and more heterogeneous ones, while on the other hand, small maximum window sizes can miss significant regional-level clusters. Second, the completeness and representativeness of surveillance data cannot be ascertained. However, this is likely to be minimal since all malaria cases are treated only in public health facilities in Bhutan (there is no private health care sector) and they report routinely to the national malaria surveillance system. Third, sub-district populations were not available, and they were projected based on the National Housing and Population Census 2007. However, the main strength of our study is the breakdown of the monthly malaria incidence, stratified by species at a fine-resolution (sub-districts), providing a more precise space–time analysis.
In conclusion, malaria has been eliminated in most sub-districts in Bhutan and a few remaining clusters were largely confined to the southern district of Sarpang. The VDCP of Bhutan continues to plan for the sub-national elimination of malaria while awaiting national elimination. Operational research to understand the burden of malaria in daily migrant workers and transmission dynamics in hotspot sub-districts could be valuable to overcome the last mile challenge of malaria elimination in Bhutan.
Supplementary Materials
The following are available online at https://www.mdpi.com/article/10.3390/ijerph18115553/s1, Figure S1: Annual parasite incidence of Plasmodium falciparum, Figure S2: Annual parasite incidence of Plasmodium vivax.
Author Contributions
Conceptualization, K.W. and K.P.; methodology, K.W.; software, K.W.; validation, K.P.; formal analysis, K.W.; data curation, K.P.; writing—original draft preparation, K.W.; writing—review and editing, K.P., T., S.L., R.N.P., P.W.G., E.D.S.F., D.J.G. and A.C.A.C.; visualization, K.W.; supervision, A.C.A.C.; project administration, K.P. and T. All authors have read and agreed to the published version of the manuscript.
Funding
No funding was required for this study. R.N.P. is a Wellcome Trust Senior Fellow in Clinical Science (200909).
Institutional Review Board Statement
Data use was approved by the Ministry of Health, Bhutan (CI/(09)/2020/16/24 and 24/11/2020). Ethical approval was waived for this study because we used secondary data without individual patient information.
Informed Consent Statement
Not applicable.
Data Availability Statement
Data supporting the results can be obtained upon request from Ministry of Health, Bhutan.
Acknowledgments
The authors would like to thank the Ministry of Health, Bhutan, for providing access to the database of malaria case notification reports between 2010–2019.
Conflicts of Interest
The authors declare that they have no competing interests.
References
- Wangdi, K.; Kaewkungwal, J.; Singhasivanon, P.; Silawan, T.; Lawpoolsri, S.; White, N.J. Spatio-temporal patterns of malaria infection in Bhutan: A country embarking on malaria elimination. Malar. J. 2011, 10, 89. [Google Scholar] [CrossRef] [Green Version]
- Wangdi, K.; Banwell, C.; Gatton, M.L.; Kelly, G.C.; Namgay, R.; Clements, A.C. Malaria burden and costs of intensified control in Bhutan, 2006–2014: An observational study and situation analysis. Lancet Glob. Health 2016, 4, e336–e343. [Google Scholar] [CrossRef] [Green Version]
- Asia Pacific Malaria Elimination Networks. Bhutan. Available online: https://apmen.org/country/Bhutan/ (accessed on 20 May 2021).
- Vector-Borne Disease Program. Strategic Plan for Elimination of Malaria and Prevention of Re-Introcution in Bhutan 2020–2025; Vector-borne Disease Program, Department of Public Health, Ministry of Health, Royal Government of Bhutan: Gelephu, Bhutan, 2020.
- Wangdi, K.; Gatton, M.L.; Kelly, G.C.; Clements, A.C. Cross-border malaria: A major obstacle for malaria elimination. Adv. Parasitol. 2015, 89, 79–107. [Google Scholar] [PubMed] [Green Version]
- Wangdi, K.; Clements, A.C.A. Ending Malaria Transmission in the Asia Pacific Malaria Elimination Network (APMEN) Countries: Challenges and the Way Forward. In Towards Malaria Elimination: A Leap Forward; IntechOpen: London, UK, 2018. [Google Scholar]
- Bi, Y.; Hu, W.; Yang, H.; Zhou, X.N.; Yu, W.; Guo, Y.; Tong, S. Spatial patterns of malaria reported deaths in Yunnan Province, China. Am. J. Trop. Med. Hyg. 2013, 88, 526–535. [Google Scholar] [CrossRef] [Green Version]
- Soto-Calle, V.; Rosas-Aguirre, A.; Llanos-Cuentas, A.; Abatih, E.; DeDeken, R.; Rodriguez, H.; Rosanas-Urgell, A.; Gamboa, D.; D’Alessandro, U.; Erhart, A.; et al. Spatio-temporal analysis of malaria incidence in the Peruvian Amazon Region between 2002 and 2013. Sci. Rep. 2017, 7, 40350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Coleman, M.; Coleman, M.; Mabuza, A.M.; Kok, G.; Coetzee, M.; Durrheim, D.N. Using the SaTScan method to detect local malaria clusters for guiding malaria control programmes. Malar. J. 2009, 8, 68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Coulibaly, D.; Rebaudet, S.; Travassos, M.; Tolo, Y.; Laurens, M.; Kone, A.K.; Traore, K.; Guindo, A.; Diarra, I.; Niangaly, A.; et al. Spatio-temporal analysis of malaria within a transmission season in Bandiagara, Mali. Malar. J. 2013, 12, 82. [Google Scholar] [CrossRef] [Green Version]
- Sweeney, A.W. The Application of GIS in Malaria Control Programs. In Proceedings of the 10th Colloquium of the Spatial Information Research Centre, University of Otago, Dunedin, New Zealand, 16–19 November 1998. [Google Scholar]
- Clements, A.C.A.; Reid, H.L.; Kelly, G.C.; Hay, S.I. Further shrinking the malaria map: How can geospatial science help to achieve malaria elimination? Lancet Infect. Dis. 2013, 13, 709–718. [Google Scholar] [CrossRef]
- Wangdi, K.; Pasaribu, A.P.; Clements, A.C.A. Addressing hard-to-reach populations for achieving malaria elimination in the Asia Pacific Malaria Elimination Network countries. Asia Pac. Policy Stud. 2021, 1–13. [Google Scholar] [CrossRef]
- Stuckey, E.M.; Smith, T.; Chitnis, N. Seasonally dependent relationships between indicators of malaria transmission and disease provided by mathematical model simulations. PLoS Comput. Biol. 2014, 10, e1003812. [Google Scholar] [CrossRef] [Green Version]
- Nguyen, M.; Howes, R.E.; Lucas, T.C.D.; Battle, K.E.; Cameron, E.; Gibson, H.S.; Rozier, J.; Keddie, S.; Collins, E.; Arambepola, R.; et al. Mapping malaria seasonality in Madagascar using health facility data. BMC Med. 2020, 18, 26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- National Statistical Bureau. 2017 Population and Housing Census of Bhutan. In National Stastical Bureau RGoB; Royal Government of Bhutan: Thimphu, Bhutan, 2017. [Google Scholar]
- National Statistics Bureau. Dzongkhag Population Projection 2006–2015; National Statistics Bureau, Royal Government of Bhutan: Thimphu, Bhutan, 2008.
- Office of the Census Commissioner. Results of Population & Housing, Census of Bhutan 2005; Royal Government of Bhutan: Thimphu, Bhutan, 2006.
- Cleveland, R.B. STL: A Seasonal-Trend decomposition Prodecures Based on Loess. J. Off. Stat. 1990, 6, 3–73. [Google Scholar]
- Sang, S.; Gu, S.; Bi, P.; Yang, W.; Yang, Z.; Xu, L.; Yang, J.; Liu, X.; Jiang, T.; Wu, H.; et al. Predicting unprecedented dengue outbreak using imported cases and climatic factors in Guangzhou, 2014. PLoS Negl. Trop. Dis. 2015, 9, e0003808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Wangdi, K.; Clements, A.C.A.; Du, T.; Nery, S.V. Spatial and temporal patterns of dengue infections in Timor-Leste, 2005–2013. Parasites Vectors 2018, 11, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Kulldorff, M. A spatial scan statistic. Commun. Stat. Theory Methods 1997, 26, 1481–1496. [Google Scholar] [CrossRef]
- Ozdenerol, E.; Williams, B.L.; Kang, S.Y.; Magsumbol, M.S. Comparison of spatial scan statistic and spatial filtering in estimating low birth weight clusters. Int. J. Health Geogr. 2005, 4, 19. [Google Scholar] [CrossRef] [Green Version]
- Kim, S.; Jung, I. Optimizing the maximum reported cluster size in the spatial scan statistic for ordinal data. PLoS ONE 2017, 12, e0182234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Wangdi, K.; Xu, Z.; Suwannatrai, A.T.; Kurscheid, J.; Lal, L.; Namgay, R.; Glass, K.; Gray, D.J.; Archie, C.A. Clements A spatio-temporal analysis to identify the drivers of malaria transmission in Bhutan. Sci. Rep. 2020, 10, 7060. [Google Scholar] [CrossRef]
- Wangdi, K.; Singhasivanon, P.; Silawan, T.; Lawpoolsri, S.; White, N.J.; Kaewkungwal, J. Development of temporal modelling for forecasting and prediction of malaria infections using time-series and ARIMAX analyses: A case study in endemic districts of Bhutan. Malar. J. 2010, 9, 251. [Google Scholar] [CrossRef] [Green Version]
- Wangdi, K.; Gatton, M.; Kelly, G.; Clements, A. Prevalence of asymptomatic malaria and bed net ownership and use in Bhutan, 2013: A country earmarked for malaria elimination. Malar. J. 2014, 13, 352. [Google Scholar] [CrossRef] [Green Version]
- Wangdi, K.; Gatton, M.L.; Kelly, G.C.; Banwell, C.; Dev, V.; Clements, A.C. Malaria elimination in India and regional implications. Lancet Infect. Dis. 2016, 16, e214–e224. [Google Scholar] [CrossRef]
- Nath, D.C.; Mwchahary, D.D. Association between climatic variables and malaria incidence: A study in Kokrajhar district of Assam, India. Glob. J. Health Sci. 2012, 5, 90–106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- World Health Organization-South-East Asia Regional Office. Meeting on cross Border Collaboration on Malaria Elimination along the India-Bhutan Border; World Health Organization: Geneva, Switzerland, 2019. [Google Scholar]
- Wangchuk, S.; Drukpa, T.; Penjor, K.; Peldon, T.; Dorjey, Y.; Dorji, K.; Chhetri, V.; Trimarsanto, H.; To, S.; Murphy, A.; et al. Where chloroquine still works: The genetic make-up and susceptibility of Plasmodium vivax to chloroquine plus primaquine in Bhutan. Malar. J. 2016, 15, 277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Douglas, N.M.; Poespoprodjo, J.R.; Patriani, D.; Malloy, M.J.; Kenangalem, E.; Sugiarto, P.; Simpson, J.A.; Soenarto, Y.; Anstey, N.M.; Price, R.N. Unsupervised primaquine for the treatment of Plasmodium vivax malaria relapses in southern Papua: A hospital-based cohort study. PLoS Med. 2017, 14, e1002379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1.
Map of Bhutan with malaria-transmitting districts and sub-districts.
Figure 2.
Decomposed (a) Plasmodium falciparum and (b) Plasmodium vivax cases of Bhutan, 2010–2019.
Figure 3.
Spatial clusters of Plasmodium falciparum malaria in Bhutan by years, 2010–2019. (Red- most likely cluster; Green- secondary cluster).
Figure 3.
Spatial clusters of Plasmodium falciparum malaria in Bhutan by years, 2010–2019. (Red- most likely cluster; Green- secondary cluster).
Figure 4.
Spatial clusters of Plasmodium vivax malaria in Bhutan by years, 2010–2019. (Red- most likely cluster; Green- secondary cluster).
Figure 4.
Spatial clusters of Plasmodium vivax malaria in Bhutan by years, 2010–2019. (Red- most likely cluster; Green- secondary cluster).
Figure 5.
Spatiotemporal clusters of (a) Plasmodium falciparum and (b) Plasmodium vivax, in Bhutan 2010–2019.
Figure 5.
Spatiotemporal clusters of (a) Plasmodium falciparum and (b) Plasmodium vivax, in Bhutan 2010–2019.
Table 1.
Socio-demographic characteristics of Plasmodium falciparum and Plasmodium vivax during the study period.
Table 1.
Socio-demographic characteristics of Plasmodium falciparum and Plasmodium vivax during the study period.
| Characteristics | Plasmodium falciparum (%) | Plasmodium vivax (%) | Total (%) |
|---|---|---|---|
| Age (years) | |||
| ≤1 | 5 (1.6) | 5 (1.1) | 10 (1.3) |
| 2–12 | 63 (20.1) | 75 (16.5) | 138 (18.0) |
| 13–20 | 61 (19.4) | 75 (16.5) | 136 (17.7) |
| 21–45 | 127 (40.5) | 203 (44.7) | 330 (43.0) |
| >45 | 58 (18.5) | 96 (21.2) | 154 (20.1) |
| Sex | |||
| Female | 107 (34.1) | 161 (35.5) | 268 (34.9) |
| Male | 207 (65.9) | 293 (64.5) | 500 (65.1) |
| Occupation | |||
| Business | 4 (1.3) | 9 (2.0) | 13 (1.7) |
| Dependent | 23 (7.3) | 23 (5.1) | 46 (6.0) |
| Farmer | 146 (46.5) | 217 (47.8) | 363 (47.3) |
| House wife | 24 (7.6) | 40 (8.8) | 64 (8.3) |
| Labour | 5 (1.6) | 11 (2.4) | 16 (2.1) |
| Monk | 2 (0.6) | 4 (0.9) | 6 (0.8) |
| Public servant | 7 (2.2) | 26 (8.1) | 33 (4.3) |
| Armed forces | 17 (5.4) | 18 (2.4) | 35 (4.6) |
| Student | 86 (27.4) | 106 (23.4) | 192 (25.0) |
| Year | |||
| 2010 | 152 (48.4) | 249 (54.9) | 401 (52.2) |
| 2011 | 84 (26.8) | 74 (16.3) | 158 (20.6) |
| 2012 | 32 (10.2) | 30 (6.6) | 62 (8.1) |
| 2013 | 6 (1.9) | 6 (1.3) | 12 (1.6) |
| 2014 | 11 (3.5) | 8 (1.8) | 19 (2.5) |
| 2015 | 13 (4.1) | 22 (4.9) | 35 (4.6) |
| 2016 | 2 (0.6) | 17 (3.7) | 19 (2.5) |
| 2017 | 1 (0.3) | 23 (5.1) | 24 (3.1) |
| 2018 | 8 (2.6) | 12 (2.6) | 20 (2.6) |
| 2019 | 5 (1.6) | 13 (2.9) | 18 (2.3) |
| District | |||
| Chukha | 5 (1.6) | 22 (4.9) | 27 (3.5) |
| Dagana | 26 (8.3) | 17 (3.7) | 43 (5.6) |
| Pemagatshel | 6 (1.9) | 13 (2.9) | 19 (2.5) |
| Samdrup Jongkhar | 43 (13.7) | 36 (7.9) | 79 (10.3) |
| Samtse | 13 (4.1) | 52 (11.5) | 65 (8.5) |
| Sarpang | 217 (69.1) | 306 (67.4) | 523 (68.1) |
| Zhemgang | 4 (1.3) | 8 (1.8) | 12 (1.6) |
Table 2.
Spatial clusters of Plasmodium falciparum in Bhutan, 2010–2019.
| Year | SaTScan Statistics | Long | Lat | Radius (km) | Population | Number of Cases | Expected Cases | No of Sub-Districts | RR | LLR | p-Value |
|---|---|---|---|---|---|---|---|---|---|---|---|
| 2010 | Most likely cluster | 90.2587 | 26.9402 | 50.36 | 70,171 | 120 | 36.51 | 26 | 11.86 | 101.71 | <0.001 |
| Secondary cluster | 91.9991 | 26.9079 | 0 | 2067 | 13 | 1.08 | 1 | 13.12 | 20.96 | <0.001 | |
| 2011 | Most likely cluster | 90.3536 | 26.9405 | 29.51 | 41,701 | 60 | 11.79 | 11 | 15.31 | 71.18 | <0.001 |
| 2012 | Most likely cluster | 90.3536 | 26.9405 | 29.51 | 42,535 | 26 | 4.51 | 11 | 26.44 | 36.44 | <0.001 |
| Secondary cluster | 91.9991 | 26.9079 | 0 | 2143 | 4 | 0.23 | 1 | 19.99 | 7.94 | <0.001 | |
| 2013 | Most likely cluster | 90.3435 | 27.0545 | 25.79 | 37,861 | 6 | 0.74 | 8 | Infinity | 12.56 | <0.001 |
| 2014 | Most likely cluster | 91.9991 | 26.9079 | 0 | 2220 | 6 | 0.078 | 1 | 167.64 | 22.14 | <0.001 |
| Secondary cluster | 90.3536 | 26.9405 | 9.4 | 14,771 | 5 | 0.52 | 3 | 16.79 | 7.97 | <0.016 | |
| 2015 | Most likely cluster | 90.4328 | 26.9283 | 15.71 | 32,603 | 11 | 1.33 | 6 | 48.10 | 19.68 | <0.001 |
| 2016 | No cluster | ||||||||||
| 2017 | No cluster | ||||||||||
| 2018 | Most likely cluster | 90.4845 | 26.9197 | 20.26 | 410,140 | 7 | 0.98 | 8 | 50.01 | 11.80 | <0.001 |
| 2019 | Most likely cluster | 90.5214 | 26.8637 | 0 | 2428 | 3 | 0.036 | 1 | 208.35 | 11.47 | <0.001 |
Lat—latitude; long—longitude; RR—relative risk; LLR—log-likelihood ratio.
Table 3.
Spatial clusters of Plasmodium vivax at sub-district level in Bhutan, 2006–2019.
| Year | SaTScan Statistics | Long | Lat | Radius (km) | Population | Number of Cases | Expected Cases | No of Sub-Districts | RR | LLR | p-Value |
|---|---|---|---|---|---|---|---|---|---|---|---|
| 2010 | Most likely cluster | 90.3536 | 26.9405 | 33.61 | 45,511 | 1166 | 38.79 | 14 | 10.84 | 164.19 | <0.001 |
| Secondary cluster | 91.6777 | 26.8583 | 9.49 | 4652 | 19 | 3.97 | 2 | 5.10 | 15.21 | <0.001 | |
| 2011 | Most likely cluster | 90.4328 | 26.9283 | 26.02 | 40,946 | 45 | 10.2 | 11 | 9.71 | 43.93 | <0.001 |
| 2012 | Most likely cluster | 90.6599 | 26.8108 | 53.39 | 65,886 | 28 | 6.54 | 26 | 50.19 | 35.78 | <0.001 |
| 2013 | Most likely cluster | 90.4845 | 26.9197 | 20.26 | 37,147 | 6 | 0.73 | 8 | Infinity | 12.78 | <0.001 |
| 2014 | Most likely cluster | 90.4845 | 26.99197 | 20.26 | 37,913 | 7 | 0.97 | 8 | 50.68 | 11.88 | <0.001 |
| 2015 | Most likely cluster | 90.2869 | 26.8729 | 28.05 | 42,318 | 20 | 2.93 | 10 | 65.08 | 33.90 | <0.001 |
| 2016 | Most likely cluster | 90.2869 | 26.8729 | 28.05 | 43,161 | 14 | 2.27 | 10 | 30.26 | 20.69 | <0.001 |
| 2017 | Most likely cluster | 89.1271 | 26.9157 | 0 | 9974 | 10 | 0.70 | 1 | 24.56 | 19.60 | <0.001 |
| 2018 | Most likely cluster | 90.5628 | 26.9293 | 8.99 | 24,825 | 8 | 0.89 | 4 | 24.91 | 13.47 | <0.001 |
| 2019 | Most likely cluster | 90.1291 | 26.7959 | 18.46 | 8678 | 7 | 0.33 | 4 | 44.50 | 16.85 | <0.001 |
Lat—latitude; long—longitude; RR—relative risk; LLR—log-likelihood ratio.
Table 4.
Significant Plasmodium falciparum space–time clusters at sub-district level in Bhutan, 2010–2019.
Table 4.
Significant Plasmodium falciparum space–time clusters at sub-district level in Bhutan, 2010–2019.
| Time Period (Month Year) | SaTScan Statistics | Long | Lat | Radius (km) | Population | Number of Cases | Expected Cases | No of Sub-Districts | RR | LLR | p-Value |
|---|---|---|---|---|---|---|---|---|---|---|---|
| January 2010–June 2012 | Most likely cluster | 90.2587 | 26.9402 | 50.36 | 76,572 | 207 | 17.71 | 26 | 32.66 | 399.91 | <0.001 |
| January 2010–December 2014 | Secondary cluster | 91.9991 | 26.9079 | 0 | 2243 | 25 | 1.07 | 1 | 25.37 | 55.86 | <0.001 |
| April 2010 | Secondary cluster | 91.3444 | 27.0405 | 56.81 | 72,597 | 9 | 0.55 | 26 | 16.69 | 16.76 | <0.001 |
| April 2010 | Secondary cluster | 88.8911 | 26.9823 | 8.33 | 11,685 | 5 | 0.089 | 3 | 57.15 | 15.28 | 0.0057 |
Lat—latitude; long—longitude; RR—relative risk; LLR—log-likelihood ratio.
Table 5.
Significant Plasmodium vivax space–time clusters at sub-district level in Bhutan, 2010–2019.
Table 5.
Significant Plasmodium vivax space–time clusters at sub-district level in Bhutan, 2010–2019.
| Time Period (Month Year) | SaTScan Statistics | Long | Lat | Radius (km) | Population | Number of Cases | Expected Cases | No of Sub-Districts | RR | LLR | p-Value |
|---|---|---|---|---|---|---|---|---|---|---|---|
| January 2010–June 2012 | Most likely cluster | 90.3536 | 26.9405 | 33.61 | 49,858 | 233 | 16.62 | 11 | 27.74 | 464.34 | <0.001 |
| March–May 2010 | Secondary cluster | 91.6777 | 26.8583 | 9.49 | 5048 | 17 | 0.17 | 2 | 104.65 | 61.91 | <0.001 |
| March–May 2010 | Secondary cluster | 88.9128 | 27.0651 | 13.09 | 25,715 | 15 | 0.87 | 7 | 17.85 | 28.85 | <0.001 |
| August–November 2017 | Secondary cluster | 89.1271 | 26.9157 | 0 | 9596 | 10 | 0.48 | 1 | 21.30 | 20.95 | <0.001 |
Lat—latitude; long—longitude; RR—relative risk; LLR—log-likelihood ratio.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Wangdi, K.; Penjor, K.; Tobgyal; Lawpoolsri, S.; Price, R.N.; Gething, P.W.; Gray, D.J.; Da Silva Fonseca, E.; Clements, A.C.A. Space–Time Clustering Characteristics of Malaria in Bhutan at the End Stages of Elimination. Int. J. Environ. Res. Public Health 2021, 18, 5553. https://doi.org/10.3390/ijerph18115553
AMA Style
Wangdi K, Penjor K, Tobgyal, Lawpoolsri S, Price RN, Gething PW, Gray DJ, Da Silva Fonseca E, Clements ACA. Space–Time Clustering Characteristics of Malaria in Bhutan at the End Stages of Elimination. International Journal of Environmental Research and Public Health. 2021; 18(11):5553. https://doi.org/10.3390/ijerph18115553
Chicago/Turabian StyleWangdi, Kinley, Kinley Penjor, Tobgyal, Saranath Lawpoolsri, Ric N. Price, Peter W. Gething, Darren J. Gray, Elivelton Da Silva Fonseca, and Archie C. A. Clements. 2021. "Space–Time Clustering Characteristics of Malaria in Bhutan at the End Stages of Elimination" International Journal of Environmental Research and Public Health 18, no. 11: 5553. https://doi.org/10.3390/ijerph18115553
APA StyleWangdi, K., Penjor, K., Tobgyal, Lawpoolsri, S., Price, R. N., Gething, P. W., Gray, D. J., Da Silva Fonseca, E., & Clements, A. C. A. (2021). Space–Time Clustering Characteristics of Malaria in Bhutan at the End Stages of Elimination. International Journal of Environmental Research and Public Health, 18(11), 5553. https://doi.org/10.3390/ijerph18115553
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
