Abstract
India has made tremendous progress in reducing malaria mortality and morbidity in the last decade. Mizoram State in North-East India is one of the few malaria-endemic regions where malaria transmission has continued to remain high. As Mizoram shares international borders with Bangladesh and Myanmar, malaria control in this region is critical for malaria elimination efforts in all the three countries. For identifying hotspots for targeted intervention, malaria data from 385 public health sub-centers across Mizoram were analyzed in the Geographic Information System. Almost all the sub-centers reporting high Annual Parasite Index (>â10) are located in Mizoramâs districts that border Bangladesh. Getis-Ord Gi* statistic shows most of the sub-centers located along the Bangladesh border in the Lawngtlai and Lunglei districts to be the malaria hotspots. The hotspots also extended into the Mamit and Siaha districts, especially along the borders of Lawngtlai and Lunglei. Analysis of terrain, climatic, and land use/land cover datasets obtained from the Global Modelling and Assimilation Office and satellite images show Mizoramâs western part (Lawngtlai, Lunglei, and Mamit districts) to experience similar topographic and climatic conditions as the bordering Rangamati district in the Chittagong division of Bangladesh. Climatic trends in this region from 1981 to 2021, estimated by the MannâKendall test and Sen's slope estimates, show an increasing trend in minimum temperature, relative humidity, rainfall, and the associated shift of climatic pattern (temperate to tropical monsoon) could facilitate malaria transmission. The quasi-Poisson regression model estimates a strong association (pâ<â0.001) between total malaria cases, temperature range, and elevation. The KruskalâWallis H test shows a statistically significant association between malaria cases and forest classes (pâ<â0.001). A regional coordination and strategic plan are required to eliminate malaria from this hyper-endemic malaria region of North-East India.
Similar content being viewed by others
Introduction
Globally, in 2020, malaria cases and deaths were estimated to be 241 million and 627,000 respectively1. Compared to 2019, this represents a 5.8% increase in cases and an 11% increase in deaths;â~â2/3rd of the additional deaths (47,000) in 2020 have been attributed to the disruption in malaria prevention and control efforts by the COVID-19 pandemic2. Even though the pandemic was predicted to increase malaria incidence and deaths in 2020, the phenomenal progress made since 2000 (27% and 51% decrease in case incidence and mortality, respectively) has started plateauing around 20172.
In India, in the last decade, malaria cases and deaths have declined significantly; 1018 deaths in 2010 have steeply decreased to 93 in 20201. Among the 11 most malaria-endemic countries, India is the only one to report a reduction in cases in 2020, compared to 20191. Despite the steep decrease in malaria incidence across India, there are a few endemic pockets where the decline in malaria incidence has been much lower than the national average and continues to pose a stiff challenge to Indiaâs malaria elimination efforts. The state of Mizoram in North-East (NE) India is one such hyper-endemic malaria region, where malaria remains a major public health challenge.
The Indo-Myanmar region where Mizoram is located in one of the top 10 biodiversity hotspots in the world3. A majority of Mizoramâs population belongs to the Kuki Chin group, with an ancestral link to the Tibeto-Burma branch of the Indo-Mongloid race4. Mizoram is dominated by tribals (~â94%)5, receives heavy rainfall (~â2000 mm annually), and an estimated 89% of Mizoram (18,748 sq. km.) is under forest cover6 with Jhum cultivation practiced in 63% of the total crop area7. Malaria in Mizoram is predominantly (>â90%) transmitted by Plasmodium falciparum8. Anopheles minimus and An. dirus have been incriminated as malaria vectors in Mizoram9. While in the rest of India, Artemisinin combination therapy (ACT) is Artesunate plus Sulfadoxine-Pyrimethamine (AS-SP), in NE India (including Mizoram), SP resistance has led to the replacement of SP with the lumefantrine in the ACT (artemether-lumefantrine)10. Mizoram shares historical ties and porous international borders with the bordering districts of Myanmar and Bangladesh4,11. Widespread artemisinin resistance has been reported in Myanmar, and a recent study has indicated the emergence of artemisinin resistance independent of kelch13 mutations in Bangladesh12. Mizoram is considered to be one of the major routes for the entry of drug-resistant parasites from Southeast Asia (SEA) to NE India8,13. India, Bangladesh, and Myanmar have all committed to malaria eradication by 2030; the porous borders shared by Mizoram with these two countries could allow the movement of parasites and vectors in either-directions jeopardizing the malaria elimination efforts of all three countries.
In many parts of the world, malaria is seasonal and is influenced by short term changes in the meteorological variablesâtemperature, rainfall and humidity14. Temperature is a key factor that influences both the parasite and vector development; vectorâs life cycle from hatching of egg to progression of various developmental stages is influenced by temperature15,16. Malaria parasiteâs sporogonic cycle that determines the vectorial capacity and disease incidence is affected by temperature17 , and so is the vectorâs gonotrophic cycle that defines the frequency of mosquito-human interactions16. Rainfall (provides vector breeding sites) and suboptimal relative humidity (reduces mosquitoâs activity and life span) are other key environmental factors affecting malaria incidence16,18,19,20,21. The elevation is also associated with the intensity of disease transmission; low temperatures in elevated regions make it hard for the vectors to survive, resulting in lower vector densities and reduced/no disease transmission22,23,24,25,26. Given the hyperendemicity of malaria in Mizoram, it would be pertinent to study the role of environmental factors affecting malaria transmission in this region.
Geographic Information System (GIS) helps in identifying the geographical distribution of cases, spatio-temporal links of the disease, spatial clusters of the disease, and the population at risk27. Worldwide, GIS has been used as an effective tool to combat infectious diseases, including malaria28,29,30,31. As Mizoram shares porous international borders with Myanmar and Bangladesh, it is important to understand the spatio-temporal distribution of malaria cases for devising regional intervention strategies. Here, using GIS, we have mapped the clusters of malaria cases at the sub-health centre level across Mizoram from 2015 to 21, and have identified the hotspots. The environmental variables (temperature, humidity, rainfall, and land use/land cover) were also mapped to study their association with malaria. Furthermore, climate change over the last 4 decades and how it might influence malaria transmission in this region are discussed.
Methods
Data source and mapping
A malaria case is a person in whom the presence of malaria parasite in the blood has been confirmed by diagnostic tests. Malaria surveillance is carried out by the public health sub-center, which is the smallest unit, and the data will be reported to the district health headquarters on a daily basis. The data is collected mainly by passive surveillance, however, at the village level, active surveillance using Rapid Diagnostic Test (RDT) is also carried out in fever cases by Accredited Social Health Activist (ASHA) workers, who report to public health sub-centers. The public health sub-center level malaria datasets from 2015 to 2021 were collected from the Directorate of Health Services, Mizoram. The datasets represent 385 sub-centers (after data cleaning) across nine districts (Lawngtlai, Lunglei, Kolasib, Mamit, Aizawl East, Aizawl West, Siaha, Champhai, and Serchhip) of Mizoram. The datasets were compiled and prepared as a geodatabase comprising the key variables: total malaria cases, percentage of Plasmodium falciparum, Plasmodium vivax and Annual Parasite Index (API). The datasets used in the study is provided in Supplementary Table 1. API is calculated based on the total number of positive slides for the parasite in a yearâÃâ1000/Total population of a particular area. The population size under the surveillance of each health sub-center is surveyed and recorded by the respective centers. As the sub-center datasets are maintained by the district health headquarters, collection of the sub-center datasets for the entire state is tedious. Some of the districts provided only the district-level datasets to the Directorate of Health Services (Mizoram), and data were missing for Lawngtlai (2015â2018), Siaha (2018), and Aizawl East (2018) districts. Hence, in addition to the sub-center malaria datasets, district-wise datasets were also prepared for mapping the spatio-temporal distributions. The total malaria cases at district-level were normalized using the projected population of 2020, based on the population estimates of Wang et.al.32. All available API data of each health sub-center during 2015â2021 was used to estimate the average API.
To characterize the regional settings of Mizoram, climatic variables, digital terrain model, and land use/land cover (LU/LC) datasets of the Global Modelling and Assimilation Office (GMAO)33 and NASA satellite missions were retrieved from NASA Power34, LP DAAC35, and USGS Earth Explorer36 tools, respectively. The long-term gridded climatic datasets representing the annual average values covering the border districts of Lunglei and Lawngtlai (Mizoram), and Rangamati (Bangladesh) were extracted from 1981 to 2021 (41Â years) and used for long-term trend analysis of temperature (minimum, maximum and range), relative humidity, and rainfall. In addition, present and future climate classification maps were used to understand the shift in regional climatic pattern37. The LU/LC map of Mizoram and adjacent regions for the year 2022 was prepared using Landsat 8 satellite images. The details of satellite images used in the study are given in Supplementary Table S2. Supervised classification logic was adopted for preparing eight LU/LC classes by selecting training samples using ERDAS Imagine software. The parametric rule-based maximum likelihood classification method was used to extract the LU/LC information38. Accuracy assessment was carried out with 160 stratified random sample locations (a minimum of 10 per LU/LC class) by comparing high-resolution datasets (Google Earth images) and published LU/LC maps (International Geospatial-Biosphere Programmeâhttps://lpdaac.usgs.gov/products/mcd12q1v006/ and Bhuvan- https://bhuvan-app1.nrsc.gov.in/thematic/thematic/index.php). The overall classification accuracy is 85.63%, and the Kappa coefficient is 0.8450, which ensures an acceptable level of accuracy.
To extract environmental variables for statistical analysis, a one km buffer was used around the public health sub-centers. The values of elevation, climate, and LU/LC parameters within the buffer zone were extracted using the Zonal Statistics tool in ArcGIS 10.4 software (https://desktop.arcgis.com), and the mean value for each parameter was obtained. All the layers were prepared as spatial layers and mapped in the ArcGIS 10.4 software using dot, choropleth, and isopleth mapping techniques.
Statistical analysis
Out of the eight LU/LC classes identified in the study area, three classes (water, wetland, and cropland) were not considered for the analysis as they were not present in the buffer zone of any of the public health sub-center. The rest of the LU/LC classes are grouped into one land use class (built-up/Jhum), and two land cover classes (mixed forest/shrub, and dense forest). The KruskalâWallis H test was used to examine if the total number of malaria cases was different for (a) built-up/Jhum (nâ=â14), (b) mixed forest/shrub (nâ=â110), and (c) dense forest (nâ=â261) classes.
To estimate the association between the total malaria cases in all districts of Mizoram and the independent variables (elevation, minimum temperature, maximum temperature, temperature range, rainfall, and relative humidity), a generalized linear regression model was used. The dependent variable (total malaria cases) is highly skewed (Skewnessâ=â9.1 and Kurtosisâ=â116.7) (Supplementary Table S3) and multicollinearity issue was observed among the independent variables in the generalized linear regression analysis (Supplementary File S4). In addition, pair-wise correlation/association of all variables was determined using a correlation matrix (Supplementary Table S5), which showed a high correlation between the temperature minimum, temperature range, and temperature maximum. A high correlation was also observed between elevation, relative humidity, and rainfall. Therefore, temperature minimum and rainfall were not considered for the multiple regression analysis using a quasi-Poisson regression model. The statistical significance was determined using the adjusted incidence rate ratio (aIRR), which is obtained by exponentiating the coefficients of a quasi-Poisson regression model. The statistical analysis was performed using STATA 17 (StataCorp, USA) and ArcGIS Pro 3.0 (ESRI, USA) software.
Optimized hotspots analysis
To understand the sub-center level hotspots of malaria in Mizoram, the Optimized Hotspot Analysis tool in ArcGIS Pro 3.0 software was used. We used all available year-wise Annual Parasite Index (API) of the 385 public health sub-centers spread across the state (Supplementary File S6). The hotspots were identified based on the Getis-Ord Gi* statistic, using z-scores (standard deviations) and p-values (probability), of each sub-centers within the context of the neighboring sub-centers39. Extreme z-scores (very high or positive and very low or negative) with significant p-values, populating in the tails of the normal distribution curve are termed hot and cold spots, respectively. In our study, we have used multiple significant levels (90%, 95%, and 99% confidence intervals) for better understanding.
The Getis-Ord Gi* statistics is given as,
where xj is the API for the particular sub-center j in the selected year, wi,j is the spatial weight between neighboring sub-centers i and j, and n is the total number of sub-centers involved in our analysis40.
We could decide whether to accept the null hypothesis or not using the p-value from the significance test obtained for the Gi* statistic41.
Here, E(G) represents the expected G, and V(G) is the variance of G, which can be calculated as
When the absolute value of the Z-score is large and the p-value is less than 0.05, the null hypothesis can be rejected i.e., the sub-centers representing API exhibit statistically significant clustering.
Trend analysis
The MAKESENS template developed for time series analysis by the Finnish Meteorological Institute was used to estimate the trends in the climatic (temperature, rainfall, humidity) variables over 41Â years. The MAKESENS estimates the presence of a trend (either increasing or decreasing) using the nonparametric MannâKendall test, while its magnitude (true slope of the linear trend) is estimated by the Sensâs slope estimate42,43.
Ethical approval
This study is a secondary data analysis and does not involve experiments with humans/or the use of human tissue samples, therefore, human ethical approval is not applicable.
Results
Characteristics of study site
Mizoram, located in NE India, shares the majority of its western borders with the Chittagong Hill Tracts of Bangladesh, and its eastern and southern borders with the Chin Hills of Myanmar (Fig. 1). Among the eight major districts (Mamit, Kolasib, Aizawl, Champhai, Serchhip, Lunglei, Lawngtlai, and Siaha) in Mizoram, Lunglei (4536 sq. km) and Aizawl (3576 sq. km) are the largest. Aizawl, the capital, has the highest (78.6%) urban population, and Mamit has the highest (82.7%) rural population. Across Mizoram, the maximum and minimum mean temperature ranges from 33.9 to 35.7 °C and 11.3 to 14.0 °C respectively. The annual average rainfall is highest in Mamit (1807 mm), and lowest in Champhai (925 mm). Champhai and Lawngtlai have the highest and lowest mean elevation, respectively. The stateâs LU/LC is dominated by evergreen forests and shrublands. At least 90% of the land area is covered by different categories forests in each of the districts44. Detailed district-wise demographic and geographic characteristics of Mizoram are given in Supplementary Table S7, and the LU/LC distribution (in sq. km) is given in Supplementary Table S8.
Descriptive epidemiology
During the study period (2015â2021), 55,778 malaria cases were reported across the nine districts of Mizoram. Figure 1 shows the cumulative malaria cases from 2015 to 21, and the increase in caseloads from the Bangladesh border (lower elevation) to the Myanmar border (higher elevation). The highest number of cases were reported from Lawngtlai (20,074), followed by Mamit (15,631), and Lunglei (13,580); these districts share international borders with Bangladesh. A total of 47,287 (84.78%) Pf cases and 8,491 (15.2%) Pv cases were recorded. Both Pf and Pv cases were highest in Lawngtlai district; 37.78% (17,865 cases) and 26.02% (2209 cases), respectively. However, the proportions of Pv cases were highest in Siaha and Aizawl districts (>â50%). The average API throughout the study period across districts was 8.40. The highest average API was reported in 2015 (27.14), and the lowest in 2018 (4.33). From 23,849 cases in 2015, a steady decrease in malaria prevalence was observed till 2018 (4221). In the last three years (2019â21), cases have gradually increased; 4870 cases (2019), 6951 cases (2020), and 5166 cases (2021) (Supplementary Table S1).
The monthly distribution of total malaria cases in the Lunglei, Serchhip, and Mamit districts is shown in Fig. 2. The cases peak from June to September, which is the monsoon season in Mizoram45. Supplementary Fig. 9 shows the distribution of the API across the sub-centers, along with the total malaria cases reported in each districts, from 2015 to 2021. Throughout the study period, the sub-centers located in the western districts of Mizoram, especially Mamit, Lunglei, and Lawngtlai, which share international borders with Bangladesh, exhibited high API. As we move towards the eastern parts of Mizoram, a significant reduction in malaria cases was observed. Till 2019, API was reported to be below 1 in the sub-centers along the Myanmar border, but a surge in API can be seen in 2020 and 2021, where four sub-centers of Champhai reported an increase in API (1.1â10) and one center each in Lunglei and Serchhip reported API above 100 in 2021. Overall, the average API of Mizoram showed a decreasing trend, from 27.22 in 2015 to 5.55 in 2021. The number of sub-centers with APIâ>â10 reduced from 2015 to 2018; however, the numbers have started to rise from 2019 to 2021. Except for two, all the sub-centers reporting high API (>â10) are located in the districts that border Bangladesh. Even though sub-center-wise data is not available for Lawngtlai (2015â2018), Siaha (2018), and Aizawl East (2018), the district-level aggregated data indicate the cases were endemic along the Bangladesh border throughout the study period.
The distribution of Pf and Pv cases across the sub-centers of Mizoram from 2015 to 2021, with the terrain as a base map, is shown in Fig. 3. Throughout the study period, a high number of cases were reported in the sub-centers in the plain regions and the least number of cases in the high terrains. In 2015, 87.11% of the malaria cases were Pf, and it consistently remained the majority over the study period (77.50% in 2016, 85.93% in 2017, 86.93% in 2018, 92.36% in 2019, 80.85% in 2020, and 78.59% in 2021). The total Pv cases were 12.89% in 2015, and it increased to 22.50% in 2016; then, there was a gradual decline until 2019 (14.07% in 2017, 13.31% in 2018, and 7.64% in 2019). A steep increase was observed again in 2020 and 2021, with 19.15% and 21.41% of the total malaria cases being Pv, respectively. Most of the sub-centers with Pv cases are located in the urban districts (northern parts of Mizoram) of Aizawl East, Aizawl West, and Kolasib, whereas the sub-centers in the rural districts (mid and southern parts of Mizoram) of Mamit, Lunglei, and Lawngtlai, predominantly reported Pf infections. In recent years (2020â21), the proportion of Pv cases has increased considerably in the high-burden southern rural districts of Lawngtlai and Siaha.
Spatial clustering analysis
The average API reported in the sub-centers of Mizoram from 2015 to 2021, with a base map of elevation, is shown in Fig. 4a. The sub-centers located in the western districts of Mizoram, especially Mamit, Lawngtlai, and Lunglei, reported a high API during the study period, where the elevation was relatively low (plains). The API gradually decreased from the Bangladesh-Mizoram border to the Mizoram-Myanmar border, as the elevation increased. The majority of the sub-centers with an APIâ>â10.1 are located in the plain regions of Mizoram, while an API of 1.1â10.1 was observed in the transition zones, and towards the highly elevated areas, the majority of the sub-centers reported an API below 1.0. The malaria hotspots identified through the location and API of sub-centers in Mizoram are shown in Fig. 4b. Most of the sub-centers located along the Bangladesh border in the Lawngtlai and Lunglei districts are malaria hotspots. The hotspots also extended into the Mamit and Siaha districts, especially along the borders of Lawngtlai and Lunglei. Most sub-centers in Aizawl East, Aizawl West, Champhai, and Serchhip are malaria coldspots with an average API of less than one. The terrain differences determine the distribution of hot and cold spots; cold spots are in the higher elevated regions, while hotspots are at lower elevation.
The evolution of the malaria hotspots from 2015 to 2021 in Mizoram is shown in Supplementary Figure S10. The sub-centers in the Western border of Mizoram, primarily in Mamit, Lunglei, and Lawngtlai, are consistently the hotspots for malaria over the study period. In 2020 and 2021, the hotspots extended eastwards, with the emergence of malaria hotspots in the sub-centers of Siaha. A similar pattern was observed with the evolution of Pf hotspots in Mizoram (Supplementary Figure S11), as the majority of the total cases are Pf.
Pv hotspots were primarily concentrated in the sub-centers of Aizawl and Kolasib cities from 2015 to 2017. From 2018 to 2019, a shift in the Pv hotspots towards the south and western parts of Aizawl was observed. The majority of the sub-centers surrounding the city of Siaha emerged as hotspots for Pv in 2020 (Fig. 5).
Association with environmental variables
The spatial distribution of LU/LC of Mizoram and the bordering areas is shown in Fig. 6. The shrubland is the dominant LU/LC type along the Western border of Mizoram, especially in Mamit, Lunglei, and Lawngtlai. The LU/LC was dominated by dense forests towards the Eastern border. A Jhum type of cultivation is dominant in the Eastern hilly terrains, especially in the Champhai district.
The spatial distribution of maximum temperature, minimum temperature, temperature range, rainfall, and relative humidity from 1981 to 2021 are shown in Fig. 7. Analysis of the annual meteorological data from 1981 to 2021 shows the western districts of Mizoram experience similar climatic conditions as the bordering districts located in the Chittagong division of Bangladesh. The terrain settings of Mizoram, along with Eastern and Western border regions, are also presented in Fig. 7.
The spatial clusters of malaria are overlaid in the maps of environmental variables for a better understanding of malaria and its environmental settings (Figs. 6 and 7). The LU/LC distribution shows the malaria hotspots of Mizoram are primarily located in the shrublands that border the mixed or dense forests. Such a pattern of LU/LC is also observed in the Chittagong Hill Tract of Bangladesh, with a transition of deciduous forests to shrublands to croplands. Analysis showed a statistically significant difference in the number of malaria cases between the mean ranks of at least one pair of groups, Ï2 (2)â=â19.652, pâ<â0.001. The pairwise comparison matrix estimated by the KruskalâWallis H test also shows a significant positive association (pâ<â0.0001) of forest categories (dense, mixed, and shrub forests) and malaria cases (Table 1 and Supplementary Table 12).
The spatial overlay of malaria clusters and climatic parameters shows the maximum temperature experienced in the malaria hotspot regions of Mizoram is 35â36 °C (Figure 7a), the minimum temperature is 13â15 °C (Fig. 7b), and the average diurnal range is 15â17 °C (Fig. 7c). Malaria hotspots were primarily clustered in Lawngtlai, Lunglei, Mamit, and Siaha districts which are with relatively higher temperatures (maximum is above 35 °C and the minimum is 13â14 °C), whereas coldspots were observed in the elevated districts (Champhai, Aizawl, and Serchhip) with relatively lower temperatures (maximum is less than 34 °C and the minimum is less than 13 °C). A similar pattern is observed with the average rainfall and relative humidity (Fig. 7d and e). The average annual rainfall and relative humidity in this hotspot region are 1500â2100 mm and 70â75%, respectively. All these climatic parameters are identical in the bordering regions of Bangladesh and the neighboring Indian state of Tripura. The annual average temperature (maximum and minimum), rainfall, and relative humidity displayed a decreasing pattern from the Chittagong Hill Tract of Bangladesh towards the Myanmar border of Mizoram (Fig. 7a,b,d, and e). All these LU/LC and climatic parameters are highly interrelated with the regional elevation (Fig. 7f). All the hotspots fall within 500 m (above MSL) elevation, and the coldspots are located predominantly above 1000 m (above MSL).
The multivariate association between the total malaria cases and environmental variablesâtemperature (min, max, and range), elevation, relative humidity, and LU/LC estimated by the quasi-Poisson regression model is presented in Table 2. A statistically significant association was observed between the total malaria cases and the environmental variables (temperature range, elevation, dense forest, and mixed forest/scrub). Total malaria cases significantly increased with a decrease in the temperature range (aIRR 0.220; 95%CI 0.099â0.490) and elevation (aIRR 0.997; 95%CI 0.995â0.999). The dense forest (aIRR 3.843; 95%CI 1.881â7.852) and mixed forest/shrub (aIRR 5.395; 95%CI 2.652â10.975) have a higher rate of malaria cases.
Climatic trend analysis
The climatic trend in the malaria hotspot (Western border) and cold-spot regions (Eastern border) of Mizoram from 1981 to 2021 are shown in Fig. 8, Supplementary Figure S13, and Supplementary Table S14. A decreasing trend in the maximum temperature and diurnal temperature range (Q: 0.05 °C/year and ââ0.07 °C/year, at Pâ<â0.001, respectively) is observed in the Western region. The minimum temperature increased at a rate of 0.015 °C/year (Pâ<â0.01). Also, the relative humidity (Qâ=â0.3%/year; Pâ<â0.001) and rainfall (Qâ=â34.28 mm; Pâ<â0.001) of the region showed an increasing trend over the past four decades. Till 2006, the annual average rainfall recorded was below 2000 mm per year, and a significant increasing trend was observed in the past decade, with the maximum rainfall recorded in 2017 (5204 mm).
In the Eastern Border of Mizoram, the climatic trend (1981â2021) follows a similar pattern as the western part. The maximum (Q: 0.03 °C/year; Pâ<â0.05) and minimum temperature, and the diurnal temperature range show a decreasing trend, while the average rainfall (Q: 19.8 °C/year; Pâ<â0.001) and relative humidity (Q: 0.38 °C/year; Pâ<â0.05) follow an increasing trend.
Using the Köppen-Geiger climate classification maps37, the malaria hotspots are compared with the present (1980â2016) climatic types (Fig. 9a). The overlay analysis shows the malaria hotspots to fall predominately under the Tropical Monsoon (Am) type climate, which is characterized by seasonal rainfall and driest months. The future climatic map (Fig. 9b) was derived based on an ensemble of climate projections from the 32 CMIP5 models (2071â2100) corresponding to the Representative Concentration Pathway (RCP) 8.5 climate scenario. It shows the climatic types of eastern parts of Mizoram will shift from the present temperate (hot summer) type to the tropical monsoon type, which is also evidenced by the climatic trend analysis (Supplementary Figure S13).
Discussion
In the last decade, India has seen a steep decline in malaria morbidity and mortality46. In most parts of the country, malaria prevalence and incidence have been consistently declining46. However, Mizoram is one of the few malaria-endemic pockets where the decrease in malaria prevalence has been relatively low and has even increased in the last two years, despite the COVID-19 pandemic. As per National Framework for Malaria Elimination (NFME), Mizoram falls under category 3 (intensified control phase); these are States/Union Territories having APIâ>â147. Malaria control in Mizoram is critical for Indiaâs malaria elimination efforts, as its western and eastern districts share porous international borders with Bangladesh and Myanmar, countries that are endemic to malaria.
Mamit, Lunglei, and Lawngtlai, the western districts of Mizoram, share international borders with Chittangong Hill Tract (CHT) of Bangladesh; the CHT comprises Rangamati, Bandarban and Khagrachhari districts. Of the 64 districts in Bangladesh, malaria is endemic in 13, and 19 million people are at risk48. In Bangladesh, from 2008 (84,690 cases) to 2020 (6,130 cases), malaria incidence has steeply declined byâ~â93%49. The remaining malaria cases are from 13 of the 64 districts, and of these 13 districts, 3 (Rangamati, Bandarban and Khagrachhari) in the Chittagong Hill Tracts (CHT) contribute toâ~â90% of the malaria cases50. GIS analysis of public health sub-centre-wise malaria data of Mizoram from 2015 to 21 shows malaria is highly endemic in the western districts, and as we move towards the east, the malaria cases decreases steeply. It is best illustrated by the steep decrease in malaria cases as we move from Mamit district (>â10,000 cases) in the west, traversing the central districts of Kolasib/Aizawl and/or Serchhip to Champhai district (<â200) in the east. The high endemicity of malaria in the CHT and western districts of Mizoram could be attributed to the favorable meteorological parameters-temperature, humidity, and rainfall. Studies predict malaria transmission is bounded by the thermal optimum of 17 °C and 34 °C51,52,53. In addition, the relatively plain terrain of this region would also be conducive for the breeding of vectors54. As we move towards the East, the temperature and RH drops, and the altitude from the mean sea level increases (Fig. 7), resulting in less favorable breeding conditions for mosquitoes and malaria transmission. Similar to Aizawl (Mizoram)21, in CHT, the minimum temperature, RH and rainfall have shown a consistent upward trend in the last 4 decades, thus favoring mosquito breeding and malaria transmission. Multiple global climate models show an overall net increase in populations at risk and climate suitability for malaria55. The trend analysis shows there is a gradual shift in average maximum temperature towards less than 34 °C, minimum temperature towards 18 °C and reduction in diurnal temperature, thus favoring malaria transmission in the CHT and western districts of Mizoram. In Western Mizoram, over the last 3 decades, an estimatedâ~â34 mm annual increase in rainfall is observed, and the increase is marked from 2015 onwards. If the increase in rainfall is also characterized by decreased spacing, it could lead to rapid expansion of mosquitoes56. Analysis of Eastern (current study) and Central Mizoram21 shows similar overall climatic trends as observed in the western region. Elevation in the eastern region of Mizoram appears to be a major factor limiting the spread of malaria. However, the predicted future climate change scenario indicates a tropical monsoon climate throughout Mizoram, paving the way for vector breeding and sustained malaria transmission. In the highland regions of Ethiopia and Colombia, 0.2 °C increase in temperature per decade has led to malaria transmission in these countriesâ high-altitude regions57,58. Therefore, unless malaria transmission is stringently controlled in the CHT and western districts of Mizoram, the effect of climate change could result in a more conducive ecological niche favoring malaria transmission in the central and eastern parts of Mizoram.
After a gradual decline in malaria prevalence from 2015 to 18, Mizoram has seen an uptick in cases from 2019 onwards, especially in its western districts. This trend is worrisome, especially with escalating malaria and the recent emergence of K13 independent ART resistance in the CHT of Bangladesh12. Molecular screening of Pfkelch13 in NE carried out in 2014â15 has identified non-synonymous mutations (K189T and A578S) in the Lunglei district, Mizoram59. ACT efficacy studies carried out in the Lunglei district from 2011 to 2013 showed no correlation between K13 propeller mutation (A578S) and treatment failures60. In the last decade, studies have reported delayed clearance phenotypes from parasites in Central (Chhattisgarh), South-West (Goa), and NE (Assam, Arunachal Pradesh, and Tripura) India, independent of K13 gene polymorphisms61,62. Decreased ART sensitivity across India and CHT and the recent increase in malaria cases in Mizoramâs international borders necessitate the need for continuous surveillance of the parasite phenotypes, especially in Mizoramâs western districts. Unlike molecular surveillance, which can be carried out in filter paper spots, assessing ART resistance through ring-stage survival assays (RSA) will be challenging as it will involve the collection and quick transport of blood samples from these arduous terrains to a malaria culture facility. Effective malaria control in the Indo-Bangladesh border districts would require coordinated efforts by both countries. Unless malaria is effectively controlled in CHT, the porous Indo-Bangladesh borders will allow the movement of parasites and vectors to Mizoramâs western districts and vice versa.
The eastern districts of Mizoram share international borders with the Chin State of Myanmar. In addition to the geographical proximity, Mizoram shares deep cultural and religious ties with Myanmar, and the people are bonded by similar customs, traditions, and beliefs11. In these porous borders, people from Mizoram and Myanmar continue to move in both directions, leading to the free movement of parasites in either direction, complicating malaria control. Historically, SEA has been a hotspot for the generation of acquired resistance to pyrimethamine, chloroquine, sulphadoxine, quinine, mefloquine, and most worryingly, artemisinins8. The first instance of drug resistance to chloroquine and SP in India has been from the NE region; these drug-resistant parasites are suspected to have entered the NE from SEA through the porous international borders63. In a state-wise study carried out across Mizoram from 2015 to 17, all the mutations identified in the pfcrt gene conferring chloroquine resistance belonged to the South-East Asian CVIET haplotype, suggesting an evolutionary link between the parasites in Mizoram and SEA8. In the Chin state of Myanmar that borders Mizoram, the K13 propeller mutations have been very low, and previous studies have not reported any decrease in ART sensitivity in this region64,65. However, in Homalin from the Sagaing region of Myanmar, 25Â km from Manipur (India),â~â47% of the parasites carried propeller mutations64. Mizoram shares domestic borders with Manipur State in the north64. Unlike the western districts bordering Bangladesh, the elevated altitude and climatic variables are less conducive to mosquito breeding and malaria transmission in the Mizoram-Myanmar border. The current absence of ART-resistant parasites does not preclude the future entry of drug-resistant parasites through Mizoram-Myanmarâs porous international borders. GIS-based analysis shows the emergence of malaria (>â10 API) in 47 sub-centers in 2021, and out of these, the highest API (>â100 API) is recorded in two sub-centers of Serchhip and Lunglei districts, bordering Myanmar. These clusters have to be monitored to assess the emergence of drug-resistant parasites.
In these high transmission malaria settings, it is also pertinent to track the status of Pv drug resistance. Pv parasites resistant to chloroquine have emerged in most Pv endemic countries66. Literature on Pv drug resistance is sparse in NE, and there are no reports on the status of Pv drug resistance in Mizoram. In Indiaâs urban settings, Pv predominates67, and GIS analysis shows Pv cases to be concentrated in the urban centers of Aizawl East, Aizawl West, Kolasib and Siaha. Therapeutic efficacy studies, as assessed by response to treatment on day 2868, have to be carried out in Lungtian, Lungpher, Sangau, and Vawmbuk sub-centers in Lawngtlai, where Pv predominates.
In addition to parasite surveillance within these international borders, it is imperative to continuously assess the vector population. A district-wise mosquito survey carried out by our team at RMRC (ICMR) Dibrugarh, from late 2000 to early 2010 shows a high density of Anopheles minimus in Western Mizoram and South Tripura, bordering Bangladesh. As observed here, urban malaria in India is mostly transmitted by Pv69, and the reasons are not clear. It is possible specific vectors well suited for survival in urban settings might transmit Pv more efficiently. Also, Pv appears to dominate at the border between plain and elevated areas. It is possible, the ecological niche and the vectors determine the distribution of Pv and Pf. Vector epidemiology and vector-parasite relationship can change spatio-temporally70,71, and is critical to devise vector-control strategies based on the existing malaria vectors.
In Mizoram,â~â54% of the rural households are engaged in shifting cultivation, also known as Jhum cultivation, which involves slashing and burning of selected forest area, followed by a short period of growing mixed crops, harvesting, following, and subsequently, the site will be abandoned72,73. In Mizoram, for Jhum cultivation, an estimated 0.2 million acre of forest land is burnt annually, resulting in significant greenhouse emissions7, which could be a contributor to the climate change trends. Furthermore, Jhum cultivation continuously exposes the farmers to malaria vectors as they stay and sleep in the forests. By opening up the forests, Jhum cultivation leads to associated streams, creating favorable breeding grounds for An. minimus74. In recent years, the Jhum cycle has reduced from 10â20 years to 2â3 years, and encroachment of virgin forested areas has led to new perennial malaria transmission zones75. As Jhum cultivation is widely practiced in Champhai and other districts bordering Myanmar (Fig. 7), there are possibilities for surge in malaria cases. In recent years, the sub-centers withâ>â1.0 API are increasing in the districts bordering Myanmar. In 2021, two sub-centers along the Myanmar border, where Jhum cultivation is widely practiced, recorded an API ofâ>â100. Thus, malaria control strategies should continuously be implemented in this region and should track, test, and treat farmers involved in Jhum cultivation.
As our analysis shows, 5 of the 8 districts in Mizoram are high malaria transmission settings, and continuous exposure to Plasmodium in these endemic areas may result in low-density infections and/or asymptomatic carriers. The low-density infections from asymptomatic carriers can lead to malaria transmission76,77. For effective malaria control, strategies should include mass screening of the community with molecular tools. Continuous molecular surveillance and prompt treatment will greatly help lessen the malaria burden and break the transmission cycle in these high-transmission settings.
One of the important limitations is the non-availability of sub-center level data month-wise and year-wise for some of the districts. Also, absence of socio-economic status datasets at the sub-center level meant appropriate confounding factors were not included in the quasi-Poisson analysis. Even though we detail the role of conducive environmental factors, it will be hard to conclude malaria endemicity in the western districts is solely due to favorable climate unless confounding factors such as socioeconomic determinants are considered and controlled. It is also possible COVID-19 pandemic might have restricted the resources, which in turn could have displaced Pv hotspots in 2020 and 2021.
Overall, an integrated management approach involving regional cooperation, molecular surveillance in sentinel sites, laboratory investigation of ART resistance using RSA, community-level screening for detection of asymptomatic carriers, appropriate vector control strategies tailor-made for existing vectors, developing skilled manpower and infrastructure for molecular studies of malaria are required to effectively control malaria in this region.
Data availability
All data analysed during this study are included in the manuscript (and the Supplementary Information files).
References
WHO. World malaria report 2021, https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2021 (2021).
WHO. More malaria cases and deaths in 2020 linked to COVID-19 disruptions, https://www.who.int/news/item/06-12-2021-more-malaria-cases-and-deaths-in-2020-linked-to-covid-19-disruptions (2022).
Barman, A. S. et al. DNA Barcoding of freshwater fishes of Indo-Myanmar biodiversity hotspot. Sci. Rep. 8, 8579. https://doi.org/10.1038/s41598-018-26976-3 (2018).
The Land Systems of Mizoram, https://landrevenue.mizoram.gov.in/uploads/files/historical-evolution-of-mizoram.pdf (1967).
Roluahpuia,. Unsettled autonomy: Ethnicity, tribes and subnational politics in Mizoram, North-east India. Nations Natl. 27, 412â426. https://doi.org/10.1111/nana.12681 (2021).
Environment, Forests & Climate Change Department, https://forest.mizoram.gov.in/page/mizoram-forest-at-a-glance (2022).
Zomuanpuii, R. et al. Epidemiology of malaria and chloroquine resistance in Mizoram, Northeastern India, a malaria-endemic region bordering Myanmar. Malar. J. 19, 95. https://doi.org/10.1186/s12936-020-03170-3 (2020).
Das, S. C. & Baruah, I. Incrimination of Anopheles minimus Theobald and Anopheles balabacensis balabacensis Baisas (A. dirus) as malaria vectors in Mizoram. Indian J Malariol. 22, 53â55 (1985).
Anvikar, A. R. et al. Antimalarial drug policy in India: Past, present & future. Indian J. Med. Res. 139, 205â215 (2014).
Vanlalruata, H. C. The Mizos and Myanmarese share deep cultural bonds that go beyond the citizen-refugee divide, https://www.thehindu.com/society/the-mizos-and-myanmarese-share-deep-cultural-bonds-that-go-beyond-the-citizen-refugee-divide/article34391179.ece (2021).
Nima, M. K. et al. Assessment of Plasmodium falciparum artemisinin resistance independent of kelch13 polymorphisms and with escalating malaria in Bangladesh. MBio https://doi.org/10.1128/mbio.03444-21 (2022).
Sarma, D. K. et al. Malaria in North-East India: Importance and implications in the era of elimination. Microorganisms https://doi.org/10.3390/microorganisms7120673 (2019).
Thomson, M. C. & Stanberry, L. R. Climate change and vectorborne diseases. N. Engl. J. Med. 387, 1969â1978. https://doi.org/10.1056/NEJMra2200092 (2022).
Depinay, J. M. et al. A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar. J. 3, 29. https://doi.org/10.1186/1475-2875-3-29 (2004).
Ye, Y., Louis, V. R., Simboro, S. & Sauerborn, R. Effect of meteorological factors on clinical malaria risk among children: An assessment using village-based meteorological stations and community-based parasitological survey. BMC Public Health 7, 101. https://doi.org/10.1186/1471-2458-7-101 (2007).
Okuneye, K., Eikenberry, S. E. & Gumel, A. B. Weather-driven malaria transmission model with gonotrophic and sporogonic cycles. J. Biol. Dyn. 13, 288â324. https://doi.org/10.1080/17513758.2019.1570363 (2019).
Yamana, T. K. & Eltahir, E. A. B. Incorporating the effects of humidity in a mechanistic model of Anopheles gambiae mosquito population dynamics in the Sahel region of Africa. Parasit. Vectors 6, 235. https://doi.org/10.1186/1756-3305-6-235 (2013).
Lindsay, S. W. & Martens, W. J. Malaria in the African highlands: Past, present and future. Bull. World Health Organ. 76, 33â45 (1998).
Lindsay, S. W. & Birley, M. H. Climate change and malaria transmission. Ann. Trop. Med. Parasitol. 90, 573â588. https://doi.org/10.1080/00034983.1996.11813087 (1996).
Karuppusamy, B. et al. Effect of climate change and deforestation on vector borne diseases in the North-Eastern Indian State of Mizoram bordering Myanmar. J. Clim. Change Health 2, 100015. https://doi.org/10.1016/j.joclim.2021.100015 (2021).
Global Health. Where Malaria Occurs, Centers for Disease Control and Prevention, https://www.cdc.gov/malaria/about/distribution.html (2020).
Omumbo, J. A., Hay, S. I., Goetz, S. J., Snow, R. W. & Rogers, D. J. Updating historical maps of malaria transmission intensity in East Africa using remote sensing. Photogramm. Eng. Remote Sens. 68, 161â166 (2002).
Drakeley, C. J. et al. Altitude-dependent and -independent variations in Plasmodium falciparum prevalence in northeastern Tanzania. J. Infect. Dis. 191, 1589â1598. https://doi.org/10.1086/429669 (2005).
Dhiman, R. C., Yadav, Y. K., Saraswat, S. & Singh, P. Altitude, temperature, and malaria vectors in Nainital and Udham Singh Nagar districts of Uttarakhand, India: An evidence-based study. J. Vector Borne Dis. 50, 220â224 (2013).
Kassam, N. A. et al. Ten years of monitoring malaria trend and factors associated with malaria test positivity rates in Lower Moshi. Malar. J. 20, 193. https://doi.org/10.1186/s12936-021-03730-1 (2021).
Saran, S., Singh, P., Kumar, V. & Chauhan, P. Review of geospatial technology for infectious disease surveillance: Use case on COVID-19. J. Indian Soc. Remote Sens. 48, 1121â1138. https://doi.org/10.1007/s12524-020-01140-5 (2020).
Eisen, R. J. & Eisen, L. Concepts and Methods in Infectious Disease Surveillance 219â229 (Wiley, 2014).
Hay, S. I. et al. A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 6, e1000048. https://doi.org/10.1371/journal.pmed.1000048 (2009).
Hay, S. I., Guerra, C. A., Tatem, A. J., Noor, A. M. & Snow, R. W. The global distribution and population at risk of malaria: Past, present, and future. Lancet. Infect. Dis 4, 327â336. https://doi.org/10.1016/S1473-3099(04)01043-6 (2004).
Snow, R. W., Guerra, C. A., Noor, A. M., Myint, H. Y. & Hay, S. I. The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature 434, 214â217. https://doi.org/10.1038/nature03342 (2005).
Wang, W., Kim, R. & Subramanian, S. V. Population estimates for districts and parliamentary constituencies in india, 2020. Harvard Dataverse V1 https://doi.org/10.7910/DVN/RXYJR6 (2021).
NASA. Global Modelling and Assimilation Office, https://gmao.gsfc.nasa.gov/reanalysis/MERRA-2/ (2022).
NASA. The Power Data Access Viewer, https://power.larc.nasa.gov/data-access-viewer/ (2023).
NASA. NASA Shuttle Radar Topography Mission Global 30 arc second, https://lpdaac.usgs.gov/products/srtmgl30v021/ (2022).
USGS. EarthExplorer, https://earthexplorer.usgs.gov.
Beck, H. E. et al. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 5, 180214. https://doi.org/10.1038/sdata.2018.214 (2018).
Maximum Likelihood, L3Harris Geospatial documentation center., <https://www.l3harrisgeospatial.com/docs/maximumlikelihood.html> (2022).
How Hot Spot Analysis (Getis-Ord Gi*) works, https://pro.arcgis.com/en/pro-app/latest/tool-reference/spatial-statistics/h-how-hot-spot-analysis-getis-ord-gi-spatial-stati.htm.
Ord, J. K. & Getis, A. Local spatial autocorrelation statistics: Distributional issues and an application. Geogr. Anal. 27, 286â306. https://doi.org/10.1111/j.1538-4632.1995.tb00912.x (1995).
Zhang, K. & Zhang, S. in 2015 23rd International Conference on Geoinformatics. 1â4 (IEEE).
Gilbert, R. O. Statistical Methods for Environmental Pollution Monitoring (Pacific Northwest National Lab. (PNNL), 1987).
Finnish Meteorological Institute, MAKESENS-application for trend calculation, https://en-beta.ilmatieteenlaitos.fi/makesens.
Environment, Forests & Climate Change Department, Government of Mizoram, https://forest.mizoram.gov.in/page/mizoram-forest-at-a-glance (2023).
CENTRAL GROUND WATER BOARD, Government of India, http://cgwb.gov.in/GW-Assessment/GWR-2020-Reports%20State/Mizoram_State_Report_Resource_2020.pdf (2022).
WHO World Malaria Report: India continues to make Impressive Gains in reduction of Malaria Burden, https://www.pib.gov.in/PressReleasePage.aspx?PRID=1677601 (2020).
NVBDCP. National framework for malaria elimination in india, https://nvbdcp.gov.in/WriteReadData/l892s/National-framework-for-malaria-elimination-in-India-2016%E2%80%932030.pdf (2016).
Haque, U. et al. Malaria burden and control in Bangladesh and prospects for elimination: An epidemiological and economic assessment. Lancet Glob. Health 2, e98-105. https://doi.org/10.1016/S2214-109X(13)70176-1 (2014).
Naher, K. M. Liberating Bangladesh from malaria: How far have we come in the last 50 years?, https://blog.brac.net/liberating-bangladesh-from-malaria-how-far-have-we-come-in-the-last-50-years/ (2021).
Clancy, E. Is Bangladesh entering the Era of Malaria Elimination?, https://speakingofmedicine.plos.org/2022/04/25/is-bangladesh-entering-the-era-of-malaria-elimination/ (2022).
Ryan, S. J. et al. Mapping physiological suitability limits for malaria in Africa under climate change. Vector Borne Zoonotic Dis. 15, 718â725. https://doi.org/10.1089/vbz.2015.1822 (2015).
Shah, M. M. et al. Malaria smear positivity among Kenyan children peaks at intermediate temperatures as predicted by ecological models. Parasit. Vectors 12, 288. https://doi.org/10.1186/s13071-019-3547-z (2019).
Mordecai, E. A. et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol. Lett. 16, 22â30. https://doi.org/10.1111/ele.12015 (2013).
Asigau, S., Hartman, D. A., Higashiguchi, J. M. & Parker, P. G. The distribution of mosquitoes across an altitudinal gradient in the Galapagos Islands. J. Vector Ecol. 42, 243â253. https://doi.org/10.1111/jvec.12264 (2017).
Caminade, C. et al. Impact of climate change on global malaria distribution. PNAS 111, 3286â3291. https://doi.org/10.1073/pnas.1302089111 (2014).
Shaman, J. & Day, J. F. Reproductive phase locking of mosquito populations in response to rainfall frequency. PLoS ONE 2, e331. https://doi.org/10.1371/journal.pone.0000331 (2007).
Lyon, B., Dinku, T., Raman, A. & Thomson, M. C. Temperature suitability for malaria climbing the Ethiopian Highlands. Environ. Res. Lett. 12, 064015. https://doi.org/10.1088/1748-9326/aa64e6 (2017).
Siraj, A. S. et al. Altitudinal changes in malaria incidence in highlands of Ethiopia and Colombia. Science 343, 1154â1158. https://doi.org/10.1126/science.1244325 (2014).
Mishra, N. et al. Emerging polymorphisms in falciparum Kelch 13 gene in Northeastern region of India. Malar. J. 15, 583. https://doi.org/10.1186/s12936-016-1636-4 (2016).
Mishra, N. et al. Surveillance of artemisinin resistance in Plasmodium falciparum in India using the kelch13 molecular marker. Antimicrob. Agents Chemother. 59, 2548â2553. https://doi.org/10.1128/AAC.04632-14 (2015).
Das, S. et al. Artemisinin combination therapy fails even in the absence of Plasmodium falciparum kelch13 gene polymorphism in Central India. Sci. Rep. 11, 9946. https://doi.org/10.1038/s41598-021-89295-0 (2021).
Chakrabarti, R. et al. Decreased in vitro artemisinin sensitivity of Plasmodium falciparum across India. Antimicrob. Agents Chemother. 63, e00101-00119. https://doi.org/10.1128/AAC.00101-19 (2019).
Mohapatra, P. K. et al. Molecular evidence of increased resistance to anti-folate drugs in Plasmodium falciparum in North-East India: A signal for potential failure of artemisinin plus Sulphadoxine-Pyrimethamine combination therapy. PLoS ONE 9, e105562. https://doi.org/10.1371/journal.pone.0105562 (2014).
Tun, K. M. et al. Spread of artemisinin-resistant Plasmodium falciparum in Myanmar: A cross-sectional survey of the K13 molecular marker. Lancet. Infect. Dis 15, 415â421. https://doi.org/10.1016/S1473-3099(15)70032-0 (2015).
Nyunt, M. H. et al. Molecular assessment of artemisinin resistance markers, polymorphisms in the K13 propeller, and a multidrug-resistance gene in the Eastern and Western border areas of Myanmar. Clin. Infect. Dis. 60, 1208â1215. https://doi.org/10.1093/cid/ciu1160 (2015).
Price, R. N. et al. Global extent of chloroquine-resistant Plasmodium vivax: A systematic review and meta-analysis. Lancet. Infect. Dis 14, 982â991. https://doi.org/10.1016/S1473-3099(14)70855-2 (2014).
Wilson, M. L. et al. Urban malaria: Understanding its epidemiology, ecology, and transmission across seven diverse ICEMR network sites. Am. J. Trop. Med. Hyg. 93, 110â123. https://doi.org/10.4269/ajtmh.14-0834 (2015).
WHO. Methods for surveillance of antimalarial drug efficacy, https://apps.who.int/iris/bitstream/handle/10665/44048/9789241597531_eng.pdf (2009).
Anvikar, A. R. et al. Epidemiology of plasmodium vivax malaria in India. Am. J. Trop. Med. Hyg. 95, 108â120. https://doi.org/10.4269/ajtmh.16-0163 (2016).
Alam, M. S. et al. Prevalence of anopheline species and their Plasmodium infection status in epidemic-prone border areas of Bangladesh. Malar. J. 9, 15. https://doi.org/10.1186/1475-2875-9-15 (2010).
Kumar, A. et al. Anopheles subpictus carry human malaria parasites in an urban area of Western India and may facilitate perennial malaria transmission. Malar. J. 15, 124. https://doi.org/10.1186/s12936-016-1177-x (2016).
Sati, V. P. Shifting cultivation in Mizoram, India: An empirical study of its economic implications. J. Mt. Sci. 16, 2136â2149. https://doi.org/10.1007/s11629-019-5416-9 (2019).
Soni, J. K., Hota, S. & Shakuntala, I. Jhum Cultivation: The Way of Livelihood Security for Mizoram (Indian Council of Agricutural Research, 2020).
Rao, T. R. The Anophelines of India (Malaria Research Centre, 1984).
Sarmah, N. P. et al. Role of Anopheles baimaii: Potential vector of epidemic outbreak in Tripura, North-east India. J. Glob. Health Rep. 3, e2019036. https://doi.org/10.29392/joghr.3.e2019036 (2019).
Lindblade, K. A., Steinhardt, L., Samuels, A., Kachur, S. P. & Slutsker, L. The silent threat: Asymptomatic parasitemia and malaria transmission. Expert Rev. Anti Infect. Ther. 11, 623â639. https://doi.org/10.1586/eri.13.45 (2013).
Sturrock, H. J. W. et al. Targeting asymptomatic malaria infections: Active surveillance in control and elimination. PLoS Med. 10, e1001467. https://doi.org/10.1371/journal.pmed.1001467 (2013).
Acknowledgements
The team would like to express their sincere gratitude to the following officers and staff of State Vector Borne Disease Control Programme (SVBDCP), Mizoram for their untiring supportâDr. Hmingthanmawii, Joint Director (Malaria), Dr. Lalthlamuana , Deputy Director (Malaria) and State Consultants of SVBDCâDr. Lalthlannunga, Zonundangi, Dr. K. Vanlalhruaia, Deborah Lalpuii, Malsawmthangi Pachuau, Thangniang Lunzosang Vaiphei, B. Lalkrawsengmawia, and District Consultants of SVBDCPâJosephine Lalramthari, K. Lalmalsawmi, Lalruatkimi Chenkual, Lalhmingmawii, Lalchhuanawmi Khiangte, Bibian Lalchhuanliana, Vanchamliana Chhangte, C. Kaipotha and Vabeimosa Syhlo.
Funding
This study was supported by the Department of Health Research grant (R.11013/11/2021-GIA/HR), funded to PBN, KB and DKS.
Author information
Authors and Affiliations
Contributions
P.B.N. conceptualized and drafted the original manuscript; P.L., K.B., M.M.J. and K.A.P. performed data curation, and data analysis; P.L., K.B., D.K.S., M.M.J., L.P. and K.A.P., reviewed and edited the manuscript; All authors have read and agreed to the final version of the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lalmalsawma, P., Balasubramani, K., James, M.M. et al. Malaria hotspots and climate change trends in the hyper-endemic malaria settings of Mizoram along the IndiaâBangladesh borders. Sci Rep 13, 4538 (2023). https://doi.org/10.1038/s41598-023-31632-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-31632-6
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.