Periodontitis affects up to 50% of individuals worldwide, and 8.5% are diagnosed with diabetes. The high-comorbidity rate of these diseases may suggest, at least in part, a shared etiology and pathophysiology. Changes in oral microbial communities have been documented in the context of severe periodontitis and diabetes, both independently and together. However, much less is known about the early oral microbial markers of these diseases. We used a subset of the ORIGINS project dataset, which collected detailed periodontal and cardiometabolic information from 787 healthy individuals, to identify early microbial markers of periodontitis and its association with markers of cardiometabolic health. Using state-of-the-art compositional data analysis tools, we identified the log-ratio of Treponema to Corynebacterium bacteria to be a novel Microbial Indicator of Periodontitis (MIP), and found that this MIP correlates with poor periodontal health and cardiometabolic markers early in disease pathogenesis in both subgingival plaque and saliva.
The human oral cavity hosts hundreds of unique microbial taxa. Within the oral cavity, there are multiple distinct niches that contain different compositions of microbial taxa. For example the supra- and subgingival tooth surface, tongue, cheek, and roof of mouth each harbor consistently distinct microbial communities1. In the context of severe periodontal disease, the composition of microbial taxa in the supra and subgingival plaque undergo dramatic changes2,3,4. Periodontitis-associated subgingival biofilms often contain a climax community dominated by Porphyromonas gingivalis, Treponema denticola, and Tannerella forsythia, referred to as the ‘red complex’5. The red complex has been studied in depth for its ability to negatively affect host physiology through virulence factors and expedite gingival deterioration in severe diseases. However, less is known about the early microbial markers of periodontitis and when compositional changes in plaque biofilms occur relative to disease onset.
Periodontal disease affects nearly half the global population6 and is one of the leading causes of tooth loss7. This disease has two primary stages, gingivitis and periodontitis. Gingivitis is the early stage of gum disease where the gingival soft tissue becomes swollen and may bleed upon provocation due to the presence of vasculitis and a local inflammatory infiltrate secondary to bacterial challenge8. With time, inflammation results in a deepening of the periodontal pocket, further expanding the ecological niche enabling greater plaque accumulation and host response to pathogens in the biofilm. Left unresolved, the early stages of periodontitis occur when the local inflammatory processes induce a progressive loss of tooth-supporting tissues and eventually tooth loss.
There is strong evidence for a link between oral health and cardiometabolic health9,10,11,12,13. While it has long been recognized that individuals with diabetes are at higher risk for periodontitis, recent evidence suggests that adverse oral microbial exposures that underlie periodontitis might also contribute to the etiology of cardiometabolic diseases, including diabetes. Few studies have investigated the relationship between oral microbial communities and markers of both periodontal disease and diabetes risk. Identifying microbial signatures that emerge early in disease could help elucidate potential shared microbial etiology of periodontitis and cardiometabolic disease.
Next-generation sequencing has enabled the high-throughput collection of detailed microbial information from thousands of samples. However, this data is inherently compositional, meaning that the results provide information only on microbial relative abundance rather than absolute abundance14. Therefore, care must be taken to avoid false positive or negative findings in these datasets, particularly because total microbial load can vary largely among subgingival plaque samples. In this analysis, we use a suite of compositional data analysis tools to identify differentially abundant bacteria in relation to measures of periodontal disease and early biomarkers of cardiometabolic disease.
We applied these tools to a subset of data collected in the Oral Infections, Glucose Intolerance and Insulin Resistance Study (ORIGINS)15. This dataset contains information from a large group of individuals free of clinical cardiometabolic disease, including a comprehensive periodontal examination, quantitative cardiometabolic markers, as well as 16 S rRNA gene amplicon sequencing from saliva and subgingival plaque from both healthy and diseased sites, processed separately, which allowed us to look for site-specific microbial markers that could be early indicators of disease.
Overview of the cohort
Wave 2 of the ORIGINS project recruited 800 participants, 787 of which had complete data collection for the current analyses. Each participant underwent extensive periodontal examination and metabolic measurements as previously described15. The subgingival plaque was collected from healthy (‘shallow’ periodontal pocket <3 mm) and diseased sites (‘deep’ periodontal pocket ≥ 4 mm) as applicable following a standardized protocol totaling 1107 plaque samples (Fig. 1). Saliva samples were collected in parallel and processed for a subset of 282 participants. Both saliva and subgingival plaque samples were processed for 16 S rRNA gene amplicon sequencing. The average full-mouth pocket depth (PD) was 1.90 mm, the average attachment loss was 0.59 mm, and 0.31% of sites exhibited bleeding on probing (BOP). 71%, 26%, and 2% of participants had no or mild periodontitis, moderate periodontitis or severe periodontitis, respectively, according to the CDC/AAP definition. As expected, the prevalence and relative abundance of canonical subgingival plaque pathogens from the red and orange complex5 were higher in deep versus shallow periodontal pockets (Table S1).
Periodontal pocket depth is strongly associated with microbial diversity
We used the robust Aitchison Principal Components Analysis (RPCA) method16, which accounts for the inherent sparsity and compositionality of next-generation sequencing experiments14, to assess beta-diversity across subgingival plaque samples. The first axis of separation showed distinct clusters of subgingival plaque collected from shallow versus deep periodontal pockets (Fig. 2A). Interestingly, the microbial composition of subgingival plaque samples from shallow periodontal pockets from different individuals were more similar to each other than the microbial composition of subgingival plaque samples from a shallow and deep periodontal pocket in the same person (Fig. 2B). This finding was only identified using RPCA and not other beta-diversity metrics (including weighted and unweighted UniFrac17, Jaccard, and Bray–Curtis), highlighting the ability of this tool to identify novel beta-diversity patterns that accord with expected patterns from prior work reporting consistent patterns of microbial dysbiosis in periodontal subgingival plaque samples16,18,19.
We performed an effect size redundancy analysis (RDA) to determine which factors explained the variation in microbial composition across samples20. Thirteen non-redundant factors were included in the RDA including, six periodontal metrics (shallow or deep periodontal pocket depth, percent of sites bleeding on probing, full-mouth classification of periodontal status using the CDC/AAP periodontitis definition, average whole-mouth pocket depth, average whole-mouth attachment loss, and percent of sites with attachment loss >3), three demographic factors (participant, sex, and age), three metabolic factors (fasting insulin, prediabetes status, and average systolic blood pressure) and one lifestyle factor (tobacco-smoking status). Of these thirteen factors, eight were found to have a significant effect size (p-value < 0.05).
The factor which explained the most variation was the disease status of the site adjacent to plaque sampling (i.e., a deep or shallow periodontal pocket (Fig. 2C)). Overall, periodontal metrics accounted for most of the explained variance (20%) followed by demographic factors (1.7%) and lastly metabolic factors (0.6%). Tobacco smoking did not have a significant effect on subgingival plaque microbial composition in this analysis (87% of the participants reported never smoking).
The multinomial regression tool Songbird21 was used to identify differentially abundant microbes in shallow versus deep periodontal pockets. Each amplicon sequence variant (ASV) was assigned a differential where higher scores reflect relative enrichment in shallow periodontal pockets, and low scores reflect relative enrichment in deep periodontal pockets (Table S2). The phylogenetic relationship among these ASVs and their associated Songbird differentials were visualized with EMPress18 (Fig. 2D). ASVs from the genus Treponema tended to have low (disease-associated) differentials, while ASVs from the genus Corynebacterium tended to have high (health-associated) differentials.
Identifying an early MIP in subgingival plaque
To further characterize these health- and disease- associated ASVs, we plotted the Songbird differentials from each ASV using the interactive tool Qurro22. As observed in the phylogenetic analysis, Treponema ASV differentials were more associated with deep periodontal pockets, whereas ASVs aligned to the genus Corynebacterium were more associated with shallow periodontal pockets (Fig. 3A). To generate a microbial indicator of periodontitis (MIP), we used Corynebacterium as a ‘reference frame’ and calculated the log-ratio of all Treponema counts to all Corynebacterium counts. This MIP was significantly higher in deep compared to shallow periodontal pockets (paired T-test < 0.0001) and revealed that the ratio of Treponema to Corynebacterium is roughly even in deep periodontal pockets, whereas the ratio in shallow periodontal pockets is heavily skewed towards Corynebacterium (Fig. 3B). Importantly, because these are relative abundance data, we cannot conclude whether this finding is due to an increase in Treponema or a decrease in Corynebacterium, but the ratio of these two organisms is a consistent biomarker of periodontal status.
To validate the robustness of the MIP, we tested the ability of the MIP to classify samples from shallow versus deep periodontal pockets. When using the entire dataset (1832 ASVs), samples were classified with an accuracy of 0.88 ± 0.04. When using just the subset of data used to generate the MIP (164 ASVs, or less than 10% of all ASVs), samples were classified with an accuracy of 0.83 ± 0.04% (Fig. 3C). To determine whether the MIP had significantly different accuracy compared to the whole dataset, stratified k-fold cross-validation with a 50:50 train and test set split was repeated 10-fold with random shuffling. On each split, a Random Forests classifier was trained and tested using the whole table versus the table filtered only for Treponema and Corynebacterium taxa, and the difference in the model performance was tested for significance on the contingency table between classifiers using a McNemar’s test23 (Table S3). Additionally, each fold classifier was evaluated by the ROC AUC between classifiers (Fig. 3C). In both evaluations, only a one-fold split showed a significant difference by the classifier, and the mean ROC AUC differed by only 5%. This demonstrates that the classification of the disease state is almost as accurate for Treponema and Corynebacterium measurements alone as all taxa.
Other microbial indicators of periodontitis have been proposed based on the analysis of chronic periodontitis24,25. We assessed these alternative microbial ratios from Chen et al.24, (Treponema denticola, Mogibacterium timidum, Fretibacterium spp., and Tannerella forsythia vs. Actinomyces naeslundii and Streptococcus sanguinis; representing 158 vs 238 ASVs, respectively), and Meuric et al.25, (Eubacterium, Campylobacter, Treponema, and Tannerella vs. Veillonella, Neisseria, Rothia, Corynebacterium, and Actinomyces; representing 27 vs. 9 ASVs, respectively) in this dataset. Both of these microbial ratios significantly discriminated healthy from diseased subgingival plaque samples. Random Forests classification with McNemar significance testing as described above showed no significant improvement of these alternative taxa ratios for discriminating between healthy versus diseased samples over the more focused Treponema:Corynebacterium MIP. We focused on the Treponema:Corynebacterium MIP as a proxy for early microbial dysbiosis to determine its association with physiological parameters of early periodontitis and cardiovascular disease.
Interestingly, even in subgingival plaque samples from shallow sites (n = 779), the MIP was significantly associated with the percent of sites bleeding on probing across the whole mouth (Pearson correlation = 0.243, p-value = 8.06e−12) (Fig. 3D). This indicates that even before clinically meaningful pocketing develops; there are microbial changes in the subgingival plaque related to gingival inflammation.
Evaluating the MIP in saliva
Collection of subgingival plaque is not trivial and requires clinically trained professionals. Saliva requires less time and participant burden to collect and can be done remotely. We performed 16 S rRNA gene amplicon sequencing on a subset of saliva samples (n = 282) collected in parallel to the subgingival plaque samples from the same individuals. Because the saliva and subgingival plaque samples were processed with different sequencing strategies at different institutions, we first assigned taxonomy to the Human Oral Microbiome Database (HOMD), a well-curated database of full-length 16 S rRNA gene amplicon sequences found in the human oral cavity26. We collapsed the subgingival plaque and saliva datasets to the species level and merged the tables together. Beta-diversity analysis of the merged table revealed that, as suspected from previous research, saliva and subgingival plaque had compositionally distinct microbial communities (Fig. 4A). The majority of microbial taxa in the merged table were found in both subgingival plaque and saliva, although each niche also contained distinct microbiota, with subgingival plaque being more diverse than saliva (Fig. 4B).
Effect size analysis using RDA of just the saliva table revealed different factors that drive microbial diversity in saliva compared to subgingival plaque (Fig. 4C). Eight factors were included in the RDA, including two demographic factors (participant and age), three metabolic factors (average systolic blood pressure, BMI, and prediabetes status), two periodontal factors (average whole-mouth attachment loss and average whole-mouth periodontal pocket depth), and one lifestyle factor (tobacco-smoking status). Overall, the percent explained was much lower in saliva compared to subgingival plaque (5.8% versus 22.3%, respectively). The only significant factors in the saliva RDA were smoking status and participant (Fig. 4C). Despite clear differences in microbial community composition between saliva and subgingival plaque, the MIP was significantly correlated (Pearson R = 0.387, p-value = 3.97E−11) (Fig. 4D).
MIP and periodontal measures
We assessed the correlation of the MIP with various whole-mouth periodontal measures. Subgingival plaque MIP was positively, significantly correlated with the percent of sites bleeding on probing, mean full-mouth pocket depth, and mean full-mouth attachment loss (Table 1). Subgingival plaque MIP was also very highly correlated with Faith’s phylogenetic diversity, a measurement of alpha diversity that accounts for phylogenetic relatedness27. All of these correlations held true when looking at only samples from healthy sites, again suggesting that microbial community shifts precede detectable disease (Table S4).
Saliva MIP was significantly correlated with average full-mouth pocket depth, but not percent of sites bleeding on probing or attachment loss (Table 2). This held true when looking only at participants with moderate or severe periodontitis but not healthy participants, suggesting that microbial changes in the subgingival pocket precede microbial changes in saliva (Table S5). Saliva MIP was also strongly correlated with Faith’s phylogenetic diversity in all patients.
MIP and metabolic metrics
The overarching goal of the ORIGINS project is to identify associations between oral microbes, oral health, and cardiometabolic health. To this end, we evaluated the correlation of the MIP with various cardiometabolic health metrics. Subgingival plaque MIP was positively, significantly correlated with average systolic and diastolic blood pressure, fasting insulin levels and HOMA-IR, after adjusting for age, sex, race, BMI, and smoking status (Table 3). Together these results suggest that the microbial changes underlying periodontal health are also associated with cardiometabolic health.
Saliva MIP was not significantly correlated with cardiometabolic metrics after adjusting for age, sex, race, BMI, and smoking status (Table 4).
Quartiles of MIP and metabolic metrics and periodontal measures
We evaluated the correlation of quartiles of subgingival plaque MIP with various periodontal measures and cardiometabolic health metrics (Table 5). As quartiles of subgingival plaque MIP increase, mean full-mouth attachment loss, mean full-mouth pocket depth, percent of sites bleeding on probing, and Faith’s phylogenetic diversity all increased significantly even after adjusting for age, sex, race, BMI, and smoking status. Similarly, average systolic and diastolic blood pressure, fasting insulin levels, and HOMA-IR are positively, significantly correlated with quartiles of MIP. In addition, we evaluated the correlation of quartiles of saliva MIP with various periodontal measures and cardiometabolic health metrics and (Table 6). Quartiles of saliva MIP were only positively, significantly correlated with mean full-mouth pocket depth and Faith’s phylogenetic diversity after adjusting for age, sex, race, BMI, and smoking status.
In a cohort of 787 healthy individuals, we were able to identify early microbial markers of periodontal disease. Microbial composition in subgingival plaque was most strongly explained by periodontal metrics such as subgingival pocket depth and percent of sites bleeding on probing. Both RPCA beta-diversity analysis and redundancy analysis (RDA) showed that the microbial composition in periodontal plaque was more similar between different individuals with the same periodontal phenotype (i.e., shallow or deep pockets), compared to plaque from a shallow versus deep periodontal site within the same person. This finding replicates prior publications showing similar patterns. Since previous studies have shown increased microbial burden in subgingival pockets with periodontitis28, it is likely that microbial load varied greatly across the samples in this dataset, and therefore, it is crucial to use scale-invariant analyses. This point is underscored by the fact that this result was not observed using metrics that do not account for compositionality (e.g., UniFrac, Bray–Curtis), which can be greatly affected by microbial load21.
RDA further revealed that saliva microbial communities were influenced by different factors compared to subgingival plaque. For instance, while tobacco smoking did not have a significant effect size in subgingival plaque microbial composition, it had the biggest effect size in saliva microbial composition. This is in line with previous reports showing that microbial composition in oral washes was affected by smoking status29, while subgingival plaque is not greatly affected by smoking status30. Together, these results demonstrate that saliva and subgingival plaque microbial communities are driven by different environmental factors.
We used the factor with the highest effect size on microbial diversity in subgingival plaque, whether the sample came from a deep or shallow pocket, to identify a microbial indicator of early periodontal disease. Using reference frames, we calculated the log-ratio of Treponema:Corynebacterium and found that it significantly differentiated healthy from diseased periodontal pocket sites. This log-ratio was used as a MIP. Previous studies have proposed alternative MIPs24,25, which were also significantly different between health and disease in this dataset, but were not statistically more predictive than the Treponema:Corynebacterium ratio for discriminating disease type. We must acknowledge that identifying microbial biomarkers in next-generation sequencing datasets always carries the risk for false positives, but in this study, we focused on the Treponema and Corynebacterium ASVs because the ability of all the species and strains of these genera to significantly discriminate sites of early disease before clinically meaningful deep pockets formed suggests that these phylogenetic branches have been evolutionarily conserved in subgingival plaque biofilm formation.
Subgingival plaque MIP was significantly correlated with poor periodontal health across a wide range of metrics when only looking at healthy plaque samples, suggesting that microbial communities change before the disease is clinically detectable. Red complex organisms canonically associated with periodontitis in the literature were also positively correlated with periodontal disease status, but they were not as widely prevalent across samples which complicates scale-invariant analyses. This is likely due to the fact that periodontitis prevalence was low and the extent and severity of disease was relatively mild and in the early stages compared to many other studies of periodontitis.
Remarkably, phylogenetically informed alpha diversity was also strongly correlated with MIP in both saliva and subgingival plaque across all periodontal status categories. Both Treponema and Corynebacterium species have been identified as microbial scaffolds in plaque biofilms. In the context of healthy periodontal plaque, reproducible biofilms with a specific taxonomic organization, referred to as ‘hedgehog’ biofilms, are widely prevalent31. In the context of severe periodontitis, Treponema taxa have been found in the deepest sections of the periodontal pocket, as they are especially sensitive to oxygen compared to Corynebacterium, and form close associations with diverse rod-like bacteria32,33,34. In light of our finding that the ratio of Treponema to Corynebacterium increases in periodontal disease, this suggests that the biofilm structure shifts from being scaffolded primarily by Corynebacterium to Treponema, where Treponema biofilms are more phylogenetically diverse than Corynebacterium biofilms.
One outstanding question that emerges from this analysis is whether microbial dysbiosis is the cause or symptom of periodontitis; does the subgingival pocket first get deeper, providing a microbial niche which Treponema species are particularly adept at colonizing? Or does Treponema colonization and protease release cause subgingival pocket deepening? The data here are suggestive that Treponema first colonizes the periodontal pocket and drives disease, because we see the MIP ratio correlating with pocket depth even before the site would be labeled as diseased (i.e., <3 mm depth). Further, Treponema can move through viscous environments, including between and through tissue35. However, further longitudinal studies are required to validate this hypothesis.
Despite the fact that saliva has a compositionally distinct microbiome compared to subgingival plaque, is driven by different metadata variables, and was sequenced independently with different parameters, we found that the MIP was significantly correlated between plaque and saliva. Saliva MIP was also correlated with poor periodontal health, although only in participants with moderate to severe periodontitis. Previous work has found that saliva flow shapes microbial organization in plaque36, providing a potential explanation for the correlation in microbial dysbiosis in these distinct microbial communities. One limitation of this analysis was that only a subset of samples had paired saliva sequencing, and increasing the number of saliva samples analyzed could further clarify the potential for saliva microbial composition as a readout for periodontal status, as has been previously suggested by smaller-scale studies37,38,39.
Periodontal disease and cardiometabolic health have been found to co-associate across diverse populations40,41. This association is potentially bidirectional. On the one hand, poor cardiovascular health has been suggested to increase the risk of periodontitis42, where factors like dysglycemia, receptor for advanced glycation end-product activity and immunological response could contribute to disease. On the other hand, microbial dysbiosis common in periodontitis has been suggested as a risk factor for cardiometabolic disease, since it can evoke persistent inflammation41,43. To further explore these associations, we evaluated the correlation of our MIP with various cardiometabolic measurements. We show that specific taxa associated with periodontal disease (the MIP) in subgingival plaque are also significantly correlated with blood pressure, fasting insulin, and HOMA-IR. There was a correlation between subgingival plaque MIP and cardiovascular health markers even in participants with no periodontitis. This suggests that the early origins of these diseases are intricately linked. One potential mechanism could be changed in the enterosalivary nitrate metabolism pathway, which affects systemically available nitric oxide and directly influences cardiometabolic outcomes44. For example, Treponema colonization and subsequently increased alpha diversity could be leading to depletion of oral nitrate reducers and nitrite depletion. However, any potential cause-and-effect relationship remains to be determined.
Through this large-scale analysis of wave 2 of the ORIGINS cohort15, we identified a simple microbial signature based on two common oral taxa associated with periodontal status that also correlates with biomarkers of cardiometabolic disease. Importantly, these microbial community composition transitions appear to occur early in disease before severe periodontitis is evident or clinical cardiovascular disease develops. The results from this analysis suggest that these microbial changes occur first in plaque, and as disease progresses can be identified in saliva. Future longitudinal sampling will allow for more definitive determinations of how the Treponema:Corynebacterium ratio is potentially involved in the pathogenesis of periodontal and cardiometabolic diseases.
ORIGINS is an occupation-based cohort study among members of the Service Employees International Union 1199 designed to investigate the relationship between oral microbial community composition and glucose metabolism. Periodontal examination, subgingival plaque, and saliva collection were performed as previously described15. In summary, 1188 subgingival plaque samples (4 samples from 297 participants) were collected from the most posterior tooth per quadrant (excluding third molars) via sterile curettes after removal of the supragingival plaque. Unstimulated saliva was collected from each participant in parallel.
The Institutional Review Board at Columbia University and the University of Minnesota approved the study protocol. All participants provided written informed consent.
DNA extraction and 16 S rRNA gene sequencing
DNA was extracted from subgingival plaque and saliva samples at The Forsyth Institute. 16 S rRNA gene amplicon sequencing was performed on subgingival plaque samples by The Forsyth Institute using primers targeting variable regions 3 and 4; Forward- CCTACGGGAGGCAGCAG (341 f) and Reverse- GGACTACHVGGGTWTCTAAT (806r). Sequencing was performed on a MiSeq using a Paired End 250 cycle kit.
16 S rRNA gene amplicon sequencing libraries on DNA extracted from saliva was performed at UC San Diego using the Earth Microbiome Project protocol45,46 with primers targeting the v4 region; Forward- GTGYCAGCMGCCGCGGTAA (515 f) and Reverse- GGACTACNVGGGTWTCTAAT (806r). Sequencing was performed on a MiSeq using a Paired End 150 cycle kit.
Raw reads were analyzed with QIIME2. Demultiplexed sequences were quality-filtered with default parameters in qiime quality-filter q-score; namely, reads were trimmed after the first appearance of three basecalls with a PHRED score of four or less, and the entire read was removed if the read was truncated to <75% of the input sequence. Quality filtered forward-read sequences were denoised using Deblur47 with the default parameters. Samples with less than 1000 quality-filtered reads were removed from downstream analysis. In order to remove reads aligned to chloroplast or mitochondrial genes, sequences were aligned using a classifier pretrained on the GreenGenes v13_8 database with 99% sequence homology using sklearn48. Sequences aligned to mitochondria or chloroplast were removed using filter-table --p-exclude (0.005% of the entire dataset). A phylogenetic tree was created using fragment insertion via SEPP49. Taxonomy was assigned using sklearn48 against the HOMD database version 15.126. All features not present in at least 1% of samples were excluded from downstream analysis.
The final quality-filtered subgingival plaque table contained 43,709,128 reads (97.6% of the raw dataset) across 1107 samples (99.2% of total collected samples) with a total of 1832 amplicon sequence variants (ASV) (26.9% of all identified ASVs), of which 99.5% were assigned to at least the phylum level. The final quality-filtered saliva table contained 4,892,251 (99.6% of the raw dataset) reads across 282 samples (100% of sequenced samples) with a total of 859 ASVs (73.0% of all identified ASVs), of which 99.5% were assigned to at least the phylum level.
Differential abundance testing
To determine which taxa are associated with which phenotypes in our dataset, we used the concept of Reference Frames21. This tool accounts for the compositional nature of next-generation sequencing experiments14. In brief, comparing relative abundances among sample groups can be misleading when the total microbial load is unknown, as is the case in this dataset. To avoid these pitfalls, we used the machine learning tool Songbird (https://github.com/biocore/songbird) to perform multinomial regression and then ranked each ASV by its coefficient in the regression model to determine each taxon’s relative differential across a given phenotype. Periodontal pocket depth was used as the formula in the model. The number of random test samples held back for validation in the model was 111 (10%). We used a batch size of ten with 500 epochs (number of passes through the entire dataset to train the model), a learning rate of 0.001 and a differential prior of 10. The resulting ranks (differentials.qza) were visualized with Qurro22 and allowed us to prioritize which taxa were most associated with a given phenotype.
To identify taxa associated with shallow versus deep periodontal pockets, we browsed the highest and lowest-ranked microbes in this category using Qurro. ASVs assigned to the genus Corynebacterium were mostly associated with shallow pockets, whereas ASVs assigned to the genus Treponema were mostly associated with deep pockets. To generate a microbial indicator of periodontitis (MIP), we used Corynebacterium as a ‘reference frame’ and calculated the log ratio of all Treponema counts to all Corynebacterium counts.
A Random Forests (Breiman 2001) (RF) model was trained to predict disease status based on shallow (pocket depth <3 mm) versus deep (pocket depth > 4 mm) periodontal pockets. The RF model was trained using a Stratified K-Folds cross-validation (CV) with 10-Fold CV splits. On each CV split, an RF model with 500 estimators was trained, and RF probability predictions were compared to the test set using the Receiver Operating Characteristic (ROC). The mean and standard deviation from the mean were calculated for the Area Under the Curve (AUC) across the 10-fold CV. This classification was performed on the whole ASV level data table and compared to the table filtered for only members of Treponema and Corynebacterium concatenated with the log-ratio of Treponema to Corynebacterium. All classification was performed through Scikit-learn (v. 0.22.2)48. To determine statistical significance, stratified k-fold cross-validation with a 50:50 train and test set split was repeated 10-fold with random shuffling. On each split, a Random Forests classifier was trained and tested using the whole table, versus the table filtered only for Treponema and Corynebacterium taxa, and the difference in the model performance was tested for significance on the contingency table between classifiers using a McNemar’s test23.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
All sequencing data and sample metadata are available through Qiita50 under study ID 11808 for saliva data and study ID 14375 for subgingival plaque samples. Raw sequence data is also available through EBI accession PRJEB50306 for saliva samples and accession PRJEB50261 for subgingival plaque samples.
Xu, X. et al. Oral cavity contains distinct niches with dynamic microbial communities. Environ. Microbiol. 17, 699–710 (2015).
Abiko, Y., Sato, T., Mayanagi, G. & Takahashi, N. Profiling of subgingival plaque biofilm microflora from periodontally healthy subjects and from subjects with periodontitis using quantitative real-time PCR. J. Periodontal Res. 45, 389–395 (2010).
Heller, D., Silva-Boghossian, C. M., do Souto, R. M. & Colombo, A. P. V. Subgingival microbial profiles of generalized aggressive and chronic periodontal diseases. Arch. Oral. Biol. 57, 973–980 (2012).
Kageyama, S. et al. Relative abundance of total subgingival plaque-specific bacteria in salivary microbiota reflects the overall periodontal condition in patients with periodontitis. PloS One 12, e0174782 (2017).
Socransky, S. S., Haffajee, A. D., Cugini, M. A., Smith, C. & Kent, R. L. Microbial complexes in subgingival plaque. J. Clin. Periodontol. 25, 134–144 (1998).
Nazir, M. A. Prevalence of periodontal disease, its association with systemic diseases and prevention. Int J. Health Sci. 11, 72–80 (2017).
Meyer, M. S., Joshipura, K., Giovannucci, E. & Michaud, D. S. A review of the relationship between tooth loss, periodontal disease, and cancer. Cancer Causes Control19, 895–907 (2008).
Könönen, E., Gursoy, M. & Gursoy, U. K. Periodontitis: a multifaceted disease of tooth-supporting tissues. J. Clin. Med. 8, E1135 (2019).
Preshaw, P. M. et al. Periodontitis and diabetes: a two-way relationship. Diabetologia 55, 21–31 (2012).
Kebschull, M., Demmer, R. T. & Papapanou, P. N. “Gum bug, leave my heart alone!”—epidemiologic and mechanistic evidence linking periodontal infections and atherosclerosis. J. Dent. Res 89, 879 (2010).
Lalla, E. & Papapanou, P. N. Diabetes mellitus and periodontitis: a tale of two common interrelated diseases. Nat. Rev. Endocrinol. 7, 738–748 (2011).
Lockhart, P. B. et al. Periodontal disease and atherosclerotic vascular disease: does the evidence support an independent association?: a scientific statement from the American Heart Association. Circulation 125, 2520–2544 (2012).
Emerging Risk Factors Collaboration. Diabetes mellitus, fasting blood glucose concentration, and risk of vascular disease: a collaborative meta-analysis of 102 prospective studies. Lancet Lond. Engl. 375, 2215–2222 (2010).
Gloor, G. B., Macklaim, J. M., Pawlowsky-Glahn, V. & Egozcue, J. J. Microbiome datasets are compositional: and this is not optional. Front. Microbiol. 8, 2224 (2017).
Demmer, R. T. et al. Periodontal bacteria and prediabetes prevalence in ORIGINS: the oral infections, glucose intolerance, and insulin resistance study. J. Dent. Res 94, 201S–211SS (2015).
Martino C. et al. A novel sparse compositional technique reveals microbial perturbations. mSystems https://doi.org/10.1128/mSystems.00016-19 (2019).
Lozupone, C. & Knight, R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ. Microbiol. 71, 8228–8235 (2005).
Chen, C. et al. Oral microbiota of periodontal health and disease and their changes after nonsurgical periodontal therapy. ISME J. 12, 1210–1224 (2018).
Shi, M. et al. The subgingival microbiome of periodontal pockets with different probing depths in chronic and aggressive periodontitis: a pilot study. Front Cell Infect. Microbiol. 8, 124 (2018).
Falony, G. et al. Population-level analysis of gut microbiome variation. Science 352, 560–564 (2016).
Morton, J. T. et al. Establishing microbial composition measurement standards with reference frames. Nat. Commun. 10, 2719 (2019).
Fedarko, M. W. et al. Visualizing’omic feature rankings and log-ratios using Qurro. NAR Genom. Bioinform. 2, lqaa023 (2020).
McNemar, Q. Note on the sampling error of the difference between correlated proportions or percentages. Psychometrika 12, 153–157 (1947).
Chen, T., Marsh, P. D. & Al-Hebshi, N. N. SMDI: an index for measuring subgingival microbial dysbiosis. J. Dent. Res. 101, 331–338 (2021).
Meuric, V. et al. Signature of microbial dysbiosis in periodontitis. Appl. Environ. Microbiol. 83, e00462–17 (2017).
Escapa, I. F. et al. New insights into human nostril microbiome from the expanded human oral microbiome database (eHOMD): a resource for the microbiome of the human aerodigestive tract. mSystems 3, e00187–18 (2018).
Faith, D. P. & Baker, A. M. Phylogenetic diversity (PD) and biodiversity conservation: some bioinformatics challenges. Evol. Bioinform. 2, 117693430600200 (2006).
Nastych, O. et al. Comparison of bacterial load parameters in subgingival plaque during peri-implantitis and periodontitis using the RT-PCR method. Acta Stomatol Croat. 54, 32–43 (2020).
Wu, J. et al. Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 10, 2435–2446 (2016).
Lanza, E. et al. Complementary clinical effects of red complex bacteria on generalized periodontitis in a caucasian population. Oral. Dis. 22, 430–437 (2016).
Mark Welch, J. L., Rossetti, B. J., Rieken, C. W., Dewhirst, F. E. & Borisy, G. G. Biogeography of a human oral microbiome at the micron scale. Proc. Natl Acad. Sci. USA 113, E791–E800 (2016).
Wecke, J. et al. A novel technique for monitoring the development of bacterial biofilms in human periodontal pockets. FEMS Microbiol Lett. 191, 95–101 (2000).
Chan, E. C. S., Vries, J., de, Harvey, R. F. & Tam, Y.-C. Use of fluorescence microscopy for monitoring periodontal disease state. Can. J. Microbiol. 27, 675–678 (2011).
Armitage, G. C., Dickinson, W. R., Jenderseck, R. S., Levine, S. M. & Chambers, D. W. Relationship between the percentage of subgingival spirochetes and the severity of periodontal disease. J. Periodontol. 53, 550–556 (1982).
Radolf, J. D. et al. Treponema pallidum, the syphilis spirochete: making a living as a stealth pathogen. Nat. Rev. Microbiol 14, 744–759 (2016).
Proctor, D. M. et al. A spatial gradient of bacterial diversity in the human oral cavity shaped by salivary flow. Nat. Commun. 9, 681 (2018).
Sun, X. et al. Alteration of salivary microbiome in periodontitis with or without type-2 diabetes mellitus and metformin treatment. Sci. Rep. 10, 15363 (2020).
Sabharwal, A. et al. The salivary microbiome of diabetic and non-diabetic adults with periodontal disease. J. Periodontol. 90, 26–34 (2019).
Lin, M. et al. Saliva microbiome changes in patients with periodontitis with and without chronic obstructive pulmonary disease. Front Cell Infect. Microbiol. 10, 124 (2020).
Zhao, M.-J. et al. Periodontal disease is associated with increased risk of hypertension: a cross-sectional study. Front. Physiol. 10, 440 (2019).
Muñoz Aguilera, E. et al. Periodontitis is associated with hypertension: a systematic review and meta-analysis. Cardiovasc. Res. 116, 28–39 (2020).
Koo, H. S. & Hong, S. M. Prevalence and risk factors for periodontitis among patients with metabolic syndrome. Metab. Syndr. Relat. Disord. 16, 375–381 (2018).
Czesnikiewicz-Guzik, M. et al. Causal association between periodontitis and hypertension: evidence from Mendelian randomization and a randomized controlled trial of non-surgical periodontal therapy. Eur. Heart J. 40, 3459–3470 (2019).
Goh, C. E. et al. Association between nitrate-reducing oral bacteria and cardiometabolic outcomes: results from ORIGINS. J. Am. Heart Assoc. 8, e013324 (2019).
Marotz, C. et al. DNA extraction for streamlined metagenomics of diverse environmental samples. BioTechniques 62, 290–293 (2017).
Thompson, L. R. et al. A communal catalogue reveals Earth’s multiscale microbial diversity. Nature 551, 457–463 (2017).
Amir, A. et al. Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems 2, e00191–e00216 (2017).
Pedregosa F. et al. Scikit-learn: Machine learning in Python. arXiv https://doi.org/10.48550/arXiv.1201.0490 (2018).
Janssen, S. et al. Phylogenetic placement of exact amplicon sequences improves associations with clinical information. mSystems 3, e00021–18 (2018).
Gonzalez, A. et al. Qiita: rapid, web-enabled microbiome meta-analysis. Nat. Methods 15, 796–798 (2018).
This research was supported by NIH grants R00 DE018739, R21 DE022422, and R01 DK 102932 (to R.T.D). R.T.D also received funding from a Calderone Research Award, Mailman School of Public Health, and a Pilot & Feasibility Award from the Diabetes and Endocrinology Research Center, College of Physicians and Surgeons (DK-63608). This publication was also supported by the National Center for Advancing Translational Sciences, National Institutes of Health, through Grant Number UL1TR001873. Clarisse Marotz was funded by NIDCR NRSA F31 Fellowship 1F31DE028478-01. Rebecca Molinsky was supported by institutional training grant T32HL007779 from the National Institutes of Health. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Marotz, C., Molinsky, R., Martino, C. et al. Early microbial markers of periodontal and cardiometabolic diseases in ORIGINS. npj Biofilms Microbiomes 8, 30 (2022). https://doi.org/10.1038/s41522-022-00289-w