Spatio-temporal trends of malaria incidence from 2011 to 2017 and environmental predictors of malaria transmission in Myanmar
Infectious Diseases of Poverty volume 12, Article number: 2 (2023)
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.
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.
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).
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.
The global malaria burden is still enormous; there were about 247 million cases in 2021, of which 5.4 million cases occurred in Southeast Asia . The Greater Mekong Subregion (GMS) consists of Cambodia, Laos, Myanmar, Thailand, Vietnam, and two provinces of China (Yunan and Guangxi), where Myanmar (87.7%) accounted for most of the indigenous cases of malaria in 2021 . Myanmar accounted for most of the Plasmodium falciparum malaria cases within the region. As falciparum malaria cases declined over time, P. vivax has become the dominant species . Despite the overall decline of malaria incidence over the last several years, Myanmar still reports the majority of malaria cases and deaths in the GMS [3,4,5].
The Mekong Malaria Elimination (MME) program is an initiative to support the GMS countries in achieving the goal of malaria elimination by 2030. However, there are several potential barriers to eliminating malaria in the GMS. First, the emergence and spread of parasites resistant to antimalarial drugs and mosquitoes resistant to insecticides hinder the elimination course. Previous reports of asymptomatic P. falciparum and P. vivax isolates carrying genes potentially associated with drug resistance suggest a silent spread of drug-resistant parasites in Myanmar [6,7,8]. The vector species display tremendous spatial heterogeneity in distribution in the GMS, many of which have developed resistance to insecticides used in public health interventions [9,10,11]. Second, there is a trend for higher levels of transmission and malaria-related mortality near international borders, making ‘border malaria’ a concern for malaria prevention. Third, because of internal military conflicts in Myanmar, refugees and internally displaced people rushed to and settled down along the international borders, increasing the risk of infectious diseases. Border regions frequently exhibit cultural and linguistic heterogeneity, complicating healthcare education, disease prevention, and administration.
Monitoring border malaria is important, but an overall analysis at the national level is essential for governments or organizations to rationally allocate material and human resources for malaria control. By 2016, 291 out of 330 townships in Myanmar were endemic, with more than 40 million people at risk for malaria . According to the data on malaria cases from October to December 2021, P. falciparum and P. vivax cases were distributed in the western, southern and northeast regions of Myanmar, whereas the central regions were essentially devoid of malaria . This spatial heterogeneity, or uneven distribution of malaria cases across the landscape, occurs at multiple spatial scales. Some nations have a greater burden than others, some states or provinces within nations have greater burdens than others, and communities likewise vary in their malaria burdens. Spatial analyses that describe and help explain the spatial and spatio-temporal patterns of malaria are important for public health surveillance, implementation of interventions, and for general policies [13,14,15]. For monthly cases, May to August each year is the peak season of malaria infection in Myanmar . However, the temporal trends of P. falciparum and P. vivax malaria in Myanmar are unclear.
Thus, the present study explored the spatiotemporal dynamics of malaria incidence in Myanmar, at both township and state/region levels, to obtain a granular view of the heterogeneous malaria distribution, detect malaria clusters for targeted control, and identify environmental predictors of P. falciparum and P. vivax malaria. Such information will be essential to informing regional malaria elimination efforts.
Myanmar lies between latitudes 9° and 29° N and longitudes 92° and 102° E; has a total area of 678,500 km2; and is bordered by China, Laos, Thailand, Bangladesh and India. The nation is divided into five physiographic regions: the northern mountains, the western ranges, the eastern plateau, the central basin and lowlands, and the coastal plains. There are three seasons: the cool season (late October to mid-February), the hot-dry season (mid-February to mid-May), and the rainy season (mid-May to late October). Agriculture, forestry, and fishing constitute the most significant contributors to Myanmar’s economy. Annual rainfall in the delta region is approximately 2500 mm, while the average annual rainfall in the dry zone in central Myanmar is less than 1000 mm . Myanmar has seven states and seven regions, including 63 districts and 330 townships, with a total population of about 50,279,900 people in the 2014 census.
Data for this study were obtained from the National Malaria Control Program (NMCP) and other partners such as the Myanmar Medical Association, Medical Action Myanmar, Myanmar Council of Churches, and Myanmar Health Assistant Association. The data include yearly counts of malaria cases and tests at the township level from 2011 to 2017 and monthly malaria cases at the state or regional level from 2011 to 2016. Population counts likewise came from the NMCP. Townships are smaller administrative units than states and regions. Areas defined as either states or regions have been classified based on political and socio-cultural factors but at the same administrative level.
Environmental predictor variables for vegetation and surface flooding were derived from Moderate Resolution Imaging Spectroradiometer (MODIS) products (MOD13Q1/MYD13Q1 250 m AQUA/TERRA 16-day composites). Since many infectious diseases, especially vector-borne diseases, are strongly influenced by environmental factors, we hypothesized that indicators of vegetation and surface flooding would correlate with malaria cases. Three environmental indices were downloaded and considered in these analyses: the normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), and a normalized flooding index (NFI) . NFI is indicative of surface water, NDVI is indicative of green surface vegetation, and EVI is an improved measure of green vegetation that is intended to account for dense forest canopies and atmospheric conditions that can lead to errors in NDVI measurements. Data were downloaded for each of these environmental indices (EI) within each 16-day period from 2011 to 2017. NDVI and EVI were strongly colinear, so we retained only EVI as an indicator of vegetation in our models. These variables were summarized at the township level by calculating the mean NFI and mean EVI for each year and for each township.
Elevation data were also downloaded from the Shuttle Radar Topography Mission 30 m dataset (https://srtm.csi.cgiar.org/srtmdata/), accessed through DIVA-GIS (https://www.diva-gis.org/gdata). Summary statistics (mean elevation and variance in elevation) were calculated for each administrative unit using QGIS 3.10.1 software (Open Source Geospatial Foundation, https://www.qgis.org).
Dynamics of malaria incidence and test positivity
Two primary epidemiological metrics were used in this research: the annual parasite incidence (API, also referred to as the case incidence) and the test positivity. API is the ratio of the number of symptomatic clinical cases of malaria, microscopically or rapid diagnostic test (RDT) confirmed, in a population in a given year to the total population of the region in that year (reported as the number of cases per 1000 people). Test positivity is the number of positive cases divided by the total number of tests. We used both in this research because testing intensity can influence case incidence. For example, some places with little-to-no testing may have high malaria burdens that are not obvious when mapping or otherwise presenting case incidence alone. Malaria cases from the two most prevalent human Plasmodium parasites (P. falciparum and P. vivax) were also shown separately.
We plotted API and test positivity by species, across years for the entire nation, monthly for the entire nation, and annually by state or region. National and regional malaria API and test positivity for P. falciparum and P. vivax were calculated by year and graphed to show annual fluctuations. Species incidences for each month during 2011–2016 were calculated to observe seasonal fluctuations in malaria transmission .
We also mapped estimated crude API, smoothed API, and test positivity at the township level. Empirical Bayes smoothing (SEB) was used with the API estimates to account for potential variance instability from differences in population estimates across geographic units (which can lead to spurious outliers that appear to be hotspots or coldspots of disease). SEB smoothing improves the ability to identify overarching spatial patterns in diseases and other phenomena [19, 20]. These smoothed APIs were calculated using GeoDa 1.14.0 software (Open Source, https://geodacenter.github.io), all maps were generated using QGIS version 3.10.1 (https://www.qgis.org).
Spatial patterns in API and test positivity
We used three approaches for analyzing spatial autocorrelation of API and test positivity: spatial correlograms, the Moran’s I statistic (both global and local), and scan statistics. Spatial correlograms illustrate the magnitude of spatial autocorrelation that exists (on the y-axis, with 1 indicating perfect clustering, − 1 indicating perfect dispersal, and 0 indicating no spatial pattern) between pairs of administrative units at different distances away from each other (distances along the x-axis).
We then tested for global and local spatial autocorrelation using Moran’s I statistic and local indicator of spatial autocorrelation (LISA). Moran’s I ranges from 1 to − 1, with a score of zero indicating the null hypothesis of spatial randomness. Positive values indicate clusters of malaria cases, while negative values indicate that neighboring areas are characterized by different malaria cases (i.e. areas with high cases neighboring areas with low cases). The LISA statistics were investigated and mapped to identify four types of clusters (high-high, low-low, high-low, and low-high) of smoothed P. falciparum and P. vivax incidence at the township level. High-high and low-low clusters present hotspots and coldspots, respectively. High-low and low-high categories represent outliers . The statistical significance was tested using 999 Monte Carlo permutations, and a P-value of 0.05 or less was considered statistically significant . Both Moran’s I and LISA statistics are based on adjacency matrices, and the spatial weight matrix is based on a Queen adjacency matrix, which establishes connections between all neighbors that share a common point or length on the boundary. Moran’s I and LISA statistics were calculated using GeoDa 1.14.0 (https://geodacenter.github.io).
Finally, we used Kulldorf’s retrospective space-time scan statistics to test for likely clusters of cases in space and time . Briefly, the process involves scanning data windows across space and time and recording the number of observations and expectations within the window for each location in comparison to expected observations outside of the window. The risk of malaria within and outside the window was tested using a likelihood ratio test with the null hypothesis of equal risk across space. The window with the largest log likelihood ratio (LLR,the most likely cluster) was considered the cluster with the highest malaria risk. The window with the next to maximum LLR represented the secondary likely cluster and was considered the area with the second highest malaria risk. The most likely cluster (primary cluster) was identified based on the maximum log likelihood ratio, other clusters with statistically significant log likelihood values were defined as secondary clusters . We used the discrete space-time Poisson model to look for malaria clusters in space and time, with 1-month temporal aggregations. The scan statistics were calculated using SaTScan™ version 9.3 (Kuldorff M. and Information Management Services, Inc. https://www.satscan.org/).
Negative binomial regressions for predictors of reported P. falciparum and P. vivax malaria cases
We used negative binomial generalized additive models to look for predictors of malaria cases at the township level and annually for P. falciparum and P. vivax malaria, respectively. The models control for repeated measures within administrative units using a random intercept, and for the relative geographic locations of the administrative units using a smoothed spline function interaction term for the geographic locations of the mean centers of the administrative units. Rather than assuming the shape of the potential association between a given variable and the outcome (counts of malaria cases), smoothing splines were fit to the continuous covariates. Covariates in the models included: mean elevation, variance in elevation, mean annual EVI, mean annual NFI, the number of exams in a given location, and the year. Regression analyses were performed separately for falciparum and vivax malaria. The model outputs are presented as plots of the spline functions. The generalized additive negative binomial regressions were calculated using the statistical software R 4.0 (R Core Team, R Foundation for statistical computing, Vienna, Austria).
Trends of countrywide malaria incidence
From 2011 to 2017, a total of 1,426,737 malaria cases were reported in Myanmar. Malaria API and test positivity in Myanmar showed an overall decline, especially since 2012 (Fig. 1A). From 2012 to 2016, P. falciparum cases had a relatively sharp decline, with a 9-fold reduction in API from 8.5 per 1000 in 2012 to 0.9 per 1000 in 2016. Test positivity was similarly reduced from 32.6% in 2012 to 0.7% in 2016. In contrast, the decrease in P. vivax cases was gradual, with a 10-fold reduction in API from 3.0 per 1000 in 2012 to 0.3 per 1000 in 2016 and test positivity decreasing from 8.8% in 2012 to 1.6% in 2016. Compared to 2016, both P. falciparum and P. vivax API experienced a rebound in 2017, although test positivity for both species continually declined (Fig. 1A). The annual API and test positivity in individual regions and states were generally consistent with the overall trend (Fig. 1B, C), while the incidence in the western border state Chin and southern border state Kayin showed a significant upward trend in 2016‒2017, responsible for the overall increase of malaria API in Myanmar in 2016‒2017 (Fig. 1B). Monthly incidence was generally high in the rainy season (May‒September), with the peak occurring in June or July each year (Fig. 1D). Some areas have a second peak in the cool/dry season (November–January, especially Rakhine and Chin states) (Fig. 2). It is noteworthy that Chin, Rakhine, and Kachin states ranked as the top malaria burden states in all years except 2012 (Fig. 2).
Township-level and region-level patterns of malaria incidence
At the township level, P. falciparum and P. vivax malaria exhibited similar spatial patterns. Correlograms of township P. falciparum and P. vivax malaria incidence showed significant spatial autocorrelation up to 100 km (Fig. 3), indicating that malaria incidence in Myanmar was not randomly distributed but rather occurred as clusters among adjacent townships. Also, in the autocorrelation analysis of the monthly incidence from 2011 to 2016 (based on region-level data), we found that the incidence was positively correlated within the range of 100 km except in 2014, further corroborating the above result (Additional file 1: Figs. S1, S2, S3).
Spatial distribution and cluster of P. falciparum and P. vivax incidence
Analysis of the smoothed APIs showed that areas with higher incidence of P. falciparum malaria were concentrated along international borders (Fig. 4). Despite the substantial reduction in P. falciparum malaria incidence over the study period, the distribution pattern remained. The number of townships reporting no P. falciparum malaria increased from 58 (17.6%) in 2011 to 69 (20.9%) in 2016, but decreased to 52 (15.8%) in 2017 (Fig. 4; Table 1). The number of townships with an API of 0‒1 per 1000 increased over the study period. In 2016 and 2017, the number of townships with P. falciparum API of 0‒1 per 1000 was 100 (30.3%) and 118 (35.8%), respectively (Fig. 4; Table 1). The number of townships with an API of > 5 per 1000 sharply decreased from 231 (70.0%) in 2011 to 71 (21.5%) in 2017 (Fig. 4; Table 1).
From the LISA statistics, clustered high API townships (red) were confined to the northern and western parts of Chin State, Kachin State, and Saging Region, while the clustered low API townships (blue) were gathered in the central and southern parts. The spatial autocorrelation (Global Moran’s I) ranked between 0.105 and 0.502 (P < 0.016) during 2011‒2014 and 2016, revealing a significant and positive spatial autocorrelation in Myanmar (Fig. 4, also Additional file 1: Figs. S4, S5). Over the seven years, the largest malaria clusters were concentrated in western areas of Myanmar. The number of townships covered by the primary clusters decreased from 19 in 2011 to 1 in 2017. Secondary clusters were characterized in the northwest areas (Fig. 4; Table 2).
The P. vivax high API areas were also concentrated along international borders, and the number declined from 2012 to 2017. The number of townships reporting no P. vivax malaria cases decreased from 69 (20.9%) in 2011 to 44 (13.3%) in 2017 (Fig. 5; Table 1). The number of townships with smoothed P. vivax API of 0‒1 per 1000 increased from 97 (29.4%) to 130 (39.4%) (Fig. 5; Table 1). In 2016 and 2017, the number of townships with a P. vivax API of > 5 per 1000 was 48 (14.6%) and 62 (18.8%), respectively (Fig. 5; Table 1). The clustered high API townships (red) were confined to the western, northern, and southeastern parts of Chin, Kachin, and Kayah states. For P. falciparum, clustered low API townships (blue) were gathered in central and southern parts. The spatial autocorrelation ranked between 0.083 and 0.221 (P < 0.035) during 2011‒2016, revealing a significant and positive spatial autocorrelation in Myanmar during this period. In 2017, there was no significantly spatial autocorrelation (Fig. 5, also Additional file 1: Figs. S6, S7). Unlike P. falciparum malaria, primary clusters of P. vivax malaria changed over the seven years. The number of townships covered by the primary clusters decreased from 37 in 2013 to 1 in 2017. Secondary clusters were characterized in the northern areas in the last two years (Fig. 5; Table 3).
The temporal clustering of townships per year was analyzed to identify periods with a higher than expected number of cases, controlling for the population in the respective administrative units. For each stratum, a temporal cluster was identified if there was the same seasonal pattern of high malaria incidence between states/regions in a given year. The temporal cluster decreased in length during the study period (June‒November in 2011 to June‒July in 2017) (Table 4). The temporal cluster analysis also showed that malaria incidence was concentrated in the rainy months, ranging from June to November. The range of months with higher temporal clustering varied slightly from year to year. In 2016, malaria risk was higher in June‒July. From January 2011 to October 2013, there was a temporal cluster of malaria cases in the study area (LLR = 207,689.62, P = 0.001). A total of 1,001,309 malaria cases occurred during the same period, and the risk of malaria infection was 3.27 times that of other periods (Table 4).
Spatio-temporal clusters of malaria
The space-time Poisson model showed two spatio-temporal malaria clusters from 2011 to 2016 (Table 5). The primary cluster was detected from January 2011 to December 2013. The primary cluster covered two states (Sagaing and Kachin, Fig. 5) and persisted from 2011 to 2016.
Relationships between reported malaria cases and environmental variables
Potential associations between variables in our negative binomial regression and malaria cases are assessed through plots of the spline functions for each variable (Fig. 6). Annual reported P. falciparum cases were highest in townships with a mean elevation of approximately 500–600 m (Fig. 6) and high variance (i.e., townships with both high and low elevations). There was an association between reported vivax malaria cases and elevation as well, but the peak association was at a slightly higher elevation (700–800 m). There was an apparent linear relationship between mean EVI and annual P. falciparum cases (Fig. 6), but not with P. vivax cases at the township level. There was a decrease in overall P. falciparum and P. vivax cases over the years, though the decrease was steeper for falciparum malaria and curvilinear for vivax malaria [first increasing until 2013 and then decreasing each year after (Fig. 6)].
Our findings showed that malaria incidence reduced in Myanmar from 2011 to 2017. There was significant spatial autocorrelation in malaria cases and test positivity across the nation. Areas with higher malaria incidence were concentrated along international borders. Primary clusters of P. falciparum malaria persisted in western townships, while clusters of P. vivax shifted geographically over the study period. Malaria cases were higher in townships with a mean elevation of 500‒600 m and falciparum cases were associated with higher levels of vegetation.
In this study, we used an array of analytic approaches to explore spatial and temporal patterns of malaria API and test positivity based on P. falciparum and P. vivax data aggregated at the township or state/region level from 2011 to 2017 in Myanmar. The NMCP categorizes geographic locations based on reported case incidence levels (API), with three overarching categories that have corresponding public health approaches: less than 1 per 1000 people per year (pre-elimination), 1 to less than 10 per 1000 people per year (moderate transmission), and above 10 per 1000 people per year (high transmission) . Our findings showed that more than 50% of townships reached malaria pre-elimination by 2017.
The reported P. falciparum API maintained a continuous decline throughout the study period, which may be mainly related to increased access to diagnoses and treatment with artemisinin-based combination therapies (ACT) . However, the reported vivax malaria API showed a more gradual decline, probably associated with the intrinsic biological features of P. vivax, including relapses induced by hypnozoites , missed diagnosis because of the low density of infection, lower accuracy from widely used rapid diagnostic tests , the early production of gametocytes favoring continuous transmission, and the high proportion of asymptomatic infections .
Some patterns with vivax malaria may be related to the increased roll-out of rapid diagnostic tests that detect both falciparum and vivax malaria (the first rapid diagnostic tests were only capable of detecting falciparum malaria). This roll-out corresponds to the initial increase in vivax malaria cases up to 2013 when they began to decrease (apparent in Fig. 1) and could be related to an increased capacity for diagnosing vivax malaria. A few studies have also shown that vivax malaria relapses following the treatment of falciparum malaria [30, 31]. It is, therefore, possible that widespread increased treatment of falciparum malaria partially drove the short increase of reported vivax malaria cases prior to their decline in 2013. Complete cure of vivax malaria is hindered by the risk of hemolysis in glucose-6-phosphate dehydrogenase deficient (G6PDd) individuals post primaquine and other 8-aminoquinoline antimalarials. This region has a high prevalence of G6PD deficiency along the China-Myanmar border (16.9%) and Thailand–Myanmar border (13.7%) [32, 33], as well as in western and central Myanmar (10% and 6.8%, respectively) [34, 35]. The need for primaquine to radically cure vivax malaria, coupled with the risk of treating G6PDd patients and the difficulties in diagnosing G6PD deficiency in field settings, remains a challenge for adequately addressing this malaria species.
The apparent declines in the API and test positivity of P. falciparum and P. vivax malaria from 2011 to 2017 in Myanmar may be illustrative of the impact of combined efforts by governmental and non-governmental organizations to eliminate malaria in the nation [36, 37]. Previous studies showed significant associations between declining reported cases and the following factors: long-lasting insecticidal net/insecticide-treated bed net distribution, indoor residual spraying, the concentration of village health workers, amount of health worker training, development of volunteers, socioeconomic status, and improved ACT availability [36, 38,39,40,41,42,43].
Our analyses also showed malaria to peak and cluster temporally from approximately May and into November. This corresponds with the rainy season in Myanmar, and the pattern has been well-described in several other studies [14, 44]. A second peak was apparent in some high burden states/regions (especially Rakhine and Chin states). This pattern has likewise been described in other locations but has not been fully explained. It may correspond to differences in human-mosquito exposure that are related to seasonal agriculture-the second peak often corresponds with rice harvesting season [45, 46], with changes in mosquito abundance , or other factors such as decreased mosquito net use after the rains have ceased .
Malaria exhibited significant spatial clustering during the study period as well. We found clustering of P. falciparum in the western and northern parts of Chin State, Kachin State, and Sagaing Region. From the Poisson model, we found that primary and secondary clusters of P. falciparum also persisted in the northern and western regions of Myanmar. Although the sizes and locations of P. falciparum clusters became gradually smaller over time, the cluster locations were relatively stable in the western and northern regions [36, 49]. A recent study found P. falciparum was the predominant species accounting for more than 80% of all infections in Paletwa Township of Chin State . Conversely, primary and secondary clusters of P. vivax malaria changed over the seven-year study period. The maps of incidence and test positivity demonstrate that high-burden areas of P. vivax malaria tend to migrate west to east. In the last three years of the study period, P. vivax clustered in eastern Kachin State bordering China and southern Shan State and Kayah State bordering Thailand. Previous studies showed P. vivax is the predominant species along many of these international borders [29, 50, 51].
Spatial patterns in malaria can be largely driven by environmental factors that vary across landscapes [52, 53]. Forested areas have long been associated with malaria in Southeast Asia [54, 55]. The hilly and mountainous areas along the international borders have had less economic development, at least partially, as a result of long-standing conflict in these areas. Changes in environmental attributes (such as forest cover) and other socioeconomic factors can lead to changes in the burden of malaria . In this study, we did not limit our measures of EVI or NFI to specific months. Instead, we took the mean value for each year of the study period. If there is a positive association in the model, the interpretation would then be that a higher mean annual value of EVI in a given year is associated with a higher number of malaria cases, etc. It is true that most cases occur during the rainy season, however in several locations there are double peaks—with the second peak occurring outside of the normal rainy season. Partially for this reason, and likewise because the model is analyzing data that have been aggregated yearly, we prefer to have environmental measures that include but are not limited to the rainy season.
Our analyses showed that malaria cases were highest in townships with a mid-level elevation (mean elevation of approximately 500‒600 m) and that falciparum cases were associated with high levels of vegetation (measured using EVI). The latter corresponds to the well-known macro associations between forests and falciparum malaria. While malaria elimination efforts have increased in the last decades in Myanmar, deforestation has also increased , and both may impact malaria transmission.
This study has a few limitations. The data might not be accurate for some years since the data were mainly derived from government health centers and the program’s village malaria volunteers. With regard to the temporal trends, it is noteworthy that the coverage of health care facilities has improved drastically in Myanmar over the last several years. For example, many community-based health clinics were set up in Kayin State (beginning around 2014 and 2015). The increase in malaria diagnosis and surveillance at first might give an impression of increased cases when in fact it is the result of increased diagnosis. We control for the influence of testing through the test positivity metric and by including the number of tests in our regressions. However, we cannot control for missing data (either through problems with surveillance systems or from a sparsity of clinics in some regions). Likewise, we do not have data on the proportions of tests that were done using RDTs or microscopy and the data are at aggregate levels meaning that we cannot perform detailed analyses of demographic risk factors. Furthermore, since our data and analyses are limited to Myanmar, it is possible that some of the patterns we see are obscured along the edges of the map. It is possible that we are missing hotspots and coldspots of malaria that would be apparent if we had corresponding data from neighboring nations (this is a type of ‘edge effect’).
Another potential limitation is related to the environmental variables included in our models. While we are confident that our model results indicate true associations between environmental variables and either falciparum or vivax malaria cases, other environmental variables could be included in similar analyses. For example, while we chose to use EVI this measure of vegetation would not necessarily differentiate forests from other kinds of vegetation. Since our a priori goal was not to look for associations between forests and malaria, we chose this more inclusive variable instead of measures that would be specific to forests.
Another limitation is the deficiency of monthly data, which were only available at the state/region level from 2011 to 2016. More granular data (i.e., at the village or village tract level) would be superior. Lastly, the coup in 2021 and COVID-19 pandemic (beginning in 2020) have disrupted many malaria control and elimination efforts. From passive case surveillance in Laiza town in Kachin State along the China-Myanmar border, malaria cases declined from 2016 to 2019 but increased rapidly beginning in 2020 . This suggests that local governments or organizations should resume malaria surveillance and implementation of interventions as soon as possible. Meanwhile, the bordering countries should pay attention to the importation of malaria cases.
This study describes spatial and temporal patterns of malaria incidence across the entire nation of Myanmar over a 7-year period (2011–2017). We describe an overall reduction in both falciparum and vivax malaria, which is likely driven at least partially by enhanced malaria elimination and control efforts during this same period. Deforestation and socioeconomic factors may likewise explain part of this reduction. However, malaria continued to cluster in some locations (falciparum malaria in the northwest, vivax malaria in the northeast), and these clusters could act as reservoirs for rebounds in other regions if vigilance is not maintained. Unfortunately, the SARS-CoV-2 (beginning in 2020) pandemic and the military coup in 2021 have disrupted many of these efforts. To achieve elimination, it will be necessary to resume heightened control and elimination efforts, and likely will require tailored approaches to address the complications of eliminating vivax malaria.
Availability of data and materials
All data analyzed for this study are included within the article.
Artemisinin-based combination therapies
Annual parasite incidence
Greater Mekong Subregion
Mekong Malaria Elimination
Moderate Resolution Imaging Spectroradiometer
Normalized difference vegetative index
Enhanced vegetation index
Normalized flooding index
Rapid diagnostic test
Local indicator of spatial autocorrelation
National malaria control program
WHO. World Malaria Report 2022. Geneva: World Health Organization; 2022.
WHO. Mekong Malaria Elimination Programme. Geneva: World Health Organization; 2021.
Inkochasan M, Gopinath D, Vicario E, Lee A, Duigan P. Access to health care for migrants in the Greater Mekong Subregion: policies and legal frameworks and their impact on malaria control in the context of malaria elimination. WHO South East Asia J Public Health. 2019;8:26–34.
Hewitt SE. Let’s ‘cut to the chase’ on malaria elimination in the Greater Mekong Subregion. Trans R Soc Trop Med Hyg. 2019;113:161–2.
Hewitt S, Delacollette C, Poirot E. Malaria control in the Greater Mekong Subregion: an overview of the current response and its limitations. Southeast Asian J Trop Med Public Health. 2013;44(Suppl 1):249–305. discussion 6–7.
Zhao Y, Liu Z, Soe MT, Wang L, Soe TN, Wei H, et al. Genetic variations associated with drug resistance markers in asymptomatic Plasmodium falciparum infections in Myanmar. Genes. 2019;10:692.
Nyunt MH, Shein T, Zaw NN, Han SS, Muh F, Lee SK, et al. Molecular evidence of drug resistance in asymptomatic malaria infections, Myanmar, 2015. Emerg Infect Dis. 2017;23:517–20.
Zhao Y, Wang L, Soe MT, Aung PL, Wei H, Liu Z, et al. Molecular surveillance for drug resistance markers in Plasmodium vivax isolates from symptomatic and asymptomatic infections at the China-Myanmar border. Malar J. 2020;19:281.
Kwansomboon N, Chaumeau V, Kittiphanakun P, Cerqueira D, Corbel V, Chareonviriyaphap T. Vector bionomics and malaria transmission along the Thailand-Myanmar border: a baseline entomological survey. J Vector Ecol. 2017;42:84–93.
Zhong D, Aung PL, Mya MM, Wang X, Qin Q, Soe MT, et al. Community structure and insecticide resistance of malaria vectors in northern-central Myanmar. Parasit Vectors. 2022;15:155.
Chaumeau V, Cerqueira D, Zadrozny J, Kittiphanakun P, Andolina C, Chareonviriyaphap T, et al. Insecticide resistance in malaria vectors along the Thailand-Myanmar border. Parasit Vectors. 2017;10:165.
MoHS. National Strategic Plan: Intensifying Malaria Control and Accelerating Progress Towards Malaria Elimination (2016–2020). Myanmar: Department of Public Health; 2016.
Childs DZ, Cattadori IM, Suwonkerd W, Prajakwong S, Boots M. Spatiotemporal patterns of malaria incidence in northern Thailand. Trans R Soc Trop Med Hyg. 2006;100:623–31.
Hu Y, Zhou G, Ruan Y, Lee MC, Xu X, Deng S, et al. Seasonal dynamics and microgeographical spatial heterogeneity of malaria along the China–Myanmar border. Acta Trop. 2016;157:12–9.
Saita S, Silawan T, Parker DM, Sriwichai P, Phuanukoonnon S, Sudathip P, et al. Spatial heterogeneity and temporal trends in malaria on the Thai-Myanmar border (2012–2017): a retrospective observational study. Trop Med Infect Dis. 2019;4(2):62.
Wikipedia, Myanmar. https://www.enwikipediaorg/wiki/Myanmar. Accessed 18 Aug 2022.
Boschetti M, Nutini F, Manfron G, Brivio PA, Nelson A. Comparative analysis of normalised difference spectral indices derived from MODIS for detecting surface water in flooded rice cropping systems. PLoS ONE. 2014;9:e88741.
Xiao D, Long Y, Wang S, Wu K, Xu D, Li H, et al. Epidemic distribution and variation of Plasmodium falciparum and Plasmodium vivax malaria in Hainan, China during 1995–2008. Am J Trop Med Hyg. 2012;87:646–54.
Chuang TW, Wimberly MC. Remote sensing of climatic anomalies and West Nile virus incidence in the northern Great Plains of the United States. PLoS ONE. 2012;7:e46882.
Marshall RJ. Mapping disease and mortality rates using empirical Bayes estimators. J R Stat Soc Ser C Appl Stat. 1991;40:283–94.
Anselin L. Local indicators of spatial Association—LISA. Geogr Anal. 1995;27:93–115.
Wangdi K, Kaewkungwal J, Singhasivanon P, Silawan T, Lawpoolsri S, White NJ. Spatio-temporal patterns of malaria infection in Bhutan: a country embarking on malaria elimination. Malar J. 2011;10:89.
Kulldorff M. A spatial scan statistic. Commun Stat Theory Methods. 1997;26:1481–96.
Liu Y, Wang X, Liu Y, Sun D, Ding S, Zhang B, et al. Detecting spatial-temporal clusters of HFMD from 2007 to 2011 in Shandong Province, China. PLoS ONE. 2013;8:e63447.
WHO. Strategy for malaria elimination in the Greater Mekong Subregion (2015–2030). Geneva: World Health Organization; 2015.
McLean ARD, Indrasuta C, Khant ZS, Phyo AK, Maung SM, Heaton J, et al. Mass drug administration for the acceleration of malaria elimination in a region of Myanmar with artemisinin-resistant falciparum malaria: a cluster-randomised trial. Lancet Infect Dis. 2021;21:1579–89.
Yuan L, Wang Y, Parker DM, Gupta B, Yang Z, Liu H, et al. Therapeutic responses of Plasmodium vivax malaria to chloroquine and primaquine treatment in northeastern Myanmar. Antimicrob Agents Chemother. 2015;59:1230–5.
Agarwal R, Choi L, Johnson S, Takwoingi Y. Rapid diagnostic tests for Plasmodium vivax malaria in endemic countries. Cochrane Database Syst Rev. 2020;11:CD013218.
Zaw MT, Thant M, Hlaing TM, Aung NZ, Thu M, Phumchuea K, et al. Asymptomatic and sub-microscopic malaria infection in Kayah State, eastern Myanmar. Malar J. 2017;16:138.
Ashley EA, Phyo AP, Carrara VI, Tun KM, Nosten F, Smithuis F, et al. Plasmodium vivax relapse rates following Plasmodium falciparum malaria reflect previous transmission intensity. J Infect Dis. 2019;220:100–4.
Hossain MS, Commons RJ, Douglas NM, Thriemer K, Alemayehu BH, Amaratunga C, et al. The risk of Plasmodium vivax parasitaemia after P. falciparum malaria: an individual patient data meta-analysis from the WorldWide Antimalarial Resistance Network. PLoS Med. 2020;17:e1003393.
Zeng W, Malla P, Xu X, Pi L, Zhao L, He X, et al. Associations among soil-transmitted helminths, G6PD deficiency and asymptomatic malaria parasitemia, and anemia in schoolchildren from a conflict zone of northeast Myanmar. Am J Trop Med Hyg. 2020;102:851–6.
Bancone G, Chu CS, Somsakchaicharoen R, Chowwiwat N, Parker DM, Charunwatthana P, et al. Characterization of G6PD genotypes and phenotypes on the northwestern Thailand-Myanmar border. PLoS ONE. 2014;9:e116063.
Han KT, Han ZY, Aye KH, Wai KT, Thi A, Cui L, et al. G6PD deficiency among malaria-infected national groups at the western part of Myanmar with implications for primaquine use in malaria elimination. Trop Med Health. 2021;49:47.
Oo NN, Bancone G, Maw LZ, Chowwiwat N, Bansil P, Domingo GJ, et al. Validation of G6PD point-of-care tests among healthy volunteers in Yangon, Myanmar. PLoS ONE. 2016;11:e0152304.
Mu TT, Sein AA, Kyi TT, Min M, Aung NM, Anstey NM, et al. Malaria incidence in Myanmar 2005–2014: steady but fragile progress towards elimination. Malar J. 2016;15:503.
Program NMC. National Plan for Malaria Elimination in Myanmar 2016–2030. 2017.
Lu G, Liu Y, Wang J, Li X, Liu X, Beiersmann C, et al. Malaria training for community health workers in the setting of elimination: a qualitative study from China. Malar J. 2018;17:95.
Linn NYY, Kathirvel S, Das M, Thapa B, Rahman MM, Maung TM, et al. Are village health volunteers as good as basic health staffs in providing malaria care? A country wide analysis from Myanmar, 2015. Malar J. 2018;17:242.
Namuganga JF, Epstein A, Nankabirwa JI, Mpimbaza A, Kiggundu M, Sserwanga A, et al. The impact of stopping and starting indoor residual spraying on malaria burden in Uganda. Nat Commun. 2021;12:2635.
Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.
Lengeler C. Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004. https://doi.org/10.1002/14651858.CD000363.pub2.
Nosten F, White NJ. Artemisinin-based combination treatment of falciparum malaria. Am J Trop Med Hyg. 2007;77:181–92.
Khine SK, Kyaw NTT, Thekkur P, Lin Z, Thi A. Malaria hot spot along the foothills of Rakhine state, Myanmar: geospatial distribution of malaria cases in townships targeted for malaria elimination. Trop Med Health. 2019;47:60.
Aung PL, Soe MT, Oo TL, Aung KT, Lin KK, Thi A, et al. Spatiotemporal dynamics of malaria in Banmauk Township, Sagaing region of Northern Myanmar: characteristics, trends, and risk factors. BMC Infect Dis. 2022;22:653.
Aung PL, Soe MT, Oo TL, Khin A, Thi A, Zhao Y, et al. Predictors of malaria rapid diagnostic test positivity in a high burden area of Paletwa Township, Chin State in Western Myanmar. Infect Dis Poverty. 2021;10:6.
Oo TT. Biology and Vector Competence of the Anopheline Mosquitoes of Myanmar with Special Consideration of Anopheles dirus. Combined Faculties for the Natural Sciences and for Mathematics of the Ruperto-Carola University of Heidelberg; 2003.
Moiroux N, Boussari O, Djenontin A, Damien G, Cottrell G, Henry MC, et al. Dry season determinants of malaria disease and net use in Benin, West Africa. PLoS ONE. 2012;7:e30558.
WHO. Malaria country Profiles-Myanmar. Geneva: World Health Organization; 2017.
Huang F, Zhang L, Xue JB, Zhou HN, Thi A, Zhang J, et al. From control to elimination: a spatial-temporal analysis of malaria along the China-Myanmar border. Infect Dis Poverty. 2020;9:158.
Geng J, Malla P, Zhang J, Xu S, Li C, Zhao Y, et al. Increasing trends of malaria in a border area of the Greater Mekong Subregion. Malar J. 2019;18:309.
Stresman GH. Beyond temperature and precipitation: ecological risk factors that modify malaria transmission. Acta Trop. 2010;116:167–72.
Dabaro D, Birhanu Z, Negash A, Hawaria D, Yewhalaw D. Effects of rainfall, temperature and topography on malaria incidence in elimination targeted district of Ethiopia. Malar J. 2021;20:104.
Mercado CEG, Lawpoolsri S, Sudathip P, Kaewkungwal J, Khamsiriwatchara A, Pan-Ngum W, et al. Spatiotemporal epidemiology, environmental correlates, and demography of malaria in Tak Province, Thailand (2012–2015). Malar J. 2019;18:240.
Prothero RM. Malaria, forests and people in Southeast Asia. Singap J Trop Geo. 1999;20:76–85.
Gallalee S, Ward AV, Aye MM, Aung NKZ, Dunn JC, Lavenberg S, et al. Factors associated with the decline of malaria in Myanmar’s Ayeyarwady Region between 2013 and 2017. Sci Rep. 2021;11:20470.
Global Forest Watch. https://www.globalforestwatchorg/. Accessed 18 Aug 2022.
Cui L, Sattabongkot J, Aung PL, Brashear A, Cao Y, Kaewkungwal J, et al. Multidisciplinary investigations of sustained malaria transmission in the Greater Mekong Subregion. Am J Trop Med Hyg. 2022;107:138–51.
We thank officials from the Myanmar National Malaria Control Program and representatives of other malaria partners for granting access to the data.
This research was supported by a Grant (U19AI089672) from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, USA.
Ethics approval and consent to participate
This study was reviewed and approved by the Institutional Review Boards of the Department of Medical Research, Myanmar (Ethics/DMR/2017/077AE/2018E/2019E/2020), and the University of South Florida. This is a secondary analysis of anonymized data; the authors had no contact with the study participants and therefore no informed consent was requested. All methods were carried out in accordance with relevant guidelines and regulations.
Consent for publication
The authors declare that they have no conflict of interests.
Figure S1. Monthly spatial correlograms with test positivity for combined malaria cases/tests. Figure S2. Monthly spatial correlograms with the API (annual parasite incidence) for combined malaria. API is calculated as the number of cases per 1000 people per year. Figure S3. Monthly spatial correlograms with SEB (smoothed empirical Bayesian) API for combined malaria. Figure S4. P. falciparum malaria distribution and spatial autocorrelation statistic using the testpositivity (TP). Figure S5. P. falciparum malaria distribution and spatial autocorrelation statistic using the annual parasite incidence (API), calculated as the number of cases per 1000 people per year. Figure S6. P. vivax malaria distribution and spatial autocorrelation statistic using the test positivity (TP). Figure S7. P. vivax malaria distribution and spatial autocorrelation statistic using the annual parasite incidence (API) calculated as the number of cases per 1000 people per year.
About this article
Cite this article
Zhao, Y., Aung, P.L., Ruan, S. et al. Spatio-temporal trends of malaria incidence from 2011 to 2017 and environmental predictors of malaria transmission in Myanmar. Infect Dis Poverty 12, 2 (2023). https://doi.org/10.1186/s40249-023-01055-6