Optimizing qPCR for the Quantification of Periodontal Pathogens in a Complex Plaque Biofilm
S.S. Kirakodu*, M. Govindaswami, M.J. Novak, J.L. Ebersole , K.F. Novak
Identifiers and Pagination:Year: 2008
First Page: 49
Last Page: 55
Publisher ID: TODENTJ-2-49
Article History:Received Date: 29/11/2007
Acceptance Date: 4/2/2008
Electronic publication date: 11/3/2008
Collection year: 2008
open-access license: This is an open access article licensed under the terms of the Creative Commons Attribution Non-Commercial License(http://creativecommons.org/licenses/by-nc/3.0/),which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.
Quantitative PCR (qPCR) has recently been used to quantify microorganisms in complex communities, including dental plaque biofilms. However, there is variability in the qPCR protocols being used. This study was designed to evaluate the validity of two of these variables with the intent of developing a more standardized qPCR protocol. The two variables evaluated were (1) the use of DNA content versus actual cell counts to estimate bacterial numbers in mixed plaque samples and (2) the effectiveness of three different universal primers versus species specific primers in amplifying specific target pathogens in these samples. Results lead to the development of a standardized protocol that was shown to be highly reproducible as demonstrated by low coefficients of variation. The results also confirmed that this standardized qPCR protocol can be used as a sensitive method for quantifying specific bacterial species in human plaque samples.
Periodontitis is initiated by polymicrobial infections and characterized by inflammatory changes in the periodontal tissues leading to loss of periodontal ligament and alveolar bone . It is widely accepted that the disease occurs as the result of an inflammatory response to a subgingival, bacterial biofilm, and is particularly related to an increase in Gram-negative anaerobes at the disease site . Numerous reports have demonstrated a close association between periodontal tissue destruction and the development of a complex of microbial species within the biofilm that includes Porphyromonas gingivalis, Tannerella forsythia, Aggregatibacter (formerly Actinobacillus) actinomycetemcomitans, Treponema denticola, Prevotella intermedia, and Fusobacterium nucleatum [3-9].
Quantitative detection of these pathogens in clinical specimens has contributed to a better understanding of the crucial host-bacterial interactions that occur in periodontitis. In the past, cultivation methods have been considered the gold standard for describing the predominant microorganisms in dental biofilms [5,10]. Limitations of this method including difficulty in growing many of the fastidious microorganisms, inherent inaccuracy of counting methodologies, and high costs in terms of supplies and technical personnel, have pushed investigators to develop more sensitive, accurate, quantitative, and cost effective molecular techniques for enumerating bacteria in complex biofilms . Approaches currently being used include DNA-DNA checkerboard hybridization [11,12], standard polymerase chain reaction (PCR) [13-15], and, more recently, real-time quantitative polymerase chain reaction (qPCR) [8,16-23]. The checkerboard hybridization technique routinely can detect 104 cells but conditions can be altered to detect 103 cells . It is generally used to quantify ranges of bacteria (<105,105-106 and >106 copies) . Standard PCR lacks the ability for precise quantification since only an end-point determination can be analyzed. In contrast, qPCR provides precise counts through direct monitoring of the increasing amount of PCR product throughout the enzymatic assay. Thus, it appears that qPCR is the most sensitive of these available methods, with detection limits estimated at 102 genome copies .
Initial evaluation of various reports using qPCR for the enumeration of oral bacteria revealed substantial variability in methodologies. These include variations in protocols used to generate the standard curves for qPCR and differences in the universal primer sets used to amplify 16S rDNA [8,16, 17,19,20,22,23,26]. There are no reports in the literature that evaluate and compare how these different methodologies affect the outcomes of bacterial quantification. Given that qPCR is becoming an accepted technology for the quantitative analysis of bacteria from mixed biofilm samples, this study was designed to compare these variables with the intent of developing a more standardized qPCR protocol. Specifically, the aims of this study were: (i) to compare the use of actual cell counts versus a standard formula based on total DNA to estimate the total number of bacteria present in bacterial cultures that are being used to generate qPCR standard curves; (ii) to determine the effectiveness of three different universal primer sequences versus species specific primers in amplifying 16S rDNA obtained from cultures of the seven target bacteria; (iii) to determine the reproducibility of qPCR using LightCycler technology in amplifying total bacteria and seven individual bacteria, and, (iv) to use a standardized qPCR procedure to quantify the total bacterial load in samples of oral biofilms, and the numbers of the seven target periodontal bacteria in the complex microbial samples.
MATERIALS AND METHODS
Bacterial Strains and Growth Conditions
Bacterial strains used in this study were purchased from ATCC (ATCC, Rockville, MD) and cultured in a Coy anaerobic chamber in an atmosphere of 85% N2, 5% CO2 and 10% H2 at 37°C. F. nucleatum ATCC 25586 and A. actinomycetemcomitans ATCC 700685 were cultured in 3% Trypticase Soy broth supplemented with 0.6% yeast extract. P. gingivalis ATCC BAA-308 and P. intermedia ATCC25611 were cultured in 2.1% Mycoplasma broth supplemented with 0.5% hemin and 0.1% menadione. T. denticola ATCC 35405 was cultured in GM-1 broth . T. forsythia ATCC 43037 was cultured in 4 % Tryptic Soy broth supplemented with 0.5% hemin and 0.001% N-acetylmuramic acid. Campylobacter rectus ATCC 33238 was cultured in 2% Mycoplasma broth supplemented with 0.2% ammonium formate, 0.3% disodium fumarate, 0.5% hemin, pH adjusted to 7.3 and 2% horse serum (not heat inactivated). All growth media were pre-reduced in the anaerobic chamber for 24-36 hrs prior to inoculation with bacteria. Following growth to late-log phase an aliquot from each culture was removed, diluted 10-fold, and used to estimate the total bacterial count (Petroff-Hausser counter). Each culture was also Gram stained to check for purity and cell morphology. Duplicate, 1 ml aliquots of each culture were dispensed in eppendorf tubes and the cells pelleted by centrifugation at 10,000g for 10 min using a microcentrifuge. The spent medium was discarded, leaving 25µl for resuspension of the cell pellet. These aliquots, containing a known number of bacteria, were then stored at -80°C until they were used for DNA isolation and generation of standard curves for qPCR.
Collection and Storage of Bacterial Plaque Samples
After obtaining appropriate institutional review board approval, two individual human subgingival bacterial plaque samples were collected at least four months apart from five patients with periodontitis (defined as the presence of 4 or more teeth with a probing depth of at least 4mm and clinical attachment loss of at least 2mm, and bleeding on probing at 35% or more of tooth sites). Samples were taken from the four deepest sites of the dentition that bled on probing and were pooled into sterile phosphate buffered saline (PBS, pH 7.2). The samples were stored at -80oC until processed for qPCR.
DNA was isolated from the frozen, pure cultures of each individual target bacterium using the MasterpureTM kit from Epicenter (Epicenter, Madison WI) and quantified by LightCycler green dye (Idaho Technologies) using the LightCycler 2.0 instrument (Roche Applied Science).
Bacterial plaque samples in sterile PBS were centrifuged at 10,000g to pellet the bacterial cells. DNA was extracted from the pelleted cells using the MasterpureTM kit and quantified as described above. DNA from each sample was resuspended in 100 μl of sterile water.
Comparing Actual Cell Counts to Determination by a Standard Formula
The generation of standard curves for the enumeration of test bacteria is an essential component of qPCR. Two methods are generally used to accomplish this – actual cell counting or calculation based on a standard formula. To compare these two methods, the seven bacterial strains were grown in pure culture to late log phase, actual cell counts evaluated using a Petroff-Hausser counting chamber, and DNA isolated. The cell counts for each test organism were then compared with the estimated number of copies of each microorganism in the isolated DNA using the following formula: (Avogadro constant X amount of DNA in µg/µl) / (genome size X mw), where mw is the molecular weight per base pair or nucleotide which is 660 Da [28,29].
Quantitative (Real-Time) Polymerase Chain Reaction
Real-time PCR was performed using a LightCycler 2.0. Each PCR was performed in a total volume of 20 µl containing 2 µl of 10X LightCycler FastStart DNA Master SYBR Green I, 0.5 µM each of HPLC purified forward and reverse primers (Table1), 4mM MgCl2 and 1μl of template DNA. Real-time PCR was carried out with an initial incubation of 10 min at 95ºC followed by 45 cycles consisting of denaturing at 95ºC for 10 sec; annealing (for annealing temperature see Table1) for 5 sec followed by amplification at 72ºC for X sec, where X is the amplification time in seconds (product size/25). At the end of each cycle, fluorescent products were detected by the LightCycler 2.0. After amplification, a melting curve analysis was performed to determine the specificity of the PCR products. For melting curve analysis PCR products were incubated for 15 sec at 5ºC below the annealing temperature for the respective primers and the temperature was increased to 95ºC with a ramp rate of 0.1/sec. DNA strands separate or melt depending on the sequences, the length of the strand and the GC content of the strand. Because melting temperatures vary according to these differences, melting temperature profiles were used to identify and genotype the PCR products.
16S rDNA Bacterial Primers and Annealing Temperatures Used in Real Time PCR and Expected Product Size
|Organism||Sequence||Size (bp)||Annealing Temperature||Ref|
Comparison of Actual Cell Counts Versus a Standard Formula* in Estimating Total Bacteria Present in Pure Cultures
|Bacteria||Cell Counts (x 107)(CT)a||Total DNA (ng)b||Genome Size (MB)c||Cell Counts (x 107)(SF)d|
* Standard formula: (Avogadro constant X amount of DNA in µg/µl) / (genome size X mw), where mw is the molecular weight per base pairs or nucleotide which is 660 Da); a: actual microscopic bacterial counts using a Petroff-Hausser counting chamber; b: total DNA isolated from each bacterial culture; c: Genome size based on published genome se-quence for each oral bacterial species; ±: estimated genome size based on published genome sequences for Campylobacter jejuni RM1221 and Campylobacter jejuni subsp. jejuni NCTC 11168 (http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi);d: calculated cell numbers using the standard formula (SF; 5,11).
Crossing Point Values Using Species Specific and Universal Primers
|Species||DNA (ng)||Species Specific Primer||Universal Primer 1||Universal Primer 2||Universal Primer 3|
The Coefficient of Variation (CV) from qPCR Runs for Different Primer Sets
* Aggregatibacter actinomycetemcomitans (Aa); Porphyromonas gingivalis (Pg), Tannerella forsythia (Tf); Fusobacterium nucleatum (Fn); Prevotella intermedia (Pi); Campylobacter rectus (Cr); Treponema denticola (Td).
Mean denotes mean CV% from 22 samples tested in triplicate for all primers except Aa, which was derived from the 12 samples that contained measureable Aa. Median denotes the median CV% for the same samples. Correl denotes the rank correlation of the CV% with the total number of bacteria (Universal) or the number of the individual species. A minimum correlation coefficient of ±0.4227 is required for a p<0.05 level of significance.
qPCR Enumeration of Total Bacteria and Specific Bacterial Species in 100 µl of Plaque DNA
|Total Bacteria (x108)||Aa (x104)||Pg (x105)||Td (x106)||Tf (x106)||Fn (x106)||Pi (x105)||Cr (x105)||Sum of 7 (x107)|
* Values are expressed as absolute counts and percentage of total bacteria (in parentheses), Universal primer 1 was used for enumeration of total bacteria and species-specific primers for A. actinomycetemcomitans (Aa), P. gingivalis (Pg), T. denticola (Td), T. forsythia (Tf), F. nucleatum (Fn), P. intermedia (Pi) and C. rectus (Cr).
Quantification analysis was performed using LightCycler 4.0 software. Quantification analysis uses sample “crossing points” (CP) to determine the presence and the concentration of the target DNA in known and unknown samples after amplification. The quantification module only considers fluorescence values measured in the exponentially growing phase of the PCR amplification process. Standard curves were generated for all the target bacteria using DNA from pure cultures and species specific DNA primers. DNA isolated from a pure culture of P. gingivalis was used in generating a standard curve for the universal primers [3,26]. Ten fold serial dilutions of genomic DNA (DNA from 107, 106, 105, 104, 103, 102 and 101 cells), from pure bacterial cultures provided seven data points for standard curve generation. In each standard curve, the concentration of standard sample (known amount of DNA/DNA from known number of bacteria) was plotted against its crossing points. Quantification of the individual target bacteria and total bacteria from the experimental samples were calculated using the standard curves.
Comparison of Actual Cell Counts to Estimation by Standard Formula
Results of experiments designedto compare the use of actual cell counts versus a theoretical formula in estimating the total number of bacteria present in pure cultures for development of standard curves revealed that in every case, estimated counts were higher based on use of the formula than those obtained from actual microscopic counts (Table2). Results using the two methods were most similar for T. forsythia and P. gingivalis although both methods yielded results in the same log value. It should be noted that the species with the greatest percentage deviation, C. rectus, required an estimation of genome size based upon the average genome sizes for sequenced Campylobacter species (Campylobacter jejuni RM1221 and Campylobacter jejuni subsp. jejuni NCTC 11168; http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi).
Comparing the Effectiveness of Different Universal Primers
Crossing point (CP) values were used to evaluate the effectiveness of three different, published, universal primer sequences in amplifying DNA obtained from pure cultures of the seven target periodontal microorganisms by comparing them with species specific primers. These species specific primers were specifically selected from published articles and have been validated and used extensively in amplifying oral bacteria from plaque samples [8, 15, 30,31,32].
Results from these experiments revealed that Universal primer set 1 (Table1) consistently yielded CP values most closely aligned with the species-specific primers across all seven target species (Table3). Universal primer set 2 yielded CP values that were lower than that of Universal primer set 1 for all species except T. denticola. In general, Universal primer set 3 did not perform well, yielding CP values that were significantly higher than those obtained with the species-specific primers. One exception was A. actinomycetemcomitans for which the CP value was lower than that found with the species-specific primer.
Reproducibility of qPCR
Reproducibility of the qPCR for the assessment of total microbial counts as well as counts of each of the seven target species of bacteria was determined using 22 individual plaque samples. Each sample, with each primer set, was evaluated in triplicate. The coefficient of variation (CV) for each sample was calculated and the median for the 22 samples was determined. The results showed a 6.9% median for the Universal primer and the median CV for the seven bacteria ranged from 4.8% to 10.4% (Table4). Thus, the qPCR quantification was very reproducible regardless of the characteristics of the individual plaque samples. Additionally, we determined if there was a relationship between the CV and the number of any particular bacterial species in the sample. The results showed no significant correlation for the CV, irrespective of whether the sample had high or low levels of a particular species.
Enumeration of Bacterial Species in a Complex Biofilm
Using the standardized protocol described above we quantified both total bacteria and individual counts of the seven target microorganisms in 10 individual mixed plaque biofilm samples (Table5). Total bacteria and specific bacteria from the plaque samples were enumerated by using the standard curves generated from the actual bacterial cell counts. After factoring in the total volume of each sample, we determined that the total bacterial load in all samples, except sample 7, averaged 3.97 X 108 cells. Total bacterial load in sample 7 was a full log higher, at 1.07 X 109 cells. In contrast, counts of the seven specific bacteria where highly variable among the samples. This variability was most notable in A. actinomycetemcomitans, where counts ranged from undetectable (samples 1, 2, 5 and 6) to 7.97 X 106. This highest value was found in sample 7, which also had the highest total bacterial count. Levels of P. gingivalis also were highly variable, ranging from 6.1 X 102 (sample 2) to 8.67 X 107 (sample 7). There was less variability among counts of other species, with all samples remaining within one log of each other. The sum of all seven target microorganisms averaged 4.65 X 107 in the ten samples. Sample 7 had higher total counts of these seven bacteria, with a total of 1.50 X 108. The data also were expressed as the percentages of these seven species in the total mixed plaque samples. These ranged from 0% (A. actinomycetemcomitans and P. gingivalis) to 8.6% (F. nucleatum). Assessment of the sum of all seven as a percentage of total counts revealed that these specific microorganisms comprised 7.5%-28.1% of the total microbial flora.
This study focused on details of a reproducible methodology to utilize modern qPCR technology to quantify the numbers and distribution of individual bacterial species within a complex microbial ecology. In addition to the differences in methodology described above, there are different systems available for carrying out qPCR reactions. For example, both TaqMan probes and SYBR Green have been used in qPCR for quantification and identification of oral bacteria. Maeda et al.  have reported that there are no significant differences between the TaqMan and SYBR Green chemistry in their specificity and sensitivity. However, the SYBR Green assay involves fewer manipulations than the TaqMan assay and therefore appears to be the less complex of the two. We therefore chose to use the SYBR Green technology in the LightCycler 2.0 system in this study.
A unique feature of this study was to compare the use of actual bacterial cell counts to the use of a standard formula to estimate the total number of a particular bacterial species present in a microbial sample.Kuboniwaet al.  reported that the yield of purified genomic DNA from pure cultures of P. gingivalis, where the cell number had been determined by colony counting, was considerably less than the calculated genomic weight. This suggested that, as one might anticipate, estimating total bacterial counts by cell counts may not be suitable to provide accurate enumeration of bacterial numbers in samples.However, both methods have been used in the development of standard curves for qPCR [5,7,33,35], thus, a more detailed comparison was desirable to evaluate validity/accuracy of these methods. The results demonstrated that, with the seven oral bacteria used in this study, differences were routinely noted between the two methods, with absolute counts being consistently lower than those obtained with the standard formula. However, the magnitude of these differences may not be significant given that in all cases, the log values of the counts were identical using either method (Table2). Given the technical limitations in microscopic counting and the likely differences in absolute genome size between bacterial strains used in various studies versus the genome sequencing strains, 20% variability between the methods would be expected to be within the margin of error.
Moreover, in clinical applications, these techniques will be used to assess differences in counts of periodontal microorganisms either longitudinally within the same patients or in cross-sectional designs between groups of individuals. Consequently, log-fold changes would generally be considered the basis for describing changes over time in patients and/or among groups. This magnitude of change is consistent with other methodologies designed to assess differences in microbial counts in clinical studies. Therefore, our data suggests that either method can be used to assess the total numbers of individual bacterial species within a complex sample. In those cases where the genome size may not be available for the bacterium of interest, eg.C. rectus, absolute bacterial counts may provide the optimal method.
Each of the universal primer sets used in this study has been used in qPCR analysis of human bacterial samples [16,26,35]. Universal Primer 1 has been used in studies of oral treponemes and P. gingivalis . Nadakarni et.al;  used Universal Primer 2 and compared its efficiency in amplifying DNA from 34 species of bacteria. They reported that the anaerobic bacteria estimated by this universal primer set in carious dentine was 40-fold greater than the total bacterial load detected by culture methods, demonstrating the utility of real-time PCR. Finally, Universal Primer 3 has been reported as a useful tool in studies of a rapid method for the quantification A. actinomycetemcomitans and P. gingivalis in saliva and subgingival plaque . However, this is the first individual study to evaluate the effectiveness of these different universal primers, as compared with species specific primers, to amplify 16S rDNA from oral bacteria. The results clearly showed variability in the consistency of different primer sets to accurately estimate individual bacteria that would be components of complex oral microbial samples. An extreme example was Universal primer 3 which demonstrated exceptionally variable CP values across the 7 oral bacteria when compared to the species-specific primers. However, it must be noted that Universal primer 3 was originally designed for use in TaqMan probe analysis and its small amplicon size may not be suitable for the SYBR method used in this study. As importantly, we determined that Universal primer 1 provided consistent CP values compared to the specific primer sets for these seven oral bacteria.
Reproducibility of qPCR reactions and methodology and the need to provide replicate sample analyses may be of concern when evaluating large numbers of samples from multiple patients. In an evaluation, universal and species specific primers were tested by triplicate qPCR runs for 22 plaque samples and we demonstrated that qPCR is highly reproducible (Table 4). These results are in accordance with the published reports of Betzl et al.  and indicate that replicate qPCR runs are not necessary to accurately assess microbial populations in mixed bacterial samples.
Using the standardized methods described in this study, we applied a qPCR approach towards quantifying the seven target bacteria in human subgingival biofilm samples from patients with periodontitis (Table5). Total bacteria in the plaque sample were enumerated using universal primer 1 and individual bacterial counts using the species specific primers (Table1). As expected, there were differences between samples in individual bacterial counts and the individual proportions of each species within the samples. These differences are reflective of the inherent variability in the microbial ecology among individuals. However, because samples were obtained from periodontitis sites, it was anticipated that one or more of the target pathogens would be present in each sample. The results demonstrated that all of the test samples harbored a minimum of five of these pathogens. In addition, the sum of the target bacteria represented a consistent proportion of the total plaque sample ranging from 5% to 28%. All of these results are consistent with previous studies showing that periodontal pockets typically harbor these bacteria within the ranges found in our study . The general overall agreement in the distribution and quantities of the bacteria supported the validity and utility of this qPCR approach for evaluating microbial species and changes in subgingival biofilm samples as a function of disease progression or interventional therapy.
This work was supported by U01 AG021406 from the National Institute on Aging and U01 DE014338 from the National Institute for Dental and Craniofacial Research.