- 1West China Medical School, West China Hospital, Sichuan University, Chengdu, Sichuan, China
- 2Department of Endocrinology and Andrology Laboratory, West China Hospital, Sichuan University, Sichuan, Chengdu, China
Introduction: The dysbiosis of the oral microbiome is associated with the progression of various systemic diseases, including diabetes. However, the precise causal relationships remain elusive. This study aims to investigate the potential causal associations between oral microbiome and type 2 diabetes (T2D) using Mendelian randomization (MR) analyses.
Methods: We conducted bidirectional two-sample MR analyses to investigate the impact of oral microbiome from saliva and the tongue T2D. This analysis was based on metagenome-genome-wide association studies (mgGWAS) summary statistics of the oral microbiome and a large meta-analysis of GWAS of T2D in East Asian populations. Additionally, we utilized the T2D GWAS summary statistics from the Biobank Japan (BBJ) project for replication. The MR methods employed included Wald ratio, inverse variance weighting (IVW), weighted median, MR-Egger, contamination mixture (ConMix), and robust adjusted profile score (RAPS).
Results: Our MR analyses revealed genetic associations between specific bacterial species in the oral microbiome of saliva and tongue with T2D in East Asian populations. The MR results indicated that nine genera were shared by both saliva and tongue. Among these, the genera Aggregatibacter, Pauljensenia, and Prevotella were identified as risk factors for T2D. Conversely, the genera Granulicatella and Haemophilus D were found to be protective elements against T2D. However, different species within the genera Catonella, Lachnoanaerobaculum, Streptococcus, and Saccharimonadaceae TM7x exhibited multifaceted influences; some species were positively correlated with the risk of developing T2D, while others were negatively correlated.
Discussion: This study utilized genetic variation tools to confirm the causal effect of specific oral microbiomes on T2D in East Asian populations. These findings provide valuable insights for the treatment and early screening of T2D, potentially informing more targeted and effective therapeutic strategies.
1 Introductions
Diabetes mellitus encompasses a range of metabolic disorders related to carbohydrate metabolism, marked by inadequate glucose utilization as an energy source and excessive glucose production due to abnormal gluconeogenesis and/or glycogenolysis, resulting in hyperglycemia (1). Approximately 537 million adults worldwide suffer from diabetes, with over 90% of these cases being type 2 diabetes (T2D), and this number is projected to increase to 783 million by 2045 (2). Recently, large population-based studies have shown that periodontal disease adversely affects glycemic control, diabetes complications, and the progression of T2D (3, 4). Treatment of oral diseases has been demonstrated to improve glycemic control and reduce HbA1c levels (5).
The oral microbiome is the second largest microbial community in the human body, following the gut microbiome. It comprises over 700 species of bacteria, fungi, viruses, and protozoa (6). These microorganisms colonize the teeth, prosthodontic surfaces, mucosal surfaces, and are abundant in saliva (7). The oral microbiome is integral not only to oral diseases such as caries and periodontitis but also to systemic health (8, 9). This influence is mediated through complex interactions with the host immune system, the gut microbiome, and various small molecule metabolites (10), and these interactions can impact systemic health by inhibiting pathogens, modulating the immune response, and affecting nutrient absorption and metabolism (11). Previous studies have indicated that dysbiosis of the oral microbiome is associated with the progression of various systemic diseases, including diabetes (12). For instance, certain periodontal pathogens, such as Porphyromonas gingivalis and Aggregatibacter actinomycetemcomitans have been directly linked to glycemic control and the risk of developing diabetes (13, 14).
However, compared to the extensive research on the gut microbiome, studies on the oral microbiome were relatively limited and often involved small sample sizes. Traditional observational studies may also suffer from inadequate control of confounding variables and the possibility of reverse causation. Therefore, we employed Mendelian randomization (MR) as an epidemiological tool (15). MR leverages naturally randomized genetic variants at conception as a form of natural experiment to uncover causal relationships between exposures and outcomes, thereby minimizing the potential for reverse causation and confounding biases (16, 17). In this study, we conducted bidirectional two-sample MR analyses to investigate the impact of oral microbiome from saliva and the tongue T2D using single nucleotide polymorphisms (SNPs) as instrumental variables.
2 Materials and methods
2.1 Study design and population
The research workflow is presented in Figure 1. The summary statistics for the metagenome-genome-wide association studies (mgGWAS) of the oral microbiome in East Asian populations were obtained from the research by Liu et al. (18), This study comprised 2017 tongue samples and 1915 salivary samples, derived from a cohort of 2984 healthy Chinese individuals with provided high-depth whole genome sequencing data. In this study, the lowest taxonomic level is species. The composition of the microbiome used in the study was determined by aligning it against 56,213 metagenome-assembled genomes (MAGs), which were organized into 3,589 species-level clusters (SGBs). The criteria for sample inclusion were as follows: a) a variant calling rate exceeding 98%; b) an average sequencing depth greater than 20×; c) no population stratification as evidenced by principal component analysis (PCA); d) the exclusion of related individuals based on pairwise identity by descent calculations. Additionally, a stringent inclusion threshold was applied for variants, requiring a mean depth greater than 8×, Hardy-Weinberg equilibrium (HWE) P-value greater than 10−5, and a genotype calling rate exceeding 98%.
Summary statistics of T2D were obtained from meta-analyses of GWAS conducted by Spracklen et al (19), which encompassed 77418 T2D cases and 356122 controls (effective sample size, Neff=211793) of East Asian individuals from 23 GWAS, For each study involved in this meta-analysis, variants were excluded according to the following criteria: a) mismatched chromosomal positions or alleles compared to the reference panel; b) ambiguous alleles (AT/CG) with a minor allele frequency (MAF) greater than 40% in the reference panel; or c) discrepancies in allele frequencies exceeding 20% when compared to East Asian-specific allele frequencies. This study was designated as the discovery cohort. Furthermore, we utilized the T2D GWAS summary statistics from the Biobank Japan (BBJ) project as a replication cohort (20), comprising 45383 T2D cases and 132032 controls of East Asian individuals. There is no sample overlap between the exposure and the outcome. The details of these datasets are shown in Supplementary Table 1.
2.2 Instrumental variables selection
In MR analysis, we employed the PLINK clumping (21) function to identify a sufficient number of independent instrumental variables (IV), which reduces linkage disequilibrium (LD) among associated genetic variants, avoids multicollinearity due to LD, and mitigates biases caused by weak instruments. Initially, we selected instrumental variables based on a genome-wide significance threshold of P <5×10-8 to test our hypothesis. However, due to the limited number of instruments available at this threshold, we adopted a more lenient threshold to ensure the robustness of our analysis: p1 = 5×10-6, p2 = 1×10-5, kb=10000Kb, and r2 = 0.001, to identify top loci. SNPs with a minor allele frequency (MAF) < 0.01 are generally considered rare SNPs, which have a limited impact on traits. Hence, only SNPs with MAF ≥ 0.01 were retained. Additionally, we applied Steiger filtering (22) to the instrumental variables and excluded instruments with F-statistics (F= (beta/se) 2) < 10 to mitigate the impact of weak instrumental variables (23).
2.3 Bidirectional MR analysis
We utilized six MR analysis methods to investigate the impact of oral microbiome from saliva and tongue dorsum on T2D, including Wald ratio, inverse variance weighting (IVW) (17), weighted median (WM) (24), MR-Egger (25), Contamination mixture (ConMix) (26), robust adjusted profile score (RAPS) (27). For MR analyses with only one instrumental variable, we employed the Wald ratio as the primary analysis method. Additionally, as the RAPS method could produce consistent results in the presence of weak and pleiotropic SNPs, it was utilized for supplementary validation in MR analyses with only one instrumental variable. IVW was chosen as the primary analysis method for MR analyses with multiple instrumental variables, given its robustness. The ConMix method explicitly models multiple potential causal estimates and infers various causal mechanisms linked to the same risk factor, each impacting the outcome to different extents. MR-Egger regression offers estimates corrected for pleiotropy. The WM estimator, which calculates the median of the weighted estimates, provides a consistent effect even when up to half of the instrumental variables are pleiotropic. We applied FDR for multiple testing correction, with FDR<0.05 indicating significance and P < 0.05 suggestive significance.
We conducted various heterogeneity and pleiotropy analyses to evaluate the robustness of our results against potential violations of multiple MR assumptions. a) Heterogeneity was assessed through the Cochran Q test of IVW and MR-Egger methods; b) Horizontal pleiotropy was evaluated using MR-Egger’s intercept. The same approach was applied to reverse MR analyses to mitigate spurious results arising from reverse causation. Additionally, a stricter threshold was applied for instrumental variables with T2D as the exposure (p1 = 5×10-8, p2 = 1×10-5, kb=10000Kb, r2 = 0.001) to enhance result reliability.
All the analyses were conducted using R software 4.2.0. The IVW, MR–Egger, WM, MR-RAPS, and ConMix methods were performed using the “TwoSampleMR” package.
3 Results
3.1 Causal effects of oral microbiome in the saliva on the development of T2D
All genetic instruments used in the MR analyses passed the Steiger test. Additionally, the F-statistics for these instruments were greater than 10, indicating strong instrument validity. (Supplementary Tables 2, 3). In the discovery cohort of T2D, a total of 89 bacterial species in saliva (46 genera, 28 families, 20 orders, 11 classes, and 8 phyla) had statistically significant relationships (P < 0.05) with T2D under either the IVW or Wald ratio MR methods. Among these, 48 bacterial species were confirmed to be significant under the RAP method, with consistent effect directions (Figure 2A; Supplementary Table 4). In the replication cohort of T2D, a total of 50 bacterial species in saliva (29 genera, 19 families, 17 orders, 11 classes, and 8 phyla) had statistically significant relationships (P < 0.05) with T2D under either the IVW or Wald ratio MR analysis methods. Among these, 26 bacterial species were confirmed to be significant under the RAP MR analysis method, with consistent effect directions (Figure 2B; Supplementary Table 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1452999/fendo-15-1452999-HTML-r1/image_m/fendo-15-1452999-g002.jpg)
Figure 2. Significant mendelian randomization (MR) results of microbiome in saliva and tongue coating for type 2 diabetes (T2D) in discovery and replication Cohorts. (A) Oral microbiome in saliva showing significant MR results in the discovery cohort. (B) Oral microbiome in saliva showing significant MR results in the replication cohort. (C) Oral microbiome on the tongue showing significant MR results in the discovery cohort. (D) Oral microbiome on the tongue showing significant MR results in the replication cohort. * indicates MR P-value < 0.05.
In both the discovery and replication cohorts, 23 bacterial species were consistently significant under either the IVW or Wald ratio MR methods. Of these, 12 species were positively correlated with the risk of developing T2D (OR > 1), while 11 species were negatively correlated (OR < 1) (Supplementary Table 6). We further examined whether the instrumental variables corresponding to these 23 bacteria were associated with other confounding factors. We found that rs10421891, the instrumental variable for Prevotella unclassified metagenome species (uMGS) 2017, has been reported to exhibit genome-wide significance about heart failure and left ventricular systolic function phenotypes. No associations with other confounding factors were identified for the remaining instrumental variables (Supplementary Table 7).
Notably, in the discovery and replication cohorts, several genera showed significant relationships with T2D risk in at least one MR analysis method other than IVW and Wald ratio. Except for the MR Egger analysis of Haemophilus D, which indicated some heterogeneity (QEgger=4.239, Pheterogeneity=3.95×10-2), the remaining sensitivity analyses of the MR showed no significant heterogeneity or pleiotropy (Figure 3; Supplementary Table 6). These included Pauljensenia, Streptococcus, Lachnoanaerobaculum, Saccharimonadaceae TM7x, Saccharimonadaceae UBA2866, and Saccharimonadaceae uMGS 1898, which were positively correlated with T2D risk, as well as Alloprevotella, Granulicatella, Streptococcus, Saccharimonadaceae uMGS 1251, and Haemophilus D, which were negatively correlated with T2D risk. Reverse MR analysis did not reveal any causal relationships (all P > 0.05, Supplementary Table 8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1452999/fendo-15-1452999-HTML-r1/image_m/fendo-15-1452999-g003.jpg)
Figure 3. Significant MR results of oral microbiome in saliva for T2D in both discovery and replication Cohorts. † indicates a P-value < 0.05 in at least one MR method other than IVW or Wald ratio.
3.2 Causal effects of oral microbiome in the tongue on the development of T2D
All genetic instruments employed for MR analyses successfully passed the Steiger test. Furthermore, the F-statistics of these genetic instruments exceeded 10, demonstrating robust instrument strength (Supplementary Tables 9, 10). In the discovery cohort of T2D, we identified 114 bacterial species in saliva (55 genera, 30 families, 21 orders, 11 classes, and 8 phyla) that exhibited statistically significant associations (P < 0.05) with T2D when analyzed using either the IVW or Wald ratio MR methods. Of these, 66 bacterial species were further validated as significant under the RAP method, maintaining consistent effect directions (Figure 2C; Supplementary Table 11). Similarly, in the replication cohort of T2D, 73 bacterial species in saliva (spanning 40 genera, 21 families, 18 orders, 11 classes, and 8 phyla) showed statistically significant associations (P < 0.05) with T2D using IVW or Wald Ratio. Out of these, 44 bacterial species were confirmed as significant via the RAP method, with effect directions remaining consistent (Figure 2D; Supplementary Table 12).
In both the discovery and replication cohorts, 33 bacterial species were consistently significant under either the IVW or Wald ratio MR methods. Within this group, 13 species were positively associated with the risk of developing T2D (OR > 1), whereas 20 species were negatively associated (OR < 1) (Supplementary Table 13). To assess potential associations with other confounding factors, we analyzed the instrumental variables for these 33 bacteria. Our findings indicated that rs4566929, the instrumental variable representing Streptococcus uMGS 2424, has shown genome-wide significance for the body weight phenotype (Supplementary Table 7). For the other instrumental variables, no significant associations with additional confounding factors were detected. Notably, six species remained significant in the MR analysis after FDR multiple corrections (Figure 4; Supplementary Table 13). These include Actinomycetaceae F0332 (genera), Streptococcus uMGS 988, Streptococcus uMGS 2424, and Prevotella uMGS 1437, all of which are positively associated with T2D risk (OR >1), and Catonella uMGS 2059 and Lachnoanaerobaculum sp000287675 MGS 1966, both of which are negatively associated with T2D risk (OR < 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1452999/fendo-15-1452999-HTML-r1/image_m/fendo-15-1452999-g004.jpg)
Figure 4. Significant MR results of oral microbiome on the tongue for T2D in both discovery and replication cohorts. † indicates a P-value < 0.05 in at least one MR method other than IVW or Wald ratio. * indicates that the FDR of the MR is less than 0.05.
Importantly, several genera demonstrated significant associations with T2D risk in at least one MR method other than IVW and Wald ratio in both cohorts. These included Actinomycetaceae F0332, Bacteroidaceae F0040, Prevotella, CAG-917 uMGS 1705, Aggregatibacter, Saccharimonadaceae TM7x, and Streptococcus, which were positively correlated with T2D risk. Conversely, Catonella, Lachnoanaerobaculum, Neisseria, Haemophilus D, Saccharimonadaceae TM7x, and Veillonellaceae F0422 were negatively correlated with T2D risk (Figure 4; Supplementary Table 13). Reverse MR analysis did not reveal any causal relationships (all P > 0.05). The sensitivity analysis of the MR indicated no significant heterogeneity or pleiotropy (Supplementary Table 14).
Notably, nine genera were shared by both the saliva and tongue, including Aggregatibacter, Catonella, Granulicatella, Haemophilus D, Lachnoanaerobaculum, Pauljensenia, Prevotella, Streptococcus and Saccharimonadaceae TM7x (Supplementary Table 15).
4 Discussion
Oral microorganisms are integral to the development and progression of oral diseases such as dental caries and periodontitis through mechanisms of pathogen inhibition and immune modulation (1, 2). Additionally, the oral microbiome contributes to overall systemic health by regulating immune responses, aiding nutrient absorption, and metabolism, and has been closely linked to systemic diseases such as T2D (28, 29). Despite this, the precise mechanisms and genetic causal relationships remain poorly understood. This study employs MR analysis to elucidate the genetic connections between the oral microbiome and T2D.
Aemaimanan et al. indicated that poor glycemic control is associated with increasing cell numbers of the red complex bacteria (Porphyromonas gingivalis, Treponema denticola, and Tannerella forsythia) within the subgingival biofilm (30). Li et al. also identified that periodontal pathogens, including Porphyromonas gingivalis, Treponema denticola, and Fusobacterium nucleatum, are significantly more abundant in T2D patients compared to normal controls. Furthermore, the Firmicutes/Bacteroidetes (F/B) ratio was higher in T2D patients than in healthy individuals (11), which aligns with the observed trend in newly diagnosed diabetes patients, where there is a decrease in the abundance of Bacteroidetes and an increase in the abundance of Firmicutes in the gut microbiome (31). Previous studies have found that an increase in Firmicutes or a higher F/B ratio is associated with obesity, a risk factor for T2D, as Firmicutes are more efficient than Bacteroidetes at extracting energy from food (32). Chen et al. found that the F/B ratio increased in patients with T2D, suggesting that this ratio may serve as a specific microbial biomarker in Chinese patients with T2D (33). Our findings also revealed that oral microbiomes genetically linked to T2D were predominantly Firmicutes.
Lu et al. observed that Treponema, Prevotella oralis, and Catonella were more abundant in the group with periodontitis and diabetes compared to the systemically healthy group (34). Additionally, Prevotella is more commonly found in the gut microbiome of T2D patients (35). We also confirmed that Prevotella in the oral microbiome is positively correlated with T2D. Previous studies have shown that Prevotella is associated with increased production of branched-chain amino acids (BCAAs). Elevated levels of BCAAs in the blood over the long term are linked to a higher risk of obesity and T2D (36). While Treponema in saliva showed a significant negative correlation. Regarding Catonella, certain species have shown a positive correlation with T2D, while more species have been found to have a negative correlation with the risk of developing T2D. Aggregatibacter actinomycetemcomitans (A.a) has been found to disrupt host mucosal defenses and was identified as one of the pathogenic bacteria involved in periodontitis (37). Castrillon et al. demonstrated that A. a detection was higher in patients with diabetes and periodontitis than in systemically healthy patients without periodontitis, with A.a being associated with periodontitis in diabetic patients (14). Our research confirmed a positive correlation between Aggregatibacter and T2D risk. A.a possesses some putative virulence factors, including leukotoxin that targets and destroys host immune cells. Previous reports have associated A.a with adverse events such as cerebral infarction in diabetic nephropathy patients undergoing hemodialysis (38). Furthermore, compared to control mice, mice infected with A.a exhibit impaired glucose tolerance and insulin resistance, along with alterations in the composition of their gut microbiota (39). Pauljensenia, a Gram-positive, strictly anaerobic, non-spore-forming bacterium from the family Actinomycetaceae. Previous studies have identified Actinomycetaceae as an oral biomarker for T2D (10). We further substantiated the positive correlation between Actinomycetaceae presence in both saliva and tongue and the risk of developing T2D. Its pathogenic role may be associated with glucose metabolism, participating in glycolysis for energy production, and accumulation of intracellular polysaccharides (40, 41), potentially increasing diabetes risk.
A study indicated a reduction in the presence of Haemophilus in the gut of patients diagnosed with T2D (31). The genus Haemophilus in the gut was identified as a defensive element against T2D (42). This finding is consistent with our study, which demonstrated a negative correlation between oral Haemophilus D and the incidence of T2D. Neri Rosario et al. employed machine learning techniques to identify Granulicatella and Prevotella as relevant genera in patients with prediabetes when compared to normoglycemic subjects (43). Previous studies have identified a negative correlation between the Plaque Index (PLI) and TM7x in patients with T2D and periodontitis (44). This suggests that TM7x may inhibit plaque formation or promote plaque clearance, potentially due to the defensive properties of certain TM7x species against T2D. Streptococcus has been identified as a marker bacterium in the oral and gut microbiome of patients with T2D (10). A study conducted on T2D patients in southern Thailand reported significantly higher total counts of salivary and plaque streptococci in diabetics compared to non-diabetics (45). Streptococcus is one of the earliest colonizers of the human body, particularly abundant in the oral cavity (46). It has been reported that pregnant individuals with pregestational diabetes with worse glycemic control were at an increased risk of group B streptococcus (GBS) colonization (46). In our study, various Streptococcus species exhibited mixed genetic causal effects on T2D, highlighting the complexity of the relationship between these bacteria and diabetes.
This study has several limitations. First, our MR analysis is concentrated on populations of East Asian ancestry, and additional validation is necessary to extend these findings to other ethnic groups. Second, factors beyond genetics, such as lifestyle, diet, and environmental influences, can also affect the oral microbiome (47, 48). The instrumental variables may explain only a small fraction of the observed variability, highlighting the need for further multidimensional research to fully comprehend the complex dynamics of the oral microbiome. Third, to ensure a sufficient number of SNPs as instrumental variables for the oral microbiome, we adopted a relatively lenient clumping threshold. Although various MR methods were used for sensitivity analysis and multiple corrections were applied to the results, the potential for some false positives cannot be entirely excluded.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
All data used were publicly available and approved by their respective institutions’ ethics committees. No additional ethical approval from our university was required as the study involved secondary data analysis. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
XL: Writing – review & editing, Writing – original draft, Methodology, Formal analysis, Data curation, Conceptualization. XX: Writing – review & editing, Writing – original draft, Software, Methodology, Data curation, Conceptualization. SS: Writing – original draft, Supervision, Formal analysis, Data curation, Conceptualization. FQ: Writing – review & editing, Writing – original draft, Resources, Project administration, Funding acquisition, Conceptualization.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Sichuan Science and Technology Program (No. 2022YFS0028) and 1•3•5 Project for Disciplines of Excellence, West China Hospital, Sichuan University (No. ZYGD23011).
Acknowledgments
The authors sincerely thank the authors who shared the original dataset in this study.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1452999/full#supplementary-material
References
1. Sacks DB, Arnold M, Bakris GL, Bruns DE, Horvath AR, Kirkman MS, et al. Guidelines and recommendations for laboratory analysis in the diagnosis and management of diabetes mellitus. Diabetes Care. (2023) 46:e151–99. doi: 10.2337/dci23-0036
2. Ahmad E, Lim S, Lamptey R, Webb DR, Davies MJ. Type 2 diabetes. Lancet. (2022) 400:1803–20. doi: 10.1016/S0140-6736(22)01655-5
3. Borgnakke WS, Yl€ostalo PV, Taylor GW, Genco RJ. Effect of periodontal disease on diabetes: systematic review of epidemiologic observational evidence. J Periodontology. (2013) 84:S135–52. doi: 10.1902/jop.2013.1340013
4. Jung Y-S, Shin MH, Kweon SS, Lee YH, Kim OJ, Kim YJ, et al. Periodontal disease associated with blood glucose levels in urban Koreans aged 50 years and older: the Dong-gu study. Gerodontology. (2015) 32:267–73. doi: 10.1111/ger.12107
5. Bharti P, Katagiri S, Nitta H, Nagasawa T, Kobayashi H, Takeuchi Y, et al. Periodontal treatment with topical antibiotics improves glycemic control in association with elevated serum adiponectin in patients with type 2 diabetes mellitus. Obes Res Clin Pract. (2013) 7:e129–38. doi: 10.1016/j.orcp.2011.11.005
6. Tan X, Wang Y, Gong T. The interplay between oral microbiota, gut microbiota and systematic diseases. J Oral Microbiol. (2023) 15:2213112. doi: 10.1080/20002297.2023.2213112
7. Tierney BT, Yang Z, Luber JM, Beaudin M, Wibowo MC, Baek C, et al. The landscape of genetic content in the gut and oral human microbiome. Cell Host Microbe. (2019) 26:283–295.e8. doi: 10.1016/j.chom.2019.07.008
8. Sharma N, Bhatia S, Sodhi AS, Batra N. Oral microbiome and health. AIMS Microbiol. (2018) 4:42–66. doi: 10.3934/microbiol.2018.1.42
9. Mosaddad SA, Mahootchi P, Safari S, Rahimi H, Aghili SS. Interactions between systemic diseases and oral microbiota shifts in the aging community: A narrative review. J Basic Microbiol. (2023) 63:831–54. doi: 10.1002/jobm.202300141
10. Guo X-J, Dai SX, Lou JD, Ma XX, Hu XJ, Tu LP, et al. Distribution characteristics of oral microbiota and its relationship with intestinal microbiota in patients with type 2 diabetes mellitus. Front Endocrinol (Lausanne). (2023) 14:1119201. doi: 10.3389/fendo.2023.1119201
11. Li Y, Qian F, Cheng X, Wang D, Wang Y, Pan Y, et al. Dysbiosis of oral microbiota and metabolite profiles associated with type 2 diabetes mellitus. Microbiol Spectr. (2023) 11:e03796–22. doi: 10.1128/spectrum.03796-22
12. Long J, Cai Q, Steinwandel M, Hargreaves MK, Bordenstein SR, Blot WJ, et al. Association of oral microbiome with type 2 diabetes risk. J Periodontal Res. (2017) 52:636–43. doi: 10.1111/jre.12432
13. Makiura N, Ojima M, Kou Y, Furuta N, Okahashi N, Shizukuishi S, et al. Relationship of Porphyromonas gingivalis with glycemic level in patients with type 2 diabetes following periodontal treatment. Oral Microbiol Immunol. (2008) 23:348–51. doi: 10.1111/j.1399-302X.2007.00426.x
14. Castrillon CA, Hincapie JP, Yepes FL, Roldan N, Moreno SM, Contreras A, et al. Occurrence of red complex microorganisms and Aggregatibacter actinomycetemcomitans in patients with diabetes. J Investig Clin Dent. (2015) 6:25–31. doi: 10.1111/jicd.12051
15. Xiong Y, Zhang F, Zhang Y, Wang W, Ran Y, Wu C, et al. Insights into modifiable risk factors of erectile dysfunction, a wide-angled Mendelian Randomization study. J Adv Res. (2023) 58:149–61. doi: 10.1016/j.jare.2023.05.008. S2090-1232(23)00147–9.
16. Bowden J, Holmes MV. Meta-analysis and Mendelian randomization: A review. Res Synth Methods. (2019) 10:486–96. doi: 10.1002/jrsm.1346
17. Slob EAW, Burgess S. A comparison of robust Mendelian randomization methods using summary data. Genet Epidemiol. (2020) 44:313–29. doi: 10.1002/gepi.22295
18. Liu X, Tong X, Zhu J, Tian L, Jie Z, Zou Y, et al. Metagenome-genome-wide association studies reveal human genetic impact on the oral microbiome. Cell Discovery. (2021) 7:117. doi: 10.1038/s41421-021-00356-0
19. Spracklen CN, Horikoshi M, Kim YJ, Lin K, Bragg F, Moon S, et al. Identification of type 2 diabetes loci in 433,540 East Asian individuals. Nature. (2020) 582:240–5. doi: 10.1038/s41586-020-2263-3
20. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. (2021) 53:1415–24. doi: 10.1038/s41588-021-00931-x
21. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. (2007) 81:559–75. doi: 10.1086/519795
22. Hemani G, Tilling K, Smith GD. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PloS Genet. (2017) 13:e1007081. doi: 10.1101/117101
23. Burgess S, Thompson SG. & CRP CHD Genetics Collaboration. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. (2011) 40:755–64. doi: 10.1093/ije/dyr036
24. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
25. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
26. Burgess S, Foley CN, Allara E, Staley JR, Howson JMM. A robust and efficient method for Mendelian randomization with hundreds of genetic variants. Nat Commun. (2020) 11:376. doi: 10.1038/s41467-019-14156-4
27. Zhao Q, Wang J, Hemani G, Bowden J, Small DS. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann Statist. (2020) 48:1742–69. doi: 10.1214/19-AOS1866
28. Wade WG. The oral microbiome in health and disease. Pharmacol Res. (2013) 69:137–43. doi: 10.1016/j.phrs.2012.11.006
29. Baker JL, Mark Welch JL, Kauffman KM, McLean JS, He X. The oral microbiome: diversity, biogeography and human health. Nat Rev Microbiol. (2024) 22:89–104. doi: 10.1038/s41579-023-00963-6
30. Aemaimanan P, Amimanan P, Taweechaisupapong S. Quantification of key periodontal pathogens in insulin-dependent type 2 diabetic and non-diabetic patients with generalized chronic periodontitis. Anaerobe. (2013) 22:64–8. doi: 10.1016/j.anaerobe.2013.06.010
31. Letchumanan G, Abdullah N, Marlini M, Baharom N, Lawley B, Omar MR, et al. Gut microbiota composition in prediabetes and newly diagnosed type 2 diabetes: A systematic review of observational studies. Front Cell Infect Microbiol. (2022) 12:943427. doi: 10.3389/fcimb.2022.943427
32. Magne F, Gotteland M, Gauthier L, Zazueta A, Pesoa S, Navarrete P, et al. The firmicutes/bacteroidetes ratio: A relevant marker of gut dysbiosis in obese patients? Nutrients. (2020) 12:1474. doi: 10.3390/nu12051474
33. Chen B, Wang Z, Wang J, Su X, Yang J, Zhang Q, et al. The oral microbiome profile and biomarker in Chinese type 2 diabetes mellitus patients. Endocrine. (2020) 68:564–72. doi: 10.1007/s12020-020-02269-6
34. Lu C, Zhao Q, Deng J, Chen K, Jiang X, Ma F, et al. Salivary microbiome profile of diabetes and periodontitis in a chinese population. Front Cell Infect Microbiol. (2022) 12:933833. doi: 10.3389/fcimb.2022.933833
35. Nagpal R, Newman TM, Wang S, Jain S, Lovato JF, Yadav H. Obesity-linked gut microbiome dysbiosis associated with derangements in gut permeability and intestinal cellular homeostasis independent of diet. J Diabetes Res. (2018) 2018:3462092. doi: 10.1155/2018/3462092
36. Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, et al. Metabolite profiles and the risk of developing diabetes. Nat Med. (2011) 17:448–53. doi: 10.1038/nm.2307
37. Fine DH, Patil AG, Velusamy SK. Aggregatibacter actinomycetemcomitans (Aa) under the radar: myths and misunderstandings of aa and its role in aggressive periodontitis. Front Immunol. (2019) 10:728. doi: 10.3389/fimmu.2019.00728
38. Murakami M, Suzuki J, Yamazaki S, Ikezoe M, Matsushima R, Ashigaki N, et al. High incidence of Aggregatibacter actinomycetemcomitans infection in patients with cerebral infarction and diabetic renal failure: a cross-sectional study. BMC Infect Dis. (2013) 13:557. doi: 10.1186/1471-2334-13-557
39. Komazaki R, Katagiri S, Takahashi H, Maekawa S, Shiba T, Takeuchi Y, et al. Periodontal pathogenic bacteria, Aggregatibacter actinomycetemcomitans affect non-alcoholic fatty liver disease by altering gut microbiota and glucose metabolism. Sci Rep. (2017) 7:13950. doi: 10.1038/s41598-017-14260-9
40. Bronzato JD, Bomfim RA, Hayasida GZP, Cúri M, Estrela C, Paster BJ, et al. Analysis of microorganisms in periapical lesions: A systematic review and meta-analysis. Arch Oral Biol. (2021) 124:105055. doi: 10.1016/j.archoralbio.2021.105055
41. Li J, Li Y, Zhou Y, Wang C, Wu B, Wan J. Actinomyces and alimentary tract diseases: A review of its biological functions and pathology. BioMed Res Int. (2018) 2018:3820215. doi: 10.1155/2018/3820215
42. Sun K, Gao Y, Wu H, Huang X. The causal relationship between gut microbiota and type 2 diabetes: a two-sample Mendelian randomized study. Front Public Health. (2023) 11:1255059. doi: 10.3389/fpubh.2023.1255059
43. Neri-Rosario D, Martínez-López YE, Esquivel-Hernández DA, Sánchez-Castañeda JP, Padron-Manrique C, Vázquez-Jiménez A, et al. Dysbiosis signatures of gut microbiota and the progression of type 2 diabetes: a machine learning approach in a Mexican cohort. Front Endocrinol (Lausanne). (2023) 14:1170459. doi: 10.3389/fendo.2023.1170459
44. Gu M, Wang P, Xiang S, Xu D, Jin C, Jiang Z, et al. Effects of type 2 diabetes and metformin on salivary microbiota in patients with chronic periodontitis. Microb Pathog. (2021) 161:105277. doi: 10.1016/j.micpath.2021.105277
45. Kampoo K, Teanpaisan R, Ledder RG, McBain AJ. Oral bacterial communities in individuals with type 2 diabetes who live in southern Thailand. Appl Environ Microbiol. (2014) 80:662–71. doi: 10.1128/AEM.02821-13
46. Field C, Bank TC, Spees CK, Germann K, Landon MB, Gabbe S, et al. Association between glycemic control and group B streptococcus colonization among pregnant individuals with pregestational diabetes. Am J Reprod Immunol. (2023) 90:e13779. doi: 10.1111/aji.13779
47. Verma D, Garg PK, Dubey AK. Insights into the human oral microbiome. Arch Microbiol. (2018) 200:525–40. doi: 10.1007/s00203-018-1505-3
Keywords: oral microbiome, type 2 diabetes, Mendelian randomization, genetic variation, causal inference
Citation: Lyu X, Xu X, Shen S and Qin F (2024) Genetics causal analysis of oral microbiome on type 2 diabetes in East Asian populations: a bidirectional two-sample Mendelian randomized study. Front. Endocrinol. 15:1452999. doi: 10.3389/fendo.2024.1452999
Received: 22 June 2024; Accepted: 06 August 2024;
Published: 23 August 2024.
Edited by:
Weiwei Liu, Chongqing Medical University, ChinaReviewed by:
Xiaodong Liu, The Chinese University of Hong Kong, ChinaMeng Jia, University of Pennsylvania, United States
Ping Zheng, The University of Melbourne, Australia
Cosmin Mihai Vesa, University of Oradea, Romania
Copyright © 2024 Lyu, Xu, Shen and Qin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Feng Qin, qinfeng@scu.edu.cn
†These authors have contributed equally to this work and share first authorship