Subscriber access provided by BUFFALO STATE
Characterization of Natural and Affected Environments
Predictors of Enteric Pathogens in the Domestic Environment from Human and Animal Sources in Rural Bangladesh Erica R. Fuhrmeister, Ayse Ercumen, Amy J Pickering, Kaitlyn M. Jeanis, Mahaa Ahmed, Sara Brown, Benjamin F. Arnold, Alan E. Hubbard, Mahfuja Alam, Debashis Sen, Sharmin Islam, Mir Himayet Kabir, Laura H. Kwong, Mahfuza Islam, Leanne Unicomb, Mahbubur Rahman, Alexandria B Boehm, Stephen P. Luby, John M. Colford, and Kara L. Nelson Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b07192 • Publication Date (Web): 29 Jul 2019 Downloaded from pubs.acs.org on July 30, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 35
Environmental Science & Technology
ACS Paragon Plus Environment
Environmental Science & Technology
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
Page 2 of 35
Predictors of Enteric Pathogens in the Domestic Environment from Human and Animal Sources in Rural Bangladesh Erica R. Fuhrmeister1, Ayse Ercumen2,3, Amy J. Pickering4, Kaitlyn M. Jeanis 1, Mahaa Ahmed 1, Sara Brown1, Benjamin F. Arnold2, Alan E. Hubbard2, Mahfuja Alam5, Debashis Sen5, Sharmin Islam5, Mir Himayet Kabir5, Laura H. Kwong6, Mahfuza Islam 5, Leanne Unicomb5, Mahbubur Rahman5, Alexandria B. Boehm6, Stephen P. Luby7, John M. Colford Jr. 2, Kara L. Nelson1* 1Department
of Civil and Environmental Engineering, University of California, Berkeley, California, 94720, United States 2School of Public Health, University of California, Berkeley, California 94720, Unites States 3Department of Forestry and Environmental Resources, North Carolina State University, Raleigh, North Carolina, United States 4Civil and Environmental Engineering, Tufts University, Medford, Massachusetts, 02153, United States 5Infectious Disease Division, International Centre for Diarrhoeal Disease Research Bangladesh, Dhaka, 1212, Bangladesh 6Department of Civil and Environmental Engineering, Stanford University, Stanford, California, 94305, United States 7Woods Institute for the Environment, Stanford University, Stanford, California 94305, United States Correspondence to:
[email protected] For submission as an original research article to Environmental Science & Technology
Abstract: 223 Main Text: 5833 Figures: 900 Tables: 600 Total: 7556
ACS Paragon Plus Environment
1
Page 3 of 35
Environmental Science & Technology
45
Abstract
46
Fecal indicator organisms are measured to indicate the presence of fecal pollution, yet the
47
association between indicators and pathogens varies by context. The goal of this study was to
48
empirically evaluate the relationships between indicator Escherichia coli, microbial source
49
tracking markers, select enteric pathogen genes and potential sources of enteric pathogens in
50
600 rural Bangladeshi households. We measured indicators and pathogen genes in stored
51
drinking water, soil and on mother and child hands. Additionally, survey and observational data
52
on sanitation and domestic hygiene practices were collected. Log10 concentrations of indicator E.
53
coli were positively associated with the prevalence of pathogenic E. coli genes in all sample types.
54
Given the current need to rely on indicators to assess fecal contamination in the field, it is
55
significant that in this study context indicator E. coli concentrations, measured by IDEXX Colilert-
56
18, provided quantitative information on the presence of pathogenic E. coli in different sample
57
types. There were no significant associations between the human fecal marker (HumM2) and
58
human-specific pathogens in any environmental sample type. There was an increase in the
59
prevalence of Giardia lamblia genes, any E. coli virulence gene, and the specific E. coli virulence
60
genes stx1/2 with every log10 increase in the concentration of the animal fecal marker (BacCow)
61
on mothers’ hands. Thus, domestic animals were important contributors to enteric pathogens in
62
these households.
63 64 65 66
ACS Paragon Plus Environment
2
Environmental Science & Technology
Page 4 of 35
67
Introduction
68
In low- and middle-income countries diarrheal illness is a leading cause of morbidity and
69
mortality.1 In Bangladesh, 6% of the 129,000 deaths in children under five in 2013 were attributed
70
to diarrheal diseases.2 Diarrheal illness results from exposure to fecal pathogens which can be
71
transmitted from feces to a new human host through a variety of environmental pathways,
72
including fingers, fields, flies, fluids, and food, described in the F-diagram.3 Recent additions to
73
the F-diagram stress the importance of animal hosts by expanding fecal sources to include feces
74
from livestock, free-roaming animals and synanthropic rodents.4 While it is well known that
75
enteric pathogens are transmitted through these pathways, few studies have measured
76
pathogens to characterize exposure from different animal reservoirs.
77 78
There is a high potential for zoonotic enteric disease transmission in low- and middle- income
79
countries where animal husbandry is a primary source of income.5 For example, in Bangladesh
80
raising livestock such as cows, goats and chickens results in animals roaming freely within the
81
home environment.6 Close proximity to domesticated animals can lead to human exposure to
82
livestock feces. Many pathogens can be transmitted from animal feces to human hosts and result
83
in diarrheal illnesses.5 Of the four most common etiological agents of moderate to severe
84
diarrhea in children 0-11 months in Bangladesh, two (Cryptosporidium and Campylobacter) are
85
known to have important animal reservoirs.7,8
86 87
Measuring enteric pathogens in the environment can help identify reservoirs and potential
88
exposure pathways, which can better inform design of strategies to reduce human exposure.
ACS Paragon Plus Environment
3
Page 5 of 35
Environmental Science & Technology
89
However, there are many different fecal pathogens from humans and animals capable of causing
90
disease, making it infeasible to measure them all, especially given that most are difficult to
91
measure in the environment due to their low concentrations in environmental matrices and
92
costly and complex methods of detection.9 Therefore, fecal indicators are commonly used to
93
indicate the presence of fecal contamination, which may contain pathogens.9 Thermotolerant
94
coliforms and Escherichia coli are recommended indicators of drinking water quality by the World
95
Health Organization.10 Other fecal indicators such as Bacteroidales are used for identifying
96
sources of fecal contamination through microbial source tracking (MST).11 However, the degree
97
to which indicators and pathogens co-occur in environmental samples varies depending on their
98
concentrations in the original fecal source (pathogen concentrations depend on the infection
99
status of the human or animal whereas indicator organism concentrations are expected to
100
remain more stable) and the relative transport and die-off/growth rates of organisms once they
101
are in the environment. Additionally, correlations between indicator organisms and pathogens
102
in natural environments, such as soil and surface waters, are impacted by the concentrations of
103
naturally occurring E. coli in the environment, often referred to as “naturalized” E. coli. Thus, the
104
relationship of indicator organisms and specific pathogens in or on different environmental
105
matrices (e.g., water, soil, hands, fomites, food) varies spatially and temporally depending on
106
these factors.12 Nonetheless, a better empirical understanding of the relationship between
107
indicators and pathogens in specific environmental reservoirs and contexts in low- and middle-
108
income countries may improve the ability to estimate human health risk and identify fecal
109
sources and exposure pathways of greatest concern. Furthermore, pinpointing dominant sources
ACS Paragon Plus Environment
4
Environmental Science & Technology
Page 6 of 35
110
and reservoirs will help improve the design of targeted interventions that can reduce exposure
111
to fecal pathogens, mitigating the burden of diarrheal illnesses.13,14
112 113
The aim of this study was to measure fecal indicators and select human pathogen genes in
114
different reservoirs in the domestic environment to: 1) determine the association between
115
concentration of indicator E. coli and presence of pathogen genes; and 2) use microbial source
116
tracking and observations of animal feces in compounds to investigate potential human and
117
animal sources of detected pathogen genes.
118
Methods
119
We selected multiple types of microorganisms that are leading causes of diarrheal illness in
120
developing countries and have a range of host specificity between humans and animals. The
121
selected targets were indicators (E. coli, human-associated Bacteroidales-like gene, and non-
122
human animal-associated Bacteroidales gene) and pathogen genes (enteropathogenic E. coli
123
(EPEC), enteroaggregative E. coli (EAEC), shiga-toxin producing E. coli (STEC), enteroinvasive E.
124
coli (EIEC), enterotoxigenic E. coli (ETEC), norovirus GII, Giardia lamblia, and Cryptosporidium
125
spp.).We measured these targets in different reservoirs in the domestic environment within rural
126
Bangladeshi compounds - stored water, soil, and mother and child hands (Figure 1).
127 128
We reviewed current literature to determine which of the selected pathogens had the potential
129
to originate from both human and animal reservoirs or only human reservoirs. Pathogens
130
associated with both human and animal feces include Cryptosporidium, Giardia lamblia, EPEC,
131
STEC and ETEC.15–17 Humans serve as a reservoir of typical EPEC, while animals such as pigs and
ACS Paragon Plus Environment
5
Page 7 of 35
Environmental Science & Technology
132
chickens can serve as reservoirs of atypical EPEC.17,18 Though typical and atypical EPEC can be
133
differentiated by the presence of the gene bfpA, this study only measured eaeA and will therefore
134
not distinguish between these strains.19 STEC is characterized by the presence of shiga toxin
135
genes (stx1 or stx2) and is known to come from ruminants.17 Enterohemorrhagic E. coli (EHEC) is
136
a common type of STEC that can also carry the eaeA gene found in EPEC.17,20 ETEC is found in
137
both humans and animals but species-specific adhesion factors confer host specificity.8 We were
138
unable to differentiate ETEC with human adhesion factors in this study because the enterotoxins
139
we measured, st1b and lt1, can be found in both humans and animals.21–24 Previous studies have
140
isolated Giardia lamblia and Cryptosporidium spp. in many fecal sources and identified closely
141
related isolates from humans and animals.15,25 On the other hand, norovirus GII, EAEC and EIEC
142
are associated with primarily human sources, largely based on genetic analysis of pathogens in
143
human and animal fecal samples.17,26,27 While EAEC is also found in animals, the gene we
144
detected, aggR, has been identified in EAEC isolated from humans but not animals.27,28 EIEC has
145
also been found in primates, which are unlikely to be important contributors to fecal pollution in
146
the study communities.17 It should be noted that the associations between pathogens and
147
possible hosts in Figure 1 are based on current scientific understanding and evidence. Zoonotic
148
pathogen transmission is difficult to demonstrate and we are limited by the methods that have
149
been employed and the pathogens, virulence genes, locations, and animal hosts that have been
150
investigated.
151
Study Setting and Design
152
600 households were randomly selected from those enrolled in the sanitation and control arms
153
(300 households per arm) of the WASH Benefits randomized controlled trial in rural Bangladesh,
ACS Paragon Plus Environment
6
Environmental Science & Technology
Page 8 of 35
154
described elsewhere.29 One study household was selected per compound and compounds were
155
usually relatives living in adjacent households surrounding a central courtyard. Compounds were
156
clustered in groups of eight to reduce spillover effects and for the ease of community promoters
157
delivering the interventions. Clusters were geographically pair-matched and then randomly
158
allocated to different intervention arms.30 Study households were from the rural Gazipur,
159
Mymensingh, Tangail and Kishoreganj districts. WASH Benefits enrolled and followed households
160
with a pregnant woman. All study households contained a young child (aged 9-44 months) at the
161
time of our visit. In this same subset of households, another analysis in preparation was
162
conducted to investigate the impact of the sanitation intervention on the prevalence of pathogen
163
genes and microbial source tracker markers.31 Another related study investigated how the
164
presence of domestic animals impacted concentrations of indicator E. coli in reservoirs in a subset
165
of households from the control arm.32
166 167
Data Collection
168
During household visits between March and October 2015, trained field staff from the
169
International Centre for Diarrhoeal Disease Research, Bangladesh (icddr,b) collected
170
environmental samples, interviewed the female caregiver of young children regarding household
171
practices on handling of animal feces, and observed the number and type of animal fecal piles in
172
compound courtyards. Field surveyors also recorded the number and type of animals owned by
173
the household as reported by the caregiver.
174
Sample Collection
ACS Paragon Plus Environment
7
Page 9 of 35
Environmental Science & Technology
175
We collected approximately 720 stored drinking water samples, 720 soil samples, 720 mother
176
hand rinses, and 360 child hand rinses from 600 study households. 176 of households in the
177
control and 156 of households in the sanitation arms were sampled twice, approximately four
178
months apart to assess the impact of season. Hand rinse samples were collected by mothers
179
placing their left hand into a sterile Whirlpak bag (Nasco, Modesto, CA) filled with 250 mL of
180
distilled water. The hand was massaged from the outside of the bag for 15 seconds, followed by
181
15 seconds of shaking. The same procedure was repeated with the right hand in the same bag.
182
To collect child hand rinse samples, respondents placed their child’s hand into a separate
183
Whirlpak bag and followed the same procedure. Soil samples were collected from a 30 x 30 cm2
184
area as close to the house entrance as possible by scraping the top layer of soil within a stencil
185
into a sterile Whirlpak bag using a sterile disposable plastic scoop. The sample area was scraped
186
both vertically and horizontally. Stored water samples were collected by asking mothers to
187
provide a glass of water as they would give it to their child under five. The provided water was
188
poured into a sterile Whirlpak bag. All samples were transported to the icddr,b field laboratory
189
on ice and processed within 12 hours of collection.
190
Indicator E. coli Enumeration and Detection of Pathogenic Genes
191
E. coli were enumerated (100 mL sample volumes) using IDEXX Colilert-18 (IDEXX Laboratories,
192
Westbrook, Maine). Trays were incubated at 44.5 °C for 18 hours.33 To prepare for IDEXX analysis,
193
hand rinse samples were diluted 1:2 with sterile distilled water in Whirpak bags while stored
194
water was analyzed undiluted. 20 g of soil was homogenized with 200 mL of distilled water and
195
diluted 1:104. Soil was dried at 110° C for 24 hours to determine moisture content. Each lab
196
technician processed one blank (sterile distilled water) for E. coli enumeration per day.
ACS Paragon Plus Environment
8
Environmental Science & Technology
Page 10 of 35
197 198
Pathogenic E. coli genes were identified using previously published methods.20,34 All samples
199
positive for E. coli as determined by IDEXX Colilert-18 were archived for subsequent pathogenic
200
E. coli gene analysis. Broth from up to 20 positive large wells was aseptically extracted from IDEXX
201
trays, composited and centrifuged. Pellets were treated with 10x the pellet volume (≈ 0.1 mL) of
202
RNAlater (Qiagen, Germantown, MD), stored at -80 °C and transported to UC Berkeley at room
203
temperature. Each lab technician analyzed one lab blank per week for the archiving process by
204
archiving wells from IDEXX trays incubated with sterile distilled water. Pellets were stored at -80
205
°C upon arrival at UC Berkeley. DNA was extracted from bacteria pellets using the DNeasy Blood
206
and Tissue Kit (Qiagen, Germantown, MD). Multiplex reactions were used to detect seven E. coli
207
virulence genes indicative of five possible pathotypes of E. coli: EAEC (aggR), EPEC or EHEC
208
(eaeA), STEC (stx1/ stx2), EIEC (ipaH), ETEC (lt1/st1b) (Table S1). PCR products were run on 2%
209
agarose gels at 110 V for 30 min. Additional information on DNA extraction and PCR protocols is
210
available in the supporting information.
211
Filtration and Nucleic Acid Extraction for qPCR Targets
212
Laboratory technicians pre-processed samples in the field laboratory by filtering 50 mL of hand
213
rinse samples and up to 500 mL of stored water (range: 100-500 mL) through a 0.45 m HA filter
214
(Millipore, Burlington, MA) to capture bacteria and viruses. 0.5 mL of 2.5M MgCl2 was added to
215
every 50 mL of sample to increase retention of viruses.35 Filters were treated with 0.5 mL of
216
RNAlater, stored at -80 °C and transported to UC Berkeley at room temperature. Once per day,
217
each lab technician processed a lab blank by filtering 5 mL of sterile distilled water with 0.5 mL
218
of MgCl2, followed by the addition of RNAlater. Five grams of soil were weighed and stored at -
ACS Paragon Plus Environment
9
Page 11 of 35
Environmental Science & Technology
219
80 °C until transport at room temperature to UC Berkeley in accordance with a USDA soil import
220
permit (PPQ 525). Filter and soil samples were stored at -80 °C at UC Berkeley. DNA and RNA
221
were extracted from filters and 0.25 grams of soil samples using modified Mobio PowerWater
222
and PowerViral (both now Qiagen, Germantown, MD) protocols, described in previous studies
223
and in the supporting information.36,37 An extraction blank was included in each batch of 23
224
samples. Extraction efficiencies for DNA and RNA were determined in a subset of samples using
225
Pseudomonas syringae pv. phaseolicola (pph6) and MS2 (Table S2).
226
qPCR Assays
227
Samples were analyzed for norovirus GII, Giardia lamblia, Cryptosporidium spp. genes, and
228
microbial source tracking markers using quantitative PCR (Table S3). Specifically, to detect human
229
fecal contamination we used the HumM2 assay which targets a gene for a hypothetical protein
230
in human-associated Bacteroidales-like microorganisms.38,39 We obtained a research license
231
(#864-15) from the U.S. Environmental Protection Agency to use the patented HumM2 assay. To
232
detect non-human, animal-associated fecal contamination we used the BacCow assay which,
233
although designed to target the 16S rRNA of Bacteroidales associated with ruminants, has been
234
shown to detect Bacteroidales in cows, ducks, goats, and chickens but not humans in rural
235
Bangladesh.36,40 The microbial source tracking markers were previously evaluated for sensitivity
236
and specificity in rural Bangladesh.36 BacCow was sensitive but not specific to ruminant feces and
237
HumM2 performed the best out of all tested human-associated assays (HumM2, HF183 and
238
BacHum), although it also amplified in the feces of chickens and goats.36
239
ACS Paragon Plus Environment
10
Environmental Science & Technology
Page 12 of 35
240
Approximately 100 stored water, soil and mother and child hand rinse samples were processed
241
for norovirus GII, Giardia and Cryptosporidium genes. Among this subset, 0) compared to those with no fecal piles. For
300
chicken/non-chicken poultry feces the prevalence ratio represents the prevalence of pathogen
301
genes in households with > 5 fecal piles compared to those with ≤ 5 fecal piles. Cow patties for
302
cooking were differentiated from cow feces because they were formed into cakes and dried in
303
the sun. Exposures (feces of different animals) were screened against outcomes (pathogen
304
genes) in bivariate models and those with a p-value