Skip to main content

The oral bacterial microbiota facilitates the stratification for ulcerative colitis patients with oral ulcers

Abstract

Background

Clinically, a large part of inflammatory bowel disease (IBD) patients is complicated by oral lesions. Although previous studies proved oral microbial dysbiosis in IBD patients, the bacterial community in the gastrointestinal (GI) tract of those IBD patients combined with oral ulcers has not been profiled yet.

Methods

In this study, we enrolled four groups of subjects, including healthy controls (CON), oral ulcer patients (OU), and ulcerative colitis patients with (UC_OU) and without (UC) oral ulcers. Bio-samples from three GI niches containing salivary, buccal, and fecal samples, were collected for 16S rRNA V3-V4 region sequencing. Bacterial abundance and related bio-functions were compared, and data showed that the fecal microbiota was more potent than salivary and buccal microbes in shaping the host immune system. ~ 22 UC and 10 UC_OU 5-aminosalicylate (5-ASA) routine treated patients were followed-up for six months; according to their treatment response (a decrease in the endoscopic Mayo score), they were further sub-grouped as responding and non-responding patients.

Results

We found those UC patients complicated with oral ulcers presented weaker treatment response, and three oral bacterial genera, i.e., Fusobacterium, Oribacterium, and Campylobacter, might be connected with treatment responding. Additionally, the salivary microbiome could be an indicator of treatment responding in 5-ASA routine treatment rather than buccal or fecal ones.

Conclusions

The fecal microbiota had a strong effect on the host’s immune indices, while the oral bacterial microbiota could help stratification for ulcerative colitis patients with oral ulcers. Additionally, the oral microbiota had the potential role in reflecting the treatment response of UC patients. Three oral bacteria genera (Fusobacterium, Oribacterium, and Campylobacter) might be involved in UC patients with oral ulcers lacking treatment responses, and monitoring oral microbiota may be meaningful in assessing the therapeutic response in UC patients.

Background

Ulcerative colitis (UC) is a chronic inflammatory disease mainly involving the colon; its pathogenesis is multifactorial, including genomic risks, immunological dysfunction, environmental factors, and gut microbial dysbiosis [1]. Increasing evidence has shown that microbial-host interaction is pivotal in homeostasis and pathogenesis. The gut microbiome educates the immune system during host development. Those germ-free animal models have immunological defects in the intestine, including smaller lymph nodes and Peyer’s patches decreased the number of the helper T cells [2]. Besides, gut microbial dysbiosis is also involved in the pathogenesis of various gut and even extra-intestinal diseases. The gut-liver [3], gut-brain [4], gut-lung [5], gut-bone marrow [6], and gut-skin axes [7] based etiologies have explained the underlying mechanisms of multi-organ diseases, such as metabolic associated fatty liver disease (MAFLD), autism spectrum disorder, chronic obstructive pulmonary disease (COPD), arthritis, and atopic dermatitis, etc.

Notably, ~ 2–34% of UC patients are complicated with oral manifestations, such as ulceration [8]. Interestingly, the oral microbiome has also been reported to be correlated with IBD [9,10,11], which indicates that the oral microbial community has the potential to trigger in situ lesions in IBD patients. Nevertheless, the microbial community in UC patients with oral ulcers has not been profiled.

To recover the microbial community in disorder, fecal microbiota transplantation (FMT) tends to be an efficient therapeutic strategy for some diseases [12]. Our previous studies have also reported that the therapeutic factors in 5-ASA routine treatment drive microbial alteration, whether in bacterial or fungal communities [13, 14]. These studies showed that the microbes-based clinical management is reasonable; however, it is unclear whether it would be meaningful for monitoring microbial profiles in discriminating disease activities and treatment efficacy. Our previous study showed that the Escherichia-Shigella richness in inflamed mucosa positively correlates with the UC activity [13]. Also, microbial contents have been used to predict activity in Crohn’s disease [15], irritable bowel syndrome [16], and other diseases [17, 18].

Based on these studies, to explore the bacterial community in the GI tract of IBD patients combined with oral ulcers, we analyzed the microbial profiles by 16S rRNA sequencing at three GI niches, including salivary, buccal, and fecal samples of UC patients with or without oral ulcers. We found that complicating with oral ulcers made UC patients a weak treatment response; and three oral bacterial genera, including Fusobacterium, Oribacterium, and Campylobacter, might be involved in the process. These results highlighted the significance of stratification of UC patients with oral diseases. Additionally, monitoring oral microbiota may be meaningful in assessing the therapeutic response in UC patients.

Methods

Study design and sample collection

This study was approved by the Institutional Medical Ethics Review Board of Peking University People’s Hospital. All subjects were enrolled at Peking University People’s Hospital from January 2017 to January 2020. To test the bacterial difference in UC patients with or without oral ulcers, we enrolled four groups of subjects, including healthy controls (CON), patients with only oral ulcers (OU), UC patients without oral ulcers (UC), and UC patients with oral ulcers (UC_OU) (Fig. 1A and Table 1). The healthy subjects were recruited from volunteers taking 5-ASA routine health examinations in Peking University People’s Hospital. The exclusion criteria of healthy participants included periodontitis, OU, UC, Crohn’s disease, Behçet’s disease, and other chronic or acute inflammatory statuses with endoscopically confirmation. The diagnosis of UC was established according to the World Gastroenterology Organization Global Guidelines [19]. Patients in the OU group were confirmed to suffer only oral ulcers without other oral and systemic diseases. In addition, those UC patients enrolled in the UC_OU group were historically recurrent oral ulcers, and suffering oral ulcers during sampling. All subjects were informed to avoid taking antibiotics, yogurt, and other probiotics within 4 weeks; hereafter, all subjects’ written informed consents were obtained before sampling. In addition, we documented the demographic and clinical information of all subjects.

Fig. 1
figure 1

Bacterial profiles at different gastrointestinal tract niches. A Study design. Three types of samples from differed GI niches were collected in this study, including salivary, buccal, and fecal samples. B Alpha diversity indices of the microbiota, including the richness, Simpson’s, Shannon’s, and Chao1 indices. Horizontal bars within boxes represent medians. The tops and bottoms of the boxes represent the 75th and 25th percentiles, respectively. The upper and lower whiskers cover 1.5 × the interquartile range from the upper and lower edges of the box, respectively. P-values were obtained using the one-way ANOVA test (comparisons among four groups). C The constrained principal coordinate analysis based on the Bray–Curtis distance. The R software (v 4.0.1) with the vegan (v 2.5–7) package were used and P-values were obtained using permutational multivariate analysis of variance (PERMANOVA). D and E The upset plot shows the bacterial family (D) and genus (E) count in each GI niche. F and G Relative abundance of the top 20 bacterial families (F) and genera (G). Visualization was performed using Circos (http://circos.ca/). The right circle in the outer part shows the groups and relative proportions of bacterial species. The left outer ring and inner bands indicate the relative proportions (%) of bacterial genera in the different groups. The left inner circle represents the relative abundances of all bacteria. H The ternary plot shows the distribution of the specific bacterial genus in each GI niche. The point color represents the phylum classification of the bacterial genus. The point size presents the mean percentage of a specific genus

Table 1 Demographic and clinical profiles of patients from different groups and healthy controls

Sample collection and DNA extraction

Bio-samples from three different spatial locations in the digestive tract, including salivary, buccal, and fecal samples, were collected for the following test. Briefly, The patient rinses his mouth with clean water before sampling and collects the unstimulated saliva into a sterile DNase-free and RNase-free centrifuge tube for about 5–8 min. After that, buccal mucosa was swabbed 50 times on the left and right side walls of mouth with a sterile cotton swab without touching teeth and gums. Fecal samples were collected in a Stool Collection Tube with Stool Stabilizer (German, Stratec Molecular) and stored with a – 80 °C freezer as in previous reports[13, 14]. According to the manufacturer’s instructions, the bacterial genomic DNA was extracted with the PSP Spin Stool DNA Kit (Stratec Molecular). Hereinafter, the bacterial 16S V3-V4 rRNA was amplified using 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R primers (5′-GGACTACHVGGGTWTCTAAT-3′) following the protocol described in our previous study [20]. In brief, a 20-μL reaction system containing the FastPfu Polymerase (TransGen Biotech Co., Beijing, China) was used to amplify the 16S rRNA [20]. All PCR products were purified and added with sample-specific barcodes.

16S rRNA sequencing and bioinformatic analysis

The Illumina Miseq platform (Illumina, San Diego) was used for 16S rRNA sequencing. The raw data were produced with the Vsearch v2.8.1 [21] and Usearch v11 (bit 64) software [22]. Vsearch was performed for merging the original data, quality control, primer, and barcode sequences excision. In total, 2,317,709 sequences were removed, leaving 55,571,009 sequences for further analysis. Redundant sequences were then filtered, and 15,157,810 unique sequences were retained. Vsearch was used for discarding low-occurrence sequences. 33,920,292 sequences were removed, and 12,254 amplicon sequences were left. Remove 4951 chimers in the 12,254 sequences and 7,303 high-quality sequences were acquired. Subsequently, chimeras filtering was conducted with the amplicon sequence variants (ASVs) method, followed by sequencing error checking using Usearch as in our previous report [23, 24]. As a result, 4,121 high-qualified ASVs were obtained and aligned with Ribosomal Database Project’s (RDP) training Set 16 [25].

Functional prediction with the PICRUSt2

The Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) software was used for functional prediction [26]. In investigating the etiology of IBD, we focused on the immune-mediated pathways in these studies [27].

The stratification of UC patients

The clinical parameters, including complement C3 (C3), complement C4 (C4), C-reactive protein (CRP), erythrocyte sedimentation rate (ESR), immunoglobulin M (IgM), immunoglobulin A (IgA), immunoglobulin G (IgG), monocyte (Mo), lymphocyte (LY), monocytes’ percentage (MO_P), white blood cell (WBC), lymphocyte’ percentage (LY_P), eosinophil (EOS), Hemoglobin(Hb), platelet (PLT), basophils(BAS), basophils’ percentage (BAS_P), neutrophils (NE), neutrophils’ percentage (NE_P), were recorded for all UC patients at the sampling baseline.

Except for a small group of UC patients losing touch, ~ 22 UC patients and 10 UC_OU patients underwent a six-month following-up. In comparing clinical parameters, especially in the Mayo clinical score, the patients were stratified into responding and non-responding groups (Fig. 1A and Table 2).

Table 2 Demographic and clinical profiles of subjects followed-up

Statistical analysis and data visualization

The STAMP software (v 2.1.3) [28] and R software (v 4.0.1; R Foundation for Statistical Computing, Vienna, Austria) with the ggplot2 (v 3.3.2) package were used for data visualization [29]. Permutational multivariate analysis of variance (PERMANOVA, Adonis test of vegan v 2.5–6) was performed for statistical analysis of beta diversity. The independent sample t-test and nonparametric Mann–Whitney U test compared two groups. One-way analysis of variance (ANOVA) and Kruskal–Wallis H nonparametric tests compared the three groups. Spearman’s correlation analysis was performed, and the P-value was corrected with the false discovery rate (FDR). Significant correlations were visualized using the pheatmap package (v1.0.12). A P-value or FDR ≤ 0.05 was considered statistically significant.

Results

The microbial community composition of salivary, buccal, and fecal samples

To investigate the characteristic microbial profiles in the GI tract, we conducted the 16S rRNA sequencing and compared the bacterial diversity in different GI locations first. The alpha diversity analysis revealed that richness and chao1 indices were gradually higher in fecal samples than in salivary and buccal samples (Fig. 1B, Additional file 5: Table S1), indicating a higher microbial community richness. In addition, the Simpson index was lower, and the Shannon index was higher in salivary and fecal samples than in buccal samples (Fig. 1B, Additional file 5: Table S1), which suggested higher evenness of microbial community composition in salivary and fecal samples. Furthermore, beta diversity analysis with constrained principal coordinates analysis (CPCoA) showed an adjacency between salivary and buccal samples while a separation between fecal and oral (salivary and buccal) samples in Bray–Curtis distance. Nonetheless, there were significant differences among the three groups (Adonis P < 0.05, Fig. 1C, Additional file 5: Table S2). These data suggested that salivary and buccal samples share similar microbial community composition to a certain degree, whereas fecal samples did not.

Further analysis showed that thirty-eight bacterial families were shared in three different niches and oral samples shared seven other families (Fig. 1D and Additional file 5: Table S3). In addition, there were twenty-nine special bacterial families in fecal samples (e.g., Acholeplasmataceae, Acidaminococcaceae, and Brucellaceae), four in buccal samples (Hicobacteraceae, Nocardiaceae, Nocardiodaceae, and Staphylococcaceae) and one unique bacterial family in salivary samples (Bacillaceae_1) (Fig. 1D and Additional file 5: Table S3). At the genus level, data showed that forty-two bacterial genera were shared among the three sites, and seventy-three in oral samples. As to unique genera, there were seventy-five bacteria specific in fecal samples (e.g., Akkermansia, Bilophila, and Citrobacter, etc.), seven in buccal samples (e.g., Aeromicrobium, Helicobacter, and Nocardia, etc.), while four in salivary samples (including Alloscardovia, Bacillus, Mobiluncus, and Mogibacterium) (Fig. 1E and Additional file 5: Table S3). These data indicated a higher bacterial richness in the feces, and some bacterial species were shared in oral and fecal samples.

We further analyzed the twenty most abundant bacterial families and genera and found that Streptococcaceae and Neisseriaceae mainly colonize the oral cavity at the family level, while Lachnospiraceae, Bacteroidaceae, Ruminococcaceae, and Enterobacteriaceae colonized with a high richness in feces (Fig. 1F and Additional file 5: Table S5). At the genus level, Streptococcus, Neisseria, Haemophilus, and Rothia occupied the oral cavity with high abundance; while Faecalibacterium, Bacteroides, Bifidobacterium, Escherichia/Shigella, and Lachnospiraceae_incertae_sedis were highly abundant in feces (Fig. 1G, H and Additional file 5: Table S5-6).

These results indicate that the bacteria in salivary and buccal samples are highly similar at the family and genus levels; they share many common and similar dominant bacteria but fewer unique bacteria. In contrast, there are many specific bacterial families and genera in feces but not in oral samples, indicating that it is usually difficult for most bacteria in the oral cavity to colonize in the gut. Thus, we speculate that this may be related to the different microenvironments in the oral cavity and gut.

The salivary bacterial microbes in UC patients complicated by oral ulcers

To uncover the microbial features in salivary samples of UC_OU patients, we further profiled the salivary bacterial community in all subjects. Compared with the CON group, the richness increased in the UC group, and the diversity and evenness elevated in the UC_OU group (Fig. 2A and Additional file 5: Table S1). Moreover, the beta diversity analysis revealed that the salivary microbial community composition was significantly different between the groups of patients and the CON (Adonis P = 0.046 for OU, Adonis P = 0.002 for UC, P = 0.004 for UC_OU vs. CON, respectively, Fig. 2B, S1B, and Additional file 5: Table S2). Furthermore, the analysis of the twenty most abundant bacterial families and genera showed that there were remarkable differences in abundance between the groups of patients and the CON at the family and genus level (Fig. 2C and Additional file 5: Table S5-6). Further investigation of UC_OU and OU manifested that some bacteria decreased (e.g., Blautia, Clostridium_XIII, and Faecalibacterium), and some increased (e.g., Abiotrophia) (Fig. 2D and Additional file 5: Table S7). Surprisingly, there were up to fifty different genera between the UC_OU and UC groups, including the increasing Klebsiella, Arthrobacter, Barnesiella, and the decreasing Blautia, Clostridium_XIII, and Faecalibacterium (Fig. 2E and Additional file 5: Table S7). We also compared the groups of patients with the CON respectively, and found that some genera (e.g., Arthrobacter, Barnesiella, Alistipes) increased and some (e.g., Mogibacterium) decreased consistently (Additional file 1: Figure S1C-E and Additional file 5: Table S7).

Fig. 2
figure 2

The bacterial profile in salivary samples of UC patients with or without oral ulcers. A Alpha diversity indices of the microbiota, including the richness, Simpson’s, Shannon’s, and Chao1 indices. Horizontal bars within boxes represent medians. The tops and bottoms of the boxes represent the 75th and 25th percentiles, respectively. The upper and lower whiskers cover 1.5 × the interquartile range from the upper and lower edges of the box, respectively. P-values were obtained using the one-way ANOVA test (comparisons among four groups). B The constrained principal coordinate analysis based on the Bray–Curtis distance. The R software (v 4.0.1) with the vegan (v 2.5–7) package were used, and P-values were obtained using permutational multivariate analysis of variance (PERMANOVA). C Relative abundance of the top 20 bacterial families (the left panel) and genera (the right panel). Visualization was performed using Circos (http://circos.ca/). The right circle in the outer part shows the groups and relative proportions of bacterial species. The left outer circle and inner bands show the relative proportions (%) of bacterial genera in the different groups. The left inner circle represents the relative abundances of all bacteria. D and E Comparative analysis of bacterial genus abundance between two groups (D UC_OU vs. OU; E UC_OU vs. UC). The EdgeR package was used for comparative analysis. The difference between the two groups is shown as a Manhattan diagram. Point shape indicates the genus enriched, depleted, or not significant in the former group compared with the latter. Point size indicates the counts of a specific genus. CPM, count per million. F and G Comparative analysis of bacterial function between two groups (F UC_OU vs. OU; G UC_OU vs. UC). Phylogenetic Investigation of Communities annotated the pathway information by Reconstruction of Unobserved States (PICRUSt2) software by referring to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The STAMP software was used for data visualization. CON, healthy controls; OU, patients with only oral ulcers; UC, UC patients without oral ulcers; UC_OU, UC patients with oral ulcers; ns, not significant; *P-value < 0.05; **P-value < 0.01

Furthermore, the PICRUSt2 analysis showed that in the group of UC_OU, the immune-related signaling pathways (including antigen processing and presentation, IL-17 signaling, and Th17 cell differentiation, etc.) were active compared with the OU group. The anti-inflammatory pathway (e.g., thiamine metabolism [30]) was restrained compared with the UC group (Fig. 2F, G and Additional file 5: Table S8), indicating a higher immune activation and inflammatory states in UC_OU patients.

In general, the increment of Arthrobacte and Barnesiella and the reduction of Blautia, Clostridium_XIII, and Faecalibacterium in UC_OU patients may be essential factors in mediating the occurrence and development of the disease by regulating immune responses and inflammatory pathways.

The buccal bacterial community in UC patients complicated by oral ulcers

There were no significant differences in alpha diversity among the buccal samples of the CON, OU, UC, and UC_OU. The beta diversity analysis showed that there were significant differences when the OU/UC, but not the UC_OU group, compared with the CON group (Adonis P = 0.008 for OU, P = 0.001 for UC and P = 0.102 for UC_OU vs. CON, respectively, Fig. 3B, S2A, and Additional file 5: Table S2).

Fig. 3
figure 3

The bacterial profile in buccal samples of UC patients with or without oral ulcers. A Alpha diversity indices of the microbiota, including the richness, Simpson’s, Shannon’s, and Chao1 indices. Horizontal bars within boxes represent medians. The tops and bottoms of the boxes represent the 75th and 25th percentiles, respectively. The upper and lower whiskers cover 1.5 × the interquartile range from the upper and lower edges of the box, respectively. P-values were obtained using the one-way ANOVA test (comparisons among four groups). B The constrained principal coordinate analysis based on the Bray–Curtis distance. The R software (v 4.0.1) with the vegan (v 2.5–7) package were used, and P-values were obtained using permutational multivariate analysis of variance (PERMANOVA). C Relative abundance of the top 20 bacterial families (the left panel) and genera (the right panel). Visualization was performed using Circos (http://circos.ca/). The right circle in the outer part shows the groups and relative proportions of bacterial species. The left outer circle and inner bands show the relative proportions (%) of bacterial genera in the different groups. The left inner circle represents the relative abundances of all bacteria. D and E, Comparative analysis of bacterial genus abundance between two groups (D UC_OU vs. OU; E UC_OU vs. UC). The EdgeR package was used for comparative analysis. The difference between the two groups is shown as a Manhattan diagram. Point shape indicates the genus enriched, depleted, or not significant in the former group compared with the latter. Point size indicates the counts of a specific genus. CPM, count per million. F and G, Comparative analysis of bacterial function between two groups (F UC_OU vs. OU; G UC_OU vs. UC). Phylogenetic Investigation of Communities annotated the pathway information by Reconstruction of Unobserved States (PICRUSt2) software by referring to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The STAMP software was used for data visualization. CON, healthy controls; OU, patients with only oral ulcers; UC, UC patients without oral ulcers; UC_OU, UC patients with oral ulcers; ns, not significant; *P-value < 0.05; **P-value < 0.01

We further investigated the buccal bacterial composition in different groups of subjects. At the family level, the top twenty bacteria in abundance levels differed not so markedly among all groups of subjects (Fig. 3C and Additional file 5: Table S5). While at the genus level, the abundance of some bacteria, such as Neisseria and Rothia, declined slightly in the UC_OU group. As for Neisseria, a previous study found a falling abundance of which at the inflamed site of UC patients compared with the corresponding area of non-IBD controls [31]; in addition, the abundance of Actinobacillus and Fusobacterium varied primarily among these groups (Fig. 3C and Additional file 5: Table S6).

Using further differential analysis, we found that, compared with the CON group, the variation of bacteria was considerably conspicuous among OU, UC, and UC_OU group at the genus level; for instance, the abundance of Barnesiella, Alistipes, and Rhodopseudomonas ascended, while Actinobacillus descended (Additional file 2: Figure S2B-D and Additional file 5: Table S9). Then, we compared oral bacterial richness between UC_OU and OU groups and found very few differential bacteria, which manifested the increase of Abiotrophia, Cardiobacterium, and Klebsiella (Fig. 3D and Additional file 5: Table S9). Interestingly, Abiotrophia defectiva was related to pro-inflammatory response in the oral cavity [32]; moreover, Klebsiella pneumoniae was reported to aggravate chronic intestinal inflammation by destructing the intestinal epithelial barrier [33]. In contrast with UC patients, there was a marked increase in the abundance of Abiotrophia, Cardiobacterium, and Klebsiella in the UC_OU patients; on the contrary, the abundance of Actinobacillus decreased notably (Fig. 3E and Additional file 5: Table S9). Hence, we inferred that these bacteria probably contributed to oral ulcers in the UC_OU patients.

A previous study found that decreased sphingolipids correlated with gut inflammation in IBD subjects [34]. Our PICRUSt2 analysis also observed suppressed sphingolipid signaling pathway in the UC_OU patients compared with OU patients. In addition, there was an enrichment of protein processing in the endoplasmic reticulum in the UC_OU patients (vs. OU, Fig. 3F and Additional file 5: Table S10). It was found multiple immune cells could activate that unfolded protein response at distinct levels [35]. Furthermore, the differences in signal pathways between UC and UC_OU patients showed that the expression level of some signal pathways declined, such as dorso-ventral axis formation and glycosphingolipid biosynthesis (Fig. 3G and Additional file 5: Table S10); among which glycosphingolipid was reported to correlate with regulating immune signaling with facilitating bacterial entering host cells [36]. The comparison results between OU and CON patients unveiled that some pathways diminished in the OU patients, such as pyrimidine metabolism, RNA polymerase, protein phosphatases and associated proteins, and secondary bile acid biosynthesis (Additional file 2: Figure S2E and Additional file 5: Table S10), which were all associated with immunity [37,38,39]. Moreover, there was an enriched apoptosis pathway in the UC group compared with the CON one (Additional file 2: Figure S2F and Additional file 5: Table S10), which was paralleled with a previous review that showed the pathogenic function of caspase-mediated intestinal epithelial cell apoptosis in the IBD [40].

In brief, buccal bacterial features in UC_OU subjects differed from UC alone or OU. and the abundance of some bacteria (such as Abiotrophia, and Klebsiella) in the buccal mucosa may help to distinguish between UC_OU subjects and UC alone or OU. Furthermore, the alterations of some signaling pathways related to immune cells or processes implied that immune factors might participate in the occurrence and development of oral ulcers in UC_OU patients.

The fecal bacterial microbiota in UC_OU patients

We further analyzed the fecal bacterial composition of all subjects. We found that, compared with the CON, the alpha diversity of fecal bacteria in UC_OU patients significantly decreased, while the OU or UC patients did not (Fig. 4A and Additional file 5: Table S1); beta diversity analysis showed a marked difference in the fecal microbiota composition between patients and CON (Adonis P = 0.011 for OU, 0.002 for UC, and 0.003 for UC_OU vs. CON, respectively; Fig. 4B, S3A and Additional file 5: Table S2). In the top 20 bacterial families, Enterobacteriaceae increased in the groups of patients, while the abundances of Prevotellaceae in the UC and Veillonellaceae in the OU were decreased than in the CON group (Fig. 4C and Additional file 5: Table S5). At the genus level, compared with the CON, Prevotella and Roseburia declined in the UC group, while Escherichia/Shigella raised in the UC and UC_OU group (Fig. 4C and Additional file 5: Table S6).

Fig. 4
figure 4

The bacterial profile in fecal samples of UC patients with or without oral ulcers. A Alpha diversity indices of the microbiota, including the richness, Simpson’s, Shannon’s, and Chao1 indices. Horizontal bars within boxes represent medians. The tops and bottoms of the boxes represent the 75th and 25th percentiles, respectively. The upper and lower whiskers cover 1.5 × the interquartile range from the upper and lower edges of the box, respectively. P-values were obtained using the one-way ANOVA test (comparisons among four groups). B The constrained principal coordinate analysis based on the Bray–Curtis distance. The R software (v 4.0.1) with the vegan (v 2.5–7) package were used, and P-values were obtained using permutational multivariate analysis of variance (PERMANOVA). C, Relative abundance of the top 20 bacterial families (the left panel) and genera (the right panel). Visualization was performed using Circos (http://circos.ca/). The right circle in the outer part shows the groups and relative proportions of bacterial species. The left outer circle and inner bands show the relative proportions (%) of bacterial genera in the different groups. The left inner circle represents the relative abundances of all bacteria. D and E, Comparative analysis of bacterial genus abundance between two groups (D UC_OU vs. OU; E UC_OU vs. UC). The EdgeR package was used for comparative analysis. The difference between the two groups is shown as a Manhattan diagram. Point shape indicates the genus enriched, depleted, or not significant in the former group compared with the latter. Point size indicates the counts of a specific genus. CPM, count per million. F and G, Comparative analysis of bacterial function between two groups [F UC_OU vs. OU; G UC_OU vs. UC]. Phylogenetic Investigation of Communities annotated the pathway information by Reconstruction of Unobserved States (PICRUSt2) software by referring to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The STAMP software was used for data visualization. CON, healthy controls; OU, patients with only oral ulcers; UC, UC patients without oral ulcers; UC_OU, UC patients with oral ulcers; ns, not significant; *P-value < 0.05; **P-value < 0.01

We further performed a comparative analysis at the genus level. The results showed that, compared with the CON group, some bacteria (e.g., Arthrobacter, Knoellia, Bacillus, Peptostreptococcus, etc.) increased in the groups of patients (OU/UC/UC_OU) and some decreased (e.g., Lautropia) (Additional file 3: Figure S3B-D and Additional file 5: Table S11). Interestingly, Bacillus and Peptostreptococcus have been reported to correlate with severe infections [41] and IBD progression [42]. In addition, no differential genera could be found between UC_OU and OU groups (Fig. 4D); however, Alloprevotella, Granulicatella, Lactobacillus, and Holdemanella ascended strikingly in the UC_OU group compared with the UC group (Fig. 4E and Additional file 5: Table S11).

We then had the PICRUSt2 analysis. The results showed that some pathways, such as aminobenzoate degradation, xylene degradation, and bile secretion, lessened in the UC_OU patients compared with OU patients (Fig. 4F and Additional file 5: Table S12), among which xylene can exacerbate allergic inflammation [43]. Furthermore, the comparison between UC_OU patients and UC patients showed that metabolism-related pathways (including aminobenzoate degradation, platinum drug resistance, styrene degradation, etc.) were repressed (Fig. 4G and Additional file 5: Table S12). In addition, nucleotide metabolism and naphthalene degradation pathways in OU patients, and amino acid metabolism and fatty acid degradation in UC patients were activated compared with the CON (Additional file 3: Figure S3E-G and Additional file 5: Table S12).

To conclude, the fecal flora characteristics in the UC_OU group were distinguished from UC alone, and Alloprevotella, Granulicatella, Lactobacillus, and Holdemanella may contribute to differentiating between UC_OU and UC alone. In addition, signal pathways related to metabolization and immunity can involve the pathogenic process in the UC_OU.

The correlation of GI spatial microbiome with clinical, immunological parameters

To further assess whether the disease activities and clinical parameters correlated with the spatial microbiome alteration, we collected twenty clinical parameters responding to inflammatory and immunological statuses, including C-reaction protein (CRP), ESR, the percentage of monocytes (MO_P), etc. We performed a redundancy analysis (RDA) between bacterial beta diversity and clinical parameters and Spearman’s correlation between alpha diversity and clinical parameters (Fig. 5A–F and Additional file 5: Table S13). The inflammatory indices were statistically correlated with the salivary microbiome in alpha diversity and beta diversity (Fig. 5A). Of note, the salivary microbial alpha diversity correlated with the lymphocyte’s (rho = 0.281, FDR = 0.022 for chao1, and rho = 0.282, FDR = 0.022 for richness, respectively), and neutrophil’s percentage (rho = − 0.372, FDR = 0.002 for chao1, and rho = − 0.372, FDR = 0.002 for richness, respectively; Fig. 5B and Additional file 5: Table S13). The buccal bacterial Shannon index correlated with subjects’ inflammatory markers, such as C4 (rho = 0.567, FDR = 0.002) and CRP (rho = 0.399, FDR = 0.039) (Fig. 5C and Additional file 5: Table S13). Like the salivary microbiome, buccal bacterial beta diversity showed negative results correlating with clinical parameters (Figs. 5C, D and Additional file 5: Table S13). The fecal microbiota showed a significant correlation between beta diversity and immunological indices, such as the monocytes’ (pseudo F = 2.218, P = 0.039) and basophils’ percentage (pseudo F = 1.984, P = 0.035), rather than inflammatory markers (Fig. 5E, F, and Additional file 5: Table S13). Based on these results, we inferred that the fecal bacteria play a more important role in shaping the host immune system.

Fig. 5
figure 5

The correlation between clinical parameters (inflammatory and immunological) and bacterial diversity (alpha and beta). A and B The correlation between inflammatory (A) or immunological (B) indices and salivary bacterial diversity. The length of the black bar shows the R square value of the indices-beta diversity correlation. The results were acquired by the redundancy analysis (RDA) using the R software (v 4.0.1) with the vegan (v 2.5–7) package. The color box shows the rho value of indices-alpha diversity correlation. The results were analyzed from Spearman’s correlation using the R software (v 4.0.1) with the vegan (v 2.5–7) package. The star symbols behind the left words show the P-values acquired from the RDA. The box’s star symbols show the adjusted P-values (false discovery rate, FDR) received from the Spearman’s correlation. (same methodology as in C–F). C and D The correlation between inflammatory (C) or immunological (D) indices and buccal bacterial diversity. E and F The correlation between inflammatory (E) or immunological (F) indices and buccal bacterial diversity

Treatment responses in UC patients

We then analyzed the treatment response for UC patients with a six-month following-up. Thirty-two participants were followed-up, including twenty-two UC and ten UC_OU patients, for which we documented the Mayo clinic score to differentiate patients with or without treatment response. There were no significant difference of some clinical parameters (Mayo score, ESR, CRP, IgA, IgG, IgM, C3, and C4) before and after treatments in UC and UC_DU groups (Additional file 4: Figure S4). Interestingly, only one out of ten UC_OU patients responded to 5-ASA routine treatment, significantly lower than that in UC patients (Chi-squared test, X2 = 19.09, P < 0.001, Fig. 6A and Additional file 5: Table S14). We further analyzed the GI spatial bacterial profiles in UC patients with or without treatment responses. Only one subject had responded to the treatment in UC_OU patients, so this group of patients was excluded from this part of the investigation. Notably, there were eighteen genera with a remarkable difference in these two subgroups of patients in the salivary bacteria, including Prevotella, Alloprevotella, Fusobacteria, Oribacterium, Campylobacter, and Rothia, etc. (Fig. 6B and Additional file 5: Table S14). Compared with salivary samples, the buccal biopsies showed fewer bacterial contents, with a significant difference between UC patients with and without treatment responses (Fig. 6C and Additional file 5: Table S14). It could be found that there were three same genera also enriched in the non-responding UC patients’ oral mucosae, which were Fusobacterium, Oribacterium, and Campylobacter, respectively. Additionally, only one content, i.e., Blautia, was represented in fecal microbiota in non-responding patients (Fig. 6D and Additional file 5: Table S14). Based on these data, we deduced that the richness of Fusobacterium, Oribacterium, and Campylobacter might be involved in non-response; the salivary microbiome also had a potential for indicating treatment response in UC patients.

Fig. 6
figure 6

Microbial differentiation between UC patients with and without treatment response. A The treatment response rate in UC patients with or without oral ulcers. The Chi-squared test was performed to test the response difference between the two groups. UC, UC patients without oral ulcers; UC_OU, UC patients with oral ulcers. B–D The comparative analysis of bacterial genus abundance of the samples from salivary (B), buccal (C), and fecal (D) niches. The STAMP software was used for data visualization. E The graphic summary

Discussion

Previous studies indicate up to 34% of UC patients are accompanied by oral ulcers [8], while the microbial community composition of oral bacteria in patients with or without oral ulcers and the relationship between oral microbiota and IBD prognosis are still unclear.

This study found that each GI region had a different bacterial community, which consisting with previous study reported by Vasapolli [44]. Salivary samples and buccal samples shared more common and less unique bacteria, while fecal samples were richer in bacterial diversity and individual genera. These may be related to the distinct microenvironments [45], gastrointestinal acid [46], bile acid chemical barriers [47], intestinal microbial colonization resistance [48], host pattern recognition receptors [49], and loss of oxygen [50]. Moreover, compared with oral bacterial community, the fecal one was closely correlated with patient clinical parameters, especially the proportion of immune cells rather than inflammation-related indicators, suggesting that the fecal microbiota may modulate the function of the immune cells and participate in the pathogenesis of UC_OU patients. The abundance and part of the fecal bacteria may partially account for this phenomenon.

The variation in bacterial richness at different GI niches facilitates to differentiating UC_OU individuals from UC patients. Previous research has been proved that the Klebsiella isolated from the saliva of IBD patients colonizes in the intestine ectopically to elicit colitis development [51] and Klebsiella-reactive Th17 may migrate to the gut to exacerbate intestinal inflammation [52]. The Blautia and Faecalibacterium were treated as the functional genera with potential probiotic properties, for their metabolic regulation and butyrate-producing capacity [53, 54]. In our study, compared with UC patients, the Klebsiella, Arthrobacter, Barnesiella were enriched, and the Blautia and Faecalibacterium were depleted in the UC_OU patients’ oral samples, which indicates a more robust state of inflammation and immune activation in UC_OU patients.

The PICRUSt2 analysis indicated that the UC_OU group differs from the UC group and may be a new subclass of UC patients. Thiamine (also known as vitamin B1) is synthesized by bacteria and fungi and regulates the immune system by activating immune cells and proteins [55]. It is also an anti-inflammatory factor regulating inflammatory agent expression and preventing recurring inflammation [30]. Interestingly, Klaassen et al. have found that Crohn’s Disease exacerbations were associated with a decrease in microbial genes involved in the biosynthesis of the anti-inflammatory mediators, including thiamine [56]. Similarly, in our study, the reduction of thiamine metabolism in the salivary samples of UC_OU patients indicate a higher pro-inflammatory state than in UC patients. Furthermore, previous research found that sphingolipids are significantly decreased in IBD subjects and negatively correlated with gut inflammation [34]. The decreased glycosphingolipid biosynthesis in UC_OU patients’ buccal samples may imply more severe intestinal inflammation compared with the UC patients. Moreover, there was a positive correlation between styrene levels and disease activity in UC [57], and the suppressed styrene degradation pathway of fecal samples suggested higher styrene levels in UC_OU patients compared with UC. The functional analysis of the microbiota at the three sites suggests that the UC_OU group is different from the UC and may have a more severe inflammatory status.

We also compare the bacterial profiles at different GI spatial niches of treatment-non-responding UC patients. The genera, including Fusobacterium, Oribacterium, and Campylobacter, were enriched in the oral cavity of non-responding UC patients. Of note, these bacteria were also represented in patients with systemic lupus erythematosus (SLE) [58] and were involved in the intestinal inflammation of patients with distal gastrectomy [59]. A previous study reported that a Fusobacterium species, i.e., F. nucleatum, in gut mucosa triggered gut inflammation and positively correlated with host IBD status [60, 61]. Campylobacter spp. was reported to recruit neutrophils to lead to gastroenteritis [62,63,64]. These data indicate that these oral bacteria might be related to therapeutic failure in UC patients. Of note, the salivary microbiome presented a markable difference in UC patients without response rather than the buccal or fecal microbiome. Additionally, a previous report confirmed that the oral microbes could recover to their initial state from collapsing in a few days and have sufficient robustness to serve as biomarkers [65, 66]. These results highlight the potential role of salivary microbiota as an indicator of IBD treatment efficiency.

Conclusion

In conclusion, we analyzed the microbial community composition of the salivary, buccal and fecal samples of UC patients with or without oral ulcers by 16S rRNA sequencing and their correlation with clinical indicators. We found that UC patients with oral ulcers lacked treatment responses, and three oral bacterial genera might be involved. The fecal microbiota had the most significant impact on the host immune indices, and the salivary microbiota had the potential role in reflecting the treatment response of UC patients and is promising to be an indicator for UC treatment efficiency. Our data throw light on the significance of stratification of UC patients with oral diseases, and monitoring oral microbiota may be meaningful in assessing the therapeutic response in UC patients (Fig. 6E).

Availability of data and materials

The 16S rRNA-seq data have been deposited in the Genome Sequence Archive in National Genomics Data Center, China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences (GSA: CRA006292) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa/. The data supporting the findings of this study are available within the article and its Supplementary Information files.

Abbreviations

C3:

Complement C3

C4:

Complement C4

CRP:

C-reactive protein

ESR:

Erythrocyte sedimentation rate

IgM:

Immunoglobulin M

IgA:

Immunoglobulin A

IgG:

Immunoglobulin G

MO:

Monocyte

LY:

Lymphocyte

MO_P:

Monocytes’ percentage

WBC:

White blood cell

LY_P:

Lymphocyte’ percentage

EOS:

Eosinophil

Hb:

Hemoglobin

PLT:

Platelet

BAS:

Basophils

BAS_P:

Basophils’ percentage

NE:

Neutrophils

NE_P:

Neutrophils’ percentage

ASV:

Amplicon sequence variants

COPD:

Chronic obstructive pulmonary disease

CPCoA:

Constrained Principal Co-ordinates Analysis

FDR:

False discovery rate

FMT:

Fecal microbiota transplantation

IBD:

Inflammatory bowel disease

MAFLD:

Metabolic associated fatty liver disease

PERMANOVA:

Permutational multivariate analysis of variance

RDP:

Ribosomal Database Project

PICRUSt2:

Phylogenetic Investigation of Communities annotated the pathway information by Reconstruction of Unobserved States

KEGG:

Kyoto Encyclopedia of Genes and Genomes

References

  1. Ungaro R, Mehandru S, Allen PB, Peyrin-Biroulet L, Colombel JF. Ulcerative colitis. Lancet. 2017;389(10080):1756–70. https://doi.org/10.1016/S0140-6736(16)32126-2.

    Article  PubMed  Google Scholar 

  2. Round JL, Mazmanian SK. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol. 2009;9(5):313–23. https://doi.org/10.1038/nri2515.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Scott A. Gut-liver axis: menace in the microbiota. Nature. 2017;551(7681):S94–5. https://doi.org/10.1038/d41586-017-06924-3.

    Article  CAS  PubMed  Google Scholar 

  4. Kotsiliti E. Gut microbiome and autism spectrum disorder. Nat Rev Gastroenterol Hepatol. 2022;19(1):6. https://doi.org/10.1038/s41575-021-00564-9.

    Article  CAS  PubMed  Google Scholar 

  5. Dang AT, Marsland BJ. Microbes, metabolites, and the gut-lung axis. Mucosal Immunol. 2019;12(4):843–50. https://doi.org/10.1038/s41385-019-0160-6.

    Article  CAS  PubMed  Google Scholar 

  6. Doonan J, Tarafdar A, Pineda MA, Lumb FE, Crowe J, Khan AM, Hoskisson PA, Harnett MM, Harnett W. The parasitic worm product ES-62 normalises the gut microbiota bone marrow axis in inflammatory arthritis. Nat Commun. 2019;10(1):1554. https://doi.org/10.1038/s41467-019-09361-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Fang Z, Pan T, Li L, Wang H, Zhu J, Zhang H, Zhao J, Chen W, Lu W. Bifidobacterium longum mediated tryptophan metabolism to improve atopic dermatitis via the gut-skin axis. Gut Microbes. 2022;14(1):2044723. https://doi.org/10.1080/19490976.2022.2044723.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Tan CXW, Brand HS, de Boer NKH, Forouzanfar T. Gastrointestinal diseases and their oro-dental manifestations: part 2: ulcerative colitis. Brit Dent J. 2017;222(1):53–7. https://doi.org/10.1038/sj.bdj.2017.37.

    Article  CAS  PubMed  Google Scholar 

  9. Read E, Curtis MA, Neves JF. The role of oral bacteria in inflammatory bowel disease. Nat Rev Gastroenterol Hepatol. 2021;18(10):731–42. https://doi.org/10.1038/s41575-021-00488-4.

    Article  PubMed  Google Scholar 

  10. Said HS, Suda W, Nakagome S, Chinen H, Oshima K, Kim S, Kimura R, Iraha A, Ishida H, Fujita J, et al. Dysbiosis of salivary microbiota in inflammatory bowel disease and its association with oral immunological biomarkers. DNA Res. 2014;21(1):15–25. https://doi.org/10.1093/dnares/dst037.

    Article  CAS  PubMed  Google Scholar 

  11. Xun Z, Zhang Q, Xu T, Chen N, Chen F. Dysbiosis and ecotypes of the salivary microbiome associated with inflammatory bowel diseases and the assistance in diagnosis of diseases using oral bacterial profiles. Front Microbiol. 2018;9:1136. https://doi.org/10.3389/fmicb.2018.01136.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Bafeta A, Yavchitz A, Riveros C, Batista R, Ravaud P. Methods and reporting studies assessing fecal microbiota transplantation: a systematic review. Ann Intern Med. 2017;167(1):34–9. https://doi.org/10.7326/M16-2810.

    Article  PubMed  Google Scholar 

  13. Xu J, Chen N, Wu Z, Song Y, Zhang Y, Wu N, Zhang F, Ren X, Liu Y. 5-Aminosalicylic acid alters the gut bacterial microbiota in patients with ulcerative colitis. Front Microbiol. 2018;9:1274. https://doi.org/10.3389/fmicb.2018.01274.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Jun X, Ning C, Yang S, Zhe W, Na W, Yifan Z, Xinhua R, Yulan L. Alteration of fungal microbiota after 5-ASA treatment in UC patients. Inflamm Bowel Dis. 2020;26(3):380–90. https://doi.org/10.1093/ibd/izz207.

    Article  PubMed  Google Scholar 

  15. Tedjo DI, Smolinska A, Savelkoul PH, Masclee AA, van Schooten FJ, Pierik MJ, Penders J, Jonkers DM. The fecal microbiota as a biomarker for disease activity in Crohn’s disease. Sci Rep. 2016;6:35216. https://doi.org/10.1038/srep35216.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Mei L, Zhou J, Su Y, Mao K, Wu J, Zhu C, He L, Cui Y. Gut microbiota composition and functional prediction in diarrhea-predominant irritable bowel syndrome. BMC Gastroenterol. 2021;21(1):105. https://doi.org/10.1186/s12876-021-01693-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Peled JU, Gomes ALC, Devlin SM, Littmann ER, Taur Y, Sung AD, Weber D, Hashimoto D, Slingerland AE, Slingerland JB, et al. Microbiota as predictor of mortality in allogeneic hematopoietic-cell transplantation. N Engl J Med. 2020;382(9):822–34. https://doi.org/10.1056/NEJMoa1900623.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Ingham AC, Kielsen K, Mordhorst H, Ifversen M, Aarestrup FM, Muller KG, Pamp SJ. Microbiota long-term dynamics and prediction of acute graft-versus-host disease in pediatric allogeneic stem cell transplantation. Microbiome. 2021;9(1):148. https://doi.org/10.1186/s40168-021-01100-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Bernstein CN, Eliakim A, Fedail S, Fried M, Gearry R, Goh KL, Hamid S, Khan AG, Khalif I, Ng SC, et al. World gastroenterology organisation global guidelines inflammatory bowel disease: update august 2015. J Clin Gastroenterol. 2016;50(10):803–18. https://doi.org/10.1097/MCG.0000000000000660.

    Article  PubMed  Google Scholar 

  20. Xu J, Zhang JN, Sun BH, Liu Q, Ma J, Zhang Q, Liu YX, Chen N, Chen F. The role of genotype and diet in shaping gut microbiome in a genetic vitamin A deficient mouse model. J Genet Genomics. 2022;49(2):155–64. https://doi.org/10.1016/j.jgg.2021.08.015.

    Article  CAS  PubMed  Google Scholar 

  21. Rognes T, Flouri T, Nichols B, Quince C, Mahe F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4: e2584. https://doi.org/10.7717/peerj.2584.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8. https://doi.org/10.1038/nmeth.2604.

    Article  CAS  PubMed  Google Scholar 

  23. Ren X, Xu J, Zhang Y, Chen G, Zhang Y, Huang Q, Liu Y. Bacterial alterations in post-cholecystectomy patients are associated with colorectal cancer. Front Oncol. 2020;10:1418. https://doi.org/10.3389/fonc.2020.01418.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Zhang Y, Xu J, Wang X, Ren X, Liu Y. Changes of intestinal bacterial microbiota in coronary heart disease complicated with nonalcoholic fatty liver disease. BMC Genomics. 2019;20(1):862. https://doi.org/10.1186/s12864-019-6251-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, Brown CT, Porras-Alfaro A, Kuske CR, Tiedje JM. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014;42:D633-642. https://doi.org/10.1093/nar/gkt1244.

    Article  CAS  PubMed  Google Scholar 

  26. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21. https://doi.org/10.1038/nbt.2676.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51. https://doi.org/10.1093/nar/gkaa970.

    Article  CAS  PubMed  Google Scholar 

  28. Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4. https://doi.org/10.1093/bioinformatics/btu494.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Villanueva RAM, Chen ZJ. ggplot2: elegant graphics for data analysis, 2nd edition. Meas-Interdiscip Res. 2019;17(3):160–7. https://doi.org/10.1080/15366367.2019.1565254.

    Article  Google Scholar 

  30. Shoeb M, Ramana KV. Anti-inflammatory effects of benfotiamine are mediated through the regulation of the arachidonic acid pathway in macrophages. Free Radical Bio Med. 2012;52(1):182–90. https://doi.org/10.1016/j.freeradbiomed.2011.10.444.

    Article  CAS  Google Scholar 

  31. Hirano A, Umeno J, Okamoto Y, Shibata H, Ogura Y, Moriyama T, Torisu T, Fujioka S, Fuyuno Y, Kawarabayasi Y, et al. Comparison of the microbial community structure between inflamed and non-inflamed sites in patients with ulcerative colitis. J Gastroen Hepatol. 2018;33(9):1590–7. https://doi.org/10.1111/jgh.14129.

    Article  CAS  Google Scholar 

  32. Sasaki M, Shimoyama Y, Kodama Y, Ishikawa T. Abiotrophia defectiva DnaK promotes fibronectin-mediated adherence to HUVECs and induces a proinflammatory response. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22168528.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Qi Y, Wu HM, Yang Z, Zhou YF, Jin L, Yang MF, Wang FY. New insights into the role of oral microbiota dysbiosis in the pathogenesis of inflammatory bowel disease. Dig Dis Sci. 2022;67(1):42–55. https://doi.org/10.1007/s10620-021-06837-2.

    Article  PubMed  Google Scholar 

  34. Brown EM, Ke XB, Hitchcock D, Jeanfavre S, Avila-Pacheco J, Nakata T, Arthur TD, Fornelos N, Heim C, Franzosa EA, et al. Bacteroides-derived sphingolipids are critical for maintaining intestinal homeostasis and symbiosis. Cell Host Microbe. 2019;25(5):668-680.e667. https://doi.org/10.1016/j.chom.2019.04.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Shi Y, Lu Y, You J. Unfolded protein response in the activation-induced biological processes of CD8(+) T cells. Pharmacol Res. 2021;169: 105654. https://doi.org/10.1016/j.phrs.2021.105654.

    Article  CAS  PubMed  Google Scholar 

  36. Yokoyama N, Hanafusa K, Hotta T, Oshima E, Iwabuchi K, Nakayama H. Multiplicity of glycosphingolipid-enriched microdomain-driven immune signaling. Int J Mol Sci. 2021;22(17):9565. https://doi.org/10.3390/ijms22179565.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Kushnareva Y, Mathews IT, Andreyev AY, Altay G, Lindestam Arlehamn CS, Pandurangan V, Nilsson R, Jain M, Sette A, Peters B, et al. Functional analysis of immune signature genes in Th1* memory cells links ISOC1 and pyrimidine metabolism to IFN-gamma and IL-17 production. J Immunol. 2021;206(6):1181–93. https://doi.org/10.4049/jimmunol.2000672.

    Article  CAS  PubMed  Google Scholar 

  38. Graczyk D, White RJ, Ryan KM. Involvement of RNA polymerase III in immune responses. Mol Cell Biol. 2015;35(10):1848–59. https://doi.org/10.1128/mcb.00990-14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Campbell C, McKenney PT, Konstantinovsky D, Isaeva OI, Schizas M, Verter J, Mai C, Jin WB, Guo CJ, Violante S, et al. Bacterial metabolism of bile acids promotes generation of peripheral regulatory T cells. Nature. 2020;581(7809):475–9. https://doi.org/10.1038/s41586-020-2193-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Subramanian S, Geng H, Tan XD. Cell death of intestinal epithelial cells in intestinal diseases. Sheng Li Xue Bao. 2020;72(3):308–24.

    PubMed  PubMed Central  Google Scholar 

  41. Koop L, Garg R, Nguyen T, Gujjula NR, Velagapudi M. Bacillus cereus: beyond gastroenteritis. WMJ. 2021;120(2):145–7.

    PubMed  Google Scholar 

  42. Yang Y, Jobin C. Novel insights into microbiome in colitis and colorectal cancer. Curr Opin Gastroenterol. 2017;33(6):422–7. https://doi.org/10.1097/mog.0000000000000399.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Satou N, Ishihara K, Hiratsuka M, Tanaka H, Endo Y, Saito S, Iwakura Y, Leonard WJ, Hirasawa N. Induction of thymic stromal lymphopoietin production by xylene and exacerbation of picryl chloride-induced allergic inflammation in mice. Int Arch Allergy Immunol. 2012;157(2):194–201. https://doi.org/10.1159/000327545.

    Article  CAS  PubMed  Google Scholar 

  44. Vasapolli R, Schutte K, Schulz C, Vital M, Schomburg D, Pieper DH, Vilchez-Vargas R, Malfertheiner P. Analysis of transcriptionally active bacteria throughout the gastrointestinal tract of healthy individuals. Gastroenterology. 2019;157(4):1081-1092.e1083. https://doi.org/10.1053/j.gastro.2019.05.068.

    Article  CAS  PubMed  Google Scholar 

  45. Lloyd-Price J, Mahurkar A, Rahnavard G, Crabtree J, Orvis J, Hall AB, Brady A, Creasy HH, McCracken C, Giglio MG, et al. Strains, functions and dynamics in the expanded Human Microbiome Project. Nature. 2017;551(7679):61–6. https://doi.org/10.1038/nature24485.

    Article  CAS  Google Scholar 

  46. Hunt RH, Camilleri M, Crowe SE, El-Omar EM, Fox JG, Kuipers EJ, Malfertheiner P, McColl KEL, Pritchard DM, Rugge M, et al. The stomach in health and disease. Gut. 2015;64(10):1650–68. https://doi.org/10.1136/gutjnl-2014-307595.

    Article  CAS  PubMed  Google Scholar 

  47. Jia W, Xie G, Jia W. Bile acid-microbiota crosstalk in gastrointestinal inflammation and carcinogenesis. Nat Rev Gastroenterol Hepatol. 2018;15(2):111–28. https://doi.org/10.1038/nrgastro.2017.119.

    Article  CAS  PubMed  Google Scholar 

  48. Ducarmon QR, Zwittink RD, Hornung BVH, van Schaik W, Young VB, Kuijper EJ. Gut microbiota and colonization resistance against bacterial enteric infection. Microbiol Mol Biol Rev. 2019;83(3):e00007-00019. https://doi.org/10.1128/MMBR.00007-19.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Soderholm AT, Pedicord VA. Intestinal epithelial cells: at the interface of the microbiota and mucosal immunity. Immunology. 2019;158(4):267–80. https://doi.org/10.1111/imm.13117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Donaldson GP, Lee SM, Mazmanian SK. Gut biogeography of the bacterial microbiota. Nat Rev Microbiol. 2016;14(1):20–32. https://doi.org/10.1038/nrmicro3552.

    Article  CAS  PubMed  Google Scholar 

  51. Atarashi K, Suda W, Luo C, Kawaguchi T, Motoo I, Narushima S, Kiguchi Y, Yasuma K, Watanabe E, Tanoue T, et al. Ectopic colonization of oral bacteria in the intestine drives TH1 cell induction and inflammation. Science. 2017;358(6361):359–65. https://doi.org/10.1126/science.aan4526.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kitamoto S, Nagao-Kitamoto H, Jiao YZ, Gillilland MG, Hayashi A, Imai J, Sugihara K, Miyoshi M, Brazil JC, Kuffa P, et al. The intermucosal connection between the mouth and gut in commensal pathobiont-driven colitis. Cell. 2020;182(2):447-462.e414. https://doi.org/10.1016/j.cell.2020.05.048.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Liu X, Mao B, Gu J, Wu J, Cui S, Wang G, Zhao J, Zhang H, Chen W. Blautia-a new functional genus with potential probiotic properties? Gut Microbes. 2021;13(1):1–21. https://doi.org/10.1080/19490976.2021.1875796.

    Article  CAS  PubMed  Google Scholar 

  54. Lenoir M, Martin R, Torres-Maravilla E, Chadi S, Gonzalez-Davila P, Sokol H, Langella P, Chain F, Bermudez-Humaran LG. Butyrate mediates anti-inflammatory effects of Faecalibacterium prausnitzii in intestinal epithelial cells through Dact3. Gut Microbes. 2020;12(1):1–16. https://doi.org/10.1080/19490976.2020.1826748.

    Article  CAS  PubMed  Google Scholar 

  55. Manzetti S, Zhang J, van der Spoel D. Thiamin function, metabolism, uptake, and transport. Biochemistry-Us. 2014;53(5):821–35. https://doi.org/10.1021/bi401618y.

    Article  CAS  Google Scholar 

  56. Klaassen MAY, Imhann F, Collij V, Fu J, Wijmenga C, Zhernakova A, Dijkstra G, Festen EAM, Gacesa R, Vich Vila A, et al. Anti-inflammatory Gut microbial pathways are decreased during Crohn’s disease exacerbations. J Crohns & Colitis. 2019;13(11):1439–49. https://doi.org/10.1093/ecco-jcc/jjz077.

    Article  Google Scholar 

  57. De Preter V, Machiels K, Joossens M, Arijs I, Matthys C, Vermeire S, Rutgeerts P, Verbeke K. Faecal metabolite profiling identifies medium-chain fatty acids as discriminating compounds in IBD. Gut. 2015;64(3):447–58. https://doi.org/10.1136/gutjnl-2013-306423.

    Article  CAS  PubMed  Google Scholar 

  58. Vieira JRP, Rezende ATO, Fernandes MR, da Silva NA. Intestinal microbiota and active systemic lupus erythematosus: a systematic review. Adv Rheumatol. 2021;61(1):42. https://doi.org/10.1186/s42358-021-00201-8.

    Article  PubMed  Google Scholar 

  59. Yang M, Sun B, Li J, Yang B, Xu J, Zhou X, Yu J, Zhang X, Zhang Q, Zhou S, et al. Alteration of the intestinal flora may participate in the development of Graves’ disease: a study conducted among the Han population in southwest China. Endocr Connect. 2019;8(7):822–8. https://doi.org/10.1530/EC-19-0001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Strauss J, Kaplan GG, Beck PL, Rioux K, Panaccione R, Devinney R, Lynch T, Allen-Vercoe E. Invasive potential of gut mucosa-derived Fusobacterium nucleatum positively correlates with IBD status of the host. Inflamm Bowel Dis. 2011;17(9):1971–8. https://doi.org/10.1002/ibd.21606.

    Article  PubMed  Google Scholar 

  61. Engevik MA, Danhof HA, Ruan W, Engevik AC, Chang-Graham AL, Engevik KA, Shi Z, Zhao Y, Brand CK, Krystofiak ES, et al. Fusobacterium nucleatum secretes outer membrane vesicles and promotes intestinal inflammation. MBio. 2021;12(2):e02706-02720. https://doi.org/10.1128/mBio.02706-20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Callahan SM, Hancock TJ, Johnson JG. Characterization of Campylobacter jejuni-neutrophil interactions. Curr Protoc. 2021;1(11): e294. https://doi.org/10.1002/cpz1.294.

    Article  CAS  PubMed  Google Scholar 

  63. Jing X, Korchagina AA, Shein SA, Muraoka WT, Koroleva E, Tumanov AV. IL-23 contributes to Campylobacter jejuni-induced intestinal pathology via promoting IL-17 and IFNgamma responses by innate lymphoid cells. Front Immunol. 2020;11: 579615. https://doi.org/10.3389/fimmu.2020.579615.

    Article  CAS  PubMed  Google Scholar 

  64. Callahan S, Doster RS, Jackson JW, Kelley BR, Gaddy JA, Johnson JG. Induction of neutrophil extracellular traps by Campylobacter jejuni. Cell Microbiol. 2020;22(8): e13210. https://doi.org/10.1111/cmi.13210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Wang J, Jia Z, Zhang B, Peng L, Zhao F. Tracing the accumulation of in vivo human oral microbiota elucidates microbial community dynamics at the gateway to the GI tract. Gut. 2020;69(7):1355–6. https://doi.org/10.1136/gutjnl-2019-318977.

    Article  PubMed  Google Scholar 

  66. Duran-Pinedo A, Solbiati J, Teles F, Teles R, Zang Y, Frias-Lopez J. Long-term dynamics of the human oral microbiome during clinical disease progression. BMC Biol. 2021;19(1):240. https://doi.org/10.1186/s12915-021-01169-z.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This work was supported by the Beijing Science and technology planning project (No. Z171100000417022 to LYL), the Beijing Municipal Natural Science Foundation (No. 7232196 to LYL, No. 7214267 to XJ), the National Natural Science Foundation of China (No. 82070539, 81873549 to LYL, and 82000496 to XJ), and Peking University People’s Hospital Scientific Research Development Funds (No. PTU2021-03 to LYL, No. RDY2020-21, RS2021-09, and RDL2021-11 to XJ).

Author information

Authors and Affiliations

Authors

Contributions

LYL and XJ designed this study. XJ, HYB, WY and ZYF performed subject enrollment and sample collection. XJ, YL, ZY, FXH, KZL, ZY, ZJH, SHT and CN documented and analyzed clinical parameters. XJ performed microbial bioinformatics analysis and visualization. XJ, ZY, and FXH wrote, and LYL and XJ revised this manuscript.

Corresponding author

Correspondence to Yu-Lan Liu.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Medical Ethics Review Board of Peking University People’s Hospital.

Consent for publication

Obtained. The authors thank all of the subjects who participated in this study.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

The comparative analysis of the salivary bacterial community. A, The constrained principal coordinate analysis (CPCoA) for three GI site samples based on the Bray-Curtis distance. The R software (v 4.0.1) with the vegan (v 2.5-7) package were used and P-values were obtained using permutational multivariate analysis of variance (PERMANOVA) (same methodology as B). B, The CPCoA for salivary samples based on the Bray-Curtis distance. C to E, The comparative analysis of salivary bacterial genus abundance between two groups ((C) OU vs. CON; (D) UC vs. CON; (E) UC_OU vs. CON). The EdgeR package was used for comparative analysis. The difference between the two groups is shown as a Manhattan diagram. Point shape indicates the genus enriched, depleted, or not significant in the former group compared with the latter. Point size indicates the counts of a specific genus. CPM, count per million. F to H, The comparative analysis of salivary bacterial function between two groups ((F) OU vs. CON; (G) UC vs. CON; (H) UC_OU vs. CON). Phylogenetic Investigation of Communities annotated the pathway information by Reconstruction of Unobserved States (PICRUSt2) software by referring to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The STAMP software was used for data visualization. Abbreviations: CON, healthy controls; OU, patients with only oral ulcers; UC, UC patients without oral ulcers; UC_OU, UC patients with oral ulcers.

Additional file 2: Figure S2.

The comparative analysis of the buccal bacterial community. A, The constrained principal coordinate analysis (CPCoA) for buccal samples based on the Bray-Curtis distance. The R software (v 4.0.1) with the vegan (v 2.5-7) package were used and P-values were obtained using permutational multivariate analysis of variance (PERMANOVA). B to D, The comparative analysis of buccal bacterial genus abundance between two groups ((B) OU vs. CON; (C) UC vs. CON; (D) UC_OU vs. CON). The EdgeR package was used for comparative analysis. The difference between the two groups is shown as a Manhattan diagram. Point shape indicates the genus enriched, depleted, or not significant in the former group compared with the latter. Point size indicates the counts of a specific genus. CPM, count per million. E to G, The comparative analysis of buccal bacterial function between two groups ((E) OU vs. CON; (F) UC vs. CON; (G) UC_OU vs. CON). Phylogenetic Investigation of Communities annotated the pathway information by Reconstruction of Unobserved States (PICRUSt2) software by referring to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The STAMP software was used for data visualization. Abbreviations: CON, healthy controls; OU, patients with only oral ulcers; UC, UC patients without oral ulcers; UC_OU, UC patients with oral ulcers.

Additional file 3: Figure S3.

The comparative analysis of the fecal bacterial community. A, The constrained principal coordinate analysis (CPCoA) for fecal samples based on the Bray-Curtis distance. The R software (v 4.0.1) with the vegan (v 2.5-7) package were used and P-values were obtained using permutational multivariate analysis of variance (PERMANOVA). B to D, The comparative analysis of fecal bacterial genus abundance between two groups ((B) OU vs. CON; (C) UC vs. CON; (D) UC_OU vs. CON). The EdgeR package was used for comparative analysis. The difference between the two groups is shown as a Manhattan diagram. Point shape indicates the genus enriched, depleted, or not significant in the former group compared with the latter. Point size indicates the counts of a specific genus. CPM, count per million. E to G, The comparative analysis of bacterial function between two groups ((E) OU vs. CON; (F) UC vs. CON; (G) UC_OU vs. CON). Phylogenetic Investigation of Communities annotated the pathway information by Reconstruction of Unobserved States (PICRUSt2) software by referring to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The STAMP software was used for data visualization. Abbreviations: CON, healthy controls; OU, patients with only oral ulcers; UC, UC patients without oral ulcers; UC_OU, UC patients with oral ulcers.

Additional file 4: Figure S4.

The alteration of clinical parameters from baseline (Before) to the end of follow-up in UC patients with or without oral ulcers. A, The mayo clinic score. B, Erythrocyte sedimentation rate (ESR); C, C-reaction protein (CRP); D, Immunoglobulin A (IgA); E, Immunoglobulin G (IgG); F, Immunoglobulin G (IgG); G, Complement 3 (C3); H, Complement 4 (C4).

Additional file 5: Supplementary tables.

Table S1. Bacterial alpha diversity of each sample. Table S2. Bray Curtis distance of each sample. Table S3. Shared bacterial families and genera in three different GI niches. Table S4. Specific bacterial abundance for the ternary plot. Table S5. Relative bacterial family abundance of each sample. Table S6. Relative bacterial genus abundance of each sample. Table S7. Comparative analysis for bacterial community in the salivary samples (all results presented). Table S8. Comparative analysis for bio-function in the salivary samples (only significant results presented). Table S9. Comparative analysis for bacterial community in the buccal samples (all results presented). Table S10. Comparative analysis for bio-function in the buccal samples (only significant results presented). Table S11. Comparative analysis for bacterial community in the fecal samples (all results presented). Table S12. Comparative analysis for bio-function in the fecal samples (only significant results presented). Table S13. Correlative analysis between clinical parameters (inflammatory and immunological) and bacterial diversity (alpha and beta). Table S14. Bacterial comparative analysis of UC patients with or without treatment response at three different GI niches.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, J., Zhang, Y., Fang, XH. et al. The oral bacterial microbiota facilitates the stratification for ulcerative colitis patients with oral ulcers. Ann Clin Microbiol Antimicrob 22, 99 (2023). https://doi.org/10.1186/s12941-023-00646-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12941-023-00646-3

Keywords