Geospatial Patterns of Antimicrobial Resistance Genes in the US EPA National Rivers and Streams Assessment Survey.

Antimicrobial resistance (AR) is a serious global problem due to the overuse of antimicrobials in human, animal, and agriculture sectors. There is intense research to control the dissemination of AR, but little is known regarding the environmental drivers influencing its spread. Although AR genes (ARGs) are detected in many different environments, the risk associated with the spread of these genes to microbial pathogens is unknown. Recreational microbial exposure risks are likely to be greater in water bodies receiving discharge from human and animal waste in comparison to less disturbed aquatic environments. Given this scenario, research practitioners are encouraged to consider an ecological context to assess the effect of environmental ARGs on public health. Here, we use a stratified, probabilistic survey of nearly 2000 sites to determine national patterns of the anthropogenic indicator class I integron Integrase gene (intI1) and several ARGs in 1.2 million kilometers of United States (US) rivers and streams. Gene concentrations were greater in eastern than in western regions and in rivers and streams in poor condition. These first of their kind findings on the national distribution of intI1 and ARGs provide new information to aid risk assessment and implement mitigation strategies to protect public health.

Exposure to human-associated fecal indicators and self-reported illness among swimmers at recreational beaches: a cohort study.

BACKGROUND: Fecal indicator bacteria used to assess illness risks in recreational waters (e.g., Escherichia coli, Enterococci) cannot discriminate among pollution sources. To address this limitation, human-associated Bacteroides markers have been proposed, but the risk of illness associated with the presence of these markers in recreational waters is unclear. Our objective was to estimate associations between human-associated Bacteroides markers in water and self-reported illness among swimmers at 6 U.S. beaches spanning 2003-2007. METHODS: We used data from a prospectively-enrolled cohort of 12,060 swimmers surveyed about beach activities and water exposure on the day of their beach visit. Ten to twelve days later, participants reported gastroinestinal, diarrheal, and respiratory illnesses experienced since the visit. Daily water samples were analyzed for the presence of human-associated Bacteroides genetic markers: HF183, BsteriF1, BuniF2, HumM2. We used model-based standardization to estimate risk differences (RD) and 95% confidence intervals (CI). We assessed whether the presence of Bacteroides markers were modifiers of the association between general Enterococcus and illness among swimmers using interaction contrast. RESULTS: Overall we observed inconsistent associations between the presence of Bacteroides markers and illness. There was a pattern of increased risks of gastrointestinal (RD = 1.9%; 95% CI: 0.1%, 3.7%), diarrheal (RD = 1.3%; 95% CI: -0.2%, 2.7%), and respiratory illnesses (RD = 1.1%; 95% CI: -0.2%, 2.5%) associated with BsteriF1. There was no evidence that Bacteroides markers acted as modifiers of Enterococcus and illness. Patterns were similar when stratified by water matrix. CONCLUSIONS: Quantitative measures of fecal pollution using Bacteroides, rather than presence-absence indicators, may be necessary to accurately assess human risk specific to the presence of human fecal pollution.

Improved HF183 quantitative real-time PCR assay for characterization of human fecal pollution in ambient surface water samples.

Quantitative real-time PCR (qPCR) assays that target the human-associated HF183 bacterial cluster within members of the genus Bacteroides are among the most widely used methods for the characterization of human fecal pollution in ambient surface waters. In this study, we show that a current TaqMan HF183 qPCR assay (HF183/BFDrev) routinely forms nonspecific amplification products and introduce a modified TaqMan assay (HF183/BacR287) that alleviates this problem. The performance of each qPCR assay was compared in head-to-head experiments investigating limits of detection, analytical precision, predicted hybridization to 16S rRNA gene sequences from a reference database, and relative marker concentrations in fecal and sewage samples. The performance of the modified HF183/BacR287 assay is equal to or improves upon that of the original HF183/BFDrev assay. In addition, a qPCR chemistry designed to combat amplification inhibition and a multiplexed internal amplification control are included. In light of the expanding use of PCR-based methods that rely on the detection of extremely low concentrations of DNA template, such as qPCR and digital PCR, the new TaqMan HF183/BacR287 assay should provide more accurate estimations of human-derived fecal contaminants in ambient surface waters.

Distribution of genetic marker concentrations for fecal indicator bacteria in sewage and animal feces.

Very little is known about the density and distribution of fecal indicator bacteria (FIB) genetic markers measured by quantitative real-time PCR (qPCR) in fecal pollution sources. Before qPCR-based FIB technologies can be applied to waste management and public health risk applications, it is vital to characterize the concentrations of these genetic markers in pollution sources (i.e., untreated wastewater and animal feces). We report the distribution of rRNA genetic markers for several general FIB groups, including Clostridium spp., Escherichia coli, enterococci, and Bacteroidales, as determined by qPCR on reference collections consisting of 54 primary influent sewage samples collected from treatment facilities across the United States and fecal samples representing 20 different animal species. Based on raw sewage sample collection data, individual FIB genetic markers exhibited a remarkable similarity in concentration estimates from locations across the United States ranging from Hawaii to Florida. However, there was no significant correlation between genetic markers for most FIB combinations (P > 0.05). In addition, large differences (up to 5 log(10) copies) in the abundance of FIB genetic markers were observed between animal species, emphasizing the importance of indicator microorganism selection and animal source contribution for future FIB applications.

Comparison of Enterococcus quantitative polymerase chain reaction analysis results from fresh and marine waters on two real-time instruments.

The U.S. Environmental Protection Agency will be recommending a quantitative polymerase chain reaction (qPCR) method targeting Enterococcus spp. as an option for monitoring recreational beach water quality. A practical consideration for implementation of this and other qPCR methods is whether the results are comparable on different PCR instruments. In this study, quantitative estimates of Enterococcus densities from marine and freshwater samples were determined by the qPCR method from cycle threshold (Ct) measurements obtained on Applied Biosystems StepOnePlus and Cepheid SmartCycler instruments. Three variations of a comparative Ct model, differing in their sources of calibration data, were used in the estimations. Both traditional and Bayesian statistical modeling approaches were examined in the instrument comparisons. The traditional analysis of variance (ANOVA) approach indicated no significant differences (p>0.05) between mean density estimates from the instruments in two of the three model variations. The Bayesian approach indicated that the 95% Bayesian credible intervals of density estimates from the instruments overlapped in all models; however, the uncertainty of the estimates varied depending on the model. These results support the interchangeable use of the two instruments in the method and also illustrate the importance of defining the source of calibration data used in the comparative Ct model.

Interlaboratory comparison of real-time PCR protocols for quantification of general fecal indicator bacteria.

The application of quantitative real-time PCR (qPCR) technologies for the rapid identification of fecal bacteria in environmental waters is being considered for use as a national water quality metric in the United States. The transition from research tool to a standardized protocol requires information on the reproducibility and sources of variation associated with qPCR methodology across laboratories. This study examines interlaboratory variability in the measurement of enterococci and Bacteroidales concentrations from standardized, spiked, and environmental sources of DNA using the Entero1a and GenBac3 qPCR methods, respectively. Comparisons are based on data generated from eight different research facilities. Special attention was placed on the influence of the DNA isolation step and effect of simplex and multiplex amplification approaches on interlaboratory variability. Results suggest that a crude lysate is sufficient for DNA isolation unless environmental samples contain substances that can inhibit qPCR amplification. No appreciable difference was observed between simplex and multiplex amplification approaches. Overall, interlaboratory variability levels remained low (<10% coefficient of variation) regardless of qPCR protocol.

Large-scale comparison of E. coli levels determined by culture and a qPCR method (EPA Draft Method C) in Michigan towards the implementation of rapid, multi-site beach testing.

Fecal pollution remains a challenge for water quality managers at Great Lakes and inland recreational beaches. The fecal indicator of choice at these beaches is typically Escherichia coli (E. coli), determined by culture-based methods that require over 18 h to obtain results. Researchers at the United States Environmental Protection Agency (EPA) have developed a rapid E. coli qPCR methodology (EPA Draft Method C) that can provide same-day results for improving public health protection with demonstrated sensitivity, specificity, and data acceptance criteria. However, limited information is currently available to compare the occurrence of E. coli determined by cultivation and by EPA Draft Method C (Method C). This study provides a large-scale data collection effort to compare the occurrence of E. coli determined by these alternative methods at more than 100 Michigan recreational beach and other sites using the complete set of quantitative data pairings and selected subsets of the data and sites meeting various eligibility requirements. Simple linear regression analyses of composite (pooled) data indicated a correlation between results of the E. coli monitoring approaches for each of the multi-site datasets as evidenced by Pearson R-squared values ranging from 0.452 to 0.641. Theoretical Method C threshold values, expressed as mean log10 target gene copies per reaction, that corresponded to an established E. coli culture method water quality standard of 300 MPN or CFU /100 mL varied only from 1.817 to 1.908 for the different datasets using this model. Different modeling and derivation approaches that incorporated within and between-site variability in the estimates also gave Method C threshold values in this range but only when relatively well-correlated datasets were used to minimize the error. A hypothetical exercise to evaluate the frequency of water impairments based on theoretical qPCR thresholds corresponding to the E. coli water quality standard for culture methods suggested that the methods may provide the same beach notification outcomes over 90% of the time with Method C results differing from culture method results that indicated acceptable and unacceptable water quality at overall rates of 1.9% and 6.6%, respectively. Results from this study provide useful information about the relationships between E. coli determined by culture and qPCR methods across many diverse freshwater sites and should facilitate efforts to implement qPCR-based E. coli detection for rapid recreational water quality monitoring on a large scale in the State of Michigan.

Performance assessment PCR-based assays targeting bacteroidales genetic markers of bovine fecal pollution.

There are numerous PCR-based assays available to characterize bovine fecal pollution in ambient waters. The determination of which approaches are most suitable for field applications can be difficult because each assay targets a different gene, in many cases from different microorganisms, leading to variation in assay performance. We describe a performance evaluation of seven end-point PCR and real-time quantitative PCR (qPCR) assays reported to be associated with either ruminant or bovine feces. Each assay was tested against a reference collection of DNA extracts from 247 individual bovine fecal samples representing 11 different populations and 175 fecal DNA extracts from 24 different animal species. Bovine-associated genetic markers were broadly distributed among individual bovine samples ranging from 39 to 93%. Specificity levels of the assays spanned 47.4% to 100%. End-point PCR sensitivity also varied between assays and among different bovine populations. For qPCR assays, the abundance of each host-associated genetic marker was measured within each bovine population and compared to results of a qPCR assay targeting 16S rRNA gene sequences from Bacteroidales. Experiments indicate large discrepancies in the performance of bovine-associated assays across different bovine populations. Variability in assay performance between host populations suggests that the use of bovine microbial source-tracking applications will require a priori characterization at each watershed of interest.

Performance of PCR-based assays targeting Bacteroidales genetic markers of human fecal pollution in sewage and fecal samples.

There are numerous PCR-based assays available to characterize human fecal pollution in ambient waters. Each assay employs distinct oligonucleotides and many target different genes and microorganisms leading to potential variations in assay performance. Performance comparisons utilizing feces and raw sewage samples are needed to determine which assays are best suited for expensive and time-consuming field validation, fate, transport, and epidemiology studies. We report the assessment of five end-point PCR and 10 real-time quantitative PCR (qPCR) assays that target genes from presumptive Bacteroidales microorganisms reported to be associated with human feces. Each assay was tested against a reference collection of 54 primary influent sewage samples collected from different geographical locations across the United States and 174 fecal DNA extracts from 23 different animal sources. Experiments indicate that human-associated genetic markers are distributed across a broad range of human populations but show substantial differences in specificity for human feces suggesting that particular assays may be more suitable than others depending on the abundance of genetic marker required for detection and the animal sources impacting a particular watershed or beach of interest.

Quantitative PCR for genetic markers of human fecal pollution.

Assessment of health risk and fecal bacterial loads associated with human fecal pollution requires reliable host-specific analytical methods and a rapid quantification approach. We report the development of quantitative PCR assays for quantification of two recently described human-specific genetic markers targeting Bacteroidales-like cell surface-associated genes. Each assay exhibited a range of quantification from 10 to 1 x 10(6) copies of target DNA. For each assay, internal amplification controls were developed to detect the presence or absence of amplification inhibitors. The assays predominantly detected human fecal specimens and exhibited specificity levels greater than 97% when tested against 265 fecal DNA extracts from 22 different animal species. The abundance of each human-specific genetic marker in primary effluent wastewater samples collected from 20 geographically distinct locations was measured and compared to quantities estimated by real-time PCR assays specific for rRNA gene sequences from total Bacteroidales and enterococcal fecal microorganisms. Assay performances combined with the prevalence of DNA targets in sewage samples provide experimental evidence supporting the potential application of these quantitative methods for monitoring fecal pollution in ambient environmental waters.

Sticks and stones: racism as experienced by adolescents in New Brunswick.

The purpose of the study was to describe both the nature of racism as experienced by adolescent self-described victims in the province of New Brunswick, Canada, and their response to the perceived racist incidents. A qualitative methodology based on the constructivist paradigm was used. In-depth interviews were conducted with non-White adolescent victims of racism and with parents of victims. Although the study was initiated in response to an eruption of publicity about teenage racial violence, the findings indicate that racist incidents were not a new phenomenon for the participants. They described a low-key but long-term problem that had begun when they entered the public school system. Name-calling was by far the most common form of racism identified, and it played a part in most of the other incidents described; dismissed as harmless by authority figures, it appeared to have long-term consequences for its targets.The participants' response to racism was found to have 3 phases: splintered universe, spiralling resistance, and disengagement. The results suggest that nurses working in the field of school health should address issues of racism among children and adolescents.

A Bayesian method for calculating real-time quantitative PCR calibration curves using absolute plasmid DNA standards.

BACKGROUND: In real-time quantitative PCR studies using absolute plasmid DNA standards, a calibration curve is developed to estimate an unknown DNA concentration. However, potential differences in the amplification performance of plasmid DNA compared to genomic DNA standards are often ignored in calibration calculations and in some cases impossible to characterize. A flexible statistical method that can account for uncertainty between plasmid and genomic DNA targets, replicate testing, and experiment-to-experiment variability is needed to estimate calibration curve parameters such as intercept and slope. Here we report the use of a Bayesian approach to generate calibration curves for the enumeration of target DNA from genomic DNA samples using absolute plasmid DNA standards. RESULTS: Instead of the two traditional methods (classical and inverse), a Monte Carlo Markov Chain (MCMC) estimation was used to generate single, master, and modified calibration curves. The mean and the percentiles of the posterior distribution were used as point and interval estimates of unknown parameters such as intercepts, slopes and DNA concentrations. The software WinBUGS was used to perform all simulations and to generate the posterior distributions of all the unknown parameters of interest. CONCLUSION: The Bayesian approach defined in this study allowed for the estimation of DNA concentrations from environmental samples using absolute standard curves generated by real-time qPCR. The approach accounted for uncertainty from multiple sources such as experiment-to-experiment variation, variability between replicate measurements, as well as uncertainty introduced when employing calibration curves generated from absolute plasmid DNA standards.

Quantification of plasmid DNA standards for U.S. EPA fecal indicator bacteria qPCR methods by droplet digital PCR analysis.

An obstacle to establishing widely useful data acceptance criteria for U.S. Environmental Protection Agency (EPA) qPCR methods has been the unavailability of standardized reference materials. Earlier versions of EPA Methods 1609 and 1611 for enterococci used cellular reference materials for quantifying enterococci in unknown test samples, however, EPA updates to these fundamentally DNA-based analysis methods have shifted toward the use of DNA standards. This report describes the application of droplet digital PCR (ddPCR) analysis for the quantification of a set of synthetic plasmid DNA standards that have been made available for updated EPA Methods 1609.1 and 1611.1 as well as for EPA Draft Method C for Escherichia coli. To obtain the most accurate concentration estimates possible, part of this effort was to develop a data analysis model for determining the fluorescence thresholds that distinguish positive from negative droplets produced by the ddPCR reactions. Versions of this model are described for applications with individual reactions, multiple reactions within a ddPCR system run, and multiple reactions within and across different system runs. The latter version was applied toward determinations of error in the concentration estimates of the standards from replicate analyses of each standard in multiple ddPCR system runs. Mean concentration estimates for the five standards from the ddPCR analyses were 4.356, 3.381, 2.371, 1.641 and 1.071 log10 copies/5 µL with associated standard deviations of 0.074, 0.082, 0.108, 0.131 and 0.188, respectively. These estimates contrasted with expected log10 concentrations of 4.6, 3.6, 2.6, 1.9 and 1.3 copies/5 µL, respectively, based on the yield of the plasmid reported by the vendor and spectrophotometric analysis of the initial stock solution of this material. These results illustrate how the analyses of original stocks may lead to potential bias(es) in the concentration estimates of final DNA standards and subsequently in the estimates of unknown test samples determined from these standards in qPCR analyses.

Multi-laboratory survey of qPCR enterococci analysis method performance in U.S. coastal and inland surface waters.

Quantitative polymerase chain reaction (qPCR) has become a frequently used technique for quantifying enterococci in recreational surface waters, but there are several methodological options. Here we evaluated how three method permutations, type of mastermix, sample extract dilution and use of controls in results calculation, affect method reliability among multiple laboratories with respect to sample interference. Multiple samples from each of 22 sites representing an array of habitat types were analyzed using EPA Method 1611 and 1609 reagents with full strength and five-fold diluted extracts. The presence of interference was assessed three ways: using sample processing and PCR amplifications controls; consistency of results across extract dilutions; and relative recovery of target genes from spiked enterococci in water sample compared to control matrices with acceptable recovery defined as 50 to 200%. Method 1609, which is based on an environmental mastermix, was found to be superior to Method 1611, which is based on a universal mastermix. Method 1611 had over a 40% control assay failure rate with undiluted extracts and a 6% failure rate with diluted extracts. Method 1609 failed in only 11% and 3% of undiluted and diluted extracts analyses. Use of sample processing control assay results in the delta-delta Ct method for calculating relative target gene recoveries increased the number of acceptable recovery results. Delta-delta tended to bias recoveries from apparent partially inhibitory samples on the high side which could help in avoiding potential underestimates of enterococci--an important consideration in a public health context. Control assay and delta-delta recovery results were largely consistent across the range of habitats sampled, and among laboratories. The methodological option that best balanced acceptable estimated target gene recoveries with method sensitivity and avoidance of underestimated enterococci densities was Method 1609 without extract dilution and using the delta-delta calculation method. The applicability of this method can be extended by the analysis of diluted extracts to sites where interference is indicated but, particularly in these instances, should be confirmed by augmenting the control assays with analyses for target gene recoveries from spiked target organisms.

Development and evaluation of a culture-independent method for source determination of fecal wastes in surface and storm waters using reverse transcriptase-PCR detection of FRNA coliphage genogroup gene sequences.

A method, incorporating recently improved reverse transcriptase-PCR primer/probe assays and including controls for detecting interferences in RNA recovery and analysis, was developed for the direct, culture-independent detection of genetic markers from FRNA coliphage genogroups I, II & IV in water samples. Results were obtained from an initial evaluation of the performance of this method in analyses of waste water, ambient surface water and stormwater drain and outfall samples from predominantly urban locations. The evaluation also included a comparison of the occurrence of the FRNA genetic markers with genetic markers from general and human-related bacterial fecal indicators determined by current or pending EPA-validated qPCR methods. Strong associations were observed between the occurrence of the putatively human related FRNA genogroup II marker and the densities of the bacterial markers in the stormwater drain and outfall samples. However fewer samples were positive for FRNA coliphage compared to either the general bacterial fecal indicator or the human-related bacterial fecal indicator markers particularly for ambient water samples. Together, these methods show promise as complementary tools for the identification of contaminated storm water drainage systems as well as the determination of human and non-human sources of contamination.

Detection of Cyclospora cayetanensis using a quantitative real-time PCR assay.

Cyclospora cayetanensis, a coccidian parasite, with a fecal-oral life cycle, has become recognized worldwide as an emerging human pathogen. Clinical manifestations include prolonged gastroenteritis. While most cases of infection with C. cayetanensis in the United States have been associated with foodborne transmission, waterborne transmission has also been implicated. We report on the development and application of a real-time, quantitative polymerase chain reaction assay for the detection of C. cayetanensis oocysts, which is the first reported use of this technique for this organism. Both a species-specific primer set and dual fluorescent-labeled C. cayetanensis hybridization probe were designed using the inherent genetic uniqueness of the 18S ribosomal gene sequence of C. cayetanensis. The real-time polymerase chain reaction assay has been optimized to specifically detect the DNA from as few as 1 oocyst of C. cayetanensis per 5 microl reaction volume.

Quantitative PCR analysis of selected Aspergillus, Penicillium and Paecilomyces species.

A total of 65 quantitative PCR (QPCR) assays, incorporating fluorigenic 5' nuclease (TaqMan) chemistry and directed at the nuclear ribosomal RNA operon, internal transcribed spacer regions (ITS1 or ITS2) was developed and tested for the detection of selected Aspergillus, Penicillium and Paecilomyces species. The assays varied in specificity from species or subspecies to closely related species groups, subject to the amount of nucleotide sequence variation in the different organisms. A generic assay for all target species of Aspergillus, Penicillium and Paecilomyces was also developed and tested. Using a previously reported DNA extraction method, estimated conidia detection limits for target species ranged from less than one to several hundred per sample for the different assays. Conidia detection limits for non-target species were at least 1,000 fold higher in nearly all instances. The assays were used to analyze ten HVAC dust samples from different sources around the US. Total quantities of Aspergillus, Penicillium and Paecilomyces conidia in the samples, determined by the generic assay and the summed totals from the specific assays, were in general agreement, suggesting that all of the numerically dominant species in the samples were accounted for by the specific assays. QPCR analyses of these samples after spiking them with selected target organisms indicated that the enumeration results were within approximately a one-half log range of the expected values 95% of the time. Evidence is provided that the commonly used practices of enumerating Aspergillus and Penicillium as a single group or only by genus can be misleading in understanding the indoor populations of these organisms and their potential health risks.

Quantitative polymerase chain reaction analysis of fungi in dust from homes of infants who developed idiopathic pulmonary hemorrhaging.

Fungal concentrations were measured in the dust of 6 homes in Cleveland, Ohio, where an infant developed pulmonary hemorrhage (pulmonary hemorrhage homes [PHH]) and 26 reference homes (RH) with no known fungal contamination. Quantitative polymerase chain reaction assays for 82 species (or assay groups) were used to identify and quantify fungal concentrations. The ratios of the geometric means of PHH to RH were >1 for 26 species (group I). However, the same ratios were <1 for 10 species (group II). Probit analysis of the sum of the logs of the concentrations of these 2 groups resulted in a 95% probability range for separating PHH from RH homes. The same 82 fungal species were also tested for hemolysin production on sheep's blood agar (incubated at 37 degree C). Hemolysins were more commonly produced by group I species (42%) compared with group II species (10%).

Comparison of Enterococcus quantitative polymerase chain reaction analysis results from Midwest U.S. river samples using EPA Method 1611 and Method 1609 PCR reagents.

Enterococci target sequence density estimates from analyses of diluted river water DNA extracts by EPA Methods 1611 and 1609 and estimates with lower detection limits from undiluted DNA extracts by Method 1609 were indistinguishable. These methods should be equally suitable for comparison with U.S. EPA 2012 Recreational Water Quality Criteria values.

Standardization of enterococci density estimates by EPA qPCR methods and comparison of beach action value exceedances in river waters with culture methods.

The U.S. EPA has published recommendations for calibrator cell equivalent (CCE) densities of enterococci in recreational waters determined by a qPCR method in its 2012 recreational water quality criteria (RWQC). The CCE quantification unit stems from the calibration model used to estimate enterococci densities in recreational beach waters in the EPA National Epidemiological and Environmental Assessment of Recreational (NEEAR) Water Study and directly informed the derivation of the RWQC recommendations. Recent studies have demonstrated that CCE estimates from the method can vary when using different cultured Enterococcus cell preparations in calibrator samples. These differences have been attributed to differences in the quantities of targeted gene copies (target sequences) that are recovered per nominal calibrator cell by DNA extraction. Standardization of results from the calibration model will require the estimation of target sequence recoveries from the calibrator and water samples. In addition, comparisons of water sample results with the RWQC values will require a knowledge of target sequence recoveries from the NEEAR study calibrator samples. In this study recoveries of target sequences and the mean target sequence/cell ratio for the NEEAR study calibrator samples were retrospectively estimated with a corroborated standard curve. A modification of the calibration model was then used to estimate enterococci target sequence quantities in water samples from eight midwestern U.S. rivers. CCE estimates were obtained by dividing these target sequence quantities by the mean NEEAR study target sequence/cell ratio. This target sequence-based quantification approach resulted in a high degree of agreement in beach action decisions (determinations of whether bacterial fecal indicator densities are above or below RWQC-recommended values) from CCE results of the qPCR method and from culture dependent enumeration of both enterococci and Eschericia coli in the corresponding water samples.