ARTICLE |
a Department of International Health, Johns Hopkins School of Public Health, Baltimore, Maryland
b Biomedical Research, AB PRISMA Morona 452, Iquitos, Peru
c Naval Medical Research Center Detachment, Lima, Peru
d Regional Health Department, Loreto, Peru
| ABSTRACT |
|---|
|
|
|---|
PATIENTS AND METHODS. We conducted a 43-month-long prospective, community-based diarrheal disease surveillance in 442 children <72 months of age in the Peruvian Amazon between October 1, 2002, and April 15, 2006.
RESULTS. The incidence of diarrheal disease was 4.38 episodes per child-year. The incidence rate for shigellosis was 0.34 episodes per child-year in children <72 months of age and peaked in children between 12 and 23 months at 0.43 episodes per child-year. Maternal education at or beyond the primary grade level, piped water supply, weight-for-age z score, and improved water-storage practices were the most significant determinants of disease in this community with living conditions comparable to many rural areas in the developing world.
CONCLUSIONS. Children living in this region had a 20-fold higher rate of disease incidence detected by active surveillance as those recently estimated by passive detection. Most symptomatic disease was caused by S flexneri, although the diversity of serotypes will require a multivalent vaccine to have a significant impact on the burden of disease caused by shigellosis. Several other public health disease-control interventions targeted at water source and improved storage, nutritional interventions, and improved maternal education seem to have a greater impact than a univalent S flexneri 2a vaccine.
Key Words: antibiotic resistance bacterial infections diarrhea nutrition vaccines
Abbreviations: PCR—polymerase chain reaction HAZ—height-for-age z score WAZ—weight-for-age z score WHZ—weight-for-age z score GEE—generalized estimating equation CI—confidence interval QIC—Quasi-likelihood Information Criterion IR—incidence rate IRR—Incidence Rate Ratio
Diarrhea causes 1 billion illness episodes and 2.5 million deaths per year of children under the age of 5 in the developing world and China.1 Improved case management of acute dehydrating diarrhea through the use of oral rehydration therapy2 has increased the relative importance of dysentery and persistent diarrhea as the cause of severe and fatal disease.3,4 Shigella is a major cause of diarrhea and mortality in the developing world5 and is the enteric infection that has been most consistently associated with the clinical dysentery syndrome,6,7 prolonged episodes of acute diarrhea, and the development of persistent diarrhea.8 Episodes of shigellosis may result in a protein-losing enteropathy9 and have been shown to be a risk factor for subsequent malnutrition and growth shortfalls.10 Antibiotic resistance is becoming a progressive problem in southeast Asia,11,12 Africa,13 and South America.14,15 For these reasons, the World Health Organization has targeted the development of a vaccine for Shigella as a high priority.
Three key factors are elemental in the adequate control of the morbidity and mortality resulting from shigellosis. The first factor is regular periodic antibiotic-susceptibility testing of isolates to guide local empiric therapy for dysentery. The second factor is intermittent regional monitoring of disease burden with serotype breakdown to guide vaccine development and implementation and identify sites for vaccine trials. The final factor is the determination of risk factors to guide disease-control interventions. We sought to measure the disease burden, the prevalence of antibiotic resistance in Shigella isolates, the distribution of serotypes, and risk factors for shigellosis in the context of rural South America.
| METHODS |
|---|
|
|
|---|
Children were weighed on Salter scales (Salter Housewares Ltd, Tonbridge, England), and their height was measured by using a marked platform with a sliding footboard at 1-month intervals, with each child being measured on the day of their birth. Participating families were visited 3 times weekly by a trained health promoter to document the number, consistency, and characteristics of stool passed over the previous 24-hour period, as well as other gastrointestinal symptoms, fever, anorexia, and treatments given. Reference stool cultures were obtained quarterly to determine asymptomatic carriage rates of Shigella and other enteropathogens.
Diarrhea was defined as an overt change in a child's normal stool pattern characterized by an increase in the frequency to at least 3 unformed stools in a 24-hour period. Dysentery was defined as a diarrheal episode in which a technician observed blood in an unformed stool that coincided with a diarrheal episode on active surveillance. Episodes of diarrhea and dysentery were considered to be separate if separated by 3 diarrhea-free days,16 and the episode was considered attributable to the enteric pathogen if the pathogen was isolated during or within 2 days of the episode. Shigellosis was defined as a diarrheal episode from which Shigella was isolated on stool culture.
Specimens were obtained by field workers as soon as possible after meeting the case definition for diarrhea. One sample was sought for all episodes; however, children who were culture-negative for Campylobacter and Shigella but continued to have diarrhea were asked to provide a sample every fourth day until the episode remitted.
Microbiology
Stool
Fresh voided stools were examined for the presence of blood, mucus, and the state of formation. Testing for occult blood was performed by using hemoccult (SmithKline Diagnostics, Palo Alto, CA) cards, and a fecal smear was stained with methylene blue to measure intestinal inflammation. An aliquot was placed in Cary-Blair medium, refrigerated, transported in a cooler from the field site, and plated on the day of collection. Samples were directly plated on MacConkey, Hektoen, and Blaser's medium. Samples were plated on thiosulfate-citrate bile salts after 6 to 8 hours of pre-enrichment in alkaline peptone water. Colonies displaying an appearance consistent with Shigella by standard biochemical testing were confirmed by agglutination with serogroup specific antisera (Denka-Seiken, Tokyo, Japan). Shigella flexneri isolates were further typed by agglutination with type-specific antisera (Denka-Seiken). Antibiotic sensitivity was tested by using the disk-diffusion method.
Polymerase Chain Reaction
Selected stools frozen at –20°C were processed for invasion plasmid antigen H (IpaH) as per previously described methods.17,18 Polymerase chain reaction (PCR) was run with an initial 5 minutes at 94°C followed by 35 cycles of 94°C for 30 seconds, 55°C for 1 minute, and 72°C for 1 minute, and a final 5 minutes at 72°C. PCR products were electrophoresed on a 2% agarose gel and stained with ethidium bromide and visualized with UV light.
Risk Factors
Risk factors examined in the analysis included demographic and socioeconomic factors, the source and quantity of water available, the use of latrines and the presence of drainage ditches/open sewers in and around the household, and nutritional status. All data except those related to nutritional status and breastfeeding data were derived from the census. Demographic factors included age and gender of the child, maternal age and education (maximum grade level obtained), and people per area of roofed area (people per square meter). Socioeconomic indicators included family and per-capita monthly income, and proxy economic variables regarding the cost of the materials used for housing. Six water-use variables were assessed: (1) the source of water used for drinking (categorized as a piped chlorinated water supply to the household, piped chlorinated water available from a public tap, a deep sealed well, shallow improvised wells, and river/inlet); (2) the volume of the smallest container used for water storage; (3) the number of hours in the day in which piped water was available; (4) total volume of water stored per capita; (5) per-capita volume of water stored with properly sealing tops; and (6) distance from the water source to the kitchen. For sanitation, latrine presence and sharing of a latrine (with >1 household) were recorded. Additional sanitation variables examined included the presence of drainage ditches/open sewers in and around the household and whether they were covered or were reported to overflow into the living area of the family. Nutritional factors examined included the birth weight and the height-for-age z score (HAZ), weight-for-age z score (WAZ), and weight-for-height z score (WHZ). The relationship of breastfeeding to disease incidence was determined by classifying time periods according to months of exclusive breastfeeding (defined as breast milk only with no additional intake of liquids and solids other than medicines), partial breastfeeding (breastfeeding with additional liquids given), mixed breastfeeding (breastfeeding with additional liquids and solids), and weaned.
Analysis
Raw incidence rates (IRs) for diarrheal illness, shigellosis, and dysentery were calculated by dividing the number of incident episodes by the observed time at risk. Generalized estimating equations (GEEs) were used to estimate robust 95% confidence intervals (CIs) for IRs and to conduct risk-factor analysis.19,20 Ninety-five percent CIs were computed by modeling the number of incident episodes per child-year in the study, whereas risk factor analysis modeled the odds of experiencing an episode of Shigella given a hypothesized risk factor.
The multivariate model used GEEs to estimate a logistic regression to model the probability of shigellosis. The model was selected by identifying a maximum set of variables on the basis of examination of descriptive results. Variables were excluded in a stepwise fashion if P values from Wald statistics were >.40 and if the quasi-likelihood information criterion (QIC) statistic indicated improved model fit on variable exclusion.21 Variables examined descriptively that did not enter the maximum model were evaluated for entry after a final full model was chosen. Because of curvilinear relationships, higher power terms for age, normalized anthropometric measures, and time (as continuous variables) were explored and compared with categorical definitions of these variables. The QIC statistic was used to evaluate fit between categorical and continuous definitions. In addition to QIC, goodness of model fit was evaluated by using Horton et al's goodness-of-fit test for logistic regression estimated by GEE.22 Odds ratios are reported, which we interpret in terms of relative risk given our large sample and satisfaction of the rare outcome assumption. Anthropometric data were converted to z scores by using the World Health Organization standards from 1978. Statistical analyses were performed by using SAS 9.2 (SAS Institute, Inc, Cary, NC).
| RESULTS |
|---|
|
|
|---|
Four hundred forty-two children were enrolled in the study and contributed to a total of 914.3 child-years of surveillance. The number of child-years of surveillance is less than the product of participants and the study period because of staggered enrollment, exclusion on the completion of 72 months of life, and because of temporary intermittent out-migrations of the study children from the surveillance area. One child was withdrawn by the mother, and 3 children were lost to external migration. There were no deaths in the cohort during the study period.
Incidence of Diarrhea and Dysentery
In the 334030 child-days of surveillance, there were 4023 episodes of diarrhea detected for an overall incidence of 4.38 episodes per child-year (95% CI: 4.09–4.68). The mean duration of observed episodes was 3.6 days, and the study population spent 4.2% of the total surveillance period with diarrhea. Only 1.7% of episodes persisted beyond 14 days. Of the total number of episodes of diarrhea, samples were available for etiologic analysis in 3756 episodes (93% of all observed episodes). Ninety percent (89.6%) of the study children experienced diarrhea during the period of their active surveillance. Table 1 shows the incidence of all-cause diarrhea, dysentery, shigellosis, and S flexneri according to age group.
|
Incidence and Clinical and Laboratory Characteristics of Shigellosis
Shigella was detected in 8.3% of all diarrheal episodes for which stool specimens were available and had an annual incidence of 0.34 episodes per child-year (95% CI: 0.30–0.38). Shigella incidence was lowest among 0- to 5-month-olds and 6- to 11-month-olds at 0.23 and 0.16 episodes per child-year, respectively. Shigellosis resulted in overt dysentery in only 20.2%, whereas the remainder of cases appeared as acute watery diarrhea. Shigellosis was associated with fecal occult blood in 30.8% of documented episodes and fecal leukocyte detection in 23.4% of episodes, significantly more than in all-cause diarrhea (7.9% and 10.9%; P < .001). The mean duration of an episode of shigellosis was 3.68 days (range: 1–31 days). Shigellosis presented with fever in 73.4% of all episodes and was more likely to be associated with fever than all-cause diarrhea (58.6%; P < .001).
Shigella isolates were also obtained from 3.2% of surveillance stool cultures from children in the absence of diarrheal illness. One hundred three isolates of the 414 isolates analyzed were from asymptomatic children. S flexneri accounted for 67.1% (278) of isolates, Shigella sonnei for 11.8% (49), Shigella boydii for 11.4% (47), Shigella dysenteriae for 2.4%, and other Shigella species that exhibited a biochemical profile of Shigella but did not agglutinate with type specific sera, 7.2% (30). Serotyping of S flexneri isolates revealed 33.1% of all S flexneri and 1.8% of all diarrheal episodes at the community level were caused by S flexneri 2a. S flexneri serotype 3a was present in 19.4%, S flexneri serotype 6 was present in 16.5%, and S flexneri 4a was present in 10.1% of all flexneri isolates. The remaining S flexneri serotypes accounted for <10% of S flexneri isolates and are listed in Table 2.
|
Risk-Factor Analysis
Univariate Analysis
Variation in risk across age categories is shown in Table 1. No differences in incidence of shigellosis were noted according to the gender of the child (P = .58), the age of the child's mother (P = .25), or the presence of a younger child born within 2 years of the sentinel child in the household (P = .49) (see Table 3). Children whose mother's received primary school or less education (<6 years) had a 52% higher (IR: 0.41 [95% CI: 0.36–0.48]) IR than those with more education (IR: 0.27 [95% CI: 0.22–0.32]; P < .005). Greater household crowding trended toward a higher shigellosis incidence, but the differences were not significant (P = .7). Lower monthly per-capita income and poor material used for flooring in the household were both correlated with shigellosis, and children living in houses with wood or bark floors had almost twice the rate of shigellosis as children living on cement floors (P = .01 and < .005, respectively).
|
There was no increased incidence of shigellosis in low birth weight children. The HAZ exhibited a negative trend in respect to the risk of developing shigellosis; however, this trend did not reach statistical significance (P = .25). However, WAZ and WHZ scores were negatively correlated with the incidence of shigellosis. Children in the lowest WAZ quartile had a 77% increase in shigellosis incidence relative to children in the highest WAZ quartile (P = .003), and children in the lowest WHZ quartile had a 61% increase in incidence relative to children in the highest WHZ quartile (P = .01). For children who entered the study under the age of 9 months, breastfeeding status was associated with a 54.3% decrease in the shigellosis incidence in exclusively breastfed children relative to that of weaned children (P = .11) that was not statistically significant. When all children were considered independent of their time of enrollment, breastfeeding was not significantly correlated with the development of shigellosis despite a trend toward a decrease in incidence in exclusively or predominantly breastfed children.
The incidence of disease was evaluated temporally according to season. Transmission intensity was greater in the months of March to October than in November to February. The IRs in the high-transmission seasons were 2.23 times that of the low-transmission seasons (P < .0001). The incidence of shigellosis during low-transmission seasons was 0.21 episodes per child-year (95% CI: 0.16–0.26) and was not significantly different between years in the 3-year study period. The incidence during peak transmission seasons in 2003 (IR: 0.56 [95% CI: 0.46–0.68]) and 2004 (IR: 0.44 [95% CI: 0.34–0.55]) were not significantly different from each other (P = .10) but were significantly higher than the peak in 2005 (IR: 0.29 [95% CI: 0.22–0.38]), largely because of an early decrease in cases in September and October 2005. If the high-transmission season in 2005 was truncated to March through August, the IRs during this season increase to 0.36 episodes per child-year (95% CI: 0.27–0.47), slightly lower than in the 2 previous seasons (P = .04).
Multivariate Model
Descriptive results identified the following maximum set of variables for model selection: age of the child, maternal education, per-capita income, materials used for flooring, water source, sealed water storage, latrine type, coverage of sewers around household, anthropometrics (added separately), breastfeeding status, and seasonality. Model selection procedures resulted in the selection of the subset of variables seen in Table 4. The final model included linear, quadratic, and cubic terms for child's age and WAZ; however, a categorical definition for time defined as March through October of each year was used to capture the seasonal effect. Child's age was centered at 24 months and divided by 12 to ease interpretation of estimates (ie, age = [child's age in months – 12]/24). The final model has a QIC of 2507.62 and Horton's
2 statistic of 8.34 (P = .50), both indicating good fit of the data.
|
6 years, in low-transmission season, who lives in a home with cement floors, uncovered sewers, uses the river as a drinking-water source, and stores water in a properly sealed container. There is a clear curvilinear relationship between child's age and shigellosis risk. Children 12 months of age have approximately half the risk of shigellosis when compared with children over 18 months of age. Figure 1 shows the predicted risk of shigellosis for all ages in children with mothers who did not complete primary schooling versus those who did. Low maternal education is associated with a 34% increase in the odds of developing Shigella (P = .03).
|
|
Household handling of water and sewage also influenced shigellosis risk. Households with an improved water source (household faucet, community faucet, or well) had a 50% decreased risk of shigellosis compared with children living in households that obtained their drinking water from the river (P = .013). Households that did not store water in properly sealed containers had a 30% increased risk of disease (P = .04). The improvised covering of sewers was associated with a 15% reduction in risk, and although the relationship was not significant, the inclusion of this variable in the model improved the model fit according to the QIC statistic.
Seasonality was important in determining risk. Shigellosis risk more than doubled from March through October of each year compared with November through February (P < .0001), even after controlling for all other variables in the model.
| DISCUSSION |
|---|
|
|
|---|
Recently, disease-burden estimates for shigellosis have been increasingly derived from passive detection.23,29,30 Although it is true that this allows for the inclusion of greater areas of study and is more economical, it is clear that significant underreporting occurs. A recent review from Vietnam, Bangladesh, Thailand, Pakistan, and Indonesia estimated an average of 0.013 episodes of shigellosis per child-year in children <5 years old, 20 times lower than seen under active surveillance in this study. Other similar passive surveillance systems have added correction factors (of 5 times the detected number of cases) to estimate true disease burdens.30 Because underreporting of disease is likely to be variable and depends on factors such as transport, health infrastructure, distance to health center, intensity of investigator presence, and local perceptions of disease severity, improved modeling of health care uses are needed before passive surveillance with estimated correction factors can replace community-based longitudinal studies to provide reliable estimates of disease burden.
Despite the fact that Shigella was the most common agent isolated in all children <5 years old with dysentery, the majority of episodes of shigellosis (79.8%) were nondysenteric. In the first year of life, rates of dysentery attributable to Campylobacter were nearly as frequent. Our rates of dysentery caused by Shigella are likely to be decreased by our interventions: all children were treated as soon as Shigella was isolated from a diarrheal stool and on the same day plating was done, which meant that antibiotics were generally being given on illness day 3. This likely prevented the evolution of more severe disease in this population. In addition, our definition of dysentery as visible blood in an unformed stool is more strict that that of other groups, which allow for microscopic criteria to be used.31
Although rates of diarrhea, dysentery, and shigellosis are high, it should be noted that the level of hygiene in this community is similar to that of many populations in the developing world. The community has completed all of the World Health Organization's management of childhood illness guidelines for diarrhea. Standard case management is used and community knowledge and use of oral rehydration solutions are high. Vaccination coverage rates, including measles vaccination, surpass 95%. Furthermore, nearly all children receive near-exclusive breastfeeding until 6 months of age and are not completely weaned until 1 year of age, and an improved water source is available to >94.3% of the population. Despite these interventions, the burden of diarrheal illness, and shigellosis in particular, remain extremely high.
In this cohort, low maternal education, an untreated water source, poor water-storage practices, and being underweight were the key determinants of risk for shigellosis. Improved water storage in particular is an extremely affordable intervention that would seem to be an excellent disease-control strategy in this setting, because children living in homes where water was stored with perfectly fitted lids had a 31.8% decreased incidence of disease and local containers are available at a cost of $3 to $4 US dollars. Maternal education, targeted nutritional interventions to decrease the prevalence of underweight and wasted children, and improved water sources are other interventions that would be expected to have an important impact on the disease burden caused by shigellosis.
The Shigella isolated from children was resistant to the tetracycline, ampicillin, erythromycin, sulfamethoxazole/trimethoprim, and chloramphenicol, the treatments they were most likely to receive for dysentery, in >60% of all cases. Before the antibiotic-resistance data were available, the majority of children received ineffective drug therapy. Although resistance to ciprofloxacin, azithromycin, ceftriaxone, and nalidixic acid was uncommon, these medications are less available in this and other endemic areas. A policy change was necessary to make these antimicrobial agents available on an ongoing basis in this community, a change that has not yet been made on a regional or national level. Empiric strategies for the treatment of the dysentery syndrome need centrally planned periodic microbiologic evaluation to determine their adequacy. Current recommendations that more broadly specify an agent to which local strains are susceptible are not adequate given the lack of information regarding the susceptibility of isolates in most clinical contexts in endemic areas. Regional and national programs should recommend specific agents with microbiologic evidence of in vitro efficacy. Now that generic forms of ciprofloxacin and azithromycin are available and the cost of these more effective treatments has decreased significantly, policy makers should be able to recommend these more efficacious treatments with less reservation.
The use of a vaccine that offered complete protection against S flexneri 2a, 3a, and 6 would protect children against 72.6% of the incident cases of shigellosis caused by S flexneri but only 47.5% of all the cases of shigellosis. The most developed Shigella vaccines concentrate on S flexneri 2a alone, which in this population would have only a limited effect, because even a vaccine with 100% efficacy would prevent 34.1% of the cases of S flexneri but only 22.3% of the cases of Shigella and only 2% of all-cause diarrhea experienced by this cohort. In addition, because 24.2% of the episodes of shigellosis occurred in children under the age of 2 years, unless the vaccine is a conjugate or live-attenuated vaccine administered early in life, the efficacy is likely to be significantly hampered in this and other similar highly endemic settings.
| ACKNOWLEDGMENTS |
|---|
We thank M.A. Single for critical review of an earlier version of the manuscript and S.T. Unt for logistic and field support.
| FOOTNOTES |
|---|
Address correspondence to Margaret Kosek, MD, Johns Hopkins School of Public Health, Department of International Health, 615 N Wolfe St, W 5515, Baltimore, MD 21205. E-mail: mkosek{at}jhsph.edu
The authors have indicated they have no financial relationships relevant to this article to disclose.
The views expressed in this article are those of the authors and do not necessarily reflect the official policy or position of the Department of the Navy, Department of Defense, or US government.
| What's Known on This Subject Since the 1980's the number of prospective cohort studies on diarrheal diseases has fallen markedly. Current disease burden estimates for shigellosis are based on active surveillance data that is often decades old or passive surveillance estimates.
|
| What This Study Adds This study demonstrates high rates of shigellosis and brings into question the accuracy of recently used passive surveillance estimates, provides serotype-specific disease burdens to measure potential disease burden reduction achievable with Shigella vaccines, and compares this with alternative disease-control measures.
|
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||