Correlating Transport Behavior with Cell Properties ... - ACS Publications

Jun 11, 2010 - environment affect variability of cell properties and transport behavior of eight porcine E. ... in bacterial deposition rates to Pyrax...
2 downloads 0 Views 520KB Size
Environ. Sci. Technol. 2010, 44, 5008–5014

Correlating Transport Behavior with Cell Properties for Eight Porcine Escherichia coli Isolates C A R L H . B O L S T E R , * ,† KIMBERLY L. COOK,† IAN M. MARCUS,‡ BERAT Z. HAZNEDAROGLU,‡ AND SHARON L. WALKER‡ U.S. Department of AgriculturesAgricultural Research Service, 230 Bennett Ln., Bowling Green, Kentucky 42104, and Department of Chemical & Environmental Engineering, University of California, Riverside, Riverside, California 92521

Received March 31, 2010. Revised manuscript received May 25, 2010. Accepted May 28, 2010.

In this study we investigate how growth stage and depositional environment affect variability of cell properties and transport behavior of eight porcine E. coli isolates. We compared the surface properties for cells harvested during exponential and stationary growth phase and their transport behavior through columns packed with either uncoated or Fe-coated quartz sand. We then investigated correlations between measured cell properties and fitted bacterial attachment efficiencies. For both growth stages we found that bacterial attachment efficiencies in the uncoated quartz sand varied among the eight different isolates by over an order of magnitude whereas attachment efficiencies in the Fe-coated sands varied by a factor of less than two. With the exception of one isolate, growth condition had minimal impact on attachment efficiencies to the uncoated sands. A strong and statistically significant inverse relationship was observed between bacterial attachment efficiencies in the uncoated quartz sand columns and log-transformed zeta potential, whereas a mild yet statistically significant relationship between bacterial attachment efficiencies in the Fe-coated sandsandcellwidthwasobserved.Fortheexperimentalconditions used in our study, we found that variability in E. coli transport was more dependent on the depositional environment than on growth conditions.

Introduction Infiltration of fecal material into the subsurface can result in the contamination of groundwater supplies by pathogenic microorganisms such as bacteria, viruses, and protozoa, thereby posing a threat to public health. Studies have shown that fecal contamination is routinely detected in drinking water wells and groundwater is often reported as the source of waterborne disease outbreaks (1). To assess whether a groundwater source is at risk for fecal contamination, agencies responsible for monitoring water supplies generally test for the presence of nonpathogenic microorganisms commonly found in fecal materialsreferred to as indicator organismssdue to the difficulties and limitations of testing for individual pathogens. Indicator organisms are also used * Corresponding author phone: (270) 781-2632 ext. 244; fax: (270) 781-7994; e-mail: [email protected]. † U.S. Department of Agriculture. ‡ University of California, Riverside. 5008

9

ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 44, NO. 13, 2010

in watershed fate and transport models to predict pathogen movement through the landscape (2). While the use of indicator organisms for assessing the microbiological quality of drinking water has notable limitations, their use will likely continue in the foreseeable future (3). Therefore, a good understanding of the factors controlling their fate and transport in the environment is needed. One of the most commonly used indicator organisms in groundwater systems is E. coli. (Though some strains of E. coli are pathogenic, the majority are not) Under the new Ground Water Rule, E. coli is one of only three indicators EPA has established for determining whether a groundwater source used as a public water supply is fecally contaminated (1). Given the importance of E. coli as an indicator of fecal contamination in groundwater systems, there have been a number of recent studies focusing on understanding the various biological, geochemical, and physical factors affecting E. coli transport through porous media (4). These factors include cell acclimation time (5, 6), growth conditions (7, 8), composition of the extracellular polymeric substances (EPS) (5, 9), surface charge and hydrophobicity (9-11), cell motility (11, 12), temperature (6), ionic strength (13), sediment size (14), hydraulic conductivity (15), presence of metal-oxyhydroxide coated sand grains (13), and soil moisture content (16). Of particular relevance concerning E. coli transport through porous media are recent studies that have shown that a large diversity exists in both cell properties and transport behavior among different E. coli isolates. For instance, Morrow et al. (17) reported a nearly 2-fold difference in bacterial deposition rates to Pyrax for three E. coli isolates. Bolster et al. (9) observed deposition rates to vary by a factor of 10 among 12 E. coli isolates obtained from six different sources. Lutterodt et al. (11) also observed bacterial attachment efficiencies to vary by a factor of 10 among six E. coli isolates obtained from a soil used for cattle grazing. Similarly, Foppen et al. (10) measured breakthrough of 54 different E. coli isolates from various sources and observed bacterial attachment efficiencies to range by nearly a factor of 10 resulting in a two-log variation in effluent bacteria concentration following passage through 7 cm columns packed with quartz sand. The cause of this variability in E. coli transport, however, remains unclear and warrants further research. Our objective in this study was to further our understanding of the factors controlling variability in E. coli transport by conducting cell characterization and transport experiments using eight different E. coli isolates obtained from a swine lagoon. We compared the cell properties and transport behavior for cells harvested at two different growth stages (exponential and stationary) and two depositional environments (clean and Fe-coated quartz sand). We then investigated correlations between measured cell properties and fitted bacterial attachment efficiencies.

Materials and Methods E. coli Isolation and Isolate Selection. Five 1 L samples of liquid swine effluent were collected from a swine lagoon located on the farm at Western Kentucky University. The samples were collected at five randomly chosen locations within the lagoon at a depth of ∼1 m using 1 L sterilized Nalgene bottles attached to an adjustable sampling pole. The samples were returned to the laboratory and processed immediately. Ten mL from each 1 L bottle were serially diluted to 10-6 in phosphate buffer. One hundred µL from each dilution was plated onto m-Tech plates and incubated at 44.5 °C for 24 h. Following incubation, each presumptive E. 10.1021/es1010253

 2010 American Chemical Society

Published on Web 06/11/2010

coli colony (i.e., purple colony) was plated on EMB and mFC agar. This yielded a total of 339 presumptive E. coli colonies. E. coli isolates were fingerprinted by BOX PCR analysis using the BOX AIR primer as previously described (18). Eight isolates were chosen from the major fingerprint groups to be used in the column studies. Each of the eight isolates was further confirmed to be E. coli by chemical analysis using Enterotube following the manufacturer’s directions and by PCR analysis of the uidA gene as previously described (19). All isolates were stored at -80 °C with 15% glycerol. Isolate Typing and Evaluation of fimH and agn43. The eight isolates were placed in one of four phylogenetic groups (A, B1, B2, or D) using the triplex PCR method of Clermont et al. (20) as modified by Higgins et al. (21). The PCR-based lipopolysaccharide core typing method of Amor et al. (22) was used to identify the outer core oligosaccharide type of each E. coli isolate. PCR analyses for the presence of agn43 and fimH gene clusters were carried out as previously described (23) except the initial denaturation step was extended to 15 min. Assays were carried out in Qiagen HotStart Taq Master Mix (Qiagen, Valencia, CA) in a total volume of 25 µL with 2 µL of genomic DNA extract and 0.5 µM primer in a PTC-200 model thermal cycler (MJ Research/ BioRad, Hercules, CA). PCR mixtures (10 µL) were electrophoresed on 1% agarose gel 1 h at 75 V. Ethidium bromide stained gel images were captured with a FOTO/Analyst Investigator/Eclipse system (Fotodyne Inc., Hartland, WI). Primer sequences, product size, and targeted genes are shown in Supporting Information (SI) Table S1. Cell Preparation. Prior to transport or surface characterization experiments, 40 µL of an overnight culture of E. coli was inoculated into 40 mL of Luria-Bertani (LB) broth and incubated at 37 °C until reaching either midexponential growth phase (3.5 h) or stationary phase (16 h). A refrigerated centrifuge equipped with fixed angle rotor was used to pellet the cells with an applied 3700g force for 15 min at 4 °C. The cell pellet was resuspended in 1 mM KCl solution prepared with deionized (DI) water and reagent-grade KCl (Fisher Scientific) with no pH adjustment (pH 5.6-5.8). This process was repeated twice in order to ensure complete removal of the growth medium. Cell Characterization. The electrophoretic mobility, hydrophobicity, surface charge density, extracellular polymeric substance (EPS) composition, and cell size were measured for each isolate using previously published methods (9). Greater details of the cell characterization methods can be found in the SI. Transport Experiments. Transport experiments were conducted through water-saturated columns packed with either clean quartz sand (Unimin, New Caanan, CT) or Fecoated quartz sand ranging in size from 710 to 850 µm in diameter. This sand size was chosen to minimize the impact of physical straining on the removal of bacteria from the aqueous phase. Sand was cleaned and portions of the acidwashed sand were coated with Fe-oxyhydroxides using previously published methods (13). Streaming potential of the clean and Fe-coated sand was measured using an electro kinetic analyzer (Anton Paar GmbH, Graz, Austria) following the procedure described in Kim et al. (24). Columns were packed by slowly pouring water-saturated autoclaved sand into 2.5 cm diameter Chromaflex Chromatography Columns (Kontes Glass Co., Vineland, NJ) filled with a 1 mM KCl electrolyte solution while vibrating the columns. Porosity ranged from 0.36 to 0.39 as determined gravimetrically. Column lengths ranged from 9.9 to 10.1 cm. After packing was completed, columns were operated in an upward direction using syringe pumps (model 200 syringe pump, KD Scientific Inc., New Hope, PA), and approximately 10 pore volumes of the electrolyte solution were passed through each column to equilibrate the sand pack. Columns

were operated at a flow rate of 0.67 mL/min for a Darcian velocity of ∼0.36 cm/min. Bacterial solutions with an optical density of ∼0.20 (∼1 × 108 cells/mL) were injected into each column for 38 min (∼1.2 pore volumes) followed by 2 pore volumes of bacteria-free electrolyte solution. Effluent was collected every 3 min using a Spectra/Chrom CF-1 fraction collector (Spectrum Chromatography, Houston, TX). Effluent concentrations of E. coli were determined by measuring the optical density of the samples at a wavelength of 546 nm. All experiments were duplicated on a separate day to ensure true replication. In several columns bromide was added as a conservative tracer so that a dispersion coefficient for the columns could be determined. Data Analysis. The breakthrough curves (BTCs) for each treatment were fit with both first-order and second-order deposition models using weighted least-squares regression. Weights were calculated assuming a coefficient of variation in the observed BTC concentrations of 5%; a value equivalent to the median coefficient of variation observed between duplicate data points for all our column experiments combined. Bacterial attachment efficiencies (R), which describe a cell’s ability to become associated with the sand surface, were calculated using the hemispheres-in-cell model of Ma et al. (25) from the fitted deposition parameters obtained from the best-fit model for each treatment. (See SI for further details on the transport equations and modeling methodology used.) Correlation and regression analyses were used to determine whether fitted R values were significantly correlated with any of the measured cell properties. To account for the effects of multiple parameters, multiple regression analysis was performed. All regression analyses were performed using PROC REG in SAS 9.2 (SAS Institute, 2003).

Results and Discussion Isolate Selection and Characterization. Using a threshold of 80% similarity in BOX fingerprint patterns, 25 unique profiles were differentiated from the 339 E. coli isolates. Eight isolates, representing 81% of the diversity of isolates, were selected for transport studies and cell surface characterization. Each of these eight isolates was confirmed to be E. coli by PCR amplification of the uidA gene and by the Enterotube II method (SI Tables S4 and S5). Based on the genotyping triplex assay of Clermont et al. (20), six of the isolates grouped with the B1 genotype and two with group A (SI Table S5). Clermont et al. (20) also found that isolates from swine feces were dominated by the B1 genotype with remaining isolates being genotype A. The outer core oligosaccharide types of the eight E. coli isolates were identified as R1, R3, K12, and R1/R3 mixed type (SI Table S5). The identification of isolate SP2C04 as core type K12 was unexpected since this group is infrequent among natural populations. R1 core type dominates in natural populations and three of our isolates fit this core type. Isolate SP3B03 was positive for both R1 and R3 core groups. Appelmelk et al. (26) described an R1/R3 mixed type which reacted equally well with polyclonal R-R1 and with R-R3 antiserum. They suggest that the R3 reaction might be due to cross-reactions with non-LPS antigens. Cell Properties. Measured cell properties of the eight porcine E. coli isolates analyzed in this study are shown in Table 1. In general, we observed a large range in cell properties among the different E. coli isolates. Zeta potential varied by a factor of nearly 25 for exponential-phase cells with values ranging from -2.0 to -46 mV and varied by a factor of nearly 100 for stationary-phase cells with values ranging from -0.8 to -70.5 mV. Hydrophobicity also varied by over an order of magnitude ranging from 6 to 81% for the exponentialphase cells and from 5 to 81% for the stationary-phase cells. Surface charge density varied by a factor of 6 and 11 for VOL. 44, NO. 13, 2010 / ENVIRONMENTAL SCIENCE & TECHNOLOGY

9

5009

TABLE 1. Measured Values of Percent Recovery and Selected Cell Properties (Standard Deviations Are in Parentheses) for Exponential- (exp) and Stationary-Phase (stat) Cells hydrophobicity (%) isolate SP1C10 SP1H01 SP2A12 SP2B07 SP2C04 SP3B03 SP4C08 SP4H03

exp 26 (9.3) 35 (2.3) 13 (1.5) 11 (1.8) 7.0 (0.3) 11 (5.5) 81 (2.1) 6.0 (0.2)

zeta potential (mV)

surface charge density (µC/cm2)

stat

exp

stat

exp

15 (11.6) 46 (6.5) 33 (8.3) 5.0 (2.9) 81 (1.2) 36 (10.) 45 (1.5) 19 (3.8)

-4.5 (0.04) -6.2 (1.8) -26 (4.5) -45 (0.76) -25 (1.6) -0.89 (2.19) -7.4 (0.49) -46 (0.78)

-4.1 (0.16) -4.4 (0.95) -31 (0.09) -38 (4.8) -26 (1.0) -0.80 (0.51) -6.5 (1.3) -71 (2.3)

1.0 × 10

9

4

1.3 × 10

3

1.2 × 104

4.4 × 102

1.7 × 103

1.5 × 103

2.4 × 103

1.0 × 103

1.1 × 104

7.3 × 102

1.8 × 103

4.1 × 102

1.7 × 103

1.6 × 102

8.6 × 103

7.6 × 102

exponential- and stationary-phase cells, respectively, whereas cell size was relatively uniform among the different isolates with cell length only ranging from 1.5 to 1.8 µm and cell width ranging from 0.76 to 0.97 µm. The sugar and protein contents of the EPS varied by factors ranging from 4 (sugar content of stationary-phase cells) to 9 (sugar content of exponential-phase cells) (SI Table S6). Such large diversity in cell properties among different E. coli isolates is consistent with what has been reported in the literature under a variety of experimental conditions (5, 8-11, 17). Increasing incubation time in LB broth from 3.5 to 16 h had varying effects on measured cell properties (Table 1). For example, increasing incubation time increased hydrophobicity for five of the isolates while decreasing it for three, with the largest change observed for SP2C04 which went from hydrophilic (7%) to very hydrophobic (81%). Similar to earlier work (8, 27, 28), changes in zeta potential were generally minimal with the exception of SP4H03 in which increasing incubation time resulted in zeta potential becoming significantly more negative (-46 to -71 mV). For seven of the isolates surface charge density values were lower, in most cases by a factor of 10 or greater, for the stationaryphase cells compared with the exponential-phase cells. This is likely due to the fact that the bacteria enter a starvation state by 16 h, and therefore are conserving their resources by producing less proteins, sugars, and lipids resulting in less functional groups to titrate on the surface (7, 29, 30). Cell length remained relatively constant with growth stage whereas cell width notably increased for seven of the isolates. Total EPS increased in stationary-phase cells for six of the isolates (SI Table S6). Transport Experiments. Similar to cell properties, transport behavior through the uncoated quartz sand varied among the different E. coli isolates (Figure 1A and B). The percent recovery of cells eluted from the columns ranged from 48% to nearly 98% for both the exponential- and stationary-phase cells. Five of the eight isolates had recoveries of 89% or greater whereas two of the eight isolates had recoveries of 60% or less with the lowest recovery (∼48%) observed for SP3B03. Bacterial attachment efficiencies (R) varied among the eight isolates by over an order of magnitude for both growth conditions; for the exponential-phase cells R ranged from 0.032 to 0.47 whereas for the stationary-phase cells R ranged from 0.024 to 0.55 (Figure 1C). Attachment efficiencies measured under the two growth conditions were similar for all isolates with the exception of SP4H03, where R decreased from 0.088 to 0.024 when increasing incubation time from 3.5 to 16 h. This is likely due 5010

stat

ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 44, NO. 13, 2010

cell length (µm)

cell width (µm)

exp

stat

exp

stat

1.53 (0.033) 1.61 (0.12) 1.61 (0.092) 1.54 (0.055) 1.73 (0.008) 1.74 (0.024) 1.59 (0.098) 1.66 (0.024)

1.60 (0.035) 1.77 (0.071) 1.65 (0.15) 1.50 (0.078) 1.62 (0.089) 1.73 (0.193) 1.74 (0.083) 1.60 (0.022)

0.86 (0.023) 0.81 (0.037) 0.76 (0.013) 0.76 (0.030) 0.88 (0.010) 0.83 (0.029) 0.86 (0.021) 0.92 (0.023)

0.88 (0.038) 0.96 (0.019) 0.91 (0.054) 0.87 (0.016) 0.91 (0.037) 0.94 (0.069) 0.97 (0.023) 0.89 (0.003)

to the fact that zeta potential for this isolate became significantly more negative with increasing incubation time (Table 1), thereby increasing the electrostatic repulsive barrier to deposition (see discussion below). Walker et al. (28) observed stationary-phase cells of E. coli K12 to be more adhesive to quartz sand than exponential-phase cells. The authors attributed this behavior to the decrease in surface charge density they observed for the stationary-phase cells. Similar to Walker et al. (28) we also observed a large decrease in surface charge density (Table 1) with increasing incubation time yet we did not find any differences in transport behavior through the uncoated quartz sand for seven of the eight isolates tested. For the one isolate that we did see an effect, the effect was opposite to what Walker et al. (28) observed thus highlighting the importance of not extrapolating transport behavior observed for one strain of E. coli to all strains of E. coli. Because studies have shown that both surface properties and transport behavior of E. coli are dependent on growth conditions such as growth media and starvation status of the cells (5, 12), additional studies using a wider range of environmentally relevant growth conditions are needed to fully determine the importance of growth stage on E. coli transport behavior. The presence of the Fe-oxyhydroxide coatings on the quartz sand grains resulted in decreased effluent concentrations for seven of the eight isolates (Figure 2A and B), results consistent with previously published studies (13, 31-33). The one isolate, SP3B03, which did not experience decreased transport through the Fe-coated sand had a near-neutral charge (Table 1). The range in transport behavior among the different E. coli isolates through the Fe-coated sands was much less than what was observed for the uncoated sands with percent recoveries ranging from ∼35% for SP2C04 to ∼54% for SP3B03. Bacterial attachment efficiencies ranged by a factor of less than two for the Fe-coated sands with values ranging from 0.45 to 0.82 for the exponential-phase cells and from 0.52 to 0.87 for the stationary-phase cells (Figure 2C). For all isolates the R values for the stationaryphase cells were higher than those observed for the exponential-phase cells, though these differences were generally small and within the uncertainties of R. The presence of metal-oxyhydroxide coatings on quartz sand is often assumed to result in a net positive charge on the sand surface thereby removing any electrostatic barrier to bacterial deposition. Under such favorable conditions, all of the cells that strike a collector are expected to attach to the collector and therefore variability in transport is expected to be reduced compared with unfavorable conditions. This

FIGURE 1. Observed (symbols) and fitted (solid line) breakthrough curves for eight porcine E. coli isolates through uncoated quartz sand columns for cells harvested in the exponential (A) and stationary (B) phases. Also included are the fitted values and 95% confidence intervals (C) of the bacterial attachment efficiencies (r) for exponential- (closed symbols) and stationary-phase (open symbols) cells. Because we assumed that the only source of uncertainty in r was the uncertainty in k obtained from the weighted least-squares fit, the uncertainties associated with our fitted r values may underestimate their true uncertainties.

FIGURE 2. Observed (symbols) and fitted (solid line) breakthrough data for eight porcine E. coli isolates through Fe-coated quartz sand columns for cells harvested in the exponential (A) and stationary (B) phases. Also included are the fitted values and 95% confidence intervals (C) of the bacterial attachment efficiencies (r) for exponential- (closed symbols) and stationary-phase (open symbols) cells. Because we assumed that the only source of uncertainty in r was the uncertainty in k obtained from the weighted least-squares fit, the uncertainties associated with our fitted r values may underestimate their true uncertainties.

is indeed what we observed with our E. coli isolates and what has been reported for nonbiological colloids (34). However, all of our fitted R values were significantly lower than the expected value of one. In measuring the zeta potential of our Fe-coated sands we found that the surface charge was negative rather than positive (-6.1 mV), and DLVO calculations show that a barrier to deposition existed in the Fecoated sands for seven of the isolates (SI Table S7), likely

explaining, at least in part, why our observed values of R were less than one. The presence of the Fe-coatings did, however, reduce the energy barrier by 50% or more compared to the uncoated sand thereby resulting in less unfavorable conditions and increased deposition (SI Table S7). The one exception was with isolate SP3B03 in which no energy barrier existed for either sand type and explains why the R values were similar for the uncoated and Fe-coated sands for this one organism. The fact that R values for this organism were VOL. 44, NO. 13, 2010 / ENVIRONMENTAL SCIENCE & TECHNOLOGY

9

5011

FIGURE 3. Relationship between zeta potential and r for uncoated sands (A) and cell width and r for Fe-coated sands (B). Included are the best-fit linear regressions. below one, even though no energy barrier existed, indicates that additional factors such as steric repulsion, surface charge heterogeneities, or surface structures were limiting bacterial deposition. Alternatively, our lower-than-expected values of R may be a result of limitations in using colloid-filtration theory for rod-shaped bacteria. Correlating r with Cell Properties. We observed a strong and statistically significant inverse correlation (r ) -0.89; P < 0.001) between R and the log-transform of the absolute value of zeta potential in the uncoated sands (Figure 3A; SI Table S8). No other parameter was found to be significantly (P < 0.05) correlated with R for the uncoated sands (SI Table S8). Multiple linear regression analysis indicated that the best overall model (r2 ) 0.85; P < 0.001) for describing the variability in R was a 2-parameter model that included sugar content and log-transformed zeta potential; however, zeta potential was clearly the most important factor with a partial r2 of 0.78 compared to 0.07 for sugar content. The observation that zeta potential is the dominant factor controlling bacterial transport through columns packed with like-charged collectors is in qualitative agreement with Derjaguin-Landau-Verwey-Overbeek (DLVO) theory which predicts that bacterial attachment efficiencies to like-charged collectors will be higher for cells with a lower net negative charge due to a reduction in the electrostatic double-layer repulsion between the negatively charged bacteria and sand grains. This behavior is commonly observed for individual isolates when the surface charge has been modified by changes in ionic strength (7, 13, 35); however, when comparing multiple E. coli isolates, bacterial deposition rates have not been shown to be correlated with zeta potential (9-11). We cannot explain why this correlation was not 5012

9

ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 44, NO. 13, 2010

observed in these studies but note that two of the three studies compared E. coli isolates obtained from multiple sources (9, 10). It is possible that additional confounding factors may be present when correlating bacterial transport to surface properties for E. coli isolates obtained from multiple sources. For the Fe-coated sands a small yet statistically significant correlation (r ) 0.54, P ) 0.03) was observed between cell width and R (Figure 3B; SI Table S8). Using multiple linear regression analysis, no other parameters were found to add significant (P < 0.05) improvements for describing the observed variability in R. Because we used cell width in our calculation of η, this correlation suggests that cell width may play an additional role in bacterial deposition not currently incorporated into η. However, additional studies are needed to determine if this observed correlation is of both statistical and practical importance. In addition to measuring various cell surface properties, we also evaluated the presence of two genes that encode for phase-variable surface proteins (fimH which encodes for type 1 fimbriae and agn43 which encodes for Antigen 43) that have been associated with bacterial adhesion to surfaces (36, 37). Consistent with Yang et al. (23), we found that all isolates were positive for fimH suggesting that this gene may be ubiquitous in livestock fecal E. coli strains and therefore its occurrence is not a useful tool for predicting transport behavior (SI Table S5). Conversely, only three of the eight isolates tested in our study were positive for the agn43 gene (SI Table S5) and two of these isolates, SP3B03 and SP1H01, had high attachment efficiencies to the quartz sand grains. Although the presence of the agn43 gene does not necessarily mean that the gene is expressed and that the Ag43 outer membrane protein is present (36), our results are consistent with both Yang et al. (12) and Lutterodt et al. (11) who found that the expression of this gene decreased E. coli transport through porous media. Implications for Modeling E. coli Transport. Consistent with other studies, we observed a large range in the transport behavior of different E. coli isolates through porous media. Morrow et al. (17) attributed their observed differences among the three different isolates used in their study to differences in cell hydrophobicity. Bolster et al. (9) observed a mild inverse relationship between cell width and bacterial deposition rates for their 12 E. coli isolates. Both Lutterodt et al. (11) and Yang et al. (12) found that cell motility and Ag43 expression were the dominant factors controlling transport of various E. coli isolates through columns packed with a variety of materials. On the other hand, Foppen et al. (10) did not find any correlation between transport behavior and measured cell properties for the 54 E. coli isolates used in their study. In this study we found that cell retention in columns packed with clean quartz sand was inversely correlated with the zeta potential of the cells whereas a mild relationship between cell retention and cell width in Fe-coated sands was observed. These results indicate that the factors controlling the large diversity in transport behavior among different E. coli isolates likely vary depending on the experimental conditions and source of the E. coli isolates. This large diversity in cell properties and transport behavior among different E. coli isolates clearly demonstrates the inadequacy of using a single strain of E. coli to assess the general transport behavior of all E. coli strains. Furthermore, this variability shows the insufficiency of using a single deposition rate to model E. coli transport, especially for conditions unfavorable to deposition. As an example, we conducted an additional set of transport experiments with bacterial solutions containing both SP3B03 (the isolate with the highest R) and SP2C04 (the isolate with the lowest R) at three different concentration ratios (90:10, 50:50, and 10:90). The combined transport behavior of the cells was then predicted using a simple bimodal transport model where

FIGURE 4. Observed (symbols) and predicted (lines) BTCs of three different bacterial solution compositions. The solutions included 90% SP3B03 and 10% SP2C04 (90/10), 50% SP3B03 and 50% SP2C04 (50/50), and 10% SP3B03 and 90% SP2C04 (10/90). the attachment efficiency for both isolates was weighted by the fraction of that isolate in solution (38). Figure 4 shows that by properly accounting for the bimodal distribution in R, we were able to adequately predict the BTC behavior of the combined solutions through the uncoated quartz sand columns. Though this is an admittedly oversimplified example, it does illustrate the importance of including a distribution of bacterial deposition rates when modeling E. coli transport from a single fecal source. Continued research is required to assess the impacts of biological variability (both inter- and intrastrain variability) on modeling the transport of E. coli in the environment.

Acknowledgments We thank Stacy Antle for assistance with the transport studies. This research was part of USDA-ARS National Program 206: Manure and By-product Utilization. Mention of trade names or commercial products is solely for the description of experimental procedures and does not imply recommendation or endorsement by the USDA.

Supporting Information Available Additional information on cell characterization methods and Tables S1-S8. This material is available free of charge via the Internet at http://pubs.acs.org.

Literature Cited (1) U.S. Environmental Protection Agency. National Primary Drinking Water Regulations; Ground Water Rule; Final Rule; EPA: Washington, DC, 2006. (2) Pachepsky, Y. A.; Sadeghi, A. M.; Bradford, S. A.; Shelton, D. R.; Guber, A. K.; Dao, T. Transport and fate of manure-borne pathogens: Modeling perspective. Agric. Water Manage. 2006, 86, 81–92. (3) Indicators for Waterborne Pathogens; National Research Council: Washington, DC, 2004. (4) Foppen, J. W. A.; Schijven, J. F. Evaluation of data from the literature on the transport and survival of Escherichia coli and thermotolerant coliforms in aquifers under saturated conditions. Water Res. 2006, 40, 401–426. (5) Haznedaroglu, B. Z.; Bolster, C. H.; Walker, S. L. The role of starvation on Escherichia coli adhesion and transport in saturated porous media. Water Res. 2008, 42, 1547–1554. (6) Castro, F. D.; Tufenkji, N. Relevance of nontoxigenic strains as surrogates for Escherichia coli O157:H7 in groundwater contamination potential: Role of temperature and cell acclimation time. Environ. Sci. Technol. 2007, 41, 4332–4338. (7) Walker, S. L.; Hill, J. E.; Redman, J. A.; Elimelech, M. Influence of growth phase on adhesion kinetics of Escherichia coli D21g. Appl. Environ. Microbiol. 2005, 71, 3093–3099.

(8) Yang, H.-H.; Morrow, J. B.; Grasso, D.; Vinopal, R. T.; Smets, B. F. Intestinal versus external growth conditions change the surficial properties in a collection of environmental Escherichia coli isolates. Environ. Sci. Technol. 2006, 40, 6976–6982. (9) Bolster, C. H.; Haznedaroglu, B. Z.; Walker, S. L. Diversity in cell properties and transport behavior among 12 different environmental Escherichia coli isolates. J. Environ. Qual. 2009, 38, 465–472. (10) Foppen, J. W.; Lutterodt, G.; Roling, W. F. M.; Uhlenbrook, S. Towards understanding inter-strain attachment variations of Escherichia coli during transport in saturated quartz sand. Water Res. 2009. (11) Lutterodt, G.; Basnet, M.; Foppen, J. W. A.; Uhlenbrook, S. The effect of surface characteristics on the transport of multiple Escherichia coli isolates in large scale columns of quartz sand. Water Res. 2009, 43, 595–604. (12) Yang, H.-H.; Morrow, J. B.; Grasso, D.; Vinopal, R. T.; Dechesne, A.; Smets, B. F. Antecedent growth conditions alter retention of environmental Eschericia coli isolates in transiently wetted porous media. Environ. Sci. Technol. 2008, 42, 9310–9316. (13) Bolster, C. H.; Walker, S. L.; Cook, K. L. Comparison of Escherichia coli and Campylobacter jejuni transport in saturated porous media. J. Environ. Qual. 2006, 35, 1018–1025. (14) Bradford, S. A.; Simunek, J.; Walker, S. L. Transport and straining of E. coli O157:H7 in saturated porous media. Water Resour. Res. 2006, 42, W12S12. (15) Levy, J.; Sun, K.; Findlay, R. H.; Farruggia, F. T.; Porter, J.; Mumy, K. L.; Tomaras, J.; Tomaras, A. Transport of Escherichia coli bacteria through laboratory columns of glacial-outwash sediments: Estimating model parameter values based on sediment characteristics. J. Contam. Hydrol. 2007, 89, 71–106. (16) Jiang, G.; Noonan, M. J.; Buchan, G. D.; Smith, N. Transport of Escherichia coli through variably saturated sand columns and modeling approaches. J. Contam. Hydrol. 2007, 93, 2–20. (17) Morrow, J. B.; Stratton, R.; Yang, H.-H.; Smets, B. F.; Grasso, D. Macro- and nanoscale observations of adhesive behavior for several E. coli strains (O157:H7 and environmental isolates) on mineral surfaces. Environ. Sci. Technol. 2005, 39, 6395–6404. (18) Dombek, P. E.; Johnson, L. K.; Zimmerley, S. T.; Sadowsky, M. J. Use of repetitive DNA sequences and the PCR to differentiate Escherichia coli isolates from human and animal sources. Appl. Environ. Microbiol. 2000, 66, 2572–2577. (19) Ram, J. L.; Ritchie, R. P.; Fang, J.; Gonzales, F. S.; Selegean, J. P. Sequence-based source tracking of Escherichia coli based on genetic diversity of Β-Glucuronidase. J. Environ. Qual. 2004, 33, 1024–1032. (20) Clermont, O.; Bonacorsi, S.; Bingen, E. Rapid and simple determination of the Escherichia coli phylogenetic group. Appl. Environ. Microbiol. 2000, 66, 4555–4558. (21) Higgins, J.; Hohn, C.; Hornor, S.; Frana, M.; Denver, M.; Joerger, R. Genotyping Escherichia coli from enviornmental and animal samples. J. Microbiol. Methods 2007, 70, 227–235. (22) Amor, K.; Heinrichs, D. E.; Frirdich, E.; Ziebell, K.; Johnson, R. P.; Whitfield, C. Distribution of core oligosaccharide types in lipopolysaccharides from Escherichia coli. Infect. Immun. 2000, 68, 1116–1124. (23) Yang, H.-H.; Vinopal, R. T.; Grasso, D.; Smets, B. F. High diversity among environmental Escherichia coli isolates from a bovine feedlot. Appl. Environ. Microbiol. 2004, 73, 1528–1536. (24) Kim, H. N.; Bradford, S. A.; Walker, S. L. Escherichia coli O157: H7 transport in saturated porous media: Role of solution chemistry and surface macromolecules. Environ. Sci. Technol. 2009, 43, 4340–4347. (25) Ma, H.; Pedel, J.; Fife, P.; Johnson, W. P. Hemispheres-in-cell geometry to predict colloid deposition in porous media. Environ. Sci. Technol. 2009, 43, 8573–8579. (26) Appelmelk, B. J.; An, Y.-Q.; Hekker, T. A. M.; Thijs, L. G.; MacLaren, D. M.; deGraaf, J. Frequencies of lipopolysaccharide core types in Escherichia coli strains from bacteraemic patients. Microbiology 1994, 140, 1119–1124. (27) Li, J.; McLandsborough, L. A. The effects of the surface charge and hydrophobicity of Escherichia coli on its adhesion to beef muscle. Int. J. Food Microbiol. 1999, 53, 185–193. (28) Walker, S. L.; Redman, J. A.; Elimelech, M. Influence of growth phase on bacterial deposition: Interaction mechanisms in packed-bed column and radial stagnation point flow systems. Environ. Sci. Technol. 2005, 39, 6405–6411. (29) Hong,Y.;Brown,D.B.Cellsurfaceacid-basepropertiesofEscherichia coliandBacillusbrevisandvariationofgrowthphase,nitrogensource and C:N ratio. Colloids Surf., B 2006, 50, 112–119. (30) Huisman, G. W.; Siegele, D. A.; Zambrano, M. M.; Kolter, R. Morphological and physiological changes during stationary phase. Escherichia coli and Salmonella. Cell. Mol. Biol. 1996, 1672–1682. VOL. 44, NO. 13, 2010 / ENVIRONMENTAL SCIENCE & TECHNOLOGY

9

5013

(31) Bolster, C. H.; Mills, A. L.; Hornberger, G. M.; Herman, J. S. Effect of surface coatings, grain size, and ionic strength on the maximum attainable coverage of bacteria on sand surfaces. J. Contam. Hydrol. 2001, 50, 287–305. (32) Mills, A. L.; Herman, J. S.; Hornberger, G. M.; DeJesus, T. H. Effect of solution ionic strength and iron coatings on mineral grains on the sorption of bacterial cells to quartz sand. Appl. Environ. Microbiol. 1994, 60, 3300–3306. (33) Scholl, M. A.; Mills, A. L.; Herman, J. S.; Hornberger, G. M. The influence of mineralogy and solution chemistry on the attachment of bacteria to representative aquifer materials. J. Contam. Hydrol. 1990, 6, 321–336. (34) Li, X.; Scheibe, T. D.; Johnson, W. P. Apparent decreases in colloid deposition rate coefficients with distance of transport under unfavorable deposition conditions: A general phenomenon. Environ. Sci. Technol. 2004, 38, 5616–5625.

5014

9

ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 44, NO. 13, 2010

(35) Walker, S. L.; Redman, J. A.; Elimelech, M. Role of cell surface lipopolysaccharides in Escherichia coli K12 adhesion and transport. Langmuir 2004, 20, 7736–7746. (36) Danese, P. N.; Pratt, L. A.; Dove, S. L.; Kolter, R. The outer membrane protein antigen 43 mediates cell-to-cell interactions within Escherichia coli biofilms. Mol. Microbiol. 2000, 37, 424– 432. (37) Otto, K.; Elwing, H.; Hermansson, M. The role of type 1 fimbriae in adhesion of Escherichia coli to hydrophilic and hydrophobic surfaces. Colloids Surf. 1999, 15, 99–111. (38) Bolster, C. H.; Mills, A. L.; Hornberger, G. M.; Herman, J. S. Effect of intra-population variability on the long distance transport of bacteria. Ground Water 2000, 38, 370–375.

ES1010253