Subgingival Microbiome and Specialized Pro-Resolving Lipid Mediator Pathway Profiles Are Correlated in Periodontal Inflammation

Failure of resolution pathways in periodontitis is reflected in levels of specialized pro-resolving lipid mediators (SPMs) and SPM pathway markers but their relationship with the subgingival microbiome is unclear. This study aimed to analyze and integrate lipid mediator level, SPM receptor gene expression and subgingival microbiome data in subjects with periodontitis vs. healthy controls. The study included 13 periodontally healthy and 15 periodontitis subjects that were evaluated prior to or after non-surgical periodontal therapy. Samples of gingival tissue and subgingival plaque were collected prior to and 8 weeks after non-surgical treatment; only once in the healthy group. Metabololipidomic analysis was performed to measure levels of SPMs and other relevant lipid mediators in gingiva. qRT-PCR assessed relative gene expression (2^-ΔΔCT) of known SPM receptors. 16S rRNA sequencing evaluated the relative abundance of bacterial species in subgingival plaque. Correlations between lipid mediator levels, receptor gene expression and bacterial abundance were analyzed using the Data Integration Analysis for Biomarker discovery using Latent cOmponents (DIABLO) and Sparse Partial Least Squares (SPLS) methods. Profiles of lipid mediators, receptor genes and the subgingival microbiome were distinct in the three groups. The strongest correlation existed between lipid mediator profile and subgingival microbiome profile. Multiple lipid mediators and bacterial species were highly correlated (correlation coefficient ≥0.6) in different periodontal conditions. Comparing individual correlated lipid mediators and bacterial species in periodontitis before treatment to healthy controls revealed that one bacterial species, Corynebacterium durum, and five lipid mediators, 5(S)6(R)-DiHETE, 15(S)-HEPE, 7-HDHA, 13-HDHA and 14-HDHA, were identified in both conditions. Comparing individual correlated lipid mediators and bacterial species in periodontitis before treatment to after treatment revealed that one bacterial species, Anaeroglobus geminatus, and four lipid mediators, 5(S)12(S)-DiHETE, RvD1, Maresin 1 and LTB4, were identified in both conditions. Four Selenomonas species were highly correlated with RvD1, RvE3, 5(S)12(S)-DiHETE and proinflammatory mediators in the periodontitis after treatment group. Profiles of lipid mediators, receptor gene and subgingival microbiome are associated with periodontal inflammation and correlated with each other, suggesting inflammation mediated by lipid mediators influences microbial composition in periodontitis. The role of correlated individual lipid mediators and bacterial species in periodontal inflammation have to be further studied.

Introduction

In the United States, around 60 million adults over 30 years of age have periodontitis (42.2%) with 7.8% having severe periodontitis (1). Periodontitis is a biofilm induced chronic inflammatory disease characterized by gingival inflammation and destruction of alveolar bone. Disproportionate host responses and microbiota dysbiosis are the two major etiological factors (2). Periodontal tissue damage is primarily mediated by bacterially induced immune responses. The composition of the subgingival microbiota is associated with the inflammatory status in periodontal tissue (3). It is hypothesized that the change of local environment induced by inflammation results in the shifts of the subgingival microbiota (4, 5). Many studies have shown that the composition of the subgingival microbiota shifts following periodontal therapy performed to control inflammation (6, 7). In periodontitis, the inflamed tissues with deep periodontal pocket provide an anaerobic environment with breakdown products of tissue destruction, plasma proteins and hemoglobin as nutrients for the growth of several anaerobic gram-negative bacteria, resulting in microbiota shifts. Once inflammation is controlled, the environment is not suitable for periodontal pathogens anymore and commensal microbiota can be re-established in a homeostatic relationship with the host.

Resolution of inflammation is a proactive process induced by specialized pro-resolving lipid mediators (SPMs), including lipoxins, resolvins, protectins and maresins, that bind to specific G protein-coupled receptors on a variety of cells. In the resolution phase of inflammation, there is decreased infiltration of neutrophils, reduced levels of pro-inflammatory cytokines and lipid mediators, and increased recruitment of resolving macrophages that clear the lesion by efferocytosis (8, 9). It has been demonstrated that SPMs control inflammatory diseases, such as inflammatory bowel disease (10), diabetes (11) and periodontitis (12) in the preclinical models. Specifically, in experimental periodontitis, SPMs, such as resolvin E1, topically applied on gingiva can prevent bone loss, regenerate the lost bone, change gene expression patterns in gingiva and result in shifts of the oral microbiota (13, 14). Resolution of inflammation induced by SPMs can influence the composition of the subgingival microbiota in periodontal inflammation.

In humans, SPMs have been found in milk (15), serum, lymphoid tissue (16), saliva and gingival crevicular fluid (17, 18). The levels of SPMs and related lipid mediators in various specimens are associated with the inflammatory status of mammary glands (15), the stability of atherosclerotic plaques (19), the severity of tuberculous meningitis (20) and the disease status of periodontitis (17, 18). Recently, SPMs, SPM pathway markers and SPM corresponding receptor genes are identified in gingival tissues, the periodontal inflammation site (21). Profiles of these lipid mediators and receptor genes are associated with the severity of inflammation. However, the associations between these lipid mediators and the oral microbiome have not been explored. Since preclinical studies demonstrated that resolution of inflammation induced by SPMs is associated with shifts in the taxonomic composition of the oral microbiota (13, 14), there is a need to investigate the clinical relationship between profiles of SPM relevant lipid mediators and subgingival microbiome in periodontal inflammation to clarify how inflammatory reactions mediated by SPMs affect subgingival microbial composition in humans and vice versa.

The aim of this study was to assess correlations between levels of lipid mediators including SPMs and SPM pathway markers, expression of SPM receptor genes in human gingiva and the relative abundance of bacterial species in subgingival plaque. The results of this study indicate the potential interactions between lipid mediators and bacterial species in different periodontal inflammatory conditions.

Materials and Methods

Clinical Study Design

The study was conducted in accordance with the guidelines of the World Medical Association’s Declaration of Helsinki and approved by the University of Texas Health Science Center at Houston (UTHealth) Committee for the Protection of Human Subjects (HSC-DB-16-0167). All participants provided written informed consent. The inclusion criteria were subjects aged ≥18 with ≥24 teeth and no history of systematic periodontal therapy within the past two years. These subjects should not have received anti-inflammatory drugs for more than one week and systemic antibiotics within three months before sample collection, did not routinely take fish oil supplements, had no presence of diabetes mellitus or any systemic condition that entails a diagnosis of “systemic disorders that have a major impact on the loss of periodontal tissues by influencing periodontal inflammation” (22), were not pregnant, and were not current users of tobacco products or nicotine replacement medication.

A comprehensive periodontal charting and examination were performed for all subjects. Probing depth (PD), level of free gingival margin (FGM), clinical attachment level (CAL) and presence of bleeding of probing (BOP) were measured at six sites per tooth. Full mouth series (FMS) radiographs were obtained prior to all treatments for all subjects to confirm alveolar bone level. The subjects in the healthy group were required to meet all inclusion criteria, as well as have all teeth with a probing depth of ≤3mm, clinical attachment loss of ≤2mm (except teeth with mid-buccal or lingual gingival recession), and radiographic bone levels ≤2mm from the cementoenamel junction (CEJ) (no radiographic bone loss due to periodontitis). The subjects in the periodontitis group were required to present with ≥8 teeth with a probing depth of ≥5mm, clinical attachment loss of ≥3mm, and radiographic bone levels >2mm from the CEJ.

Clinical Sample Collection

Subjects in the healthy group were seen for one clinical research visit, prior to prophylaxis. During the visit, four gingival biopsy samples were obtained from interproximal sites (mesiobuccal, distobuccal, mesiolingual, or distolingual sites of each tooth) of two representative posterior teeth and subgingival plaque samples were obtained from eight interproximal sites of the same two teeth. Subjects in the periodontitis group were seen for two clinical visits; the first visit being prior to the first scaling and root planing (SRP) appointment, and the second visit eight weeks following completion of the second SRP appointment. In these two visits, gingival biopsies and subgingival plaque samples were collected from interproximal sites of the same two representative posterior teeth. The most severely affected teeth with more bone loss than other teeth were selected for gingival biopsy and plaque sample collection. All first-visit samples collected in the periodontitis subjects were from sites with deep probing depths (≥ 5mm).

Gingival biopsies were performed as previously described (21). At each site for gingival sample collection, intrasulcular incisions were made from the papilla zenith toward the base of the interproximal papillae, but not to exceed the nearest line angle of each tooth. A subsequent horizontal incision was made at the base of the interproximal papillae to connect the intrasulcular incisions. The incisions extended to the alveolar bone, and split thickness papillae were partially elevated to permit the harvesting of the interproximal gingival tissue. The gingival sample was required to include epithelium and the underlying supracrestal connective tissue. Of the four gingival tissue samples collected, the sample selected for Lipid Mediator-SPM metabololipidomic analysis was from the site with the deepest probing depth and the sample selected for real-time quantitative reverse transcription assay (qRT-PCR) was from the site with the second-deepest probing depth. The other two gingival tissue samples were stored for future use.

Prior to subgingival plaque collection, supragingival plaque was removed with a separate sterilized scaler. Each subgingival plaque sample was collected using a sterile end of a Gracey mini-curette (Hu-Friedy, Chicago, IL, USA), which was inserted into the depth of the sulcus, with a standardized 20 single strokes. The plaque sample collected from one single site was placed in a microcentrifuge tube with 150µl of Tris-EDTA (TE) buffer, using an up and down motion to ensure deposition of the plaque sample into the TE buffer. The tube was immediately placed in liquid nitrogen and then transferred and stored at -80°C. Only one of the eight plaque samples was processed for 16S rRNA sequencing where possible to avoid affecting microbiome data by pooling samples from multiple sites with different periodontal conditions. Eight plaque samples were collected because sometimes the DNA concentration of a single site could be low. For one subject, the DNA concentration of the single site was too low for PCR amplification, thus samples from two sites were combined. Generally, the selected plaque samples were from sites with the deepest probing depth matched to the sites for Lipid Mediator-SPM metabololipidomics when possible.

Lipid Mediator-SPM Metabololipidomics for Gingival Tissue Samples

The gingival tissue samples were quantitatively analyzed for levels of SPMs, SPM pathway markers and proinflammatory mediators using Lipid Mediator-SPM metabololipidomics as described earlier (21, 23). The quantities of these lipid mediators in gingival tissue samples are presented as mean ± standard deviation ng per mg of the total protein in each sample. The lower limit of quantitation is 0.015ng.

Expression of SPM Receptor Genes in Gingival Tissue Samples via qRT-PCR

Expression of SPM corresponding receptor genes was analyzed with qRT-PCR, utilizing the TaqMan Gene Expression Assays Protocol (Applied Biosystems, Foster City, CA, USA). Details of the experiment were previously described (21). The TaqMan Assay primers included GAPDH (assay ID: Hs99999905_m1, Applied Biosystems) as the housekeeping gene compared to ALX (FPR2) (assay ID: Hs00265954_m1, Applied Biosystems), BLT1 (LTB4R) (assay ID: Hs00175124_m1, Applied Biosystems), ChemR23 (ERV1) (assay ID: Hs01386063_m1, Applied Biosystems), GPR18 (DRV2) (assay ID: Hs00245542_m1, Applied Biosystems), GPR32 (DRV1) (assay ID: Hs00265986_s1, Applied Biosystems), GPR37 (assay ID: Hs0017374_m1, Applied Biosystems),and LGR6 (assay ID: Hs00663887_m1, Applied Biosystems), which were labeled with FAM dye. The calculated cycle threshold (Ct) values from each sample were obtained and relative expression (2^-ΔΔCt) was calculated to determine expression of specific SPM receptor genes between groups.

16S rRNA Sequencing for Subgingival Plaque Samples

DNA extraction of one representative plaque sample was performed utilizing the protocol of a MO BIO PowerSoil DNA Isolation Kit (MO BIO Laboratories, Carlsbad, CA USA) (24). Eukaryotic and prokaryotic cells from the plaque sample were lysed with garnet beads to release DNA. The DNA was isolated on a PowerSoil spin filter and subsequently washed and collected as an eluate in a separate collection tube. The DNA eluate was measured on the UV-Vis spectrophotometer (Nanodrop, Thermo-Fisher Scientific, Waltham, MA, USA) to determine the amount of DNA in the sample. Then DNA samples were processed for 16S rRNA Sequencing (LC Biosciences, Houston, TX). The 16S rRNA V3-V4 regions were amplified via PCR and sequenced on the Illumina MiSeq platform (Illumina Inc., San Diego, CA, USA). The amplification primers (forward primer (338F): 5’-ACTCCTACGGGAGGCAGCAG-3’; reverse primer (806R): 5’-GGACTACHVGGGTWTCTAAT-3’) contained adapters for MiSeq sequencing and single-index barcodes that resulted in PCR products that were pooled and sequenced directly. Read pairs were de-multiplexed based on barcodes and merged. All of the samples were then processed and sequenced together, and the cleaned, merged data was imported into CLC Genomics Workbench with the Microbial Genomics module. The 16S rRNA gene sequences were allocated to specific operational taxonomic units (OTUs) at 98% identity using the Human Oral Microbiome Database (HOMD). OTUs without a match to the HOMD database were given an OTU number, and genus/species were identified by BLAST search against the bacterial 16S rRNA database (25).

Statistical Analysis

The Lipid Mediator-SPM metabololipidomics data and qRT-PCR data were first log transformed, and Wilcoxon rank sum test was used to detect differentially expressed lipid mediators and receptor genes between two conditions except for comparison between the periodontitis prior to SRP group and the periodontitis after SRP group, where Wilcoxon signed rank test was applied. The p-values were adjusted using the False Discovery Rate (FDR) multi-test correction method. Data were considered statistically significant when p<0.05 (21).

For the subgingival microbiome, community diversity (alpha and beta diversity) was assessed using the Microbial Genomics Diversity module of CLC Genomics Workbench. OTUs from the abundance table were aligned using MUSCLE with a required minimum abundance of 100. Aligned OTUs were used to construct a phylogenetic tree using Maximum Likelihood Phylogeny using the Neighbor Joining method and the Jukes Cantor substitution model. Rarefication analysis was done by sub-sampling the OTU abundances in the different samples at a range of depths from 1 to 100,000; the number of different depths sampled was 20, with 100 replicates at each depth. The OTUs with relative abundance less than 0.01% were excluded. Alpha diversity measures were calculated for observed OTUs and Shannon index. Statistical significance in alpha diversity between groups was calculated with Student’s t-test or paired t-test. PERMANOVA (Permutational Multivariate Analysis of Variance) analysis was used to detect significant differences in Beta diversity between groups, and comparisons were visualized using Principal Coordinate Analysis (PCoA). Diversity measures were calculated using the Bray-Curtis formula.

Differences in taxa abundance of subgingival plaque samples between the three groups were identified using DESeq2 (26). OTUs were selected through DESeq2 testing. The OTUs were color-coded according to the phyla they belong to and plotted according to their log2 fold change. The relative fold change was calculated between the three subject groups.

R random Forest package (27) was used to calculate the relative abundance of OTUs. Those OTUs with mean decrease Gini larger than 0.15 were selected for subsequent integration analysis. The lipid mediator levels, SPM receptor gene expression levels and bacterial abundance were first normalized to zero mean and unit variance, and then integrated and analyzed using Data Integration Analysis for Biomarker discovery using a Latent cOmponents (DIABLO) method implemented in R mixOmics package (28). The resulting profile correlations were visualized in DIABLO. Correlations between individual lipid mediators and bacterial species in each subject group were further analyzed using the sparse partial least squares (SPLS) method (29) implemented in the R mixOmics package. The data at the genus level were also analyzed using the same methods.

Results

Composition of the Subgingival Microbiome

Thirteen periodontally healthy subjects and 15 periodontitis subjects before and after treatment were included in this study. All periodontitis subjects were diagnosed with generalized/localized periodontitis, Stage II or III, Grade B or C (Supplementary Table 1). In periodontitis subjects, the mean probing depth of sample collection sites was significantly reduced following non-surgical periodontal therapy (SRP); however, the average PD after treatment remained significantly greater than in health (Table 1).

TABLE 1

Table 1 Characteristics of subjects and sample collection sites.

Three hundred and sixty one OTUs were identified across all samples. The beta-diversity, which represents the variation in microbial composition between groups, was significantly different between the healthy group and the periodontitis prior to SRP group as well as between the healthy group and the periodontitis after SRP group (p=<0.01,<0.01, respectively; Figure 1). Significance was not seen between the periodontitis prior to SRP group and the periodontitis after SRP group (p= 0.57), suggesting periodontal treatment did not entirely change the microbial composition to a healthy condition. The bacterial species richness (the number of OTUs) in the periodontitis prior to SRP group was significantly higher than in the healthy and periodontitis after SRP groups (p=<0.01, 0.03 respectively, Figure 2A). The mean Shannon index (evenness) was the highest in the periodontitis prior to SRP group and lowest in the healthy group. There was only a significant difference in bacterial species evenness between the periodontitis prior to SRP and healthy groups (p=0.04, Figure 2B). These results indicate that species were the most diverse in periodontitis before treatment.

FIGURE 1

Figure 1 Bray-Curtis principal coordinate analysis (PCoA) plot for subgingival microbiome profiles. Axis one accounts for 13% of sample variance, while axis two and three account for 9% and 8% of variance, respectively. One dot represents one sample in each group. This plot demonstrates clusters of samples based on their similarity of microbial composition. These three groups display distinct microbial compositions. (H (green): healthy; P (red): periodontitis before non-surgical therapy; A (blue): periodontitis after non-surgical therapy).

FIGURE 2

Figure 2 Alpha diversity of the subgingival microbiome. (A) Richness of the subgingival microbiome. Richness is represented by the number of operational taxonomic units (OTUs). (B) Evenness of the subgingival microbiome. Evenness is represented by Shannon index. (*p < 0.05; H (green): healthy; P (red): periodontitis before non-surgical therapy; A (blue): periodontitis after non-surgical therapy).

Considering the OTU species assignments, Corynebacterium matruchotii, Fusobacterium nucleatum_subsp._animalis, Fusobacterium nucleatum_subsp._polymorphum, Fusobacterium nucleatum_subsp._vincentii, and Veillonella dispar, were seen in all three groups with varying abundance (Table 2). The top three species that had the highest relative abundance in the healthy group were Neisseria oralis (10%), Corynebacterium matruchotii (6%), and Actinomyces sp._oral_taxon_169 (5%). Periodontitis sites before and after therapy were dominated by Fusobacterium nucleatum subsp. vincentii at 8% and 6%, respectively. While considering species with relative abundance ≥1%, three putative red complex periodontal pathogens, Porphyromonas gingivalis, Tannerella forsythia and Treponema denticola were present in the periodontitis prior to SRP group and Porphyromonas gingivalis and Treponema denticola were present in the periodontitis after SRP group. These species had very low abundance and were infrequently present in the healthy group.

TABLE 2

Table 2 Relative abundance of bacterial species.

Among two-group comparisons, there were 61 bacterial species with significant differences in relative fold change between the healthy and periodontitis prior to SRP groups, 32 bacterial species with significant differences in relative fold change between the healthy and periodontitis after SRP groups and 28 bacterial species with significant differences in relative fold change between the periodontitis prior to SRP and after SRP groups (Figures 3A–C and Supplementary Table 2). The relative amounts of P. gingivalis, T. denticola, and T. forsythia in both the periodontitis prior to SRP and after SRP groups were significantly higher than in the healthy group (Figures 3A, B), but there was no significant difference between the two periodontitis groups. Periodontal health related species, Rothia aeria and Corynebacterium durum, in the healthy group had significantly higher relative fold values than the two periodontitis groups.

FIGURE 3

Figure 3 Differentially abundant operational taxonomic units (OTUs) in three comparisons. Each circle represents one OTU (bacterial species). The OTUs are color-coded based on the phylum to which they belong and plotted based on their log2 fold change in each comparison. The size of the circle is proportional to the mean count of each species. All differentially abundant OTUs in each comparison are listed. (A) Healthy (H) vs. Periodontitis before non-surgical therapy (P); (B) Healthy (H) vs. Periodontitis after non-surgical therapy (A); (C) Periodontitis before non-surgical therapy (P) vs. Periodontitis after non-surgical therapy (A). In these comparisons, the second group is the reference group. The fold change is calculated by dividing the bacterial abundance in the first group by the bacterial abundance in the second group (H/P, H/A, P/A respectively). These graphs demonstrate different microbial compositions between the three groups.

Correlations Between Lipid Mediator, SPM Receptor Gene and Subgingival Microbiome

Previously, we demonstrated that the level of six lipid mediators (5-HETE, 15-HETE, 15(S)-HEPE, 4-HDHA, 7-HDHA, 17-HDHA) and expression of D-series resolvins corresponding receptor genes (GPR18, GPR32) were significantly altered by periodontitis treatment (21). In the present study, data of lipid mediator levels, receptor gene expression and bacterial species abundance were integrated and analyzed by the DIABLO method to demonstrate profile correlations. When all samples from the three subject groups are included, the correlation coefficient between subgingival microbiome and lipid mediator profiles was 0.61, the correlation coefficient between subgingival microbiome and receptor gene profiles was 0.45 and the correlation coefficient between lipid mediator and receptor gene profiles was 0.34. (Figure 4A). As the results at the species level, the correlation between subgingival microbiome and lipid mediator profiles was also the highest at the genus level. (correlation coefficient= 0.5) (Supplementary Figure 1A). Since the subgingival microbiome and lipid mediator profiles had the highest correlation, the subgingival microbiome-lipid mediator profile correlation was further assessed in each group and the correlation coefficients were 0.65, 0.58 and 0.43 in the healthy, periodontitis prior to SRP and periodontitis after SRP groups respectively (Figure 4B). The correlation patterns at the genus level in each group were also similar to those at the species level (Supplementary Figure 1B).

FIGURE 4

Figure 4 Correlations between lipid mediator profile, specialized pro-resolving lipid mediator (SPM) receptor gene profile and subgingival microbiome profile. (A) The component correlation plots represent subgingival microbiome-lipid mediator profile (upper-middle box), subgingival microbiome-receptor gene profile (upper-right box), and lipid mediator-receptor gene profile (middle-right box). The correlation coefficient of these profiles: subgingival microbiome-lipid mediator (middle-left box)= 0.61; subgingival microbiome-receptor gene (lower-left box)= 0.45; lipid mediator-receptor gene (lower-middle box)= 0.34. (B) The correlation plot demonstrates the subgingival microbiome-lipid mediator correlation patterns in each subject group. The correlation coefficients for the H, P and A groups are 0.65, 0.58 and 0.43, respectively. Most of the subject’s subgingival microbiome-lipid mediator profiles (subjects are labeled in the plot) in periodontitis before treatment move toward the healthy patterns after treatment. (H (green): healthy; P (red): periodontitis before non-surgical therapy; A (blue): periodontitis after non-surgical therapy; one dot represents one subject in each group; numbers on axes represent relative levels of bacterial species, lipid mediators or receptor gene expression; the analyses are performed using the DIABLO method.)

Specific bacterial species and lipid mediators were correlated (correlation coefficient ≥0.5) while including all samples of the three groups (Figure 5). Then, SPLS analysis was conducted to identify individual bacterial species and lipid mediators contributing to the subgingival microbiome-lipid mediator profile correlation in each group. Using a correlation coefficient cut off of ≥0.6, five bacterial species and six lipid mediators were correlated in the healthy group, five bacterial species and 11 lipid mediators were correlated in the periodontitis prior to SRP group and nine bacterial species and nine lipid mediators were correlated in the periodontitis after SRP group (Supplementary Table 3 and Figure 6). These correlations indicate that the relative abundance of these bacteria species and the level of these lipid mediators are highly associated.

FIGURE 5

Figure 5 Correlations between individual lipid mediators and bacterial species. The circos plot shows correlations between lipid mediators (yellow) and bacterial species (pink). Three lines outside the circle represent the relative levels of lipid mediators or bacterial species (H (green): healthy; P (red): periodontitis before non-surgical therapy; A (blue): periodontitis after non-surgical therapy). Purple lines inside the circle represent positive correlations between lipid mediators and bacterial species with correlation coefficient ≥0.5. None of the bacterial species and lipid mediators has a significant negative correlation coefficient ≤-0.5. The analysis is performed using the DIABLO method.

FIGURE 6

Figure 6 Correlations between bacterial species and lipid mediators in the three groups. (A) correlated lipid mediators and bacterial species in the healthy group. (B) correlated lipid mediators and bacterial species in the periodontitis before non-surgical therapy group. (C) correlated lipid mediators and bacterial species in the periodontitis after non-surgical therapy group. The analysis is performed using the Sparse Partial Least Squares (SPLS) method. Bacterial species are listed in alphabetical order. Lipid mediators are grouped by biosynthetic pathways; leukotriene pathway marker: 5(S),6(R)-DiHETE; lipoxygenase pathway marker: 5(S),12(S)-DiHETE; prostaglandin pathway marker: 12(S)-HHTrE; E-series resolvins and pathway markers (derived from omega-3 EPA): 15(S)-HEPE, 18-HEPE, RvE3; D-series resolvins and pathway markers (derived from omega-3 DHA or DPA): RvD1, 7-HDHA, 13-HDHA; Maresins and pathway markers (derived from omega-3 DHA or DPA): 14-HDHA, MaR1, 7(S)-MaR1; pro-inflammatory lipid mediators: LTB4, PGD2, PGE2, PGF2a. Some pathway markers are involved in multiple pathways but only one major pathway is listed. All correlated lipid mediators-bacterial species have absolute correlation coefficients ≥0.60.

When comparing these correlated species in periodontitis before treatment to healthy controls, only Corynebacterium durum, a periodontal health related species, was identified in both conditions, but correlated with different lipids. In health, Corynebacterium durum was positively correlated with the leukotriene pathway marker, 5(S)6(R)-DiHETE, and one of the E-series resolvins, RvE3. In periodontitis before treatment, Corynebacterium durum was negatively correlated with the resolvin pathway marker, 15(S)-HEPE, and positively correlated with the proinflammatory mediator, LTB4, and SPMs, including 7(S)-Maresin1, Maresin1 and RvD1. Both conditions had five lipid mediators, 5(S)6(R)-DiHETE, 15(S)-HEPE, 7-HDHA, 13-HDHA and 14-HDHA, correlated with different bacterial species.

When comparing these correlated species in periodontitis after treatment to healthy controls, Treponema socranskii, a periodontal pathogen, was identified in both conditions. This bacterial species was negatively correlated to 5(S)6(R)-DiHETE in the healthy group and positively correlated with different lipids in the periodontitis after SRP group. The E-series resolvin, RvE3, was the only lipid mediator identified in both conditions.

When comparing these correlated species in periodontitis before treatment to after treatment, a bacterial species, Anaeroglobus geminatus, was identified in both conditions that positively correlated with different lipid mediators. Both conditions had four lipid mediators, 5(S)12(S)-DiHETE, RvD1, Maresin 1 and LTB4, correlated with different bacteria species. Among the nine bacterial species identified in the periodontitis after SRP group, four Selenomonas species (Selenomonas sp._oral_taxon_136, Selenomonas sp._oral_taxon_137, Selenomonas sp._oral_taxon_138, Selenomonas sp._oral_taxon_479) are highly correlated with multiple lipid mediators.

Discussion

Many studies demonstrate that shifts in the bacterial load and bacterial composition within the subgingival microbiota are correlated with periodontal conditions (30, 31). Associations between periodontal inflammation and molecular-level host responses (32, 33) or metabolite profiles (34, 35) have also been investigated. However, there is a remarkable lack of studies assessing clinical host-microbiota relationship in periodontitis. The current study showed different correlations between levels of lipid mediators, expression of SPM receptor genes and abundance of subgingival bacterial species in health vs. disease. To the best of our knowledge, this is the first study demonstrating that the molecular profiles for lipid-mediated resolution of inflammation activity are associated with microbial composition shifts in human periodontitis. The correlations between specific lipid mediators and bacterial species indicate their potential interactions, which are important for deciphering the mechanisms of resolution of inflammation in the pathogenesis of periodontitis.

In this study, most of the periodontitis subjects have severe periodontitis (Stage III, Grade C). In addition to the lack of maintenance or poor self-care, these patients may have abnormal host responses causing severe, rapid periodontal destruction and dysbiotic subgingival microbiota. Current results indicate that failure of resolution associated with the imbalanced lipid mediator profiles could be the potential mechanisms (21). This study compared samples from the diseases sites of the periodontitis subjects to the healthy sites of the subjects without periodontitis because lipid mediator levels and profiles are distinct at the subject level (16, 17). However, it would be interesting to analyze samples from both diseases sites and healthy sites of the same periodontitis subjects to evaluate the subgingival microbiome-lipid mediator correlations at different local environments.

According to the profile correlation analysis, the correlation between lipid mediator profile and receptor gene profile was lower than both the subgingival microbiome-lipid mediator profile correlation and subgingival microbiome-receptor gene profile correlation. The low lipid mediator-receptor gene profile correlation supports the hypothesis that an imbalance between lipid mediator levels and receptor expression results in failure of resolution of inflammation in periodontitis (21). It is known that expression of SPM corresponding receptors is important for resolution of inflammation. Deficient expression of these receptors has been associated with increased inflammation in peritonitis (36, 37), paw edema (38), and microbial sepsis (39) in pre-clinical models. Overexpression of SPM receptors, such as ChemR23 (ERV1), was associated with immune responses favorable for inflammation resolution in a dorsal air pouch model (40).

As with the clustering in both lipid mediators and SPM receptor gene expression (21), the clusters of subgingival microbiome were associated with inflammatory conditions in gingiva (Figure 1). Current results of 16S rRNA sequencing show that periodontitis subjects had a higher abundance of periodontal pathogens and microbial diversity than healthy subjects. Healthy subjects also have a higher abundance of periodontal health related bacterial species. These findings were similar to published results in the literature (41–45).

The highest subgingival microbiome-lipid mediator profile correlation (Figure 4A) suggests that the composition of the subgingival microbiota could be significantly affected by periodontal inflammation mediated by SPMs and other relevant lipid mediators (13, 14). The correlation pattern in the periodontitis after SRP group was closer to the healthy group than the periodontitis prior to SRP group suggesting periodontal treatment moves profiles of lipid mediators and subgingival microbiota toward non-diseased profiles (Figure 4B). Non-surgical periodontal therapy not only improves clinical parameters but also changes molecular profiles closer to the homeostatic condition. It is possible that local debridement stimulates SPM-mediated resolution of inflammation influencing microbiota shifts in the re-establishment of the subgingival biofilm. In the correlation analysis between individual lipid mediators and bacterial species, most of the correlated lipid and bacterial species are specific for the subject group and not shared by other groups. These results also support that the subgingival microbiome and lipid mediator profiles are specifically associated with periodontal inflammatory conditions.

The correlation analysis identified some bacterial species infrequently discussed in the literature and several lipid mediators potentially important for resolution of periodontal inflammation (Figure 6). A periodontal health related bacterial species, Corynebacterium durum, had relatively low abundance in periodontitis, but its positive correlation with several SPMs indicate that gingival tissues produce SPMs to regulate periodontal inflammation resulting in the presence of this periodontal health related bacterial species. In periodontitis before treatment, one bacterial species, Actinobaculum sp._oral_taxon_848 (46), had the greatest number of correlated resolution related lipid mediators, 5(S)12(S)-DiHETE, 5(S)6(R)-DiHETE, 18-HEPE, 7-HDHA, 13-HDHA, 14-HDHA, Maresin1, and higher correlations than other species. A recently identified bacterial species, Anaeroglobus geminatus (47), related to periodontitis and rheumatoid arthritis, was identified in periodontitis before and after treatment, but was highly correlated with different lipid mediators in the two inflammatory conditions. These bacterial species could be important for periodontal inflammation and deserve further investigation.

Selenomonas species were dominant in the identified bacterial species correlated with multiple lipid mediators in the periodontitis after SRP group. Although Selenomonas species are not frequently discussed in the literature, some studies show associations between the presence, as well as abundance, of Selenomonas species and periodontitis (48–50), and the fundamental role of Selenomonas species in the structure of subgingival biofilm (51). Tanner and coworkers found that Selenomonas noxia was associated with the progression of periodontitis (52). In the current study, most of the identified Selenomonas species had the highest relative abundance in the periodontitis after SRP group although the relative abundance was generally low (<1%, Supplementary Table 4). The abundance of Selenomonas species was highly and positively associated with levels of two SPMs, RvE3 and RvD1 in periodontitis after treatment. An in vitro study investigating the impact of dietary lipid supplements on the microbiota in the rumen ecosystem showed that fish oil in the diet is associated with an increased amount of Selenomonas species in fermenters (53) suggesting omega-3 fatty acids influence the composition of Selenomonas species and other bacterial species. The current findings support the association between the abundance of Selenomonas species and levels of lipid mediators in resolution of periodontal inflammation following treatment. Although the exact role of Selenomonas species in the resolution of periodontal inflammation has to be further studied, Selenomonas species appear to be sensitive to the resolution of inflammation after treatment and the changes of their abundance could be critically associated with the changes of periodontal inflammatory conditions.

As with Selenomonas species, two periodontitis pathogens, Treponema socranskii (50, 54, 55) and Tannerella sp._oral_taxon_808 (56), and an oral biofilm early colonizer species, Actinomyces gerencseriae, were also correlated with RvE3, RvD1 and other lipid mediators, including proinflammatory mediators, LTB4, PGD2, PGE2 and PGF2a, in periodontitis after treatment. These results indicate that the resolution phase of periodontal inflammation induced by periodontal treatment includes a combination of proresolving and proinflammatory mediators highly correlated with non-periodontitis or periodontitis associated bacterial species, explaining the specific phase of inflammation resolution between untreated periodontitis and healthy phases.

We previously reported that levels of six lipid mediators in periodontitis before treatment were significantly higher than in periodontitis after treatment (21). In this study, among the six lipid mediators, only 15(S)-HEPE, a pathway marker for E-series resolvins, was negatively correlated with Corynebacterium durum, and 7-HDHA, a pathway marker for D-series resolvins, was positively correlated with Actinobaculum sp._oral_taxon_848 in periodontitis before treatment. These two SPM pathway markers appear to be more associated with the subgingival microbiota in periodontitis than the others. Two lipids derived from arachidonic acid, 5(S),12(S)-DiHETE and 5(S)6(R)-DiHETE, were correlated with multiple bacterial species in the three subject groups. 5(S),12(S)-DiHETE is an epimer of leukotriene B₄ weakly chemotactic for neutrophils (57). 5(S)6(R)-DiHETE is a hydrolysis product of leukotriene A₄ and can bind to the leukotriene receptor (58). The identification of these two lipids in the correlation analysis suggests their relevant synthetic pathways for leukotrienes could be actively associated with changes in the subgingival microbial composition.

Resolution of inflammation induced by SPMs is critical for periodontal regeneration and healing given chronic and excessive inflammation is detrimental to homeostasis (9). Several SPMs, including RvD1, RvE3, Maresin1 and 7(S)-Maresin 1, were identified in the current correlation analysis. Resolvins are known to resolve inflammation, stimulate periodontal tissue healing and regenerate lost alveolar bone. RvE1 induces periodontal regeneration in pre-clinical periodontitis models (12–14, 59, 60), promotes regenerative properties of periodontal ligament (PDL) stem cells (61), and regulates osteoclasts as well as osteoblasts favorably for bone preservation (62–64). RvE3 demonstrates potent resolution properties by limiting neutrophil infiltration in the peritonitis model and inhibiting neutrophil chemotaxis (65–67). RvD1 inhibits osteoclast differentiation and reduces bone destruction in an arthritis preclinical model (68), inhibits production of pro-inflammatory cytokines from gingival fibroblasts (69) and enhances PDL fibroblast proliferation as well as wound closure (70). Maresins can promote periodontal tissue healing and regulate periodontal inflammation. MaR1 promotes survival and inhibits apoptosis in PDL cells (71), induces regenerative properties of PDL stem cells (61), and restores phagocytic capacity of macrophages from periodontitis patients (72). Both RvD1 and RvE3 were highly correlated with multiple bacterial species in periodontitis after treatment suggesting their important role in resolution of periodontal inflammation associated with subgingival microbiota shift. Their mean levels also increased in periodontitis after treatment as compared to periodontitis before treatment (21). According to the current findings and molecular actions described in the literature, RvD1 and RvE3 potentially could be utilized to treat periodontitis.

In general, results of the correlation analyses indicate significant associations between lipid mediator levels and subgingival bacterial abundance in periodontal inflammation. It is possible that immune responses mediated by SPMs and other lipid mediators drive the shifts of the subgingival microbiota, which can be supported by preclinical models showing oral microbiota shifts following the local administration of SPMs (13, 14). According to in vitro experiments, SPMs do not have direct anti-bacterial properties (73). However, the other direct effects of these lipid mediators on bacterial growth or other activities could not be excluded (74). On the other hand, oral bacteria might affect lipid synthesis and metabolism. It is known that the metabolites from gut microbiota affects host lipid metabolism and lipid composition (75). Oral microbiota could do the same. Additionally, periodontal tissues could produce more lipid mediators to regulate excessive inflammation induced by the increased amount of bacteria in periodontitis (73). These potential mechanisms for subgingival microbiome-lipid mediator correlations have to be further studied. All of these possibilities indicate the important role of lipid mediator profiles in dysbiosis of the subgingival microbiota in periodontitis.

With some limitations, current results should be carefully interpreted. Profiles of lipid mediators and receptor genes were derived from gingival tissues. Since lipid mediators can be produced by and receptor genes can be expressed on a variety of cells, including neutrophils, macrophages, fibroblasts, osteoblasts and osteoclasts (14, 62–64, 69, 70, 76), these molecular profiles should be interpreted as a snapshot of inflammatory reactions in gingiva, but cannot fully explain cellular activities. Additionally, results of the correlation analysis could be biased by the limited number of subjects and subject variability. It is necessary to conduct in vitro or in vivo experiments to further investigate these correlations and potential interactions between bacterial species and lipid mediators.

In conclusion, profiles of lipid mediators, receptor genes and the subgingival microbiome were distinct in different periodontal inflammatory conditions. The highly correlated lipid mediator and subgingival microbiome profiles and specific correlations between individual lipid mediators and bacterial species indicate that periodontal inflammation regulated by lipid mediators drives the shifts of the subgingival microbiota. Elucidation of these correlations facilitates the understanding of resolution of inflammation in periodontitis and identification of potential biological targets for the development of novel therapies.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI BioProject ID: PRJNA726254.

Ethics Statement

The studies involving human participants were reviewed and approved by UTHealth Committee for the Protection of Human Subjects. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

C-TL and TD contributed conception and design of the study. C-TL and BF collected clinical samples. C-TL, GT and BF processed clinical samples and performed experiments. C-TL, RL, LZ, GT, WZ and KM performed the data analysis. C-TL, RL, LZ, GT, WZ, KM, NA, and TD contributed to data interpretation. C-TL wrote the first draft of the manuscript. All authors wrote sections of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported in part by UTHealth School of Dentistry Grants Program in Dental Research (C-TL), UTHealth the Center for Clinical and Translational Sciences (CCTS) Translational Technologies Core Laboratories Awards (C-TL), USPHS grant DE025020 from the National Institute of Dental and Craniofacial Research (TD), CPRIT RP170668 (WZ) and NIH 1 UL1 TR003167 01 (WZ).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank Dr. Srinivas Ayilavarapu, Department of Periodontics and Oral Hygiene, The University of Texas Health Science Center at Houston School of Dentistry, for his professional opinions on this project.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.691216/full#supplementary-material

Abbreviations

5(S),6(R)-DiHETE, 5(S),6(R)-dihydroxyeicosatetraenoic acid; 5(S),12(S)-DiHETE, 5(S),12(S)-dihydroxy-6E,8Z,11E,14Z-eicosatetraenoic acid; 7-HDHA, 7-hydroxydocosahexaenoic acid; 12(S)-HHTrE, 12(S)-hydroxyheptadecatrienoic acid; 13-HDHA, 13-hydroxydocosahexaenoic acid; 14-HDHA, 14-hydroxydocosahexaenoic acid; 15(S)-HEPE, 15(S)-hydroxyeicosapentaenoic acid; 18-HEPE, 18-hydroxyeicosapentaenoic acid; DHA, docosahexaenoic acid; DPA, docosapentaenoic acid; EDTA, ethylenediaminetetraacetic acid; EPA, eicosapentaenoic acid; LTB4, leukotriene B4; MaR1, maresin 1, 7(S)-MaR1, 7-epi maresin 1; PGD2, prostaglandin D2; PGE2, prostaglandin E2; PGF2a, prostaglandin F2a; RvD1, resolvin D1; RvE3, resolvin E3.

References

1. Eke PI, Thornton-Evans GO, Wei L, Borgnakke WS, Dye BA, Genco RJ. Periodontitis in US Adults: National Health and Nutrition Examination Survey 2009-2014. J Am Dent Assoc (2018) 149(7):576–88 e6. doi: 10.1016/j.adaj.2018.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Meyle J, Chapple I. Molecular Aspects of the Pathogenesis of Periodontitis. Periodontol 2000 (2015) 69(1):7–17. doi: 10.1111/prd.12104

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Socransky SS, Haffajee AD, Cugini MA, Smith C, Kent RL Jr. Microbial Complexes in Subgingival Plaque. J Clin Periodontol (1998) 25(2):134–44. doi: 10.1111/j.1600-051x.1998.tb02419.x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Marsh PD, Devine DA. How is the Development of Dental Biofilms Influenced by the Host? J Clin Periodontol (2011) 38(Suppl 11):28–35. doi: 10.1111/j.1600-051X.2010.01673.x

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Van Dyke TE, Bartold PM, Reynolds EC. The Nexus Between Periodontal Inflammation and Dysbiosis. Front Immunol (2020) 11:511. doi: 10.3389/fimmu.2020.00511

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Chen C, Hemme C, Beleno J, Shi ZJ, Ning D, Qin Y, et al. Oral Microbiota of Periodontal Health and Disease and Their Changes After Nonsurgical Periodontal Therapy. ISME J (2018) 12(5):1210–24. doi: 10.1038/s41396-017-0037-1

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Belstrom D, Grande MA, Sembler-Moller ML, Kirkby N, Cotton SL, Paster BJ, et al. Influence of Periodontal Treatment on Subgingival and Salivary Microbiotas. J Periodontol (2018) 89(5):531–9. doi: 10.1002/JPER.17-0377

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Serhan CN, Chiang N, Dalli J. The Resolution Code of Acute Inflammation: Novel Pro-Resolving Lipid Mediators in Resolution. Semin Immunol (2015) 27(3):200–15. doi: 10.1016/j.smim.2015.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Van Dyke TE. Pro-Resolving Mediators in the Regulation of Periodontal Disease. Mol Aspects Med (2017) 58:21–36. doi: 10.1016/j.mam.2017.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Gewirtz AT, Collier-Hyams LS, Young AN, Kucharzik T, Guilford WJ, Parkinson JF, et al. Lipoxin A4 Analogs Attenuate Induction of Intestinal Epithelial Proinflammatory Gene Expression and Reduce the Severity of Dextran Sodium Sulfate-Induced Colitis. J Immunol (2002) 168(10):5260–7. doi: 10.4049/jimmunol.168.10.5260

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tang Y, Zhang MJ, Hellmann J, Kosuri M, Bhatnagar A, Spite M. Proresolution Therapy for the Treatment of Delayed Healing of Diabetic Wounds. Diabetes (2013) 62(2):618–27. doi: 10.2337/db12-0684

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Osorio Parra MM, Elangovan S, Lee CT. Specialized Pro-Resolving Lipid Mediators in Experimental Periodontitis: A Systematic Review. Oral Dis (2019) 25(5):1265–76. doi: 10.1111/odi.12979

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Hasturk H, Kantarci A, Goguet-Surmenian E, Blackwood A, Andry C, Serhan CN, et al. Resolvin E1 Regulates Inflammation At the Cellular and Tissue Level and Restores Tissue Homeostasis In Vivo. J Immunol (2007) 179(10):7021–9. doi: 10.4049/jimmunol.179.10.7021

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Lee CT, Teles R, Kantarci A, Chen T, McCafferty J, Starr JR, et al. Resolvin E1 Reverses Experimental Periodontitis and Dysbiosis. J Immunol (2016) 197(7):2796–806. doi: 10.4049/jimmunol.1600859

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Arnardottir H, Orr SK, Dalli J, Serhan CN. Human Milk Proresolving Mediators Stimulate Resolution of Acute Inflammation. Mucosal Immunol (2016) 9(3):757–66. doi: 10.1038/mi.2015.99

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Colas RA, Shinohara M, Dalli J, Chiang N, Serhan CN. Identification and Signature Profiles for Pro-Resolving and Inflammatory Lipid Mediators in Human Tissue. Am J Physiol Cell Physiol (2014) 307(1):C39–54. doi: 10.1152/ajpcell.00024.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Elabdeen HR, Mustafa M, Szklenar M, Ruhl R, Ali R, Bolstad AI. Ratio of Pro-Resolving and Pro-Inflammatory Lipid Mediator Precursors as Potential Markers for Aggressive Periodontitis. PloS One (2013) 8(8):e70838. doi: 10.1371/journal.pone.0070838

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tobon-Arroyave SI, Isaza-Guzman DM, Gomez-Ortega J, Florez-Alzate AA. Salivary Levels of Specialized Pro-Resolving Lipid Mediators as Indicators of Periodontal Health/Disease Status. J Clin Periodontol (2019) 46(10):978–90. doi: 10.1111/jcpe.13173

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Fredman G, Hellmann J, Proto JD, Kuriakose G, Colas RA, Dorweiler B, et al. An Imbalance Between Specialized Pro-Resolving Lipid Mediators and Pro-Inflammatory Leukotrienes Promotes Instability of Atherosclerotic Plaques. Nat Commun (2016) 7:12859. doi: 10.1038/ncomms12859

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Colas RA, Nhat LTH, Thuong NTT, Gomez EA, Ly L, Thanh HH, et al. Proresolving Mediator Profiles in Cerebrospinal Fluid are Linked With Disease Severity and Outcome in Adults With Tuberculous Meningitis. FASEB J (2019) 33(11):13028–39. doi: 10.1096/fj.201901590R

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ferguson B, Bokka NR, Maddipati KR, Ayilavarapu S, Weltman R, Zhu L, et al. Distinct Profiles of Specialized Pro-Resolving Lipid Mediators and Corresponding Receptor Gene Expression in Periodontal Inflammation. Front Immunol (2020) 11:1307. doi: 10.3389/fimmu.2020.01307

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Jepsen S, Caton JG, Albandar JM, Bissada NF, Bouchard P, Cortellini P, et al. Periodontal Manifestations of Systemic Diseases and Developmental and Acquired Conditions: Consensus Report of Workgroup 3 of the 2017 World Workshop on the Classification of Periodontal and Peri-Implant Diseases and Conditions. J Periodontol (2018) 89 Suppl 1:S237–S48. doi: 10.1002/JPER.17-0733

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Norris PC, Skulas-Ray AC, Riley I, Richter CK, Kris-Etherton PM, Jensen GL, et al. Identification of Specialized Pro-Resolving Mediator Clusters From Healthy Adults After Intravenous Low-Dose Endotoxin and Omega-3 Supplementation: A Methodological Validation. Sci Rep (2018) 8(1):18050. doi: 10.1038/s41598-018-36679-4

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lee CT, Tribble G, Osborne A, Wu E. Dna Extraction From Low Biomass Plaque Using Two Homogenization Methods for 16S rRNA Sequencing: A Comparative Study. J Periodontics Implant Dent (2018) 1(2):74–80. doi: 10.3966/261634032018120102004

CrossRef Full Text | Google Scholar

25. Belstrom D, Sembler-Moller ML, Grande MA, Kirkby N, Cotton SL, Paster BJ, et al. Microbial Profile Comparisons of Saliva, Pooled and Site-Specific Subgingival Samples in Periodontitis Patients. PloS One (2017) 12(8):e0182992. doi: 10.1371/journal.pone.0182992

PubMed Abstract | CrossRef Full Text | Google Scholar

26. 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. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Liaw A, Wiener M. Classification and Regression by Randomforest. R News (2002) 2(3):18–22.

Google Scholar

28. Rohart F, Gautier B, Singh A, Le Cao KA. mixOmics: An R Package for ‘Omics Feature Selection and Multiple Data Integration. PloS Comput Biol (2017) 13(11):e1005752. doi: 10.1371/journal.pcbi.1005752

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Le Cao KA, Martin PG, Robert-Granie C, Besse P. Sparse Canonical Methods for Biological Data Integration: Application to a Cross-Platform Study. BMC Bioinf (2009) 10:34. doi: 10.1186/1471-2105-10-34

CrossRef Full Text | Google Scholar

30. Teles R, Teles F, Frias-Lopez J, Paster B, Haffajee A. Lessons Learned and Unlearned in Periodontal Microbiology. Periodontology 2000 (2013) 62(1):95–162. doi: 10.1111/prd.12010

PubMed Abstract | CrossRef Full Text | Google Scholar

31. 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. doi: 10.1128/JB.183.12.3770-3783.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Davanian H, Stranneheim H, Bage T, Lagervall M, Jansson L, Lundeberg J, et al. Gene Expression Profiles in Paired Gingival Biopsies From Periodontitis-Affected and Healthy Tissues Revealed by Massively Parallel Sequencing. PloS One (2012) 7(9):e46440. doi: 10.1371/journal.pone.0046440

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Kebschull M, Guarnieri P, Demmer RT, Boulesteix AL, Pavlidis P, Papapanou PN. Molecular Differences Between Chronic and Aggressive Periodontitis. J Dental Res (2013) 92(12):1081–8. doi: 10.1177/0022034513506011

CrossRef Full Text | Google Scholar

34. Barnes VM, Teles R, Trivedi HM, Devizio W, Xu T, Mitchell MW, et al. Acceleration of Purine Degradation by Periodontal Diseases. J Dental Res (2009) 88(9):851–5. doi: 10.1177/0022034509341967

CrossRef Full Text | Google Scholar

35. Barnes VM, Ciancio SG, Shibly O, Xu T, Devizio W, Trivedi HM, et al. Metabolomics Reveals Elevated Macromolecular Degradation in Periodontal Disease. J Dental Res (2011) 90(11):1293–7. doi: 10.1177/0022034511416240

CrossRef Full Text | Google Scholar

36. Chiang N, Dalli J, Colas RA, Serhan CN. Identification of Resolvin D2 Receptor Mediating Resolution of Infections and Organ Protection. J Exp Med (2015) 212(8):1203–17. doi: 10.1084/jem.20150225

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Norling LV, Dalli J, Flower RJ, Serhan CN, Perretti M. Resolvin D1 Limits Polymorphonuclear Leukocyte Recruitment to Inflammatory Loci: Receptor-Dependent Actions. Arterioscler Thromb Vasc Biol (2012) 32(8):1970–8. doi: 10.1161/ATVBAHA.112.249508

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Dufton N, Hannon R, Brancaleone V, Dalli J, Patel HB, Gray M, et al. Anti-Inflammatory Role of the Murine Formyl-Peptide Receptor 2: Ligand-Specific Effects on Leukocyte Responses and Experimental Inflammation. J Immunol (2010) 184(5):2611–9. doi: 10.4049/jimmunol.0903526

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Gobbetti T, Coldewey SM, Chen J, McArthur S, le Faouder P, Cenac N, et al. Nonredundant Protective Properties of FPR2/ALX in Polymicrobial Murine Sepsis. Proc Natl Acad Sci USA (2014) 111(52):18685–90. doi: 10.1073/pnas.1410938111

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Herrera BS, Hasturk H, Kantarci A, Freire MO, Nguyen O, Kansal S, et al. Impact of Resolvin E1 on Murine Neutrophil Phagocytosis in Type 2 Diabetes. Infect Immun (2015) 83(2):792–801. doi: 10.1128/IAI.02444-14

PubMed Abstract | CrossRef Full Text | Google Scholar

41. 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. doi: 10.1038/ismej.2012.174

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Szafranski SP, Wos-Oxley ML, Vilchez-Vargas R, Jauregui R, Plumeier I, Klawonn F, et al. High-Resolution Taxonomic Profiling of the Subgingival Microbiome for Biomarker Discovery and Periodontitis Diagnosis. Appl Environ Microbiol (2015) 81(3):1047–58. doi: 10.1128/AEM.03534-14

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Griffen AL, Beall CJ, Campbell JH, Firestone ND, Kumar PS, Yang ZK, et al. Distinct and Complex Bacterial Profiles in Human Periodontitis and Health Revealed by 16S Pyrosequencing. ISME J (2012) 6(6):1176–85. doi: 10.1038/ismej.2011.191

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Duran-Pinedo AE, Chen T, Teles R, Starr JR, Wang X, Krishnan K, et al. Community-Wide Transcriptome of the Oral Microbiome in Subjects With and Without Periodontitis. ISME J (2014) 8(8):1659–72. doi: 10.1038/ismej.2014.23

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Mark Welch JL, Rossetti BJ, Rieken CW, Dewhirst FE, Borisy GG. Biogeography of a Human Oral Microbiome At the Micron Scale. Proc Natl Acad Sci USA (2016) 113(6):E791–800. doi: 10.1073/pnas.1522149113

PubMed Abstract | CrossRef Full Text | Google Scholar

46. LaMonte MJ, Genco RJ, Buck MJ, McSkimming DI, Li L, Hovey KM, et al. Composition and Diversity of the Subgingival Microbiome and Its Relationship With Age in Postmenopausal Women: An Epidemiologic Investigation. BMC Oral Health (2019) 19(1):246. doi: 10.1186/s12903-019-0906-2

PubMed Abstract | CrossRef Full Text | Google Scholar

47. 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. doi: 10.1002/art.34539

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Kumar PS, Griffen AL, Moeschberger ML, Leys EJ. Identification of Candidate Periodontal Pathogens and Beneficial Species by Quantitative 16s Clonal Analysis. J Clin Microbiol (2005) 43(8):3944–55. doi: 10.1128/JCM.43.8.3944-3955.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Goncalves LF, Fermiano D, Feres M, Figueiredo LC, Teles FR, Mayer MP, et al. Levels of Selenomonas Species in Generalized Aggressive Periodontitis. J Periodontal Res (2012) 47(6):711–8. doi: 10.1111/j.1600-0765.2012.01485.x

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Nibali L, Sousa V, Davrandi M, Spratt D, Alyahya Q, Dopico J, et al. Differences in the Periodontal Microbiome of Successfully Treated and Persistent Aggressive Periodontitis. J Clin Periodontol (2020) 47(8):980–90. doi: 10.1111/jcpe.13330

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Drescher J, Schlafer S, Schaudinn C, Riep B, Neumann K, Friedmann A, et al. Molecular Epidemiology and Spatial Distribution of Selenomonas Spp. in Subgingival Biofilms. Eur J Oral Sci (2010) 118(5):466–74. doi: 10.1111/j.1600-0722.2010.00765.x

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Tanner A, Maiden MF, Macuch PJ, Murray LL, Kent RL Jr. Microbiota of Health, Gingivitis, and Initial Periodontitis. J Clin Periodontol (1998) 25(2):85–98. doi: 10.1111/j.1600-051x.1998.tb02414.x

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Potu RB, AbuGhazaleh AA, Hastings D, Jones K, Ibrahim SA. The Effect of Lipid Supplements on Ruminal Bacteria in Continuous Culture Fermenters Varies With the Fatty Acid Composition. J Microbiol (2011) 49(2):216–23. doi: 10.1007/s12275-011-0365-1

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Visser MB, Ellen RP. New Insights Into the Emerging Role of Oral Spirochaetes in Periodontal Disease. Clin Microbiol Infect (2011) 17(4):502–12. doi: 10.1111/j.1469-0691.2011.03460.x

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Moter A, Riep B, Haban V, Heuner K, Siebert G, Berning M, et al. Molecular Epidemiology of Oral Treponemes in Patients With Periodontitis and in Periodontitis-Resistant Subjects. J Clin Microbiol (2006) 44(9):3078–85. doi: 10.1128/JCM.00322-06

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Beall CJ, Campbell AG, Griffen AL, Podar M, Leys EJ. Genomics of the Uncultivated, Periodontitis-Associated Bacterium Tannerella Sp. Bu045 (Oral Taxon 808). mSystems (2018) 3(3) 3(3):e000018–18. doi: 10.1128/mSystems.00018-18

CrossRef Full Text | Google Scholar

57. Lee TH, Menica-Huerta JM, Shih C, Corey EJ, Lewis RA, Austen KF. Characterization and Biologic Properties of 5,12-Dihydroxy Derivatives of Eicosapentaenoic Acid, Including Leukotriene B5 and the Double Lipoxygenase Product. J Biol Chem (1984) 259(4):2383–9. doi: 10.1016/S0021-9258(17)43363-1

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Muller A, Rechencq E, Kugel C, Lellouche JP, Beaucourt JP, Niel G, et al. Comparative Biological Activities of the Four Synthetic (5,6)-Dihete Isomers. Prostaglandins (1989) 38(6):635–44. doi: 10.1016/0090-6980(89)90046-4

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Hasturk H, Abdallah R, Kantarci A, Nguyen D, Giordano N, Hamilton J, et al. Resolvin E1 (Rve1) Attenuates Atherosclerotic Plaque Formation in Diet and Inflammation-Induced Atherogenesis. Arteriosclerosis Thrombosis Vasc Biol (2015) 35(5):1123–33. doi: 10.1161/ATVBAHA.115.305324

CrossRef Full Text | Google Scholar

60. Hasturk H, Kantarci A, Ohira T, Arita M, Ebrahimi N, Chiang N, et al. Rve1 Protects From Local Inflammation and Osteoclast- Mediated Bone Destruction in Periodontitis. FASEB J Off Publ Fed Am Soc Exp Biol (2006) 20(2):401–3. doi: 10.1096/fj.05-4724fje

CrossRef Full Text | Google Scholar

61. Albuquerque-Souza E, Schulte F, Chen T, Hardt M, Hasturk H, Van Dyke TE, et al. Maresin-1 and Resolvin E1 Promote Regenerative Properties of Periodontal Ligament Stem Cells Under Inflammatory Conditions. Front Immunol (2020) 11:585530. doi: 10.3389/fimmu.2020.585530

PubMed Abstract | CrossRef Full Text | Google Scholar

62. El Kholy K, Freire M, Chen T, Van Dyke TE. Resolvin E1 Promotes Bone Preservation Under Inflammatory Conditions. Front Immunol (2018) 9:1300. doi: 10.3389/fimmu.2018.01300

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Gao L, Faibish D, Fredman G, Herrera BS, Chiang N, Serhan CN, et al. Resolvin E1 and Chemokine-Like Receptor 1 Mediate Bone Preservation. J Immunol (2013) 190(2):689–94. doi: 10.4049/jimmunol.1103688

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Zhu M, Van Dyke TE, Gyurko R. Resolvin E1 Regulates Osteoclast Fusion Via Dc-STAMP and Nfatc1. FASEB J Off Publ Fed Am Soc Exp Biol (2013) 27(8):3344–53. doi: 10.1096/fj.12-220228

CrossRef Full Text | Google Scholar

65. Fukuda H, Ikeda H, Muromoto R, Hirashima K, Ishimura K, Fujiwara K, et al. Synthesis of Resolvin E3, A Proresolving Lipid Mediator, and Its Deoxy Derivatives: Identification of 18-Deoxy-Resolvin E3 As A Potent Anti-Inflammatory Agent. J Org Chem (2020) 85(21):14190–200. doi: 10.1021/acs.joc.0c01701

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Isobe Y, Arita M, Iwamoto R, Urabe D, Todoroki H, Masuda K, et al. Stereochemical Assignment and Anti-Inflammatory Properties of the Omega-3 Lipid Mediator Resolvin E3. J Biochem (2013) 153(4):355–60. doi: 10.1093/jb/mvs151

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Isobe Y, Arita M, Matsueda S, Iwamoto R, Fujihara T, Nakanishi H, et al. Identification and Structure Determination of Novel Anti-Inflammatory Mediator Resolvin E3, 17,18-Dihydroxyeicosapentaenoic Acid. J Biol Chem (2012) 287(13):10525–34. doi: 10.1074/jbc.M112.340612

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Benabdoun HA, Kulbay M, Rondon EP, Vallieres F, Shi Q, Fernandes J, et al. In Vitro and In Vivo Assessment of the Proresolutive and Antiresorptive Actions of Resolvin D1: Relevance to Arthritis. Arthritis Res Ther (2019) 21(1):72. doi: 10.1186/s13075-019-1852-8

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Khaled M, Shibani NA, Labban N, Batarseh G, Song F, Ruby J, et al. Effects of Resolvin D1 on Cell Survival and Cytokine Expression of Human Gingival Fibroblasts. J Periodontol (2013) 84(12):1838–46. doi: 10.1902/jop.2013.120388

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Mustafa M, Zarrough A, Bolstad AI, Lygre H, Mustafa K, Hasturk H, et al. Resolvin D1 Protects Periodontal Ligament. Am J Physiol Cell Physiol (2013) 305(6):C673–9. doi: 10.1152/ajpcell.00242.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Du L, Li Y, Liu W. Maresin 1 Regulates Autophagy and Inflammation in Human Periodontal Ligament Cells Through Glycogen Synthase Kinase-3beta/Beta-Catenin Pathway Under Inflammatory Conditions. Arch Oral Biol (2018) 87:242–7. doi: 10.1016/j.archoralbio.2017.12.023

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Wang CW, Colas RA, Dalli J, Arnardottir HH, Nguyen D, Hasturk H, et al. Maresin 1 Biosynthesis and Proresolving Anti-Infective Functions With Human-Localized Aggressive Periodontitis Leukocytes. Infect Immun (2015) 84(3):658–65. doi: 10.1128/IAI.01131-15

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Chiang N, Fredman G, Backhed F, Oh SF, Vickery T, Schmidt BA, et al. Infection Regulates Pro-Resolving Mediators That Lower Antibiotic Requirements. Nature (2012) 484(7395):524–8. doi: 10.1038/nature11042

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Agans R, Gordon A, Kramer DL, Perez-Burillo S, Rufian-Henares JA, Paliy O. Dietary Fatty Acids Sustain the Growth of the Human Gut Microbiota. Appl Environ Microbiol (2018) 84(21) 84(21):e01525–18. doi: 10.1128/AEM.01525-18

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Schoeler M, Caesar R. Dietary Lipids, Gut Microbiota and Lipid Metabolism. Rev Endocr Metab Disord (2019) 20(4):461–72. doi: 10.1007/s11154-019-09512-0

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Herrera BS, Ohira T, Gao L, Omori K, Yang R, Zhu M, et al. An Endogenous Regulator of Inflammation, Resolvin E1, Modulates Osteoclast Differentiation and Bone Resorption. Br J Pharmacol (2008) 155(8):1214–23. doi: 10.1038/bjp.2008.367

PubMed Abstract | CrossRef Full Text | Google Scholar