Terminal restriction fragment length polymorphism (t-RFLP) analysis is a powerful comparative fingerprinting molecular technique that is widely used to describe microbial community structure for different temporal and spatial variations, as well as geochemical alterations caused by natural and/or engineered interventions (
1,
8,
13,
32,
38-
40,
46). In t-RFLP analysis, the phylogenetic signatures (such as that of the gene encoding 16S rRNA) of total community DNA are first PCR amplified with fluorescence-labeled primers, and the PCR products are subsequently digested with a restriction enzyme to generate a series of peaks corresponding to various fluorescence-labeled terminal restriction fragment (T-RF) lengths. The profile of a series of T-RFs provides an estimate of the number of phylotypes in a microbial community, and the fluorescence intensity of each peak reflects the relative abundance of each phylotype (theoretically) in the microbial community (
39). Because many phylogenetically similar species can contribute to the same T-RF, T-RFs are commonly referred to as operational taxonomic units. Nevertheless, the t-RFLP assay is still considered a useful fingerprinting method, particularly for comparative analysis of microbial communities.
Real-time PCR analysis has been accepted as a very sensitive quantitative molecular method that has applications ranging from clinical microbiology (
4,
5,
24,
44) to molecular ecology (
25,
53,
54) and environmental microbiology (
3,
10,
18,
22,
28,
33,
34,
37). Unlike other quantitative molecular techniques, real-time PCR is a fluorescence-based kinetic PCR method, and it allows rapid quantification of the initial number of copies of the target gene (
29). In TaqMan real-time PCR, an extra fluorescence-labeled probe complementary to the target DNA is initially bound to the target DNA. As PCR progresses, the fluorescent label is cleaved to release a fluorescent signal that is recorded at the end of each PCR cycle (
23,
31). The detected fluorescent signals are compared to a set of standard curves for quantification of the target gene in the samples. In addition to its sensitivity for quantification, the real-time PCR assay can be designed for a specific strain, a phylogenetic group, or even a group of microorganisms exhibiting a similar function (
49).
Currently, researchers are able to quantify a specific species of interest by using real-time PCR or are able to acquire qualitative information for microbial community structure by using molecular techniques, like t-RFLP analysis. However, integrating results obtained from various molecular assays in order to describe microbial diversity and the absolute abundance of species populations in a given community is quite a challenge and may not even be possible. Therefore, a single quantitative molecular method that allows effective fingerprinting of a given microbial community will be particularly valuable in microbe-related studies.
In this study, a single molecular assay, the real-time-t-RFLP assay, which allowed simultaneous quantification of all microbial ribotypes in a microbial community, was developed. This assay was applied to environmental samples (wastewater) and enrichment cultures (estrogen-degrading consortia). The estrogen-degrading consortia enriched from activated sludge were capable of degrading three different estrogens, 17α-estradiol, 17β-estradiol, and estrone. The limitations and potential applications of the real-time-t-RFLP assay are also discussed.
MATERIALS AND METHODS
Chemicals.
Three natural estrogens (17α-estradiol, 17β-estradiol, and estrone) were used in this study. 17β-Estradiol (>99% pure), 17α-estradiol (>98% pure), and estrone (>98% pure) were purchased from Sigma-Aldrich Inc. (St. Louis, Mo.). Stock solutions of estrogens were prepared in acetone. Dimethylformamide and N,O-bis(trimethylsilyl)trifluoroacetamide were purchased from Pierce Chemicals, Dallas, Tex. Hi-Di formamide and the GeneScan 500 ROX size standard were purchased from Applied Biosystems, Warrington, United Kingdom.
Bacterial cultures.
Three bacterial strains were used in this study.
Escherichia coli (gram negative),
Pseudomonas fluorescens (gram negative), and
Bacillus thuringiensis (gram positive) were obtained from the culture collection at the Center for Environmental Biotechnology, University of Tennessee, Knoxville. The 16S rRNA gene sequences of
E. coli,
P. fluorescens, and
B. thuringiensis exhibit 99% identity to the sequences deposited under GenBank accession numbers
AF233451,
AF228366, and
AF155954, respectively. The bacterial strains were grown in Luria-Bertani medium at 30°C overnight before they were harvested by centrifugation for experimental use.
Estrogen-degrading consortia.
Estrogen-degrading consortia were enriched from activated sludge from two wastewater treatment plants (WWTPs) (designated WWTP1 and WWTP2) near Knoxville, Tenn. WWTP1 treats 1.9 × 104 m3 of wastewater per day and consists of an oxidation ditch activated sludge system. Wastewater samples were collected from the following three locations: primary treatment (i.e., after the primary treatment serving as the influent for the ditch), the center of the ditch, and the clarifier. WWPT2 treats 1.1 × 105 m3 of wastewater per day and consists of a one-sludge nitrification system. A sludge sample was collected from the nitrified activated sludge tank of WWTP2.
Enrichment of the estrogen-degradation consortia was performed in 40-ml glass vials containing 10 ml of nitrate mineral salts medium (
14), 1 ml of activated sludge as the inoculum, and 1 mg of estrogen (17α-estradiol, 17β-estradiol, or estrone) per liter. The inoculated vials and abiotic controls (without activated sludge) were incubated at 30°C at 150 rpm. The abiotic loss of estrogen was less than 5%. The experiment was performed in duplicate. Weekly subculturing was performed by transferring 1 ml of the liquid suspension to a new vial containing nitrate mineral salts medium with an estrogenic compound added as described above. After five transfers, enrichment cultures were collected for microbial community structure analysis by the real-time-t-RFLP assay. The microbial community structure after enrichment with three different estrogenic compounds was compared to the structure in the original activated sludge samples.
Analytical methods. (i) DNA extraction.
Genomic DNA of each bacterium was extracted with a FastDNA kit (Q-Biogene, Carlsbad, Calif.) in accordance with the manufacturer's instructions. For environmental samples (wastewater and activated sludge) and enrichment cultures, a FastDNA SPIN kit for soil (Q-Biogene) was used with the following minor modifications: the silica binding matrix-DNA complex was washed twice with 80% (vol/vol) ethanol after the recommended salt-ethanol wash step (
18). DNA concentrations were determined with a Hoefer DyNa Quant 200 fluorometer (Pharmacia Biotech, San Francisco, Calif.).
(ii) Real-time-t-RFLP assay.
A region of the 16S rRNA gene sequence (length, 352 bp) was selected for development of the real-time-t-RFLP assay in this study, because primers (16S1055f and 16S1392r) and a TaqMan probe (16STaq1115f) for this region have been designed and successfully tested with wastewater samples (
28). With no base pair mismatches, the primers and probe capture more than 9,000 bacterial 16S rRNA gene sequences in the GenBank database (
28) and generate a longer PCR amplicon (352 bp versus 263 bp) than the other primer set and probe for total bacterial 16S rRNA genes (
2). The target region in extracted DNA was amplified by using the primers and the TaqMan probe, with the following modifications. The PCR was performed with fluorescence-labeled forward primer 16SHex1055f (5′-hexachlorofluorescein-ATGGCTGTCGTCAGCT-3′), reverse primer 16S1392r (ACGGGCGGTGTGTAC-3′), and TaqMan probe 16STaq1115f (5′-[6-carboxyfluorescein]-CAACGAGCGCAACCC-[6-carboxytetramethylrhodamine]-3′). The PCR mixture (total volume, 25 μl) contained 2.5 μl of 10× PCR buffer A, 1.25 U of
Taq polymerase (Fisher Scientific, Fair Lawn, N.J.), 0.5 μg of bovine serum albumin, each of four deoxynucleoside triphosphates at a concentration of 200 μM, 4.5 mM MgCl
2, each primer at a concentration of 600 nM, 250 nM TaqMan probe 16STaq1115f, and 10 ng of DNA template. PCR amplification was performed with a DNA Engine Opticon continuous fluorescence detection system (MJ Research, Waltham, Mass.). The PCR conditions were as follows: initial denaturation at 95°C for 10 min, followed by 30 cycles of 95°C for 30 s, 54°C for 1 min, a plate read step, and 72°C for 2 min, and a final incubation at 4°C. Standard curves for quantifying 16S rRNA gene copies in samples were constructed by using plasmid #931, which carries a
Nitrospira partial 16S rRNA gene (GenBank accession number
AF420301) (
18,
28). Standard curves based on amounts ranging from 2.3 × 10
3 to 2.3 × 10
8 copies of the 16S rRNA gene were generated parallel to the sample analysis.
PCR products were separated on a 1.5% agarose gel in 1× Tris-acetate-EDTA buffer, and the DNA bands of the expected size (352 bp) were excised. The PCR products were then recovered and purified by using MicroSpin columns (Amersham Biosciences, Piscataway, N.J.), followed by ethanol precipitation and two desalinization steps with 80% ethanol.
Purified PCR products were then digested with restriction enzyme MspI at 37°C for 8 h following 20 min of incubation at 70°C. The digestion reaction mixture (30 μl) contained 60 to 90 ng of the purified PCR products and 10 U of MspI (or HhaI). Digested PCR products were precipitated with ethanol and resuspended in 40 μl of high-performance liquid chromatography water.
The lengths of T-RFs of digested PCR products were determined with an ABI Prism 310 genetic analyzer (Applied Biosystems, Foster City, Calif.) in GeneScan mode. Each sample was prepared by adding 5 μl of desalted digested PCR product, 10 μl of Hi-Di formamide, and 0.5 μl of Genescan ROX 500 size standards. Samples were denatured at 95°C for 5 min, followed by rapid chilling on ice. By using an injection time of 10 s, a final DNA concentration in the range from 0.05 to 0.7 ng/μl was introduced electrokinetically into a capillary containing polymer (GS STR POP4). Electrophoresis was performed at 15 kV at 60°C for 25 min. The lengths of T-RFs were automatically determined by comparison with the internal standards by using the GeneScan software, version 3.1.
(iii) Estrogen analysis.
Concentrations of estrogens in the growth medium were determined as described by Raman et al. (
50). The estrogen in a liquid suspension was extracted with an equal volume of ethyl ether (pesticide grade; Fisher Scientific). The ether layer (500 μl) was carefully transferred to 2-ml autosampling gas chromatography vials and evaporated under a gentle stream of nitrogen. The dried extracts were resuspended with 450 μl of dimethylformamide and then derivatized with 50 μl of
N,
O-bis(trimethylsilyl)trifluoroacetamide. The derivatized samples (1 μl) were injected into an HP 6890 series gas chromatograph equipped with a DB-5MS capillary column (length, 30 m; film thickness, 0.25 μm; J&W Scientific) and an HP 5973 mass selective detector. The analysis was performed in the selective ion monitoring mode. Standard curves for estrogens, containing at least five data points ranging from 0.001 to 5 mg/liter, were used for calculation of estrogen concentrations.
Data acquisition and analysis.
The number of 16S rRNA gene copies of the
ith T-RF per unit volume (
Ci) in a microbial community can be determined as follows:
where
Co is the initial number of gene copies in the sample, expressed in number of 16S rRNA gene copies per milliliter;
jAi is the peak area of the
ith T-RF of
j length (in bases) that is measured from the electropherogram; and
n is the total number of T-RFs in the electropherogram.
Co was determined by comparing the threshold cycle (
Ct) mesured from samples to standard curves (
28). The
Ct was defined as the number of PCR cycles when the generated fluorescence exceeded 10 times the standard deviation of background fluorescence. The
Ct value was determined automatically with computer software (Opticon Monitor, version 1.4) installed in a DNA Engine Opticon continuous fluorescence detection system (MJ Research). Standard curves were constructed by plotting the
Ct values over a range of known numbers of gene copies per unit volume. The ratio of each peak area to the total peak area,
\( \(\frac{^{j}A_{i}}{{{\sum}_{i{=}1}^{n}}(^{j}A_{i})}\) \), can be considered the relative abundance of the
ith T-RF in a microbial community.
Since T-RFs less than 30 bp long were observed in all samples and their intensities were independent of the input DNA, the peaks of these T-RFs were excluded from the data analyses. A threshold peak area value of 1% has been used in other studies (
11,
43), and use of a 5% relative peak area threshold has been shown to minimize the effects of the amount of labeled products (one of the significant factors) affecting t-RFLP profiles (
51). In this study, peak area that contributed less than 1% of the total area was regarded as background noise and was excluded from the analysis.
DISCUSSION
In studies performed with other fingerprinting methods, such as t-RFLP analysis and denaturing gradient gel electrophoresis, workers observed fallacious microbial community structure due to the presence of artificial peaks that commonly resulted from incomplete digestion of PCR products (
15,
19,
38,
45,
48). For real-time-t-RFLP analysis, complete digestion of the PCR product is essential for obtaining accurate quantification of each ribotype since only complete digestion can result in a correct percentage of each ribotype (T-RF) and the percentage is subsequently used to calculate the number of 16S rRNA gene copies for each ribotype. To optimize the real-time-t-RFLP assay for complete digestion of PCR products, different mass ratios of PCR products to restriction enzyme (MspI) were examined by using two different digestion procedures: a single-step procedure (MspI was added at the beginning) and a two-step procedure (equal doses of MspI were added in the beginning and 6 h after the first dose). PCR products were digested under recommended digestion temperature and duration conditions. It was found that a single 10-U dose of MspI was sufficient to completely digest 60- to 90-ng PCR products in a 30-μl digestion reaction mixture (results not shown), as described in Materials and Methods.
In addition to incomplete digestion, many other factors (including nonspecific PCR products) can result in artificial peaks in t-RFLP profiles (
15). In this study, a factor that caused artificial peaks that has not been reported previously was noted. When a digested DNA sample was overloaded during capillary electrophoresis, artificial peaks occurred at sizes 6 to 10% shorter than the sizes of measured T-RFs. As shown in Fig.
7D, an artificial peak at 48 bp was observed along with the measured T-RF (51 bp) of HhaI-digested PCR products when plasmid #931 was used as a template. When the loaded sample sizes were decreased proportionally (2, 1, 0.5, and 0.25 ng/μl), the area under the T-RF peak at 48 bp decreased rapidly but not proportionally. As expected, the area under the T-RF peak at 51 bp decreased in proportion to the loaded sample size. Furthermore, the area ratios of the 48-bp T-RF peak to the 51-bp T-RF peak did not remain constant, indicating that the 48-bp T-RF peak was an artificial peak. To avoid generation of artificial peaks, digested samples were loaded at concentrations of less than 0.5 ng/μl. Nevertheless, it is also recommended that each sample be run with at least two different dilutions during capillary electrophoresis to avoid overestimating the diversity and/or underestimating the abundance of ribotypes in a microbial community.
Selection of a proper amplicon size for a real-time-t-RFLP assay can be a challenge due to the inherent limitations of the real-time PCR and t-RFLP techniques. For TaqMan real-time PCR, the recommended amplicon length is less than 400 bp because shorter amplicons allow effective binding of probes and primers to the complementary sequences (
12). For t-RFLP analysis, the recommended size is between 400 and 700 bp, an amplicon size that is long enough to include sequence divergence among microorganisms and at least one restriction site (
38). Based on both criteria, the real-time-t-RFLP assay with an amplicon size of 352 bp was successfully developed. According to TAP analysis of the 16S rRNA gene sequences in the Ribosomal Database Project database (allowing at most no mismatch in the last four nucleotides from the 3′ end of the forward primer and three mismatches in any other nucleotide), this amplicon should generate 188 unique T-RFs when a sample is digested with MspI and only 157 of 11,274 primed sequences (1.4%) showed no MspI restriction sites in this 352-bp amplicon. In this study, the 352-bp T-RF was observed in only two digested samples, the clarifier sample and the 17α-estradiol-degrading consortium enriched from the clarifier sample from WWTP1 (<5% of the total area under T-RF peaks), suggesting that there was a lack of restriction sites for MspI in amplicons of these two samples. As the amplicon used in this study was not located in variable region I in the 16S rRNA gene and the first 500 bp of variable region I provides better resolution for t-RFLP analysis (
39), designing different amplicons (such as amplicons targeting a different region of the 16S rRNA gene) and using different restriction enzymes to digest amplified PCR products should improve the resolution and avoid undigested ampilcons.
Although the real-time-t-RFLP assay is subject to systematic biases like any PCR-based method (
56) and has the same limitations as t-RFLP analysis (
39), the application of the real-time-t-RFLP assay to environmental samples was considered successful. The microbial ecology of activated sludge is very complex and has not yet been fully illustrated. Many predominant bacterial groups in activated sludge have been identified by using a wide range of molecular techniques, including clonal library analysis (
7), in situ hybridization (
36,
52), denaturing gradient gel electrophoresis (
47), and t-RFLP analysis (
21,
57). The predominant bacterial phyla are the
Proteobacteria,
Planctomycetes,
Bacteroidetes,
Firmicutes, and
Actinobacteria. As shown in Table
2, the microbial ecology in wastewater and activated sludge analyzed by the real-time-t-RFLP assay (Fig.
4 and
5) was consistent with the microbial ecology determined in other studies. For example, the measured T-RFs are consistent with the expected T-RFs of a wide range of known microorganisms, including foam-causing, filamentous, nitrifying-denitrifying, and presumptive phosphorus-accumulating bacteria. By developing different combinations of amplicons and restriction enzymes, workers may be able to quantify specific functional groups of microorganisms by the real-time-t-RFLP assay. Quantitative information pertaining to specific functional groups of microorganisms not only provides interesting insights into these functional communities but may also provide a fundamental basis for fine-tuning operational parameters to prevent foaming and sludge bulking, as well as to enhance the nitrogen and phosphorus removal efficiency in treatment units. The quantitative feature of the real-time-t-RFLP assay is a powerful tool for understanding the microbial ecology underlying microbially meditated processes.
Biodegradation of estrogen in activated sludge has been strongly suggested in previous studies; however, the diversity and abundance of estrogen-degrading cultures in activated sludge remain unclear. The changes in the microbial diversity and quantities in response to addition of three different types of estrogens (Fig.
6) indicated the potential presence of estrogen-degrading cultures in activated sludge. Studies which focus on isolation and characterization of estrogen-degrading cultures in activated sludge should be useful for the development of optimal operating conditions for better estrogen removal from wastewater.
The strategy used in real-time-t-RFLP assay development can be extended from eubacteria to archaea, a specific phylogenetic group, and/or functionally similar groups. In this study, the 16S rRNA gene was used to target all eubacteria. Similarly, some functional genes can be used for development of a method for quantifying a specific group of functional microorganisms. For example, the
amoA genes have shown higher diversity than the 16S rRNA gene in ammonia-oxidizing bacteria (
49) and can be further adapted into a real-time- t-RFLP assay for studying ammonia-oxidizing bacteria. For studies designed to increase the T-RF resolution at the species level, the 16S-23S intergenic region that has significant heterogeneity in both length and nucleotides (
20) can be considered during design of the amplicon.
The real-time-t-RFLP assay will be particularly useful for in situ or ex situ bioremediation in which the microbial activity changes in response to increases or decreases in specific microbial populations. The information obtained should allow better assessment of the treatment duration required for achieving target clean-up goals. In addition, this assay should also be useful for studying microbial ecology in soils and water in which the microbial communities vary due to spatial or temporal differences, pH variation, nutrient availability, and many engineered biological treatment processes (such as processes in bioreactors) in which microbial populations shift in response to environmental toxicants and/or stimulants.