Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2004, p. 7545-7549, Vol. 70, No. 12
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.12.7545-7549.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
| SHORT REPORT |
Kankyo Engineering Co., Ltd., Kisarazu,1 Institute for Biological Resources and Functions, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba,2 Department of Molecular Engineering, Tokyo Institute of Technology, Nagatsuta, Midori-ku, Yokohama,3 Food Engineering Division, National Food Research Institute, Tsukuba, Ibaraki, Japan4
Received 26 December 2003/ Accepted 21 July 2004
ABSTRACT
We reevaluated the bias toward a 1:1 ratio of products in multitemplate PCR used in ecological studies and showed that the template reannealing at the annealing step would not cause the bias; however, the preferential homoduplex formation during temperature decrease from denaturation to annealing step would cause the bias.
Multitemplate PCR has been frequently used in microbial ecology to determine the microbial community structure in environmental samples. However, it has been reported that the analysis of microbial communities based on multitemplate PCR can contain wrong information on the abundance and diversity of genes (3, 6, 21). This problem is classified into two categories: bias and artifact (6).
Artifact problems can be caused by the formation of a heteroduplex and a chimera. A heteroduplex is formed in PCR by the cross-hybridization of heterologous sequences (5, 16, 20) and can give data for nonexistent genes by further analysis using cloning and sequencing (22), denaturing gradient gel electrophoresis (8, 18, 23), terminal restriction fragment length polymorphism (T-RFLP) (12), and so on (1). Intrastrand reannealing of the PCR product can also cause additional gene diversity in T-RFLP analysis (2). A chimera is formed from an incompletely extended primer (13, 17) and template switching (11, 14) during multitemplate PCR and gives artificial gene diversity.
Bias can be caused in multitemplate PCR by differences in primer binding energy (4, 9, 15) and reannealing of templates (10, 19). Reannealing of templates has been reported to cause a strong bias toward 1:1 mixtures of genes in final PCR products, regardless of the initial ratio of the templates (10, 19). The mechanism of this bias has been explained as the amplification rate of more abundant PCR products declining faster with the reannealing of PCR products than that of less abundant products in the same tube; hence, the final product concentrations are biased toward a 1:1 ratio independent of the initial template concentrations (10, 19). However, in our preliminary experiments, we did not observe this 1:1 bias in several cases. Suzuki and Giovannoni (19) reported that the PCR with low amplification efficiencies and low PCR product yields did not show bias. This result was true in our case, but even with high PCR product yields we did not observe the bias in some cases.
In this study, we reevaluated the 1:1 bias by using a real-time quantitative PCR method to monitor amplification efficiencies and the amount of PCR product. We have found that the reannealing of products at the annealing step does not cause the bias, but the preferential homoduplex formation by reannealing of the products during the decrease in temperature from the denaturation to the annealing step does cause the bias.
We prepared fluorescently labeled PCR products of bacterial 16S rRNA genes as templates. The fluorescently labeled templates were required for T-RFLP analysis, described below. First, three bacterial strains, Microbacterium imperiale DSM 20530, Pseudomonas fluorescens DSM 50108, and Agromyces mediolanus DSM 20152, were purchased from Deutsche Sammlung von Mikroorganismen und Zellkulturen (Braunschweig, Germany), and they were cultured overnight at 30°C in test tubes containing 5 ml of a medium of the following composition: 1.0% (wt/vol) casein peptone, 0.5% (wt/vol) yeast extract, 0.5% (wt/vol) glucose, and 0.5% (wt/vol) NaCl. The genomic DNA of each microbe was extracted from the cultured cells by using a DNeasy Tissue kit (QIAGEN GmbH, Hilden, Germany) according to the manufacturer's instructions. The 16S rRNA genes were then amplified by PCR from the genomic DNA with the primer pair EC27f-BODIPY FL (5'-CAGAGTTTGATCCTGGCTCAG-3') and E1389r (5'-ACGGGCGGTGTGTACAAG-3'). All primers were supplied by Espec Oligo Service Corp (Tsukuba, Japan). The primer EC27f-BODIPY FL was designed by adding a BODIPY FL-modified cytosine to the 5' end of primer 27f. The PCR products were purified with Microcon PCR filters (Millipore Corp., Bedford, Mass.) and were used as templates in this study. The PCR products from M. imperiale DSM 20530, P. fluorescens DSM 50108, and A. mediolanus DSM 20152 were designated Mim, Pfl, and Ame, respectively. The concentrations of the purified PCR products were measured with a Pico Green dsDNA Quantification kit (Molecular Probes, Inc.). Accession numbers of 16S rDNA sequences of M. imperiale DSM 20530, P. fluorescens DSM 50108, and A. mediolanus DSM 20152 in GenBank were X77442, AF094728, and X77449, respectively.
In PCR bias experiments, a quenching primer PCR method (7) was used for real-time monitoring of the amount of PCR product. This method uses a primer containing a BODIPY FL-modified cytosine at its 5' end (Fig. 1). When such a primer is hybridized with target DNA, the fluorescence of BODIPY FL is quenched by guanine in the target complementary to modified cytosine, and the quenching rate is proportional to the amount of target DNA (7). Therefore, the amount of PCR product is calculated with the following formula: Cproduct = (R/Rmax) x Cprimer, where Cproduct (nanomolar concentration) is the concentration of PCR products, Cprimer (nanomolar concentration) is the initial concentration of the modified primer, R is the measured quenching rate, and Rmax is the maximum quenching rate when all primers are hybridized with the target.
![]() View larger version (17K): [in a new window] |
FIG. 1. Fluorescence quenching of a BODIPY FL-labeled primer.
|
In PCR bias experiments, the quenching primer PCR was carried out with a LightCycler System with an initial denaturation at 95°C for 120 s followed by a cycle of denaturation at 95°C, primer annealing at 56°C, and extension at 72°C. The duration of each step and the ramp time are shown in Table 1. Templates Pfl and Ame were mixed at a ratio of 9:1, 3:1, 1:1, 1:3, or 1:9. Primers EC27f-BODIPY FL and E350r (5'-CAATATTCCCCACTGCTGCC-3') were used in the experiments. The reaction mixture (20 µl) contained 107 molecules of mixed templates, 500 nM each primer, 2 mM MgCl2, 0.25 mg of bovine serum albumin/ml, 2 µl of reaction buffer (10x rTaq buffer), 1.6 µl of deoxynucleoside triphosphate mixture (2.5 mM each), 0.5 U of DNA polymerase (rTaq), and 0.32 µl of TaqStart antibody (Clontech Laboratories Inc., Palo Alto, Calif.). The fluorescence was measured after denaturation at 95°C and extension at 72°C in each cycle. The fluorescence quenching rate was calculated with the following formula: Rn = [1 (F72, n/F95, n)/(F72, 10/F95, 10)] x 100, where Rn is the quenching rate in cycle n, F72, n is the fluorescence intensity at 72°C in cycle n, and F95, n is the fluorescence intensity at 95°C in cycle n. The amount of PCR product at each cycle was calculated with this quenching rate. The PCR was stopped at 40 cycles if the reaction had reached a plateau phase, because the bias has been reported to be observed in late PCR cycles (19). Otherwise, the reaction was stopped at 80 cycles, at which time the reaction had already reached a plateau phase. In separate experiments, PCR was stopped at 15, 17, 20, 24, 29, and 40 cycles, and the amount of PCR product was determined by a Pico Green dsDNA Quantification kit. The value was compared to that calculated from the quenching rate, and as a result they were almost identical (data not shown), indicating that we can know the amount of product by the quenching rate.
|
View this table: [in a new window] |
TABLE 1. PCR conditions
|
Before T-RFLP analysis we examined the accuracy of this analysis, because in an ABI 310, different fragments will be differentially loaded into capillaries by electroinjection, and this may affect the quantification. We prepared an equimolar mixture of six BODIPY FL-labeled fragments having various lengths (22 to 531 bp) and applied it to the ABI 310. As a result, the longer fragment showed a greater peak area, but the peak height was nearly the same for all fragments (data not shown), as has been reported by other researchers (9). Therefore, we used the peak height for the calculation of PCR product ratios and verified the linear relationship between the molar ratio and the peak height ratio by using a mixture of the fragments at various ratios (data not shown).
To reevaluate a PCR bias toward 1:1, we performed PCR with templates Pfl and Ame in various conditions (shown in Table 1). Condition A was similar to that in which Suzuki and Giovannoni observed PCR bias (19). The temperature ramp time in condition A is the usual time that is used when PCR is carried out in a conventional Peltier cycler with an aluminum block and plastic tubes. In condition A, the amount of PCR products increased until the 26th cycle, at which time it showed stable values (Fig. 2). PCR was carried out with several replicates and was stopped after the 13th, 19th, or 40th cycle, and then the product ratio was determined by T-RFLP (Fig. 3A). The product ratio in the 13th cycle was the same as the initial ratio. In the 19th cycle, the product ratio was biased toward a 1:1 ratio, and in the 40th cycle, the product ratio was even more biased. Thus, the PCR bias to a 1:1 ratio appeared in the late PCR cycles, as has been shown by other researchers (10, 19). The same experiment was carried out with plastic tubes and a conventional Peltier cycler. As a result, almost the same bias was observed after the 40th cycle (data not shown), indicating that the difference in reaction tubes and heating-cooling systems did not affect PCR bias.
![]() View larger version (21K): [in a new window] |
FIG. 2. Real-time monitoring of the amount of PCR products in condition A. PCR was carried out with several replicates and was stopped after the 13th (open square), 19th (open circle), or 40th cycle (open triangle).
|
![]() View larger version (33K): [in a new window] |
FIG. 3. Relation between initial template ratio and product ratio after PCR under various conditions. (A to G) PCR was carried out in conditions A to G, which are described in Table 1. The template for PCR was a mixture of the templates Ami and Pfl, and the ratio of Pfl to Ami is shown. (A) PCR was stopped after the 13th (open square), 19th (open circle), and 40th (open triangle) cycle, and the product ratio was analyzed.
|
Our results indicate that the ramp time from the denaturation step to the annealing step was one of the factors behind the bias. The mechanism of this bias can be explained as follows: during the decrease in temperature from the denaturation step to the annealing step, three kinds of duplexes begin to form in the order of homoduplexes, heteroduplexes, and duplexes between primers and templates, because the thermodynamic stability of these duplexes is in this order. Here we divided the temperature area into areas 1, 2, and 3 (Fig. 4). We defined each area as follows: in area 1, homoduplexes are preferentially formed; in area 2, both homoduplexes and heteroduplexes are formed at similar rates of possibility; in area 3, three kinds of duplexes are formed. In PCR condition C (Table 1), the ramp time from the denaturation step to the annealing step was 20 times longer than that in condition B. This means that the reaction mixture stayed 20 times longer in area 1 in condition C, and therefore higher amounts of homoduplexes were formed in condition C. The homoduplex formation occurs more frequently within more abundant PCR products than within less abundant products, and therefore the decreasing rate of single strands is be higher in more abundant products, as suggested by Suzuki and Giovannoni (19). This leads to the bias toward a 1:1 ratio. In areas 2 and 3, both homoduplexes and heteroduplexes are formed at similar rates of possibility, and the decreasing rate of single strands is the same for all PCR products. Therefore, it is not likely that homoduplex and heteroduplex formation in areas 2 and 3 cause the bias.
![]() View larger version (29K): [in a new window] |
FIG. 4. Conceptional diagram of temperature areas 1, 2, and 3.
|
![]() View larger version (28K): [in a new window] |
FIG. 5. Relation between initial template ratio and product ratio when the pair of templates were Ame plus Pfl (open triangle), Min plus Pfl (open circle), and Ame plus Min (open square).
|
Our experimental results showed that reannealing of PCR products at the annealing step would not cause the bias. Rather, preferential homoduplex formation due to reannealing of PCR products during the decrease in temperature from the denaturation step to the annealing step caused this bias.
Our experiments suggest that the way to minimize PCR bias toward a 1:1 ratio is to set the fastest ramp rate from the denaturation step to the annealing step, using PCR machines with a fast ramp rate. However, this may enhance the formation of heteroduplexes when PCR reaches plateau phase (16). The heteroduplexes will not cause bias in PCR but will give artificial gene diversity by further analysis (6). In such cases, elimination of heteroduplexes after PCR will be needed (16, 20).
The elucidation of the mechanism of PCR bias and artifact would be very important for improving the accuracy of microbial community analysis based on multitemplate PCR. Our findings will contribute to reducing PCR bias toward a 1:1 ratio.
REFERENCES
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»