Previous Article | Next Article ![]()
Journal of Clinical Microbiology, November 2005, p. 5753-5759, Vol. 43, No. 11
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.11.5753-5759.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
University of Pittsburgh, School of Dental Medicine, Pittsburgh, Pennsylvania,1 Harvard University, School of Dental Medicine, Boston, Massachusetts,2 University of Pittsburgh, Hillman Cancer Center, Pittsburgh, Pennsylvania,3 National Institute of Dental and Craniofacial Research, NIH, Division of Intramural Clinical Research, Bethesda, Maryland,4 Forsyth Institute, Boston, Massachusetts,5 University of Oslo, Institute of Oral Biology, Oslo, Norway,6 Twins Institute for Genetics Research, Montes Claros, Brazil,7 UNIMONTES, Montes Claros, Brazil8
Received 5 July 2005/ Returned for modification 15 August 2005/ Accepted 29 August 2005
|
|
|---|
|
|
|---|
In previous studies, conventional culturing methods have been used to show that well-known species, such as Streptococcus mutans and Lactobacillus spp., are associated with dental caries (13). These species have been reported as potential contributors to caries onset and development. More recently, advanced molecular methods of bacterial identification, such as PCR techniques and 16S rRNA gene sequencing analysis, have become available and have revealed that the bacterial involvement in the development of dental caries is more complex than previously believed (2).
The aim of this study was to use molecular identification methods, such as 16S rRNA gene sequence and reverse-capture checkerboard hybridization, for identification of the bacteria associated with dental caries and health in infants and children. We compared the bacteria found in biofilms of caries-active children with different degrees of lesion severity with those found in biofilms of caries-free children. In addition, we compared biofilms of healthy surfaces of caries-active subjects with biofilms of healthy surfaces of caries-free subjects. A simple univariate test was used to determine the overabundance and underabundance of bacterial species in the diseased and in the healthy groups. Features identified with this univariate test were used to construct a probabilistic disease prediction model. With proper machine learning-based evaluation, we found that this model was successful in utilizing biofilm bacterial risk indicators to predict disease and health. Our modeling approach splits a data set into two groups of samples, a training and a test set. Because we evaluated the model by performing learning (classifier construction) in the training set and evaluated in the test set, the performance of the prediction model could be statistically validated and, thus, could be generalized.
|
|
|---|
Study design and sampling. Two examiners conducted dental caries examinations according to National Institute of Dental and Craniofacial Research criteria (3) modified to distinguish caries lesions with a chalky whitish/yellowish opaque appearance, without clinically detectable loss of substance (white spot lesions), from cavitated carious lesions. Interproximal surface caries were assessed using digital imaging fiber-optic transillumination (DIFOTI, Irvington, NY). About 20% of the total cohort was caries free, 25% had low to moderate levels of dental caries rates, and 55% had high to rampant caries rates.
Plaque samples were collected in the morning, and the children were asked not to brush their teeth or eat the night before the exam. Supragingival dental plaque samples were collected from the entire cohort over a period of 4 weeks. The caries-free twins had pooled plaque samples collected from three healthy surfaces, including anterior and posterior teeth. Caries-active subjects had plaque samples collected separately from four types of surfaces: surface of intact enamel (site 1), surface of white spot lesions (site 2), surface of initial enamel lesions (site 3), and excavated plaque from deep dentinal lesions (site 4) (Fig. 1). All caries-active subjects provided three to four sites of plaque collected separately from different teeth according to the severity of disease. For intact enamel and white spot lesions, plaque was collected by swiping the tooth surface with a Stimudent (Johnson & Johnson, New York, NY), whereas plaque from cavitated lesions was collected by means of a small Gracey curette (1-2; Hu-Friedy, Chicago, IL).
![]() View larger version (22K): [in a new window] |
FIG. 1. Study design. Model 1 compared site 1 (biofilm of surfaces of intact enamel) in the caries-free group with all sites combined in the caries-active group: site 1 (biofilm of surfaces of intact enamel), site 2 (white spot lesions), site 3 (dentinal lesions), and site 4 (deep dentinal lesions). Model 2 compared site 1 (biofilm present on surfaces of intact enamel) from the caries-free group to only site 1 (biofilm present on surfaces of intact enamel) from the caries-active group.
|
Isolation of bacterial DNA. Samples were placed inside a sterile microcentrifuge tube already containing 50 µl of TE buffer (50 mM Tris-HCl, pH 7.6, 1 mM EDTA) and immediately processed for DNA extraction. Five microliters of 5% Tween 20 and 1 µl of 10-mg/ml proteinase K (100 mg; V3021; Promega) were added to each sample. Samples were incubated at 55°C for 2 h. The proteinase K was inactivated by heating the samples at 95°C for 10 min. Samples were frozen at 20°C and transported in dry ice from Brazil to the United States. Samples were frozen at 80°C for further analysis.
Cloning, sequencing, and clonal analysis. Initially, four random samples from different types of lesions (intact healthy surfaces, white spot lesions, surfaces of initial enamel lesions, and surfaces of deep dentinal lesions) were selected for clonal analysis for initial characterization of microbial profiles in this cohort and to further determine the target species to be used for checkerboard analysis. Cloning procedures, which included 16S rRNA sequencing and data analysis, were performed with the TOPO TA cloning kit (Invitrogen, San Diego, CA) as previously described (10). A total of 200 clones with the insert of the correct size of approximately 1,500 bases were analyzed (100 per sample) and compared to known 16S rRNA gene sequences in the Ribosomal Database Project (7).
Amplification of bacterial DNA. The 16S rRNA genes were amplified under standard conditions using a universal forward primer and a universal reverse primer (10). PCR was performed in thin-walled tubes with a Perkin-Elmer 9700 thermocycler. One microliter of the DNA template was added to 22.5 µl of a PCR SuperMix (Invitrogen, Carlsbad, CA). In a hot start protocol, samples were preheated at 94°C for 4 min, followed by amplification under the following conditions: denaturation at 94°C for 45 s, annealing at 60°C 45 s, and elongation for 1.5 min with an additional 5 s for each cycle. A total of 30 cycles were performed and were then followed by a final elongation step at 72°C for 10 min. The results of PCR amplification were examined by electrophoresis in 1% agarose gel. DNA was stained with ethidium bromide and visualized under short-wavelength UV light. For checkerboard hybridization, the resulting PCR product from each sample was used to run three sets of PCR amplification using a forward primer labeled at the 5' end with digoxigenin and a universally conserved reverse primer. This approach was used based on the fact that for each sample we needed to run three different checkerboards (we could only fit 30 probes on each checkerboard). Using the same PCR product for the subsequent labeled PCR reduced the variation in the samples.
Checkerboard hybridization. The reverse-capture checkerboard hybridization assay was used, as previously described (11), to detect levels (abundance) of 82 oral bacterial species or groups (see the Appendix, below). Briefly described, reverse-capture DNA probes (complementary oligonucleotide DNAs of known sequence) are used to target polynucleotides of unknown sequence (16S rRNA) bacterial genes in the biological sample solution. Probes were placed on a nylon membrane in separate horizontal lanes using a Mini Slot apparatus. 16S rRNA genes from plaque samples were PCR amplified using a specific labeled primer. Hybridizations were performed in vertical channels in a Miniblotter apparatus with labeled amplicons (target 16S rRNA genes) for up to 45 samples. A total of 1,350 hybridizations were performed simultaneously using a single membrane. Standard chemifluorescence detection was performed using the Storm system (Amersham, Piscataway, NJ). For each spot on the membranes, signal levels were extracted from their background by applying spot edge detection methodology. This method locates the average intensity around the spot's outline and then applies this as the background for the spot. The background was therefore calculated independently for each spot, and signal levels (normalized to mean counts) were calculated independently for each spot (ImageQuant software; Amersham, Piscataway, NJ). Low-quality spots were also filtered for quality control, and background noise was eliminated from the analysis. Two lanes in each run had the universal probes to serve as standards, and signal levels were converted to mean counts by comparison with standards on the membrane. Signal levels were then adjusted for abundance by comparing them to the universal control probes. This approach allowed for computing the abundance of the target species individually by adjusting the DNA concentration in each sample.
Computation of microbial abundance patterns. To determine if differences in bacterial levels were present between caries-free subjects and caries-active subjects, the Wilcoxon rank sum test for nonparametric distribution was first used for mean comparisons across all 82 bacterial species. However, classification approaches employed in biological samples should be capable of accurately representing not only differences across groups but also patterns that reflect true relationships among samples and disease profiles. For that reason, the Gene Expression Data Analysis software (9) was used for the classification of global patterns of microbial colonization and abundance levels in the two groups, and data were evaluated as follows: (i) the data were cube root transformed as a means of stabilizing the variance of the data across the full range of bacterial levels combined from different checkerboards, (ii) a simple univariate test was used to determine the overabundance and underabundance of bacterial species in the diseased and in the healthy groups, (iii) significant bacterial species in the two groups were identified using the Baïve Bayes classifier, and (iv) computation validation (random resampling validation) was performed using the significant species as features. In this model, we used a classification framework named PACE (permutation-achieved classification error) (6), which allowed us to compare the performance statistic of the classifier on true data samples and to perform validation procedures in the same data with randomly reassigned class labels. Our modeling approach split the data set into two groups of samples, a training and a test set, by using a 70/30 ratio. Learning (classifier construction) was conducted on the training set and evaluated on the test set. Features (differentially expressed bacterial species) were selected using the specific test and threshold, and class prediction was performed on the test data set (Baïve Bayes posterior probability). Sensitivity and specificity of the model were measured using the area under the receiver/operator curve (ROC). A statistically significant result increases the chance of a discriminative signal in the data. The objective of performing classification analyses was to build a predictive model to determine whether, and how, the differentially abundant species could be used to determine if differences were present in the biofilm of the caries-active group relative to the caries-free group.
Two approaches were used for analysis (Fig. 1). Model 1 compared site 1 (biofilm of surfaces of intact enamel) in the caries-free group with all sites combined in the caries-active group: site 1 (biofilm of surfaces of intact enamel), site 2 (white spot lesions), site 3 (dentinal lesions), and site 4 (deep dentinal lesions) (1). In this model, a ratio was created to adjust bacterial levels within species by the total number of sites collected for each caries-active subject. Model 2 compared site 1 (biofilm present on surfaces of intact enamel) from the caries-free group to only site 1 (biofilm present on surfaces of intact enamel) from the caries-active group.
|
|
|---|
Model 1. Model 1 compared site 1 (biofilm of surfaces of intact enamel) in the caries-free group with all sites combined in the caries-active group, i.e., site 1 (biofilm of surfaces of intact enamel), site 2 (white spot lesions), site 3 (dentinal lesions), and site 4 (deep dentinal lesions) (1). The achieved classification error was statistically significant under ROC (specificity, 0.86; sensitivity, 0.84) with a permutation test and PACE analysis at the 99% level. A total of 10 species were found to be significantly overabundant in the caries-active subjects, and 10 species were found to be significantly underabundant in this group compared to caries-free subjects (Fig. 2).
![]() View larger version (25K): [in a new window] |
FIG. 2. Overabundance and underabundance of bacterial species in caries-active subjects. Using model 1, the combined sites (intact enamel [1], white spot lesions [2], enamel lesions [3], and dentin lesions [4]) in caries-active subjects were compared to the biofilm of caries-free subjects. Species that are overabundant in the caries-active group relative to the caries-free group are shown in red; species that are underabundant in the caries-active group relative to the caries-free group are shown in green.
|
![]() View larger version (23K): [in a new window] |
FIG. 3. Overabundance and underabundance of bacterial species in caries-active subjects. Using model 2, the biofilm of caries-active subjects was compared to the biofilm of caries-free subjects. Species that are overabundant in the caries-active group relative to the caries-free group are shown in red; species that are underabundant in the caries-active group relative to the caries-free group are shown in green.
|
|
View this table: [in a new window] |
TABLE 1. Detection of bacterial species in subjectsa
|
|
|
|---|
The results of our study suggest that an abundance of the cariogenic microflora, i.e., Actinomyces sp., S. mutans, and Lactobacillus sp., in concert with an underabundance of health-associated species (S. parasanguinis, Abiotrophia defectiva, Gemella hemolysans, S. mitis/oralis, and S. sanguinis) represented a profile of the caries-active group (Fig. 2 and 3). By using this model, it was possible to identify microbial profiles that distinguished caries-active subjects from caries-free subjects.
In this study we employed two approaches of analysis. The first model looked at bacterial colonization present in the biofilm of caries-free subjects compared to the different biofilms of sites, with different lesion severity, in the caries-active subjects (Fig. 1). In the second model the microflora present in the biofilms of caries-free subjects was compared with only the microflora present in biofilms of surfaces of intact enamel in the caries-active subjects (Fig. 1). The same patterns of bacterial overabundance and underabundance were found to be significant in both models (Fig. 2 and 3). Notably, by employing model 2 one could conceivably use bacterial biofilms of healthy surfaces of caries-active subjects to profile subjects who may be at increased risk for developing caries. Thus, this model predicted that biofilm bacteria abundance on enamel intact surfaces may provide useful information for early diagnosis and detection of dental caries.
The results of this study confirmed those of previous reports (1, 8) by indicating that the microflora associated with caries is dominated by gram-positive bacteria, particularly the genera Actinomyces, Lactobacillus, and Streptococcus, and more predominantly by Streptococcus mutans (1). Earlier culture-based studies (13) implicated mutans streptococci in caries initiation, with the lactobacilli as important contributory bacteria. These results corroborate our findings in that S. mutans was found to be overabundant in about 90% of caries-active subjects and underabundant in the caries-free subjects (Fig. 2 and 3). Similarly, Lactobacillus spp. were found with high frequency in caries-active subjects but found only in one subject in the caries-free group (Table 1). Most studies (13) have found that although Actinomyces spp. are commonly detected in the human mouth, their role in carious lesions is variable and inconclusive. However, in our study Actinomyces spp. were found in abundance in both caries-active and caries-free subjects but with statistically significant differences between the groups (Table 1). One particular species (Actinomyces sp. strain B19SC) was found in very high levels and was significantly overabundant in caries-active subjects relative to caries-free subjects (both models) (Fig. 2 and 3).
Caries-active subjects in this population had a distinctive microflora represented by an underabundance of beneficial health-associated species (Fig. 2 and 3). Conversely, caries-free subjects exhibited an overabundance of known beneficial species (data not shown), i.e., S. parasanguinis, A. defectiva, G. hemolysans, S. mitis/oralis, Streptococcus cristatus, and S. sanguinis. Streptococcus sanguinis has long been associated with health (13). To understand what constitutes health, before disease ensues, requires ascertaining which health-associated species play an essential role, since beneficial species may protect the host against harmful pathogens, thereby preventing the onset and development of dental caries.
Lastly, results from our study suggest that research is needed to identify ways to translate these findings into clinically useful diagnostic models. We speculate that the most effective treatment strategies and preventive measures may well depend on elucidation of specific microbial profiles in health and disease. Biofilm bacteria can be very resistant to antibiotic treatment, and the mechanisms by which the biofilm-grown bacteria attain this resistance are still unknown. Moreover, biofilm bacteria infections are rarely resolved by the host's immune system, and they are resistant to the host's defense mechanisms (2). Knowing which bacteria are beneficial can help in the development of new treatment strategies for caries, whereby one manipulates these communities so as to prevent disease and promote health as an adjunct to standard prevention treatments.
|
|
|---|
|
View this table: [in a new window] |
TABLE A1. Bacterial species and groups used in checkerboard analysis
|
|
|
|---|
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»