Spatio-temporal patterns of malaria infection in Bhutan: a country embarking on malaria elimination

Respiratory syncytial virus among hospitalized patients of severe acute respiratory infection in Bhutan: Cross-sectional study

Introduction

Respiratory syncytial virus (RSV) is a leading cause of lower respiratory tract infections worldwide, particularly in young children. In Bhutan, respiratory disease continues to be among the top 10 diseases of morbidity for several years. This study aimed to estimate the prevalence of RSV among hospitalized patients with severe acute respiratory infection (SARI) in Bhutan.

Method

Respiratory specimens were collected from SARI patients of all ages in 2016 and 2018 following influenza surveillance guidelines. Specimens were tested for influenza and RSV, human metapneumovirus, adenovirus, and human parainfluenza virus types 1, 2, and 3 using real-time reverse-transcription polymerase chain reaction assay. Descriptive statistics were used to analyze the result in STATA 16.1.

Result

Of the 1339 SARI specimens tested, 34.8% were positive for at least one viral pathogen. RSV was detected in 18.5% of SARI cases, followed by influenza in 13.4% and other respiratory viruses in 3%. The median age of SARI cases was 3 (IQR: 0.8-21 years) years. RSV detection was higher among children aged 0-6 (Adj OR: 3.03; 95% CI: 1.7-5.39) and 7-23 months (Adj OR: 3.01; 95% CI: 1.77-5.12) compared with the children aged 5-15 years. RSV was also associated with breathing difficulty (Adj OR: 1.73; 95% CI: 1.17-2.56) and pre-existing lung disease, including asthma (Adj OR: 2.78; 95% CI: 0.99-7.8).

Conclusion

Respiratory viruses were detected in a substantial proportion of SARI hospitalizations in Bhutan.

Assessment of malaria risk in Southeast Asia: a systematic review

BACKGROUND

Several countries in Southeast Asia are nearing malaria elimination, yet eradication remains elusive. This is largely due to the challenge of focusing elimination efforts, an area where risk prediction can play an essential supporting role. Despite its importance, there is no standard numerical method to quantify the risk of malaria infection. Thus, there is a need for a consolidated view of existing definitions of risk and factors considered in assessing risk to analyse the merits of risk prediction models. This systematic review examines studies of the risk of malaria in Southeast Asia with regard to their suitability in addressing the challenges of malaria elimination in low transmission areas.

METHODS

A search of four electronic databases over 2010-2020 retrieved 1297 articles, of which 25 met the inclusion and exclusion criteria. In each study, examined factors included the definition of the risk and indicators of malaria transmission used, the environmental and climatic factors associated with the risk, the statistical models used, the spatial and temporal granularity, and how the relationship between environment, climate, and risk is quantified.

RESULTS

This review found variation in the definition of risk used, as well as the environmental and climatic factors in the reviewed articles. GLM was widely adopted as the analysis technique relating environmental and climatic factors to malaria risk. Most of the studies were carried out in either a cross-sectional design or case-control studies, and most utilized the odds ratio to report the relationship between exposure to risk and malaria prevalence.

CONCLUSIONS

Adopting a standardized definition of malaria risk would help in comparing and sharing results, as would a clear description of the definition and method of collection of the environmental and climatic variables used. Further issues that need to be more fully addressed include detection of asymptomatic cases and considerations of human mobility. Many of the findings of this study are applicable to other low-transmission settings and could serve as a guideline for further studies of malaria in other regions.

Spatio-temporal trends of malaria incidence from 2011 to 2017 and environmental predictors of malaria transmission in Myanmar

BACKGROUND

Myanmar bears the heaviest malaria burden in the Greater Mekong Subregion (GMS). This study assessed the spatio-temporal dynamics and environmental predictors of Plasmodium falciparum and Plasmodium vivax malaria in Myanmar.

METHODS

Monthly reports of malaria cases at primary health centers during 2011-2017 were analyzed to describe malaria distribution across Myanmar at the township and state/region levels by spatial autocorrelation (Moran index) and spatio-temporal clustering. Negative binomial generalized additive models identified environmental predictors for falciparum and vivax malaria, respectively.

RESULTS

From 2011 to 2017, there was an apparent reduction in malaria incidence in Myanmar. Malaria incidence peaked in June each year. There were significant spatial autocorrelation and clustering with extreme spatial heterogeneity in malaria cases and test positivity across the nation (P < 0.05). Areas with higher malaria incidence were concentrated along international borders. Primary clusters of P. falciparum persisted in western townships, while clusters of P. vivax shifted geographically over the study period. The primary cluster was detected from January 2011 to December 2013 and covered two states (Sagaing and Kachin). Annual malaria incidence was highest in townships with a mean elevation of 500‒600 m and a high variance in elevation (states with both high and low elevation). There was an apparent linear relationship between the mean normalized difference vegetative index and annual P. falciparum incidence (P < 0.05).

CONCLUSION

The decreasing trends reflect the significant achievement of malaria control efforts in Myanmar. Prioritizing the allocation of resources to high-risk areas identified in this study can achieve effective disease control.

Distribution and Risk Factors of Malaria in the Greater Accra Region in Ghana

Malaria remains a serious public health challenge in Ghana including the Greater Accra Region. This study aimed to quantify the spatial, temporal and spatio-temporal patterns of malaria in the Greater Accra Region to inform targeted allocation of health resources. Malaria cases data from 2015 to 2019 were obtained from the Ghanaian District Health Information and Management System and aggregated at a district and monthly level. Spatial analysis was conducted using the Global Moran's I, Getis-Ord Gi*, and local indicators of spatial autocorrelation. Kulldorff's space-time scan statistics were used to investigate space-time clustering. A negative binomial regression was used to find correlations between climatic factors and sociodemographic characteristics and the incidence of malaria. A total of 1,105,370 malaria cases were reported between 2015 and 2019. Significant seasonal variation was observed, with June and July being the peak months of reported malaria cases. The hotspots districts were Kpone-Katamanso Municipal District, Ashaiman Municipal Districts, Tema Municipal District, and La-Nkwantanang-Madina Municipal District. While La-Nkwantanang-Madina Municipal District was high-high cluster. The Spatio-temporal clusters occurred between February 2015 and July 2017 in the districts of Ningo-Prampram, Shai-Osudoku, Ashaiman Municipal, and Kpone-Katamanso Municipal with a radius of 26.63 km and an relative risk of 4.66 (p < 0.001). Malaria cases were positively associated with monthly rainfall (adjusted odds ratio [AOR] = 1.01; 95% confidence interval [CI] = 1.005, 1.016) and the previous month's cases (AOR = 1.064; 95% CI 1.062, 1.065) and negatively correlated with minimum temperature (AOR = 0.86, 95% CI = 0.823, 0.899) and population density (AOR = 0.996, 95% CI = 0.994, 0.998). Malaria control and prevention should be strengthened in hotspot districts in the appropriate months to improve program effectiveness.

Spatial analysis to evaluate risk of malaria in Northern Sumatera, Indonesia

Background

As Indonesia aims for malaria elimination by 2030, provisional malaria epidemiology and risk factors evaluation are important in pursue of this national goal. Therefore, this study aimed to understand the risk factor of malaria in Northern Sumatera.

Methods

Malaria cases from 2019 to 2020 were obtained from the Indonesian Ministry of Health Electronic Database. Climatic variables were provided by the Center for Meteorology and Geophysics Medan branch office. Multivariable logistic regression was undertaken to understand the risk factors of imported malaria. A zero-inflated Poisson multivariable regression model was used to study the climatic drivers of indigenous malaria.

Results

A total of 2208 (indigenous: 76.0% [1679] and imported: 17.8% [392]) were reported during the study period. Risk factors of imported malaria were: ages 19–30 (adjusted odds ratio [AOR] = 3.31; 95% confidence interval [CI] 1.67, 2.56), 31–45 (AOR = 5.69; 95% CI 2.65, 12.20), and > 45 years (AOR = 5.11; 95% CI 2.41, 10.84). Military personnel and forest workers and miners were 1,154 times (AOR = 197.03; 95% CI 145.93, 9,131.56) and 44 times (AOR = 44.16; 95% CI 4.08, 477,93) more likely to be imported cases as compared to those working as employees and traders. Indigenous Plasmodium falciparum increased by 12.1% (95% CrI 5.1%, 20.1%) for 1% increase in relative humidity and by 21.0% (95% CrI 9.0%, 36.2%) for 1 °C increase in maximum temperature. Plasmodium vivax decreased by 0.8% (95% CrI 0.2%, 1.3%) and 16.7% (95% CrI 13.7%, 19.9%) for one meter and 1 °C increase of altitude and minimum temperature. Indigenous hotspot was reported by Kota Tanjung Balai city and Asahan regency, respectively. Imported malaria hotspots were reported in Batu Bara, Kota Tebing Tinggi, Serdang Bedagai and Simalungun.

Conclusion

Both indigenous and imported malaria is limited to a few regencies and cities in Northern Sumatera. The control measures should focus on these risk factors to achieve elimination in Indonesia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04262-y.

Cross-border malaria drivers and risk factors on the Brazil–Venezuela border between 2016 and 2018

Globally, cross-border importation of malaria has become a challenge to malaria elimination. The border areas between Brazil and Venezuela have experienced high numbers of imported cases due to increased population movement and migration out of Venezuela. This study aimed to identify risk factors for imported malaria and delineate imported malaria hotspots in Roraima, Brazil and Bolivar, Venezuela between 2016 and 2018. Data on malaria surveillance cases from Roraima, Brazil and Bolivar, Venezuela from 2016 to 2018 were obtained from national surveillance systems: the Brazilian Malaria Epidemiology Surveillance Information System (SIVEP-Malaria), the Venezuelan Ministry of Health and other non-government organizations. A multivariable logistic regression model was used to identify the risk factors for imported malaria. Spatial autocorrelation in malaria incidence was explored using Getis-Ord (Gi*) statistics. During the study period, there were 11,270 (24.3%) and 4072 (0.7%) imported malaria cases in Roraima, Brazil and Bolivar, Venezuela, respectively. In the multivariable logistic regression for Roraima, men were 28% less likely to be an imported case compared to women (Adjusted Odds Ratio [AOR] = 0.72; 95% confidence interval [CI] 0.665, 0.781). Ages 20–29 and 30–39 were 90% (AOR = 1.90; 95% CI 1.649, 2.181) and 54% (AOR = 1.54; 95% CI 1.331, 1.782) more likely to be an imported case compared to the 0–9 year age group, respectively. Imported cases were 197 times (AOR = 197.03; 95% CI 175.094, 221.712) more likely to occur in miners than those working in agriculture and domestic work. In Bolivar, cases aged 10–19 (AOR = 1.75; 95% CI 1.389, 2.192), 20–29 (AOR = 2.48; 95% CI 1.957, 3.144), and 30–39 (AOR = 2.29; 95% CI 1.803, 2.913) were at higher risk of being an imported case than those in the 0–9 year old group, with older age groups having a slightly higher risk compared to Roraima. Compared to agriculture and domestic workers, tourism, timber and fishing workers (AOR = 6.38; 95% CI 4.393, 9.254) and miners (AOR = 7.03; 95% CI 4.903, 10.092) were between six and seven times more likely to be an imported case. Spatial analysis showed the risk was higher along the international border in the municipalities of Roraima, Brazil. To achieve malaria elimination, cross-border populations in the hotspot municipalities will need targeted intervention strategies tailored to occupation, age and mobility status. Furthermore, all stakeholders, including implementers, policymakers, and donors, should support and explore the introduction of novel approaches to address these hard-to-reach populations with the most cost-effective interventions.

Spatio-temporal modelling of malaria mortality in India from 2004 to 2013 from the Million Death Study

BACKGROUND

India has a substantial burden of malaria, concentrated in specific areas and population groups. Spatio-temporal modelling of deaths due to malaria in India is a critical tool for identifying high-risk groups for effective resource allocation and disease control policy-making, and subsequently for the country's progress towards United Nations 2030 Sustainable Development Goals.

METHODS

In this study, a spatio-temporal model with the objective of understanding the spatial distribution of malaria mortality rates and the rate of temporal decline, across the country, has been constructed. A spatio-temporal "random slope" model was used, with malaria risk depending on a spatial relative risk surface and a linear time effect with a spatially-varying coefficient. The models were adjusted for urban/rural status (residence of the deceased) and Normalized Difference Vegetation Index (NDVI), using 2004-13 data from the Million Death Study (MDS) (the most recent data available), with nationwide geographic coverage. Previous studies based on MDS had focused only on aggregated analyses.

RESULTS

The rural population had twice the risk of death due to malaria compared to the urban population. Malaria mortality in some of the highest-risk regions, namely the states of Odisha and Jharkhand, are declining faster than other areas; however, the rate of decline was not uniformly correlated with the level of risk. The overall decline was faster after 2010.

CONCLUSION

The results suggest a need for increased attention in high-risk rural populations, which already face challenges like inadequate infrastructure, inaccessibility to health care facilities, awareness, and education around malaria mortality and prevalence. It also points to the urgent need to restart the MDS to document changes since 2013, to develop appropriate malaria control measures.

Space-Time Clustering Characteristics of Malaria in Bhutan at the End Stages of Elimination

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.

Spatial and Temporal Patterns of Malaria in Phu Yen Province, Vietnam, from 2005 to 2016

Malaria in Vietnam has become focal to a few provinces, including Phu Yen. This study aimed to assess correlations between intervention (population proportion protected by insecticide-treated nets and indoor residual spraying) and climatic variables with malaria incidence in Phu Yen Province. The Vietnam National Institute of Malariology, Parasitology, and Entomology provided incidence data for Plasmodium falciparum and Plasmodium vivax for 104 communes of Phu Yen Province from January 2005 to December 2016. A multivariable, zero-inflated Poisson regression model was developed with a conditional autoregressive prior structure to identify the underlying spatial structure of the data and quantify associations with covariates. There were a total of 2,778 P. falciparum and 1,770 P. vivax cases during the study period. Plasmodium falciparum and P. vivax incidence increased by 5.4% (95% credible interval [CrI] 5.1%, 5.7%) and 3.2% (95% CrI 2.9%, 3.5%) for a 10-mm increase in precipitation without lag, respectively. Plasmodium falciparum and P. vivax incidence decreased by 7.7% (95% CrI 5.6%, 9.7%) and 10.5% (95% CrI 8.3%, 12.6%) for a 1°C increase in minimum temperature without lag, respectively. There was a > 95% probability of a higher than provincial average trend of P. falciparum and P. vivax in Song Cau and Song Hoa districts. There was a > 95% probability of a lower than provincial average trend in Tuy Dong Xuan and Hoa districts for both species. Targeted distribution of resources, including intensified interventions, in this part of the province will be required for local malaria elimination.

A spatio-temporal analysis to identify the drivers of malaria transmission in Bhutan

At a time when Bhutan is on the verge of malaria elimination, the aim of this study was to identify malaria clusters at high geographical resolution and to determine its association with local environmental characteristics. Malaria cases from 2006–2014 were obtained from the Vector-borne Disease Control Program under the Ministry of Health, Bhutan. A Zero-Inflated Poisson multivariable regression model with a conditional autoregressive (CAR) prior structure was developed. Bayesian Markov chain Monte Carlo (MCMC) simulation with Gibbs sampling was used to estimate posterior parameters. A total of 2,062 Plasmodium falciparum and 2,284 Plasmodium vivax cases were reported during the study period. Both species of malaria showed seasonal peaks with decreasing trend. Gender and age were not associated with the transmission of either species of malaria. P. falciparum increased by 0.7% (95% CrI: 0.3%, 0.9%) for a one mm increase in rainfall, while climatic variables (temperature and rainfall) were not associated with P. vivax. Insecticide treated bed net use and residual indoor insecticide coverage were unaccounted for in this study. Hot spots and clusters of both species were isolated in the central southern part of Bhutan bordering India. There was significant residual spatial clustering after accounting for climate and demographic variables.

Malaria elimination in Bhutan: asymptomatic malaria cases in the Bhutanese population living in malaria-risk areas and in migrant workers from India

In 2018, Bhutan reported 54 cases of malaria, of which six were indigenous, 14 introduced and 34 imported. Considering the continuous reduction in the number of indigenous cases, Bhutan plans to eliminate malaria by 2025 under the Bhutan Malaria Elimination Strategy. The study was conducted to assess the presence of asymptomatic plasmodial infection in both, Bhutanese population living in malaria-risk areas and in migrant workers to guide the elimination strategies. A cross-sectional study was conducted from April to May 2016 in 750 Bhutanese people and 473 migrant workers. Plasmodium falciparum and Plasmodium vivax infections were investigated by using a rapid diagnostic test (RDT) and the polymerase chain reaction (PCR). Prevalence of asymptomatic plasmodial infection based on PCR was 0.27% (95% CI: 0.05–1.07%) among Bhutanese people with a mean age of 43 years old. The proportions of males and females were 45% and 55%, respectively. Among migrant workers, the prevalence of asymptomatic plasmodial infection was 0.42% (95% CI: 0.07–1.69%) with a mean age of 30 years old. The majority of migrant workers were from the neighboring Indian State of West Bengal (57.51%), followed by Assam (12.26%). RDT in both study groups did not detect any plasmodial infection. The presence of a low prevalence of asymptomatic plasmodial infection indicates that the current elimination strategies and interventions are effective.

Malaria Hotspots: Is There Epidemiological Evidence for Fine-Scale Spatial Targeting of Interventions?

As data at progressively granular spatial scales become available, the temptation is to target interventions to areas with higher malaria transmission - so-called hotspots - with the aim of reducing transmission in the wider community. This paper reviews literature to determine if hotspots are an intrinsic feature of malaria epidemiology and whether current evidence supports hotspot-targeted interventions. Hotspots are a consistent feature of malaria transmission at all endemicities. The smallest spatial unit capable of supporting transmission is the household, where peri-domestic transmission occurs. Whilst the value of focusing interventions to high-burden areas is evident, there is currently limited evidence that local-scale hotspots fuel transmission. As boundaries are often uncertain, there is no conclusive evidence that hotspot-targeted interventions accelerate malaria elimination.

Spatial Heterogeneity and Temporal Trends in Malaria on the Thai⁻Myanmar Border (2012⁻2017): A Retrospective Observational Study

Malaria infections remain an important public health problem for the Thai-Myanmar border population, despite a plan for the elimination by the end of 2026 (Thailand) and 2030 (Myanmar). This study aimed to explore spatiotemporal patterns in Plasmodium falciparum and Plasmodium vivax incidence along the Thai-Myanmar border. Malaria cases among Thai citizens in 161 sub-districts in Thailand's Kanchanaburi and Tak Provinces (2012-2017) were analyzed to assess the cluster areas and temporal trends. Based on reported incidence, 65.22% and 40.99% of the areas studied were seen to be at elimination levels for P. falciparum and P. vivax already, respectively. There were two clear clusters of malaria in the region: One in the northern part (Cluster I), and the other in the central part (Cluster II). In Cluster I, the malaria season exhibited two peaks, while there was only one peak seen for Cluster II. Malaria incidence decreased at a faster rate in Cluster I, with 5% and 4% reductions compared with 4% and 3% reductions in P. falciparum and P. vivax incidence per month, respectively, in Cluster II. The decreasing trends reflect the achievements of malaria control efforts on both sides of the Thai-Myanmar border. However, these clusters could act as reservoirs. Perhaps one of the main challenges facing elimination programs in this low transmission setting is maintaining a strong system for early diagnosis and treatment, even when malaria cases are very close to zero, whilst preventing re-importation of cases.

Bayesian spatiotemporal analysis of malaria infection along an international border: Hlaingbwe Township in Myanmar and Tha-Song-Yang District in Thailand

BACKGROUND

One challenge in moving towards malaria elimination is cross-border malaria infection. The implemented measures to prevent and control malaria re-introduction across the demarcation line between two countries require intensive analyses and interpretation of data from both sides, particularly in border areas, to make correct and timely decisions. Reliable maps of projected malaria distribution can help to direct intervention strategies. In this study, a Bayesian spatiotemporal analytic model was proposed for analysing and generating aggregated malaria risk maps based on the exceedance probability of malaria infection in the township-district adjacent to the border between Myanmar and Thailand. Data of individual malaria cases in Hlaingbwe Township and Tha-Song-Yang District during 2016 were extracted from routine malaria surveillance databases. Bayesian zero-inflated Poisson model was developed to identify spatial and temporal distributions and associations between malaria infections and risk factors. Maps of the descriptive statistics and posterior distribution of predicted malaria infections were also developed.

RESULTS

A similar seasonal pattern of malaria was observed in both Hlaingbwe Township and Tha-Song-Yang District during the rainy season. The analytic model indicated more cases of malaria among males and individuals aged ≥ 15 years. Mapping of aggregated risk revealed consistently high or low probabilities of malaria infection in certain village tracts or villages in interior parts of each country, with higher probability in village tracts/villages adjacent to the border in places where it could easily be crossed; some border locations with high mountains or dense forests appeared to have fewer malaria cases. The probability of becoming a hotspot cluster varied among village tracts/villages over the year, and some had close to no cases all year.

CONCLUSIONS

The analytic model developed in this study could be used for assessing the probability of hotspot cluster, which would be beneficial for setting priorities and timely preventive actions in such hotspot cluster areas. This approach might help to accelerate reaching the common goal of malaria elimination in the two countries.

Analysis of clinical malaria disease patterns and trends in Vietnam 2009-2015

BACKGROUND

Viet Nam has made tremendous progress towards reducing mortality and morbidity associated with malaria in recent years. Despite the success in malaria control, there has been a recent increase in cases in some provinces. In order to understand the changing malaria dynamics in Viet Nam and measure progress towards elimination, the aim of this study was to describe and quantify spatial and temporal trends of malaria by species at district level across the country.

METHODS

Malaria case reports at the Viet Nam National Institute of Malariology, Parasitology, and Entomology were reviewed for the period of January 2009 to December 2015. The population of each district was obtained from the Population and Housing Census-2009. A multivariate (insecticide-treated mosquito nets [ITN], indoor residual spraying [IRS], maximum temperature), zero-inflated, Poisson regression model was developed with spatial and spatiotemporal random effects modelled using a conditional autoregressive prior structure, and with posterior parameters estimated using Bayesian Markov chain Monte Carlo simulation with Gibbs sampling. Covariates included in the models were coverage of intervention (ITN and IRS) and maximum temperature.

RESULTS

There was a total of 57,713 Plasmodium falciparum and 32,386 Plasmodium vivax cases during the study period. The ratio of P. falciparum to P. vivax decreased from 4.3 (81.0% P. falciparum; 11,121 cases) in 2009 to 0.8 (45.0% P. falciparum; 3325 cases) in 2015. Coverage of ITN was associated with decreased P. falciparum incidence, with a 1.1% (95% credible interval [CrI] 0.009%, 1.2%) decrease in incidence for 1% increase in the ITN coverage, but this was not the case for P. vivax, nor was it the case for IRS coverage. Maximum temperature was associated with increased incidence of both species, with a 4% (95% CrI 3.5%, 4.3%) and 1.6% (95% CrI 0.9%, 2.0%) increase in P. falciparum and P. vivax incidence for a temperature increase of 1 °C, respectively. Temporal trends of P. falciparum and P. vivax incidence were significantly higher than the national average in Central and Central-Southern districts.

CONCLUSION

Interventions (ITN distribution) and environmental factors (increased temperature) were associated with incidence of P. falciparum and P. vivax during the study period. The factors reviewed were not exhaustive, however the data suggest distribution of resources can be targeted to areas and times of increased malaria transmission. Additionally, changing distribution of the two predominant malaria species in Viet Nam will require different programmatic approaches for control and elimination.

Travel to farms in the lowlands and inadequate malaria information significantly predict malaria in villages around Lake Tana, northwest Ethiopia: a matched case-control study

Background

In Ethiopia, malaria has declined in the last decade; only a small number of cases have been reported, primarily from hotspots. The contribution of house proximity to water bodies and the role of migration in malaria transmission has not yet been examined in detail in northwest Ethiopia. Individual and household-level environmental and socio-demographic drivers of malaria heterogeneity were explored contextually in meso-endemic villages around Lake Tana, northwest Ethiopia.

Methods

A health facility-based paired age-sex matched case–control study involving 303 matched pairs was undertaken from 10 October 2016, to 30 June 2017. Geo-referencing of case households, control households, proximate water bodies, and health centres was carried out. A pretested and structured questionnaire was used to collect data on socio-demography, household assets, housing, travel history, and malaria intervention measures. Medians (interquartile range) were computed for continuous variables. Pearson’s Chi square/Fisher’s exact test was used to detect significant differences in proportions. Principal component analysis was performed to estimate household wealth. Stratified analysis was used to confirm confounding and interaction. A multivariable conditional logistic regression model was used to detect risk factors for malaria.

Results

Of 303 malaria cases, 59 (19.5% [15.4–24.3]) were imported malaria cases whereas 244 (80.5% [75.7–84.6]) were locally acquired malaria cases. In bivariate analysis, marital status, educational status, and bed net ownership were significantly associated with malaria cases. In multivariable adjustment, travel to malarious lowlands in the preceding month (adjusted mOR = 7.32; 95% CI 2.40–22.34), household member’s travel to malarious lowlands (adjusted mOR = 2.75; 95% CI 1.02–7.44), and inadequate health information on malaria (adjusted mOR = 1.57; 95% CI 1.03–2.41) were predictors of malaria. Stratified analysis confirmed that elevation of households and travel to malarious lowlands were not effect modifiers. Travel to malarious lowlands had a confounding effect on malaria but elevation of households did not.

Conclusions

In this study, travel to farms in the lowlands and inadequate health information on malaria were risk factors for malaria in villages around Lake Tana. This evidence is critical for the design of improved strategic interventions that consider imported malaria cases and approaches for accessing health information on malaria control in northwest Ethiopia.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2434-y) contains supplementary material, which is available to authorized users.

Malaria burden and costs of intensified control in Bhutan, 2006-14: an observational study and situation analysis

INTRODUCTION

The number of malaria cases has fallen in Bhutan in the past two decades, and the country has a goal of complete elimination of malaria by 2016. The aims of this study are to ascertain the trends and burden of malaria, the costs of intensified control activities, the main donors of funding for the control activities, and the costs of different preventive measures in the pre-elimination phase (2006-14) in Bhutan.

METHODS

We undertook a descriptive analysis of malaria surveillance data from 2006 to 2014, using data from the Vector-borne Disease Control Programme (VDCP) run by the Department of Public Health of Bhutan's Ministry of Health. Malaria morbidity and mortality in local Bhutanese people and foreign nationals were analysed. The cost of different control and preventive measures were calculated, and the average numbers of long-lasting insecticidal nests per person were estimated.

FINDINGS

A total of 5491 confirmed malaria cases occurred in Bhutan between 2006 and 2014. By 2013, there was an average of one long-lasting insecticidal net for every 1·51 individuals. The cost of procuring long-lasting insecticidal nets accounted for more than 90% of the total cost of prevention measures. The Global Fund to Fight AIDS, Tuberculosis and Malaria was the main international donor, accounting for more than 80% of the total funds.

INTERPRETATION

The malaria burden in Bhutan decreased significantly during the study period with high coverage of long-lasting insecticidal nets. The foreseeable challenges that require national attention to maintain a malaria-free status after elimination are importation of malaria, especially from India; continued protection of the population in endemic districts through complete coverage with long-lasting insecticidal nets and indoor residual spraying; and exploration of local funding modalities post-elimination in the event of a reduction in international funding.

FUNDING

None.

An Evaluation of the National Malaria Surveillance System of Bhutan, 2006-2012 as It Approaches the Goal of Malaria Elimination

INTRODUCTION

Bhutan is progressing toward malaria elimination. The purpose of this evaluation was to assess the ability of the surveillance system from 2006 to 2012 to meet the objectives of the Bhutan Vector-borne Disease Control Program (VDCP) and to highlight priorities requiring attention as the nation transitions to elimination.

METHODS

The evaluation was conducted using the Center for Disease Control guidelines for evaluating public health surveillance systems. Data sources included a search of publically available literature, VDCP program data, and interviews with malaria surveillance personnel. Blood slide quality assurance and control through formal assessment of slide preparation and measures of between-reader correlation were performed.

RESULTS

Total malaria cases declined from 2006 to 2012. The average slide positivity rate decreased from 3.4% in 2006 to 0.2% in 2012. The proportion of non-residents in all cases increased to its highest value of 22.6% in 2012, and significant clustering in the border regions of India was noted, with Sarpang accounting for more cases than any other district from 2009 onward. Case detection was almost exclusively passive, but flexibility and sensitivity was demonstrated by the later addition of active case detection and specification of imported and locally acquired cases. Spatial data were limited to the village level, not allowing identification of transmission hotspots. For blood smears, statistical measures of between-reader agreement and predictive value were not computed. Blood smear quality was suboptimal by at least one criterion in over half of evaluated smears. Timeliness in reporting of cases was on a weekly to monthly basis, and did not meet the WHO goal of immediate notification.

CONCLUSION

As of 2012, the national malaria surveillance system demonstrated flexibility, representativeness, simplicity, and stability. The full potential for data analysis was not yet realized. Attaining the goal of malaria elimination will require system function enhancement through increased and more accurate case detection and rapid investigation, improved health worker training and accountability, focally targeted response measures, and, in particular, the challenge of finding re-introductions of infections from India. Many such measure have been undertaken or planned as part of the next phase of the Bhutan's National Strategic Plan.

Malaria elimination in India and regional implications

The malaria situation in India is complex as a result of diverse socio-environmental conditions. India contributes a substantial burden of malaria outside sub-Saharan Africa, with the third highest Plasmodium vivax prevalence in the world. Successful malaria control in India is likely to enhance malaria elimination efforts in the region. Despite modest gains, there are many challenges for malaria elimination in India, including: varied patterns of malaria transmission in different parts of the country demanding area-specific control measures; intense malaria transmission fuelled by favourable climatic and environment factors; varying degrees of insecticide resistance of vectors; antimalarial drug resistance; a weak surveillance system; and poor national coordination of state programmes. Prevention and protection against malaria are low as a result of a weak health-care system, as well as financial and socioeconomic constraints. Additionally, the open borders of India provide a potential route of entry for artesunate-resistant parasites from southeast Asia. This situation calls for urgent dialogue around tackling malaria across borders-between India's states and neighbouring countries-through sharing of information and coordinated control and preventive measures, if we are to achieve the aim of malaria elimination in the region.

Where chloroquine still works: the genetic make-up and susceptibility of Plasmodium vivax to chloroquine plus primaquine in Bhutan

BACKGROUND

Bhutan has made substantial progress in reducing malaria incidence. The national guidelines recommend chloroquine (CQ) and primaquine (PQ) for radical cure of uncomplicated Plasmodium vivax, but the local efficacy has not been assessed. The impact of cases imported from India on the genetic make-up of the local vivax populations is currently unknown.

METHODS

Patients over 4 years of age with uncomplicated P. vivax mono-infection were enrolled into a clinical efficacy study and molecular survey. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. On day 28 a 14-day regimen of PQ (0.25 mg/kg/day) was commenced under direct observation. After day 42, patients were followed up monthly for a year. The primary and secondary endpoints were risk of treatment failure at day 28 and at 1 year. Parasite genotyping was undertaken at nine tandem repeat markers, and standard population genetic metrics were applied to examine population diversity and structure in infections thought to be acquired inside or outside of Bhutan.

RESULTS

A total of 24 patients were enrolled in the clinical study between April 2013 and October 2015. Eight patients (33.3 %) were lost to follow-up in the first 6 months and another eight patients lost between 6 and 12 months. No (0/24) treatment failures occurred by day 28 and no (0/8) parasitaemia was detected following PQ treatment. Some 95.8 % (23/24) of patients were aparasitaemic by day 2. There were no haemolytic or serious events. Genotyping was undertaken on parasites from 12 autochthonous cases and 16 suspected imported cases. Diversity was high (H E 0.87 and 0.90) in both populations. There was no notable differentiation between the autochthonous and imported populations.

CONCLUSIONS

CQ and PQ remains effective for radical cure of P. vivax in Bhutan. The genetic analyses indicate that imported infections are sustaining the local vivax population, with concomitant risk of introducing drug-resistant strains.

Development and evaluation of a spatial decision support system for malaria elimination in Bhutan

Background

Bhutan has reduced its malaria incidence significantly in the last 5 years, and is aiming for malaria elimination by 2016. To assist with the management of the Bhutanese malaria elimination programme a spatial decision support system (SDSS) was developed. The current study aims to describe SDSS development and evaluate SDSS utility and acceptability through informant interviews.

Methods

The SDSS was developed based on the open-source Quantum geographical information system (QGIS) and piloted to support the distribution of long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) in the two sub-districts of Samdrup Jongkhar District. It was subsequently used to support reactive case detection (RACD) in the two sub-districts of Samdrup Jongkhar and two additional sub-districts in Sarpang District. Interviews were conducted to ascertain perceptions on utility and acceptability of 11 informants using the SDSS, including programme and district managers, and field workers.

Results

A total of 1502 households with a population of 7165 were enumerated in the four sub-districts, and a total of 3491 LLINs were distributed with one LLIN per 1.7 persons. A total of 279 households representing 728 residents were involved with RACD. Informants considered that the SDSS was an improvement on previous methods for organizing LLIN distribution, IRS and RACD, and could be easily integrated into routine malaria and other vector-borne disease surveillance systems. Informants identified some challenges at the programme and field level, including the need for more skilled personnel to manage the SDSS, and more training to improve the effectiveness of SDSS implementation and use of hardware.

Conclusions

The SDSS was well accepted and informants expected its use to be extended to other malaria reporting districts and other vector-borne diseases. Challenges associated with efficient SDSS use included adequate skills and knowledge, access to training and support, and availability of hardware including computers and global positioning system receivers.

Spatio-Temporal Variations in Malaria Incidence in Children Less than 10 Years Old, Health District of Sokone, Senegal, 2010-2013

Introduction

Malaria is a leading cause of morbidity and mortality in sub-Saharan Africa. Detailed characterization of the risks for malaria, among populations living in areas where the disease is endemic, is an important priority, especially for planning and evaluating future malaria-control tools. A prospective cohort study was implemented in children under ten years living in rural areas with high Plasmodium falciparum transmission in Senegal.

Methods

Malaria incidence was prospectively evaluated over three year follow-up among a cohort of children aged less than 10 years old living in eight villages of the Sokone health district. The parents of 1316 children comprising a passive case detection cohort were encouraged to seek care from the study health centers at any time their child felt sick. In the event of reported history of fever within 24 hours or measured axillary temperature ≥ 37.5°C, a Rapid Diagnostic Test (RDT) was performed.

Results

From November 2010 to October 2013, among the 1468 reported febrile episodes, 264 were confirmed malaria episodes. Over the 3 years, 218 (16.9%) children experienced at least one clinical malaria episode. Cumulative malaria incidence was 7.3 episodes per 100 children-year at risk, with remarkably heterogeneous rates from 2.5 to 10.5 episodes per 100 children-year at risk. Clinical malaria prevalence ranged from 11.5 to 28.4% in the high transmission season versus from 9.6 to 21.2% in the low transmission season.

Conclusion

This longitudinal community-based study shows that occurrence of clinical malaria was not evenly distributed among all the cohort children in the eight villages. It demonstrates the complexity of spatial distribution of malaria incidence at a local level, even in a region of vegetation and altitudinal homogeneity.

Epidemiological and Molecular Characterization of Dengue Virus Circulating in Bhutan, 2013-2014

Dengue is one of the most significant public health problems in tropical and subtropical countries, and is increasingly being detected in traditionally non-endemic areas. In Bhutan, dengue virus (DENV) has only recently been detected and limited information is available. In this study, we analyzed the epidemiological and molecular characteristics of DENV in two southern districts in Bhutan from 2013–2014. During this period, 379 patients were clinically diagnosed with suspected dengue, of whom 119 (31.4%) were positive for DENV infection by NS1 ELISA and/or nested RT-PCR. DENV serotypes 1, 2 and 3 were detected with DENV-1 being predominant. Phylogenetic analysis of DENV-1 using envelope gene demonstrated genotype V, closely related to strains from northern India.

We describe the epidemiological and molecular features of DENV currently circulating in the two southwestern districts of Bhutan, demonstrating a shift in serotype dominance from previous DENV-3 (2004–2006) to current DENV-1 (2013–2014). The presence of the dengue virus in Bhutan is a relatively recent one. Unfortunately, dengue epidemiological and molecular data in this country is scarce. A fever outbreak in 2013 and 2014 saw patients seeking care at medical facilities in two district of southwestern Bhutan bordering with India. Analyses of serum specimens collected from these patients indicated that dengue virus was at least a major source of this outbreak. These specimens were analyzed in the Public Health Laboratory in Bhutan and in AFRIMS, Thailand. With a combination of three different assays, we established that 31% of all cases captured were caused by dengue virus, although the proportion was higher in 2013 than in 2014. Three different serotypes of dengue virus were found: DENV-1, -2 and -3. No DENV-4 was found. We successfully isolated DENV-1, from which was sequenced the E gene for further analyses. Our analyses revealed that the current DENV-1 in Bhutan probably originated from India.

Cross-border malaria: a major obstacle for malaria elimination

Movement of malaria across international borders poses a major obstacle to achieving malaria elimination in the 34 countries that have committed to this goal. In border areas, malaria prevalence is often higher than in other areas due to lower access to health services, treatment-seeking behaviour of marginalized populations that typically inhabit border areas, difficulties in deploying prevention programmes to hard-to-reach communities, often in difficult terrain, and constant movement of people across porous national boundaries. Malaria elimination in border areas will be challenging and key to addressing the challenges is strengthening of surveillance activities for rapid identification of any importation or reintroduction of malaria. This could involve taking advantage of technological advances, such as spatial decision support systems, which can be deployed to assist programme managers to carry out preventive and reactive measures, and mobile phone technology, which can be used to capture the movement of people in the border areas and likely sources of malaria importation. Additionally, joint collaboration in the prevention and control of cross-border malaria by neighbouring countries, and reinforcement of early diagnosis and prompt treatment are ways forward in addressing the problem of cross-border malaria.

Plasmodium vivax malaria at households: spatial clustering and risk factors in a low endemicity urban area of the northwestern Peruvian coast

BACKGROUND

Peru has presented a decreasing malaria trend during the last decade, particularly in areas on northwestern coast; however, a limited number of cases continues to be reported yearly mainly in malaria hotspots.

METHODS

A two-phase study was conducted to identify spatial and temporal clusters of incident Plasmodium vivax malaria, as well as to determine risk factors associated with households (HH) presenting P. vivax malaria episodes in an urban area of the northwestern Peruvian Coast from June 2008 to May 2010. In the first stage, a full census of the study population was conducted, including geo-referencing of reported P. vivax episodes. In the second stage, a population-based case-control study allowed the identification of risk factors associated with HHs reporting episodes. A total of 117 case HHs with reported P. vivax and 117 control HHs without malaria episodes were assessed. A semi-structured questionnaire was used to interview the head of households and to collect data on HH location and structure, availability of public services, preventive malaria measures, family member with outdoor occupation (farmer, moto-taxi driver), and other HH characteristics. Univariate and multivariate logistic regression analyses were performed to determine case-HH risk factors. SaTScan was used to detect spatial and temporal P. vivax malaria clusters.

RESULTS

The most likely spatial cluster of malaria incidence included 1,040 people (22.4% of total population) in 245 HHs (24.6% of total HHs) accounting for 283 malaria episodes (40.1% of total episodes) during the study period (RR = 2.3, p < 0.001). A temporal cluster was also identified from April 12, 2009 to July 4, 2009 accounting for 355 malaria episodes (50.4% of total episodes) (RR = 7.2, p = 0.001). Factors significantly associated with case HHs compared with control HHs were: proximity to water drain < 200 metres (OR = 2.3, 95% CI: 1.3, 4.0); HH size >5 individuals (OR = 1.8, 95% CI: 1.0, 3.2); lack of potable water (OR = 1.8, 95% CI: 1.1, 3.2); and having domestic and peridomestic animals (OR = 3.6, 95% CI: 1.3, 9.5).

CONCLUSION

Plasmodium vivax malaria incidence is highly heterogeneous in space and time in the urban study area with important geographical and housing risk factors associated with symptomatic episodes.

Prevalence of asymptomatic malaria and bed net ownership and use in Bhutan, 2013: a country earmarked for malaria elimination

BACKGROUND

With dwindling malaria cases in Bhutan in recent years, the government of Bhutan has made plans for malaria elimination by 2016. This study aimed to determine coverage, use and ownership of LLINs, as well as the prevalence of asymptomatic malaria at a single time-point, in four sub-districts of Bhutan.

METHODS

A cross-sectional study was carried out in August 2013. Structured questionnaires were administered to a single respondent in each household (HH) in four sub-districts. Four members from 25 HH, randomly selected from each sub-district, were tested using rapid diagnostic tests (RDT) for asymptomatic Plasmodium falciparum and Plasmodium vivax infection. Multivariable logistic regression models were used to identify factors associated with LLIN use and maintenance.

RESULTS

All blood samples from 380 participants tested negative for Plasmodium infections. A total of 1,223 HH (92.5% of total HH) were surveyed for LLIN coverage and use. Coverage of LLINs was 99.0% (1,203/1,223 HH). Factors associated with decreased odds of sleeping under a LLIN included: washing LLINs <six months and >nine months compared to washing LLINs every six months; HH in the least poor compared to the most poor socio-economic quintile; a HH income of Nu 5,001-10,000 (US$1 = Nu 59.55), and Nu >10,000, compared to HH with income of <Nu 1,500; HH located one to three hours walking distance to a health centre compared to being located closer to a health centre; a reported lack of knowledge as to what to do in event of LLINs being torn; and keeping LLINs in a box compared to keeping them hanging in the place of use. Factors associated with use of LLINs for purposes other than the intended use included: income group Nu 1,501-3,000 and HH located one to three hours walking distance from a health centre.

CONCLUSIONS

There was high coverage of LLINs in the study area with regular use of LLINs throughout the year. LLIN use for purposes other than malaria prevention was low. With high coverage and regular use of LLINs, and a zero prevalence of malaria infection found in historically high-risk communities during the peak malaria season, it appears Bhutan is on course to achieve malaria elimination.

Further shrinking the malaria map: how can geospatial science help to achieve malaria elimination?

Malaria is one of the biggest contributors to deaths caused by infectious disease. More than 30 countries have planned or started programmes to target malaria elimination, often with explicit support from international donors. The spatial distribution of malaria, at all levels of endemicity, is heterogeneous. Moreover, populations living in low-endemic settings where elimination efforts might be targeted are often spatially heterogeneous. Geospatial methods, therefore, can help design, target, monitor, and assess malaria elimination programmes. Rapid advances in technology and analytical methods have allowed the spatial prediction of malaria risk and the development of spatial decision support systems, which can enhance elimination programmes by enabling accurate and timely resource allocation. However, no framework exists for assessment of geospatial instruments. Research is needed to identify measurable indicators of elimination progress and to quantify the effect of geospatial methods in achievement of elimination outcomes.

The changing epidemiology of malaria elimination: new strategies for new challenges

Malaria-eliminating countries achieved remarkable success in reducing their malaria burdens between 2000 and 2010. As a result, the epidemiology of malaria in these settings has become more complex. Malaria is increasingly imported, caused by Plasmodium vivax in settings outside sub-Saharan Africa, and clustered in small geographical areas or clustered demographically into subpopulations, which are often predominantly adult men, with shared social, behavioural, and geographical risk characteristics. The shift in the populations most at risk of malaria raises important questions for malaria-eliminating countries, since traditional control interventions are likely to be less effective. Approaches to elimination need to be aligned with these changes through the development and adoption of novel strategies and methods. Knowledge of the changing epidemiological trends of malaria in the eliminating countries will ensure improved targeting of interventions to continue to shrink the malaria map.

Effect of bednets and indoor residual spraying on spatio-temporal clustering of malaria in a village in south Ethiopia: a longitudinal study

BACKGROUND

Understanding the spatio-temporal pattern of malaria transmission where prevention and control measures are in place will help to fine-tune strategies. The objective of this study was to assess the effect of mass distribution of bednets and indoor residual spraying (IRS) with insecticides on the spatio-temporal clustering of malaria in one malaria endemic village in south Ethiopia.

METHODS

A longitudinal study was conducted from April 2009 to April 2011. The average population was 6631 in 1346 locations. We used active and passive searches for malaria cases for 101 weeks. SatScan v9.1.1 was used to identify statistically significant retrospective space-time clusters. A discrete Poisson based model was applied with the aim of identifying areas with high rates. PASW Statistics 18 was used to build generalized Poisson loglinear model.

RESULTS

The total number of both types of malaria episodes was 622, giving 45.1 episodes per 1000 persons per year; among these, episodes of Plasmodium falciparum and vivax infection numbered 316 (22.9 per 1000 per year) and 306 (22.2 per 1000 per year), respectively. IRS with Dichlorodiphenyltrichloroethane (DDT) and later with Deltamethrin and free mass distribution of insecticide-treated nets (ITNs) were carried out during the study period. There was space-time clustering of malaria episodes at a household level. The spatio-temporal clustering of malaria was not influenced by free mass distribution of ITNs; however, the time-span of the spatio-temporal clustering of malaria cases ended after IRS with Deltamethrin. The presence of clusters on the south-east edge of the village was consistent with the finding of an increasing risk of acquiring malaria infection for individuals who lived closer to the identified vector breeding site.

CONCLUSION

The risk of getting malaria infection varied significantly within one village. Free mass distribution of ITNs did not influence the spatio-temporal clustering of malaria, but IRS might have eliminated malaria clustering.

Malaria control in Bhutan: case study of a country embarking on elimination

Background

Bhutan has achieved a major reduction in malaria incidence amid multiple challenges. This case study seeks to characterize the Bhutan malaria control programme over the last 10 years.

Methods

A review of the malaria epidemiology, control strategies, and elimination strategies employed in Bhutan was carried out through a literature review of peer-reviewed and grey national and international literature with the addition of reviewing the surveillance and vector control records of the Bhutan Vector-Borne Disease Control Programme (VDCP). Data triangulation was used to identify trends in epidemiology and key strategies and interventions through analysis of the VDCP surveillance and programme records and the literature review. Enabling and challenging factors were identified through analysis of socio-economic and health indicators, corroborated through a review of national and international reports and peer-review articles.

Findings

Confirmed malaria cases in Bhutan declined by 98.7% from 1994 to 2010. The majority of indigenous cases were due to Plasmodium vivax (59.9%) and adult males are most at-risk of malaria. Imported cases, or those in foreign nationals, varied over the years, reaching 21.8% of all confirmed cases in 2006.

Strategies implemented by the VDCP are likely to be related to the decline in cases over the last 10 years. Access to malaria diagnosis in treatment was expanded throughout the country and evidence-based case management, including the introduction of artemisinin-based combination therapy (ACT) for P. falciparum, increasing coverage of high risk areas with Indoor Residual Spraying, insecticide-treated bed nets, and long-lasting insecticidal nets are likely to have contributed to the decline alongside enabling factors such as economic development and increasing access to health services.

Conclusion

Bhutan has made significant strides towards elimination and has adopted a goal of national elimination. A major challenge in the future will be prevention and management of imported malaria infections from neighbouring Indian states. Bhutan plans to implement screening at border points to prevent importation of malaria and to targeted prevention and surveillance efforts towards at-risk Bhutanese and migrant workers in construction sites.