|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
The geographic distribution of filarial infection is nonuniform.3,4 Since the factors influencing this distribution have not been identified, elimination programs initially conduct nationwide blood surveys to construct national LF prevalence maps.5 However, certification of elimination requires intensive screening for infection in areas at high risk for LF transmission, and the cost of refining the initial national surveys to precisely demarcate at-risk areas and identify regions for certification sampling may be prohibitive.6 Consequently, there is a need to ascertain which environmental, geographic, and demographic factors determine the distribution of filarial infection and develop a tool for prediction of community risk of filarial infection using easily-obtainable information.
In 1999, the Centers for Disease Control and Prevention (CDC) initiated an LF elimination demonstration program in Leogane Commune, Haiti, in collaboration with Hôpital Sainte Croix, the main health facility in the Commune, and Notre Dame University (Notre Dame, IN). Leogane Commune is located on Haitis southern peninsula and covers an area of approximately 400 km2. Home to 150,000200,000 people, the Commune is divided into 13 sections (Hôpital Sainte Croix, unpublished data). The three sections in the coastal plain (1, 2, and 3) are the most populous; the remaining 10 sections are mountainous and more sparsely populated. Initial program activities included an assessment of the geographic distribution of W. bancrofti infection in the Commune. We describe this assessment and evaluate the relationships between infection and potential geographic, socioeconomic, and demographic risk factors.
| METHODS |
|
|
|---|
The survey protocol and consent form were reviewed and approved by both the CDC and Hôpital Sainte Croix institutional review boards (IRBs). The IRB of Emory University reviewed and approved the proposed analyses as secondary analyses of pre-existing data.
Grid sampling.
Twenty-five schools were selected using a grid sampling scheme. Using ArcView software version 3.2 (Environmental Systems Research Institute, Inc., Redlands, CA), a grid of rectangles of equal size was overlaid on a map of the Commune (Figure 1
), excluding remote mountain areas. A point within each grid rectangle was then chosen at random and its longitude and latitude were read from the map. Using a handheld global positioning system (GPS) unit (Magellan unit; Thales Navigation, Inc., Santa Clara, CA) to find each random point in the field, an advance team located the primary school(s) nearest that point (Figure 1
). We aimed to enroll at least 100 children per grid rectangle. If a single school with 100 or more eligible children did not exist near the random point, we selected multiple schools whose combined eligible children numbered more than 100. On occasion, in sparsely populated areas with few schools, enrollment fell short of 100 students. Once schools were identified, the team met with each headmaster to discuss the purpose of the school prevalence survey and encourage him to allow his school to participate. No school declined to participate.
|
Mountain school selection.
Many mountainous areas were remote and rugged, and could only be accessed on foot. The grid sampling strategy could not be executed in these areas because locating random points and nearby schools on foot proved unfeasible. To obtain good coverage of the mountainous portion of the Commune, the advance team met with community health workers (CHWs) from sections 610 and 1315 to enlist their help in identifying schools in their sections. All 17 schools identified by the CHWs present at the meeting participated in the prevalence survey (Figure 1
).
Leogane town and immediate surroundings.
Prior to establishing the LF elimination program in Leogane Commune, CDC conducted a school-based randomized control trial to examine the efficacy of diethylcarbamazine (DEC) and albendazole in treating W. bancrofti and intestinal helminth infections in 511-year-old children living in and around the town of Leogane (Beach M, unpublished data). Ten schools, five from the town of Leogane proper and five from surrounding communities (within 5 km of town), participated in the trial; 1,329 children were enrolled. Of these, 291 children from nine schools were evaluated for W. bancrofti infection at baseline using the same test used in our prevalence survey. Their test results were used as the school prevalence data for the town of Leogane and neighboring communities (Figure 1
).
Antigen testing and individual-level data collection. Each childs birth date (if known), age, sex, and home community were recorded during the consent process before infection status assessment. W. bancrofti infection status was evaluated using an immunochromatographic test (ICT) (AMRAD, Sydney, New South Wales, Australia), a qualitative test for the presence of circulating W. bancrofti adult worm antigen. Capillary (fingerstick) blood (100 µL) was collected from each child and applied to an ICT card. The test result was read as antigen positive (infected) or antigen negative (uninfected) after 15 minutes. Infected children were treated with a single tablet (50 mg) of DEC and all children were treated with albendazole.
School-level data collection. The longitude and latitude of each school were recorded using GPS (TSC1 Asset Surveyor with Pro XRS receiver and Pathfinder Office software, Trimble Navigation, Ltd., Sunnyvale, CA; and a Garmin eMap unit with Map Source software, Garmin International, Inc., Olathe, KS). Using these GPS-derived school coordinates, we read the altitude of each school in meters above sea level from a detailed topographic map of the Commune (scale = 1:50,000). We also used the geographic coordinates to categorize the schools by administrative section and topographic zone; in the latter case, schools were characterized as being located on the plains, in the foothills, or in the mountains based on school location and knowledge of the surrounding terrain.
Each headmaster was asked to provide information about 1) the schools tuition, 2) the in-school nutrition program for students (if such a program existed), and 3) the communities served by the school.
Analytic methods. All analyses were performed using SAS version 8.2 (SAS Institute, Inc., Cary, NC) and all maps were constructed using ArcView. Prior to modeling, descriptive analyses were performed on the independent and dependent (infection status) variables. Logistic generalized linear mixed modeling was conducted using SASs GLIMMIX macro.79 All models included random school effects to account for within-school correlation (correlation between students attending the same school). We first fit single-variable models to identify possible associations between infection status and individual independent variables.
All independent variables were included in the initial multivariate model. However, topographic zone, section, and altitude were highly collinear; including two (or all three) of these variables in the same model caused parameter estimate variances to become large and rendered the model unstable. As a result, although each variable provided important, different, and complimentary information, only one of the three could be included at a time in subsequent multivariate models.
We considered four multivariate mixed models containing topographic zone, section, quartiles of altitude, or alternate altitude ranges, respectively. All four models contained age, sex, tuition, nutrition program and random schools effects; all models also initially contained an age x sex interaction term. The age x sex interaction was non-significant in all models (all P > 0.70) and was therefore dropped, leaving only main effects models of the form
![]() |
where j = 1, . . . , ni, with ni being the number of children in the ith school; i = 1, . . . , 57;
ij denotes the probability of infection for the jth child from the ith school; geography denotes topographic zone (k = 1,2), section (k = 1,2 . . . 8), or altitude (k = 1,2,3 or 1,2,3,4); and
i is the random effect of the ith school. The level of significance (
) was set at 0.05 for all multivariate models. The results from these mixed models should be interpreted given the random effects in the model; in the interests of brevity, this condition is implied rather than explicitly stated in the sections that follow.
| RESULTS |
|
|
|---|
|
25 meters above sea level; more than half were located
75 meters above sea level. Thirteen (22.8%) schools were located > 400 meters above sea level. The distributions of age and sex were similar in all three topographic zones. In contrast, the distribution of tuition varied by topographic zone. Plains schools accounted for 70% of schools in the highest quartile of tuition; 9 of 22 plains schools (40.9%), compared with 2 of 21 mountain schools (9.5%) and 2 of 13 foothills schools (15.4%), were in the highest tuition quartile. Of the 14 schools in the lowest quartile of tuition, seven were mountain schools and six were plains schools.
The distribution of nutrition programs also varied by tuition and topographic zone. Of 40 schools in the middle and highest quartiles of tuition, 30 (75%) did not have in-school nutrition programs. In contrast, more than 70% (10 of 14) of lowest quartile schools had such programs. Forty-eight percent (10 of 21) of plains schools had in-school nutrition programs, whereas only 16.7% (2 of 12) of foothills schools and 31.8% (7 of 22) of mountain schools had nutrition programs.
Distribution of W. bancrofti infection and single-variable model results.
Five hundred thirty-three children (16.1%) tested positive for W. bancrofti infection. W. bancrofti infection prevalence by student home community is shown in Figure 2
. Table 2
summarizes the distribution of infection with respect to student and school characteristics and presents the single-variable model results. We observed modest associations between age, sex, and infection status. Children attending schools in the middle quartiles of tuition were less likely to be infected than children attending lowest quartile schools, whereas children attending highest quartile schools were almost three times as likely to be infected as children from lowest quartile schools. Seventy-nine percent (421 of 533) of infected children attended plains schools; children attending schools in the foothills and plains were 4 and 29 times as likely to be infected, respectively, as their counterparts in the mountains. Sections 1 through 3, which encompass the entire plains portion of the Commune as well as some foothills areas, had a combined infection prevalence of 30.9%. These three sections accounted for 92.2% of infections. In contrast, predominantly mountainous sections 615 had a combined prevalence of 2.4%. Mountain sections 7 and 14 and mountain/foothills section 15 accounted for 78.6% of the infections in this group. Children attending schools in sections 1, 2, 3, 7, 14, and 15 were more likely to be infected than children attending schools in sections 10 and 13, two of the most remote and rugged sections in the Commune (prevalence odds ratio [POR] range = 3.2169.69). In contrast, children attending schools in sections 6, 8, 9, 11, and 12 were no more likely to be infected than children in the reference group. The prevalence of infection decreased dramatically with increasing altitude; we observed similar patterns with quartiles of altitude and absolute altitude levels.
|
|
|
The association of infection with altitude also remained unattenuated in the multivariate analyses, for both quartiles of altitude and absolute altitude levels (Table 3
).
Age continued to be associated with infection status in the multivariate analyses, regardless of the geographic variable included in the model. Seven and eight-year-old and 911-year-old children were 3031% (P value range = 0.0390.048) and 2528% (P value range = 0.0650.095) more likely, respectively, to be infected than five and six-year-old children. In contrast, sex was not associated with infection status after controlling for the other covariates.
The strong associations between infection status and tuition suggested by the single-variable models were weaker after controlling for other covariates and all confidence intervals included the possibility of no association. Depending on the geographic variable included in the model, children attending schools in the middle quartiles of tuition were as likely or slightly less likely to be infected than children attending schools in the lowest quartile of tuition. In contrast, children attending highest tuition quartile schools were up to 70% more likely to be infected than children attending lowest tuition quartile schools.
The association between nutrition program and infection status was also weakened after controlling for the other covariates and regardless of the geographic variable included in the model, confidence intervals included the null (POR = 1). Children attending schools with nutrition programs were up to 86% more likely to be infected than children from schools with no programs.
| DISCUSSION |
|
|
|---|
Ideally, we would have constructed a single model containing all three geographic variables, but the near collinearity between them prevented this. Alone, each variable provided important information, but no single variable provided a complete picture of the geographic distribution of W. bancrofti infection in Leogane Commune. The altitude findings allow for coarse identification of areas at risk for infection (areas lying below 70 meters above sea level), but some low-prevalence communities were located at lower altitudes, while the intermediate- and high-prevalence communities in sections 7 and 14 were located between 360 and 450 meters above sea level. Without the section variable, the foothill and mountain "hot spots" in sections 7, 14, and 15 might have gone unnoticed, masked by the overall low prevalence of infection observed in the mountains. Without topographic zone, the marked plains-foothills-mountains gradient of infection prevalence would not have been identified. Therefore, multiple geographic descriptors should be considered when identifying areas at risk for W. bancrofti infection.
The geographic variable results contribute in several ways to our understanding of the geographic distribution of W. bancrofti infection in Haiti, and potentially in other countries where Culex quinquefasciatus transmits LF, and may help to guide LF elimination program activities. In areas where anophelines are the predominant vectors, transmission of W. bancrofti infection at higher altitudes is assumed. Onapa and others have found W. bancrofti infection prevalences ranging from 18.3% to 30.1% in Ugandan communities located 1,0001,100 meters above sea level.10 In areas with non-anopheline vectors, however, it is commonly accepted that W. bancrofti infection prevalence decreases sharply with increasing altitude. However, a recent study in Nepal, where Cx. quinquefasciatus is the main LF vector, found infection prevalences of 13.816.5% at various altitudes up to 900 meters above sea level and 9.6% between 900 meters and 1,500 meters.11 Our results indicate that in Haiti, the force of infection is greatest (and most relevant to filariasis elimination program activities) at much lower altitudes; the odds of infection drop by at least one order of magnitude above 70 meters. That said, the section results offer a caveat: although individuals living on floodplains at low altitude are much more likely to be infected than individuals living at higher altitudes, higher altitude communities are not universally exempt from infection. Filariasis elimination programs in Haiti and other areas where Cx. quinquefasciatus predominates should recognize that there may be substantial pockets of infection and potentially even transmission of infection at higher altitudes, and extend infection mapping activities into what have been assumed to be non-LF altitudes. The current results also suggest that foothills communities serve as transition areas between areas of high and low infection risk.
The characteristics, if any, shared by at-risk plains, foothills, and mountain communities remain unclear. They may share, or lack, certain landscape traits, perhaps sharing microenvironments that favor mosquito breeding, survival, and transmission of infection. Alternatively, at-risk communities may be similar with respect to the exposure opportunities of their populations; active transmission of infection might not be occurring in "unusually located" high prevalence communities themselves, but substantial proportions of their populations might be exposed elsewhere. Section 7, although in the mountains, overlooks a suburb of Port au Prince with a W. bancrofti prevalence of 7.5% (Ministère de la Santé Publique et de la Population, Haiti, unpublished data). Due to its geographic location, there is more movement between this section and a source of infection than might be expected for a mountain community. Similarly, the presence of good roads into section 14 and through section 15 results in more population movement between these sections and the plains than exists in more remote mountain settings.
Lymphatic filariasis has traditionally been considered a disease of poverty, inadequate sanitation, and underdevelopment,1214 although the association between socioeconomic status (SES) and W. bancrofti infection has never been formally evaluated. We included tuition in the models as an indicator of the relative SES of a childs family. School tuition may not be a precise indicator of family SES, since families value education differently; for example, some very poor families make financial sacrifices to send their children to good schools. However, to the extent that tuition is a valid marker for SES and our point estimates are stable, our results suggest that higher SES may not confer protection against infection with W. bancrofti. In this setting, children with the highest SES were, at best, as likely as children with the lowest SES to be infected. Higher SES residents of Leogane Commune may be more mobile than the less affluent, traveling frequently to urban and periurban plains areas (areas in which the force of infection and population density are highest) or settling in these more desirable areas, thus increasing their opportunity for exposure. Alternatively, all residents of a given area share the same environmental risk profile, and in this part of Haiti, higher SES does not appear to be associated with increased use of physical means of protecting oneself from vector-borne infections, such as window screens or bed nets, as it might be in other settings. Therefore, there may be no difference between the highest and lowest SES groups in terms of physical protection from W. bancrofti infection. Regardless, the tuition findings are of value to elimination programs in terms of their implications for community education and motivation and drug coverage. If SES (as defined locally) does not confer protection against W. bancrofti infection, then everyone, regardless of SES, must participate in LF elimination programs to ensure program success.
Since data collection involved only a brief encounter in the school setting with each parent-child pair, very little personal data could be collected. Consequently, there may be confounding by unmeasured covariates, particularly by individual-level characteristics such as mosquito bite prevention behaviors (e.g., bed net use, personal and household insecticide use), location of residence (urban versus rural), and home community population density. In addition, because there was no way to obtain geographic coordinates for each childs place of residence, school location was used as a proxy for home location and values for the geographic location variables were assigned based on the location of the childs school. However, risk of infection probably depends primarily on home location, since Cx. quinquefasciatus, the mosquito vector for W. bancrofti in Haiti, feeds at night. A review of available home community data suggested that most children lived near their schools. Nevertheless, if a childs home and school geographies differed, there could have been non-differential misclassification of one or more geographic variables, which would bias the relevant PORs towards the null.
Our sampling methodology and the participation of all selected schools provided high resolution geographic coverage of Leogane Commune, yielding good estimates of prevalence and allowing us to discern local variations in the distribution of W. bancrofti infection. Given the clear spatial pattern suggested by the prevalence data, future work will focus on formally describing the spatial correlation between outcomes and accounting for this correlation when modeling infection status. In addition, we aim to determine whether at-risk locations share common landscape characteristics by assessing the relationship between landscape traits, region, and W. bancrofti infection using remote sensing data and these school data.
Received September 22, 2003. Accepted for publication November 25, 2003.
Acknowledgments: We thank Joyanna Wendt and Jeanne Radday (Leogane Commune LF elimination program coordinators at the time this work was undertaken), Dr. Jack Lafontant (Director, Hôpital Sainte Croix), and Dr. Thomas Streit for their support during the school survey data collection. Dr. David Molyneuxs support and enthusiasm for the project were also much appreciated. We are very grateful to Allen Hightower for his invaluable advice and technical assistance, without which the GPS mapping work could not have been accomplished. We also thank Dr. Brian Stephens for the two months of hard work required to create the GPS base map of the Commune, and Renn Doyle for digitizing the only existing map of the Commune. We are indebted to the following elimination program staff members for their dedication to the program and for their willingness to travel the length and breadth of the Commune during data collection activities: Conus Brevaus, Jean Marc Brissau, Dardith Desire, Erubens Elaus, Shiler Emile, Fred Francis, Jean Steven Hilaire, Jean Willy Hilaire, Jesnel Jorseme, Michelet Leriche, Jean Clausel Louis, Rose Guertha Louis-Charles, Wesly Pierre, Osnel Sagaille, and Kethlie Toussaint.
Financial support: This work was supported by the Department for International Development, United Kingdom, through the Liverpool School of Tropical Medicine Lymphatic Filariasis Support Centre and by the CDC Emerging Infections Program.
Authors addresses: Heather A. Boyd, Department of Epidemiology, Rollins School of Public Health, Emory University, Room 468, 1518 Clifton Road, Atlanta, GA 30322 and Division of Parasitic Diseases, Centers for Disease Control and Prevention, Chamblee Facility, Mailstop F-22, Building 24, Room 1006, 4770 Buford Highway, Atlanta, GA 30341-3717, Telephone: 770-488-7595, Fax: 770-488-4465, E-mail: heb6{at}cdc.gov. Lance A. Waller, Department of Biostatistics, Rollins School of Public Health, Emory University, Room 326, 1518 Clifton Road, Atlanta, GA 30322, Telephone: 404-727-1057, Fax: 404-727-1370, E-mail: lwaller{at}sph.emory.edu. W. Dana Flanders, Department of Epidemiology, Rollins School of Public Health, Emory University, Room 464, 1518 Clifton Road, Atlanta, GA 30322, Telephone: 404-727-8716, Fax: 404-727-8737, E-mail: wflande{at}sph.emory.edu. Michael J. Beach, Centers for Disease Control and Prevention, Chamblee Facility, Building 102, Room 1411, Mailstop F-22, 4770 Buford Highway, Atlanta, GA 30341-3717, Telephone: 770-488-7763, Fax: 770-488-7761, E-mail: mjb3{at}cdc.gov. J. Sony Sivilus and Rodrigue Lovince, Hôpital Sainte Croix, Leogâne, Haiti, Telephone: 509-235-1044; E-mails: jeanssivilus{at}yahoo.fr and Rodrigue.Lovince.1{at}nd.edu. Patrick J. Lammie, Centers for Disease Control and Prevention, Chamblee Facility, Building 23, Room 1021, Mailstop F-22, 4770 Buford Highway, Atlanta, GA 30341-3717, Telephone: 770-488-4054, Fax: 770- 488-4108, E-mail: pjl1{at}cdc.gov. David G. Addiss, Centers for Disease Control and Prevention, Chamblee Facility, Building 24, Room 1007, Mailstop F-22, 4770 Buford Highway, Atlanta, GA 30341-3717, Telephone: 770-488-7770, Fax: 770-488-4465, E-mail: dga1{at}cdc.gov.
Reprint requests: David G. Addiss, Centers for Disease Control and Prevention, Chamblee Facility, Mailstop F-22, 4770 Buford Highway, Atlanta, GA 30341-3717.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
M. B. DE ROCHARS, S. KANJILAL, A. N. DIRENY, J. RADDAY, J. G. LAFONTANT, E. MATHIEU, R. D. RHEINGANS, A. C. HADDIX, T. G. STREIT, M. J. BEACH, et al. THE LEOGANE, HAITI DEMONSTRATION PROJECT: DECREASED MICROFILAREMIA AND PROGRAM COSTS AFTER THREE YEARS OF MASS DRUG ADMINISTRATION Am J Trop Med Hyg, November 1, 2005; 73(5): 888 - 894. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. V. E. B. DE ROCHARS, M. D. MILORD, Y. ST. JEAN, A. M. DESORMEAUX, J. J. DORVIL, J. G. LAFONTANT, D. G. ADDISS, and T. G. STREIT GEOGRAPHIC DISTRIBUTION OF LYMPHATIC FILARIASIS IN HAITI Am J Trop Med Hyg, November 1, 2004; 71(5): 598 - 601. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |