Next Article in Journal
Effect of High-Intensity Interval Training and Intermittent Fasting on Body Composition and Physical Performance in Active Women
Next Article in Special Issue
Periodontal Conditions and Pathogens Associated with Pre-Eclampsia: A Scoping Review
Previous Article in Journal
Computational Intelligence-Based Model for Mortality Rate Prediction in COVID-19 Patients
Previous Article in Special Issue
Prevalence and Associated Factors of Self-Reported Gingival Bleeding: A Multicenter Study in France
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Analysis of Risk Factors for Periodontitis Focusing on the Saliva Microbiome and Polymorphism

1
Department of Preventive Dentistry, Okayama University Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, 2-5-1 Shikata-cho, Kita-ku, Okayama 700-8558, Japan
2
Department of Epidemiology for Community Health and Medicine, Kyoto Prefectural University of Medicine, 465 Kajii-cho, Kamigyo-ku, Kyoto 602-8566, Japan
3
Public Health Informatics Unit, Department of Integrated Health Sciences, Nagoya University Graduate School of Medicine, Nagoya 461-8673, Japan
4
Laboratory for Genotyping Development, RIKEN Center for Integrative Medical Sciences, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama City 230-0045, Japan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(12), 6430; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126430
Submission received: 20 May 2021 / Revised: 9 June 2021 / Accepted: 10 June 2021 / Published: 14 June 2021
(This article belongs to the Special Issue Oral Microbiota and Oral Health)

Abstract

:
Few studies have exhaustively assessed relationships among polymorphisms, the microbiome, and periodontitis. The objective of the present study was to assess associations simultaneously among polymorphisms, the microbiome, and periodontitis. We used propensity score matching with a 1:1 ratio to select subjects, and then 22 individuals (mean age ± standard deviation, 60.7 ± 9.9 years) were analyzed. After saliva collection, V3-4 regions of the 16S rRNA gene were sequenced to investigate microbiome composition, alpha diversity (Shannon index, Simpson index, Chao1, and abundance-based coverage estimator) and beta diversity using principal coordinate analysis (PCoA) based on weighted and unweighted UniFrac distances. A total of 51 single-nucleotide polymorphisms (SNPs) related to periodontitis were identified. The frequencies of SNPs were collected from Genome-Wide Association Study data. The PCoA of unweighted UniFrac distance showed a significant difference between periodontitis and control groups (p < 0.05). There were no significant differences in alpha diversity and PCoA of weighted UniFrac distance (p > 0.05). Two families (Lactobacillaceae and Desulfobulbaceae) and one species (Porphyromonas gingivalis) were observed only in the periodontitis group. No SNPs showed significant expression. These results suggest that periodontitis was related to the presence of P. gingivalis and the families Lactobacillaceae and Desulfobulbaceae but not SNPs.

1. Introduction

Periodontitis, which induces alveolar bone loss, is inflammation of periodontal tissue caused by bacterial infection. The World Health Organization reported that periodontitis is one of the main causes of tooth loss and can subsequently worsen an individual’s quality of life due to tooth loss [1]. Periodontitis affects some diseases, such as diabetes mellitus, kidney disease, premature birth, aspiration pneumonia, and arteriosclerosis [2]. Therefore, prevention of periodontitis is important to maintain both systemic and oral health.
Periodontitis is caused by a biofilm, which is a collection of any types of bacteria [3,4]. A biofilm of 1 mg includes over 1 billion bacteria. Some bacteria attach on teeth and make microcolonies. Those colonies gather and make biofilms. In biofilms, bacteria interact and obtain some benefits, including a broader habitat range, efficient metabolism, increased tolerance, and enhanced virulence [5]. The concept of microbiome analysis in periodontal science is to investigate the diversity of communities of bacteria.
It is important to control risk factors for periodontitis to prevent disease. Previous studies reported several risk factors, including smoking, diabetes mellitus, stress, and obesity [2]. Recently, studies reported the microbiome [6] and polymorphisms [7] as risk factors for periodontitis. Polymorphisms are the most common types of genetic variation among individuals [8]. They occur once in 500–1000 nucleotides. Most polymorphisms have no effect on periodontal health, but some polymorphisms are associated with risks of disease. Munz et al. showed that MTND1P5 (re16870060-G) and SHISA9 (rs729876-T) were associated with the risk of periodontitis [7].
Dentists have been struggling to prevent worsening periodontitis. Because the mechanism of worsening periodontitis is a host–parasite interaction, we have to pay attention to both the oral microbiome and host factors. However, few studies showed associations between the microbiome (parasite) and polymorphisms (host) at the same time. We hypothesized that there are associations between polymorphisms, the microbiome, and periodontitis. The objective of the present study was to assess associations simultaneously between polymorphisms, the microbiome, and periodontitis.

2. Materials and Methods

2.1. Study Population

The present study was a cross-sectional study including individuals who were enrolled in the Japan Multi-Institutional Collaborative Cohort Study (J-MICC study) second survey in the Kyoto area from 2013 to 2017. Participants underwent oral examinations, were measured for height and weight, provided saliva and blood samples, and answered self-reported questionnaires. The inclusion criterion was participants who provided saliva samples. The exclusion criteria were participants who had no teeth and could not be evaluated for the community periodontal index (CPI) and who had CPI = 1–3.

2.2. Ethical Procedures and Informed Consent

The study was approved by the Institutional Ethics Committee of Kyoto Prefectural University of Medicine (approval number: RBMR-E-36-8, in 2013) and was conducted in accordance with the principles of the Declaration of Helsinki. All participants provided written, informed consent before participation.

2.3. Oral Examinations

One general dentist checked dental status using a dental mirror and explorer under artificial light. Periodontal status was assessed using the CPI (World Health Organization, 4th edition). The CPI was measured for 10 teeth (maxilla: right first and second molars, right central incisor, left first and second molars; mandible: right first and second molars, left central incisor, left first and second molars). A CPI of 0 was defined as no periodontitis (control group), and a CPI of 4 was defined as periodontitis (periodontitis group).

2.4. Questionnaires

A self-administered questionnaire was used to determine whether participants had risk factors for periodontitis. Based on the Japanese Association of Periodontology Clinical Practice Guideline for the Periodontal Treatment, 2015 [2], we checked for risk factors such as age, sex, stress, smoking, and diabetes mellitus as follows:
  • Have you felt stress during the past year? (Severe/Moderate/Mild/None)
  • Are you a smoker? (Yes/Past/No)
  • Do you take at least one medication per week to lower blood sugar levels? (Yes/No)
We combined “Severe” and “Moderate” responses as positive awareness, and “Mild” and “None” responses as negative awareness to evaluate stress. Participants who took medicines to lower blood sugar levels were defined as having diabetes mellitus.

2.5. Salivary Microbiome Analysis

To investigate microbiome composition, saliva samples were collected. Saliva was collected by passive drool through a 1 inch straw into a vial. Approximately 1.0 mL of saliva was collected into a 10 mL centrifuge tube (Salivette®, Sarstedt AG & Co., Nümbrecht, Germany) from all examinees. Saliva samples in centrifuge tubes were stored at −20 °C until assayed. To prevent contamination, sterilized tubes were used for saliva collection.
DNA sequencing was performed based on the 16S metagenomics sequencing library preparation protocol [9] at the Oral Microbiome Center (Taniguchi Dental Clinic, Kagawa, Japan). First, DNA was extracted from saliva. Saliva was suspended using lysis buffer and crushed for 3 min. After centrifugation, DNA was extracted using GenFind 2.0 (Beckman Coulter Inc., Brea, CA, USA). During the first PCR, the V3 and V4 regions of the 16S rRNA gene were amplified using primers 341F (forward primer: 5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGNNNCCTACGGGNGGCWGCAG-3′) and 806R (reverse primer: 5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGNNNGACTACHVGGGTATCTAATCC-3′), using Kapa HiFi HotStart 2× ReadyMix DNA polymerase (Kapa Biosystems Ltd., London, UK). Cycle conditions were 95 °C (3 min), 28 cycles of 95 °C (30 s), 55 °C (30 s), and 72 °C (30 s), followed by a final extension of 72 °C (5 min). After extension, primers were removed using 50 μL AMPure XP. Thereafter, the second PCR was performed using primers including the 8 nt identifying index. Cycle conditions were 95 °C (3 min), 8 cycles of 95 °C (30 s), 55 °C (30 s), and 72 °C (30 s), followed by a final extension of 72 °C (5 min). Next-generation sequencing was performed using the obtained sequences and the MiSeq platform (MiSeq Reagent V3 600 cycles, Illumina, San Diego, CA, USA). Obtained reads showed that the mean Phred quality score 30 was 81.4%. Low-quality reads were discarded. After sequencing, UPARSE [10] was used to generate operational taxonomic units (OTUs), UCLUST [11] to divide the sequence into a cluster, and BLAST to search for homology (reference database, Greengenes 13.5).

2.6. Genotyping

The 14,539 study participants from the 12 areas of the J-MICC study were genotyped at the RIKEN Center for Integrative Medical Sciences using the Human OmniExpressExome-8 v1.2 BeadChip array (San Diego, CA, USA) [12]. The quality control of samples and SNPs was performed as in a previous study [8]. Briefly, data were subjected to quality control procedures, by which SNPs with a call rate of <0.98 or a Hardy–Weinberg equilibrium p-value <1 × 10−6 or a low minor allele frequency <0.01 were filtered out. Genotype imputation was performed using SHAPIT [13] and Minimac3 [14] software based on the 1000 Genomes Project cosmopolitan reference panel (phase 3). After genotype imputation, variants with an R2 < 0.3 were excluded, resulting in 12,617,547 variants. Of the SNPs that passed quality control, 51 SNPs related to periodontitis were identified [15,16,17,18,19,20,21,22].

2.7. Assessment of Other Factors

Data were collected from participants’ medical records for triglyceride, blood glucose, and HbA1c levels, and the body mass index (BMI) was calculated as weight (kg) divided by the square of height (m2).

2.8. Statistical Analysis

Participants were assigned to the periodontitis group or control group using propensity score matching with a 1:1 ratio [23] to adjust for confounders of periodontitis. Propensity scores were calculated using logistic regression models adjusted for age, sex, smoking status, stress, diabetes mellitus, and BMI, which were confounders. To check balanced covariates between the periodontitis and control groups, Fisher’s exact test and two-sample t-tests were performed after propensity score matching. Then, the variance inflation factor (VIF) was used to assess the multicollinearity of parameters used for propensity score matching. No multicollinearity of parameters was defined as VIF < 10. Model fitness using propensity score matching was assessed by the Hosmer–Lemeshow test.
Alpha (richness and evenness of bacterial taxa within a community) and beta (ecological distances between samples) diversities were analyzed to show the difference in microbiome composition between the periodontitis and control groups. Alpha diversities were assessed by the Shannon index, Simpson’s index, Chao1, and the abundance-based coverage estimator (ACE) [24]. Beta diversities were assessed by principal coordinate analysis (PCoA) based on weighted and unweighted UniFrac distances [25]. The weighted UniFrac distance reflects the difference in abundant lineages. Unweighted UniFrac distances show the difference in rare lineages [26]. The differences in microbial composition between the two groups were assessed by Adonis and ANOSIM, and the difference in dispersion was judged by Permdisp2 using “Calypso” (http://cgenome.net:8080/calypso-8.84/faces/uploadFiles.xhtml/; accessed on 6 March 2020).
Linear discriminant analysis effect size (LEfSe) showed differences in OTU richness between periodontitis and control groups using “Galaxy” (https://huttenhower.sph.harvard.edu/galaxy/; accessed on 6 March 2020). A cutoff value of 2.0 was used. For LDA scores >2.0, Mann–Whitney U tests were used to determine significant differences between the two groups. The difference in SNP expression between two groups was assessed using Fisher’s exact test.
Statistical analyses were conducted using SPSS version 22 (IBM, Tokyo, Japan) for Fisher’s exact test, the chi-squared test, the two-sample t-test, and the Mann–Whitney U test. All p-values <0.05 were considered significant. To adjust for multiple comparisons, the q-value was calculated using the Benjamini and Hochberg false discovery rate [27] based on the results of the p-value. All q-values <0.05 were considered significant.

2.9. Data Availability

Requests to access the dataset should be addressed to the J-MICC Study Central Office via the corresponding author of this paper.

3. Results

3.1. Participants’ Characteristics

Figure 1 shows the study flowchart. Overall, 385 of 3917 participants provided saliva samples in the J-MICC study. Then, 31 participants who met the inclusion criteria and with CPI = 0 or 4 were selected, and 22 participants (11 males and 11 females; mean age ± standard deviation (SD), 60.7 ± 9.9 years) were finally included after propensity score matching. No significant differences in risk factors for periodontitis (age, sex, smoking status, stress, diabetes mellitus, and BMI) or blood test items were observed between the periodontitis and control groups (Table 1). All parameters showed VIF < 10 (data not shown). The model showed good fit (Hosmer–Lemeshow test; p = 0.554).

3.2. Salivary Microbiome Analysis

A total of 2,314,918 reads were obtained from 22 saliva samples, and 1,993,572 quality-filtering reads (mean ± SD: 90,617 ± 19,112) were used for analysis. A total of 471 OTUs were obtained based on 97% sequence similarity, and 12 phyla, 23 classes, 36 orders, 58 families, 112 genera, and 349 species were identified. In alpha diversity, there were no indices that showed significance (Shannon index, p = 0.59; Simpson index, p = 0.24; Chao1, p = 0.18; ACE, p = 0.18) (Supplementary Figure S1).
Figure 2 shows the results of the PCoA of weighted and unweighted UniFrac distances. For unweighted UniFrac distances, there were significant differences between the periodontitis and control groups (ANOSIM, p = 0.004, R = 0.18; Adonis, p = 0.006, R2 = 0.105). However, there was no significant difference in the homogeneity of dispersions between the two groups (Permdisp2, p = 0.589). In contrast, for weighted UniFrac distances, there were no significant differences (ANOSIM, p = 0.136, R = 0.066; Adonis, p = 0.347, R2 = 0.049).
Figure 3 shows LEfSe at the species level. LDA scores >2.0 were observed for 1 phylum, 5 classes, 6 orders, 10 families, 14 genera, and 26 species in the periodontitis group, and 2 genera and 10 species in the control group. Significant differences were observed in one phylum (Synergistetes), five classes (Coriobacteriia, Bacteroidetes [C-1], Mollicutes, Deltaproteobacteria, Synergistia), five orders (Coriobacteriales, Bacteroidetes [O-1], Mycoplasmatales, Desulfobacterales, Synergistales), seven families (Atopobiaceae, Bacteroidetes [F-1], Lactobacillaceae, Mycoplasmataceae, Desulfobulbaceae, TM7 [F-2], Synergistaceae), two genera (Olsenella, Fretibacterium), and one species (Porphyromonas gingivalis) (Table 2).

3.3. Genotyping

No SNPs showed significantly different expressions between the two groups (q > 0.05, data not shown). Associations between periodontitis and SNPs are shown in Supplementary Table S1.

4. Discussion

The microbiome and polymorphisms for periodontitis were investigated in 22 individuals. Two families and one species were specifically confirmed in a periodontitis group. However, single-nucleotide polymorphisms were not related to periodontitis. To the best of our knowledge, the present study is the first to assess the associations among periodontitis, the microbiome, and genotyping. Beta diversity based on unweighted UniFrac distance showed a significant difference between the periodontitis and control groups, but the weighted UniFrac distance did not. Alpha diversity and SNPs were not associated with periodontitis. Unweighted UniFrac distance considers only species presence and is most efficient in detecting abundance change in rare lineages [26]. The present study showed that two families (Lactobacillaceae and Desulfobulbaceae) and one species (P. gingivalis) existed only in the periodontitis group (Table 2). Thus, the presence of the families Lactobacillaceae and Desulfobulbaceae and of P. gingivalis could be related to periodontitis. These results also suggested that the microbiome affects periodontitis more than SNPs do clinically.
P. gingivalis is the most abundant species in the subgingival microbiota of patients with periodontitis and is related to chronic periodontitis [28]. Damgaard et al. reported that the presence of P. gingivalis in saliva is significantly associated with periodontitis [28]. Guerra et al. reported that P. gingivalis in the salivary microbiome is a risk factor for periodontitis [6]. The present results support the existence of P. gingivalis as a risk factor for periodontitis.
The families Lactobacillaceae and Desulfobulbaceae were significantly more frequent in the periodontitis group than in the control group. Lactobacillaceae is isolated in the mouth, but its contribution to oral health has not been determined [29]. Desulfobulbaceae is a sulfate-reducing bacterium that is isolated from periodontal pockets [30]. Because few studies showed periodontitis associated with the families Lactobacillaceae and Desulfobulbaceae, further investigations are needed to clarify the associations.
Beta diversity in unweighted UniFrac distance was found to be significantly different between the periodontitis and control groups, whereas alpha and beta diversities in weighted UniFrac distance were not. Schulz et al. reported a significant difference in beta diversity, but not in alpha diversity, for periodontitis [31]. Acharya et al. also showed that alpha diversity was not significantly different between periodontitis and healthy gingiva [32]. The present findings support these findings. In contrast, Takeshita et al. showed that alpha diversity was significantly different between periodontitis and healthy gingiva [33]. One potential explanation for discrepancies in the results for diversity might be the effect of the circadian rhythm. Takayasu et al. reported that genera accounting for 79.3% of the oral microbiome change roughly every 24 h [34]. The circadian rhythm resulted in a greater number of bacteria at the phylum level (Firmicutes and Bacteroidetes) [34]. Therefore, the circadian rhythm might affect alpha and beta diversities in weighted UniFrac distance, which is based on the number of bacteria.
The results of the present study might have clinical relevance. The prevalence of periodontitis was significantly associated with the salivary microbiome, but not with SNPs. Some studies showed the association between SNPs and periodontitis [15,16,17,18,19,20,21,22], but few studies assessed the relationships among polymorphisms, the microbiome, and periodontitis simultaneously. The results showed that the microbiome, rather than SNPs, affected the prevalence of periodontitis. The effect of individual SNPs on periodontitis might be small. A previous study showed that many SNPs related to periodontitis had small odds ratios (<1.5) [35]. Therefore, clinicians should pay attention to microbiome composition to prevent periodontitis.
The present study had some strengths. First, associations among periodontitis, the microbiome, and SNPs were examined simultaneously. Second, covariates were reduced because propensity score matching was used [23]. Since periodontitis is a multifactorial disease, propensity score matching is useful to investigate the associations among periodontitis, the microbiome, and SNPs.
Some limitations must be considered when interpreting the results. First, the conditions of saliva collection were not controlled. The oral microbiome is affected by diurnal rhythm, oral health behavior, and food intake [36]. For generalization, it might be important to control the timing of saliva collection, food consumption, and oral health behaviors. The interval between saliva collection and eating was checked immediately before saliva collection, and it was not significantly different between the two groups (periodontitis group, 264.6 ± 79.8 min; control group, 241.2 ± 72.7 min; two-sample t-test, p = 0.585). However, when participants brushed their teeth before saliva collection was not known. Second, the causal association is unclear because the present study was a cross-sectional study. Therefore, it is unclear whether the present results are related to the onset or the progression of periodontitis. Third, there might be sampling bias. Saliva samples were provided from only a portion of participants in the J-MICC study (9.8%). Therefore, caution is needed before generalizing the present findings.

5. Conclusions

In conclusion, periodontitis was related to the presence of P. gingivalis and the families Lactobacillaceae and Desulfobulbaceae, but not to SNPs. The present results suggest that the oral microbiome, not polymorphism, is a risk factor for periodontitis, and that clinicians should pay more attention to microbiome composition than to host factors in the routine work of periodontal examination and diagnosis. Further investigations are needed to clarify the role of the families Lactobacillaceae and Desulfobulbaceae in periodontitis.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijerph18126430/s1, Figure S1: Species richness estimated by the Shannon index (a), Simpson’s index (b), Chao1 (c), and abundance-based coverage estimator (ACE) (d) based on operational taxonomic unit (OTU) abundance for the periodontitis and control groups. Table S1: Difference in genotypic frequencies between the periodontitis and control groups.

Author Contributions

Conceptualization, N.T. and D.E.; methodology, N.T. and D.E.; formal analysis, N.T.; investigation, N.T. and D.E.; resources, D.M., T.K., M.N., Y.M., and M.K.; data curation, N.T, D.M., and T.K.; writing—original draft preparation, N.T., D.E., and M.M.; writing—review and editing, N.T., D.E., D.M., T.K., M.N., Y.M., M.K., and M.M.; visualization, N.T., D.E., and M.M.; supervision, D.E. and M.M.; project administration, D.E. and M.M.; funding acquisition, N.T., D.E., D.M., T.K., and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the 8020 Promotion Foundation, grant number 19-5-14.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Ethics Committee of Kyoto Prefectural University of Medicine (protocol code ERB-E-36-9 and 6 June 2019).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Requests to access the dataset should be addressed to the J-MICC Study Central Office via the corresponding author of this paper.

Acknowledgments

This study was supported by Grants-in-Aid for Scientific Research for Priority Areas of Cancer (No. 17015018) and Innovative Areas (No. 221S0001) and by JSPS KAKENHI Grants (No. 16H06277 and 15H02524) from the Japanese Ministry of Education, Culture, Sports, Science and Technology. This work was also supported in part by a Grant-in-Aid from the Japan Society for the Promotion of Science (JSPS) KAKENHI Grant (Grant-in-Aid for Young Scientists) Number 15K19236 and 17K15840. This study was supported in part by funding for the BioBank Japan Project from the Japan Agency for Medical Research and development since April 2015, and the Ministry of Education, Culture, Sports, Science, and Technology from April 2003 to March 2015.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Petersen, P.E. The World Oral Health Report 2003: Continuous Improvement of Oral Health in the 21st Century—The Approach of the WHO Global Oral Health Programme. Commun. Dent. Oral Epidemiol. 2003, 31, 3–24. [Google Scholar] [CrossRef] [PubMed]
  2. JSP Clinical Practice Guideline for the Periodontal Treatment. 2015. Available online: file:///C:/Users/Toyama/Zotero/storage/9KLV5QTF/guideline_perio_plan2015.pdf (accessed on 8 June 2020).
  3. Vieira Colombo, A.P.; Magalhães, C.B.; Hartenbach, F.A.R.R.; Martins do Souto, R.; Maciel da Silva-Boghossian, C. Periodontal-Disease-Associated Biofilm: A Reservoir for Pathogens of Medical Importance. Microb. Pathog. 2016, 94, 27–34. [Google Scholar] [CrossRef] [PubMed]
  4. Kinane, D.F.; Stathopoulou, P.G.; Papapanou, P.N. Periodontal Diseases. Nat. Rev. Dis. Primers 2017, 3, 1–14. [Google Scholar] [CrossRef]
  5. Marsh, P.D.; Moter, A.; Devine, D.A. Dental Plaque Biofilms: Communities, Conflict and Control. Periodontology 2000 2011, 55, 16–35. [Google Scholar] [CrossRef]
  6. Guerra, F.; Mazur, M.; Ndokaj, A.; Corridore, D.; La Torre, G.; Polimeni, A.; Ottolenghi, L. Periodontitis and the Microbiome: A Systematic Review and Meta-Analysis. Minerva Stomatol. 2018, 67, 250–258. [Google Scholar] [CrossRef]
  7. Munz, M.; Richter, G.M.; Loos, B.G.; Jepsen, S.; Divaris, K.; Offenbacher, S.; Teumer, A.; Holtfreter, B.; Kocher, T.; Bruckmann, C.; et al. Meta-Analysis of Genome-Wide Association Studies of Aggressive and Chronic Periodontitis Identifies Two Novel Risk Loci. Eur. J. Hum. Genet. 2019, 27, 102–113. [Google Scholar] [CrossRef] [Green Version]
  8. Shastry, B.S. SNPs: Impact on Gene Function and Phenotype. In Single Nucleotide Polymorphisms: Methods and Protocols; Komar, A.A., Ed.; Methods in Molecular BiologyTM; Humana Press: Totowa, NJ, USA, 2009; pp. 3–22. ISBN 978-1-60327-411-1. [Google Scholar]
  9. 16S Metagenomic Sequencing Library Preparation. Available online: https://support.illumina.com/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf (accessed on 8 June 2020).
  10. Edgar, R.C. UPARSE: Highly Accurate OTU Sequences from Microbial Amplicon Reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  11. Edgar, R.C. Search and Clustering Orders of Magnitude Faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  12. Koyama, T.; Kuriyama, N.; Ozaki, E.; Matsui, D.; Watanabe, I.; Takeshita, W.; Iwai, K.; Watanabe, Y.; Nakatochi, M.; Shimanoe, C.; et al. Genetic Variants of RAMP2 and CLR Are Associated with Stroke. J. Atheroscler. Thromb. 2017, 24, 1267–1281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Delaneau, O.; Marchini, J.; Zagury, J.-F. A Linear Complexity Phasing Method for Thousands of Genomes. Nat. Methods 2012, 9, 179–181. [Google Scholar] [CrossRef]
  14. Das, S.; Forer, L.; Schönherr, S.; Sidore, C.; Locke, A.E.; Kwong, A.; Vrieze, S.I.; Chew, E.Y.; Levy, S.; McGue, M.; et al. Next-Generation Genotype Imputation Service and Methods. Nat. Genet. 2016, 48, 1284–1287. [Google Scholar] [CrossRef] [Green Version]
  15. Brodzikowska, A.; Górska, R.; Kowalski, J. Interleukin-1 Genotype in Periodontitis. Arch. Immunol. Ther. Exp. 2019, 67, 367–373. [Google Scholar] [CrossRef] [Green Version]
  16. Lin, W.; Xu, H.; Wu, Y.; Wang, J.; Yuan, Q. In Silico Genome-Wide Identification of M6A-Associated SNPs as Potential Functional Variants for Periodontitis. J. Cell. Physiol. 2020, 235, 900–908. [Google Scholar] [CrossRef] [PubMed]
  17. Deng, H.; Liu, F.; Pan, Y.; Jin, X.; Wang, H.; Cao, J. BsmI, TaqI, ApaI, and FokI Polymorphisms in the Vitamin D Receptor Gene and Periodontitis: A Meta-Analysis of 15 Studies Including 1338 Cases and 1302 Controls. J. Clin. Periodontol. 2011, 38, 199–207. [Google Scholar] [CrossRef] [PubMed]
  18. Schaefer, A.S.; Bochenek, G.; Manke, T.; Nothnagel, M.; Graetz, C.; Thien, A.; Jockel-Schneider, Y.; Harks, I.; Staufenbiel, I.; Wijmenga, C.; et al. Validation of Reported Genetic Risk Factors for Periodontitis in a Large-Scale Replication Study. J. Clin. Periodontol. 2013, 40, 563–572. [Google Scholar] [CrossRef]
  19. Taiete, T.; Casati, M.Z.; Stolf, C.S.; Corrêa, M.G.; Santamaria, M.P.; Andere, N.M.R.B.; Coletta, R.D.; Sallum, E.A.; Júnior, F.H.N.; Silvério, K.G.; et al. Validation of Reported GLT6D1 (Rs1537415), IL10 (Rs6667202), and ANRIL (Rs1333048) Single Nucleotide Polymorphisms for Aggressive Periodontitis in a Brazilian Population. J. Periodontol. 2019, 90, 44–51. [Google Scholar] [CrossRef] [Green Version]
  20. Leite, F.R.M.; Enevold, C.; Bendtzen, K.; Baelum, V.; López, R. Pattern Recognition Receptor Polymorphisms in Early Periodontitis. J. Periodontol. 2019, 90, 647–654. [Google Scholar] [CrossRef]
  21. Marchesan, J.T.; Jiao, Y.; Moss, K.; Divaris, K.; Seaman, W.; Webster-Cyriaque, J.; Zhang, S.; Yu, N.; Song, C.; Bencharit, S.; et al. Common Polymorphisms in IFI16 and AIM2 Genes Are Associated With Periodontal Disease. J. Periodontol. 2017, 88, 663–672. [Google Scholar] [CrossRef] [PubMed]
  22. Zupin, L.; Robino, A.; Navarra, C.O.; Pirastu, N.; Lenarda, R.D.; Gasparini, P.; Crovella, S.; Bevilacqua, L. LTF and DEFB1 Polymorphisms Are Associated with Susceptibility toward Chronic Periodontitis Development. Oral Dis. 2017, 23, 1001–1008. [Google Scholar] [CrossRef]
  23. Fernandez, R.; Tufanaru, C. Understanding Propensity Scores. Int. J. Evid. Based Healthc. 2017, 15, 142–143. [Google Scholar] [CrossRef] [Green Version]
  24. Kim, B.-R.; Shin, J.; Guevarra, R.B.; Lee, J.H.; Kim, D.W.; Seol, K.-H.; Lee, J.-H.; Kim, H.B.; Isaacson, R.E. Deciphering Diversity Indices for a Better Understanding of Microbial Communities. J. Microbiol. Biotechnol. 2017, 27, 2089–2093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Goodrich, J.K.; Di Rienzi, S.C.; Poole, A.C.; Koren, O.; Walters, W.A.; Caporaso, J.G.; Knight, R.; Ley, R.E. Conducting a Microbiome Study. Cell 2014, 158, 250–262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chen, J.; Bittinger, K.; Charlson, E.S.; Hoffmann, C.; Lewis, J.; Wu, G.D.; Collman, R.G.; Bushman, F.D.; Li, H. Associating Microbiome Composition with Environmental Covariates Using Generalized UniFrac Distances. Bioinformatics 2012, 28, 2106–2113. [Google Scholar] [CrossRef] [PubMed]
  27. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  28. Damgaard, C.; Danielsen, A.K.; Enevold, C.; Massarenti, L.; Nielsen, C.H.; Holmstrup, P.; Belstrøm, D. Porphyromonas Gingivalis in Saliva Associates with Chronic and Aggressive Periodontitis. J. Oral Microbiol. 2019, 11, 1653123. [Google Scholar] [CrossRef] [Green Version]
  29. Spinler, J.K. Human Microbiome, Lactobacillaceae in the. In Encyclopedia of Metagenomics; Nelson, K.E., Ed.; Springer: New York, NY, USA, 2014; pp. 1–8. ISBN 978-1-4614-6418-1. [Google Scholar]
  30. van der Hoeven, J.S.; van den Kieboom, C.W.; Schaeken, M.J. Sulfate-Reducing Bacteria in the Periodontal Pocket. Oral Microbiol. Immunol. 1995, 10, 288–290. [Google Scholar] [CrossRef] [PubMed]
  31. Schulz, S.; Porsch, M.; Grosse, I.; Hoffmann, K.; Schaller, H.-G.; Reichert, S. Comparison of the Oral Microbiome of Patients with Generalized Aggressive Periodontitis and Periodontitis-Free Subjects. Arch. Oral Biol. 2019, 99, 169–176. [Google Scholar] [CrossRef] [PubMed]
  32. Acharya, A.; Chen, T.; Chan, Y.; Watt, R.M.; Jin, L.; Mattheos, N. Species-Level Salivary Microbial Indicators of Well-Resolved Periodontitis: A Preliminary Investigation. Front. Cell. Infect. Microbiol. 2019, 9. [Google Scholar] [CrossRef] [Green Version]
  33. Takeshita, T.; Kageyama, S.; Furuta, M.; Tsuboi, H.; Takeuchi, K.; Shibata, Y.; Shimazaki, Y.; Akifusa, S.; Ninomiya, T.; Kiyohara, Y.; et al. Bacterial Diversity in Saliva and Oral Health-Related Conditions: The Hisayama Study. Sci. Rep. 2016, 6, 22164. [Google Scholar] [CrossRef] [Green Version]
  34. Takayasu, L.; Suda, W.; Takanashi, K.; Iioka, E.; Kurokawa, R.; Shindo, C.; Hattori, Y.; Yamashita, N.; Nishijima, S.; Oshima, K.; et al. Circadian Oscillations of Microbial and Functional Composition in the Human Salivary Microbiome. DNA Res. 2017, 24, 261–270. [Google Scholar] [CrossRef]
  35. Shimizu, S.; Nagasawa, T.; Furuichi, Y. A genome-wide association study of periodontitis. Nihon Shishubyo Gakkai Kaishi (J. Jpn. Soc. Periodontol.) 2019, 61, 127–135. [Google Scholar] [CrossRef]
  36. Cornejo Ulloa, P.; van der Veen, M.H.; Krom, B.P. Review: Modulation of the Oral Microbiome by the Host to Promote Ecological Balance. Odontology 2019, 107, 437–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Study flowchart. CPI, community periodontal index.
Figure 1. Study flowchart. CPI, community periodontal index.
Ijerph 18 06430 g001
Figure 2. Principal coordinate analysis (PCoA) to assess the difference in beta diversity between the periodontitis and control groups. Plots were based on weighted UniFrac distance (a), unweighted UniFrac distance (b), and relative abundances of bacterial operational taxonomic units.
Figure 2. Principal coordinate analysis (PCoA) to assess the difference in beta diversity between the periodontitis and control groups. Plots were based on weighted UniFrac distance (a), unweighted UniFrac distance (b), and relative abundances of bacterial operational taxonomic units.
Ijerph 18 06430 g002
Figure 3. Graphics of linear discriminant analysis (LDA) effect size (LEfSe) for the periodontitis and control groups. Horizontal bars show the effect size for each taxon. The length of the bar indicates the log10 transformed LDA score, indicated by vertical dotted lines. The threshold on the logarithmic LDA score for discriminative features was set to 2.0. The taxon of bacteria with significant change in relative abundance (n = 74, linear discriminant analysis, p < 0.05) is written alongside the horizontal lines. The name of the taxon level was abbreviated as follows: p, phylum; c, class; o, order; f, family; g, genus; and s, species.
Figure 3. Graphics of linear discriminant analysis (LDA) effect size (LEfSe) for the periodontitis and control groups. Horizontal bars show the effect size for each taxon. The length of the bar indicates the log10 transformed LDA score, indicated by vertical dotted lines. The threshold on the logarithmic LDA score for discriminative features was set to 2.0. The taxon of bacteria with significant change in relative abundance (n = 74, linear discriminant analysis, p < 0.05) is written alongside the horizontal lines. The name of the taxon level was abbreviated as follows: p, phylum; c, class; o, order; f, family; g, genus; and s, species.
Ijerph 18 06430 g003
Table 1. Differences in parameters between the periodontal disease and control groups.
Table 1. Differences in parameters between the periodontal disease and control groups.
Periodontal Disease GroupControl Groupp
N = 11N = 11
GenderMale6 (54.5) a5 (45.5)1.000 b
Female5 (45.5)6 (54.5)
Passive stressNegative5 (45.5)2 (18.2)0.361
Positive6 (54.5)9 (81.8)
Smoking statusNo3 (27.3)7 (63.6)0.061 c
Past4 (36.4)4 (36.4)
Current4 (36.4)0 (0)
DiabetesNo9 (81.8)11 (100)0.333
Past1 (9.1)0 (0)
Current1 (9.1)0 (0)
Age 64.3 ± 7.6 d57.1 ± 11.00.091 e
Body mass index 23.9 ± 3.321.3 ± 2.70.051
Triglyceride (mg/dL) 102.5 ± 34.789.2 ± 48.80.470
HbA1c (%) 5.8 ± 0.35.0 ± 0.30.065
Blood glucose (mg/dL) 100.8 ± 19.392.7 ± 13.00.263
a N (%), b Fisher’s exact test, c chi-square test, d mean ± SD, e two sample t test. HbA1c, hemoglobin A1c.
Table 2. Differences in operational taxonomic units (OTUs) between the periodontitis and control groups.
Table 2. Differences in operational taxonomic units (OTUs) between the periodontitis and control groups.
OTUPeriodontitisControl paqb
SpeciesMedian25th Percentiles75th PercentilesMedian25th Percentiles75th Percentiles
Actinomyces cardiffensis1.31 × 10−402.48 × 10−40000.0470.074
Rothia aeria1.20 × 10−39.45 × 10−51.15 × 10−28.36 × 10−32.28 × 10−32.16 × 10−20.0470.071
Propionibacterium propionicum002.74 × 10−54.58 × 10−52.51 × 10−58.32 × 10−50.0230.075
Olsenella sp. oral taxon 8074.87 × 10−501.29 × 10−40000.0650.087
Olsenella uli5.47 × 10−501.31 × 10−4002.13 × 10−50.0400.076
Slackia exigua2.25 × 10−42.75 × 10−56.58 × 10−43.02 × 10−505.12 × 10−50.0400.072
Bacteroidaceae [G-1] sp. oral taxon 2723.56 × 10−42.19 × 10−57.54 × 10−42.48 × 10−501.07 × 10−40.0470.068
Porphyromonas endodontalis7.06 × 10−31.14 × 10−32.50 × 10−21.12 × 10−303.69 × 10−30.0340.068
Porphyromonas gingivalis6.61 × 10−44.96 × 10−42.79 × 10−3000<0.0010.007
Tannerella forsythia1.73 × 10−31.05 × 10−38.62 × 10−34.29 × 10−44.22 × 10−51.29 × 10−30.0280.072
Prevotella marshii2.36 × 10−508.26 × 10−50000.0280.067
Prevotella nanceiensis8.27 × 10−41.26 × 10−46.66 × 10−38.12 × 10−32.73 × 10−31.60 × 10−20.0230.069
Prevotella sp. oral taxon 472002.75 × 10−51.02 × 10−46.45 × 10−51.49 × 10−40.0070.084
Prevotella sp. oral taxon 5261.31 × 10−406.76 × 10−40000.0190.086
Lactobacillus crispatus001.08 × 10−40000.1510.160
Lactobacillus ultunensis5.90 × 10−501.44 × 10−40000.0650.084
Streptococcus oligofermentans8.05 × 10−52.53 × 10−59.74 × 10−52.13 × 10−502.77 × 10−50.0050.090
Streptococcus sp. oral taxon 0663.28 × 10−21.34 × 10−25.78 × 10−21.28 × 10−29.79 × 10−32.09 × 10−20.0130.078
Catonella sp. oral taxon 1642.74 × 10−501.31 × 10−40000.0760.094
Oribacterium parvum0002.29 × 10−504.30 × 10−50.1010.121
Peptostreptococcaceae [XI][G-1] sp. oral taxon 383007.18 × 10−50000.1510.151
Mollicutes [G-2] sp. oral taxon 9061.61 × 10−501.10 × 10−40000.1010.117
Mycoplasma faucium6.68 × 10−401.26 × 10−3002.15 × 10−50.0130.067
Anaeroglobus geminatus3.59 × 10−503.14 × 10−40000.0190.076
Leptotrichia sp. oral taxon 2214.87 × 10−507.67 × 10−47.46 × 10−41.37 × 10−45.99 × 10−30.0400.069
Ottowia sp. oral taxon 894002.74 × 10−51.07 × 10−42.56 × 10−51.60 × 10−40.0470.065
Kingella sp. oral taxon 459000001.51 × 10−40.1510.151
Neisseria subflava1.89 × 10−41.61 × 10−53.02 × 10−35.72 × 10−22.60 × 10−27.43 × 10−20.0080.058
Desulfobulbus sp. oral taxon 0415.51 × 10−502.51 × 10−40000.0280.063
Aggregatibacter sp. oral taxon 5122.36 × 10−501.45 × 10−41.54 × 10−42.51 × 10−51.26 × 10−30.0400.065
TM7 [G-2] sp. oral taxon 350002.20 × 10−40000.1330.150
TM7 [G-3] sp. oral taxon 3511.04 × 10−22.38 × 10−31.90 × 10−22.18 × 10−21.36 × 10−23.79 × 10−20.0230.064
TM7 [G-5] sp. oral taxon 3563.58 × 10−41.77 × 10−41.13 × 10−31.26 × 10−403.05 × 10−40.0190.068
Treponema lecithinolyticum001.09 × 10−40000.1330.145
Fretibacterium fastidiosum8.75 × 10−54.87 × 10−51.89 × 10−41.79 × 10−506.04 × 10−50.0280.059
Fretibacterium sp. oral taxon 3602.92 × 10−45.90 × 10−57.66 × 10−42.56 × 10−509.16 × 10−50.0070.063
Genus
Bifidobacterium003.03 × 10−40000.1330.142
Olsenella1.75 × 10−42.75 × 10−53.29 × 10−4003.57 × 10−50.0020.032
Slackia2.25 × 10−42.75 × 10−56.58 × 10−43.02 × 10−505.12 × 10−50.0400.058
Bacteroidetes [G-3]2.87 × 10−44.38 × 10−51.02 × 10−3001.83 × 10−40.0340.060
Bacteroidaceae [G-3]3.56 × 10−42.19 × 10−57.54 × 10−42.48 × 10−501.07 × 10−40.0470.063
Lactobacillus6.56 × 10−502.51 × 10−40000.0340.054
Peptostreptococcaceae [XI][G-6]1.64 × 10−41.09 × 10−41.47 × 10−34.26 × 10−501.25 × 10−40.0100.053
Mollicutes [G-2]1.61 × 10−501.10 × 10−40000.1010.115
Mycoplasma1.05 × 10−31.31 × 10−41.51 × 10−31.24 × 10−42.77 × 10−55.63 × 10−40.0230.061
Anaeroglobus3.59 × 10−503.14 × 10−40000.0190.076
Ottowia002.74 × 10−51.07 × 10−42.56 × 10−51.60 × 10−40.0470.058
Desulfobulbus5.51 × 10−502.51 × 10−40000.0280.056
TM7 [G-2]002.20 × 10−40000.1330.133
TM7 [G-3]1.04 × 10−22.38 × 10−31.90 × 10−22.18 × 10−21.36 × 10−23.79 × 10−20.0230.053
TM7 [G-5]3.58 × 10−41.77 × 10−41.13 × 10−31.26 × 10−403.05 × 10−40.0190.061
Fretibacterium4.14 × 10−41.97 × 10−41.65 × 10−35.02 × 10−51.79 × 10−51.81 × 10−40.0040.032
Family
Atopobiaceae1.75 × 10−42.75 × 10−53.29 × 10−4003.57 × 10−50.0020.020
Coriobacteriaceae2.46 × 10−42.75 × 10−56.58 × 10−43.02 × 10−51.79 × 10−55.12 × 10−50.0470.052
Bacteroidetes [F-1]8.77 × 10−41.61 × 10−42.55 × 10−31.07 × 10−402.77 × 10−40.0340.049
Bacteroidaceae7.09 × 10−41.38 × 10−41.21 × 10−39.14 × 10−501.83 × 10−40.0190.063
Lactobacillaceae6.56 × 10−502.51 × 10−40000.0340.043
Mollicutes [F-1]1.61 × 10−501.10 × 10−40000.1010.101
Mycoplasmataceae1.05 × 10−31.31 × 10−41.51 × 10−31.24 × 10−42.77 × 10−55.63 × 10−40.0230.046
Desulfobulbaceae5.51 × 10−502.51 × 10−40000.0280.047
TM7 [F-2]3.58 × 10−41.77 × 10−41.13 × 10−31.26 × 10−403.05 × 10−40.0190.048
Synergistaceae4.14 × 10−41.97 × 10−41.65 × 10−35.02 × 10−51.79 × 10−51.81 × 10−40.0040.020
Order
Coriobacteriales4.18 × 10−41.65 × 10−47.51 × 10−45.02 × 10−52.15 × 10−57.45 × 10−50.0050.015
Bacteroidetes [O-1]8.77 × 10−41.61 × 10−42.55 × 10−31.07 × 10−402.77 × 10−40.0340.041
Mollicutes [O-2]1.61 × 10−501.10 × 10−40000.1010.101
Mycoplasmatales1.05 × 10−31.31 × 10−41.51 × 10−31.24 × 10−42.77 × 10−55.63 × 10−40.0230.046
Desulfobacterales5.51 × 10−502.51 × 10−40000.0280.042
Synergistales4.14 × 10−41.97 × 10−42.23 × 10−35.02 × 10−51.79 × 10−51.81 × 10−40.0020.012
Class
Coriobacteriia9.98 × 10−41.97 × 10−45.06 × 10−35.02 × 10−505.02 × 10−50.0050.013
Bacteroidetes [C-1]8.77 × 10−41.61 × 10−42.55 × 10−33.02 × 10−501.81 × 10−40.0340.034
Mollicutes1.17 × 10−31.65 × 10−41.51 × 10−34.96 × 10−501.74 × 10−40.0160.020
Deltaproteobacteria8.26 × 10−502.51 × 10−45.54 × 10−42.13 × 10−51.13 × 10−30.0070.012
Synergistia4.14 × 10−41.97 × 10−42.23 × 10−39.14 × 10−43.05 × 10−59.90 × 10−40.0020.010
Phylum
Synergistetes4.14 × 10−41.97 × 10−42.23 × 10−35.02 × 10−51.79 × 10−51.81 × 10−40.002
a Mann–Whitney U test, b adjusted p-values by Benjamini and Hochberg’s false discovery rate (FDR).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Toyama, N.; Ekuni, D.; Matsui, D.; Koyama, T.; Nakatochi, M.; Momozawa, Y.; Kubo, M.; Morita, M. Comprehensive Analysis of Risk Factors for Periodontitis Focusing on the Saliva Microbiome and Polymorphism. Int. J. Environ. Res. Public Health 2021, 18, 6430. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126430

AMA Style

Toyama N, Ekuni D, Matsui D, Koyama T, Nakatochi M, Momozawa Y, Kubo M, Morita M. Comprehensive Analysis of Risk Factors for Periodontitis Focusing on the Saliva Microbiome and Polymorphism. International Journal of Environmental Research and Public Health. 2021; 18(12):6430. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126430

Chicago/Turabian Style

Toyama, Naoki, Daisuke Ekuni, Daisuke Matsui, Teruhide Koyama, Masahiro Nakatochi, Yukihide Momozawa, Michiaki Kubo, and Manabu Morita. 2021. "Comprehensive Analysis of Risk Factors for Periodontitis Focusing on the Saliva Microbiome and Polymorphism" International Journal of Environmental Research and Public Health 18, no. 12: 6430. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18126430

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop