Skip to main content

Anna Karenina and the subgingival microbiome associated with periodontitis

Abstract

Background

Although localized aggressive periodontitis (LAP), generalized aggressive periodontitis (GAP), and chronic periodontitis (CP) are microbially driven diseases, our inability to separate disease-specific associations from those common to all three forms of periodontitis has hampered biomarker discovery. Therefore, we aimed to map the genomic content of, and the biological pathways encoded by, the microbiomes associated with these clinical phenotypes. We also estimated the extent to which these biomes are governed by the Anna Karenina principle (AKP), which states that eubiotic communities are similar between individuals while disease-associated communities are highly individualized.

Methods

We collected subgingival plaque from 25 periodontally healthy individuals and diseased sites of 59 subjects with stage 3 periodontitis and used shotgun metagenomics to characterize the aggregate of bacterial genes.

Results

Beta-dispersion metrics demonstrated that AKP was most evident in CP, followed by GAP and LAP. We discovered broad dysbiotic signatures spanning the three phenotypes, with over-representation of pathways that facilitate life in an oxygen-poor, protein- and heme-rich, pro-oxidant environment and enhance capacity for attachment and biofilm formation. Phenotype-specific indicators were more readily evident in LAP microbiome than GAP or CP. Genes that enable acetate-scavenging lifestyle, utilization of alternative nutritional sources, oxidative and nitrosative stress responses, and siderophore production were unique to LAP. An attenuation of virulence-related functionalities and stress response from LAP to GAP to CP was apparent. We also discovered that clinical phenotypes of disease resolved variance in the microbiome with greater clarity than the newly established grading system. Importantly, we observed that one third of the metagenome of LAP is unique to this phenotype while GAP shares significant functional and taxonomic features with both LAP and CP, suggesting either attenuation of an aggressive disease or an early-onset chronic disease.

Conclusion

Within the limitations of a small sample size and a cross-sectional study design, the distinctive features of the microbiomes associated with LAP and CP strongly persuade us that these are discrete disease entities, while calling into question whether GAP is a separate disease, or an artifact induced by cross-sectional study designs. Further studies on phenotype-specific microbial genes are warranted to explicate their role in disease etiology.

Video Abstract

Introduction

Periodontitis, an infection-mediated disease that destroys tooth-supporting structures, is the sixth most prevalent disease in the world, affecting over 700 million adults worldwide [1]. The consequences of untreated disease are tooth loss, poor nutritional status, loss of speech, and masticatory function. With the annual cost of periodontal treatment exceeding 15 billion dollars in the USA alone, this disease poses a significant health burden that is comparable to outpatient treatment of cardiac and metabolic diseases [2,3,4]. Additionally, emerging evidence implicates periodontitis in the pathogenetic pathways of several potentially life-threatening diseases including coronary heart disease, pre-term births, diabetes, and cerebrovascular accidents [5], and therefore, the consequences of untreated periodontitis extend beyond the oral cavity.

Three clinical phenotypes have been most commonly described based on the rate of progression and the pattern of disease within the dentition: chronic periodontitis (CP), localized aggressive periodontitis (LAP), and generalized aggressive periodontitis (GAP) [6]. Although it is established that all three phenotypes have a microbially driven etiology, our current understanding of the microbiota associated with these phenotypes is not sufficient to explain the clinical differences. For example, the prevailing paradigm that similar bacteria are found in CP and GAP [7] does not explain the rapid disease progression seen in GAP, nor does the presence of specific bacteria in GAP and LAP [8,9,10] explain the involvement of specific teeth in LAP.

Our knowledge of the microbiome associated with chronic and aggressive periodontitis is largely based on phylogenetic characterizations of subgingival communities or investigations of specific bacteria (notably, Porphyromonas gingivalis, Treponema denticola, and Tannerella forsythia, (popularly known as the “red complex” bacteria) and Aggregatibacter actinomycetemcomitans) within these ecosystems [10,11,12]. However, these individual species are part of complex communities and their role in disease causation or perpetuation can be fully understood only when studied in an ecological context. While phylogenetic approaches such as 16S sequencing place these organisms in an ecological framework [13,14,15], their functional roles can only be inferred. Moreover, most studies are comparisons of only one disease phenotype to healthy controls [7, 16], limiting our ability to separate phenotype-specific associations from those that are common to all clinical presentations of periodontitis. It is not surprising that we have not been able to use microbiomics to develop a robust understanding of how the microbiome contributes to a specific disease phenotype.

High-throughput whole genome sequencing has provided us with an unprecedented view of the genetic composition and functional behaviors of complex biomes that can never be fully characterized by cultivation alone [17]. For example, metagenomic studies of the gut microbiome have revealed that eubiotic communities are similar between individuals but respond in a stochastic or random manner to stresses, resulting in dysbiotic communities that vary from person to person. This is called the Anna Karenina principle (AKP) after the opening line from Tolstoy’s book: “All happy families look alike; each unhappy family is unhappy in its own way” [18]. In the present study, we aimed to investigate what the microbiome is capable of doing and how this functional capacity relates to periodontal health status, in order to develop testable hypotheses about the role of the subgingival microbiome in maintaining health and causing disease. Here, we present the first functional catalog of the subgingival microbiome in the three most common phenotypes of periodontitis. Using a metagenomic approach to characterize the aggregate of bacterial genes in the subgingival microenvironment, we also highlight how each differs from periodontal health. As a tertiary aim, we examined the 1999 (based on disease phenotype) and the 2017 (based on disease extent, severity, and risk) classifications of periodontitis [6, 19] in the context of the subgingival microbiome.

Methods

Study population

This study was approved by the institutional review boards of The Ohio State University and Louisiana State University (OSU IRB 2014H0020, LSUHS-NO 8796). Thirty-four nonsmoking, normoglycemic individuals with stage 3 periodontitis as defined by the 2017 classification [19] were recruited and informed consent or assent with parental approval was obtained as appropriate. Additionally, sequences from 25 stage 3 periodontitis and 25 periodontally healthy controls from a previous study [20] were reanalyzed. Periodontal health was defined as clinical attachment loss (CAL) ≤ 1 mm, probing pocket depths (PD) ≤ 3 mm, and mean gingival index < 1. Periodontitis was classified based on both phenotype [6] and disease stage and grade [19] for comparison. Stage 3 grade A (S3gA) was defined as bone loss to age ratio (BL/age) < 0.25, mean plaque index (PI) > 1.5, and mean gingival index (GI) > 1.5; stage 3 grade B (S3gB) as 0.25 < BL/age < 1, PI > 1.5, and GI > 1.5; and stage 3 grade C (S3gC) BL/age > 1, PI ≤ 1, and GI ≤ 1. Chronic periodontitis (CP) was defined as interproximal attachment loss affecting at least 30% of the sites in individuals between 40 and 80 years of age. CP was diagnosed based on at least two periodontal assessments over at least 2 years. Localized aggressive periodontitis (LAP) was defined as interproximal attachment loss affecting first molars and incisors (and no more than 2 other teeth) in systemically healthy individuals below 40 years of age. An additional requirement was a contributory family history (of early tooth loss). Generalized aggressive periodontitis (GAP) was defined as rapid progressing interproximal attachment loss in systemically healthy individuals below 40 years of age that did not follow the patterns of distribution in the dentition described for LAP. Exclusion criteria included age below 8 years, current pregnancy, requirement for antibiotic prophylaxis prior to dental therapy, HIV infection, long-term (greater than 3 months) use of medications known to cause gingival changes (e.g., immunosuppressants, phenytoin, calcium channel blockers, aspirin, NSAIDS, bisphosphonates, or steroids), antibiotic therapy within 3 months of sample collection, and history of previous subgingival periodontal therapy.

Sample size was estimated based on the probability of least an 80% chance of detecting individual genes that differed in abundance by at least 0.01% between any two groups.

Clinical procedure

Patients were recruited after an initial screening and diagnosis visit when periodontal health related metrics (clinical attachment loss (CAL), probe depths (PD), bleeding on probing (BOP), gingival index (GI), and plaque index (PI)) were recorded. Since a common clinical presentation in all three disease phenotypes is that only some sites are affected by disease, 3 non-contiguous sites with clinical attachment loss (CAL) > 5 mm and probing depth (PD) ≥ 6 mm on 3 different teeth were selected in subjects with periodontitis. Samples were collected 1 week following the initial screening visit. This strategy was adopted to minimize the effects of bacterial translocation due to periodontal probing. Samples collected by inserting sterile endodontic paper points (Caulk-Dentsply, Milford, DE, USA) to the depth of the periodontal pocket for 30 s. From periodontally healthy subjects, samples were similarly collected and pooled from 15 non-contiguous interproximal sites. Paper points were immediately placed in 100 μL of RNALater, temporarily stored in ice, and were transferred to − 20 °C until analysis.

DNA isolation and sequencing

Bacterial DNA was isolated from paper points, using Qiagen DNA MiniAmp kit (Qiagen, Valencia, CA, USA) and quantified using Qubit fluorometer. Fifty nanograms of DNA was used to generate libraries with an Illumina TruSeq kit according to the manufacturer’s instructions. Quantified and pooled libraries were clustered on the Illumina HiSeq 4000 system (Illumina, Inc., San Diego, CA, USA) and 150 bp paired-end sequencing performed. Sequences for all 59 samples are deposited in the Sequence Read Archives under the project ID PRJNA552294 and PRJNA508385.

Statistical analysis

Trimmed and filtered sequences were uploaded to the MG-RAST metagenomics analysis pipeline (version 3.3.6) [21, 22] (Argonne National Laboratory) for quality processing and basic functional analysis. The MG-RAST API [23] and the custom Python library we developed to access it and analyze/visualize results were used throughout the analysis process to download relevant data and pipeline results (available for download at http://github.com/smdabdoub/PyMGRAST).

The phylogenetic profile of each group was determined using Kraken v1.1 [24] with a database constructed from a list of complete genomes from the Human Oral Microbiome Database [HOMD], as of September 19, 2017 (GenBank IDs available in Supplement). Although Kraken provides data at all levels of taxonomy, analysis was performed at the genus and species levels, showing high levels of specificity, accuracy, and coverage. Only taxa present in at least 20% of subjects in each group and constituting ≥ 0.000005% of the relative abundance were retained. Annotations with Kraken were done using the Ohio Supercomputer facilities (www.osc.edu). Knowledge-based assignments were made for gram staining and oxygen tolerance characteristics of the OTUs (Python script available on https://github.com/akshayparopkari/kadambari). As-yet-uncultivated phylotypes were annotated based on the characteristics of their cultivated phylogenetic neighbors. The taxonomic attribution of each function was resolved by matching the unique ID of each annotated sequence using the SQLite3 database (Version 3.28.0). Briefly, functionally annotated sequences were matched with the primary output of Kraken using SQLite3. The full SEED functional ontology of these sequences was obtained by matching their accession IDs to a reference database of the SEED subsystem.

Species diversity and richness were interrogated using the Chao 1 diversity index and abundance coverage estimator (ACE), and group-wise significance interrogated with rank-based nonparametric tests (Wilcoxon rank-sum test and Kruskall-Wallis test (using PhyloToAST diversity.py [25])). Dissimilarity between samples was calculated using Bray-Curtis metrics to estimate beta diversity. Nonmetric multidimensional scaling (NMDS) of Bray-Curtis distances was performed using the R package Vegan [24], and the R package Plotly used for visualization [26]. Significance of clustering was interrogated using adonis function, (QIIME implementation of vegan package in R) with 999 permutations. Similarity percentage analysis (SIMPER) of Bray-Curtis distances was used to determine the drivers of separation [27].

k-means clustering was used to estimate the variance in the disease-associated microbiome. Bray-Curtis dissimilarity (computed from the relative abundances of genes and species in subjects with disease) was used as input and silhouette width used to estimate the number of clusters [27]. A silhouette was created for each cluster based on the closeness of the members and separation from others. The silhouettes were then combined onto a single plot, providing an evaluation of the validity of the number of clusters. Factors that contributed to variance in the microbiome were identified using a between-class analysis method. Briefly, relative abundances of genes and species in subjects with disease were input into a principal coordinate analysis, and the eigenvectors that explained 75% of the variance in the microbiome were inducted into a linear discriminant analysis. Linear discriminant analysis (LDA) for dimensionality reduction of CSS-normalized s-OTU counts was performed using scikit-learn v0.18.0 [28] and plots were visualized using PhyloToAST. MANOVA/Wilks lambda was used to test for significance of LDA clustering.

Between-group differences in abundances and prevalence were estimated for both phylogenetic and functional data. Statistical testing for differences in relative abundances was performed using the DESeq2 package with Bayesian shrinkage of estimators (R package apeglm) [29, 30]. p-values were adjusted for multiple testing (FDR < 0.1, FDR-adjusted Wald test). Barycentric plots based on the functions that passed DESeq2 adjustments were graphed using the ggtern package in R [31]. Comparisons of functional potential between groups were made in the context of the KEGG (Kyoto Encyclopedia of Genes and Genomes) [32] and the SEED [33] ontological hierarchies. Presence/absence of the features was interrogated using Fisher’s exact test and odds ratios calculated with the R packages Questionr and DescTools.

Core features were calculated based on presence of the feature in at least 80% of patients with the particular condition. The SparCC package in python was used to construct network dependencies between features [34]. To decrease the occurrence of spurious associations due to rare taxa, co-occurrence networks were computed only on the core taxa [35]. Correlation was estimated on log-transformed abundances of the core features, and p-values computed following 99 bootstraps. Significant co-occurrences (rho ≥ 0.6 and p < 0.05) were used to create the graph structures. Gephi v0.9.1 [36] was used to visualize the resultant networks. The Zi-Pi plot calculations were done using the formula of Guimerà and Nunes Amaral and graphed using the simplified method of Olesen [37, 38].

The ability of genes to discriminate between groups was examined using a machine-learning algorithm (RandomForest package in R). Two thirds of the dataset was used to train the classifier, which was tested on the remaining data (1000 trees/10-fold cross validation). For all iterations of the test a “confusion table” was created for each of the exposures based on the number of correctly classified and misclassified samples, and this data was used to compute sensitivity and specificity. The robustness of the classifier was evaluated using ROC curves (ROCR package in R).

Results

Our primary aim was to characterize the functional and phylogenetic profiles of the subgingival microbiome in individuals with periodontitis. To do this, we obtained 72 million classifiable sequences from deep periodontal pockets of 59 systemically healthy subjects with periodontitis, and 31 million sequences from 25 periodontally healthy individuals. All subjects with periodontitis were classified as stage 3 based on disease severity and complexity (Table 1), with 17 patients demonstrating the molar-incisor phenotype (equivalent to LAP phenotype). Four subjects were classified as stage 3 grade A (S3gA), 22 as S3gB, and 33 as S3gC. Twenty-five subjects were classified as chronic periodontitis (CP; age range, 56–61 years), 17 as generalized aggressive periodontitis (GAP; age range, 24–32 years), and 17 as localized aggressive periodontitis (LAP; age range, 15–19 years). Their sequences represented 8336 functionally annotated microbial genes and 454 taxa.

Table 1 Clinical and demographic characteristics of the study population

Subgingival microbiomes follow the Anna Karenina principle

We began our analysis by creating a catalog of disease- and health-associated genes. In all three disease phenotypes, 73.8% of 454 taxa and 60% of 8336 genes were identified in 2 or more individuals within each group. Disease-associated microbiomes demonstrated greater beta-dispersion [39], with only 47% of disease-associated metagenome being shared by 80% or more of individuals with periodontitis (common core metagenome). On the other hand, over 73% of genes were identified in the common core of periodontal health, indicating that periodontitis-associated microbiomes follow the Anna Karenina principle (AKP). However, AKP was most evident in chronic periodontitis (34% of transcripts contributing to common core metagenome), followed by generalized aggressive periodontitis (51%) and localized aggressive periodontitis (62%).

We then queried whether AKP would preclude robust differentiation between health and disease. To do this, we computed a “consensus disease profile” based on genes and taxa found in all 59 patients and compared this to health using nonmetric multidimensional scaling of Bray-Curtis pairwise dissimilarities. The difference between health and disease was significantly greater than differences between any two individuals with disease (Fig. 1a, b, p < 0.001, PERMANOVA). A Random Forest Classifier (RandomForest package in R) classifier was able to predict disease with 87% sensitivity and 91% specificity based on functional profiles and with 72% sensitivity and 78% specificity when using phylogenetic metrics (Fig. 1c, d). Overall, 28% of genes were uniquely observed in disease, and 12% were unique to health, while 26% exhibited significant differential abundances in health and disease.

Fig. 1
figure 1

Between-class analysis of phylogenetic and functional profiles in health and disease. Nonmetric multidimensional scaling (NMDS) and receiver operating characteristic (ROC) curves shown. The first three dimensions of species-level (a) and gene-level (b) Bray-Curtis distances are shown (p-value < 0.001 for a and b). Each purple circle represents one of 59 subjects with periodontitis and each yellow circle represents one of 25 periodontally healthy individuals. The ability of disease-specific and health-specific indicators to predict each state is shown in c (phylogenetic metrics) and d (functional metrics)

Disease phenotype explains microbiome variance better than disease grade

Having established that periodontitis differed significantly from health both taxonomically and functionally, we next investigated if differences could be discerned within the periodontitis-associated microbiome using unsupervised cluster analysis (k-means clustering). Bray-Curtis dissimilarity distances (computed from the relative abundances of genes and species in subjects with disease) were used as input and silhouette width used to estimate number of clusters. We identified three distinct clusters taxonomically and functionally (p = 0.0008 and 0.001 respectively, ADONIS test of Bray-Curtis Dissimilarity Index, Fig. 2a, b). We then investigated the factors that drove these differences using a between-class analysis method that combines principal coordinates analysis with linear discriminant analysis (Fig. 2c–j). Relative abundances of genes and species in subjects with disease were used as input. Disease phenotype, ethnicity, and age emerged as the strongest drivers of clustering. Disease phenotype yielded the lowest degree of misclassification both taxonomically and functionally while significant misclassification was evident when using disease grade as a discriminant. Furthermore, disease phenotype explained the strong ethnicity and age-based clustering, since most of the younger individuals and those of African American ethnicity belonged to the LAP group (Fig. 2i).

Fig. 2
figure 2

Factors that explain variance in the subgingival microbiome. k-means clustering of Bray-Curtis distances revealed three clusters (LAP blue, CP red, GAP green) based on function (a, p-value < 0.001) and taxa (b, p-value < 0.001). The ellipses represent the centroids of each cluster. Linear discriminant analysis of relative abundances of function (top panels) and phylogeny (bottom panels) revealed that disease phenotype (c, d), but not disease grade (e, f), was able to discriminate between subjects with disease. The differences based on age-decade (g, h) and ethnicity (i, j) were explained by the disease phenotype. In each cluster, the larger ellipse indicates the 95% confidence region to contain the true mean for the group, and the smaller (inner) ellipse represents the region estimated to contain 50% of the population for the group. The misclassification rates of each variable are shown within each panel

Same players, different teams

Since targeted microbial investigations have previously suggested that localized aggressive periodontitis has a distinct microbial profile while the microbiota of chronic and generalized aggressive periodontitis are similar [40], we tested the hypothesis that GAP and CP are microbially similar while LAP is a taxonomically distinct entity. NMDS revealed significant class separation between the 3 diseases (p < 0.001, PERMANOVA of Bray-Curtis Dissimilarity Index, Fig. 2b). We then investigated whether GAP and CP are microbially more similar than GAP and LAP by computing pairwise dissimilarities (Bray-Curtis) between each GAP and CP subject, as well as each GAP and LAP subject. Taxonomically, GAP was more similar to LAP than to CP (p < 0.001, Dunn’s test for multiple comparisons).

Since the NMDS indicated that these three diseases were microbially discrete entities, we examined the taxonomic features that contributed to class separation. LAP exhibited significantly lower species richness (as measured by the ACE and Chao 1 indices) than the other groups (p < 0.05, Dunn’s test, Fig. 3a), and the three groups also demonstrated significant differences in beta diversity. While all 3 diseases were dominated by gram-negative anaerobic bacteria, (representing 56.8%, 62.5%, and 47.9% of the abundance in CP, GAP, and LAP, respectively (Fig. 3b)), the abundances of these groups were significantly greater in GAP when compared to LAP (p = 0.03, Wilcoxon nonparametric test). By contrast, gram-positive anaerobic bacteria were significantly higher in CP when compared to either GAP or LAP (p < 0.04, Wilcoxon).

Fig. 3
figure 3

Disease-specific taxonomical indicators. Density curves of alpha diversity (ACE) are shown in a. The peak indicates the median values for each group, and the x-axis shows the data range. LAP exhibited significantly lower alpha diversity than the other 2 groups (p < 0.0001, Dunn’s test). Distribution of species-level taxa by gram staining characteristics and oxygen requirements in is shown in b. GAP patients demonstrated significantly greater gram-negative anaerobic bacteria and lower gram-negative aerobic bacteria when compared to the other two groups (p < 0.01, Dunn’s test for multiple comparisons). Percent of the microbiome that is shared by 80% or more of individuals (common core microbiome) with CP, GAP, and LAP are graphically indicated by the Euler graphs c (i–iii), as well as the number of core species shared by all three diseases is shown in c (iv). Phylogenetic tree of species that were significantly different between groups (p < 0.05, FDR-adjusted Wald test) are shown in 3D. Bars represent normalized mean relative abundances, while the solid circles indicate species that belong to the common disease core. Data supporting this figure can be found in Supplemental Table 1

Interestingly, 349 out of 416 species were identified in all the three diseases, and only 28 species were unique to any one of the three diseases. Collectively, the unique species constituted less than 0.03% of the abundance in each group. Each disease condition had core taxa that constitute more than 50% of the identified taxa in the condition (Fig. 3c (i–iii)). When core taxa of each condition were compared, most species were present in the cores of all 3 conditions (Fig. 3c (iv)). 138 OTUs were found to be significantly differentially abundant between any two disease states (p < 0.05, FDR-adjusted Wald test — Fig. 3d and Supplementary Table 1). Using similarity percentages (SIMPER) analysis, we identified OTUs that explained 70% or more of the class separation. Several OTUs which significantly contributed to the separation were also common core taxa, demonstrating that the diseases differ in the ratios of their predominant shared taxa. One hundred and seven OTUs were significant contributors to the separation between GAP and LAP; of these, 65 were part of the common core species of GAP and 63 formed the common core of LAP (Fig. 3d, supplementary Table 1). Similarly, 50 species, 39 of which were members of the core microbiomes of GAP and CP, contributed to the separation between GAP and CP. The separation between CP and GAP was driven by Aggregatibacter actinomycetemcomitans, Fusobacterium nucleatum, Treponema socranskii, and several members of the genera Actinomyces, Campylobacter, Prevotella, and Capnocytophaga. The separation between LAP and GAP was mainly driven through Porphyromonas gingivalis and members of the genera Neisseria and Actinomyces.

Since inter-bacterial interactions play a large role in influencing microbial assemblages, we used graph theoretics to assess connectivity between species. The underlying rationale for this analysis is that taxa with the strongest connections demonstrate superior adaptation to their niche. As a corollary, diseases that present similar microenvironments will demonstrate greater co-dependency between member species than diseases that are different. The network topography is summarized in Supplemental Table 2. While GAP and LAP demonstrated robust hubs with 3568 and 2114 edges, CP presented a sparse topography, with only 489 connections, attesting to its phylogenetically idiosyncratic presentation (Fig. 4). Zi-Pi plots of both CP and LAP demonstrated expansive nodes with several putative keystone species in the network topography, while the node distribution in GAP was equitable (preventing us from creating a Zi-Pi plot) and did not demonstrate any candidate keystone species. Together, the data suggest that patients with GAP and LAP have a more homogeneous subgingival microenvironment than those with CP, which may explain the taxonomic heterogeneity observed in CP. Based on the clinical observation that 35% of untreated cases of LAP progress to GAP [41], we hypothesize that loss of the influential key players found in LAP creates a state of flux that, when observed cross-sectionally, gives rise to the observation that GAP is a distinct disease phenotype. This theory is further supported by observations that individuals with GAP demonstrate low serum antibody response to the microbial constituents, leading to its continuous periodontal destruction [40]. This is unlike the other two phenotypes which can self-arrest with time.

Fig. 4
figure 4

Same players, different teams. Zi-Pi plots of co-occurrence networks in generalized aggressive periodontitis (a), chronic periodontitis (b), and localized aggressive periodontitis (c). The text box illustrates the general topographical characteristics of each network. The plots show nodes with low connectivity (peripherals), and the putative keystone species, which are separated into local/module key players (module hubs), key players across different modules (connectors), and the combination of the two categories (network hubs). Network graphs of chronic periodontitis and localized aggressive periodontitis (d, e) show the ego network of species that have first-degree edges with the putative keystone species. A similar network graph was not possible for generalized aggressive periodontitis since there no keystone pathogens were identified in the Zi-Pi plot. Nodes are colored based on their modules. Thick edges represent positive correlations, and faint lines represent negative correlations. Data supporting this figure can be found in Supplemental Table 2

The LAP microbiome is functionally distinguishable from CP and GAP

Since the three disease phenotypes demonstrated several taxa in common, we tested the hypothesis that there would be significant functional overlap in their respective associated microbiomes using the SEED ontology to annotate genes and the KEGG database for pathway identification. A greater degree of class separation was evident based on functional capabilities than on taxonomic profiles (p < 0.0001, PERMANOVA of Bray-Curtis Dissimilarity, Fig. 5a). 20.61% of the LAP metagenome (1278/6200 genes) was unique, in that, these genes were not present in either GAP or CP (Fig. 5b and Supplemental Table 3). Twenty percent of these unique genes did not have functional role assignments, pointing to gaps in our knowledge of the microbiome of localized aggressive periodontitis. Forty percent of unique genes encoded enzymes for anaerobic degradation of aromatic compounds, methanogenesis, lysine and acetyl CoA fermentation, and anaerobic respiratory reductases. Twenty percent of the unique genes coded for gram-negative cell structures and 5% for gram-negative phages. 27.96% (2037/7286) of the LAP metagenome differed significantly from that of GAP (p < 0.05, FDR-adjusted Wald test, DESeq2). The LAP microbiome demonstrated greater capacity for inositol catabolism, and Lipid A, lipopolysaccharide, and peptidoglycan biosynthesis when compared to GAP (p < 0.05, FDR-adjusted Wald test, DESeq2). Additionally, the LAP biome demonstrated a 4-fold to 144-fold greater enrichment of genes encoding c-type cytochrome and molybdenum cofactor biosynthesis, iron-sulfur clusters, formate dehydrogenase, and oxidative stress response. 27.5% (2019/7340) of the LAP metagenome differed from CP. This was attributable to a higher representation of genes encoding acetyl CoA, lactate, mixed-acid and lysine fermentation, methanogenesis, anaerobic respiratory reductases, dehydrogenases, dehydratases and anaerobic toluene, and ethylbenzene degradation. Also overrepresented in the LAP metagenome were membrane transport functions (type II, III, IV, V, and VI secretions systems and ABC transporters), and functions related to quorum sensing and biofilm formation (Autoinducer-2 transport and processing, biofilm adhesins, and histidine kinase sensors).

Fig. 5
figure 5

Disease-specific functional indicators. Barycentric plots of significantly different virulence functions in the three groups (p < 0.05, FDR adjusted Wald test) are shown. Each dot represents a gene. The three groups generalized aggressive periodontitis, chronic periodontitis, and localized aggressive periodontitis (GAP, CP, and GAP) are used as vertices. Within each plot, the coordinates of each gene are determined by the weighted average of the coordinates of all genes, and the weights are given by the relative abundance of the gene in that group (LAP, GAP, and CP). Data supporting this figure can be found in Supplemental Table 3

The GAP microbiome — a functional chimera

Pairwise dissimilarity analysis revealed that the GAP metagenome was intermediate between CP and LAP (p = 0.91, Dunn’s multiple comparisons test on Bray-Curtis distances between GAP-CP and GAP-LAP). GAP shared 77% of its metagenome with CP and 64% of its metagenome with LAP. To reduce bias induced by sparse data, we next used core genes in each group to compute these distances. Not surprisingly, we observed a lower similarity within the common core metagenomes, but a more balanced difference, with 56% of genes shared by CP and GAP, and 55.1% by GAP and LAP. Since severe attachment loss in young adults in the presence of clinical inflammation and local factors could represent either early-onset chronic periodontitis, a generalized form of the molar-incisor phenotype, or true de novo aggressive periodontitis [41], we examined clustering of the GAP samples alone. NMDS did not reveal significant separation between the 17 samples, suggesting that this chimeric effect cannot be readily attributed to heterogenous diseases.

Fourteen percent of the core genes shared by LAP and GAP encoded as yet unknown functions (Supplemental Table 4). Among the characterized genes, the predominant shared functionalities were related to an anaerobic lifestyle. These included genes encoding heme- and hemin-dependent respiration, dehydrogenases, electron donors and acceptors other than oxygen (namely, nitrate, sulfate, hydrogen, and ferric iron), and fermentation. Other shared functions included polyamine metabolism, flagellar biosynthesis, and gram-negative cell wall components (including peptidoglycan biosynthesis), response to oxidative and osmotic stress, resistance to antibiotics and toxic compounds, phages, and conjugative transposons. The differences in the microbiomes of LAP and GAP were attributable to lower abundances of membrane transport functions (type II, III, IV, V VI, and ABC transporters), quorum sensing and biofilm formation in and higher levels of sporulation and dormancy, phages and transposable elements in GAP. Pathways involved in biofilm stability were also lower in GAP in comparison to LAP, and even lower in CP when compared to GAP.

The functional roles of 19% of the genes shared by CP and GAP were unknown. Both CP and GAP demonstrated a collective capacity for metabolism of amino acid, organic compounds, alcohols, and glycogen. The other shared functionalities included capsule and cell wall synthesis, response to oxidative and osmotic stress, and resistance to antibiotics and toxic compounds, phages, and conjugative transposons. Flagellar components and proteins associated with flagella biosynthesis and assembly machinery are also exclusively enriched in GAP as compared to CP, as were potent inflammatory triggers such as lipopolysaccharides and peptidoglycans. Other abundant functions of GAP include dormancy and sporulation, invasion and intracellular resistance, iron acquisition and siderophores, and multidrug antibiotics efflux pumps.

Discussion

Culture-based approaches to characterizing the subgingival microbiome spawned a slew of studies on individual bacteria, their responses to environmental shifts, and their roles in disease causation as independent operators [42,43,44,45,46]. However, cultivation independent methods have revealed that polymicrobial diseases are caused by the collective actions and interactions of the microbiome [47, 48], and that the “invisible” members of a microbial community have important contributions to these interactions [49, 50]. When we used a top-down approach to map the genomic content of and the biological pathways encoded by the microbiomes associated with chronic, localized, and generalized aggressive periodontitis, a surprising picture of the three phenotypes emerged.

Healthy microbiomes demonstrated significantly lower beta-dispersion than disease, similar to what has been reported in the gut microbiome [51]. This effect has been called the Anna Karenina principle and suggests that healthy microbiomes operate under stricter host control than do diseased ones. Subjects with periodontitis were carefully selected for clinical homogeneity: all sampled sites demonstrated similar attachment loss, pocket depths, and periodontal inflammation. Therefore, the variance in the microbiome could not be explained by disease severity. This dispersion was more apparent in taxonomy than in the functional profiles, corroborating our earlier finding that the microbiome associated with periodontitis is taxonomically heterogenous but functionally congruent [20]. This was also corroborated by identifying a set of putative periodontal pathogens and perturbed functions common to all three diseases, suggesting that certain genes and taxa are broadly associated with the disease process. This is not surprising, since it has long been recognized that periodontitis is a microbially heterogenous disease [52], this observation was largely based on examination of selected species or meta-taxonomic cross-sectional studies. However, by overlaying robust clinical metrics on the genomic content of the subgingival microbiome, we were able to identify that not all clinical phenotypes of periodontitis are equally heterogenous and that the heterogeneity does not extend to microbial functions. We realize, as we say this, that this is a cross-sectional study, and therefore, it is not designed to identify patterns in microbial community dynamics. Our findings serve to further reinforce the critical need for adequately powered longitudinal studies that combine granular clinical metadata with open-ended explorations of gene expression patterns and protein-protein networks, and rigorous and improved modeling of stochastic events in microbiome assembly.

Another observation was that most functions which discriminated between health and disease were significantly more abundant in disease, suggesting that the disease environment is associated with expanded functionality. Our study reveals that particular functional capability is required for life in the anaerobic, pro-oxidant, heme-rich environment of a pathologically deepened subgingival sulcus. For example, lipopolysaccharide biosynthesis, iron transport, stress response, fermentation, and metabolism of secondary amino acids were uniformly enriched in all three diseases when compared to health, even though some of them were more abundant in one of the three disease phenotypes.

We also found that species demonstrating a fitness for this environment are different in each individual, highlighting challenges associated with developing species-based biomarkers or vaccines for periodontitis. For example, flagellar proteins, potent inflammatory triggers, mapped to a variety of organisms (e.g., species belonging to the genera Treponema, Selenomonas, Pseudomonas, and Campylobacter), while genes encoding dormancy and sporulation, invasion and intracellular resistance, iron acquisition and siderophores, multidrug antibiotics efflux pumps were assigned to several species within the genera Aggregatibacter, Porphyromonas, Atopobium, and Prevotella, among others.

Integrating data from the three most common phenotypes of periodontitis allowed for identification of phenotype-specific indicators, elucidating their potential role in disease causation. Disease-specific indicators were more readily evident in the LAP microbiome than GAP or CP. When we examined the unique and overrepresented suite of genes in LAP, a picture emerged of a community with greater fitness for an anaerobic, proteolytic lifestyle, and a higher capacity for virulence than CP or GAP. For example, genes that encode the “acetate switch”, permitting a transition from rapid growth to a slower, acetate-scavenging lifestyle [53], as well as dehydrogenases, hydratases, and anaerobic reductases were identified only in LAP. Anaerobic glycerol phosphate-3-hydrogenase genes, which play a critical role in utilizing alternate nutritional sources [54] were consistently higher in LAP when compared to the other two. Genes encoding c-type cytochrome and molybdenum cofactor biosynthesis, iron-sulfur clusters, and formate dehydrogenase were also overrepresented in this cohort. Molybdenum cofactor is essential in bacterial respiration and energy conversion, especially in those species that do not have appreciable plasticity in their metabolic and respiratory pathways [55]. Molybdenum metabolism is also tightly connected to iron-sulfur cluster synthesis [56]. Formate is not only an important byproduct of anaerobic respiration, it is also a substrate for many sulfate reducing bacteria, for example, Campylobacter and Prevotella [57]. Its synthesis is mediated by formate dehydrogenase, an enzyme that contains molybdenum cofactor and iron-sulfur clusters. C-type cytochromes regulate several key pathogenic processes, notably heme synthesis, oxidative stress response, nitrosative stress response, and siderophore production [58]. This small set of genes that were unique to LAP may provide insights into disease etiology and development of LAP-specific microbial biomarkers and deserve further investigation.

On the other hand, the microbiome of GAP shared important taxonomic and functional features with both LAP and CP. This “chimera”-like appearance might, in part, explain the clinical observation that 35% of untreated cases of LAP progress to GAP [41]. Furthermore, the absence of the putative keystone species in GAP could indicate an ecosystem in a state of flux, which can explain the poor immune response from the patients, and the lack of the ability of the disease to self-arrest [59]. While the absence of clustering within the GAP samples argues against this, it is possible that our small sample size precluded a statistically significant clustering. Our observations lead us to question whether GAP is indeed a unique disease or represents, for some individuals at least, a “halfway house” between LAP and CP.

We noticed progressively greater beta-dispersion and the concomitantly smaller common core microbiomes when moving from LAP to GAP to CP, pointing to greater personalization of the disease-associated microbiome in older individuals with CP (CP average age in this cohort 58 ± 2 years old). Two possible explanations present themselves: one is that the microbiome naturally shifts with increasing age; and the age of the microbiome when dysbiosis sets in could determine the disease-associated profile. The other is that personalization occurs due to the chronicity of the disease, following the Anna Karenina principle (AKP).

We also noted an attenuation of virulence-related functionalities and stress response from LAP to GAP to CP as evidenced by 4- to 144-fold lower levels of gram-negative cell wall components (including LPS and Lipid A) and fermentation, and of genes in the glutathione pathway in GAP when compared to LAP and in CP when compared to GAP. This serves to explain how the microbiome mediates the differential periodontal destruction observed in LAP, GAP, and CP. Virulence genes play an important role in the creation of a dysbiotic ecosystem, since they allow the species expressing them to overcome the colonization resistance offered by the health-compatible indigenous species. However, virulence comes at an extreme fitness cost, and therefore, bacteria benefit by staying avirulent unless environmental cues or competition dictate it. Our data suggests that attenuation of virulence potential might be a factor contributing to chronicity of the disease. This might explain the random burst model of disease progression [60], where an increase in community virulence precedes bursts of inflammation and loss of tooth-supporting structures.

Conclusions

In summary, a comprehensive metagenomic analysis of the subgingival microbiomes in different disease phenotypes reveals broad patterns of shift in microbial functions that span all diseases. Many of these functions facilitate life in an oxygen-poor, protein- and heme-rich, pro-oxidant environment, as well as providing an enhanced capacity for attachment and biofilm formation. However, beta-dispersion metrics demonstrate that no two individuals with disease are alike, especially older individuals with chronic disease phenotype. Therapies focused on microbial modulation through mechanical, chemical, or other means will have to take into account patient-specific parameters for efficacy. Importantly, we observe that generalized aggressive periodontitis shares significant functional features with both localized aggressive periodontitis and chronic periodontitis, suggesting either attenuation of an aggressive disease or an early-onset chronic disease. We therefore question whether this is a separate disease entity, or an artifact induced by cross-sectional study designs. The present investigation also uncovers disease-specific indicators with varying discriminant abilities for each phenotype. These can not only serve as potential biomarkers for molecular identification of clinical phenotypes, but also clarify the role of the microbiome in disease pathogenesis.

Availability of data and materials

Sequences for all 59 samples are deposited in the Sequence Read Archives under the project ID PRJNA552294 and PRJNA508385.

References

  1. Kassebaum NJ, Bernabe E, Dahiya M, Bhandari B, Murray CJ, Marcenes W. Global burden of severe periodontitis in 1990-2010: a systematic review and meta-regression. J Dent Res. 2014;93(11):1045–53. https://doi.org/10.1177/0022034514552491.

    Article  Google Scholar 

  2. Mohd-Dom T, Ayob R, Mohd-Nur A, Abdul-Manaf MR, Ishak N, Abdul-Muttalib K, et al. Cost analysis of periodontitis management in public sector specialist dental clinics. BMC Oral Health. 2014;14(1):56. https://doi.org/10.1186/1472-6831-14-56.

    Article  Google Scholar 

  3. Fardal O, O'Neill C, Gjermo P, Fardal E, Sandvik L, Hansen BF, et al. The lifetime direct cost of periodontal treatment: a case study from a Norwegian specialist practice. J Periodontol. 2012;83(12):1455–62. https://doi.org/10.1902/jop.2012.110689.

    Article  Google Scholar 

  4. Gjermo PE, Grytten J. Cost-effectiveness of various treatment modalities for adult chronic periodontitis. Periodontol. 2009;51:269–75. https://doi.org/10.1111/j.1600-0757.2009.00313.x.

    Article  Google Scholar 

  5. Kumar PS. From focal sepsis to periodontal medicine: a century of exploring the role of the oral microbiome in systemic disease. J Physiol. 2016;595(2):465–76. https://doi.org/10.1113/JP272427.

    Article  Google Scholar 

  6. Armitage GC. Development of a classification system for periodontal diseases and conditions. Ann Periodontol. 1999;4(1):1–6. https://doi.org/10.1902/annals.1999.4.1.1.

    Article  Google Scholar 

  7. Demmer RT, Papapanou PN. Epidemiologic patterns of chronic and aggressive periodontitis. Periodontol. 2010;53(1):28–44.

    Article  Google Scholar 

  8. Faveri M, Mayer MP, Feres M, de Figueiredo LC, Dewhirst FE, Paster BJ. Microbiological diversity of generalized aggressive periodontitis by 16S rRNA clonal analysis. Oral Microbiol Immunol. 2008;23(2):112–8. https://doi.org/10.1111/j.1399-302X.2007.00397.x.

    Article  Google Scholar 

  9. Fine DH, Patil AG, Loos BG. Classification and diagnosis of aggressive periodontitis. J Clin Periodontol. 2018;45(S20):S95–S111. https://doi.org/10.1111/jcpe.12942.

    Article  Google Scholar 

  10. Haubek D, Ennibi OK, Poulsen K, Poulsen S, Benzarti N, Kilian M. Early-onset periodontitis in Morocco is associated with the highly leukotoxic clone of Actinobacillus actinomycetemcomitans. J Dent Res. 2001;80(6):1580–3. https://doi.org/10.1177/00220345010800062001.

    Article  Google Scholar 

  11. Shaddox LM, Huang H, Lin T, Hou W, Harrison PL, Aukhil I, et al. Microbiological characterization in children with aggressive periodontitis. J Dent Res. 2012;91(10):927–33. https://doi.org/10.1177/0022034512456039.

    Article  Google Scholar 

  12. Socransky SS, Haffajee AD, Cugini MA, Smith C, Kent RL Jr. Microbial complexes in subgingival plaque. J Clin Periodontol. 1998;25(2):134–44. https://doi.org/10.1111/j.1600-051X.1998.tb02419.x.

    Article  Google Scholar 

  13. Abusleme L, Dupuy AK, Dutzan N, Silva N, Burleson JA, Strausbaugh LD, et al. The subgingival microbiome in health and periodontitis and its relationship with community biomass and inflammation. ISME J. 2013;7(5):1016–25. https://doi.org/10.1038/ismej.2012.174.

    Article  Google Scholar 

  14. Li Y, He J, He Z, Zhou Y, Yuan M, Xu X, et al. Phylogenetic and functional gene structure shifts of the oral microbiomes in periodontitis patients. ISME J. 2014;8(9):1879–91. https://doi.org/10.1038/ismej.2014.28.

    Article  Google Scholar 

  15. Ge X, Rodriguez R, Trinh M, Gunsolley J, Xu P. Oral microbiome of deep and shallow dental pockets in chronic periodontitis. PLoS ONE. 2013;8(6):e65520. https://doi.org/10.1371/journal.pone.0065520.

    Article  Google Scholar 

  16. Picolos DK, Lerche-Sehm J, Abron A, Fine JB, Papapanou PN. Infection patterns in chronic and aggressive periodontitis. J Clin Periodontol. 2005;32(10):1055–61. https://doi.org/10.1111/j.1600-051X.2005.00828.x.

    Article  Google Scholar 

  17. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60. https://doi.org/10.1186/gb-2011-12-6-r60.

    Article  Google Scholar 

  18. Zaneveld JR, McMinds R, Vega TR. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat Microbiol. 2017;2(9):17121. https://doi.org/10.1038/nmicrobiol.2017.121.

    Article  Google Scholar 

  19. Papapanou PN, Sanz M, Buduneli N, Dietrich T, Feres M, Fine DH, et al. Periodontitis: consensus report of workgroup 2 of the 2017 World Workshop on the Classification of Periodontal and Peri-Implant Diseases and Conditions. J Periodontol. 2018;89(Suppl 1):S173–s82. https://doi.org/10.1002/jper.17-0721.

    Article  Google Scholar 

  20. Dabdoub SM, Ganesan SM, Kumar PS. Comparative metagenomics reveals taxonomically idiosyncratic yet functionally congruent communities in periodontitis. Sci Rep. 2016;6(1):38993. https://doi.org/10.1038/srep38993.

    Article  Google Scholar 

  21. Wilke A, Bischof J, Gerlach W, Glass E, Harrison T, Keegan KP, et al. The MG-RAST metagenomics database and portal in 2015. Nucleic Acids Res. 2016;44(D1):D590–4. https://doi.org/10.1093/nar/gkv1322.

    Article  Google Scholar 

  22. Tang W, Wilkening J, Bischof J, Gerlach W, Wilke A, Desai N, et al. Building scalable data management and analysis infrastructure for metagenomics. 2013.

    Google Scholar 

  23. Wilke A, Bischof J, Harrison T, Brettin T, D'Souza M, Gerlach W, et al. A RESTful API for accessing microbial community data for MG-RAST. PLoS Comput Biol. 2015;11(1):e1004008. https://doi.org/10.1371/journal.pcbi.1004008.

    Article  Google Scholar 

  24. Oksanen J, Kindt R, Legendre P, O'Hara B, Simpson GL, Solymos P, Stevens MHH, Wagner H. vegan: Community Ecology Package (R package version 1.15-1); 2008.

  25. Dabdoub SM, Fellows ML, Paropkari AD, Mason MR, Huja SS, Tsigarida AA, et al. PhyloToAST: bioinformatics tools for species-level analysis and visualization of complex microbial datasets. Sci Rep. 2016;6(1):29123. https://doi.org/10.1038/srep29123.

    Article  Google Scholar 

  26. Sievert C. Interactive web-based data visualization with R, plotly, and shiny. Boca Raton: CRC Press, Taylor and Francis Group; 2020. https://doi.org/10.1201/9780429447273.

    Book  Google Scholar 

  27. Rousseeuw PJ. Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. J Comput Appl Math. 1987;20:53–65 https://doi.org/10.1016/0377-0427(87)90125-7.

    Article  Google Scholar 

  28. Pedregosa F, Ga, #235, Varoquaux L, Gramfort A, Michel V, et al. Scikit-learn: machine learning in Python. J Mach Learn Res. 2011;12:2825–30.

    Google Scholar 

  29. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.

    Article  Google Scholar 

  30. Zhu A, Ibrahim JG, Love MI. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics. 2018;35(12):2084–92. https://doi.org/10.1093/bioinformatics/bty895.

    Article  Google Scholar 

  31. Hamilton NE, Ferry M. ggtern: ternary diagrams using ggplot2. J Stat Softw. 2018;87(1):1–17.

    Google Scholar 

  32. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. https://doi.org/10.1093/nar/28.1.27.

    Article  Google Scholar 

  33. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014;42(Database issue):D206–14. https://doi.org/10.1093/nar/gkt1226.

    Article  Google Scholar 

  34. Friedman J, Alm EJ. Inferring correlation networks from genomic survey data. PLoS Comput Biol. 2012;8(9):e1002687. https://doi.org/10.1371/journal.pcbi.1002687.

    Article  Google Scholar 

  35. Weiss S, Van Treuren W, Lozupone C, Faust K, Friedman J, Deng Y, et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. Isme J. 2016;10(7):1669–81. https://doi.org/10.1038/ismej.2015.235.

    Article  Google Scholar 

  36. Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. ICWSM. 2009;8:361–2.

    Article  Google Scholar 

  37. Guimera R, Nunes Amaral LA. Functional cartography of complex metabolic networks. Nature. 2005;433(7028):895–900. https://doi.org/10.1038/nature03288.

    Article  Google Scholar 

  38. Olesen JM, Bascompte J, Dupont YL, Jordano P. The modularity of pollination networks. Proc Natl Acad Sci. 2007;104(50):19891–6. https://doi.org/10.1073/pnas.0706375104.

    Article  Google Scholar 

  39. Anderson MJ. Distance-based tests for homogeneity of multivariate dispersions. Biometrics. 2006;62(1):245–53. https://doi.org/10.1111/j.1541-0420.2005.00440.x.

    Article  Google Scholar 

  40. Kononen E, Muller HP. Microbiology of aggressive periodontitis. Periodontol. 2014;65(1):46–78. https://doi.org/10.1111/prd.12016.

    Article  Google Scholar 

  41. Brown LJ, Albandar JM, Brunelle JA, Loe H. Early-onset periodontitis: progression of attachment loss during 6 years. J Periodontol. 1996;67(10):968–75. https://doi.org/10.1902/jop.1996.67.10.968.

    Article  Google Scholar 

  42. Lamont RJ, El-Sabaeny A, Park Y, Cook GS, Costerton JW, Demuth DR. Role of the Streptococcus gordonii SspB protein in the development of Porphyromonas gingivalis biofilms on streptococcal substrates. Microbiology. 2002;148(Pt 6):1627–36.

    Article  Google Scholar 

  43. Marsh PD, McDermid AS, McKee AS, Baskerville A. The effect of growth rate and haemin on the virulence and proteolytic activity of Porphyromonas gingivalis W50. Microbiology. 1994;140(Pt 4):861–5.

    Article  Google Scholar 

  44. Michalowicz BS, Ronderos M, Camara-Silva R, Contreras A, Slots J. Human herpesviruses and Porphyromonas gingivalis are associated with juvenile periodontitis. J Periodontol. 2000;71(6):981–8. https://doi.org/10.1902/jop.2000.71.6.981.

    Article  Google Scholar 

  45. Mikuls TR, Payne JB, Yu F, Thiele GM, Reynolds RJ, Cannon GW, et al. Periodontitis and Porphyromonas gingivalis in patients with rheumatoid arthritis. Arthritis Rheum. 2014;66(5):1090–100. https://doi.org/10.1002/art.38348.

    Article  Google Scholar 

  46. Pussinen PJ, Alfthan G, Jousilahti P, Paju S, Tuomilehto J. Systemic exposure to Porphyromonas gingivalis predicts incident stroke. Atherosclerosis. 2007;193(1):222–8.

  47. Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3. https://doi.org/10.1038/4441022a.

    Article  Google Scholar 

  48. Paster BJ, Boches SK, Galvin JL, Ericson RE, Lau CN, Levanos VA, et al. Bacterial diversity in human subgingival plaque. J Bacteriol. 2001;183(12):3770–83. https://doi.org/10.1128/JB.183.12.3770-3783.2001.

    Article  Google Scholar 

  49. Bernard G, Pathmanathan JS, Lannes R, Lopez P, Bapteste E. Microbial dark matter investigations: how microbial studies transform biological knowledge and empirically sketch a logic of scientific discovery. Genome Biol Evol. 2018;10(3):707–15. https://doi.org/10.1093/gbe/evy031.

    Article  Google Scholar 

  50. Solden L, Lloyd K, Wrighton K. The bright side of microbial dark matter: lessons learned from the uncultivated majority. Curr Opin Microbiol. 2016;31:217–26 https://doi.org/10.1016/j.mib.2016.04.020.

    Article  Google Scholar 

  51. Armour CR, Nayfach S, Pollard KS, Sharpton TJ. A metagenomic meta-analysis reveals functional signatures of health and disease in the human gut microbiome. mSystems. 2019;4(4):e00332–18. https://doi.org/10.1128/mSystems.00332-18.

    Article  Google Scholar 

  52. Teles R, Teles F, Frias-Lopez J, Paster B, Haffajee A. Lessons learned and unlearned in periodontal microbiology. Periodontol. 2013;62(1):95–162. https://doi.org/10.1111/prd.12010.

    Article  Google Scholar 

  53. Wolfe AJ. The acetate switch. Microbiol Mol Biol Rev. 2005;69(1):12–50. https://doi.org/10.1128/MMBR.69.1.12-50.2005.

    Article  Google Scholar 

  54. Daniels JB, Scoffield J, Woolnough JL, Silo-Suh L. Impact of glycerol-3-phosphate dehydrogenase on virulence factor production by Pseudomonas aeruginosa. Can J Microbiol. 2014;60(12):857–63. https://doi.org/10.1139/cjm-2014-0485.

    Article  Google Scholar 

  55. Dhouib R, Pg Othman DSM, Essilfie A-T, Hansbro PM, Hanson JO, McEwan AG, et al. Maturation of molybdoenzymes and its influence on the pathogenesis of non-typeable Haemophilus influenzae. Front Microbiol. 2015;6(1219). https://doi.org/10.3389/fmicb.2015.01219.

  56. Bittner F, Mendel RR. Molybdenum cofactor biosynthesis and cross talk with iron–sulfur. Encyclopedia of Inorganic and Bioinorganic Chemistry. 2013.

  57. Loesche WJ, Grossman NS. Periodontal disease as a specific, albeit chronic, infection: diagnosis and treatment. Clin Microbiol Rev. 2001;14(4):727–52, table of contents. https://doi.org/10.1128/CMR.14.4.727-752.2001.

    Article  Google Scholar 

  58. Wilson A. Bacterial c-type cytochromes and pathogenicity. Eur J BioMed Res. 2015;1(1):17–21. https://doi.org/10.18088/ejbmr.1.1.2015.pp17-21.

    Article  Google Scholar 

  59. Schenkein HA, Barbour SE, Tew JG. Cytokines and inflammatory factors regulating immunoglobulin production in aggressive periodontitis. Periodontol. 2007;45(1):113–27. https://doi.org/10.1111/j.1600-0757.2007.00214.x.

    Article  Google Scholar 

  60. Goodson JM, Tanner AC, Haffajee AD, Sornberger GC, Socransky SS. Patterns of progression and regression of advanced destructive periodontal disease. J Clin Periodontol. 1982;9(6):472–81. https://doi.org/10.1111/j.1600-051X.1982.tb02108.x.

    Article  Google Scholar 

Download references

Funding

This study was funded by R01-DE022579 awarded to Purnima Kumar by the National Institutes of Dental and Craniofacial Research.

Author information

Authors and Affiliations

Authors

Contributions

KA: sample collection, DNA processing, bioinformatics, statistical analysis, and manuscript preparation. PM: sample collection and manuscript preparation. SMG: sample collection, DNA processing, bioinformatics, statistical analysis, and manuscript preparation. SMD: bioinformatics, statistical analysis, and manuscript preparation. HNN: statistical analysis and manuscript preparation. PSK: overall study design, funding, sample collection, DNA processing, bioinformatics, statistical analysis, and manuscript preparation. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Purnima S. Kumar.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the institutional review boards of The Ohio State University and Louisiana State University (OSU IRB 2014H0020, LSUHS-NO 8796). Thirty-four nonsmoking, normoglycemic individuals with stage 3 periodontitis were recruited and informed consent or assent with parental approval was obtained as appropriate.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplemental Table 1.

Genus and species-level taxa that differed between chronic, localized and generalized aggressive periodontitis are shown. Differences in detection frequency (Fisher's exact test), differential abundance (DESeq2 along with log (2) fold differences, p-value, and FDR adjusted p-value (Wald’s Test), drivers of the class separation (SIMPER), species that were identified in ≥ 80% of individuals (common core microbiome) are shown.

Additional file 2: Supplemental Table 2.

Same players, different teams. Network statistics of co-occurrence plots are shown.

Additional file 3: Supplemental Table 3.

Periodontitis and the Anna Karenina principle. Genes that differed between chronic, localized and generalized aggressive periodontitis are shown, along with log (2) fold differences, p-value, and FDR adjusted p-value (Wald’s Test).

Additional file 4: Supplemental Table 4.

GAP-a microbial chimera. Genes that were shared by all three diseases, as well as between generalized and localized aggressive periodontitis and generalized aggressive periodontitis and chronic periodontitis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Altabtbaei, K., Maney, P., Ganesan, S.M. et al. Anna Karenina and the subgingival microbiome associated with periodontitis. Microbiome 9, 97 (2021). https://doi.org/10.1186/s40168-021-01056-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-021-01056-3

Keywords