|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Frontier agricultural settlements in the Amazon Basin favor malaria transmission not only by inducing massive environmental changes, such as deforestation, but also by favoring human clustering close to vector habitats.4 Although the epidemiology of epidemic malaria in early frontier settlements in rural Amazonia was thoroughly described in the mid-1980s,5 minimal follow-up information is available for the recent decades, during which extensive social and ecologic transformations took place.6 To understand how malaria transmission is maintained in the more stable and organized communities originating from frontier settlements, a detailed knowledge of the transmission dynamics of malaria and associated risk factors is necessary.7
Here we describe the current patterns of malaria transmission in one of the largest agricultural settlements of Brazilian Amazonia, the Pedro Peixoto Settlement Project, opened in the state of Acre in the 1980s. The combined analyses of malaria morbidity data and individual and household-level risk factors for a population-based prospective cohort are presented and discussed in relation to prospects for malaria control in this and other settings with similar endemicity.
| MATERIALS AND METHODS |
|
|
|---|
|
|
To characterize the place of residence, we created a continuous covariate, house location, which gives the linear distance in km between each dwelling and an index house located in the sector of earliest ocupation. Briefly, the location of each dwelling along Ramal do Granada and on minor secondary roads was described in relation to an index house (index house A), the first house on the left side of Ramal da Linha 14 at km 14 (Figure 2
); for subjects living along Reserva da Linha 14, the place of residence was calculated as the sum of two linear distances in km: 1) the distance between the subjects house and the index house B (the last house on the left of Ramal da Linha 14 at km 30) and 2) the distance between index houses A and B. New settlers typically occupied land plots close to Iquiri River, corresponding to house locations of 12–15.5 km.
A baseline questionnaire was applied to all study participants to collect demographic, health, and socioeconomic data. Cumulative exposure to malaria was estimated using age and the duration of residence in the study area as proxies, whereas recent malaria exposure was determined from the records of slide-confirmed malaria episodes that were diagnosed in local malaria diagnosis outposts between January 2001 and March 2004.
Information on selected household assets, land ownership, type of building material, and number of inhabitants per room was used to derive a wealth index, from which socioeconomic status was estimated. We combined the following data: 1) ownership of five household assets (gas stove, coach, bicycle, motor vehicle, and cattle); 2) land tenure (yes or no); 3) type of building material (brick walls versus others); and 4) number of inhabitants per room (
1 per room or > 1 per room). Principal component analysis, carried out using the XLSTAT software, version 7.5.2 (Addinsoft, New York, NY), was used to weigh each variable, as described by Filmer and Pritchett.9 The first principal component explained 25.6% of variability and gave the greatest weights to ownership of a couch (0.670) and a motorized vehicle (car or motorcycle) (0.641), and to having a lower number of inhabitants per room (0.574). After the standardized variables were weighted, the highest scores were given to the ownership of a brick house (2.262), a couch (1.040), and a motorized vehicle (0.742). Lowest scores were given to households with no gas stove (–1.237), no land tenure (–1.054), > 1 inhabitant per room (–0.619), and no cattle (–0.614). The scores were summed to give a wealth index for each household (range, –4.871 to 5.409), and these indices were then used to stratify households into quartiles in decreasing order (first quartile, 25% richest).
Baseline malaria prevalence survey. The baseline survey participants five years of age or older were eligible for collection of 5-mL venous blood samples for malaria diagnosis and other laboratory studies, regardless of the presence of any clinical symptoms. Laboratory diagnosis of malaria was based on microscopic examination of thick smear and nested polymerase chain reaction (PCR) amplification of a species–specific segment of the 18S rRNA gene of human malaria parasites.10,11 Thick blood smears were stained with Giemsa and examined for malaria parasites under 700x magnification (at least 100 fields) by two experienced local microscopists. Malaria infections diagnosed by onsite microscopy were treated according to the latest malaria therapy guidelines in Brazil.12 All positive slides and the slides collected from febrile subjects that were declared negative by onsite microscopy were sent for review by an expert microscopist at the National Reference Laboratory of the Ministry of Health of Brazil, Brasília. The DNA templates for PCR amplification were isolated from 200 µL of whole blood using GFX genomic blood DNA purification kits (Amersham Pharmacia Biotech, Piscataway, NJ). Baseline blood samples were additionally used: 1) to investigate human genetic polymorphisms putatively associated with resistance to malaria (manuscript in preparation), 2) to measure hemoglobin and other iron status indicators and diagnose glucose-6-phosphate dehydrogenase deficiency,13 3) to measure antibodies to several arbo-viruses,8 and 4) to measure antibodies to blood-stage malarial antigens.14,15
Malaria surveillance. Between March 2004 and May 2005, the study population was placed under clinical and laboratory surveillance for symptomatic malaria episodes, which included both active and passive case detection.16 For this purpose, the study clinician and a health worker visited the area five times a week (Monday–Friday) and examined all study participants (irrespective of their ages) reporting current or recent fever, headache, or any other malaria symptom. Venous or finger-prick blood samples for malaria diagnosis by microscopy and polymerase chain reaction (PCR) were obtained. As part of surveillance, all patients files kept at the three local malaria diagnosis outposts were daily examined to detect additional malaria episodes among study participants who reported directly to these facilities. Thick blood smears were examined by local microscopists; a portion of them (all positive slides and those collected from febrile subjects that were declared negative by onsite microscopy) were reviewed by an expert microscopist at the National Reference Laboratory of the Ministry of Health of Brazil. The DNA templates for PCR amplification were obtained from venous blood samples as explained previously; those from finger-prick blood samples spotted on FTA Micro Cards were prepared as described by the manufacturer (Whatman, Clifton, NJ). As in other studies conducted in Amazonia,17,18 our PCR often detected additional malaria species that were missed by conventional microscopy. These additional species were also considered in the analysis. A clinical case of malaria (outcome variable) was present when a symptomatic patient, enrolled through either active or passive case detection, was found to have malaria parasites on thick smear. A minimal interval of 28 days between two or more consecutive sample examinations was required to count the latter positive slide as a new malaria episode. In addition, when different species were detected in samples obtained less than 28 days apart, the subject was considered to have a single episode of mixed-species infection. On the basis of the results of onsite microscopy, the treatment was given as described previously. Two additional cross-sectional surveys of the whole study population were performed in September–October 2004 and February–March 2005 to update census data and record temporary absences from the study area that lasted more than one week. These data were considered when computing the number of person-years at risk. Subjects diagnosed with malaria were excluded from the population at risk during the 28 days after initiation of treatment.
Statistical analysis. Exploratory analysis. A database was created with SPSS 13.0 software (SPSS Inc., Chicago, IL). Incidence rates were expressed as the number of cases per 100 person-years or person-months at risk, and their exact Poisson 95% confidence intervals (CIs) were calculated using StatsDirect 2.6.1 software (StatsDirect, Altrincham, Cheshire, UK). Exploratory logistic regression analysis, using R software version 2.6.0,19 examined potential risk factors and confounders, including categorical individual-level variables (gender, migration history, main occupation, recreational fishing, bednet use, recent malaria, and work-related commutation during follow-up), continuous individual-level variables (age, years of schooling, years of residence in the current house, in the study site, and in other malaria-endemic areas, age at the start of malaria-exposure, self-reported number of slide-confirmed malaria episodes, and number of slide-positive malaria episodes recorded between 2001–2004), categorical household-level variables (type of house, type of material used in roof, walls and floor, presence of gaps in the walls, presence of screens in windows and doors, presence of a separate kitchen, presence of indoor animals, type of solid and liquid waste disposal, presence of electricity on the property, wealth index [stratified into quartiles], and sharing the household with a land clearer [land clearing was hypothesized to be an occupation of high risk for malaria]), and continuous household-level variables (number of household inhabitants, years of schooling of the household head, and house location). Separate multiple logistic regression models were fitted to individual- and household-level variables using manual inclusion. Covariates were maintained in subsequent multivariate models, together with the a priori potential confounders (age and gender), if they were associated with the outcome, in exploratory unadjusted analysis, at a level of significance of 20%. The goodness of fit was assessed by either analysis of variance, Akaikes information criterion (AIC) values,20 or odds ratio (OR) changes. Interaction terms were evaluated, and diagnostic tests (Cooks distance21 for influential points and studentized residuals22) were applied using the package car of R 2.6.0 software. A few statistically significant interaction terms were identified but did not improve the models; therefore they were not included in the final models. Important influential points were identified for the variable "occupation," which was then re-stratified and reevaluated; as for the other variables examined, influential points were maintained because overall results did not change significantly when such influential points were excluded from the analysis. Models and outliers were reassembled if needed until the best fit was achieved. Final logistic regression models for individual variables included gender, age, occupation, years of residence in the study site, years of residence outside the study site, the occurrence of recent malaria episodes, recreational fishing, and cohabiting with a land clearer, whereas those for household-level variables included wealth index, house location in the study site, number of household inhabitants, and sharing the house with a land clearer.
Mixed-effects logistic regression.
We used logistic regression models to explore the association between individual and household-level covariates and the occurrence of: 1) at least one episode of either Plasmodium vivax or Plasmodium falciparum clinical malaria during the follow-up, 2) at least one episode of clinical P. vivax malaria, and 3) at least one episode of clinical P. falciparum malaria. Malaria episodes diagnosed during the baseline survey were not considered in this analysis of incident malaria. Because of the nested structure of the data (study subjects are clustered into households), we fitted mixed-effects logistic regression models with multivariate normal random effects using restricted maximum likelihood and numerical integration via adaptive Laplace23 with the lme4 package of R 2.6.0.19 Clustering into households was treated as a random effect in the model (using a variable that identified to which household the subject belonged), whereas all other covariates were entered as fixed-effect variables. Only the variables associated with statistical significance at the 5% level and those that altered the parameter estimate of fixed-effects variables (odds ratio) or the unexplained residual effects at household level (random-effects estimation) by
10% were retained in the final models. Because of the open-cohort design of the study, duration of follow-up (in person-years) was entered as an offset. The following variables remained in final mixed-effects logistic regression models for each outcome: gender, occupation, length of residence in the study area, occurrence of recent slide-confirmed malaria, recreational fishing, wealth index, house location, number of inhabitants in the household, and sharing the house with a land clearer. Because of missing values, only 458 observations remained in the final adjusted models. Given the nonlinear relationship between house location and risk of P. falciparum malaria that was revealed by logistic additive regression (see below), a quadratic term was used to fit the corresponding parametric logistic model. The relationships between duration of residence in the study site and the expected probabilities for each outcome, derived from the final mixed-effects models, were represented graphically using a locally weighted polynomial regression.
Logistic additive models. We applied logistic additive models with integrated smoothness estimation, using penalized regression splines,24 to explore the shape of the association between house location and the risk of clinical malaria during the follow-up. These analyses were performed with the mgcv package of R 2.6.0 software. A similar approach was used to study the association among age, years of residence inside and outside the study area, and the risk of clinical malaria during the follow-up. Because there was extensive overlap between age and time of residence in the study area for many subjects, we created a third variable (years of residence outside the study area) by subtracting the number of years of residence in the study area from the subjects age. Each time-related variable was modeled separately, and its contribution to the outcome was compared using AIC values, the percentage of deviance explained and the P value of the smooth term. Only years of residence in the study area was significantly associated with the risk of clinical malaria (regardless of the species) during the follow-up, whereas age and years outside the study area explained less than 1.12% of the deviance (P > 0.05) and were thus excluded from subsequent analysis (data not shown).
Spatial analysis. The Kulldorff spatial scan statistics was used to test whether malaria episodes were randomly distributed within the study area and to identify significant spatial clusters if present.25 Analysis was made using the Bernoulli model implemented in the version 4.0.3 of the SaTScan software.26 This program creates circular windows that are moved systematically throughout the geographic space to identify significant clusters of infections. The windows are centered on each of the households; with the maximum window size set to include 50% of the households (i.e., the largest possible cluster would encompass 50% of the households). For each location and size of the scanning window, SaTScan performs a likelihood ratio test to evaluate whether infections are more prevalent within that specific circular window as compared with the outside. Separate analyses were made for: 1) at least one infection with any malaria parasite during the follow-up, 2) at least one infection with P. falciparum, and 3) at least one infection with P. vivax. Only subjects who stayed in the study site during the whole follow-up period (between March–April 2004 and May 2005) were included in the spatial analysis. The P values were determined by 100,000 Monte Carlo replications of the data set; and a level of significance of 5% was adopted.
Ethical considerations. The study protocol was approved by the Ethical Review Board of the Institute of Biomedical Sciences of the University of São Paulo, Brazil (318/2002 and 538/2004), and written informed consent was obtained from each adult participant and from the parent or legal guardian of every minor.
| RESULTS |
|
|
|---|
Malaria morbidity during the follow-up.
The study cohort comprised 264 males and 245 females (male:female ratio, 1.08:1) followed for an average of 11.5 months (10 days to 14.5 months), and thus contributing 489.7 person-years of follow-up. A total of 393 study participants (77.2%) remained free of malaria diagnosis during the follow-up. Of the remaining subjects, 71 had a single symptomatic malaria episode, 26 had two episodes, and 19 had between 3 and 6 episodes confirmed by microscopy. No episodes of severe or complicated malaria were diagnosed. Of 195 laboratory-confirmed malaria infections diagnosed during the follow-up, P. vivax accounted for 114 (58.5%), P. falciparum for 48 (24.6%), P. falciparum and P. vivax for 32 (16.4%), and P. vivax and P. malariae for 1 (0.5%). Overall, P. vivax (alone or in mixed-species infections) was found in 147 malaria episodes (incidence, 30/100 person-years at risk, 95% confidence interval [CI] of 25.4–35.3/100 person-years at risk), and P. falciparum (alone or in mixed-species infections) was found in 80 malaria episodes (incidence, 16.3/100 person-years at risk, 95% CI of 12.9–20.3/ 100 person-years at risk), with a 1.8:1 species ratio. The proportions of study participants who remained free of P. vivax and of P. falciparum malaria during the follow-up were similar (82.9% and 86.8%, respectively), although the proportions of subjects with two or more consecutive episodes due to each species differed substantially (6.7% and 2.2%, respectively), possibly because of relatively high rates of P. vivax relapses in this area.27 The temporal distribution of malaria episodes according to species is shown in Figure 3
; the monthly incidence of P. vivax malaria ranged between 0.24 (April 2005) and 6.25 (May 2004) episodes/100 person-months at risk, whereas that of P. falciparum malaria ranged between 0 (November 2004) and 4.24 (June 2004) episodes/100 person-months at risk.
|
|
|
|
|
|
|
| DISCUSSION |
|
|
|---|
At least three findings of the present study have clear implications for designing appropriate interventions against malaria in our site and other settings with similar endemicity. First, malaria morbidity is clearly associated with forest-related activities such as land clearing and, to a lesser extent farming, which together represent the primary occupation of 35.9% study participants. In addition, sharing the house with a land clearer nearly doubled the odds of having malaria, after controlling for house location and several other covariates. Because land clearing also leads to changes in relative vector abundance that may favor malaria transmission in such a modified environment,30 this activity plays a crucial role in maintaining malaria transmission in the entire community. These circumstances are reminiscent of the well-studied forest-malaria phenomenon in Southeast Asia31,32 and call for specific individual- and community-level interventions to reduce the risk of forest-related malaria in frontier settlements of the Amazon Basin.
Second, the probability of having malaria decreases after five years of residence in the settlement, but is not affected by the subjects age. The decline in malaria morbidity with prolonged residence remains statistically significant after controlling for several factors that are potentially associated with increased malaria risk in newcomers, such as the place of residence (newcomers tend to settle in forest fringes, where transmission is more intense and environmental conditions differ from those in areas of earlier settling) and occupation (newcomers are often hired as land clearers before they can earn a living from their crops and cattle ranching). These findings suggest that gradual acquisition of clinical immunity to malaria exists under conditions of low endemicity that prevail in most communities originated from frontier settlements in rural Amazonia. Despite the relatively low levels of local malaria transmission, malaria-exposed subjects often harbor subclinical infections with very low parasite loads, most of them detected only by PCR.33–36 Our data indicate that asymptomatic carriers may constitute a significant reservoir of malaria not only in isolated and sparsely populated riverine communities of native Amazonians,33 but also in more densely populated agricultural settlements.
Third, we showed a significant spatial clustering of malaria risk in the study site. These results suggest that malaria control measures should be spatially targeted for maximum effect.37 Significantly, two-thirds of subjects who had P. falciparum malaria episode(s) diagnosed during the follow-up lived in one-fifth of the dwellings, all situated in the more recently occupied areas of the study site. The spatial clustering of malaria revealed by this and other studies6,7,37 indicates that selective house spraying and other household–specific interventions are possibly efficient malaria control strategies in this and other settings with similar endemicity.
Subjects with malaria episodes shortly before the baseline survey were at an increased risk for subsequent malaria infections during the follow-up, most likely a result of the continuous exposure to some of the risk factors identified in our study. Interestingly, some predictors of malaria risk affect predominantly newcomers to the area in a synergistic way. First, the newcomers, before being able to earn the living from their own land, are often hired by local farmers to work in high-risk activities such as land clearing. Second, vacant land plots available for new settlers are often covered with native rain forest and are situated close to protected forest reserves. Extensive deforestation will thus be required before this land can be used for cultivation or cattle ranching. Finally, because most newcomers are migrants from malaria-free regions, they lack acquired immunity to malaria parasites. As a result, the continuous arrival of nonimmune settlers to occupy recently demarcated plots in the periphery of more structured colonization projects results in further environmental change that maintain malaria transmission in a number of communities in rural Amazonia.
Only two malaria parasite species, P. vivax and P. falciparum, are of public health importance in our study site. The absence of P. malariae in the baseline survey and the low incidence of P. malariae malaria during the follow-up (this parasite was found in only 1 of 195 infections) contrast with the high proportion of infections with this species (9.4–11.9%) found in other frontier areas of Brazil.17,18 The reasons for these differences are unclear given that similar diagnostic methods were used in all studies. However, the proportions of mixed-species infections diagnosed in the baseline survey (15 of 60 PCR-diagnosed infections, 25%) and during the follow-up (33 of 195 infections, 16.9%) are within the previously reported ranges (determined by PCR) for native (5.4–36.1%)33 and for migrant populations (25.7–30.2%)17,18 across the Amazon Basin of Brazil.
Assessing malaria risk factors in frontier settlements is particularly challenging because of several factors: 1) there is no clear-cut correlate of cumulative exposure to malaria, as age is in rural Africa; 2) individuals cluster into households where several risk factors are shared; and 3) rural settlements usually consist of straight roads with linear or nonlinear gradients of time- and space-dependent covariates. These features were carefully considered in the present analysis, leading to the identification of several independent predictors of malaria risk in rural Amazonia that are potential targets for control measures. Further analyses of data generated by this ongoing cohort study are expected to provide additional insights into the phenomenon of frontier malaria in the Amazon Basin.
Received March 14, 2008. Accepted for publication July 7, 2008.
Acknowledgments: The authors thank the inhabitants of Ramal do Granada for their participation in the study; Sebastião Bocalom Rodrigues, Damaris de Oliveira, and Nésio M. Carvalho (Municipal Government of Acrelândia), Raimundo A. Costa and the malaria control teams in Granada and Acrelandia, for their logistic support, Adamílson Luís de Souza and Carlos E. Cavasini for help in fieldwork, Francisco das Chagas O. Luz (Ministry of Health, Brasília, Brazil) for reviewing all malaria slides, Marília Sá Carvalho and Oswaldo Gonçalves Cruz (Program of Scientific Computation, Oswaldo Cruz Foundation) for their helpful advice in the statistical analysis, Cassiano P. Nunes, for artwork, and Tatiana Havryliuk for a critical reading of the manuscript.
Financial support: CNPq (470067/2004-7) and FAPESP (03/09719-6 and 05/51988-0) to M.U.F.; scholarships from FAPESP (M.dS.N. and C.J.), and CNPq (N.S.dS. and M.U.F.) are also acknowledged.
* Address correspondence to Mônica da Silva-Nunes, Department of Parasitology, Institute of Biomedical Sciences, University of São Paulo, Av. Prof. Lineu Prestes 1374, São Paulo, São Paulo, Brazil, 05508-900. E-mail: msnunes1{at}yahoo.com.br ![]()
Authors addresses: Mônica da Silva-Nunes, Natal S. da Silva, Camila Juncansen, and Marcelo U. Ferreira, Department of Parasitology, Institute of Biomedical Sciences, University of São Paulo, Av. Prof. Lineu Prestes 1374, São Paulo, São Paulo, Brazil 05508-900, E-mails: msnunes1{at}yahoo.com.br, natalss{at}gmail.com, camilacj{at}gmail.com, and muferrei{at}usp.br. Cláudia T. Codeço, Program of Scientific Computation, Oswaldo Cruz Foundation, Av. Brasil 4365, Rio de Janeiro, Rio de Janeiro, Brazil 21045-900, E-mail: codeco{at}fiocruz.br. Rosely S. Malafronte, Laboratory of Protozoology, Institute of Tropical Medicine of São Paulo, Av. Dr. Enéas de Carvalho Aguiar 470, Cerqueira César, São Paulo, Sã o Paulo, Brazil 05403-000, E-mail: rmalafronte{at}usp.br. Pascoal T. Muniz, Department of Health Sciences, Federal University of Acre, BR 364 km 4, Rio Branco, Acre, Brazil 69915-000, E-mail: pascoal{at}ufac.br.
Reprint requests: Mônica da Silva-Nunes or Marcelo U. Ferreira, Department of Parasitology, Institute of Biomedical Sciences, University of São Paulo, Av. Prof. Lineu Prestes 1374, 05508-900 São Paulo (SP), Brazil, E-mails: msnunes1{at}yahoo.com.br and muferrei{at}usp.br.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
F. A. Souza-Silva, M. d. Silva-Nunes, B. A. M. Sanchez, I. P. Ceravolo, R. S. Malafronte, C. F. A. Brito, M. U. Ferreira, and L. H. Carvalho Naturally Acquired Antibodies to Plasmodium vivax Duffy Binding Protein (DBP) in Rural Brazilian Amazon Am J Trop Med Hyg, February 1, 2010; 82(2): 185 - 193. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Orjuela-Sanchez, N. S. da Silva, M. da Silva-Nunes, and M. U. Ferreira Recurrent Parasitemias and Population Dynamics of Plasmodium vivax Polymorphisms in Rural Amazonia Am J Trop Med Hyg, December 1, 2009; 81(6): 961 - 968. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. U. Ferreira, R. M. Hiramoto, D. P. Aureliano, M. da Silva-Nunes, N. S. da Silva, R. S. Malafronte, and P. T. Muniz A Community-based Survey of Human Toxoplasmosis in Rural Amazonia: Seroprevalence, Seroconversion Rate, and Associated Risk Factors Am J Trop Med Hyg, July 1, 2009; 81(1): 171 - 176. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |