Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2005, p. 3179-3183, Vol. 71, No. 6
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.6.3179-3183.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Microbiology, Oregon State University, Corvallis, Oregon 97331
Received 16 July 2004/ Accepted 31 December 2004
|
|
|---|
|
|
|---|
This study addresses two specific fecal pollution issues requiring source identification. The Tillamook Bay watershed on the Oregon coast is home to a large dairy industry, but it also contains a growing population of wildlife, including elk herds. Fecal pollution has had an impact on the economy of the local shellfish industry. Appropriate assessment and abatement strategies and the cooperation of local farmers depend on the ability to distinguish host sources of fecal pollution. A second issue involves runoff from urban storm drains that carry rainwater into rivers and streams. This may include fecal pollution from leaking septic systems or the excrement of domestic pets. Efficient source identification will help to evaluate human health risk and expedite water quality management procedures.
Fecal members of the order Bacteroidales are abundant in the feces of warm-blooded animals (9, 15, 18, 19, 25), and some have host species or group-specific distributions (2, 18). Using terminal restriction fragment length polymorphism and clone library analysis, Bernhard and Field (3) developed PCR primers that distinguish between human and ruminant fecal pollution based on differences in Bacteroidales partial 16S rRNA gene sequences. The primers have been used successfully in trials and field studies in both freshwater and saltwater environments (4, 5, 12, 13). Similar analyses were used to develop markers for pig and horse fecal pollution (9a). However, the clone library analysis involving eight host sources showed that humans, cats, and dogs share very closely related Bacteroidales sequences, as do cows and elk. We were unsuccessful in developing markers to differentiate cow and elk feces based on sequence data from clone libraries. In addition, although two previously developed human-specific Bacteroidales PCR primers did not amplify dog fecal DNAs (3), clone library sequences did not support design of a dog-specific primer.
Determining microbial diversity in a complex community requires an extensive analysis of large clone libraries. In a typical library of 100 to 300 clones, often only one or two sequences are identical at the species level, indicating a significant lack of coverage (11, 16). Terminal restriction fragment length polymorphism analysis to identify unique restriction patterns can reduce the number of clones sequenced, but pattern comparison involves an additional element of subjectivity that also results in error and incomplete sequence information. These gaps in sequence data are especially detrimental when the goal is to identify unique markers among very closely related sequences. The aim of this study was to identify host-specific Bacteroidales rRNA gene markers by comparing the genes empirically rather than by relying on sequence data alone. We used a technique based on subtractive hybridization, whereby genetic differences between closely related genomes are amplified to generate unique fragments. Subtractive hybridization has been used to compare bacterial genomes for identification of virulence factors (7, 8, 24), to define regions present in a sequenced genome but absent in an unsequenced relative (1), and more recently to design primers that identify individual members of microbial communities (22). Based on previously published methods (22, 26, 27), this report describes a modified subtractive hybridization for selective enrichment of Bacteroidales rRNA gene sequences. We compared Bacteroidales rRNA genes from a source of interest for primer design (target) with those of one or more reference sources (subtracters) in a solution hybridization according to the procedure described by Zwirglmaier and colleagues (27). An experimental overview is given in Fig. 1. Each hybridization took place in a single microplate well. We used unique unhybridized sequences to design primers for fecal source identification.
![]() View larger version (34K): [in a new window] |
FIG. 1. Overview of subtractive hybridization in microplate wells, adapted from the technique of Zwirglmaier et al. (27).
|
|
|
|---|
Sample collection and DNA preparation.
Fecal samples were collected by hunters, farmers, and colleagues and were also acquired from animal shelters. Samples were stored in guanidine isothiocyanate buffer (5 M guanidine isothiocyanate, 100 mM EDTA, pH 8, 0.5% Sarkosyl) at 80°C. The FastDNA kit for soils (Q-Biogene, Carlsbad, CA) was used for DNA isolation. DNAs from 10 to 30 individual fecal samples were mixed in approximately equal concentrations (3 ng/µl) to create a pool of genomic DNA for each host species. All DNA quantifications were done by a PicoGreen assay (Molecular Probes, Inc., Eugene, OR).
PCR amplification and restriction digestion.
Target and subtracter genomic DNA pools were amplified using PCR primers extended with AciI restriction sites. The Bacteroidales-specific 16S rRNA primer AciBac32F (3) and the universal 23S rRNA primer Aci422R (21) amplified approximately 2,400 bp, including most of the Bacteroidales 16S rRNA gene, the intergenic spacer region (internal transcribed spacer), and a portion of the 23S rRNA gene. The oligonucleotides used as primers and linkers in this study are listed in Table 1. Each 50-µl PCR mixture contained 1x Taq polymerase buffer, each primer at a concentration of 10 µM, each deoxynucleoside triphosphate at a concentration of 200 µM, 0.06% bovine serum albumin, 2 mM MgCl2, and 1.25 U of the proofreading Taq polymerase TaKaRa Ex Taq (Takara Bio, Inc., Shiga, Japan). Cycling parameters were as follows: 30 cycles of 94°C for 1 min, 68°C for 1 min, and 72°C for 2.5 min, followed by a final extension at 72°C for 7 min. PCR products were purified with a GeneClean II kit (Q-Biogene) and digested with AciI (New England Biolabs, Inc., Beverly, MA).
|
View this table: [in a new window] |
TABLE 1. Primers and linkers used in this study
|
Immobilization of subtracter in microplate well.
Equal amounts of the subtracter DNA from the two host sources were mixed (1 µg total) and diluted in 50 µl phosphate-buffered saline buffer (8.0 mM Na2HPO4, 1.5 mM KH2PO4, 137 mM NaCl, 27 mM KCl, [pH 7.2]; and MgCl2 added to a final concentration of 100 mM). The mixture was heated to denature the DNAs and immediately placed into a MaxiSorp microplate well (Nalge Nunc, Naperville, IL). The plate was incubated at 37°C for 1 h, the buffer was removed, and the plate was dried for 2 h at 70°C.
Solution hybridization.
Target DNA (20 ng) was diluted in 40 µl hybridization buffer (2.5x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate], 1% blocking reagent [Roche Applied Science, Indianapolis, IN], 0.01% sodium dodecyl sulfate, 0.05% N-lauroylsarcosine, and formamide). A gradient of formamide concentrations (29%, 35%, and 41%) provided hybridizations in three wells at a range of effective temperatures for stringency optimization. Target DNA was heat denatured, iced, and added to the microplate wells. Hybridization was carried out at 70°C for 2 h.
Amplification of subtraction products, cloning, sequencing, and primer design.
Two microliters of subtracted, unhybridized target DNA was removed from the supernatant and amplified with primer T1 using the cycling parameters described above for this primer. PCR fragments were gel extracted (QiaQuick gel purification kit; QIAGEN) following dilution and reamplification and cloned into TOPO TA vectors (Invitrogen Corp.). Ten individual clones were randomly selected for bidirectional sequencing on an ABI 3100 capillary sequencer. Sequences were aligned using ARB (20), either with related GenBank sequences (NCBI BLAST) or with our clone library sequences if they matched the rRNA region for which we had sequence data (16S rRNA; E. coli positions 32 to 708). The short, subtracted (SH) target sequences were added to a neighbor-joining tree of full-length sequences using the parsimony insertion tool of ARB. New primers DF475F and EF990R were designed using the Probe Design and Probe Match functions of ARB.
Primer specificity and sensitivity.
Primers were tested for cross-reactivity against host pools of fecal DNAs representing 4 to 30 individual hosts from each species. Primer specificity was optimized by manipulation of annealing temperature, MgCl2 concentration, and cycle number on a PCR Express thermocycler (Thermo Hybaid, Middlesex, United Kingdom). Primer sensitivity was estimated in serial dilutions of plasmid-inserted templates of known copy numbers. Theoretical detection limits were determined in pure water, creek water (from Beaver Creek, Alsea, OR), and seawater (from 5 miles off the central Oregon coast). Each reaction mixture contained 3 ng total DNA from the natural water sources.
Nucleotide sequence accession numbers.
The partial 16S rRNA gene sequences have been deposited in the GenBank database under accession numbers AY70144475 through AY70144492.
|
|
|---|
The restriction digests of target and subtracter DNAs resulted in 100- to 500-bp fragments, visualized by gel electrophoresis and shown for the elk experiment in Fig. 2A. Subtracted target fragments of 300 to 500 bp were recovered in both the elk and dog experiments following hybridization and PCR amplification with primer T1. A gel image of the subtracted fragments from the elk experiment is shown in Fig. 2B. Band intensity appeared to increase with increasing stringency (formamide concentration) in fragments from the elk wells versus those from the human and cow wells but not in the human subtracter-only controls. However, since the PCR may not have been quantitative, no conclusions could be drawn from this observation. Subtracted target fragments from the 35% formamide wells were cloned and sequenced. Products in the control wells were not sequenced.
![]() View larger version (46K): [in a new window] |
FIG. 2. (A) AciI restriction digest of target (elk fecal) and subtracter (cow and human [hum] fecal) rRNA genes following PCR amplification with Bacteroidales-specific primers. (B) Elk fecal subtraction products amplified with target-specific linker T1 as PCR primer. HF, human fecal; CF, cow fecal.
|
![]() View larger version (28K): [in a new window] |
FIG. 3. Relationships of subtracted target fragments to database or host-specific clone library sequences. (A) Target elk feces-derived sequences clustered separately from related database sequences. (B) Target dog feces-derived sequences clustered with dog and cat fecal Bacteroidales sequences from an earlier clone library analysis (L. K. Dick et al., unpublished data).
|
![]() View larger version (40K): [in a new window] |
FIG. 4. Host DNA pools PCR amplified with new source-specific primers derived from products of subtractive hybridization. (A) Primer EF990R distinguishes elk fecal DNA from cow and human reference sources, as well as from pig and deer sources. (B) Primer EF990R distinguishes elk from cow feces even when total genomic DNA in the cow fecal pool is in fourfold excess to total DNA in elk fecal pool. (C) Primer DF475F distinguishes dog fecal DNA from human and cat fecal DNAs as well as from other host sources.
|
Primer DF475F was used successfully in a national comparative study of fecal-source tracking methods (12, 14). Unidentified aqueous samples containing mixtures of four fecal sources or sewage were prepared in three different matrices: distilled water, salt water, and humic acids. DF475F correctly identified the marker in all samples containing dog feces, including a sample containing humic acids.
|
|
|---|
Maxisorp microplates (Nalge Nunc) contain a polystyrene surface that binds DNA noncovalently. Since immobilized DNA is stable for at least 6 to 12 h at 68°C (27), we expected that the 70°C, 2-h hybridization would occur without significant loss of subtracter DNA to the solution. A large excess of subtracter over target DNA was used to allow for any loss of subtracter due to reannealing.
Three formamide concentrations were chosen based on Na2+ concentration, GC content, and fragment sizes. Formamide is a denaturant and has the effect of decreasing the melting temperature of the DNA hybrid (23). Overly stringent conditions result in hybridization of perfect matches only, and subtraction products could lack the variability necessary for primer design. If too relaxed, the conditions result in no product at all. Since all three formamide concentrations resulted in similar bands, products of the 35% concentration were chosen for sequencing because this concentration produced a sufficient amount of unhybridized DNA for cloning. Ultimately, the design of primers that distinguish target from subtracter host Bacteroidales indicated that an appropriate stringency range was attained. Future experiments might include a greater range of stringencies or sequencing of products from more than one stringency for comparison, but because our primer design objectives were met, we did not pursue this further.
No fragments of less than 300 bp were recovered in either of the experiments, for which there are several possible explanations. Some of the restriction fragments ligated to each other, forming chimeric sequences; some small fragments may have been lost in this way. PCR bias may have resulted in preferential amplification of the most-abundant larger restriction fragments. Alternatively, the inverted terminal repeats created by ligation of identical linkers to either end of the fragments may have resulted in the formation of panhandle structures that prevented primer binding and amplification (1). This would be more kinetically favorable with smaller fragment sizes. Inouye and Hondo found that hybridization efficiency declined as fragment size fell below 227 bp (17). These observations suggest the use of a restriction enzyme that produces fragments larger than 200 bp.
The intent in using the entire 16S rRNA gene, internal transcribed spacer, and part of the 23S rRNA gene was to obtain as much sequence information as possible for primer design. Three of 20 total subtracted fragments aligned with the 23S rRNA gene, but the small number of database sequences available for comparison made it difficult to use these fragments for primer design. Fourteen of the 20 subtraction products were from the 16S rRNA gene, and all 10 of the elk feces-specific clones were from a region known to be hypervariable (E. coli positions 940 to 1370) (6).
Primer DF475F was the first and only dog-specific primer designed and tested, and it did not require stringent optimization to exclude amplification of other host fecal DNAs. This illustrates the ability of this method to enrich for regions of variability. Primer EF990R was the third elk-specific primer tested and the first to successfully distinguish cow and elk fecal Bacteroidales sequences. It did not distinguish elk and sheep Bacteroidales sequences, which was not the intent of this experiment, but that may be accomplished with another subtractive hybridization. Alternatively, it may be possible to design a primer based on the 5' and 3' ends of a target fragment. Any cross-reacting DNA could be sequenced, allowing more directed primer design.
Primer EF990R amplified elk fecal DNAs and did not amplify cow fecal DNAs, even when four times more cow fecal DNA was added. A previous study using blind samples (12, 14) demonstrated that Bacteroidales source-specific primers identified fecal sources correctly when the sources comprised as little as 1% of the total fecal contamination in the sample.
Previous studies reported the immobilization of up to seven genomic subtracters in one well (26, 27). An initial attempt to use seven subtracters with an elk-specific target resulted in subtraction products without enough variability for primer design. Future experiments may determine what limitations exist and what adjustments can be made in the use of multiple subtracters.
Microplate subtractive hybridization was successfully employed to generate a unique source-specific marker for dog fecal pollution and a marker that differentiates elk and cow fecal pollution. The results demonstrate the method's ability to enrich for variable regions of the 16S rRNA gene. In addition, the method could readily be adapted for other gene targets. The capacity to characterize markers that distinguish sources of fecal pollution without obtaining large numbers of clones for each new host will expedite the addition of new source markers for fecal pollution.
We thank Katrin Zwirglmaier for her correspondence and advice and Sarah Walters for helpful discussions and support.
|
|
|---|
e, and S. B. Grant. 2003. Tiered approach for identification of a human fecal pollution source at a recreational beach: case study at Avalon Bay, Catalina Island, California. Environ. Sci. Technol. 37:673-680.[Medline]
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»