- Open Access
Subgingival microbiome of deep and shallow periodontal sites in patients with rheumatoid arthritis: a pilot study
BMC Oral Health volume 21, Article number: 248 (2021)
Subgingival microbiome in disease-associated subgingival sites is known to be dysbiotic and significantly altered. In patients with rheumatoid arthritis (RA), the extent of dysbiosis in disease- and health-associated subgingival sites is not clear.
8 RA and 10 non-RA subjects were recruited for this pilot study. All subjects received full oral examination and underwent collection of subgingival plaque samples from both shallow (periodontal health-associated, probing depth ≤ 3mm) and deep subgingival sites (periodontal disease-associated, probing depth ≥ 4 mm). RA subjects also had rheumatological evaluation. Plaque community profiles were analyzed using 16 S rRNA sequencing.
The phylogenetic diversity of microbial communities in both RA and non-RA controls was significantly higher in deep subgingival sites compared to shallow sites (p = 0.022), and the overall subgingival microbiome clustered primarily according to probing depth (i.e. shallow versus deep sites), and not separated by RA status. While a large number of differentially abundant taxa and gene functions was observed between deep and shallow sites as expected in non-RA controls, we found very few differentially abundant taxa and gene functions between deep and shallow sites in RA subjects. In addition, compared to non-RA controls, the UniFrac distances between deep and shallow sites in RA subjects were smaller, suggesting increased similarity between deep and shallow subgingival microbiome in RA. Streptococcus parasanguinis and Actinomyces meyeri were overabundant in RA subjects, while Gemella morbillorum, Kingella denitrificans, Prevotella melaninogenica and Leptotrichia spp. were more abundant in non-RA subjects.
The aggregate subgingival microbiome was not significantly different between individuals with and without rheumatoid arthritis. Although the differences in the overall subgingival microbiome was driven primarily by probing depth, in contrast to the substantial microbiome differences typically seen between deep and shallow sites in non-RA patients, the microbiome of deep and shallow sites in RA patients were more similar to each other. These results suggest that factors associated with RA may modulate the ecology of subgingival microbiome and its relationship to periodontal disease, the basis of which remains unknown but warrants further investigation.
The human microbiome is a complex collection of microbial communities inhabiting many sites in the human body that perform functions required for the health of a host. The interactions between our immune system and the microbiome generally mediate a symbiotic relationship, where pathogenic bacteria are reduced, and commensal or beneficial organisms are tolerated. Many studies have shown that the microbiome plays a crucial role in the training and development of our immune system . For example, intestinal epithelial cells of neonates require early exposure to bacterial endotoxin in order to establish mucosal homeostasis and subsequently, permit microbial colonization without sustaining damage . As a result, there has been an increased interest in exploring the role of microbiome in autoimmune diseases. Indeed, microbial dysbiosis has been associated with several diseases characterized by immune dysregulation, including systemic lupus erythematosus, inflammatory bowel disease, and multiple sclerosis, among others . Compared to these autoimmune conditions, the role of the human microbiome in rheumatoid arthritis (RA) has received little attention.
RA is an autoimmune disease marked by chronic joint inflammation and pain, which can progress to irreparable cartilage and bone damage . The etiology of RA remains unknown, but both genetic predispositions and environmental triggers are likely responsible . Dysbiosis of the subgingival microbiome and the resulting periodontitis have been implicated as an environmental trigger [5,6,7]. Periodontal disease, a condition often associated with changes in the subgingival microbiome, is a multifactorial disease marked by inflammation and destruction of tooth-supporting tissue . Individuals with RA have an increased incidence and severity of periodontal disease as compared to their healthy counterparts [9, 10]. Even after controlling for arthritic symptoms, patients with RA are more likely to have chronic periodontitis [9,10,11,12,13]. Furthermore, the pathophysiology of RA and chronic periodontitis shares similar features, both involving immune cell infiltration into tissue, initiation of an inflammatory environment, and subsequent tissue destruction [4, 14]. Given these shared characteristics, the relationship between rheumatoid arthritis and periodontal disease is an active area of investigation.
Several studies have examined the associations between specific subgingival microbes and RA symptoms. [9, 15,16,17,18]. Elevated levels of DNA from Porphyromonas gingivalis, a pathogen commonly associated with periodontal disease, have been found in the synovial fluid of RA patients . In vitro studies have shown that P. gingivalis is capable of infiltrating primary chondrocytes and causing cell damage [20, 21]. In addition, P. gingivalis expresses peptidylarginine deiminase and produces citrullinated human proteins, which are key molecules in the initiation and progression of RA. . Similarly, Aggregatibacter actinomycetemcomitans produces leukotoxin A, which increases the production and release of citrullinated proteins in neutrophils. Indeed, RA patients have higher rates of leukotoxin A exposure compared to their healthy counterparts . This exposure to citrullinated proteins results in the development of autoantibodies to citrullinated proteins (ACPAs) in genetically susceptible individuals . Taken together, these studies suggest that oral bacteria or their metabolites may translocate from the mouth to the joints, infiltrate chondrocytes, increase citrullinated proteins, and cause cellular damage, potentially contributing to the pathogenesis of RA.
In limiting the analysis to a few selected organisms, most prior studies did not account for the complexity and dynamics of the overall subgingival microbiome, which collectively may influence RA pathogenesis. Recent advances in sequencing technology have provided a greater resolution of the subgingival microbiome, allowing for both species- and community-level comparisons between individuals. A few studies have shown that the salivary microbiome differs between healthy and RA patients [25, 26], but data for subgingival microbial dysbiosis have been mixed [15,16,17,18]. In most studies, the subgingival microbiome of subjects with and without RA share similar diversity [15, 18]. In addition, the influence of RA disease severity on microbial diversity remains unclear as increased disease activity was both negatively  and positively  associated with different diversity metrics. Lopez-Oliva et al.  found that RA influenced the overall structure of subgingival communities, but this was not supported by other studies [15, 17, 18].
Given the similarities in pathophysiology, shared risk factors, and the increased incidence and severity of periodontal disease in RA subjects, there is significant interest in investigating the linkage between these two disease processes. Previous studies of subgingival microbiome in RA patients have shown conflicting results. Several studies did not account for the severity of periodontal disease in their sampling process, as only the deepest subgingival sites or random sites were sampled [15,16,17]. In this pilot study, we analzyed the subgingival microbiome of both shallow (health-associated) and deep (disease-associated) subgingival sites in patients with rheumatoid arthritis (RA), and then compare to non-RA household controls. We hypothesized that subgingival microbiome of diseased and healthy subgingival sites in RA patients are altered compared to non-RA controls, which may offer insights into the interplay between subgingival microbiome and RA pathogenesis.
A total of 10 subjects with rheumatoid arthritis (RA) who were naïve to biologic therapy and 10 household members of the RA patients (non-RA controls) were recruited from the Rheumatology Clinic at UF Health, Gainesville, Florida. Subjects were excluded if they were (1) younger than 18 years of age, (2) had fewer than 8 natural teeth, (3) antibiotic treatment in the three months preceding sample collection, (4) immune compromising conditions (i.e., patients with HIV/AIDS, or on immunosuppressive regimens for conditions other than rheumatoid arthritis), or (5) participation in another clinical study involving the use of dental products one week prior to sample collection or during the study period. Informed consent was obtained from all subjects. The protocol was approved by the Institutional Review Board at the University of Florida. All RA subjects were on DMARDs at the time of sampling, and 75% were in remission. No RA subjects were on glucocorticoids or had been previously treated with biologics (e.g. adalimumab, etanercept, or golimumab).
Medical and oral examinations
Each RA patient underwent extensive medical screening including evaluation of autoantibodies (RF and ACPA) and inflammatory markers (ESR and CRP). Additionally, RAPID3 score, a measure of RA severity, was determined and the treatment history of RA was collected. Two of the 10 RA patients were excluded from the final analysis, as one subject had previously received biologic therapy and the diagnosis of a second subject was subsequently felt to be more consistent with mixed connective tissue disease rather than RA. Of the 8 RA subjects included in the study, 7 were seropositive for RA. Determination of remission was based on ACR/European League Against Rheumatism (EULAR) 2011 remission criteria . Information on previous and recent dental treatments were not collected during the study visit.
All study participants underwent oral examination and collection of subgingival plaque samples in the Dental Clinical Research Unit at UF Health. Oral examinations included the assessment of dental and periodontal health status by the same dentist (MMN). Subjects were grouped according to periodontal health status as presenting periodontal health or periodontitis . Periodontal health was defined as presenting sites with both probing depth (PD) and clinical attachment level (CAL) ≤ 3 mm, and bleeding on probing (BoP) < 10 %. Periodontitis at different stages of severity (I to III) was defined as presenting: (1) more than 30 % of the tooth sites with both PD and clinical attachment level (CAL) ≥ 4 mm and BoP, and (2) minimum of 6 teeth with at least one site of PD and CAL ≥ 5 mm and BoP as described elsewhere [8, 28].
Subjects were required to refrain from tooth brushing, flossing, eating, and drinking (anything other than water) for 12 h prior to sample collection. For each subject, one subgingival plaque sample was obtained from a shallow site (probing depth ≤ 3mm) and a second sample from a deep subgingival site (periodontal disease-associated, probing depth ≥ 4 mm).
Each type of plaque sample was pooled from 1 to 6 sites on a single tooth or multiple teeth with similar health conditions using sterile periodontal paper points and transferred directly to a microcentrifuge bead-beating tube containing buffer for storage (PowerBead tubes, Qiagen, Venlo, Netherlands). Sample tubes were placed on ice, and stored at − 80 °C until processing.
Microbial profiling and microbiome analysis
Bacterial DNA was extracted using the DNeasy Powersoil Kit (Qiagen, Venlo, Netherlands) with the following modifications. An SDS-containing solution was added to the PowerBead tubes to facilitate cellular lysis and membrane breakdown. Tubes were secured to a Vortex Adapter, mixed for 10 min and heated/mixed at 70 °C for 10 min. The supernatant was transferred to a clean, 2 mL collection tube and a solution which precipitates non-DNA material was added. Samples were incubated at 4 °C for 10 min, centrifuged and the supernatant was transferred to another clean tube. This step was repeated, and then a highly concentrated salt solution was added to facilitate DNA binding to the spin column. Samples were then loaded onto spin columns and washed twice to remove contaminants. Purified DNA was eluted from the column and stored at − 20 °C until further processing.
The 16S rRNA genes (V1-V3 region) were amplified using primers 27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 534R (5′-ATTACCGCGGCTGCTGG-3′) with barcodes to allow multiplex deep sequencing as described previously . The final PCR reactions contained 0.75 U Accuprime Taq High Fidelity Polymerase (Invitrogen, Carlsbad, CA), 2 µL 10X PCR buffer II, 100 nM forward primer, 100 nM reverse primer, 2 µL purified DNA template and sterile DNA/RNA-free water in a total volume of 20 µL. Cycling conditions were as follows: denaturation at 95 °C for 2 min, 25 cycles of denaturation at 95 °C for 20 s, annealing at 56 °C for 30 s and extension at 72 °C for 5 min. PCR products for each sample were analyzed on a SYBR Safe 1 % agarose gel (Invitrogen, Carlsbad, CA) to confirm an expected size of 600 base pairs. The amplicons were excised, purified using NucleoSpin Gel and PCR Clean-up kit (Macherey-Nagel, Bethehem, PA), and then quantified using Qubit HS DNA quantification kits (Invitrogen, Carlsbad, CA). Equimolar concentrations of DNA were then pooled and purified, and qPCR was performed using the Library Quant Kit (Kapa Biosystems, Wilmington, MA). 16 S rRNA sequencing was performed on the Illumina MiSeq platform using the MiSeq Reagent kit V3 and PhiX control V3 kit (Illumina, San Diego, CA).
Raw paired-end MiSeq reads of 300 nucleotides each (covering the V1–V3 hypervariable region of the 16 S rRNA gene using primers 27 F and 534R) were processed using custom scripts in R v.3.4.2 (R Core Team, 2018). Reads were filtered based on exact matches to barcode/primer and an average quality score of 30. Samples were de-multiplexed according to the combination of their unique variable length barcodes (4 to 8 nucleotides) on each paired end. The barcodes and primers were trimmed for downstream analysis. To reconstruct the original contiguous amplicon, paired-end reads were joined using FLASh (Fast Length Adjustment of SHort reads), with a minimum overlap of 10 base pairs. USEARCH alignment was employed with a 97% identity and 80% aligned query threshold to assign OTUs with reference taxonomic information from the Human Oral Microbiome Database (16 S rRNA RefSeq Version 10.1) [30, 31] to each joined read. Reads that did not meet the filtering criteria were excluded from subsequent analysis. Alpha and beta diversity metrics were generated in QIIME2 (version 2018.8, available at https://qiime2.org/)  using the core-metrics-phylogenetic pipeline and later analyzed in R v.3.4.2 (R Core Team, 2018). Reads for each sample were rarefied to 7500 prior to alpha and beta diversity analysis.
Age, ESR, RAPID3 and probing depth of periodontal sites (deep vs. shallow) were compared between RA and non-RA controls using Student’s t-tests. The other demographic and clinical were compared using Fisher’s Exact test. All statistical analyses were conducted in R v.3.4.2 (R Core Team, 2018) unless otherwise stated.
For microbiome analysis, mixed linear models were conducted using lmer() in the lme4 package (v.1.1–19)  to determine how RA status, probing depth, and their interaction influenced each alpha diversity metric. Subject identity was included as a random effect to account for the non-independence of the deep and shallow samples coming from the same subject. Within RA patients, mixed linear models were used to determine how alpha diversity varied with RA disease activity, measured using the RAPID3 score. Again, subject identity was included as a random effect. To compare the overall community composition, permutational multivariate analysis of variance (PERMANOVA) was performed on weighted and unweighted UniFrac distances between samples using adonis() in vegan v.2.5-3 . UniFrac distances are distance metrics that incorporate phylogenetic relatedness and account for the presence (unweighted) or abundance (weighted) of OTUs. Clustering was visualized with Principal Coordinate Analysis using betadisper() in vegan v.2.5-3 . Differential abundance analysis was performed using LEfSe  to identify OTUs associated with RA status or probing depth. Differentially abundant features met the minimum LDA threshold of 2. To compare how the community differences between shallow and deep sites of a single individual differed across RA status, intra-subject distances (the UniFrac distance between shallow and deep site of the same subject) were tested using unpaired t-tests, as subjects could not be paired with their household match because two RA subjects were excluded and not all RA patients had both a shallow and a deep site. Metagenomic functions were predicted using PICRUSt  and subjected to biomarker analysis using LEfSe. Differentially abundant gene functions met a minimum LDA threshold of 2. The bioBakery Python tools package contained the Python scripts used in LEfSe and PICRUSt analyses.
All methods were carried out in accordance relevant guidelines and regulations.
Clinical characteristics and demographics of RA and non-RA controls are shown in Table 1. No significant differences in age, gender, race, number of caries, and periodontal health status were observed between the two groups.
A total of 34 samples were collected, including both shallow and deep sites from each of the 10 non-RA controls and 6 RA subjects. For the remaining 2 RA subjects, a deep site was not collected from one (G1-10) and a second had no deep sites to sample (G1-2). 16 S rRNA Illumina sequencing of the 34 samples generated 970,260 raw reads. After trimming and quality control, 882,526 reads were available for analysis with an average of 25,957 reads per sample (range: 8246−50,936 reads). Reads clustered into 499 distinct operational taxonomic units (OTUs), belonging to 10 different phyla and 97 different genera. The most abundant phyla were Firmicutes (29.9%), Bacteroidetes (27.1%), Fusobacteria (21.3%) and Proteobacteria (13.2%). The dominant genera included Fusobacterium (17.9%), Prevotella (14.6%), Veillonella (11.0%), Streptococcus (10.1%), and Porphyromonas (7.1%).
To characterize the subgingival microbiome in RA, we first calculated the alpha diversity using observed OTUs, Faith’s phylogenetic diversity index, and Shannon’s diversity index (Fig. 1). Analysis using mixed linear models showed that microbial communities in deep subgingival sites had higher phylogenetic diversity than shallow sites (b = 1.93 ± 0.73SE, p = 0.022; Table 2). However, neither RA status (b = 0.60 ± 1.10SE, p = 0.69) nor the interaction between RA and probing depth (i.e. deep versus shallow sites) was associated with phylogenetic diversity (b = − 1.08 ± 1.20SE, p = 0.38). Furthermore, neither RA status nor probing depth was associated with observed OTUs or Shannon diversity. Mixed linear models showed that RA disease activity, as measured by the RAPID3 score, was not associated with differences in any of the alpha diversity measurements for both shallow and deep sites (see Additional file 2: Table S1).
Next, we compared differences at the microbial community level by calculating weighted and unweighted UniFrac distances between samples (Fig. 2). For both weighted (Fig. 2a, b) and unweighted (Fig. 2c, d) UniFrac analysis, samples clustered primarily by probing depth along the first principal coordinate axis (PERMANOVA, weighted UniFrac: p < 0.006; unweighted UniFrac: p < 0.009). In contrast, samples did not separate according to RA status (PERMANOVA, weighted UniFrac: p = 0.826; unweighted UniFrac: p = 0.881). Additionally, we found no interaction between probing depth and RA status (PERMANOVA, weighted UniFrac: p = 0.861; unweighted UniFrac: p = 0.997).
To further compare microbial composition, we utilized Linear discriminant analysis Effect Size (LEfSe) to identify specific OTUs that were differentially abundant between groups. A total of 45 OTUs were differentially enriched between deep and shallow sites in non-RA controls (see Additional file 1: Fig. S1). OTUs abundant in deep subgingival sites of non-RA subjects included several members of the red and orange complexes , and were associated with 33 differentially enriched metagenome functions predicted by PICRUSt (see Additional file 1: Fig. S2). Surprisingly, only 11 OTUs (see Fig. 3) and one predicted metagenome function (see Fig. 4) in RA subjects were differentially abundant between deep and shallow sites. Thus, in contrast to the substantial differences in subgingival microbiome between shallow and deep sites in non-RA controls, these results suggest increased similarity in the microbiome of shallow and deep periodontal sites in RA patients.
To quantify these differences, we calculated the weighted and unweighted Unifrac distances between paired sites (deep versus shallow) within each subject, then compared between RA and non-RA controls. Using this approach, large intra-subject distances indicate a high degree of separation between shallow and deep site microbiomes, whereas small distances suggest more similiar microbiomes. The weighted and unweighted intra-subject UniFrac distances for RA patients were 31 and 24% lower than non-RA controls, respectively (Fig. 5), although these differences did not achieve statistical significance (p = 0.093 and p = 0.085, respectively, unpaired t-tests). Thus, these results support a trend toward increased similarity in the subgingival microbiome between shallow and deep sites in RA.
Finally, we examined whether there were specific OTUs in the subgingival microbiome that were unique to RA subjects. Differentially abundant OTUs between RA and non-RA controls were identified using Linear discriminant analysis Effect Size (LEfSe). For both shallow and deep sites, Actinomyces meyeri and Streptococcus parasanguinis were more abundant in RA subjects compared to non-RA controls, whereas Gemella morbillorum, Kingella denitrificans, Prevotella melaninogenica, and a Leptotrichia species were depleted in RA subjects (Fig. 6). G. morbillorum was the most abundant OTU among the 6 differentially enriched OTUs but only constituted 0.31% of the reads in non-RA controls and 0.26% in RA, indicating that the four differentially abundant OTUs associated with RA were all minority components of the subgingival microbiome. There were no differences in the abundance of predicted gene functions between RA and non-RA controls.
In this pilot study, the subgingival microbiome of shallow and deep periodontal sites in patients with rheumatoid arthritis (RA) was determined by 16 S rRNA sequencing and compared to those of non-RA household controls. Historically, methods such as organism-specific PCR/qPCR and anaerobic culture techniques have been used to examine subgingival communities in RA patients [38,39,40]. However, these methods could not compare the overall microbial community profile. Only in the past decade have rapid advances in deep-sequencing technology allowed for an unbiased and more complete picture of the oral microbiome. Prior studies have used Illumina 16 s rRNA sequencing, Roche/454 16 S rRNA sequencing, or whole-genome shotgun sequencing to characterize the oral microbiome in RA [15,16,17,18, 25, 26]. Rather than sampling only the deepest subgingival sites as previously reported , we sampled paired shallow and deep sites of each subject to define the relationship between healthy and diseased subgingival microbiome within individuals in the context of RA.
It has been reported that periodontal disease is more severe and more prevalent in patients with rheumatoid arthritis . Thus it was not surprising to observe here greater pocket depths for both clinically healthy and diseased sites in RA than the controls. Consistent with previous work , we found an increase in phylogenetic diversity for deep sites compared to shallow sites, suggesting that the dysbiosis seen in periodontal disease may be characterized by increased microbial diversity. This is a very interesting finding because it differs from a general trend observed in other body habitats where dysbiosis and disease progression often involve a loss of microbial diversity [42, 43].
Consistent with previous studies [15, 17, 18] with the exception of Ganesan et al.  and Dabdoub et al. , we observed clustering of subgingival microbiomes by probing depth in both weighted and unweighted UniFrac analysis. However, we did not observe clustering according to RA status. This is consistent with some studies [15, 18] but is in contrast with others [16, 26]. The reasons for these discrepancies are unclear but may be related to the chronicity and/or remission status of RA patients, severity of periodontal disease, use of controls within the family unit, the number of sampling sites, and/or sequencing methods.
Ganesan et al. and Dabdoub et al. found large differences in subgingival communities between subjects with periodontal disease and those who are periodontally healthy [44, 45]. However, their work did not find differences between shallow and deep sites within the same subjects, suggesting that microbial dysbiosis is a global event rather than a local shift [44, 45]. Here, we found that even in the same individual, the communities of shallow and deep sites are different, implying that dysbiosis is a local event. Interestingly, the analysis of intra-subject UniFrac distances and LEfSe suggest that the differences in subgingival microbiome between deep and shallow sites were smaller in RA patients compared to non-RA controls. One possible explanation for the increased similarity in microbiome between shallow and deep sites of RA patients is that subgingival dysbiosis may have already begun in shallow sites, thereby reducing the microbial community differences between shallow and deep sites in RA. Given the small sample size of our pilot study, further investigations using longitudinal sampling and a larger study cohort are needed to confirm these findings.
The analysis of microbiota composition revealed that two OTUs, S. parasanguinis (OTU 721) and A. meyeri (OTU 671), were more abundant in RA as compared to controls. While neither OTU has been associated with RA to date, several Actinomyces species have been shown to be enriched  and several Streptococcus species were depleted in RA patients [16, 18]. S. parasanguinis is an oral commensal capable of bacterial antagonism, inhibiting the growth of periodontal pathogens including Aggregatibacter actinomycetemcomitans, P. gingivalis and Prevotella intermedia . However, one study showed that S. parasanguinis contributed to P. gingivalis colonization . Nonetheless, this organism has been associated with refractory periodontitis  and alveolar bone loss , and its role in the colonization of periodontal pathogens may be variable. A. meyeri is an opportunistic pathogen and has been shown to cause empyema necessitans . The role of A. meyeri in periodontal disease has not been well studied, though this organism has been associated with apical periodontitis . P. gingivalis possesses mechanisms for citrullinated proteins, a significant factor in the initiation and propagation of RA . However, we did not find an association with RA in our study.
We identified 4 OTUs underrepresented in RA subjects: G. morbillorum (OTU 046), K. denitrificans (OTU 582), P. melaninogenica (OTU 469) and Leptotrichia spp. (OTU 223). G. morbillorum has been associated with oral infections, including endodontic disease , apical periodontitis , and chronic periodontitis [54, 55]. This organism has not been directly associated with RA previously, but it has been reported to cause septic arthritis , including in RA patients . Similarly to G. morbillorum, K. denitrificans is an oral and respiratory tract commensal which can cause infections such as endocarditis and peritonitis in susceptible individuals [58, 59]. This OTU was also depleted in salivary samples of RA patients compared to controls . Contrary to our findings, enrichment of both Prevotella and Leptotrichia species has been shown in RA [16, 18]. Increase abundance of P. melaninogenica in the subgingival microbiome has been associated with a number of different conditions, including inflammatory bowel disease , obesity , and rheumatoid arthritis and osteoarthritis patients . Furthermore, elevated levels of antibodies to P. melaninogenica have been reported for RA patients , thus implicating a role in RA . Given the differences between our data and previous work, the role of P. melaninogenica and Leptotrichia species in RA pathogenesis warrants further investigation.
There are several limitations in this study. First, although the probing depths for deep sites were comparable between RA and non-RA controls, the shallow sites (probing depth ≤ 3 mm) were deeper in RA patients than non-RA controls. These differences in probing depths, however, did not contribute to differences in the microbiome structure of shallow sites between RA and non-RA controls. Second, we did not collect information on smoking history or detailed medical histories, both of which could have influenced the composition of subgingival microbiome. Third, while none of the RA patients received biologic therapy, they were not all receiving the same drug regimens for RA which could have contributed to variations in the microbiome. Lastly, given that this was a pilot study, our sample size was small. Nonetheless, our data suggest that RA may modulate the subgingival microbiome and support further investigation in a larger cohort study to confirm these results.
To our knowledge, this is the first study that compares paired deep and shallow sites within RA patients. Our findings suggest that the aggregate subgingival microbiome is not significantly different between individuals with and without rheumatoid arthritis and that the probing depth primarily drives the differences in the overall profile of the subgingival microbiome. Compared to the substantial differences in subgingival microbiome typically seen between deep and shallow sites in the general non-RA population , our data suggest that these communities are more similar to each other in RA patients. These results raise the hypothesis that factors associated with the onset of RA may modulate the ecology of subgingival microbiome and its relationship to periodontal disease, the basis of which remains unknown but warrants further investigation. Additionally, we identified several OTUs, which were differentially abundant between RA patients and non-RA controls, and their role in RA and periodontal disease should be explored in future studies.
Autoantibodies to citrullinated proteins
Clinical attachment level
Bleeding upon probing
Belkaid Y, Hand TW. Role of the microbiota in immunity and inflammation. Cell. 2014 ;157(1):121–41.
Chassin C, Kocur M, Pott J, Duerr CU, Gütle D, Lotz M, et al. miR-146a mediates protective innate immune tolerance in the neonate intestine. Cell Host Microbe. 2010;8(4):358–68.
De Luca F, Shoenfeld Y. The microbiome in autoimmune diseases. Clin Exp Immunol. 2019;195(1):74–85.
Smolen JS, Aletaha D, McInnes IB. Rheumatoid arthritis. Lancet. 2016;388(10055):2023–38.
Jalil SF, Arshad M, Bhatti A, Ahmad J, Akbar F, Ali S, et al. Rheumatoid arthritis: what have we learned about the causing factors? Pak J Pharm Sci. 2016;29(2):629–45.
Deane KD, Demoruelle MK, Kelmenson LB, Kuhn KA, Norris JM, Holers VM. Genetic and environmental risk factors for rheumatoid arthritis. Best Pract Res Clin Rheumatol. 2017;31(1):3–18.
Alam J, Jantan I, Bukhari SNA. Rheumatoid arthritis: recent advances on its etiology, role of cytokines and pharmacotherapy. Biomed Pharmacother. 2017;92:615–33.
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 Clin Periodontol. 2018;45:S162-70.
Dissick A, Redman RS, Jones M, Rangan BV, Reimold A, Griffiths GR, et al. Association of periodontitis with rheumatoid arthritis: a pilot study. J Periodontol. 2010;81(2):223–30.
Mikuls TR, Payne JB, Yu F, Thiele GM, Reynolds RJ, Cannon GW, et al. Periodontitis and Porphyromonas gingivalis in patients with rheumatoid arthritis. Arthritis Rheumatol. 2014;66(5):1090–100.
Tang Q, Fu H, Qin B, Hu Z, Liu Y, Liang Y, et al. A possible link between rheumatoid arthritis and periodontitis: a systematic review and meta-analysis. Int J Periodontics Restorative Dent. 2017;37(1):79–86.
Cheng Z, Meade J, Devine DA. Periodontal disease and periodontal bacteria as triggers for rheumatoid arthritis. Best Pract Res Clin Rheumatol. 2017;31(1):19–30.
Bingham C, Moni M. Periodontal disease and rheumatoid arthritis: the evidence accumulates for complex pathobiologic interactions. Curr Opin Rheumatol. 2015;25(3):345–53.
Potempa J, Mydel P, Koziel J. The case for periodontitis in the pathogenesis of rheumatoid arthritis. Nat Rev Rheumatol. 2017;13(10):606–20.
Mikuls TR, Walker C, Qiu F, Yu F, Thiele GM, Alfant B, et al. The subgingival microbiome in patients with established rheumatoid arthritis. Rheumatology (UK). 2018;57(7):1162–72.
Lopez-Oliva I, Paropkari AD, Saraswat S, Serban S, Yonel Z, Sharma P, et al. Dysbiotic subgingival microbial communities in periodontally healthy patients with rheumatoid arthritis. Arthritis Rheumatol. 2018;70(7):1008–13.
Beyer K, Zaura E, Brandt BW, Buijs MJ, Brun JG, Crielaard W, et al. Subgingival microbiome of rheumatoid arthritis patients in relation to their disease status and periodontal health. PLoS One. 2018;13(9):e0202278.
Scher JU, Ubeda C, Equinda M, Khanin R, Buischi Y, Viale A, et al. Periodontal disease and the oral microbiota in new-onset rheumatoid arthritis. Arthritis Rheum. 2012;64(10):3083–94.
Reichert S, Haffner M, Keyßer G, Schäfer C, Stein JM, Schaller HG, et al. Detection of oral bacterial DNA in synovial fluid. J Clin Periodontol. 2013;40(6):591–8.
Pischon N, Röhner E, Hocke A, N’Guessan P, Müller HC, Matziolis G, et al. Effects of Porphyromonas gingivalis on cell cycle progression and apoptosis of primary human chondrocytes. Ann Rheum Dis. 2009;68(12):1902–7.
Röhner E, Detert J, Kolar P, Hocke A, N’Guessan P, Matziolis G, et al. Induced apoptosis of chondrocytes by Porphyromonas gingivalis as a possible pathway for cartilage loss in rheumatoid arthritis. Calcif Tissue Int. 2010;87(4):333–40.
Wegner N, Wait R, Sroka A, Eick S, Nguyen KA, Lundberg K, et al. Peptidylarginine deiminase from Porphyromonas gingivalis citrullinates human fibrinogen and α-enolase: implications for autoimmunity in rheumatoid arthritis. Arthritis Rheum. 2010;62(9):2662–72.
Konig MF, Abusleme L, Reinholdt J, Palmer RJ, Teles RP, Sampson K, et al. Aggregatibacter actinomycetemcomitans-induced hypercitrullination links periodontal infection to autoimmunity in rheumatoid arthritis. Sci Transl Med. 2016;8(369):369ra176.
Firestein GS, McInnes IB. Immunopathogenesis of rheumatoid arthritis. Immunity. 2017;46(2):183–96.
Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, et al. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat Med. 2015;21(8):895–905.
Chen B, Zhao Y, Li S, Yang L, Wang H, Wang T, et al. Variations in oral microbiome profiles in rheumatoid arthritis and osteoarthritis with potential biomarkers for arthritis screening. Sci Rep. 2018;8(1):17126.
Felson DT, Smolen JS, Wells G, Zhang B, van Tuyl LHD, Funovits J, et al. American College of Rheumatology/European league against rheumatism provisional definition of remission in rheumatoid arthritis for clinical trials. Ann Rheum Dis. 2011;70(3):404–13.
Chapple ILC, Mealey BL, Van Dyke TE, Bartold PM, Dommisch H, Eickholz P, et al. Periodontal health and gingival diseases and conditions on an intact and a reduced periodontium: consensus report of workgroup 1 of the 2017 World Workshop on the classification of periodontal and peri-implant diseases and conditions. J Clin Periodontol. 2018;89:S68–77.
Kozlov A, Bean L, Hill E V., Zhao L, Li E, Wang GP. Molecular identification of bacteria in intra-abdominal abscesses using deep sequencing. Open Forum Infect Dis. 2018;5(2):ofy025.
Escapa IF, Chen T, Huang Y, Gajare P, Dewhirst FE, Lemon KP. 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. 2018;3(6):e00187-18.
Chen T, Yu W-H, Izard J, Baranova OV, Lakshmanan A, Dewhirst FE. The human oral microbiome database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Database. 2010;2010:baq013–baq013.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
Bates D, Mächler M, Bolker BM, Walker SC. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlin D et al. vegan: community ecology package. 2018.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.
Socransky SS, Haffajee AD, Cugini MA, Smith C, Kent RL. Microbial complexes in subgingival plaque. J Clin Periodontol. 1998;25(2):134–44.
Smit MD, Westra J, Vissink A, Doornbos-van der Meer B, Brouwer E, van Winkelhoff AJ. Periodontitis in established rheumatoid arthritis patients: a cross-sectional clinical, microbiological and serological study. Arthritis Res Ther. 2012;14(5):R222.
Ziebolz D, Pabel SO, Lange K, Krohn-Grimberghe B, Hornecker E, Mausberg RF. Clinical periodontal and microbiologic parameters in patients with rheumatoid arthritis. J Periodontol. 2011;82(10):1424–32.
Wolff B, Berger T, Frese C, Max R, Blank N, Lorenz H-M, et al. Oral status in patients with early rheumatoid arthritis: a prospective, case-control study. Rheumatology. 2014;53(3):526–31.
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.
Scher JU, Ubeda C, Artacho A, Attur M, Isaac S, Reddy SM, et al. Decreased bacterial diversity characterizes the altered gut microbiota in patients with psoriatic arthritis, resembling dysbiosis in inflammatory bowel disease. Arthritis Rheumatol. 2015;67(1):128–39.
Kriss M, Hazleton KZ, Nusbacher NM, Martin CG, Lozupone CA. Low diversity gut microbiota dysbiosis: drivers, functional implications and recovery. Curr Opin Microbiol. 2018;44:34–40.
Ganesan SM, Joshi V, Fellows M, Dabdoub SM, Nagaraja HN, O’Donnell B, et al. A tale of two risks: smoking, diabetes and the subgingival microbiome. ISME J. 2017;11(9):2075–89.
Dabdoub SM, Ganesan SM, Kumar PS. Comparative metagenomics reveals taxonomically idiosyncratic yet functionally congruent communities in periodontitis. Sci Rep. 2016;6(1):38993.
Herrero ER, Slomka V, Bernaerts K, Boon N, Hernandez-Sanabria E, Passoni BB, et al. Antimicrobial effects of commensal oral species are regulated by environmental factors. J Dent. 2016;47:23–33.
Maeda K, Nagata H, Nonaka A, Kataoka K, Tanaka M. Oral streptococcal glyceraldehyde-3-phosphate dehydrogenase mediates interaction with fimbriae. Microbes Infect. 2004;6(13):1163–70.
Colombo AP V, Boches SK, Cotton SL, Goodson JM, Kent R, Haffajee AD, et al. Comparisons of subgingival microbial profiles of refractory periodontitis, severe periodontitis, and periodontal health using the human oral microbe identification microarray. J Periodontol. 2009;80(9):1421–32.
Fine DH, Markowitz K, Fairlie K, Tischio-Bereski D, Ferrendiz J, Furgang D, et al. A consortium of Aggregatibacter actinomycetemcomitans, Streptococcus parasanguinis, and Filifactor alocis is present in sites prior to bone loss in a longitudinal study of localized aggressive periodontitis. J Clin Microbiol. 2013;51(9):2850–61.
Ellebrecht DB, Pross MMF, Schierholz S, Palade E. Actinomyces meyeri empyema necessitatis—a case report and review of the literature. Surg J. 2019;5:e57–61
Johnson BT, Mayo JA, Jeansonne BG. beta-Hemolytic streptococci and other beta-hemolytic organisms in apical periodontitis and severe marginal periodontitis. Endod Dent Traumatol. 1999;15(3):102–8.
Gomes BPFA, Montagner F, Jacinto RC, Pinheiro ET, Zaia AA, Ferraz CCR, et al. Gemella morbillorum in primary and secondary/persistent endodontic infections. Oral Surgery Oral Med Oral Pathol Oral Radiol Endodontology. 2008;105(4):519–25.
Signoretti FGC, Gomes BPFA, Montagner F, Jacinto RC. Investigation of cultivable bacteria isolated from longstanding retreatment-resistant lesions of teeth with apical periodontitis. J Endod. 2013;39(10):1240–4.
Albandar JM, Khattab R, Monem F, Barbuto SM, Paster BJ. The subgingival microbiota of Papillon–Lefèvre syndrome. J Periodontol. 2012;83(7):902–8.
Kuriyama T, Karasawa T, Nakagawa K, Yamamoto E, Nakamura S. Bacteriology and antimicrobial susceptibility of gram-positive cocci isolated from pus specimens of orofacial odontogenic infections. Oral Microbiol Immunol. 2002;17(2):132–5.
Desmottes M-C, Brehier Q, Bertolini E, Monteiro I, Terreaux W. Septic arthritis of the knee due to Gemella morbillorum. Int J Rheum Dis. 2018;21(5):1146–7.
Tang F-G, Yuan F-L, Liu J-K, Long Q. Suppurative arthritis caused by Gemella morbillorum in a patient with rheumatoid arthritis of the knee. Nan Fang Yi Ke Da Xue Xue Bao. 2012;32(3):437–8.
Hassan IJ, Hayek L. Endocarditis caused by Kingella denitrificans. J Infect. 1993;27(3):291–5.
Kopyt N, Kumar A, Agrawal V. A case of suppurative peritonitis by a commensal oral organism, Kingella denitrificans, in an adult peritoneal dialysis patient. Perit Dial Int. 2015;35(1):105–7.
Brito F, Zaltman C, Carvalho ATP, Fischer RG, Persson R, Gustafsson A, et al. Subgingival microflora in inflammatory bowel disease patients with untreated periodontitis. Eur J Gastroenterol Hepatol. 2013;25(2):239–45.
Maciel SS, Feres M, Gonçalves TED, Zimmermann GS, da Silva HDP, Figueiredo LC, et al. Does obesity influence the subgingival microbiota composition in periodontal health and disease? J Clin Periodontol. 2016;43(12):1003–12.
Ogrendik M, Kokino S, Ozdemir F, Bird PS, Hamlet S. Serum antibodies to oral anaerobic bacteria in patients with rheumatoid arthritis. MedGenMed Medscape Gen Med. 2005;7(2):2.
Ogrendik M. Does periodontopathic bacterial infection contribute to the etiopathogenesis of the autoimmune disease rheumatoid arthritis? Discov Med. 2012;13(72):349–55.
Kirst ME, Li EC, Alfant B, Chi YY, Walker C, Magnusson I, et al. Dysbiosis and alterations in predicted functions of the subgingival microbiome in chronic periodontitis. Appl Environ Microbiol. 2015;81(2):783–93.
We thank members of the Wang laboratory for helpful suggestions.
This work was supported by Department of Defense grant W81XWH-15-1-0320, and in part by the Gatorade Trust through funds distributed by the University of Florida, Department of Medicine.
Ethics approval and consent to participate
Informed consent was obtained from all patients and controls under a protocol approved by the Institutional Review Board at the University of Florida.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Differentially abundant OTUs between deep and shallow sites in non-RA controls. Differentially abundant taxa were identified using LEfSe and met the minimum LDA score of 2. OTUs enriched in deep subgingival sites are shown in red, and OTUs more abundant in shallow sites are shown in dark blue. OTU: Operational taxonomic unit. Figure S2: Differentially abundant metagenome functions between deep and shallow sites in non-RA controls. Metagenome functions were predicted using 16S rRNA data and PICRUSt. Differentially abundant gene functions were identified using LEfSe and met the minimum LDA score of 2. Gene functions enriched in deep sites are shown in red, whereas functions more abundant in shallow sites are shown in dark blue.
Table S1. Association between RA disease severity and alpha diversity of the subgingival microbiome. Mixed linear models that included the RAPID3 score of RA disease severity and probe depth (shallow vs. deep) were created for each alpha diversity metric. Subject identity was included as a random effect. Model coefficient estimates, standard error (SE) of the estimates, and p values are shown. OTUs: Operational taxonomic units.
About this article
Cite this article
Lehenaff, R., Tamashiro, R., Nascimento, M.M. et al. Subgingival microbiome of deep and shallow periodontal sites in patients with rheumatoid arthritis: a pilot study. BMC Oral Health 21, 248 (2021). https://doi.org/10.1186/s12903-021-01597-x
- Subgingival microbiome
- Rheumatoid arthritis
- Microbial dysbiosis
- Periodontal disease
- 16S rRNA sequencing