Next Article in Journal
Polymorphisms in TNF Receptor Superfamily 1B (TNFRSF1B:rs3397) are Linked to Mycobacterium avium paratuberculosis Infection and Osteoporosis in Rheumatoid Arthritis
Next Article in Special Issue
Characterization of the Gut Microbiota of Individuals at Different T2D Stages Reveals a Complex Relationship with the Host
Previous Article in Journal
Isolation and Molecular Typing of Leishmania infantum from Phlebotomus perfiliewi in a Re-Emerging Focus of Leishmaniasis, Northeastern Italy
Previous Article in Special Issue
Gut Microbiome Modulation Based on Probiotic Application for Anti-Obesity: A Review on Efficacy and Validation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 7
Issue 12
10.3390/microorganisms7120645
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Profiling the Oral Microbiome and Plasma Biochemistry of Obese Hyperglycemic Subjects in Qatar

by
Muhammad U. Sohail
1,*,
Mohamed A. Elrayess
1,
Asma A. Al Thani
1,2,
Maha Al-Asmakh
1,2 and
Hadi M. Yassine
1,2,*
1
Biomedical Research Center, Qatar University, Doha 2713, Qatar
2
Department of Biomedical Sciences, College of Health Sciences, QU Health, Qatar University, Doha 2713, Qatar
*
Authors to whom correspondence should be addressed.
Microorganisms 2019, 7(12), 645; https://doi.org/10.3390/microorganisms7120645
Submission received: 18 September 2019 / Revised: 19 November 2019 / Accepted: 28 November 2019 / Published: 3 December 2019
(This article belongs to the Special Issue Gut Microbiota: Its Role in Diabetes and Obesity)
Browse Figures
Review Reports Versions Notes

Abstract

:
The present study is designed to compare demographic characteristics, plasma biochemistry, and the oral microbiome in obese (N = 37) and lean control (N = 36) subjects enrolled at Qatar Biobank, Qatar. Plasma hormones, enzymes, and lipid profiles were analyzed at Hamad Medical Cooperation Diagnostic Laboratory. Saliva microbiome characterization was carried out by 16S rRNA amplicon sequencing using Illumina MiSeq platform. Obese subjects had higher testosterone and sex hormone-binding globulin (SHBG) concentrations compared to the control group. A negative association between BMI and testosterone (p < 0.001, r = −0.64) and SHBG (p < 0.001, r = −0.34) was observed. Irrespective of the study groups, the oral microbiome was predominantly occupied by Streptococcus, Prevotella, and Veillonella species. A generalized linear model revealed that the Firmicutes/Bacteroidetes ratio (2.25 ± 1.83 vs. 1.76 ± 0.58; corrected p-value = 0.04) was higher, and phylum Fusobacteria concentration (4.5 ± 3.0 vs. 6.2 ± 4.3; corrected p-value = 0.05) was low in the obese group compared with the control group. However, no differences in microbiome diversity were observed between the two groups as evaluated by alpha (Kruskal–Wallis p ≥ 0.78) and beta (PERMANOVA p = 0.37) diversity indexes. Certain bacterial phyla (Acidobacteria, Bacteroidetes, Fusobacteria, Proteobacteria, Spirochaetes, and Firmicutes/Bacteroidetes) were positively associated (p = 0.05, r ≤ +0.5) with estradiol, fast food consumption, creatinine, breastfed during infancy, triglycerides, and thyroid-stimulating hormone concentrations. In conclusion, no differences in oral microbiome diversity were observed between the studied groups. However, the Firmicutes/Bacteroidetes ratio, a recognized obesogenic microbiome trait, was higher in the obese subjects. Further studies are warranted to confirm these findings in a larger cohort.

1. Introduction

Diabetes mellitus is a chronic metabolic disorder with devastating consequences. Genetics and environment are the two main etiological factors that partially or collectively contribute to type 2 diabetes mellitus (T2DM) pathology. In most cases, several genes and epigenetic factors are involved in the development of the disease [1]. However, the significant rise in T2DM incidences suggests that environmental factors play a relatively more important role in the disease pathology than we previously thought. Many environmental factors, such as diet and physical activity, play a crucial role in T2DM development and pathogenesis [2]. Accordingly, a better understanding of the disease pathology and etiology could improve prevention and treatment options. However, no single-time, definitive therapy is available for the treatment of diabetes, and mostly the disease is managed by continuous, long-term medications. In this regard, exploring potential preventative and therapeutic measures requires a better understanding of the long pre-diabetic phase that can last for several years [3]. Hyperglycemia in the pre-diabetic stage is correlated with obesity, blood pressure, and dyslipidemia [4]. However, like diabetes, obesity also does not have a well-defined etiology or pathogenesis and, in many cases, obesity is considered as a prodrome of insulin resistance (IR) in diabetes [5].
Evaluation and characterization of the mucosal microbiome as an environmental factor for obesity and diabetes has received much attraction in recent years. Several studies showed that alterations in the gut microbiome could be liked with IR and obesity and that the microbiome could be a hopeful therapeutic goal for prevention of these metabolic syndromes [6]. Furthermore, microbiome studies in diabetic patients revealed a strong correlation between microbiome diversity and obese/diabetic phenotypes [7,8]. Particularly, elevated Firmicutes and Proteobacteria populations and a drop in Bacteroidetes levels are reported in obese and IR subjects [9,10]. However, most of the published literature is focused on patients with established diabetes, and very few studies have explored the microbiome of the pre-diabetic phase of hyperglycemia [11,12]. Although the oral microbiome has been shown to recapitulate the gut microbiome [13], it is relatively less explored. Accordingly, the association between the oral microbiome and disease prognosis is not well defined.
The oral microbiome has been implicated in several oral cavity (dental caries, periodontitis, endodontic, alveolar osteitis, and tonsillitis) and systemic diseases (cardiovascular disease, stroke, pneumonia, and diabetes) [14,15,16,17]. Long et al. [18] reported that several bacterial taxa in the phylum Actinobacteria are negatively associated with obesity and diabetes. Similarly, many other studies have also reported phylogenetic differences in the microbiome profile of obese and lean humans [19,20]. However, the compositional changes in microbial ecology that are observed in diabetic and obese subjects are often not in agreement across these studies. This makes it challenging to interpret the role of the microbiome in disease development and pathogenicity. We, therefore, sought to characterize the oral microbial community as a biomarker of obesity and diabetes. To achieve this objective, we studied the oral microbiome in obese hyperglycemic subjects.

2. Materials and Methods

2.1. Study Cohort

This study enrolled obese pre-diabetic (hyperglycemic) adult males and age-matched healthy individuals who voluntarily enrolled at Qatar Biobank (QBB; https://www.qatarbiobank.org.qa/). The inclusion criteria for the participants were: age ≥ 30 years, body mass index (BMI) characterized as healthy (<25) or obese (≥30), and not clinically diagnosed as diabetic. All participants were Qatari nationals. Subjects who had a history of any metabolic disease or had consumed antibiotics or steroids in the preceding three months of sampling were excluded from this study. Seventy-three volunteers participated in this study, and were categorized according to their BMI into obese (N = 37) and lean control (N = 36) groups. Information on disease history, oral health (Table S4), medication, feeding habits, physical activity, socioeconomic status, smoking habit, and family disease history was collected on the prescribed questionnaire. All participants in this study signed an informed consent form for the use of their information for sample analysis as anonymous volunteers. The institutional ethics review boards of the QBB (MOPH-QBB-IRB-011) and Qatar University (QU-IRB 969-A/18) approved this study in compliance with participant anonymity, research ethical, moral, and biosafety standards.

2.2. Plasma Biochemistry

Blood samples were collected in anticoagulant-coated evacuated tubes (BD, Mississauga, ON, Canada). Plasma concentrations of the hormones, enzymes, and lipid markers were analyzed at Hamad Medical Corporation (HMC) diagnostic laboratory using Cobas 6000 analyzer (Roche Diagnostics), as described previously [21,22,23]. A complete list of instruments and reagents used for plasma biochemistry is available in the supplementary file (Table S5).

2.3. 16S rRNA Sequencing

Saliva samples were collected from the participants by spitting saliva in sterile tubes. The samples were transported on ice from QBB to the Biomedical Research Center (BRC) of Qatar University (QU). Only 69 (Obese 36 and Control 33) saliva samples were available for sequencing as three participants did not provide a saliva sample, and we lost one sample during DNA extraction. Genomic DNA was extracted from the samples using a commercially available DNA extraction kit (QIAamp DNA Mini Kit, 51306, Germantown, MD, USA). The quality and quantity of the DNA were evaluated using NanoDrop-2000 (Thermo Fisher Scientific, Waltham Massachusetts, US) and Qubit-4 (Life Technologies, Carlsbad, California, US). The DNA samples were then subjected to 16S rRNA library preparation protocol using an Illumina Nextera XT Library Prep. Kit (FC-131-1002, Illumina Inc., San Diego, CA, USA). In brief, the V3–V4 region of the 16S rRNA gene was amplified using a 337F/805R primer pair [24], followed by an Illumina two-step amplification library preparation strategy [25]. Prepared libraries were cleaned and normalized using magnetic beads (Agencourt Ampure XP, Beckman Coulter, IN, USA). Finally, all libraries were pooled together in equal volumes and denatured using 0.2 N NaOH. The sequencing was performed on Illumina MiSeq (San Diego, CA, USA) using a 600 cycles v3 kit (MS-102-3003; Illumina, San Diego, CA, USA).

2.4. Bioinformatics

The data were obtained as paired-end reads. Forward and reverse reads were merged before analysis. The data were subjected to quality filtration and chimera removal using the DADA2 plugin implemented in QIIME2 [26,27]. The first thirteen bases of the forward and reverse reads were trimmed, while truncation was performed at 255 bases to allow sufficient overlapping of the forward and reverse reads. The DADA2 plugin generated 7110 sequence features, defined as unique 16S rRNA gene sequence variants. Phylogenetic diversity analysis was performed on QIIME2 using q2-phylogeny plugin that wraps mafft-fasttree program. Taxonomic classification was performed using Greengenes 13-8 database as the reference [28,29]. Each feature sequence was assigned taxonomy for >97% identity (or < 3% divergence) at the species, >95% at the genus, >90% at the family, >85% at the order, >80% at the class, and >77% at the phylum level [30].

2.5. Statistical Analysis

Demographic data were arranged from lowest to highest possible values and categorically numbered from 1 to 6 (Table S1). The Mann–Whitney U test was performed to compare the mean differences in plasma biochemistry and demographic characteristics between the study groups. Spearman′s rank correlation coefficient was applied to measure the correlation between these variables. The Benjamini–Hochberg false discovery rate (FDR) was used to perform multiple comparisons of the p-values, and an adjusted p < 0.05 was considered statistically significant [31]. Core microbiome diversity analysis was performed at 28,153 sequencing depth using observed_OTUs, faith_PD, and Shannon indexes for alpha diversity analysis. Weighted_unifrac unweighted_unifrac, and Bray_Curtis_matrix indexes were used for beta diversity analysis to generate principal coordinates (PCoA) plots. The Wallis test was used to compare within-sample diversity (alpha diversity). Beta diversity (between samples) significance analysis was performed using the PERMANOVA test. For taxonomic assessment of microbiota, taxa were represented at a particular phylogenetic resolution (phylum, family, and genus). Only those taxa that had a relative abundance of at least 0.5% in any of the two groups were included in statistical analysis, and remaining data were discarded. The Kruskal–Wallis test was applied to compare taxonomic differences at all hierarchical levels. The association between microbiome and biochemical parameters was assessed using Spearman′s rank correlation coefficient and stepwise linear regression. A power calculation based on a previous similar study [32] indicated that a sample size of 35 per group has 95% power to detect a minimal difference of 11% in microbiome composition between diabetics and controls, with 10% deviation from the mean value (σ) at a level of α = 0.01 [33].

3. Results

The demographic and plasma biochemistry features of the study groups are presented in Table 1. The age distribution was similar across the two groups, with a mean age of 36.75 ± 8.0 years in the obese and 36.52 ± 7.9 years in the control group. The mean BMI was 35.65 ± 4.9 and 22.73 ± 1.5 for the obese and control groups, respectively. Fasting glucose and Hemoglobin A1c (HbA1c) levels were higher in the obese group compared to the control (p ≤ 0.001). Testosterone, sex hormone-binding globulin (SHBG), and insulin concentrations were lower (p ≤ 0.001), whereas estradiol hormone concentrations were higher (p = 0.03) in the obese group compared to the control group. No differences were observed for dietary habits, physical activity, smoking, and antibiotic usage between the study groups (Table 1). Spearman’s correlation analysis revealed a negative correlation between testosterone and BMI (r = −0.68, p ≤ 0.001) and insulin (r = −0.50, p ≤ 0.001). On the other hand, a positive correlation was observed between BMI and estradiol (r = 0.44, p ≤ 0.01) (Figure 1).
High-throughput 16S rRNA amplicon sequencing yielded 6,790,910 sequences for all analyzed samples (n = 69, median ± SD = 96,330 ± 20,849.7). After quality control, 4,028,561 sequences of good quality (Phred quality score ≥ ASCII 30), belonging to 7110 features, were used for further analysis. All obtained sequences were demultiplexed and deposited to the NCBI Sequence Read Archive (SRA) for future reference under study accession number PRJNA587625.
Figure 2 depicts the alpha diversity indexes of bacterial communities in the two study groups. The Faith_PD, Shannon, and observed_OTUs indexes were utilized to determine taxonomic diversity (species richness and evenness) within the samples. None of the diversity measures significantly differed between the obese and the control groups (Kruskal–Wallis p ≥ 0.05). Similarly, weighted_unifrac distance matrix-based PCoA plots showed that no distinct clustering pattern was present between the microbiome of the study samples (Figure 3). The Permutational multivariate analysis of variance (PERMANOVA) tests applied in beta diversity indexes did not show significant differences (PERMANOVA p = 0.37) among the studied groups.
The taxonomic analysis was performed using the greengenes 13-8 database as a reference. In total, 50 phyla, 124 classes, 188 orders, 412 families, 595 genera, and 676 species were identified in the 69 saliva samples. In descending order, the four most dominant phyla present in the oral microbiome of the study population were Firmicutes (43.3%), Bacteroidetes (25.2%), Proteobacteria (10.7%), and Actinobacteria (8.2%). Streptococcus spp., (22.2%) Prevotella melaninogenica (13.4%), and Veillonella dispar (6.6%) were the most abundant bacterial species that belonged to family Streptococcaceae (22.5%), Prevotellaceae (18.1%), and Veillonellaceae (10.1%), respectively. The Kruskal–Wallis test showed no differences (p ≥ 0.05) in microbiome population at any phylogenetic level between the studied groups (Figure 4). However, a generalized linear model revealed that the Firmicutes/Bacteroidetes ratio was significantly higher in obese IR subjects when compared to insulin sensitive control (2.25 ± 1.83 vs. 1.76 ± 0.58; corrected p-value = 0.04) after correcting for potential confounders including HbA1c, insulin, and triglycerides. Similarly, Fusobateria also exhibited a significant difference between the lean and obese groups (6.2 ± 4.3 vs. 4.5 ± 3.0, respectively; corrected p-value = 0.05) after correcting for potential confounders (HbA1c, insulin, and triglycerides) (Figure 5).
Figure 6 presents Spearman’s correlation analysis between microbiome, population demographic characteristics, and plasma biochemistry data. Irrespective of the treatment group, certain bacterial phyla (Acidobacteria, Bacteroidetes, Fusobacteria, Proteobacteria, Spirochaetes, and Firmicutes/Bacteroidetes) were positively associated (corrected p-value = 0.05, r ≤ +0.5) with estradiol, fast food consumption, creatinine, breastfed during infancy, triglycerides, and thyroid-stimulating hormone (TSH) concentrations. Particularly, the phyla Proteobacteria and Acidobacteria were positively associated with elevated estradiol concentrations. The Firmicutes/Bacteroidetes ratio was negatively associated with HDL cholesterol and positively associated with TSH concentrations. In order to identify the best predictors of microbial taxon association with the demographic parameters, a stepwise linear regression was carried out. The regression model indicated that regardless of BMI, estradiol and HDL cholesterol were the best predictors of Acidobacteria, whereas estradiol, HDL cholesterol and triiodothyronine were the best predictors of Firmicutes. The model also revealed that creatinine was the best predictor of Fusobacteria regardless of BMI (Table 2).

4. Discussion

To date, limited attention has been directed to studying the oral microbiome during the pre-diabetic phase, when metabolic and immune-inflammatory perturbation are silently underway. To fill this gap, we investigated the association between the oral microbiome and metabolic markers of clinically healthy adults categorized as obese hyperglycemic and healthy normoglycemic subjects in Qatar.
We observed significant differences between the two investigated groups in terms of HbA1c, glucose, insulin, Homeostatic Model Assessment of Insulin Resistance (HOMA-IR), C-peptide, triglycerides, and HDL cholesterol levels, even though both the groups had no differences in demographic characteristics. On the other hand, we observed no association between these obesogenic parameters and daily activity, breastfed during infancy, smoking, and fast food consumption. Furthermore, we observed no association between BMI/hyperglycemia and thyroid gland activity, although this gland has a principal role in the regulation of cellular metabolism. It is worth noting that obese hyperglycemic males enrolled in this study had hypogonadism (drop in testosterone and SHBG concentrations) and elevated levels of the estradiol hormone. This observed negative association between testosterone and BMI/hyperglycemia is now a well-known phenomenon that describes a greater tendency to acquire hypogonadism in T2DM males [34]. Furthermore, a low testosterone concentration has been recognized as a reliable predictor of IR and the probability of developing T2DM in the future [35]. Meyer et al. (2017) suggested that in male obese subjects, testosterone is converted to estrogen by the aromatase enzyme produced in adipose tissue, which eventually results in male hypogonadism [36]. On the other hand, the estrogen hormone is a strong regulator of body weight and insulin sensitivity through the activation of G-protein-coupled estrogen receptors and the overactivation of SHBG [37]. Narrated together, all three sex hormones have a significant implication on T2DM development that corresponds directly to host metabolic rate and thyroid hormones activities [38]. Thyroid gland hormones are involved in cellular metabolism, and the overactivity of the gland is considered as a risk factor for diabetes development, especially in pre-diabetes subjects [39].
Since the microbiome has been enormously described as a co-regulator of host metabolism and adipogenesis [40], we ran microbiome analyses on saliva samples and studied its association with metabolic markers as indicated above. We observed no difference in microbiome alpha and beta diversity indexes between the two groups. Weighted_unifrac analysis of beta diversity revealed that approximately 60% of the samples were scattered along axis 1, which indicates that no changes in microbiome diversity as per hyperglycemia were present. However, after correcting for potential confounders, including HbA1c, insulin, and triglycerides, the Firmicutes/Bacteroidetes ratio was found to be significantly higher in obese IR subjects when compared to insulin-sensitive controls. Similar observations have previously been reported by Demmer et al. (2017), who studied the subgingival microbiome of oral infections, glucose intolerance and IR in non-diabetic adults [8]. The authors reported that the increase in Firmicutes or drop in Bacteroidetes population is associated with periodontitis and systemic inflammations [8]. An increase in the Firmicutes/Bacteroidetes ratio is considered a prognostic factor for the development of T2DM [41], as it is observed that most of the bacterial genera in the phylum Firmicutes can contribute to host weight gain and obesity [42].
In general, Firmicutes, Bacteroidetes, and Proteobacteria were the most abundantly present bacterial phyla in the salivary microbiome of our study samples. Previous studies also showed that these bacterial phyla are most abundantly present in our mucosal communities and are associated with host metabolic rate and energy homeostasis [42,43]. However, there are also contradictory findings in the literature about microbial phylogenetic association with T2DM. Xiao et al. (2017) found that the oral microbiome of diabetic and pre-diabetes mice groups is different from normoglycemic mice [44]. However, as we observed, they also found higher population abundances of Firmicutes, Bacteroidetes, and Proteobacteria. Similarly, Long et al. (2017) found that although Firmicutes are the most abundantly present bacteria in oral microbiome, bacterial taxa found in the phylum Actinobacteria are associated with the risk of T2DM development [18]. In contrast, Anbalagan et al. [43] observed that the oral microbiome of T2DM patients was not different from healthy controls. These controversies in observations could be attributed to different study designs and sampling methods. In our study, the lack of changes in microbiome alpha and beta diversity may be associated with the sampling method or demographic characteristics of the participants. Two major contributors to microbiome diversity, food and physical activity, were constant in our study population [45]. Furthermore, saliva samples were obtained using the spitting technique; however, many studies found that microbiome composition could change due to different sampling procedures: saliva/spitting, dental surface, inner cheeks, and lingual swabbing were different [46,47].
We performed a correlation analysis between demographic data and the microbiome at the phylum level. It was observed that the phylum Firmicutes was negatively associated with estradiol hormone levels. The Firmicutes/Bacteroidetes ratio was positively associated with triglycerides and TSH concentrations and negatively associated with HDL cholesterol. However, no correlation between microbiome diversity and BMI, plasma glucose, or Hb1Ac was observed. Therefore, to observe the predictors of microbiome change, the subsequent analysis of microbiome association was performed without consideration of confounding factors. The regression model indicated that regardless of BMI, estradiol and HDL cholesterol were the best predictors of Acidobacteria, whereas estradiol, HDL cholesterol, and triiodothyronine were the best predictors of Firmicutes. The model also revealed that creatinine was the best predictor of Fusobacteria regardless of BMI. These phyla are associated with host metabolism and energy harvest and, therefore, may support obesogenic phenotype and IR [48]. Previously, similar observations were reported by Si et al. (2017), who researched the oral microbiome when exploring biomarkers of metabolic syndrome and reported a linear correlation of HDL cholesterol and triglycerides with members of the phyla Firmicutes and Proteobacteria [49].
In conclusion, we observed that obese pre-diabetic male subjects had significantly low testosterone and sex hormone-binding globulin that may compromise their sexual activity. Overall, oral microbial ecology was highly diverse, including 7110 features at different hierarchal levels. Certain bacterial phyla were associated with reproductive and metabolic hormones, triglycerides, and HDL cholesterol concentrations. Significant differences in the Firmicutes/Bacteroidetes ratio were observed between the pre-diabetic and control groups. The association between the oral microbiome and host metabolic health identified in this study may be advantageous as the key aid in the early prediction of T2DM. More studies with larger samples are warranted to confirm these findings in both genders and different ethnicities.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/7/12/645/s1: Table S1: Data provenance*; Table S2: Dietary habits of the study population—the consumption of dairy products; Table S3: Dietary habits of the study population—fast food consumption; Table S4: Oral health of the participants; Table S5: List of instruments and reagent kits used for serum biochemistry analyses.

Author Contributions

M.U.S. conceived the concept and performed laboratory work and bioinformatics. M.A.E. performed statistical analyses. H.M.Y. developed the concept and supervised lab work and data analysis. M.A.-A. and A.A.A.T. contributed to writing the manuscript. All authors read the manuscript and approved it before submission.

Funding

This research was funded by Qatar National Research Fund through grant number UREP23-062-3-014.

Acknowledgments

The authors acknowledge the undergraduate students (A.A.B., W.A.-A., M.I., B.S., O.T.J. and J.A.) who were part of this project and contributed to lab work. The authors also acknowledge the support of the Qatar BioBank research team for contributing to sample collection, laboratory analysis, and data arrangement.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fuchsberger, C.; Flannick, J.; Teslovich, T.M.; Mahajan, A.; Agarwala, V.; Gaulton, K.J.; Ma, C.; Fontanillas, P.; Moutsianas, L.; McCarthy, D.J. The genetic architecture of type 2 diabetes. Nature 2016, 536, 41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Colberg, S.R.; Sigal, R.J.; Yardley, J.E.; Riddell, M.C.; Dunstan, D.W.; Dempsey, P.C.; Horton, E.S.; Castorino, K.; Tate, D.F. Physical activity/exercise and diabetes: A position statement of the American Diabetes Association. Diabetes Care 2016, 39, 2065–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tabák, A.G.; Herder, C.; Rathmann, W.; Brunner, E.J.; Kivimäki, M. Prediabetes: A high-risk state for developing diabetes. Lancet 2012, 379, 2279. [Google Scholar] [CrossRef] [Green Version]
  4. Collaboration, E.R.F. Diabetes mellitus, fasting blood glucose concentration, and risk of vascular disease: A collaborative meta-analysis of 102 prospective studies. Lancet 2010, 375, 2215–2222. [Google Scholar] [CrossRef] [Green Version]
  5. Sohail, M.U.; Yassine, H.M.; Sohail, A.; Al Thani, A.A. Impact of Physical Exercise on Gut Microbiome, Inflammation, and the Pathobiology of Metabolic Disorders. Rev. Diabet. Stud. RDS 2019, 15, 35–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Wang, J.-Z.; Du, W.-T.; Xu, Y.-L.; Cheng, S.-Z.; Liu, Z.-J. Gut microbiome-based medical methodologies for early-stage disease prevention. Microb. Pathog. 2017, 105, 122–130. [Google Scholar] [CrossRef]
  7. Wang, J.; Jia, H. Metagenome-wide association studies: Fine-mining the microbiome. Nat. Rev. Microbiol. 2016, 14, 508. [Google Scholar] [CrossRef]
  8. Demmer, R.T.; Breskin, A.; Rosenbaum, M.; Zuk, A.; LeDuc, C.; Leibel, R.; Paster, B.; Desvarieux, M.; Jacobs, D.R., Jr.; Papapanou, P.N. The subgingival microbiome, systemic inflammation and insulin resistance: The oral infections, glucose intolerance and insulin resistance study. J. Clin. Periodontol. 2017, 44, 255–265. [Google Scholar] [CrossRef]
  9. Ding, S.; Chi, M.M.; Scull, B.P.; Rigby, R.; Schwerbrock, N.M.; Magness, S.; Jobin, C.; Lund, P.K. High-fat diet: Bacteria interactions promote intestinal inflammation which precedes and correlates with obesity and insulin resistance in mouse. PLoS ONE 2010, 5, e12191. [Google Scholar] [CrossRef] [Green Version]
  10. Turnbaugh, P.J.; Ley, R.E.; Mahowald, M.A.; Magrini, V.; Mardis, E.R.; Gordon, J.I. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 2006, 444, 1027. [Google Scholar] [CrossRef]
  11. Allin, K.H.; Tremaroli, V.; Caesar, R.; Jensen, B.A.; Damgaard, M.T.; Bahl, M.I.; Licht, T.R.; Hansen, T.H.; Nielsen, T.; Dantoft, T.M. Aberrant intestinal microbiota in individuals with prediabetes. Diabetologia 2018, 61, 810–820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, R.-R.; Xu, Y.-S.; Ji, M.-M.; Zhang, L.; Li, D.; Lang, Q.; Zhang, L.; Ji, G.; Liu, B.-C. Association of the oral microbiome with the progression of impaired fasting glucose in a Chinese elderly population. J. Oral Microbiol. 2019, 11, 1605789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhang, X.; Zhang, D.; Jia, H.; Feng, Q.; Wang, D.; Liang, D.; Wu, X.; Li, J.; Tang, L.; Li, Y. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat. Med. 2015, 21, 895. [Google Scholar] [CrossRef] [PubMed]
  14. Blekkenhorst, L.C.; Bondonno, N.P.; Liu, A.H.; Ward, N.C.; Prince, R.L.; Lewis, J.R.; Devine, A.; Croft, K.D.; Hodgson, J.M.; Bondonno, C.P. Nitrate, the oral microbiome, and cardiovascular health: A systematic literature review of human and animal studies. Am. J. Clin. Nutr. 2018, 107, 504–522. [Google Scholar] [CrossRef] [Green Version]
  15. Michaud, D.S.; Izard, J.; Rubin, Z.; Johansson, I.; Weiderpass, E.; Tjønneland, A.; Olsen, A.; Overvad, K.; Boutron-Ruault, M.C.; Clavel-Chapelon, F. Lifestyle, dietary factors, and antibody levels to oral bacteria in cancer-free participants of a European cohort study. Cancer Causes Control 2013, 24, 1901–1909. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Vanhatalo, A.; Blackwell, J.R.; L’Heureux, J.E.; Williams, D.W.; Smith, A.; van der Giezen, M.; Winyard, P.G.; Kelly, J.; Jones, A.M. Nitrate-responsive oral microbiome modulates nitric oxide homeostasis and blood pressure in humans. Free Radic. Biol. Med. 2018, 124, 21–30. [Google Scholar] [CrossRef] [PubMed]
  17. Zarco, M.; Vess, T.; Ginsburg, G. The oral microbiome in health and disease and the potential impact on personalized dental medicine. Oral Dis. 2012, 18, 109–120. [Google Scholar] [CrossRef] [PubMed]
  18. Long, J.; Cai, Q.; Steinwandel, M.; Hargreaves, M.K.; Bordenstein, S.R.; Blot, W.J.; Zheng, W.; Shu, X.O. Association of oral microbiome with type 2 diabetes risk. J. Periodontal Res. 2017, 52, 636–643. [Google Scholar] [CrossRef] [PubMed]
  19. Turnbaugh, P.J.; Hamady, M.; Yatsunenko, T.; Cantarel, B.L.; Duncan, A.; Ley, R.E.; Sogin, M.L.; Jones, W.J.; Roe, B.A.; Affourtit, J.P. A core gut microbiome in obese and lean twins. Nature 2009, 457, 480. [Google Scholar] [CrossRef] [Green Version]
  20. Zupancic, M.L.; Cantarel, B.L.; Liu, Z.; Drabek, E.F.; Ryan, K.A.; Cirimotich, S.; Jones, C.; Knight, R.; Walters, W.A.; Knights, D. Analysis of the gut microbiota in the old order Amish and its relation to the metabolic syndrome. PLoS ONE 2012, 7, e43052. [Google Scholar] [CrossRef]
  21. Al Kuwari, H.; Al Thani, A.; Al Marri, A.; Al Kaabi, A.; Abderrahim, H.; Afifi, N.; Qafoud, F.; Chan, Q.; Tzoulaki, I.; Downey, P. The Qatar Biobank: Background and methods. BMC Public Health 2015, 15, 1208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Al-Dabhani, K.; Tsilidis, K.; Murphy, N.; Ward, H.; Elliott, P.; Riboli, E.; Gunter, M.; Tzoulaki, I. Prevalence of vitamin D deficiency and association with metabolic syndrome in a Qatari population. Nutr. Diabetes 2017, 7, e263. [Google Scholar] [CrossRef] [PubMed]
  23. Thani, A.A.; Fthenou, E.; Paparrodopoulos, S.; Marri, A.A.; Shi, Z.; Qafoud, F.; Afifi, N. Qatar Biobank Cohort Study: Study Design and First Results. Am. J. Epidemiol. 2019. [Google Scholar] [CrossRef] [PubMed]
  24. Anwar, H.; Suchodolski, J.S.; Ullah, M.I.; Hussain, G.; Shabbir, M.Z.; Mustafa, I.; Sohail, M.U. Shiitake Culinary-Medicinal Mushroom, Lentinus edodes (Agaricomycetes), Supplementation Alters Gut Microbiome and Corrects Dyslipidemia in Rats. Int. J. Med. Mushrooms 2019, 21. [Google Scholar] [CrossRef] [PubMed]
  25. Sohail, M.U.; Hume, M.E. Evaluation of Antimicrobial Action of Chitosan and Acetic Acid on Broiler Cecal Bacterial Profiles in Anaerobic Cultures Inoculated with Salmonella Typhimurium. J. Appl. Poult. Res. 2018, 28, 176–183. [Google Scholar] [CrossRef]
  26. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581. [Google Scholar] [CrossRef] [Green Version]
  27. Wajid, B.; Sohail, M.U.; Ekti, A.R.; Serpedin, E. The A, C, G, and T of Genome Assembly. BioMed Res. Int. 2016, 2016. [Google Scholar] [CrossRef] [Green Version]
  28. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F. QIIME 2: Reproducible, interactive, scalable, and extensible microbiome data science. Peer J. Preprints 2018, 6, e27295v2. [Google Scholar] [CrossRef]
  29. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  30. Suchodolski, J.S.; Foster, M.L.; Sohail, M.U.; Leutenegger, C.; Queen, E.V.; Steiner, J.M.; Marks, S.L. The fecal microbiome in cats with diarrhea. PLoS ONE 2015, 10, e0127378. [Google Scholar] [CrossRef] [Green Version]
  31. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  32. Larsen, N.; Vogensen, F.K.; Van Den Berg, F.W.; Nielsen, D.S.; Andreasen, A.S.; Pedersen, B.K.; Al-Soud, W.A.; Sørensen, S.J.; Hansen, L.H.; Jakobsen, M. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PLoS ONE 2010, 5, e9085. [Google Scholar] [CrossRef] [PubMed]
  33. Brant, R. Hypothesis Testing: Two-Sample Inference-Estimation of Sample Size and Power for Comparing Two Means, Fundamentals in Biostatistics. 2013. Available online: https://www.stat.ubc.ca/~rollin/stats/ssize/n2.html (accessed on 17 September 2019).
  34. Al Hayek, A.A.; Robert, A.A.; Alshammari, G.; Hakami, H.; Al Dawish, M.A. Assessment of hypogonadism in men with type 2 diabetes: A cross-sectional study from Saudi Arabia. Clin. Med. Insights Endocrinol. Diabetes 2017, 10, 1179551417710209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Asare-Anane, H.; Ofori, E.; Agyemang, Y.; Oppong, S.; Tagoe, E.; Bani, S.; Ateku, R.; Bawa, T. Obesity and testosterone levels in Ghanaian men with type 2 diabetes. Clin. Diabetes 2014, 32, 61–65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Meyer, M.R.; Clegg, D.J.; Prossnitz, E.R.; Barton, M. Obesity, insulin resistance and diabetes: Sex differences and role of oestrogen receptors. Acta Physiol. 2011, 203, 259–269. [Google Scholar] [CrossRef] [PubMed]
  37. Wallace, I.R.; McKinley, M.C.; Bell, P.M.; Hunter, S.J. Sex hormone binding globulin and insulin resistance. Clin. Endocrinol. 2013, 78, 321–329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Afsar, B.; Karaca, H. The relationship between insulin, insulin resistance, parathyroid hormone, cortisol, testosterone, and thyroid function tests in the presence of nephrolithiasis: A comprehensive analysis. Cent. Eur. J. Urol. 2014, 67, 58. [Google Scholar]
  39. Chaker, L.; Ligthart, S.; Korevaar, T.I.; Hofman, A.; Franco, O.H.; Peeters, R.P.; Dehghan, A. Thyroid function and risk of type 2 diabetes: A population-based prospective cohort study. BMC Med. 2016, 14, 150. [Google Scholar] [CrossRef] [Green Version]
  40. Paun, A.; Danska, J.S. Modulation of type 1 and type 2 diabetes risk by the intestinal microbiome. Pediatric Diabetes 2016, 17, 469–477. [Google Scholar] [CrossRef]
  41. Keying, X.; Gui, Q.; Yang, Y. Study on the changes in intestinal phylum firmicutes in elderly patients with type 2 diabetes. Chin. J. Geriatr. 2017, 36, 195–198. [Google Scholar]
  42. Koliada, A.; Syzenko, G.; Moseiko, V.; Budovska, L.; Puchkov, K.; Perederiy, V.; Gavalko, Y.; Dorofeyev, A.; Romanenko, M.; Tkach, S. Association between body mass index and Firmicutes/Bacteroidetes ratio in an adult Ukrainian population. BMC Microbiol. 2017, 17, 120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Anbalagan, R.; Srikanth, P.; Mani, M.; Barani, R.; Seshadri, K.G.; Janarthanan, R. Next generation sequencing of oral microbiota in Type 2 diabetes mellitus prior to and after neem stick usage and correlation with serum monocyte chemoattractant-1. Diabetes Res. Clin. Pract. 2017, 130, 204–210. [Google Scholar] [CrossRef] [PubMed]
  44. Xiao, E.; Mattos, M.; Vieira, G.H.A.; Chen, S.; Corrêa, J.D.; Wu, Y.; Albiero, M.L.; Bittinger, K.; Graves, D.T. Diabetes enhances IL-17 expression and alters the oral microbiome to increase its pathogenicity. Cell Host Microbe 2017, 22, 120–128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Barton, W.; Penney, N.C.; Cronin, O.; Garcia-Perez, I.; Molloy, M.G.; Holmes, E.; Shanahan, F.; Cotter, P.D.; O’Sullivan, O. The microbiome of professional athletes differs from that of more sedentary subjects in composition and particularly at the functional metabolic level. Gut 2018, 67, 625–633. [Google Scholar] [CrossRef]
  46. Hu, J.; Iragavarapu, S.; Nadkarni, G.N.; Huang, R.; Erazo, M.; Bao, X.; Verghese, D.; Coca, S.; Ahmed, M.K.; Peter, I. Location-specific oral microbiome possesses features associated with CKD. Kidney Int. Rep. 2018, 3, 193–204. [Google Scholar] [CrossRef] [Green Version]
  47. Zaura, E.; Keijser, B.J.; Huse, S.M.; Crielaard, W. Defining the healthy” core microbiome” of oral microbial communities. BMC Microbiol. 2009, 9, 259. [Google Scholar] [CrossRef] [Green Version]
  48. Goodson, J.; Groppo, D.; Halem, S.; Carpino, E. Is obesity an oral bacterial disease? J. Dent. Res. 2009, 88, 519–523. [Google Scholar] [CrossRef]
  49. Si, J.; Lee, C.; Ko, G. Oral microbiota: Microbial biomarkers of metabolic syndrome independent of host genetic factors. Front. Cell. Infect. Microbiol. 2017, 7, 516. [Google Scholar] [CrossRef] [Green Version]
Microorganisms 07 00645 g001 550
Figure 1. Spearman’s correlation analysis was applied for the pairwise analysis of plasma biochemical profile and demographic characteristics. The Benjamini–Hochberg false discovery rate (FDR) was used to perform multiple comparisons of the p-values. The color intensity shows the strength of correlation. Asterisks in each box indicate the corrected p-value; *** ≤ 0.001, ** ≤ 0.01 and * ≤ 0.05. Analyses were performed using the corrplot package in RStudio version 3.5.0.
Figure 1. Spearman’s correlation analysis was applied for the pairwise analysis of plasma biochemical profile and demographic characteristics. The Benjamini–Hochberg false discovery rate (FDR) was used to perform multiple comparisons of the p-values. The color intensity shows the strength of correlation. Asterisks in each box indicate the corrected p-value; *** ≤ 0.001, ** ≤ 0.01 and * ≤ 0.05. Analyses were performed using the corrplot package in RStudio version 3.5.0.
Microorganisms 07 00645 g001
Microorganisms 07 00645 g002 550
Figure 2. Alpha diversity analysis at 28,153 sequencing depth for observed_OTUs, faith_PD, and Shannon indexes.
Figure 2. Alpha diversity analysis at 28,153 sequencing depth for observed_OTUs, faith_PD, and Shannon indexes.
Microorganisms 07 00645 g002
Microorganisms 07 00645 g003 550
Figure 3. Beta diversity analysis on a 3D principal coordinates (PCoA) plot for weighted_unifrac, unweighted_unifrac and Bray Curtis indexes.
Figure 3. Beta diversity analysis on a 3D principal coordinates (PCoA) plot for weighted_unifrac, unweighted_unifrac and Bray Curtis indexes.
Microorganisms 07 00645 g003
Microorganisms 07 00645 g004 550
Figure 4. Taxonomic distribution of the microbiome at the phylum, family and species level. Height of each bar represents the relative proportion of that phylotype. Only microbes with a relative proportion of at least 0.5% are presented in the bar plots. The Kruskal–Wallis test shows no differences (p ≤ 0.05) between the study groups.
Figure 4. Taxonomic distribution of the microbiome at the phylum, family and species level. Height of each bar represents the relative proportion of that phylotype. Only microbes with a relative proportion of at least 0.5% are presented in the bar plots. The Kruskal–Wallis test shows no differences (p ≤ 0.05) between the study groups.
Microorganisms 07 00645 g004
Microorganisms 07 00645 g005a 550Microorganisms 07 00645 g005b 550
Figure 5. A bar chart of the estimated means of the Firmicutes/Bacteroidetes ratio and Fusobacteria in the control and obese groups after adjusting for triglycerides, insulin and HBA1c. Data are presented as the mean (95%CI).
Figure 5. A bar chart of the estimated means of the Firmicutes/Bacteroidetes ratio and Fusobacteria in the control and obese groups after adjusting for triglycerides, insulin and HBA1c. Data are presented as the mean (95%CI).
Microorganisms 07 00645 g005aMicroorganisms 07 00645 g005b
Microorganisms 07 00645 g006 550
Figure 6. Spearman’s correlation analysis between bacterial phylum and demographic/plasma biochemistry values. The color intensity indicates the strength of the correlation depicted as r-value. Analyses were performed using the corrplot package in RStudio version 3.5.0. Asterisks in each box indicate the corrected p-value; ** ≤ 0.01 and * ≤ 0.05.
Figure 6. Spearman’s correlation analysis between bacterial phylum and demographic/plasma biochemistry values. The color intensity indicates the strength of the correlation depicted as r-value. Analyses were performed using the corrplot package in RStudio version 3.5.0. Asterisks in each box indicate the corrected p-value; ** ≤ 0.01 and * ≤ 0.05.
Microorganisms 07 00645 g006
Table
Table 1. Demographic characteristics * and plasma biochemistry profile of the study cohort.
Table 1. Demographic characteristics * and plasma biochemistry profile of the study cohort.
Study ParameterStudy GroupCorrected
p-Value
ObeseControl
BMI35.65 ± 4.9222.73 ± 1.52< 0.001
Hemoglobin A1c (HbA1c) (%)5.68 ± 0.515.14 ± 0.76< 0.001
Fasting glucose (mmol/L)5.96 ± 0.424.80 ± 0.35< 0.001
Insulin (umol/L)16.56 ± 7.956.42 ± 1.86< 0.001
Homeostatic Model Assessment of Insulin Resistance (HOMA-IR)4.3 ± 1.71.3 ± 0.3< 0.001
Triglyceride (mmol/L)1.74 ± 1.121.18 ± 0.740.051
Total cholesterol (mmol/L)5.09 ± 0.975.32 ± 1.030.382
HDL cholesterol (mmol/L)1.11 ± 0.281.32 ± 0.430.772
LDL cholesterol (mmol/L)3.23 ± 0.893.47 ± 0.860.215
Thyroid-stimulating hormone (mIU/L)1.91 ± 1.381.95 ± 1.480.890
Triiodothyronine (pmol/L)4.21 ± 0.703.99 ± 0.640.213
Thyroxine (pmol/L)12.62 ± 1.4113.01 ± 1.480.208
Testosterone (nmol/L)13.38 ± 3.6722.74 ± 7.350.005
Sex hormone-binding globulin (nmol/L)26.20 ± 12.0137.00 ± 16.280.034
Estradiol (pmol/L)111.43 ± 42.2692.27 ± 36.890.030
Creatinine (µmol/L)72.05 ± 8.7277.53 ± 9.190.027
C-peptide (ng/mL)3.04 ± 0.961.51 ± 0.500.048
Aspartate aminotransferase (U/L)24.78 ± 9.8421.61 ± 8.740.120
Alanine aminotransferase (U/L)39.68 ± 18.9724.75 ± 13.860.004
Alkaline phosphatase (U/L)73.92 ± 14.3568.36 ± 12.830.111
Homocysteine (umol/L)9.65 ± 2.5410.53 ± 2.850.216
* Daily activity2.39 ± 1.613.00 ± 1.660.195
* Breastfed in infancy1.12 ± 0.331.00 ± 0.000.137
* Fast food consumption2.53 ± 1.423.22 ± 1.550.212
* Smoking1.19 ± 1.241.51 ± 1.360.343
* Antibiotic usage in last one year2.39 ± 1.041.49 ± 0.500.204
Data are the mean ± standard deviation. The Mann–Whitney U test was performed to compare mean differences between the groups. The Benjamini–Hochberg false discovery rate (FDR) was used to perform multiple comparisons of the p-values. * Data provenance: The demographic data were collected as per the questionnaire designed by Qatar Biobank. For statistical analysis, the data were categorically arranged from 1 to 6, where 1 is lowest or negative value and 6 is the highest possible value. Please refer to the supplementary file for data provenance.
Table
Table 2. Predictors of microbial taxa by stepwise linear regression after correcting for all potential confounders (BMI, lipids, glucose, and insulin).
Table 2. Predictors of microbial taxa by stepwise linear regression after correcting for all potential confounders (BMI, lipids, glucose, and insulin).
BacteriumPredictorAdjusted R SquareStd. Error of the Estimatep-Value
AcidobacteriaEstradiol0.0671.7 0.023
Estradiol0.131.7 0.005
HDL cholesterol0.024
FirmicutesEstradiol0.0979.70.008
Estradiol0.1419.40.002
HDL cholesterol0.046
Estradiol0.1869.20.003
HDL cholesterol0.028
Triiodothyronine0.043
FusobacteriaCreatinine0.1033.60.006
A stepwise linear regression analysis was performed on the bacterial taxa that showed significant association with different plasma biochemistry markers.

Share and Cite

MDPI and ACS Style

Sohail, M.U.; Elrayess, M.A.; Al Thani, A.A.; Al-Asmakh, M.; Yassine, H.M. Profiling the Oral Microbiome and Plasma Biochemistry of Obese Hyperglycemic Subjects in Qatar. Microorganisms 2019, 7, 645. https://doi.org/10.3390/microorganisms7120645

AMA Style

Sohail MU, Elrayess MA, Al Thani AA, Al-Asmakh M, Yassine HM. Profiling the Oral Microbiome and Plasma Biochemistry of Obese Hyperglycemic Subjects in Qatar. Microorganisms. 2019; 7(12):645. https://doi.org/10.3390/microorganisms7120645

Chicago/Turabian Style

Sohail, Muhammad U., Mohamed A. Elrayess, Asma A. Al Thani, Maha Al-Asmakh, and Hadi M. Yassine. 2019. "Profiling the Oral Microbiome and Plasma Biochemistry of Obese Hyperglycemic Subjects in Qatar" Microorganisms 7, no. 12: 645. https://doi.org/10.3390/microorganisms7120645

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