- Open Access
The role of age, ethnicity and environmental factors in modulating malaria risk in Rajasthali, Bangladesh
Malaria Journal volume 10, Article number: 367 (2011)
Malaria is endemic in the Rajasthali region of the Chittagong Hill Tracts in Bangladesh and the Rajasthali region is the most endemic area of Bangladesh. Quantifying the role of environmental and socio-economic factors in the local spatial patterns of malaria endemicity can contribute to successful malaria control and elimination. This study aimed to investigate the role of environmental factors on malaria risk in Rajasthali and to quantify the geographical clustering in malaria risk unaccounted by these factors.
A total of 4,200 (78.9%; N = 5,322) households were targeted in Rajasthali in July, 2009, and 1,400 individuals were screened using a rapid diagnostic test (Falci-vax). These data were linked to environmental and socio-economic data in a geographical information system. To describe the association between environmental factors and malaria risk, a generalized linear mixed model approach was utilized. The study investigated the role of environmental factors on malaria risk by calculating their population-attributable fractions (PAF), and used residual semivariograms to quantify the geographical clustering in malaria risk unaccounted by these factors.
Overall malaria prevalence was 11.7%. Out of 5,322 households, 44.12% households were living in areas with malaria prevalence of ≥ 10%. The results from statistical analysis showed that age, ethnicity, proximity to forest, household density, and elevation were significantly and positively correlated with the malaria risk and PAF estimation. The highest PAF of malaria prevalence was 47.7% for third tertile (n = 467) of forest cover, 17.6% for second tertile (n = 467) of forest cover and 19.9% for household density >1,000.
Targeting of malaria health interventions at small spatial scales in Bangladesh should consider the social and socio-economic risk factors identified as well as alternative methods for improving equity of access to interventions across whole communities.
Malaria eradication is the ultimate goal of the World Health Organization (WHO) . Most malaria endemic countries are presently shifting their efforts from malaria control to eradication . The WHO has targeted eight to ten countries for elimination of malaria by 2015 and afterwards in all other endemic countries . Although the long-term goal is world-wide eradication, there remains some debate over the feasibility of elimination in Africa and eleven Asian malaria endemic countries [2–5]. Accurate maps of malaria incidence are important tools in malaria control as they can guide interventions and assess their effectiveness. Risk maps and geo-spatial data are increasingly used to support malaria elimination .
The Bangladesh national malaria survey in 2007 indicated that malaria was endemic in 13 of 64 administrative districts and that the crude prevalence is 4.0% [7, 8]. In Bangladesh, prevalence is very low across large parts of the country except in Chittagong Hill Tracts (CHT), where risk is especially high among marginalized and hard-to-reach communities . Moderate-to-high risk remains in well-defined areas across Bangladesh . Active case detection with rapid diagnosis tests (RDT) conducted in 2007 revealed that the total malaria prevalence in the CHT was 11.7% with a majority of infections (>90%) caused by Plasmodium falciparum. One particular sub-district, Rajasthali, had a prevalence of 36% and is considered the area most endemic to malaria in CHT . In 2009, another prevalence study carried out in Rajasthali found 11.5% of participants positive with malaria infection despite 2 years of an extensive intervention and control program supported by the Global Fund . Malaria is transmitted in a broad range of eco-epidemiological settings in Bangladesh due to high species diversity and the presence of species displaying different ecological behaviors . Successful malaria control and elimination in Bangladesh can only be achieved if the country can target efficiently the endemic foci. Most of these hotspots are situated in CHT, along the borders of India and Myanmar where 90% of total cases occurs [6–8, 11]. These areas are remote, forested and populated by ethnic minorities living a traditional rural life. The vectors present in this forested area are Anopheles baimai (dirus), Anopheles philippinensis, Anopheles vagus and Anopheles minimus. Malaria risk is related to environmental factors affecting these vectors including altitude, forest, household density, cultivation practices, urbanization, and distance from water bodies [9, 11–13], but the role of these factors to malaria risk has not been quantified.
At present controlling malaria in Bangladesh is based on both preventing infections and on prompt effective treatment of clinical cases. Through the National Malaria Strategic Plan, the malaria and parasitic disease control unit supports malaria prevention and treatment services in 13 endemic districts of Bangladesh.
Predicting the abundance, as well as understanding the risk factors, spatial distribution, and spread of malaria in endemic settings can significantly contribute to local malaria control strategies. A deeper understanding of the role of landscape attributes in the spatial distribution of malaria is crucial so that appropriate local elimination efforts can be developed.
Malaria prevalence in Bangladesh shows a pronounced heterogeneity of transmission . Models predict that unrecognized heterogeneity reduces the efficacy of disease control strategies . While escalating control measures in disease hotspots may be very effective at reducing overall transmission. Recent efforts to develop pragmatic global maps, such as Malaria Atlas Project (MAP), suggests a new era of using maps to define regional populations at risk of malaria, that will guide future global malaria control and the distribution of funds . The most modifiable cause of malaria risk is having two or fewer bed nets; non-modifiable risk factors are associated with sex, age, forest cover, elevation and household density . The relative importance of a risk factor can be estimated using the PAF. The odds ratio and the relative risk can only measure the level of risk associated with the exposure to a risk factor. They never reflect the impact of the factor in a population. PAF is the proportional reduction in population disease that would occur if exposure to a risk factor were reduced over a specified time interval from the population. Many diseases are caused by multiple risk factors, and individual risk factors may interact in their impact on overall risk of disease. PAF estimates can provide some insight for policymakers in planning public health interventions. The goal of this study was to estimate the PAF for different environmental factors for malaria in Rajasthali that will generate key information as a new phase of malaria control in Bangladesh begins with the support of the Global Fund.
Covering a surface area of 145.04 km2, Rajasthali lies between 22°20' to 22°26' north latitude and 92°08' to 92°22' east longitude in CHT. It is a hilly, remote area covered with dense forest. The total population of Rajasthali is 24,097 with 5,322 households. Most households belong to ethnic clans leading a traditional rural life. All households (n = 5,322) in Rajasthali were georeferenced .
Malaria prevalence, socio-economic and environmental data
The malaria prevalence data used in this analysis was collected during a 2009 malaria-prevalence survey in Rajasthali using a two-stage cluster sampling technique. A questionnaire survey was conducted at the same time to obtain demographic information, bed net usage and socio economic status. A full description of the survey has been published . Blood samples from 1,400 (> 0-102 years of age) individuals were tested using rapid diagnosis test. Out of 1,400 samples, 161 (11.5%) were positive for either P. falciparum (158, 11.3%) or Plasmodium vivax (11, 0.8%) . Geometrically corrected, cloud free Landsat 5 Thematic Mapper image was taken on January 23rd, 2010 (Path 136, Row 44) from the United States Geological Survey. Based on the review of existing literature and considering the importance of land cover types, a supervised classification was performed. The study area was classified into six categories (deep water, shallow water, brown open land, bright open land, forest, and grassland/bush) based on the maximum likelihood method for PG-Steamer (Pixoneer Geomatics Inc., Tae-jon). Ground-truth sites and known land-cover were identified from high-resolution Google Earth images. The proportion of forested area within 2 km from each of 1,400 sampled households was calculated using Arc GIS 9.3. Later the proportion of forested area was categorized into tertiles for the analysis. A Shuttle Rader Topographic Mission Digital Elevation Model (SRTM DEM) of 3 arc-second (approximately 90 m) resolution was used to estimate the altitude of each sampled household using ArcGIS 9.3. Local house density was also used for analysis. Data preparation and collation of these items has been described in detail .
This study was reviewed and approved by both the research review committee and ethical review committee of International Center for Diarrhoeal Diseases Research Bangladesh (ICDDR, B). Written consent was taken from study participants or their caretakers .
All statistical analyses were carried out using the statistical software STATA 11 (Stata Corp., College Station, TX). Data included individual-level variables (Age, sex, tribe); household level variables (Household head education and occupation, number of bed nets, long lasting insecticidal net (LLIN) use, bed net treatment and household floor materials) and environmental variables (forest, altitude, household density). The statistical analysis was carried out in two phases using malaria infection status (based on rapid diagnostic test result) as the outcome variable of interest. Firstly, all individual, household and environmental variables were screened using univariate logistic regression for statistical association with malaria infection status based on a liberal p-value of 0.20 in the likelihood-ratio test. Secondly, all variables significant in the screening phase were considered for inclusion through a manual backward stepwise variable selection process in a multivariable logistic regression analysis. The criterion for removal of risk factors was based on statistical considerations using the likelihood ratio test with a significance level of p > 0.05. The variables age, sex, ethnicity, number of bed net, forest, altitude, household density and floor were retained for the final multivariate model
Assessing spatial autocorrelation
The semivariogram was used as a graphical representation of the spatial dependence in the data. The residuals of the final model (age, forest, altitude and household density) were extracted and screened for spatial dependence using semivariograms developed with the package geoR of the statistical software R version 2.9.0 (The R foundation for Statistical Computing, Vienna). The residuals of a model was also tested without the non-significant (p > 0.1) variables (sex, tribes, number of bed net, LLIN, bed net treatment, floor, household head education and occupation), to observe changes to the semivariograms.
Estimation of the population attributable fractions of malaria
The PAF of risk factors included in the final multivariable model was estimated using the AFLOGIT procedure of STATA . PAF estimates the proportional amount of disease risk if all the risk factors were simultaneously eliminated from the population. The adjusted PAF highlights how much of the disease risk is attributed to all the factors included in the model. From the initial descriptive analysis, reference levels for each factor with the lowest risk were determined. That ensured that the PAFs were derived from positive associations with the outcome. A univariate analysis for each risk factor was performed. Then adjusted ORs, PAF and 95% CIs from multivariate models including all the statistically significant risk factors were reported.
The spatial distribution of malaria prevalence in the 1,400 locations surveyed in Rajasthali shows that the distribution of malaria prevalence is heterogeneous across communities (Figure 1) ranging from 4.26% to 45.24% (Table 1). Out of 5322 households, 44.12% were in areas with malaria prevalence of ≥ 10%.
Age, sex, ethnic association, number of bed nets, forest, altitude, household density and floor type proved statistically significant in univariate analyses (Table 2). Only age, altitude, forest and household density proved statistically significant in the multivariate model. Children aged up to 10 years and people living in higher altitude areas had much higher odds of malaria infection than their peers. The odds of malaria increased with increasing elevation. People living around the densest forest (upper third tertile) were at a significantly higher risk of being infected compared to the first and second tertiles. In addition to age, altitude and forest, the odds of malaria also increased with increasing household density. In univariate analysis, lower household density (1-200) proved as the highest risk factor for infection. However, in the multivariate model, higher household density (>1000) produced an increased risk. Adjusted malaria risk with >1,000 households was significantly higher than for (1-200) households density (OR, 3.61 [95% CI, 1.34-9.74]; P = 0.011).
After taking into account the effect of the individual- and survey-level environmental variables, there was little evidence for spatial dependence (Figure 2). For that reason, extra spatial parameters were not fitted to the regression.
Population attributable fraction
The PAFs were estimated for the environmental covariates; altitude, population density and forest (Table 3). Adjusting for all other risk factors, the factors with the highest PAFs were the third tertile of forest cover (47.7%), being a member of the Marma community (23.8%) and living in a place where house density was > 1,000/sq km (19.9%). Most of the Marma villages are located in the western and central regions of Rajasthali. Household density is highest in the central area of Rajasthali where most of the administrative buildings are located. Forest density is higher in the eastern parts of Rajasthali. The second tertile of forest cover accounted for 17.6% of the PAF, Bengali people accounted for 14.5% and household elevation (altitude 51-100 m from sea level) accounted for 11.2% (Table 3). Other communities including Khiang and Chakma (7.3%) , and living at > 100 m altitude (9.3%) with lower household density (201-500) accounted for the 6.5% PAFs (Table 3).
The purpose of this study was to estimate the PAFs for different contributors to malaria risk and to develop a local malaria risk model in a highly endemic area of Bangladesh. The results show that at local spatial scales there was little evidence to include environmental covariates (such as temperature and precipitation) in the models and only the covariates age, forest, altitude and household density were significantly associated with infection. Furthermore, there was little evidence for spatial dependence in semivariograms after the model was developed, suggesting little evidence of substantial unexplained variation that varied predictably over the region. Taken together these results suggest that locally, individual-level factors (e.g. socioeconomic factors; behavioral factors; adherence to preventative measures) are more likely to determine the spatial distribution of malaria infections. The dataset did not include additional individual-level variables describing behavior and further surveys aimed at developing spatial risk models at small spatial scales should endeavor to include this type of information. The inclusion of these individual-level variables may improve the discriminatory ability of the model but unfortunately were not available in the Rajasthali dataset. The greatest PAFs were altitude, household density and forest cover.
Presumably due to lower acquired immunity, children are more vulnerable for malaria infection [7, 18–20]. Results from this study are consistent with the age related effects. The high PAFs for forest coverage and high altitude suggest a critical role of specific malaria vectors. In this region An. dirus, An. minimus, and a diverse fauna of other anopheline species have been reported as main malaria vectors, with An. dirus spp being an efficient vector in forests habitats [21, 22]. Although the vector species have not been sampled in the study area, the higher risk of infection in hilly and forested area may implicate the primary vector species as An. dirus. Despite the small overall elevation distribution in Rajasthali (from 22 m to 359 m above the sea level), this factor may have some contribution to characterize differences of transmission risk, and malaria prevalence rate. Altitude proved one of the key factors responsible for malaria transmission in CHT . Few studies have confirmed that elevation is one of the key factors associated with malaria [23–25]. Study results from Afghanistan showed no transmission in villages at elevations >2,000 m . The altitude in all households of Rajasthali is < 360 m and, thus, not high enough to provide cooler temperatures.
Univariate analyses and PAF further revealed that lower household density was a significant risk factor for malaria infection. This pattern is similar to a study from Ghana that reported higher malaria risk in smaller villages and in outer areas of each village . However, for the multivariate model, the effect of household density was reversed and definitive conclusions could not easily be made. Strong correlations were observed between the forest coverage and house density, and may account for the results in the multivariate model.
The relative importance of risk factors for public health intervention can be estimated using PAF. PAF defines the reduction of incidence and can be achieved if the population had been entirely unexposed compared to its current exposure pattern . The results indicate that altitude, household density and forest were the most important risk factors associated with malaria prevalence in the study area.
These findings are important for targeted intervention and resource allocation. The greatest use of PAFs here has highlighted modifiable risk factors, predicting how much disease can be avoided with their elimination. The non-modifiable factors cannot be eliminated here (e.g., ethnicity). Modifiable risk factors should be considered to prioritize and target public health intervention strategies in CHT. The different ethnic community living in the very remote regions (close to deep forest and in high altitude areas) accounted for the highest PAFs. These results suggest that if living conditions and access to treatment care services can improve, then it would be possible to prevent a significant proportion of malaria incidence in this population. PAF estimates obtained in this study will be useful since the malaria control programme in Bangladesh is enhancing its control efforts and has committed to reduce malaria cases by 60% by 2015 as the baseline year/cases of 2015 .
The increasing rate of malarial infection with greater distances from the village centre could also be due to proximity to mosquito breeding sites. Residences on the village periphery may be located closer to swamps and agricultural fields that result in a greater local mosquito density, although we have no data to indicate this occurred. Socioeconomic status can be protective against malaria if it provides better access to anti-malarial medicines or bed nets, which is the ultimate priority in accordance to malaria control program in Bangladesh. In the central and western areas, it is possible that the more established families with more resources live near the town centre or market . These trends could also be explained by the proximity of health clinics and markets that sell bed nets and anti-malarial medications to the village centre. These findings are pivotal for national control programs and could help policy makers identify high-risk areas that can be targeted first (Figure 1). The calculations (Table 3) are then able to inform program managers if they need to re-distribute or recruit additional health workers to cover 100% of the population in this area. It has important implications in meeting the first key objective of the control programme, which is to effectively diagnosis and treat 100% of estimated malaria cases . Targeting those foci with higher expected prevalence within forested areas will become particularly important for Bangladesh.
This study indicates that prevalence in Rajasthali remains high. Through identification of high-risk malaria zones one may generate new hypotheses regarding the spatial distribution of malaria. It can also be used as advocacy for channeling more funds to conduct operational research in high-risk areas. Targeting of interventions at such fine spatial scales will be helpful for the malaria control programme. Given the social and socio-economic risk factors are identified, more consideration could be given to methods of improving equity of access to interventions among whole communities in the region.
WHO: World malaria report. 2010, World Health Organization, Geneva
Roberts L, Enserink M: Malaria. Did they really say... eradication?. Science. 2007, 318: 1544-1545. 10.1126/science.318.5856.1544.
Feachem R, Sabot O: A new global malaria eradication strategy. Lancet. 2008, 371: 1633-1635. 10.1016/S0140-6736(08)60424-9.
Lines J, Schapira A, Smith T: Tackling malaria today. BMJ. 2008, 337: a869-10.1136/bmj.a869.
Tanner M, Hommel M: Towards malaria elimination--a new thematic series. Malar J. 2010, 9: 24-10.1186/1475-2875-9-24.
Heidi R, Ubydul H, Archie C, Andrew T, Andrew V, Syed Masud A, Akramul I, Rashidul H: Mapping malaria risk in Bangladesh using bayesian geostatistical models. Am J Trop Med Hyg. 2010, 83: 861-867. 10.4269/ajtmh.2010.10-0154.
Haque U, Ahmed SM, Hossain S, Huda M, Hossain A, Alam MS, Mondal D, Khan WA, Khalequzzaman M, Haque R: Malaria prevalence in endemic districts of Bangladesh. PLoS One. 2009, 4: e6737-10.1371/journal.pone.0006737.
Haque U, Magalhaes RJ, Reid HL, Clements AC, Ahmed SM, Islam A, Yamamoto T, Haque R, Glass GE: Spatial prediction of malaria prevalence in an endemic area of Bangladesh. Malar J. 2010, 9: 120-
Haque U, Hashizume M, Glass GE, Dewan AM, Overgaard HJ, Yamamoto T: The role of climate variability in the spread of malaria in Bangladeshi highlands. PLoS One. 2010, 5: e14341-10.1371/journal.pone.0014341.
Lindsay SW, Wilkins HA, Zieler HA, Daly RJ, Petrarca V, Byass P: Ability of Anopheles gambiae mosquitoes to transmit malaria during the dry and wet seasons in an area of irrigated rice cultivation in The Gambia. J Trop Med Hyg. 1991, 94: 313-324.
Haque U, Sunahara T, Hashizume M, Shields T, Yamamoto T, Haque R, Glass GE: Malaria prevalence, risk factors and spatial distribution in a hilly forest area of Bangladesh. PLoS One. 2011, 6: e18908-10.1371/journal.pone.0018908.
Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW: Urbanization, malaria transmission and disease burden in Africa. Nat Rev Microbiol. 2005, 3: 81-90. 10.1038/nrmicro1069.
Bolzoni L, Real L, De Leo G: Transmission heterogeneity and control strategies for infectious disease emergence. PLoS One. 2007, 2: e747-10.1371/journal.pone.0000747.
Hay SI, Okiro EA, Gething PW, Patil AP, Tatem AJ, Guerra CA, Snow RW: Estimating the global clinical burden of Plasmodium falciparum malaria in 2007. PLoS Med. 2010, 7: e1000290-10.1371/journal.pmed.1000290.
Haque U, Hashizume M, Sunahara T, Hossain S, Masud Ahmed S, Haque R, Yamamoto T, Glass GE: Progress and challenges to control malaria in a remote area of Chittagong hill tracts, Bangladesh. Malar J. 2010, 9: 156-10.1186/1475-2875-9-156.
Greenland S, Drescher K: Maximum likelihood estimation of the attributable fraction from logistic models. Biometrics. 1993, 49: 865-872. 10.2307/2532206.
Carneiro I, Roca-Feltrer A, Griffin JT, Smith L, Tanner M, Schellenberg JA, Greenwood B, Schellenberg D: Age-patterns of malaria vary with severity, transmission intensity and seasonality in sub-Saharan Africa: a systematic review and pooled analysis. PLoS One. 2010, 5: e8988-10.1371/journal.pone.0008988.
Maude RJ, Dondorp AM, Faiz MA, Yunus EB, Samad R, Hossain A, Rahman MR: Malaria in southeast Bangladesh: a descriptive study. Bangladesh Med Res Counc Bull. 2008, 34: 87-89.
Yeshiwondim A, Gopal S, Hailemariam AT, Dengela DO, Patel HP: Spatial analysis of malaria incidence at the village level in areas with unstable transmission in Ethiopia. Int J Health Geogr. 2009, 8: 5-10.1186/1476-072X-8-5.
Alam M, Khan M, Chaudhury N, Deloer S, Nazib F, Bangali A, Haque R: Prevalence of anopheline species and their Plasmodium infection status in epidemic-prone border areas of Bangladesh. Malar J. 2010, 9: 15-10.1186/1475-2875-9-15.
Rosenberg R, Maheswary NP: Forest malaria in Bangladesh. II. Transmission by Anopheles dirus. Am J Trop Med Hyg. 1982, 31: 183-191.
Cooper RD, Waterson DG, Frances SP, Beebe NW, Sweeney AW: Speciation and distribution of the members of the Anopheles punctulatus(Diptera: Culicidae) group in Papua New Guinea. J Med Entomol. 2002, 39: 16-27. 10.1603/0022-2585-39.1.16.
Drakeley CJ, Carneiro I, Reyburn H, Malima R, Lusingu JP, Cox J, Theander TG, Nkya WM, Lemnge MM, Riley EM: Altitude-dependent and -independent variations in Plasmodium falciparum prevalence in northeastern Tanzania. J Infect Dis. 2005, 191: 1589-1598. 10.1086/429669.
Salehi M, Mohammad K, Farahani MM, Zeraati H, Nourijelyani K, Zayeri F: Spatial modeling of malaria incidence rates in Sistan and Baluchistan province, Islamic Republic of Iran. Saudi Med J. 2008, 29: 1791-1796.
Brooker S, Leslie T, Kolaczinski K, Mohsen E, Mehboob N, Saleheen S, Khudonazarov J, Freeman T, Clements A, Rowland M, Kolaczinski J: Spatial epidemiology of Plasmodium vivax, Afghanistan. Emerg Infect Dis. 2006, 12: 1600-1602. 10.3201/eid1210.060051.
Kreuels B, Kobbe R, Adjei S, Kreuzberg C, von Reden C, Bater K, Klug S, Busch W, Adjei O, May J: Spatial variation of malaria incidence in young children from a geographically homogeneous area with high endemicity. J Infect Dis. 2008, 197: 85-93. 10.1086/524066.
Rothman Kenneth: Modern Epidemiology. 1988, Los Angles, USA: Lippincott-Raven, 2
This research study was funded by ICDDR, B and its donors which provide unrestricted support to ICDDR, B for its operations and research. Current donors providing unrestricted support include: Government of the People's Republic of Bangladesh; Canadian International Development Agency (CIDA), Embassy of the Kingdom of the Netherlands (EKN), Swedish International Development Cooperation Agency (Sida), and the Department for International Development, UK (DFID). The Johns Hopkins Malaria Research Institute (JHMRI) provided support for one of the authors (GG). We gratefully acknowledge these donors for their support and commitment to ICDDR, B's research efforts.
Conflict of interest
The authors declare that they have no competing interests.
Conceived the study design: UH, RJSM. Data preparation: UH. Data analysis: UH, DM. Prepared the manuscript: UH, GG. Conception, overall scientific management, interpretation of result, and critically reviewed the final report: KK, RH, RJSM, WP, GG. All authors approved the final version of the manuscript. Principal investigator of the project: UH.