Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2006, p. 1604-1612, Vol. 72, No. 2
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.2.1604-1612.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
University of North CarolinaChapel Hill, Chapel Hill, North Carolina,1 Southern California Coastal Water Research Project, Westminster,2 University of Southern California, Los Angeles,3 Santa Monica Baykeeper, Santa Monica, California4
Received 16 June 2005/ Accepted 16 November 2005
|
|
|---|
|
|
|---|
Despite the impairment of water quality and risks to human health, identification and elimination of the sources of bacteria responsible for the beach warnings remain elusive. The difficulty in identifying and eliminating the sources of bacteria results from three important factors. First, the traditional indicators of fecal pollution on the basis of which the water quality thresholds were developed are not specific to humans. These fecal indicator bacteria can be shed from any warm-blooded organism, including wild and domesticated animals (12). Therefore, source tracking turns into a challenging scenario when these diffuse and frequently intermittent or episodic fecal releases occur. The second difficulty in identifying and eliminating sources of fecal indicator bacteria is their ubiquity in urban environments. Finally, unlike many human pathogens of concern, fecal indicator bacteria may survive and even grow in the environment (see, e.g., references 24, 39, and 44).
Viruses are one tool that could prove useful in source-tracking studies, because they include many pathogens of concern, and they are generally species specific. Viruses are known to cause a significant portion of waterborne disease, mostly from ingestion of sewage-contaminated water and seafood (10). Until recently, however, virus detection and quantification have relied on cell culture-based approaches that are much too slow to be effective source-tracking tools. Recently developed molecular techniques, such as quantitative reverse transcriptase PCR (QRT-PCR), can detect and quantify viral genetic material directly from water samples. Tests conducted previously in Southern California (11, 22, 32, 41, 42), Florida (17, 36), and Europe (35) using conventional RT-PCR or PCR have detected genetic material from human-specific viruses, including enterovirus, hepatitis A virus, rotavirus, and adenovirus, in urban runoff discharges or seawater samples.
A different approach would be to use alternative bacterial indicators for source tracking that might be much more abundant in urban discharges. For example, Bacteroides spp. make up approximately one-third of the human fecal microflora, considerably outnumbering fecal coliforms, E. coli, and Enterococcus spp. Bacteroides spp. are obligate anaerobes, so there is little concern over persistence or regrowth in the environment. More importantly, human-specific Bacteroides markers have been developed, increasing the value of this potential indicator (2, 3, 7).
Both viruses and alternative bacterial indicators such as Bacteroides spp. have been shown to be potentially useful source-tracking tools. Griffith et al. (18) and Noble et al. (33) concluded that genetics-based methods, such as PCR, consistently provided the best information in efforts to conduct source tracking on mixed-source samples. To date, however, no single method has all of the traits necessary to be the consummate source-tracking tool. Therefore, a multitiered, multiindicator approach has been recommended by some investigators (4, 40). By using multiple tools, investigators can utilize the strengths of each to ascertain inputs and track fates that will ultimately lead to successful management solutions.
The objective of this study was to identify the contributions and quantify the loading of fecal contamination affecting the SMB by using a multitiered approach. The first tier included traditional measurements of fecal indicator bacteria. The second tier included molecular assays developed and conducted for Enterococcus spp., a human-specific Bacteroides marker, and enterovirus. These methods rely on conventional PCR, quantitative PCR (QPCR), or QRT-PCR, which have not previously been applied in conjunction with one another for source-tracking studies in urban watersheds. The third tier involved sequencing of the enterovirus from the field samples with the greatest concentrations so as to determine the likely type of enterovirus amplified in the assay. The multitiered approach was applied using a mass-based design to quantify inputs and flux through an urban watershed to the beach.
|
|
|---|
![]() View larger version (28K): [in a new window] |
FIG. 1. Map of the Ballona Creek watershed in Los Angeles, Calif. Tributary and main-stem sampling sites for the water quality study are indicated. (Inset) Santa Monica Bay, in Southern California.
|
|
View this table: [in a new window] |
TABLE 1. Sampling sites along the main stem and major tributaries of Ballona Creek
|
After collection, samples were placed on ice and transported immediately to the University of Southern California for processing. For each composite sample, 100 ml of water was devoted to indicator bacteria analysis, and 200 to 600 ml of the sample volume was vacuum filtered through replicate 47-mm-diameter, 0.4-µm-pore-size polycarbonate filters (Poretics, Inc., Livermore, CA) using a filter funnel and receiver (Millipore, Inc., Bedford, MA) for Enterococcus sp. analyses by QPCR or Bacteroides sp. analysis by conventional PCR, as suggested by Haugland et al. (20). In addition, replicate filtrations were also conducted using 47-mm-diameter, 0.45-µm-pore-size type HA (Millipore, Inc., Bedford, MA) mixed cellulose ester filters for enterovirus analysis as suggested by Fuhrman et al. (11). For each filter, the volume filtered was the maximum filterable within ca. 10 min of the start of filtration. The total volume filtered was dependent on the location and turbidity of each individual sample, and the filter volumes were carefully recorded to the nearest 1 ml. All filters were placed in microcentrifuge tubes and stored at 80°C for further analysis.
Analysis of indicator bacteria.
Concentrations of E. coli and Enterococcus spp. were measured by defined substrate technology using kits supplied by IDEXX Laboratories, Inc. (Westbrook, ME) according to the manufacturer's instructions. Briefly, 10-fold and 100-fold dilutions of the water samples were made with deionized water containing the appropriate media and sodium thiosulfate, mixed to dissolve, dispensed into trays (Quanti-Tray/2000), and heat sealed. E. coli was measured using the Colilert-18 reagents, while Enterococcus spp. were measured using Enterolert reagents. Samples were incubated overnight according to the manufacturer's instructions and inspected for positive wells. Conversion of positive wells from these tests to a most probable number was done following Hurley and Roscoe (21).
Extraction of DNA.
The polycarbonate filters were processed for DNA extraction using the UltraClean fecal DNA isolation kit (MoBio Laboratories, Inc., Carlsbad, CA) according to the manufacturer's alternative protocol for maximum yield. Eluted DNA extracts were stored at 20°C until use.
The concentration of the extracted DNA was measured using the Quant-iT Picogreen double-stranded DNA reagent (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. Standard curves were generated in duplicate using Lambda DNA standards between 2.5 ng/ml and 600 ng/ml and a negative control (0 ng/ml). Fluorometric measurements were made using a Bio-Rad VersaFluor fluorometer.
Bacterial analyses using QPCR.
The total Enterococcus sp. primers and probe are described by Ludwig and Schleifer (28) and were constructed using the 23S rRNA gene regions around the target site of a well-established Enterococcus group-specific primer (ENC854R). Primer ECST748F targets Enterococcus spp., lactococci, and several clostridia. The target site of probe GPL813TQ is present in the 23S rRNA genes from a variety of representatives of gram-positive bacteria with low G+C DNA contents (28).
The master mix contained 1x Taq buffer, 4 mM MgCl2, 3 mM deoxynucleoside triphosphates (dNTPs), 2.5 U Ex Taq R polymerase (R PCR kit for quantitative PCR; TaKaRa Mirus Bio, Madison, WI), 1 µM ENC854R, 1 µM ECST748F, 0.1 µM GPL813 TQ Cy3 Probe (synthesized by MWG Biotech, High Point, NC), and nuclease-free water, yielding a final volume of 20 µl, to which 5 µl of sample (either DNA extract from an environmental sample, ranging from 1 to 76 ng genomic DNA, or 5 µl of lysed cell suspension or genomic equivalents) was added for a final volume of 25 µl. The samples were run under the following optimized assay conditions for PCR: 1 cycle consisting of an initial hold at 95°C for 2 min and 45 cycles of denaturation at 94°C for 15 s and annealing/extension at 60°C for 30 s (the optics were turned on during the annealing step). The Cepheid Smart Cycler II was set with the following specific parameters for this assay. The Dye Set was set for FCTC25. The cycle threshold analysis mode was set for growth curve (linear) analyses, with a manual threshold typically set at 5 to 15 fluorescence units. The background subtract level was set at a minimum of 12 and a maximum of 40. The BoxCar averaging feature was set at zero. The assay was previously optimized for Taq, Mg2+, and dNTP concentrations as well as all cycling parameters (data not shown). For quality control, Enterococcus faecalis ATCC 29212 and Enterococcus faecium ATCC 35667 combined were used as our calibration strains for the total-Enterococcus primer and probe set. Control bacterial preparations were prepared by boiling bacteria for 5 min, centrifuging for 2 min at 12,000 x g at 4°C, and storing immediately on ice. E. faecalis and E. faecium cells were enumerated using SYBR Green I epifluorescence microscopy (30). Serial dilutions of the standards were made in diethyl pyrocarbonate-treated sterile water, and four-point standard curves were run in duplicate in concert with the unknown samples on the Smart Cycler II instrument. Total-Enterococcus primers were tested with all 19 validly described species of the genus Enterococcus and demonstrated amplification of rRNA genes of all strains, with varying efficiencies (28).
Analysis of Bacteroides spp. using conventional PCR.
Amplification of the human-specific Bacteroides/Prevotella marker generally followed the procedure of Bernhard and Field (2), with PCR primers that amplify partial 16S rRNA from the human feces-specific group (BAC708R and HS183F). A range of extracted DNA quantities (representing 1 to 70 ng per assay, with most samples in the range of 5 to 20 ng) was tested to avoid problems with inhibition of the PCR. DNA was amplified with the Bacteroides-Provotella-specific primers described by Bernhard and Field (2). The PCR mixture was 1x Taq polymerase buffer, 1µm each primer, 200 µm dNTPs, 1.25 U Taq polymerase, 640 ng µl1 BSA, and 1.5 mM MgCl2. The PCR conditions were specifically optimized for this study and differed from those in the original publication: 2 min at 95°C, then 25 cycles of 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s, followed by a 5-min extension at 72°C. Then 1 µl of each PCR product was reamplified using the conditions given above for another 25 cycles. PCR was performed on a 3000MX thermal cycler (Stratagene). PCR products were visualized in a 2% agarose gel stained with 1x SYBR Gold (Molecular Probes, Eugene, OR) and compared to a 100-bp DNA ladder (Promega). Positive results yielded 525-bp amplicons. The positive control was a human fecal sample extracted with a QIAamp DNA stool kit (QIAGEN, Valencia, CA). Negative controls contained water instead of sample. All samples were initially run with 5 µl of extracted material. After the initial analyses, all negative samples (14 of 60) were spiked with 0.1 ng of positive-control DNA and reanalyzed to determine possible inhibition. Three of the 14 negative samples (two from Benedict Box Channel and one from Adams) were determined to be inhibited. Inhibited samples were reanalyzed using 2 µl of sample, but all three remained negative.
Enterovirus analyses.
Samples, along with negative controls, were extracted using a modified RNeasy plant minikit (QIAGEN Plant and Fungi RNA isolation protocol). The RLT homogenization buffer supplied was supplemented with polyvinyl pyrrolidone 40 (PVP-40) at a final concentration of 2%. Prior to the extraction, fresh ß-mercaptoethanol (Sigma Chemical Co.) was added to the extraction buffer in the exact concentration recommended by the manufacturer. Filters (type HA) were manually homogenized with a pipette tip, and 700 µl of the RLT/filter slurry was applied to a QiaShredder column (QIAGEN) and spun at maximum speed,
8,000 x g, for 2 min to aid in viral lysis as well as to separate filter particles from the filtrate (11). The filtrate was then carefully removed without disturbing the pelleted material and was placed into a new 1.5-ml tube. The volume of solution in each tube was estimated by pipetting, and 0.4 volume of 5 M potassium acetate (pH 6.5) was added. Tubes were mixed by inversion and incubated on ice for 15 min. The mixture was then spun (12,000 x g) at 4°C for 15 to 30 min and the supernatant transferred to a new 1.5-ml microcentrifuge tube. Subsequently, the Plant and Fungi RNA isolation protocol was followed starting at step 5. RNA was eluted with 50 µl of the supplied RNase-free water. A one-step TaqMan QRT-PCR was performed on the extracted RNA, with final reaction volumes of 25 µl, using a QIAGEN One-Step RT-PCR kit. Five microliters of extracted RNA was added to 20 µl of master mix containing 1x RT buffer, 6 mM MgCl2, 500 nM dNTPs, 700 nM EV1 reverse primer (5'-TGTCACCATA AGCAGCCA-3'), 700 nM EV1 forward primer (5'-CCCTGAATGCGGCTAAT-3'), 30 µg bovine serum albumin, 20 U of recombinant RNasin (Promega Corp.), 1.5% PVP-25 (29) (Sigma Chemical Co.), 100 genomic equivalent units of a competitive internal positive control (CIPC) developed in-house (16) by following the general approach of Kleiboeker (25), plus 300 nM CIPC probe 5'-Cy5-TGTGCTGCAAGGCGATTAAGTTGGGT-BHQ-2-3', and 300 nM EV-BHQ probe 5'-6-carboxyfluorescein-ACGGACACCCAAAGTAGTCGGTTC-BHQ-1-3', and 1 µl of enzyme mix (containing both reverse transcriptase and DNA polymerase). The probe and primers were synthesized by MWG Biotech, Inc. The Cepheid Smart Cycler II was programmed as follows: a 1-h reverse transcription step at 50°C followed by a 15-min hold at 95°C for DNA polymerase activation, then 45 cycles of 94°C for 15 s (denaturation) and 60°C for 1 min (annealing and extension, with optics on).
QRT-PCR results were available 3 h after the start of analysis, making the total RNA extraction, QRT-PCR preparation, and analysis time less than 5 h. After the first analyses, those samples that appeared to have inhibition of the QRT-PCR (as indicated by the CIPC) had RNA sample volumes reduced in half and were rerun. No inhibition was observed at this lower RNA concentration. Standard curves were generated using a synthetic enterovirus transcript that was quantified using fluorometric analysis, and sample genome concentrations were interpolated from the standard curve using the manufacturer's curve-fitting software. Quantitative results are reported per liter sample volume.
Sequencing of enterovirus QRT-PCR positives.
QRT-PCR-positive samples from the enterovirus analyses were sequenced to assess the specificity and fidelity of our enterovirus QRT-PCR and to elucidate the identities of the enterovirus genomes being amplified by the assay. Following the initial enteroviral analysis of the Ballona Creek samples, RNA samples identified as positive for enterovirus genomes were again amplified using the QRT-PCR protocol. However, after numerous repetitions of freezing and thawing of the extracted RNAs for various analyses, only Cochran Ave. samples from 9:00, 10:00, and 11:00 contained amplifiable enterovirus genetic material (11). The 144-bp enterovirus QRT-PCR product was distinguished from the 126-bp internal positive control product using a 10% polyacrylamide-1x Tris-borate-EDTA gel. The gel was stained with ethidium bromide and visualized with a 254-nm UV light. The larger 144-bp enterovirus QRT-PCR products were excised, homogenized in 50 µl of 1x QIAGEN One-Step RT-PCR buffer using a microcentrifuge pestle, and incubated overnight at 37°C with shaking. The enterovirus products were purified (Wizard SV gel and PCR Clean Up System; Promega), cloned into the 3.9-kilobase pCR 2.1 TOPOVector (TOPO TA cloning kit; Invitrogen, Carlsbad, CA), transformed into TOP10 chemically competent E. coli, and plated on Luria-Bertani agar plates containing 100 µg/ml of ampicillin. Bacterial clones were screened using QPCR, and positive colonies from each of the three sites were selected and grown individually in 3 ml of 2x medium overnight at 37°C. Plasmid DNA was isolated (PerfectPrep plasmid minikit; Eppendorf, Westbury, CT), and the DNA concentration was calculated using yeast extract-tryptone Ribogreen (Molecular Probes). Plasmid DNA was sequenced bidirectionally by MWG Biotech, and the chromatograms were inspected using Sequencher, version 4.2 (Gene Codes Corp., Ann Arbor, MI). The sequences were aligned to sequences in the NCBI GenBank using BLAST and are available by searching nucleotide accession numbers DQ196482 through DQ196487.
Calculations and statistical analyses.
Data analysis comprised four steps. First, the hydrologic budget was evaluated to determine if the majority of the flow was sampled. This evaluation was conducted by comparing the volumetric inputs from each of the tributaries to the volumetric discharges along the main stem of Ballona Creek. The second step was to examine temporal and spatial trends in the flux of fecal indicator bacteria. The flux of indicator bacterial cells per hour was calculated by multiplying the concentration of the indicator bacteria (per deciliter [100 ml]) by the flow rate (in deciliters per hour). The mean hourly flux (temporal analysis) was calculated by averaging the flux of indicator bacteria at all main-stem locations for each hourly interval (n = 6). The mean flux at each site (spatial analysis) was calculated by averaging the flux at each main-stem or tributary site for all hourly samples (n = 6). Statistical analysis of the differences in bacterial flux between hourly time periods or, alternatively, between main-stem sites was conducted using analysis of variance (46). The third data analysis step was to examine spatial and temporal patterns in the frequency of Bacteroides detection. The Bacteroides method used in this study was a presence/absence end point. This examination was conducted by tabulating the locations and time periods of Bacteroides detection to detect patterns moving downstream, adjacent to tributaries, or over time. The fourth data analysis step was to examine the spatial and temporal extent of enterovirus concentrations. Unlike that of Bacteroides, the presence of enterovirus was quantified, so the magnitude of enterovirus concentrations was tabulated among the different locations across the different time periods. As in the Bacteroides data analysis, patterns moving downstream, adjacent to tributaries, or over time were examined.
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 2. Schematic diagram depicting additive flow in the main channel of Ballona Creek and the percentage contributed by each tributary sampled.
|
![]() View larger version (8K): [in a new window] |
FIG. 3. Mean flux of E. coli and Enterococcus cells (expressed as cells per hour) at main-channel sampling sites of Ballona Creek (26 August 2004).
|
![]() View larger version (9K): [in a new window] |
FIG. 4. Mean hourly flux of Enterococcus spp. (expressed as cells per hour) along the main channel of Ballona Creek, measured by using either an IDEXX chromogenic substrate (Enterolert) or QPCR methods on 26 August 2004.
|
![]() View larger version (28K): [in a new window] |
FIG. 5. Loading of Enterococcus spp. (expressed as cells per hour) in the main channel and tributaries of Ballona Creek at 9:00 (a), 10:00 (b), 11:00 (c), 12:00 (d), 13:00 (e), and 14:00 (f) on 26 August 2004.
|
The human-specific Bacteroides marker was positively detected throughout the main stem of Ballona Creek; its presence in 31 of the 36 main-stem samples tested (86%) was confirmed (Table 2). A positive signal was observed at all sites during the early morning hours, but during periods of heightened UV radiation (midday and late afternoon), a decrease in the number of positive results for the human-specific Bacteroides marker was observed. For example, only three of six samples at 14:00 were positive for the human-specific Bacteroides marker and exceeded the water quality threshold for E. coli.
|
View this table: [in a new window] |
TABLE 2. Number of enterovirus genomes per liter along the main stem of Ballona Creek in Los Angeles, CA
|
The highest sequence homology observed for the three enterovirus-positive Cochran Ave. samples that were sequenced (>95%) was with human coxsackievirus A22 (GenBank accession no. AF499643), human coxsackievirus A19 (accession no. AF499641), and human enterovirus 90 (accession no. AY773285). The 144-bp QRT-PCR product corresponded to nucleotides 453 to 596 of human coxsackievirus A22, nucleotides 457 to 600 of human enterovirus 90, and nucleotides 454 to 597 of human coxsackievirus A19.
|
|
|---|
Our study is not the first to examine the presence of viruses in urban runoff entering shorelines in the SMB and other Southern California urban watersheds. For example, Gold and colleagues (13, 14) found viruses in repeated samples from multiple storm drains to the SMB by using both cell culture and RT-PCR techniques. Haile et al. (19) detected human-specific viruses in all three storm drains tested in their epidemiological study of the SMB. Noble and Fuhrman (32) found human enteric virus genomes in the near-shore marine waters of the SMB. Jiang et al. (22) found human adenovirus in samples collected at 12 sites between Malibu and the Mexican border, and Fuhrman et al. (11) previously found human enterovirus genomes in Ballona Creek.
The multitiered approach used in this study can assist watershed managers in determining sources and efficiently abating the most significant inputs of fecal indicator bacteria. If managers relied solely on the patterns in fecal indicator bacteria from Ballona Creek, then the only option would be to treat the entire 37-m3 s1 discharge furthest downstream at Inglewood Ave., because the flux of fecal indicator bacteria was similar at all sources. The use of multiple tools, however, allows managers to prioritize the most important sources. In this case, the presence of human enterovirus was greatest at the Cochran Ave. site, where the system daylights from the underground storm drain system beneath Los Angeles and the discharge volume is one-third of the volume at Inglewood Ave. Previous studies of Southern California storm drains have detected a human-pathogenic virus signal (22, 23, 31, 32). Since Cochran Ave. had the most frequent occurrence and highest concentrations of enterovirus, plus a consistent co-occurrence of the human-specific Bacteroides marker, this source would appear to be the most likely candidate for future management actions. The sequencing results that confirmed the presence of several potential risks to human health (human coxsackievirus and enterovirus) should provide the reassurance most managers would need before planning future management steps.
The lack of correlation between bacterial indicator levels and levels of human-pathogenic viruses has been observed in previous studies (8, 9) and demonstrates the value of the multitiered approach used here for source identification. For example, analysis of wild shellfish from the Atlantic coast of France indicated no significant correlation between fecal coliforms and enteroviruses or hepatitis A virus (26, 27), and viruses have sometimes been found in oysters without coliform contamination (15, 45). Noble and Fuhrman (32) detected enterovirus in 35% of the 50 shoreline samples they examined over a 5-year period, and no significant statistical relationship to any of the standard bacterial indicators was found. Virus and fecal indicator bacteria were measured in dry-weather urban runoff in drains along 300 km of shoreline from Santa Barbara to San Diego, CA (31). Although 40% of the storm drains contained detectable enterovirus, there was no correlation with concentrations of fecal indicator bacteria. It is also possible that differential rates of degradation of viruses and bacteria can explain much of the discordant relationship between viral pathogens and indicator bacteria (see, e.g., reference 34).
The use of QPCR to measure fecal indicator bacteria presents unique opportunities and challenges. An advantage of using QPCR for measuring fecal indicator bacteria is speed; the method potentially provides measurements in less than 3 h (18). However, culture-based methods quantify only viable bacteria, while QPCR measures the DNA from both cultivable and noncultivable microbes. This was most apparent in the temporal trends from Ballona Creek. Levels of Enterococcus spp. determined using culture-based methods generally decreased as the day progressed, most likely as a result of photoinactivation of cells by sunlight (5, 6, 34). Ballona Creek is a 40-m-wide concrete-lined channel, concentrating solar energy into the shallow creek in the channel invert. The QPCR results, however, remained steady, indicating that the bacterial DNA was still intact and detectable, even though the Enterococcus spp. were not viable.
Overall, the use of multiple approaches provided convincing evidence of the extent and types of microbial contamination in this urban watershed. We believe such studies should provide invaluable information for researchers and managers trying to balance regulatory burdens and public safety.
The contents of this document do not necessarily reflect the views and policies of the SWRCB, nor does mention of trade names or commercial products constitute endorsement or recommendation for use.
We thank Orange County Public Health Laboratories (Douglas Moore) for their high-titer poliovirus stock and Ian Hewson, Joshua Steele, Mike Schwalbach, and Sheila O'Brien for assistance with sample collection and processing. We also thank the volunteers of Santa Monica Baykeeper for assistance with sample collection and environmental parameter data collection.
|
|
|---|
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»