Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2004, p. 4971-4979, Vol. 70, No. 8
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.8.4971-4979.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Ingrid Brettar,1 Richard Christen,2 Sebastien Flavier,2 Julia Bötel,1 and Manfred G. Höfle1*
Department of Environmental Microbiology, GBF-German Research Centre for Biotechnology, 38124 Braunschweig, Germany,1 UMR 6543 CNRS and Université Nice Sophia Antipolis Centre de Biochimie, F-06108 Nice, France2
Received 13 December 2003/ Accepted 9 April 2004
|
|
|---|
|
|
|---|
Real-time PCR is a new and highly sensitive method for the quantification of microorganisms from different samples (40). Typical quantitative PCR (Q-PCR) techniques utilize approaches originally developed for clinical applications. Q-PCR has been used for the quantification of taxa and target genes from various bacteria, such as Staphylococcus aureus (15), Bacillus anthracis (24, 30), and Legionella spp. (45), and from viruses (11) or for the quantification of insulin-like growth factors (28). Recently, applications of real-time PCR have also been demonstrated with marine samples to quantify diatoms and pelagophytes (41) and numbers of cells of Synechococcus spp., Prochlorococcus spp. (37), and Pseudomonas stutzeri (12); with soil samples to quantify cells of Rhodococcus spp. (31), methanotrophs (19), and fungi (9); with sediments to quantify cells of Geobacter spp. (36); in wastewater treatment plants to quantify nitrifying bacteria (13); and with cultures to determine dissimilatory (bi)sulfite reductase gene expression in Desulfobacterium autotrophicum (26). One advantage of PCR techniques, apart from their great potential for automation, is their sensitivity. Theoretically, PCR should detect one copy of the 16S rRNA gene (39). Lockey et al. (21) described a real-time PCR method and validation strategy useful to detect a single-copy herpes simplex virus thymidine kinase gene. Moreover, a single B. anthracis spore was detected in 100 liters of air sampled by real-time PCR (24). However, for the most part, sensitivity becomes limited below 100 gene copies per reaction mixture (12, 35).
Sensitivity and specificity can be improved by adopting a "nested" or a "heminested" PCR approach, whereby initial amplification is carried out with a pair of primers with broad specificity (14, 38). This has been successfully used to detect autotrophic ammonia-oxidizing bacteria in samples of lake water (14) or in combination with species-specific PCR product melting temperatures (with SybrGreen) to identify pathogenic Campylobacter species (23).
Our main interest is the ecology of microbial catalysts involved in major biogeochemical processes of the central Baltic Sea. Previous studies (3, 4) have shown that autotrophic denitrification may be an important process for N balance in the central Baltic Sea. Thus, in the present study, bacteria responsible for autotrophic denitrification were chosen as a target for our research. In a combined experimental and molecular approach, a bacterium related to Thiomicrospira denitrificans was identified as a major catalyst. Quantification of bacterial catalysts independent from their ribosome contents or the necessity of culturing these bacteria is a tool urgently needed for elucidating their ecology. We developed a specific Q-PCR technique allowing quantification of T. denitrificans-like bacteria both in experimental enrichments and in environmental samples. The method is based on specific primers (OST 1R and OST 1F) designed from a consensus sequence obtained from community fingerprints of Baltic bacterioplankton (M. G. Höfle, S. Flavier, R. Christen, J. Bötel, M. Labrenz, and I. Brettar, submitted for publication). The method gave reliable results when it was applied to DNA as well as to RNA extracted from bacterioplankton. Our approach allowed the in situ detection of a bacterial catalyst under consideration down to 90 cells ml1 or a relative abundance of 0.01% without the need of absolute standards in real-time experiments. Based on total bacterial cell counts, the number of ribosomes per cell could also be estimated. This approach may help to improve quantification of major bacterial catalysts, especially for oligotrophic environments, and allow us to gain new insights into their in situ ecology.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Comparison of total 16S rRNA contents from seawater samples and enrichmentsa
|
Primers.
The application of 16S ribosomal complementary DNA (rcDNA) single-strand conformation polymorphism analysis and sequencing of dominant bands (29) demonstrated that an autotrophic denitrifying bacterium related to, but not identical with, T. denitrificans was the dominant autotrophic denitrifier in the above-described enrichments (Höfle et al., submitted). Subsequent development of highly specific primers against this bacterium has also been described in detail (Höfle et al., submitted). T. denitrificans-like forward primer OST 1F (5'-TCAGATGTGAAATCCAATGGCTCA-3') corresponds to 16S rRNA Escherichia coli positions 663 to 686, and reverse primer OST 1R (5'-CTTAGCGTCAGTTATGTTCCAGG-3') corresponds to E. coli positions 803 to 825. Primers used for bacterial 16S rRNA gene and 16S rcDNA amplification from environmental DNA or RNA are given by Schwieger and Tebbe (34). Primer set Com 1 (5'-CAGCAGCCGCGGTAATAC-3') and phosphorylated Com 2ph (5'-CCGTCAATTCCTTTGAGTTT-3') amplify positions 519 to 926 of the E. coli 16S rRNA gene.
16S rcDNA generation.
To retrieve 16S rcDNA, 30 ng of template RNA was reverse transcribed at 50°C with the Thermoscript reverse-transcription (RT)-PCR system (Invitrogen) and the universal reverse primer 1492R (5'-GGTTACCTTGTTACGACTT-3') (20), according to the manufacturer's instructions. In each reverse transcription reaction mixture, some RNA samples were not supplemented with reverse transcriptase to rule out DNA contamination.
Relative quantification protocol without a preamplification step.
A Rotor-Gene real-time cycler (Corbett Research), the Brilliant SybrGreen Q-PCR master mix (Stratagene), and nucleic acids isolated from Baltic seawater experiment (Bx) enrichments most likely containing large amounts of the target species were used for the establishment of Q-PCR and quantitative RT-PCR (Q-RT-PCR) protocols. Prior quantifications, minimum Com and OST primer concentrations (in a range from 50 to 900 nM for each primer), and optimal template concentrations were determined. For optimal primer binding, a temperature of 59°C was determined, where no binding to the sequence of the most closely related bacterium which could be cultured, i.e., T. denitrificans, occurred (Fig. 1). The optimized PCR mixtures (25 µl) contained 1x Brilliant SybrGreen Q-PCR master mix (Stratagene), the Com (a 900 nM concentration of each primer) or OST 1 (OST 1F, 900 nM; OST 1R, 300 nM) primer system, and 2 to 4 ng of nucleic acids extracted from the Bx enrichments. An initial denaturing step at 95°C for 15 min was followed by 40 cycles of 94°C for 30 s, 59°C for 40 s, and 72°C (a measuring step) for 50 s. Potential development of primer dimers was determined by a melting-point analysis in a range from 40 to 94°C. Each measurement was performed multiple times (at least three replicates).
![]() View larger version (53K): [in a new window] |
FIG. 1. Agarose (1.2%) gel of RT-PCR amplification products based on the OST 1F-OST 1R (OST) or Com 1-Com 2ph (Com) primer system from various template RNAs. Templates originated from the thiosulfate-plus-nitrate enrichment Bx21, T. denitrificans DSM 1251T, or E. coli K-12. The first lane contains molecular size markers (100 bp). NC, negative control.
|
Real-time data analysis.
Relative standards were prepared using serially diluted (1:10) nucleic acids isolated from exactly the same DNA or RNA extract. Each sample measured was calibrated by its own respective Com and OST standards. Due to potentially lower PCR efficiencies of nondiluted DNA or cDNA extracts, each sample to be measured was always diluted 1:10. Calibration curves were generated by Rotor-Gene software version 4.6. For each standard, the concentration was plotted against the cycle number at which the fluorescence signal increased above the background or cycle threshold (CT value). The slope of each calibration curve was incorporated into the following equation to determine the reaction efficiency: efficiency = 10(1/slope) 1. According to this formula, an efficiency of 1 means a doubling of product in each cycle. Only data resulting from measurements with comparable Com- and OST-specific PCR efficiencies were used for further calculations.
Calculation of numbers of T. denitrificans-like cells and 16S rRNA molecules.
In order to calculate numbers of T. denitrificans-like cells, total cell numbers were used and the Com product/OST product ratios were obtained from relative 16S rRNA gene quantifications, assuming that all cell counts represented bacteria. As with the calculation of T. denitrificans-like 16S rRNA molecule numbers, total 16S rRNA molecule numbers were used and the Com product/OST product ratios were obtained from relative 16S rcDNA quantification analyses. Total 16S rRNA molecule numbers were calculated as follows. Total RNA was quantified as described above. The RNA extraction efficiency of our phenol extraction protocol has been demonstrated to range from 52 to 66% and to be relatively constant for gram-negative bacteria (43). Therefore, we used a 60% extraction efficiency to calculate the real amount of total RNA based on the fluorescence measured. Neglecting tRNA and mRNA, which are minor components of the total RNA pool, the 23S and 16S rRNAs are present in a ratio of 2:1. Consequently, 16S rRNA quantities are just one-third of the total rRNA contents. Calculations for conversions between the mass and the number of copies of a 1,542-bp 16S rRNA were performed as follows (PCR Application Manual, 2nd ed., online version, Roche Molecular Biochemicals, Mannheim, Germany). To calculate molecular mass, the following formula was used: 1,542 bp x 340 Da/bp = 524,280 Da = 5.24 x 105 Da. To calculate the number of picomoles per microgram, we used the formula 106 pg/µg x [1 pmol/(5.24 x 105 pg)] = 1.91 pmol/µg. To calculate the number of molecules per microgram, we used the formula [(6.022 x 1011 molecules)/pmol] x 1.91 pmol/µg = 1.15 x 1012 molecules/µg = 1.15 x 109 molecules/ng.
Evaluation of Q-PCR and Q-RT-PCR.
In order to evaluate the quality of our quantification protocols, the number of ribosomes per T. denitrificans-like bacterial cell was calculated and compared to literature data. Calculations were based on respective 16S rRNA molecule numbers divided by the specific numbers of T. denitrificans-like bacterial cells. A flow diagram, comprising the several steps, analyses, and calculations for the assessment of the copy numbers for 16S rRNA molecules and genes as well as the numbers of ribosomes per cell, is given in Fig. 2.
![]() View larger version (40K): [in a new window] |
FIG. 2. Flow diagram of the culture-independent approach to assess the abundance of T. denitrificans-like bacteria. Solid lines depict those characteristics which have to be measured for our quantification approach. Dotted lines depict characteristics which have to be calculated in order to estimate the numbers of T. denitrificans-like cells, 16S rRNA molecules, and ribosomes. Detailed calculations are given in Materials and Methods.
|
|
|
|---|
Establishment of Q-PCR and Q-RT-PCR assays.
Using nucleic acids extracted from the Bx enrichments, quantification of T. denitrificans-like 16S rRNA molecules as well as 16S rRNA genes was successful without a preamplification step in a range from 8 to 50% relative abundance (Table 1). With 16S rRNA genes, the average PCR efficiency was slightly above 1 for the Com and the OST primer systems, respectively (Table 2). Similarly, with 16S rRNA, the average PCR efficiency was 0.8 for the Com as well as the OST primer system. However, this protocol was not transferable to environmental samples, where the T. denitrificans-like organism was probably present in smaller amounts. Despite the fact that melting-curve temperatures indicated satisfactory amplification products for the bacterial primer system, environmental 16S rRNA and 16S rRNA genes used as templates for real-time PCR sometimes resulted in impossible PCR efficiencies for the bacterial primer system or no amplification with the OST primer system (Fig. 3). Thus, our quantification protocol was further developed, resulting in a nested-PCR approach. Preamplification of 16S rRNA genes or 16S rcDNA with the 27F-1492R primer system before RT analyses with the Com and OST primer systems increased the molecule copy number 3.5 x 104-fold and resulted in improved PCR efficiencies for the Com as well as the OST primer system (Table 2). Moreover, Q-PCR and Q-RT-PCR was successfully applied to most of the environmental samples (Table 1). Regarding 16S rRNA analyses, the linearity of Com and OST amplification was given from quantities of 7.8 x 1010 to 1.7 x 103 molecules (Fig. 4), with a correlation coefficient of 0.999. An upper quantification limit could not be evaluated. A comparison of the specific abundances obtained by real-time PCR without a preamplification step and those from analyses including a preamplification step is demonstrated in Fig. 5. Statistically, there was no evidence for a significant variation between results gained by both quantification procedures (Student's t test; P < 0.05).
|
View this table: [in a new window] |
TABLE 2. Properties of PCR and RT-PCR and the different primer sets used for Q-PCR protocolsa
|
![]() View larger version (28K): [in a new window] |
FIG. 3. Insufficient quantification of bacterial 16S rRNA genes from sample G6 (depth, 120 m) by real-time PCR (without preamplification). (A) Environmental DNA was used as the template in quantities of 6.1 x 105 to 6.1 x 102 bacterial chromosomes for the Com primer set, representing CT values ranging from 20.6.0 to 25.1. Specific PCR products did not appear when environmental DNA was used as the template for the OST primer set. DNA contamination represented CT values ranging from 30.4 to 33.6. Norm. Fluoro., normalized fluorescence. (B) Relation between the CT and the apparent amount of 16S rRNA gene molecules by SybrGreen real-time PCR. Each point represents an amount of 16S rRNA gene molecules corresponding to the CT value with the Com primer set. Values for standard molecules are represented as blue points, and those for samples are shown as red points. The correlation coefficient of the regression line was 0.999. PCR efficiency was 1.78. (C) Melting-curve analysis of online PCR signals. Peak a shows DNA contaminants melting at 72.0 to 73.0°C; peak b shows Com products melting at 81.5 to 82.0°C. deg., degrees.
|
![]() View larger version (33K): [in a new window] |
FIG. 4. Detection limits of environmental T. denitrificans-like 16S rRNA in sample G3 (50-m depth) from the Baltic Sea by real-time PCR. (A) Sensitivity of the Com and OST primer sets in quantifying T. denitrificans-like 16S rRNA. After RT-PCR with primer 1492R and preamplification with primers 27F and 1492R, 16S rcDNA was used as the template in quantities of 7.8 x 1010 to 7.8 x 107 molecules for the Com primer set and 1.7 x 107 to 1.7 x 103 molecules for the OST primer set (values shown in italics), representing CT values ranging from 3.0 to 32.2. Norm. Fluoro., normal fluorescence. (B) Relation between the CT and the apparent amount of 16S rcDNA molecules obtained with SybrGreen real-time PCR. Each point represents an amount of 16S rcDNA molecules corresponding to the CT value with the Com and OST primer sets. The correlation coefficient of the regression line was 0.999. The PCR efficiency was 0.82.
|
![]() View larger version (19K): [in a new window] |
FIG. 5. Mean values of T. denitrificans-like 16S rRNA fractions (percent relative abundance) achieved from quantifications without (lane 1) and with (lane 2) a preamplification step. (A) Enrichment Bx21. The Student t test (P < 0.05) showed no evidence for a significant variation between results gained by using both quantification protocols (P value, 0.96; n = 15; degrees of freedom, 13). (B) Enrichment Bx24. The Student t test (P < 0.05) showed no evidence for a significant variation between results gained by both quantifications protocols (P value, 0.83; n = 8; degrees of freedom, 6).
|
|
|
|---|
Sensitivity and detection limits.
Contaminant genomic DNA is usually present in polymerase preparations and may lead to significant fluorescence in the no-template controls after about 35 cycles (1). In our experiments, this was indeed observed with the Com primer system, mostly with E. coli K-12 as a contaminant (data not shown), but usually not with the OST primer set (data not shown). This may be due to the high specificity of the T. denitrificans-like primers employed. In an attempt to overcome the contamination problem, Corless et al. (6) applied several treatments, such as UV light, DNase I enzyme treatment, or restriction endonuclease digestion, but all treatments were accompanied with an associated decrease in sensitivity. In order to increase sensitivity, we applied the preamplification approach, which enriched the target for nested Q-PCR by more than 4 orders of magnitude. This approach resulted in an observable PCR product, but it must be noted that the additional PCR step may also increase the PCR biases. Furthermore, it should also be noted that our enrichments consisting of at least 8% target DNA or RNA were successfully quantified without preamplification. However, quantification of these enrichments using a preamplification step did not yield significantly different results (Fig. 5), and community fingerprints with and without preamplification were also found to be highly similar (data not shown). Rodrigues et al. found that real-time PCR has a 6-log dynamic range of detection for DNA from pure cultures (102 to 107 cells ml1) but that it is less sensitive with soil samples (31). Moreover, their experiments suggested that the type of target gene sequence can also be an important factor in detection. They found different standard curves for the fcb-targeted gene and the 16S rRNA-targeted gene. Other factors affecting detection limits in PCR-based methods included the type and composition of the sample matrix, the type of target organism, the number and diversity of bacteria in the sample, the Taq polymerase used, and the DNA extraction protocol (22, 31). However, Teo et al. (38) suggested that a nested-PCR approach using a combination of cycle-limited conventional PCR and Q-PCR enables the detection of low-copy-number target DNA when it is present in a high background of nontarget DNA. In results similar to ours (Table 2), they found that the relative amplification efficiency of a nested-PCR approach equals that obtained with the target DNA and quantification standards. However, their experiments were based on human immunodeficiency virus type 1 DNA, and environmental samples can have high levels of chemical and genetic complexity not encountered in tissue or pure cultures (36). In general, detection limits for quantification by PCR range between 10 and 200 copies per reaction. For example, Baldwin et al. (2) detected aromatic oxygenase genes down to a detection limit of 200 copies with a SybrGreen assay, and Kolb et al. (19) detected 10 to 100 copies of the pmoA gene from methanotrophs in soil. However, quantification of DNA extracts from pure cultures of S. aureus in a LightCycler system using SybrGreen proved to be less sensitive (60 nuc gene copies/µl) than using a fluorigenic TaqMan probe (6 nuc gene copies/µl) (15). Considering the sample dilutions as well as preamplification, detection limits of our approach were much lower, being 20 copies for RNA and 0.2 copy for DNA ml of seawater1 (Table 2).
PCR efficiencies.
Another source of error which may be encountered is varied PCR efficiencies. A perfect amplification, in which each target sequence is replicated during every cycle, would have an efficiency of 1.0. Bach et al. (1) detected correlations between fluorescence and starting copy number for three peptidase gene classes of the same high quality independently of a SybrGreen or a TaqMan assay. However, the PCR efficiencies of a SybrGreen assay, a 5'-exonuclease assay, and a fluorescence resonance energy transfer assay have been compared by Newby et al. (27). They found that the SybrGreen assay usually showed lower efficiencies, which was explained by the presence of the intercalating dye that can be inhibitory to PCR (46). Comparing PCR efficiencies for the fcb-targeted gene with those for the 16S rRNA-targeted gene, Rodrigues et al. (31) suggested that the target sequence can also be an important factor resulting in different efficiencies. In our studies, various PCR efficiencies were also detected (Table 2). To overcome this problem, each measurement was verified and standardized by its own domain-specific as well as T. denitrificans-like-organism-specific standardization derived from a dilution series of a corresponding environmental target 16S rcDNA or 16S rRNA gene; i.e., each sample had its own set of standards and dilution series.
Conversion of numbers of gene copies to specific cell numbers.
Some general aspects regarding quantification by real-time PCR must be taken into consideration. The final determination of bacterial load by real-time PCR in a multispecies community will be influenced by the variation in the number of rRNA operons in a given species (7). Further, the number of replication forks is directly related to the generation time, which in turn depends on the metabolic status of the bacteria at the time of sampling (17). Not knowing the exact number of copies of 16S rRNA operons in a given target species at the time of sampling represents a limitation to the assessment of bacterial numbers by real-time PCR based on the 16S rRNA gene (25). Thus, a DNA standard representing those bacteria most likely to predominate in any given habitat is important for a more accurate determination of the total bacterial load. For this reason, our external relative standardizations were specifically prepared for every sample quantified and solely based on its respective original community DNA or RNA, and they were as close to natural conditions as possible. Moreover, our approach combining different methods still allowed an estimation of numbers of taxon-specific cells and 16S rRNA molecules.
Assessment of the number of ribosomes as an indicator of cell activity.
In general, the number of ribosomes per cell is considered to be a measure of the protein-synthesizing capacity of the cell and its growth potential or growth rate. Thus, the number of ribosomes per T. denitrificans-like cell was calculated to verify our quantifications. Calculation uncertainties because of various operon numbers or replication forks have already been discussed above. An average 16S rRNA gene copy number per cell of 2.3 is assumed for "Epsilonproteobacteria" (18). The average 16S rRNA gene copy number of our investigated microbial community is unknown. However, in comparison to other data, our calculations may support our quantification assay: E. coli, growing at a maximal growth rate under laboratory conditions, can contain up to 72,000 ribosomes per cell (8). Dependent on its activity level, the marine oligotrophic ultramicrobacterium Sphingomonas sp. strain RB2256 contained 200 to 2,000 ribosomes per cell (8). Investigating bacteria from four permanently cold marine stations, Sahm and Berninger (32) calculated numbers of 350 to 1,300 ribosomes per bacterial cell. Depending on the activity level, the number of ribosomes per T. denitrificans-like cell (± the standard deviation) ranged from 100 (±70) to 1,500 (±500). This appears to be reasonable for a marine taxon. These results indicate that our real-time PCR approach can be used for the assessment of both specific and relative abundances of uncultured marine bacteria.
This work was supported by funds from the European Commission for the AQUA-CHIP project (QLK4-2000-00764).
We are solely responsible for the content of this publication. It does not represent the opinion of the European Commission. The European Commission is not responsible for any use that might be made of data appearing therein.
Present address: IOW-Institute for Baltic Sea Research Warnemuende, Biology Section, 18119 Rostock-Warnemuende, Germany. ![]()
|
|
|---|
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»