Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2005, p. 5992-5998, Vol. 71, No. 10
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.10.5992-5998.2005
Agriculture and Agri-Food Canada, London, Ontario, Canada N5V 4T3,1 Agriculture and Agri-Food Canada, Central Experimental Farm, Ottawa, Ontario, Canada K1A 0C62
Received 24 December 2004/ Accepted 24 May 2005
|
|
|---|
|
|
|---|
This is particularly true when freshwater resources are in proximity to land subject to increasing agricultural activity and burgeoning human populations, increasing the risk to adjacent waters from agricultural runoff, sewage effluent, leaking rural septic systems, and stormwater discharge. The microbiological quality of surface and drinking water has traditionally been evaluated by quantifying fecal indicator bacteria, notably Escherichia coli, using standard microbiological methods. The presence of this organism is implicit evidence for fecal contamination and indicates a possible risk of contamination with viral, bacterial, or parasitic pathogens of enteric origin. On this basis, many jurisdictions mandate compliance with drinking and recreational water standards (4, 5).
A watershed management approach has been proposed for protecting surface water quality (12, 28). A key requirement of this strategy is to identify and then eliminate or abate sources of significant fecal contamination. In watersheds with mixed urban, agricultural, and industrial activities, the identification of pollution sources can pose a significant challenge. There has thus been significant interest in identifying attributes of fecal indicator bacteria that distinguish the host source. Genotypic and phenotypic approaches that have been investigated for this purpose include repetitive extragenic palindromic-type PCR fingerprinting, ribotyping, AFLP, pulsed-field gel electrophoresis, detection of source-specific marker genes, antibiotic resistance profiles, and carbon utilization profiles (6, 8, 13, 15, 16, 19-21, 23, 27). Repetitive extragenic palindromic PCR methods are a method of choice because of the relatively low cost, operational ease, and success at correctly classifying the host source (3, 20). Typically with this approach, environmental isolates are compared with reference collections of bacteria obtained from potential sources of fecal pollution in the area and on the basis of similarity are ascribed a probable host source. For any given study, the likelihood of correct classification will depend on several factors, including the size and representability of the reference collection, geographic size of the study area, choice of methods for image analysis and pattern recognition, and statistical methods for comparison of environmental isolates with the reference collection (1, 8, 13, 19-22, 27). The construction of robust fingerprint libraries is expensive and therefore requires an informed sampling strategy.
In many types of farming systems, animals or poultry are raised confined in barns, and their manure is stored for several months prior to release into the environment. Swine, dairy, and egg-laying poultry operations notably store excreta as a liquid or slurry in sometimes extremely large manure holding tanks. In Canada, for example, about 85% of swine and 43% of dairy cattle are produced on farms that use liquid manure storage systems (25). These manure holding tanks therefore represent a key sampling point for construction of reference collections in areas with confined production systems. Furthermore, the strain composition of enteric bacterial communities shed by farm animals could be altered significantly during manure storage, biasing the genetic composition of populations released into the broader environment. In this context, we (i) compared the structure and diversity of E. coli collections obtained from manure holding tanks with those of concurrent collections obtained from fresh manure shed by the herds on a dairy farm and a swine farm, (ii) examined the E. coli populations in swine and dairy manure holding tanks at different collection times to assess potential for population differences over time, and (iii) evaluated the spatial variability in the E. coli community composition within large manure holding tanks.
|
|
|---|
Approximately 500 g of freshly excreted fecal material were obtained from each animal in the dairy barn during morning feeding. Swine in-barn samples were taken by collecting a composite sample of approximately 500 g of freshly excreted fecal material from pens representing the various age groups in the barn. Sampling of individuals in the swine barn was not possible due to the number of individuals per pen and the type of pens used in the barn. Holding tank samples were collected in 1-liter sterile bottles (Systems Plus, Woodstock, Ontario, Canada) at a discrete depth of either 2.5 or 0.5 m below the surface, using a Sludge Judge Ultra sampler (NASCO Canada, Aurora, Ontario, Canada). Ten samples were taken around the circumference of the holding tank, five of which were taken at 2-foot depth and five of which were taken at 9-foot depth.
In-barn samples were pooled by thoroughly mixing each individual sample and transferring 100 g into a clean Ziploc bag. This composite sample was then mixed thoroughly prior to any subsampling for microbiological purposes.
E. coli isolation and identification.
Manure samples were kept at 4°C and processed within 24 h. Samples were serially diluted in sodium metaphosphate buffer (2 g/liter; Fisher, Mississauga, Ontario, Canada) and mixed thoroughly. Dilutions were spread plated onto mFC basal medium (Difco, Toronto, Ontario, Canada) supplemented with 100 mg/liter of 3-bromo-4-chloro-5-indolyl-ß-D-glucopyranoside (BCIG) (Med-Ox Diagnostics, Ottawa, Ontario, Canada) and incubated overnight at 44.5°C. After overnight growth, single blue colonies were picked and streaked onto LB agar (Difco, Toronto, Ontario, Canada) (four isolates/plate) and grown at 37°C overnight. The isolates were purified by restreaking twice on LB agar. The purified colonies were inoculated into sterile 96-well microtiter plates containing 100 µl fresh LB broth (Difco, Toronto, Ontario, Canada) per well and grown statically overnight at 37°C. For confirmation, the isolate cultures were replica plated (10 µl per well) into sterile 96-well microtiter plates containing 100 µl of lactose broth (containing, per liter, 10 g Proteose Peptone no. 3, 3 g yeast extract, 5 g NaCl, 10 g lactose, and 20 mg bromcresol purple) or 100 µl of tryptone broth (containing, per liter, 10 g Bacto tryptone and 5 g NaCl) and were incubated overnight at 37°C. Positive confirmation was indicated by a color change from purple to yellow in lactose broth (lactose fermentation) and by the formation of a red-pink color upon addition of 40 µl of Kovac's reagent to the tryptone broth wells, indicating indole production (24). Isolates were considered to be Escherichia coli if they grew at 44.5°C, had a positive reaction for ß-glucuronidase (blue color on mFC-BCIG agar), fermented lactose, and produced indole. Confirmed isolates were inoculated into sterile 96-well microplates containing 100 µl/well of LB broth and incubated overnight at 37°C. Sterile glycerol (Sigma, Mississauga, Ontario, Canada) was then added to each well at a final concentration of 15% (vol/vol), and the plates were stored at 70°C.
Template preparation and ERIC PCR.
Cell suspensions of E. coli were prepared by inoculating 100 µl of fresh LB broth per well in a sterile 96-well microtiter plate with frozen stock cultures. Cells were grown statically at 37°C overnight to an A600 of about 1 and centrifuged at 710 x g for 25 min (Centra CL3 microplate centrifuge; Thermo IEC, Needham Heights, MA). The pelleted cells were resuspended in 100 µl of sterile Milli-Q H2O (100 µl) and shaken at 1,000 rpm with a microplate shaker (Sarstedt, Montreal, Quebec, Canada) for 5 min. The resuspended cells were used directly as template for the PCR or frozen at 20°C until required.
Primers used for enterobacterial repetitive intergenic consensus (ERIC) PCR were the same as described by Versalovic et al. (29). The final reaction mix (25µl) consisted of 1x PCR buffer (Promega, Madison, WI), 3 mM MgCl2, 0.1 mg/ml gelatin, 200 µM of each deoxynucleoside triphosphate (Invitrogen, Burlington, Ontario, Canada), 2 µM each of forward and reverse primers ERIC-1 and ERIC, 1 U of Taq polymerase (Promega), and 2 µl of E. coli suspended cells as template. Amplification was performed in a Hybaid OmniGene thermocycler (InterSciences Inc., Markham, Ontario, Canada) as follows: after an initial denaturation at 95°C for 10 min, 34 cycles of denaturation (94°C, 3 seconds), (92°C, 30 seconds), annealing (50°C, 1 min), and extension (65°C, 1 min) were performed, followed by a final extension (65°C, 8 min).
PCR products were resolved by horizontal gel electrophoresis in a 25-cm by 50-cm gel (Gator A3-1; Owl Separations, Portsmouth, NH) prepared with 1.5% (wt/vol) agarose (Invitrogen, Mississauga, Ontario, Canada) and 1x Tris-borate-EDTA buffer. Six microliters of loading dye was added to 25 µl of PCR product, and 7 µl of this mixture was loaded into wells prepared with an 8-mm by 1-mm comb tooth size. Every eighth well received the MassRuler DNA ladder (Fermentas, Burlington, Ontario, Canada). Gels were subjected to 4 V/cm for 2.5 h in 1x Tris-borate-EDTA. The gel was stained with 1 µg/ml ethidium bromide solution for 10 min and destained in Milli-Q water for 10 min. Gel images were captured as 8-bit TIFF images, using Quantity One gel documentation software (Bio-Rad, Mississauga, Ontario, Canada) with a CCD gel documentation system (Bio-Rad, Mississauga, Ontario, Canada).
Computer-assisted image and data analysis.
Normalization of gel images and assignment of fingerprints to isolates were done with Bionumerics (version 3.5; Applied Maths, Kortrijk, Belgium). Positions of fingerprints on gels were normalized using the MassRuler DNA ladder as the external standard in the range of 300 bp to 3,000 bp. Similarity coefficients were generated using the curve-based cosine correlation coefficient. Similarity trees were generated using the unweighted-pair group method using average linkage, and a similarity cutoff of 80% was used in order to determine related fingerprint types. Fingerprint types and numbers of isolates per fingerprint type were tabulated in Microsoft Excel. The diversity captured in the E. coli collections was estimated by rarefaction analysis using the analytical approximation algorithm of Hurlbert (11) and 95% confidence intervals estimated as described by Heck et al. (9). Calculations were performed with the freeware program Analytical Rarefaction 1.3, available at http://www.uga.edu/
strata/software/. Curves were plotted using SigmaPlot (version 8.02; SPSS Inc., Chicago, IL). The asymptotes of the rarefaction curves were estimated using the Michaelis-Menten equation, which is available in SigmaPlot as the one-site saturation ligand model (10). The asymptote is a measure of richness at sampling saturation and was used to estimate the fraction of total community diversity captured within our E. coli collections. The SigmaPlot curve fitter uses the Marquardt-Levenberg algorithm to find the coefficients that give the best fit between the equation and the data (17).
Evaluation of spatial distribution variability in manure holding tanks.
Five sampling sites were located horizontally at about 70o angles from each other and about 1 m from the holding tank edge. At each site 1-liter samples were taken from depths of 0.5 m and 2.5 m. A Classification and Regression Trees (CART)-based classification tree approach was used to classify the dominant fingerprints on the basis of sampling depth and lateral sampling location (2, 26). CART is a widely accepted automated, binary recursive partitioning technique that selects predictor variables (independent variables) and their interactions that optimally predict a dependent measure. Dominant fingerprints in both the dairy and swine data sets were defined as being those with
30 observations. The classification trees were produced using a Gini splitting criterion. Each fingerprint class was also treated as if it was uniformly distributed in the population regardless of the observed sample proportions, essentially treating each fingerprint class as equally important for classification accuracy purposes. Misclassification rates were determined using a CART-based cross-validation procedure. In this study, a 10-fold approach was employed.
|
|
|---|
![]() View larger version (30K): [in a new window] |
FIG. 1. Rarefaction curves indicating the relative richness of E. coli collections obtained from fresh manure and the manure holding tank from a swine farm and a dairy farm (n = 500 in each case). All of the collections obtained in December 2003 were obtained on the same day. In order to evaluate the temporal stability of the populations in the stored manure, the holding tanks were subsequently resampled.
|
|
View this table: [in a new window] |
TABLE 1. Estimated total richness in E. coli populations
|
The distributions of genotypes in the swine collection sources were compared (Fig. 2). The 2003 swine barn, 2003 holding tank, and 2004 holding tank collections had 19, 39, and 22 unique genotypes (i.e., represented by only one isolate), respectively. In general, the genotypes that were most well represented in the barn were also dominant in the holding tank at both sampling times. Taking into account the relative abundance of each distinct genotype, 63% of the isolates from the December 2003 holding tank collection were represented in the barn collection. Eighty-four percent of isolates from the March 2004 holding tank collection were represented in the barn collection. Overall, the majority of the isolates obtained from the holding tank were also detected in the barn, but each collection had a significant number of unique genotypes.
![]() View larger version (22K): [in a new window] |
FIG. 2. Frequency of occurrence of genotypes (ERIC fingerprint types) of E. coli obtained from the barn and the manure holding tank of a swine farm (n = 500 in each case). The holding tank was sampled on two different dates.
|
![]() View larger version (18K): [in a new window] |
FIG. 3. Frequency of occurrence of genotypes (ERIC fingerprint types) of E. coli obtained from the barn and the manure holding tank of a dairy farm (n = 500 in each case).
|
![]() View larger version (25K): [in a new window] |
FIG. 4. Rarefaction curves for manure collections obtained from the manure holding tanks at a depth of 0.5 m or 2.5 m. Samples taken from these depths at five different lateral locations were pooled (n = 500).
|
![]() View larger version (31K): [in a new window] |
FIG. 5. Rarefaction curves for manure collections obtained from the manure holding tanks at five different lateral locations. Samples from depths of 0.5 and 2.5 m at each lateral location were pooled (n= 200).
|
There was excellent classification (percentage classified correct) on the basis of spatial location in the holding tank for genotypes designated 17, 66, and 67. Approximately 87% (for both testing and learning sets) of genotype 67 observations were located at one 0.5-m-depth location in the tank, while 97% (for both testing and learning sets) of genotype 66 observations were found at one 2.5-m-depth location in the tank. However, for genotype 17, where there was excellent classification prediction (93% for learning and 88% for testing), the fingerprint was found at three distinct lateral locations (two at 0.5-m depth and one at 2.5-m depth). Variable importance as defined on the basis of the primary splitting variables in the classification indicated that lateral location and depth were essentially of equal relative importance in discriminating the genotype data (2).
For the swine manure holding tank, as with the dairy tank, the optimal classification tree consisted of genotype classes representing each of the lateral and depth locales (total of 10 sites). There were a total of 739 observations. However, there was much poorer classification accuracy on the basis of the independent variables, relative to the dairy data set. The greatest prediction success occurred for the genotype designated 19, at only 50% correct classification for both the learning and test sample sets. The variable importance indicated dramatic differences between lateral location and depth as classification tree predictor variables (lateral location was scaled at 100% importance, relative to 19% relative importance for depth).
|
|
|---|
Our sample size captured from 60% (dairy holding tank, December 2003) to 85% (swine barn, December 2003) of the predicted genotypes in the manures (Table 1). The number of isolates required to capture 50% of the predicted numbers of genotypes in the manures ranged from 111 to 353 isolates. Clearly, on a watershed that has many confined livestock production farms, the number of isolates that will be required to build representative sample libraries for this potential source will be very large. This finding is in agreement with the large E. coli population diversity now having been characterized in other potential sources of fecal pollution (7, 8, 13, 18). In any MST study, the choice of source reference library size will be determined by the size and complexity of the study area, the accuracy of source identity required by the investigators, and the financial resources available (13, 19, 21, 27).
In both farms studied here, the diversity in the E. coli community was higher in the December samples than those of the following March (swine) or July (dairy) (Fig. 1). We did not simultaneously sample the barns in the later samplings and therefore do not know if the holding tank communities are reflective of changes in the community shed by the herd at those times. In southern Ontario, Canada, the stored manure is typically at a temperature of 0 to 5°C in December and in the mid-20oCs in the summer (data not shown). It is highly likely that the E. coli population is more dynamic at warmer temperatures. A cursory examination of the genotype distributions (Fig. 2 and 3) suggests that libraries constructed from the herd would miss a number of the genotypes found at low representation in the holding tank. In the swine farm, the dominant genotypes found in the barn in December were also well represented in the holding tank in both December and March. In the dairy farm, the December barn and holding tank communities were heavily dominated by one genotype, whereas the July holding tank community was heavily dominated by genotypes that were relatively underrepresented in the December samples. Overall, these data suggest that the composition of E. coli populations released into the environment following manure storage will be variable, subject to seasonal effects that will vary according to the local climate, and that libraries constructed from that farm at any time may misrepresent the composition of E. coli populations subsequently released by that farm into the environment.
We detected clear differences in dominant genotype spatial structure between the swine and dairy manure tank data. For dairy, there are distinct locations in the tanks that harbor specific genotypes, and a sampling strategy should therefore ensure lateral and vertical sampling components. However, for swine, the results here suggested that a bulk sample from any location would likely capture the dominant fingerprints. The more homogeneous population distribution in the swine manure holding tank may be due to the aeration that it received and whatever agitation this provided. Nevertheless, from a conservative perspective, it is likely best that a holding tank sampling scheme include lateral and vertical sampling components.
An intrinsic dilemma that plagues many environmentally based sampling programs is the development of a sampling design that captures spatial and temporal variability while maintaining logistical feasibility. This dilemma is underscored for watershed-scale MST studies, where it is critical that as many representative fecal sources are sampled as feasible and that the numbers of bacteria obtained be sufficiently large to be representative of what will be released into the broader environment. This study indicated that even at the scale of a single source point, diversity and genotype representation can vary with respect to spatial location in manure storage facilities as well as with respect to the time the sample was taken. Overall, the size of reference libraries required to capture the genetic diversity of E. coli from potential fecal sources on a watershed scale is one of a number of operational constraints that limit the applicability and likely accuracy of MST methods that require host source reference libraries (8, 13, 19, 27).
We sincerely thank K. and K. Nagelschmitz and C. and M. Bontja for access to their farms. L. Sabourin and P. Bastedo assisted with farm sampling. A. Lachance provided valuable advice.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»