Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2002, p. 6077-6086, Vol. 68, No. 12
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.12.6077-6086.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Colleen M. Cavanaugh,2 and Martin F. Polz1*
Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139,1 Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 021382
Received 1 April 2002/ Accepted 29 August 2002
|
|
|---|
|
|
|---|
Despite this extensive sequencing effort, the diversity of microbial communities is still not sufficiently explored, a fact that is exacerbated by potentially significant variation in community structure even between physically and chemically similar environments (21). Thus, a major limiting step in the application of the rRNA framework remains the potentially large number of sequences to be determined and probes to be synthesized for a satisfactory correlation of specific environmental factors to the presence of specific populations. In this context, there is a need for an efficient means to generate large sets of oligonucleotides with high specificity.
Here, we present a novel method for the generation of oligonucleotide probes which give good phylogenetic resolution without prior knowledge of the target sequence. The method relies on two consecutive PCR amplifications with Bacteria-specific primers bracketing a highly variable region in the 16S rDNA, each followed by detachment of one primer by alkaline hydrolysis of a 3' ribonucleotide. Ultimately, a single-stranded oligonucleotide of the variable region is obtained. Comparative analysis of sequences obtained from databases and empirical testing of representative isolates show good differentiation of targets by these PCR-generated probes. Because of the low cost and potential for production of the probes in microtiter plate format, anticipated applications include the rapid screening of clone libraries or of strain collections for potential abundance in specific environmental samples.
|
|
|---|
-, ß-, and
-Proteobacteria (7) (Fig. 1). Synechococcus and Prochlorcoccus were grown in nutrient-enriched Sargasso seawater as previously described (12), whereas all other isolates were grown in Luria broth at 20°C.
![]() View larger version (54K): [in a new window] |
FIG. 1. Neighbor-joining tree of 16S rRNAs for a collection of 44 bacterial isolates from a coastal marine environment, reference strains, and sequences retrieved from GenBank. The scale bar indicates 10% sequence difference. Organism names for the isolates are from the closest matches in GenBank. Group affiliations, GenBank accession numbers, and isolate-specific target sequences for PCR-generated 16S rRNA oligonucleotide probes are shown for each sequence.
|
Nucleic acid extractions.
DNA for PCR was extracted by bead beating of cell pellets in 2-ml tubes with 0.5 g of zirconium beads (0.2-mm diameter) in 650 µl of TE (10 mM Tris, 1 mM EDTA, pH 8) buffer for 30 s at 5,000 rpm. Lysates were incubated for 30 min at 37°C with 0.2 mg of proteinase K ml-1 and 0.06% sodium dodecyl sulfate (SDS), followed by inactivation of the enzyme at 65°C and addition of NaCl to a final concentration of 1 M. Nucleic acids were purified by standard phenol-chloroform extraction and by passing the final aqueous solution over Strataprep PCR purification columns (Stratagene, La Jolla, Calif.).
Total nucleic acids for membrane hybridizations were extracted by bead beating as described above except in the presence of 600 µl of RNase-free lysis buffer (5.25 M guanidine-HCl, 1.5% Triton X-100, 2.5 mM Tris, 0.25 mM EDTA, pH 7.5). Cell debris and beads were sedimented at 12,000 x g for 1 min, and 400 µl of the supernatant was transferred to a new tube containing 200 µl of isopropanol. Total nucleic acids were purified by using SNAP total RNA spin columns (Invitrogen, Carlsbad, Calif.).
Primer design.
The Bacteria-specific primers 64f (Table 1) and 104r (Table 2) were designed based on a complete alignment of the Ribosomal Database Project (RDP) in 1996 by using the following procedure. First, the alignment of representative 16S rRNA sequences published by the RDP was inspected for two stretches of high sequence conservation bracketing a region of high sequence variation (probe region) (Fig. 2). Second, a forward and a reverse PCR primer, each containing a single degeneracy, were designed matching the highly conserved regions. Third, additional degeneracies were added to the primer sequences based on comparison of the primers with the entire RDP alignment of several thousand 16S rRNA sequences to maximize matches of the primers against bacterial sequences. The resultant primers, 64f and 104r, were then tested for amplification success and properties against the 44 bacterial isolates and standard isolates listed in Fig. 1. In addition, the universal primer 519r (8) (Table 3) used in the first step of the probe synthesis was tested against the same alignment to test how well it matches bacterial sequences. Finally, all three primers were compared to the most recent alignment of 16,277 sequences in the RDP II (10) by using the Rose software (2). Archaeal sequences and database entries lacking sequence information in the target regions were manually excluded from the analysis.
|
View this table: [in a new window] |
TABLE 1. Design of primer 64f and similarity to 16S rRNA sequences in the RDP II databasea
|
|
View this table: [in a new window] |
TABLE 2. Design of primer 104r and similarity to 16S rRNA sequences in the RDP II databasea
|
![]() View larger version (23K): [in a new window] |
FIG. 2. Partial secondary structure of 16S rRNA from E. coli (4), with target regions for the three Bacteria universal primers (64f, 104r, and 519r) used in the probe synthesis indicated by solid lines. The dashed line indicates the variable probing region.
|
|
View this table: [in a new window] |
TABLE 3. Primer 519r similarity to 16S rRNA sequences in the RDP II database
|
PCR generation of oligonucleotide probes.
The oligonucleotide probes were produced by using the following five-step procedure outlined in Fig. 3, including two PCR amplifications with detachable primers.
![]() View larger version (34K): [in a new window] |
FIG. 3. Outline of PCR-based method to generate population-specific oligonucleotide probes without prior sequence information.
|
(ii) Binding of PCR product (step B).
Biotinylated products were bound to 0.1 mg of paramagnetic and streptavidin-coated Dynabeads M-280 (Dynal, Skøyen, Norway) according to the instructions supplied by the manufacturer after excess primer and nucleotides were removed by using the QIAquick PCR purification kit (Qiagen, Valencia, Calif.). To remove the nonbiotinylated complementary strand, beads were washed twice at 85°C for 1 min in autoclaved Q-grade water (Millipore), using an MPC magnet (Dynal) to separate beads from liquid.
(iii) Primer detachment by alkaline hydrolysis of ribonucleotide (step C).
The beads with the bound PCR product were washed once at room temperature in sterile Q-grade water before the ribonucleotide was hydrolyzed by incubation of the beads for 12 h in 100 µl of 0.2 M NaOH at 37°C with gentle mixing. After hydrolysis, primers bound to beads were removed from solution with the MPC magnet. The aqueous solution containing the single-stranded PCR product was diluted to 400 µl in 10 mM Tris (pH 8) and precipitated with sodium acetate and ethanol (15). The precipitated DNA was dissolved in 50 µl of Tris (10 mM, pH 8) and stored at -20°C.
(iv) Linear amplification (step D).
The single-stranded fragment from step C served as the template in a linear amplification of the population-specific probing region with primer 104r modified with a 5' biotin and 3' ribonucleotide. Approximately 0.5 pmol of single-stranded template was added to each tube. Amplification conditions and reagents were as described above.
(v) Recovery of oligonucleotide probe (step E).
The reaction product from the linear amplification was bound to 0.2 mg of magnetic beads, which were washed twice at 85°C in Q-grade water to remove the complementary strand. The amplified product was separated from biotinylated primers by alkaline hydrolysis and removal of the magnetic beads as described above. The final product, which is complementary to the probing region of the 16S rRNA (Fig. 2), was precipitated, dried, and dissolved in 10 mM Tris (pH 8) for long-term storage at -20°C.
Hydrolysis efficiency.
To test the efficiency of the primer detachment by alkaline hydrolysis, two permutations of the 104r primer, one with and one without the 3' ribonucleotide (Synthetic Genetics, San Diego, Calif.), were 5' end labeled with 32P by using T4 polynucleotide kinase (New England Biolabs, Beverly, Mass.) (15). Primers were separated from unreacted nucleotides in Sephadex G-50 spin columns (Amersham Pharmacia Biotech, Piscataway, N.J.), and both were used in separate PCRs with nonmodified primer 64f.
PCRs were performed in 40-µl volumes containing 100 ng of genomic DNA from E. coli K-12 by using primer pair 64f-104r or 64f-104r-3' ribonucleotide and the same conditions and reagents as described above. PCR products were agarose gel purified from unreacted nucleotides and primers by using QIAquick gel extraction kits (Qiagen). Half of the volume of each purified product (64f-104r and 64f-104r-3' ribonucleotide) was subjected to alkaline hydrolysis for 12 h (0.2 M NaOH, 37°C), whereas the second half was kept frozen in 10 mM Tris (pH 7.5) as a reference. Both hydrolyzed and reference samples were precipitated with ammonium acetate and ethanol (16), resuspended in 20 µl of Q-grade water, and separated in a 10% denaturing acrylamide gel on a Mini-Protean 3 system (Bio-Rad, Hercules, Calif.) at 9 V cm-1. Gels were vacuum dried on blot paper, and bands were identified by 32P-phosphorimaging with multipurpose plates and a Cyclone phosphorimager (Packard Instruments).
Linear amplification optimization.
The yield of the linear amplification for the probe synthesis (Fig. 3D) was determined with PCR mixtures containing 32P-labeled primer 104r. Single-stranded PCR product was generated from 40 ng of Shewanella baltica genomic DNA as outlined in Fig. 3A to C and used as a template. Individual reaction tubes were removed from the thermocycler after 3 to 40 amplification cycles. Products were separated on a 10% denaturing acrylamide gel and detected as described above. The intensity of the radioactive signal was estimated by image densitometry with background correction by using the Optiquant 3.1 software (Packard Instruments).
Determination of probe purity and Td.
The size and purity of a PCR-generated oligonucleotide from E. coli was compared to those of purchased synthetic probes that had been either desalted or gel purified. Purified oligonucleotides were labeled with 32P as described above and separated by size on 10% denaturing acrylamide gels, followed by phosphorimaging detection.
For determination of dissociation temperatures (Td), total nucleic acids extracted from E. coli were blotted on Zetaprobe nylon membranes (Bio-Rad) with a Minifold II slot blot apparatus (Schleicher & Schuell, Keene, N.H.) after a 30-s alkaline denaturation (10 mM NaOH, 1 mM Tris) at 85°C. Each slot contained 150 ng of total nucleic acids (approximately 50 ng of 16S rRNA). After each slot was washed with 200 µl of diethyl pyrocarbonate-treated water (16), nucleic acids were UV-cross-linked (120 mJ) in a Stratalinker 2400 (Stratagene) and rinsed in diethyl pyrocarbonate-treated water before drying.
The Td of the PCR-generated and the purchased synthetic probes were determined as described by Polz and Cavanaugh (15). Membrane strips containing triplicate slots of E. coli nucleic acids were prewetted in 2 ml of hybridization buffer (5x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate], 1% SDS) at 37°C for 30 min, followed by addition of 10 µl of buffer containing 1 pmol of probe per slot. Membranes were hybridized overnight at 37°C, and unbound probe was removed by two 15-min washes at the same temperature in washing buffer (1x SSC, 0.1% SDS). Subsequently, individual slots were cut out and washed sequentially at gradually increasing temperatures (37 to 87°C) by transferring them to 7-ml scintillation vials containing 2 ml of washing buffer for 5 min at a time. The radioactivity remaining in the wash solutions and activity remaining on the filter after the final wash were determined by liquid scintillation after adding 3 ml of Scintisafe Plus 50% (Packard) to each vial. The total activity in the wash solutions and the activity remaining on the membrane were summed, and the Td was estimated as the temperature at which 50% of the probe was dissociated from the membranes.
Hybridization specificity.
To test the hybridization specificity of the PCR-generated probes, nucleic acids from 30 different closely and distantly related bacterial strains (Fig. 1; see Fig. 7A) were blotted on Zetaprobe nylon membranes as described above. A eubacterial probe, EUB338 (20), and PCR-generated population-specific probes from E. coli and Vibrio anguillarum were 32P labeled at the 5' termini with bacteriophage T4 polynucleotide kinase (16). Each membrane was incubated overnight at 37°C with 4 pmol of the respective probe. Subsequently, membranes were washed twice for 20 min at 37°C in washing buffer (1x SSC, 0.1% SDS) before a final stringent wash for 30 min at a higher temperature (EUB338 at 50°C and the longer PCR-generated probes at 70°C). Hybridization was detected by phosphorimaging as described above.
![]() View larger version (107K): [in a new window] |
FIG. 7. Membrane hybridization of total nucleic acids from 30 bacterial strains with 32P-labeled oligonucleotide probes. (A) Identity and location of immobilized nucleic acids on the membrane; (B) hybridizations with universal bacterial probe EUB338; (C) PCR-generated probe specific for E. coli; (D) PCR-generated probe specific for V. anguillarum.
|
|
|
|---|
Alkaline hydrolysis of 3' ribonucleotides is a robust and efficient way to detach primers from amplification products. This is evident from the complete hydrolysis of primer 104r following incubation in 0.2 M NaOH with no apparent damage to the DNA (Fig. 4) (17). The mechanism of hydrolysis ensures that the ribonucleotide is removed with the primer. An OH--induced attack and subsequent breaking of the 3' phosphodiester bond produces an OH group at the 5' terminus of the synthesized probe (13). Experiments with primer 64f showed a similar efficiency of detachment (data not shown). In contrast to its sensitivity to alkali, the 3' ribonucleotide is stable in the presence of RNase, since addition of the enzyme did not affect amplification and alkaline detachment of the primer (data not shown).
![]() View larger version (92K): [in a new window] |
FIG. 4. (A) Phosphorimaging scan showing detachment of the 5'-32P-labeled primer from the PCR product after alkaline hydrolysis of the ribonucleotide at the 3' end of primer 104r. PCR products with primer 104r lacking the 3' ribonucleotides are included for comparison. Lanes: 1, primer 104r; 2, PCR product with 3' ribonucleotide-modified primer; 3, PCR product with 3' ribonucleotide-modified primer after alkaline hydrolysis; 4, PCR product with nonmodified primer; 5, PCR product with nonmodified primer after alkaline hydrolysis.
|
![]() View larger version (36K): [in a new window] |
FIG. 5. Linear amplification with 5'-32P-labeled primer 104r. (A) Gel electrophoresis and phosphorimaging detection of amplified fragments after 3 to 40 cycles. (B) Quantification by image analysis using local background subtraction shows a linear accumulation of amplification product with cycle number. AU, absorbance units.
|
![]() View larger version (30K): [in a new window] |
FIG. 6. (A) Gel electrophoresis of a 5'-32P-labeled PCR-produced oligonucleotide probe targeting 16S rRNA of E. coli. Two synthetic oligonucleotides of identical sequence (desalted or purified by acrylamide gel electrophoresis) are included as a reference. (B) Dissociation curves for a commercially produced, gel-purified oligonucleotide probe and a PCR-generated oligonucleotide probe with identical sequence (5'ACTCGTCAGCAAAGAGCAAGCTTCTTCCTGTTACCGTT3'). Each symbol represents the mean dissociation of probe hybridized to triplicate slot blots with total nucleic acids from E. coli as a target ( 50 ng of 16S rRNA). Error bars represent ±1 standard deviation, and dashed lines indicate the Td for the two probes.
|
20-nucleotide probing region. The only exceptions in this group are the two Sulfitobacter strains, S. pontiacus and Sulfitobacter sp. strain EE-36 which have an identical probing region and show only 0.4% variation in their 16S rRNA sequences. Other strains that were not well resolved at the isolate level include three of the Synechococcus and Prochlorococcus isolates (1 to 3.9% sequence variation in 16S rRNA), the two Zooglea isolates (1.5% sequence variation in 16S rRNA), and five Pseudomonas isolates (1.1 to 4.5% sequence variation in16S rRNA). This overall picture of probe specificity was confirmed by detailed comparison of deeply branching and distantly related bacterial groups, such as the thermophilic oxygen reducers, Thermotogales, green nonsulfur bacteria, deinococci, fusobacteria, and actinomycetes (data not shown). The PCR-generated probes also showed good ability to differentiate 30 closely and distantly related bacterial strains in slot blot hybridization of total nucleic acids (Fig. 7). Each slot contained 150 ng of total nucleic acids, and the relative abundances of 16S rRNA molecules in each slot, assessed by hybridization with the 32P-labeled Bacteria-specific probe EUB338, appeared to be similar in magnitude (Fig. 7B). Hybridizations of PCR-generated probes targeting E. coli and V. anguillarum showed high specificity for the respective target under the tested conditions (Fig. 7C and D), with the only exception of slight cross-hybridization of the E. coli probe with Pantoea agglomerans nucleic acids. This was surprising and remains unexplained at this point, since the target sequences are much more dissimilar than other sequences that were perfectly differentiated in the hybridization.
|
|
|---|
The generation of population-specific probes by PCR has previously been attempted. Ludwig et al. (9) used a 148- to 358-nucleotide-long variable region within bacterial 23S rRNA flanked by conserved primers as hybridization probes, with good success. However, variable regions within the 16S rRNA are significantly shorter, making primer removal necessary. This was recently achieved by enzymatic removal of primers carrying a 3' phosphorothioate followed by digestion of the remaining complementary overhang by mung bean nuclease (6). The resultant probe was a double-stranded polynucleotide that was used to identify 16S rDNA in temperature gradient gel electrophoresis gels for microbial community analysis. Unfortunately, this procedure did not appear to generate probes of uniform size, and the hybridization efficiency is likely impaired by the presence of a complementary strand that competes with the target. Thus, we propose the protocols presented here as a generally applicable improvement to generation of hybridization probes via PCR.
The probing region identified in the present study represents one of the most variable stem structures in the 16S rRNA (http://rrna.uia.ac.be/) (22) and is thus a suitable target for identification of bacterial populations. Nonetheless, there is no single region within the 16S rRNA that will differentiate all closely related sequences (18). This is evident in our analysis of several isolates of Pseudomonas, Roseobacter, Prochlorococcus, and Synechococcus, which display on the order of 0.4 to 4.2% sequence variation along the entire 16S rRNA but are not distinguishable in the probing region. Thus, if population-level identification of a specific microorganism were required, it would be prudent, as in any other hybridization approach, to use more than one probe. This can be accomplished by using additional PCR-generated probes, since we have recently identified at least two suitable regions within the 23S rRNA (unpublished).
Several steps within the probe generation protocol were optimized to ensure reproducible production of high-quality, uniform oligonucleotides. First, it is important to generate a single-stranded product by removing the complementary strand from the PCR product after the first amplification prior to hydrolysis of the ribonucleotide (Fig. 2B). This is necessary to avoid regeneration of the priming region in subsequent amplification when the strand shortened by the detachment of the primer can act as a primer for the complementary strand. The removal is efficiently achieved by washing the bead-bound PCR product for 1 min at 85°C. Washes at higher temperature or at 37°C in presence of 50 to 100% formamide as recommended by the manufacturer resulted in the loss of most bound product, probably due to dissociation of the biotin-streptavidin bond. Second, to maximize the final yield of probe after the second linear amplification, it is desirable to have a nearly complete incorporation of biotinylated primer into products, since free primers may compete with PCR product for available streptavidin binding sites on the beads. It is possible to compensate for this by adding a larger amount of beads, but this will increase costs for production of probes substantially.
The PCR-generated probes displayed melting behavior typical of oligonucleotides (Fig. 5B). However, the Td of PCR-generated probes were 5°C higher than those of gel-purified synthetic probes of identical sequence (Fig. 5A). This may be due to the lower purity of the synthetic probes (Fig. 5B), but the difference could also be due to some unknown structural differences between synthetic and PCR-generated nucleic acids. Polz and Cavanaugh (15) and McMahon et al. (11) observed a similar difference in Td between native and in vitro-transcribed 16S rRNAs, but the observations also remained unexplained. The substantial length heterogeneity of the probing region (Fig. 1) may result in Td variability between different probes. This may complicate the use of the PCR-generated probes in modern chip technologies, and further development of the method may be required. One possible solution to this problem is to use other variable regions with less length heterogeneity for PCR production of probes.
In summary, the PCR-based probe synthesis outlined in this study provides the means to rapidly generate oligonucleotide probes for identification of bacteria at a low cost. The region identified in this study is not perfect in differentiation of all closely related bacteria, but the specificity of the analysis can readily be improved by combined use of similarly produced probes targeting other variable regions. The procedure can easily be adapted to a 96-well microplate format with only minor modifications, facilitating high-throughput production of probes. Both genomic DNAs from pure cultures and cloned 16S rDNA fragments can be used as templates for probe synthesis, and hence the method is suitable for rapid hybridization-based screening of samples. Other possible applications are the high-intensity labeling of fluorescent in situ hybridization probes by incorporation of fluorescent nucleotides during the PCR or the production of oligonucleotide probes for microarray-based identification of bacterial populations (3).
This work was supported by the K&A Wallenberg Foundation (postdoctoral fellowship to S.B.) and by grants from the NSF to C.M.C. and M.F.P.
Present address: Department of Evolutionary Biology, Limnology, Uppsala University, SE-75236 Uppsala, Sweden. ![]()
|
|
|---|
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»