Fifty-one samples were sequenced and a total of 3,537,171 sequences were amplified, ranging from a minimum of 58,853 to a maximum of 74,594 sequences per sample, with a mean of 69,356 sequences per sample.
A total of 148 different bacterial species were detected. The most abundant species among all the samples was Porphyromonas_gingivalis (mean 17.82%), followed by Leptotrichia_buccalis (3.78%) and Neisseria_elongata (2.09%) (Fig. 1A). Data analysis of all the abundant species revealed a lower abundance of Tannerella_forsythia and Treponema_medium in the chewer without periodontitis group than in the other groups (P < 0.05)( Fig. 1B). Conversely, Haemophilus_parainfluenzae were more abundant in the chewer without periodontitis group (P < 0.05).
3.1 Species richness and diversity (alpha diversity)
Diversity analysis (Alpha diversity) can reflect the abundance and diversity of microbial communities. A statistically significant difference in PD Whole tree was observed (PD Whole tree, P=0.049) (Fig. 2). No significant differences were observed in Chao1, Observed Species,Simpsion and Shannon index among groups (Chao1, P=0.597; Simpsion, P=0.367, Observed Species, P=0.444; Shannon index, P=0.204) (Fig. 2).
3.2 PCoA and MRPP analysis (beta diversity)
To evaluate diversity difference, PCoA and MRPP analysis were analyzed. PCoA plots demonstrated clustering among groups (Fig. 3A). The unweighted UniFrac distance metric significantly differed (P=0.002) (Fig. 3B, Tab.1). Based on species and genus-level analysis, there were significant differences in distance between samples (species.bray, P=0.023; genus.bray, P= 0.04)(Fig.3CD, Tab.1).
3.3 Biomarker analysis ( LDA Effect Size)
To evaluate microbial contents among groups, linear discriminant analysis (LDA) was performed using LEfSe. From 946 biomarkers or higher taxa, 79 taxa were statistically significantly different among groups (P< 0.05) and had a linear discriminant analysis (LDA) score above 2 (Fig. 4). Taxa are ranked by the effect size in LEfSe(Fig. 4).
The most discriminative taxa in the samples of chewer with periodontitis group were Spirochaetae; in the samples of chewer without periodontitis group, the taxa were Treponema_denticola and Lentimicrobiaceae; in the samples of non-chewer with periodontitis group, the taxa were Bacilli and Lactobacillales; in the samples of non-chewer without periodontitis group, the taxa were Fretibacterium and Synergistetes.
Table 1 MRPP analysis