Alterations in the human oral microbiota in systemic lupus erythematosus
Journal of Translational Medicine volume 21, Article number: 95 (2023)
Alterations in oral microbiota in patients with systemic lupus erythematosus (SLE) is less evaluated. The aim of this study was to compare the characteristics of the oral microbiome in SLE patients and healthy controls, and construct an SLE classifier based on the oral microbiota.
We sequenced tongue-coating samples of individuals in treatment-naïve SLE (n = 182) and matched healthy controls (n = 280). We characterized the oral microbiome and constructed a microbial classifier in the derivation cohort and validated the results in the validation cohorts. Furthermore, the oral microbiome of posttreatment SLE (n = 73) was characterized.
The oral microbial diversity of SLE was increased, and the microbial community was different between SLE and healthy controls. The genera Prevotella and Veillonella were enriched, while Streptococcus and Porphyromonas were reduced in SLE. In addition, an increase was noted in 27 predicted microbial functions, while a decrease was noted in 34 other functions. Thirty-nine operational taxonomy units (OTUs) were identified to be related with seven clinical indicators. Two OTUs were identified to construct a classifier, which yielded area under the curve values of 0.9166 (95% CI 0.8848–0.9483, p < 0.0001), 0.8422 (95% CI 0.7687–0.9157, p < 0.0001), and 0.8406 (95% CI 0.7677–0.9135, p < 0.0001) in the derivation, validation, and cross-regional validation groups, respectively. Moreover, as disease activity increased, Abiotrophia and Lactobacillales increased, while Phyllobacterium and unclassified Micrococcusaceae decreased. Finally, nine OTUs were selected to construct a classifier distinguishing posttreatment SLE patients from healthy controls, which achieved a diagnostic efficacy of 0.9942 (95% CI 0.9884–1, p < 0.0001).
Our study comprehensively characterizes the oral microbiome of SLE and shows the potential of the oral microbiota as a non-invasive diagnostic biomarker in SLE.
Systemic lupus erythematosus (SLE) is a chronic autoimmune disease featured by loss of autoimmune tolerance and activation of autoimmune response, ultimately resulting in the formation of multiple autoantibodies and damage to various organ systems . It affects females of childbearing age ten times more than males . The clinical manifestations of SLE are markedly heterogeneous, ranging from mild symptoms to life-threatening conditions . To date, the etiology and pathogenesis of SLE is unclear, and it is believed that genetic, race, immune disorders, hormone, and environmental factors may together contribute to the occurrence of the disease .
All the microorganisms that inhabit the human body, including bacteria, fungi, archaea, viruses and protozoa, form the human commensal microbiota . Despite the primary habitat being the gut, thriving microbial populations can be found in areas throughout the body, such as the skin, digestive system, respiratory system, and reproductive system . As the largest microsystem in the human body [6, 7], the gut microbiota is considered to be an important environmental trigger for initiating and promoting the progression of SLE through molecular mimicry, interference with the host metabolism axis, and engaging in type I interferon pathways [8,9,10,11,12,13,14]. In SLE, the gut microbiota was disturbed with decreased bacterial diversity and altered flora composition and function [9, 15,16,17], and the gut microbiota was altered along with disease activity changes [18, 19].
The oral microbiota is also an important microecosystem in the human body that is composed of more than 700 unique bacterial species .Oral microbiota dysbiosis has been found in many diseases [21,22,23,24], and it has been shown that oral microbiota biomarkers could apply to the diagnosis for rheumatoid arthritis (RA) and COVID-19 [21, 24]. However, only a small number of researches have focused on the characteristics of oral microbiota in SLE, and the sample sizes of existing studies were small [25,26,27]. We conducted the present study to characterize the oral microbiome profiles in patients with SLE from China. We employed 16S rRNA MiSeq sequencing on tongue coating samples from 255 SLE cases and 280 healthy controls (HC) to profile the unique oral microbiota of SLE, identify the unique oral microbiota biomarkers, and establish an oral microbiota diagnostic model of SLE.
Materials and methods
Participants and ethic statement
The design of this work was based on prospective specimen collection followed by a blinded retrospective assessment. It was conducted according to Helsinki Declaration. Patients with SLE were consecutively selected from the outpatient clinic and inpatient ward of the Rheumatology Department of the First Affiliated Hospital of Zhengzhou University as well as the Rheumatology Department of Hainan General Hospital between July 2019 and July 2021. Additionally, we included healthy volunteers who went to the First Affiliated Hospital of Zhengzhou University for annual physical examination during the same period.
Included patients satisfied the following criteria: met the classification criteria of SLE in the American College of Rheumatology ; aged 18 years or old; did not receive glucocorticoids or any immunosuppressive drugs (treatment-naïve patients), or in stable disease condition under treatment with glucocorticoids (≤ 10 mg/d prednisolone or equivalent) and/or hydroxychloroquine (posttreatment stable patients). The exclusion criteria for SLE were either of the following: complicated with other autoimmune diseases, diabetes, neuropsychiatric diseases, infectious diseases, oral mucosal diseases, gingival and throat diseases; taken antibiotics, probiotics, yogurt, black tea, pickles or other fermented foods in the past 4 weeks; or lacked of clinical data. The collected data of participants included age, gender, body mass index (BMI), routine blood and urine tests, HBsAg, HBsAb, anti-HCV, TPPA, HIV antibody, liver function, kidney function. Data of antinuclear antibody, anti- extractable nuclear antigen antibody, complement C3 and C4, immunoglobulins A, G and M were collected from SLE patients. The initial protocol of the study obtained approval of the Ethics Committee at the First Affiliated Hospital of Zhengzhou University (2019-KY-200). All participants signed a written informed consent prior to participation.
All subjects gargled with sterile mouthwash twice before collecting tongue coating samples, and then, professional operators scraped the posterior to anterior middle regions of the tongue coating with swabs. The collected tongue coating specimen was immediately put into the tube and transferred to a − 80 °C refrigerator within 2 h.
DNA extraction, PCR amplification, and sequencing
Extraction of bacterial DNA was done via the E.Z.N.A. ® Stool DNA Kit (Omega Bio-tek, Inc., GA) as per the steps stated in the kit’s instructions. Primers targeting the V3-V4 region (5ʹ- CCTACGGGNGGCWGCAG -3ʹ and 5ʹ-GACTACHVGGGTATCTAATCC-3ʹ) with high variation of 16S rRNA were used for PCR amplification of the extracted DNA samples. DNA sequences were identified on an Illumina MiSeq platform (Illumina Inc., USA) by Shanghai Mobio Biomedical Technology Co., Ltd. Raw data of the collected samples were stored in the Sequence Read Archive database (PRJNA789129).
USEARCH (version 11.0.667) was used to process raw sequencing data. UPARSE (http://drive5.com/uparse/) was used to classify operational taxonomy units (OTUs) based on 97% sequence similarity, and OTUs were annotated using the SILVA reference database (SSU138). The Shannon, observed OTUs, Ace, Chao and Simpson indices were calculated by using Mothur v1.42.1 to evaluate α diversity. The R package (version 3.6.0) was used to analyze β diversity of bacteria, including nonmetric multidimensional scaling (NMDS), principal coordinate analysis (PCoA) and principal component analysis (PCA).
As per the normalized relative abundance matrix, linear discriminant analysis (LDA) and LDA effect size (LEfSe) were adopted (https://github.com/SegataLab/lefse) to analyze the characteristics between different groups of bacteria. The threshold was set to LDA score (log10) = 3. PICRUSt2 v2.4.1 (https://github.com/picrust/picrust2/wiki) was employed to predict functional abundances based on 16S rRNA gene sequences.
By verifying the reads in the representative sequence of the derivation and the validation sets, the OTU frequencies of the derivation and the validation sets are obtained. Using the Wilcoxon test (p < 0.05), significant OTU biomarkers were isolated to undergo further analyses. The identified OTUs from the derivation set was used for construct diagnostic model based on random forest model. To obtain the optimal combination of diagnostic markers, we ranked OTUs from high to low according to their importance values (Mean Decrease Accuracy), increased the number of OTUs in order, and calculated the error of constructing the diagnostic model using different numbers of biomarkers. Five-fold cross-validation was also performed in the derivation set to obtain five different cross-validation (CV) error curves. The mean value of CV error is calculated based on the data of the five curves, and a line graph drawn with the mean value is obtained. We select the minimum CV-error in the mean curve plus the standard deviation of the CV error as the cut-off value. We list all sets of OTU markers whose CV error is less than the cut-off value and select the set with the smallest number of OTUs as the optimal set. Finally, the probability of disease (POD) index of the derivation and validation sets were calculated by the determined optimal OTU set. The POD index refers to the ratio of the number of randomly generated decision trees whose predicted sample is SLE to to the number of predicted HC [24, 29]. Then, a receiver operating characteristic (ROC) curve (R 3.3.0, pROC package) was constructed to evaluated the model through the area under the curve (AUC).
Continuous variables that are normally and non-normal distributed were represented by the mean (standard deviation) and median (inter-quartile range), respectively. Categorical variables were expressed as percentages. To compare the difference between SLE patients and HC, Student's t-test and Wilcoxon rank-sum test were used for continuous variables with normal and non-normal distribution, respectively. In addition, the Chi-square test or Fisher's exact test was used for categorical variables. SPSS V.21.0 (SPSS, Chicago, Illinois, USA) was used to conduct the statistical analysis. Differences were identified as statistically significant with p < 0.05.
Study design and characteristics of the participants
An overall number of 563 tongue coating samples were collected from variable areas around China prospectively. After strict exclusion and diagnosis procedures, 535 tongue coating samples were analyzed, including 213 SLE (140 treatment-naïve patients and 73 posttreatment patients) from Zhengzhou, 280 HC from Zhengzhou, and 42 treatment-naïve SLE patients from Haikou. The treatment-naïve samples from Zhengzhou were randomly divided into a derivation cohort (100 SLE vs. 200 HC) and an independent validation cohort (40 SLE vs. 80 HC) based on a random number table. Characterization of the oral microbiome, identification of the key microbial markers, and construction of an SLE classifier were performed in the derivation cohort. The efficacy of diagnosis of the SLE classifier was validated in the independent validation and cross-regional validation cohort (42 SLE from Haikou vs. 80 HC). We further divided the 140 SLE patients from Zhengzhou into three groups according to the SLE disease activity index (SLEDAI) and characterized the oral microbiome in these three groups [30, 31]. Finally, we compared the oral microbiome of 73 posttreatment SLE patients and 146 HC (Fig. 1).
In the derivation and validation cohorts, the genders, ages, and BMI of SLE and HC were matched. Serum levels of white blood cells, red blood cells, hemoglobin, platelets, and lymphocyte counts were decreased significantly, and serum levels of globulin were increased significantly in SLE compared with HC. The detail is shown in Table 1.
Increased oral microbial diversity in SLE
During the derivation cohort, the analysis of rarefaction revealed that the count of OTU richness almost reached saturation in both groups, and it was increased significantly in SLE versus HC (Fig. 2A). The results of the rarefaction analysis were shown in Additional file 1: Fig. S1. As measured by the Shannon and Simpson indices, the oral microbial diversity increased significantly in SLE compared with HC (all p < 0.001), (Additional file 2: Data S1, Fig. 2B, C).
In addition, beta diversity was conducted to display microbiome space between each sample. As estimated by NMDS and PCoA analysis, the OTU distribution differed significantly between the two groups (Fig. 2D, E, Additional file 1: Fig. S2). A Venn diagram illustrated the overlap of OTUs between SLE and HC. Totally, 709 OTUs were shared, while 70 OTUs were specific to SLE (Fig. 2F). The key 46 OTUs between the two groups was shown in Additional file 1: Fig. S3; (Additional file 2: Data S2). Notably, 24 OTUs showed significant enrichment in SLE, and 22 OTUs showed significant enrichment in HC.
Phylogenetic profiles of the oral microbiome in SLE
The taxa composition and variation of the oral microbiota were compared in SLE and HC in the derivation cohort. At the phylum and genus levels, the relative abundance of taxa of oral microbiota in each sample were displayed in Additional file 1: Figs. S4 and S5 (Additional file 2: Data S3 and S4), and the average composition and relative abundance of taxa in both groups were shown in Fig. 3A and C, respectively. Bacteroidota, Firmicutes and Proteobacteria were the three most predominant phyla in both groups.
We further analyzed bacterial abundance to identify key differentially abundant bacteria between the two groups. As shown in Fig. 3B, SLE patients exhibited increased abundance of six phyla, including Bacteroidota, Firmicutes, Fusobacteriota, Patescibacteria, Campilobacterota and Cyanobacteria, and decreased abundance of Proteobacteria, Actinobacteriota, and Bacteria_unclassified compared with HC (p < 0.05), (Additional file 2: Data S5). Figure 3D showed that increased abundance in 28 bacteria, including Prevotella, Veillonella, Leptotrichia, Alloprevotella, and Lachnoanaerobaculum, and decreased abundance in 22 bacteria, including Streptococcus, Porphyromonas, Rothia, Haemophilus, and Granulicatella, were observed in SLE at the genus level (p < 0.05) (Additional file 2: Data S6). In addition, we analyzed the oral microbiota community in the two groups at the class, order, and family levels. The composition and abundance of taxa in each sample at the class level were shown in Additional file 1: Fig. S6, while at the order level and at the family level are shown in Additional file 1: Fig. S9 and Figure S12, respectively. At class level, the relative abundance of Bacteroidia and Gammaproteobacteria is relatively high, reaching more than 50% in both HC and SLE (Additional file 1: Fig. S7); at the order level, the relative abundance in the two groups was higher for Bacteroidales, Burkholderiales and Fusobacteriales (Additional file 1: Fig. S10); while at the family level, the relative abundance in the two groups was higher for Prevotellaceae, Neisseriaceae and Veillonellaceae (Additional file 1: Fig. S13). At the class level, compared to HC, SLE had 8 bacterial populations significantly enriched (e.g., Cyanobacteriia, Coriobacteriia, and Campylobacteriia) and 4 bacterial populations significantly reduced (e.g., Actinobacteria, Bacilli and Gammaproteobacteria) (Additional file 1: Fig. S8 and Additional file 2: Data S7); At the order level, SLE had 13 bacterial populations significantly enriched (e.g., Campylobacterales, Flavobacteriales, and Pseudomonadales) and 12 bacterial populations significantly reduced (e.g., Micrococcales, Pasteurellales and Staphylococcales) (Additional file 1: Fig. S11 and Additional file 2: Data S8); while at the family level, SLE had 17 bacterial populations significantly enriched (e.g., Campylobacteraceae, Sphingomonadaceae, and Chloroplast) and 18 bacterial populations significantly reduced (e.g., Micrococcaceae, Pasteurellaceae and Porphyromonadaceae) (Additional file 1: Figs. S14 and Additional file 2: Data S9).
Crucial bacteria and microbial functions related to SLE
To identify the special bacterial and key bacterial taxa contributing to the alteration of SLE oral microbiota, we performed LDA and LEfSe to identify the difference between SLE and HC. A cladogram revealing the phylogenetic profile of different bacterial taxa showed the maximum differences in taxa between the two groups, which implied dysbiosis of oral microbiota in SLE (Additional file 1: Fig. S15, Additional file 2: Data S10). At the genus level, considerable differences were observed between SLE and HC according to LDA selection. As shown in Fig. 4A, increased abundance in 20 genera, including Prevotella, Leptotrichia, Veillonella, Alloprevotella, and Megasphaera, and decreased abundance in 15 bacteria, including Bacteria_unclassified, Micrococcaceae_unclassified, Gemella, Bacilli_unclassified, and Pasteurellaceae_unclassified, were observed in SLE compared with HC (p < 0.01) (Additional file 2: Data S11).
The functional profiles and the key microbial functions between two groups were shown in Additional file 1: Fig. S16, which implied dysbiosis of functional patterns in SLE. Compared with HC, 27 functional modules, including ansamycin biosynthesis, arginine and ornithine metabolism, and histidine metabolism increased significantly, while 34 functions, including aminobenzoate degradation, unsaturated fatty acid biosynthesis, and taurine and hypotaurine metabolisms decreased significantly in SLE (all p < 0.05), (Fig. 4B, Additional file 2: Data S12). We further analyzed the correlations between oral microbiome and clinical data of SLE by Spearman’s correlation analysis. As shown in Fig. 4C and Additional file 1: Fig. S17. 39 OTUs were related to seven clinical indicators, including gender, age, white blood cell count, hemoglobin, platelet count, lymphocyte count, and globulin (Additional file 2: Data S13).
Diagnostic potential of the oral microbial markers for SLE
For assessment of the diagnostic potential of oral microbial markers in SLE, a random forest classifier model was established in order to distinguish SLE from HC in the derivation group. A total of two OTU markers, including OTU49 (Sphingomonas) and OTU71 (Actinomyces) were identified as the optimal marker set (Fig. 5A, B). The relative abundance of the two OTU markers in each sample was displayed in Additional file 2: Data S14. Then, the two OTUs were used in the calculation of the POD value of each sample in the derivation group (Additional file 2: Data S15). As shown in Fig. 5C and D, the POD value was higher in SLE than in HC (p < 0.05), and it reached an AUC of 0.9166 (95% CI 0.8848–0.9483, p < 0.0001), which implied that the oral microbial marker-based classifier had powerful diagnostic potential for distinguishing SLE from HC.
In addition, an independent validation and a cross-regional validation group were used for validation of the diagnostic efficiency of the microbial marker model for SLE. Additional file 2: Datas S16 and S17 showed the relative abundance of the two OTU markers in each sample in both groups. Figure 5E and G displayed that there was a significant increase in POD value in SLE than in HC in the independent validation and cross-reginal validation groups (all p < 0.05) (Additional file 2: Datas S18 and S19). The AUC value of PODs reached 0.8422 (95% CI 0.7687–0.9157, p < 0.0001) and 0.8406 (95% CI 0.7677–0.9135, p < 0.0001) in the independent validation and cross-regional validation groups (Fig. 5F and H), respectively, which validated the diagnostic potential of the oral microbial marker-based classifier model for SLE.
Oral microbiome alterations in SLE with different disease activities
To specify the characteristics of the oral microbiome in variable disease activities of SLE, the 140 SLE specimens were divided into mild (n = 68, SLEDAI ≤ 6), moderate (n = 58, SLEDAI 7–11), and severe (n = 14, SLEDAI ≥ 12) disease activity groups. Among the three groups, the analysis of rarefaction showed that the number of OTU richness was significantly increased as disease activity increased (Fig. 6A). Moreover, there were 69 OTUs, 78 OTUs, and 18 OTUs unique to the groups with mild, moderate, and severe disease activity, respectively (Fig. 6B). At the genus level, increased abundance in Abiotrophia and Lactobacillales_P5D1-392, and decreased abundance in Phyllobacterium and Micrococcaceae_unclassified were observed as disease activity increased (p < 0.05), (Fig. 6C, Additional file 2: Data S20). Based on LDA analysis, Eubacterium_nodatum_group and Anaerovoracaceae enrichment was significant in the severe disease activity group, while Solobacterium, Erysipelotrichales, and Erysipelotrichaceae enrichment was significant in the mild group (Fig. 6D, Additional file 2: Data S21).
Alterations of oral microbiota in posttreatment stable SLE
To evaluate the oral microbiota features of patients with SLE in posttreatment stable condition, we further collected tongue-coating samples of 73 posttreatment SLE with SLEDAI = 0, and compared the oral microbiota features with those of 146 HC. Oral microbial diversity showed significant increase in posttreatment stable SLE compared to HC (Fig. 7A). The beta diversity analysis showed an obvious difference in the microbiome space community between the two groups (Fig. 7B). At the genus level, the average composition and relative abundance of the oral microbiota in both groups were shown in Fig. 7D. Increased abundance in 21 bacteria and decreased abundance in 21 bacteria were observed in stable SLE versus HC (Fig. 7C, Additional file 2: Data S22). The heatmap showing the top 50 different OTUs revealed that 22 OTUs enrichment was seen in SLE, and 28 OTUs enrichment was seen in HC (Fig. 7E). A random forest classifier model was established between 73 SLE and 146 HC using nine OTUs to identify the diagnostic potential of oral microbial markers for posttreatment stable SLE (Fig. 7F). The POD estimated value in SLE was considerably higher than in HC (p < 0.05), (Fig. 7G, Additional file 2: Data S23), with an AUC of 0.9942 (95% CI 0.9884–1, p < 0.0001) (Fig. 7H), which indicated that the oral microbial marker-based classifier model had powerful diagnostic potential for distinguishing posttreatment stable SLE from HC.
During the study, we comprehensively characterized the oral microbiota of treatment-naïve SLE cases in a large sample of the Chinese population by 16S rRNA MiSeq sequencing. We identified the crucial bacterial and microbial functions related to SLE and investigated the correlativity between oral microbiota and clinically collected data of SLE. Notably, we constructed a microbial classifier that demonstrated strong diagnostic efficacy in distinguishing SLE from HC and achieved cross-regional validation. In addition, we profiled the oral microbiome in SLE with different disease activities. Finally, we analyzed the oral microbiota alterations in posttreatment SLE and identified the optimal microbiota biomarkers. These findings suggested that oral microbiota dysbiosis may be involved in SLE development and recovery, and the microbiota classifier may have potential as a non-invasive biomarker that can be valuable in SLE diagnosis.
We found that the oral microbiome was dysbiosis in SLE, which displayed increased oral microbiome diversity and an altered bacterial community. One study showed that SLE patients had higher oral microbiome diversity than HC in subjects without periodontal disease . Another study observed a difference in oral microbiome diversity between SLE and primary Sjogren’s Syndrome, and SLE had an increased oral microbiome diversity compared with primary Sjogren’s Syndrome . Li et al. found that oral microbiome diversity was decreased in SLE patients compared to HC . Yet, the study sample was small and enrolled both treatment-naïve and posttreatment SLE patients, which may weaken the reliability of the findings. Moreover, diet and environment may partly contribute to the different results between different studies. Additionally, we found that the taxa composition, crucial bacteria, and predicted microbial functions of the oral microbiome were different between SLE and HC. Zhang et al. found that the oral microbiota bacterial community and richness were different between RA and HC . Recently, we characterized the oral microbiome of patients with COVID-19 and found that the phylogenetic profiles, crucial bacteria and predicted microbial functions of COVID-19 differed significantly from those of HC . In addition, one study found that in diabetic mice, the disturbed oral microbiota possesses enhanced pathogenicity . Importantly, studies have shown that disturbance of the gut microbiota may trigger and promote the autoimmunity and inflammatory response in SLE by various pathways [8,9,10,11,12,13,14]. These findings suggested that SLE owns its unique oral microbiome and the altered oral microbiota may be implicated in the initiation and pathogenesis of SLE.
The specific oral microbiome can be used as a diagnostic marker, and this has been studied in many diseases. Zhang et al. developed a diagnostic model in RA using the oral microbiota and verified its powerful diagnostic efficacy in an independent cohort . Flemer established an oral microbiota-based diagnostic model that showed high specificity in identifying colorectal cancer from HC . In this study, we developed and validated an SLE classifier based on two optimal OTUs of the oral microbiome. We calculated the POD of each sample in the derivation, independent validation, and cross-regional validation groups, which yielded AUCs of 0.9926, 0.9809, and 0.9851, respectively. This study is the first to successfully conduct an oral microbiota-based SLE classifier and achieve independent validation and cross-regional validation, suggesting that the oral microbiota could potentially serve as a non-invasive diagnostic method for SLE.
In addition, further analyses were performed to establish the correlation between oral microbiota and clinical indices in SLE. This analysis revealed a close association between the oral microbiota and seven clinical indices. As shown in our previous study, the abundance of specific bacteria in the gut microbiota was related to clinical indicators in chronic renal disease . Ma et al. also found a strong correlation between gut microbiota and serum indices in autoimmune hepatitis . Another study demonstrated a significant correlation between the variation in gut microbiota and clinical data in SLE . Notably, one important finding of our study was that the oral microbiome characteristics differed in different SLE disease activities. As disease activity increased in SLE, the abundance of Abiotrophia and Lactobacillales_P5D1-392 increased, and the abundance of Phyllobacterium and Micrococcaceae_unclassified decreased. These results suggested that oral microbiota alteration may be involved in SLE progression and acceleration. Currently, studies have found that dysbiosis of gut microbiota could play a role in SLE progression, and disease activity can be partly relieved by probiotics, diet intervention, and antibiotics . We hypothesized that interventions targeting the oral microbiota may relieve disease activity and become a therapeutic approach for SLE in the future.
It has been observed that the microbiota altered as disease recovery. One study found that the perturbed oral microbiota partly returned to normal in RA after treatment . Similar results have been observed in periodontitis . Our previous study found that the oral microbiota was related to disease recovery in COVID-19 . During this research, the features of oral microbiota of SLE patients in stable posttreatment conditions were evaluated. The results showed that although patients with SLE achieved remission with neither clinical nor serological activity following treatment, the oral microbiota was dysbiosis and differed from that of HC. Moreover, we established a diagnostic model based on oral microbiota that could differentiate posttreatment stable SLE patients from HC. The results suggested that the oral microbiota may be involved in the disease recovery of SLE.
It is worth doing note that antioxidants may have a potential therapeutic effect in SLE patients. Studies have shown that T cells in lupus mice and SLE patients are continuously activated and therefore have hyperpolarized mitochondria compared to T cells in healthy mice or healthy individuals, which can lead to increased production of reactive oxygen species [35, 36]. In addition, the levels of glutathione and cysteine in serum and peripheral blood individual nuclei of SLE patients are lower than with healthy individuals, suggesting the presence of oxidative stress . It has been shown that antioxidants have specific pharmacological effects in autoimmune diseases [38, 39].
In summary, we firstly profiled the characteristic alterations of the oral (tongue) microbiota in SLE across a large cohort, identified the crucial bacteria and predicted the specific microbial functions related to SLE. Additionally, we revealed that the oral microbiota was correlated with clinical indicators in SLE and that the oral microbiota was altered as disease activity increased. These results suggested that oral microbiota dysbiosis may contribute to the development of SLE. In addition, we are the first to develop a microbial-based SLE classifier and validate its diagnostic efficiency in patients from multiple areas around China, providing a non-invasive diagnostic tool for SLE. Furthermore, we characterized the oral microbiota in posttreatment SLE patients under stable conditions, which suggested that the oral microbiota may has a relation to disease recovery in SLE. However, the outcomes need validation in multicenter cohorts of individuals of different ethnicities. In the end, these findings suggested that the oral microbiota may be a therapeutic target for exploring new treatment approaches for SLE in the future.
Contribution to the field
Systemic lupus erythematosus (SLE) is a chronic autoimmune disease of unknown etiology and mechanism. It has been shown that oral microbiota biomarkers could be served as a non-invasive diagnostic tool for many diseases, however, it is unclear whether oral microbiome could discriminate SLE from healthy controls. Compared with previous studies, we firstly analyzed characteristic alterations in the oral (tongue) microbiome of SLE in a large cohort, identified key bacteria associated with SLE and predicted specific microbial functions, and developed a first non-invasive microbial-based model for diagnosing SLE, validating its diagnostic efficiency for the first time in patients from multiple regions of China. Furthermore, we characterized the oral microbiome of treated SLE patients under stable conditions, which suggests that the oral microbiome may be associated with disease recovery in SLE. Our study suggests that oral microbiota dysbiosis may contribute to the development and progression of SLE. In addition, our study lays the foundation for future targeted oral microecological treatment of SLE.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the Sequence Read Archive database repository, https://www.ncbi.nlm.nih.gov/sra/.
Tsokos GC. Systemic lupus erythematosus. N Engl J Med. 2011;365:2110–21.
Rees F, Doherty M, Grainge MJ, Lanyon P, Zhang W. The worldwide incidence and prevalence of systemic lupus erythematosus: a systematic review of epidemiological studies. Rheumatology. 2017;56:1945–61.
Goldblatt F, O’Neill SG. Clinical aspects of autoimmune rheumatic diseases. Lancet. 2013;382:797–808.
Belkaid Y, Naik S. Compartmentalized and systemic control of tissue immunity by commensals. Nat Immunol. 2013;14:646–53.
Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.
Nenci A, Becker C, Wullaert A, Gareus R, van Loo G, Danese S, Huth M, Nikolaev A, Neufert C, Madison B, et al. Epithelial NEMO links innate immunity to chronic intestinal inflammation. Nature. 2007;446:557–61.
Ley RE, Peterson DA, Gordon JI. Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell. 2006;124:837–48.
Parks CG, de Santos ASE, Barbhaiya M, Costenbader KH. Understanding the role of environmental factors in the development of systemic lupus erythematosus. Best Pract Res Clin Rheumatol. 2017;31:306–20.
Guo M, Wang H, Xu S, Zhuang Y, An J, Su C, Xia Y, Chen J, Xu ZZ, Liu Q, et al. Alteration in gut microbiota is associated with dysregulation of cytokines and glucocorticoid therapy in systemic lupus erythematosus. Gut Microbes. 2020;11:1758–73.
Lippens C, Duraes FV, Dubrot J, Brighouse D, Lacroix M, Irla M, Aubry-Lachainaye JP, Reith W, Mandl JN, Hugues S. IDO-orchestrated crosstalk between pDCs and Tregs inhibits autoimmunity. J Autoimmun. 2016;75:39–49.
López P, de Paz B, Rodríguez-Carrio J, Hevia A, Sánchez B, Margolles A, Suárez A. Th17 responses and natural IgM antibodies are related to gut microbiota composition in systemic lupus erythematosus patients. Sci Rep. 2016;6:24072.
Li R, Meng X, Chen B, Zhao L, Zhang X. Gut microbiota in lupus: a butterfly effect? Curr Rheumatol Rep. 2021;23:27.
Rodríguez-Carrio J, López P, Sánchez B, González S, Gueimonde M, Margolles A, de Los Reyes-Gavilán CG, Suárez A. Intestinal dysbiosis is associated with altered short-chain fatty acids and serum-free fatty acids in systemic lupus erythematosus. Front Immunol. 2017;8:23.
Zhang L, Qing P, Yang H, Wu Y, Liu Y, Luo Y. Gut microbiome and metabolites in systemic lupus erythematosus: link mechanisms and intervention. Front Immunol. 2021;12:686501.
Hevia A, Milani C, López P, Cuervo A, Arboleya S, Duranti S, Turroni F, González S, Suárez A, Gueimonde M, et al. Intestinal dysbiosis associated with systemic lupus erythematosus. mBio. 2014;5:e01548-01514.
He Z, Shao T, Li H, Xie Z, Wen C. Alterations of the gut microbiome in Chinese patients with systemic lupus erythematosus. Gut Pathog. 2016;8:64.
Chen BD, Jia XM, Xu JY, Zhao LD, Ji JY, Wu BX, Ma Y, Li H, Zuo XX, Pan WY, et al. An autoimmunogenic and proinflammatory profile defined by the gut microbiota of patients with untreated systemic lupus erythematosus. Arthritis Rheumatol. 2021;73:232–43.
Li Y, Wang HF, Li X, Li HX, Zhang Q, Zhou HW, He Y, Li P, Fu C, Zhang XH, et al. Disordered intestinal microbes are associated with the activity of systemic lupus erythematosus. Clin Sci. 2019;133:821–38.
Azzouz D, Omarbekova A, Heguy A, Schwudke D, Gisch N, Rovin BH, Caricchio R, Buyon JP, Alekseyenko AV, Silverman GJ. Lupus nephritis is linked to disease-activity associated expansions and immunity to a gut commensal. Ann Rheum Dis. 2019;78:947–56.
Paster BJ, Boches SK, Galvin JL, Ericson RE, Lau CN, Levanos VA, Sahasrabudhe A, Dewhirst FE. Bacterial diversity in human subgingival plaque. J Bacteriol. 2001;183:3770–83.
Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, Wu X, Li J, Tang L, Li Y, et al. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat Med. 2015;21:895–905.
Xiao E, Mattos M, Vieira GHA, Chen S, Corrêa JD, Wu Y, Albiero ML, Bittinger K, Graves DT. Diabetes enhances IL-17 expression and alters the oral microbiome to increase its pathogenicity. Cell Host Microbe. 2017;22:120-128.e124.
Flemer B, Warren RD, Barrett MP, Cisek K, Das A, Jeffery IB, Hurley E, O’Riordain M, Shanahan F, O’Toole PW. The oral microbiota in colorectal cancer is distinctive and predictive. Gut. 2018;67:1454–63.
Ren Z, Wang H, Cui G, Lu H, Wang L, Luo H, Chen X, Ren H, Sun R, Liu W, et al. Alterations in the human oral and gut microbiomes and lipidomics in COVID-19. Gut. 2021;70:1253–65.
Li BZ, Zhou HY, Guo B, Chen WJ, Tao JH, Cao NW, Chu XJ, Meng X. Dysbiosis of oral microbiota is associated with systemic lupus erythematosus. Arch Oral Biol. 2020;113:104708.
van der Meulen TA, Harmsen HJM, Vila AV, Kurilshikov A, Liefers SC, Zhernakova A, Fu J, Wijmenga C, Weersma RK, de Leeuw K, et al. Shared gut, but distinct oral microbiota composition in primary Sjögren’s syndrome and systemic lupus erythematosus. J Autoimmun. 2019;97:77–87.
Corrêa JD, Calderaro DC, Ferreira GA, Mendonça SM, Fernandes GR, Xiao E, Teixeira AL, Leys EJ, Graves DT, Silva TA. Subgingival microbiota dysbiosis in systemic lupus erythematosus: association with periodontal status. Microbiome. 2017;5:34.
Hochberg MC. Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. 1997;40:1725.
Ren Z, Li A, Jiang J, Zhou L, Yu Z, Lu H, Xie H, Chen X, Shao L, Zhang R, et al. Gut microbiome analysis as a tool towards targeted non-invasive biomarkers for early hepatocellular carcinoma. Gut. 2019;68:1014–23.
Bombardier C, Gladman DD, Urowitz MB, Caron D, Chang CH. Derivation of the SLEDAI. A disease activity index for lupus patients. The Committee on Prognosis Studies in SLE. Arthritis Rheum. 1992;35:630–40.
Fanouriakis A, Kostopoulou M, Alunno A, Aringer M, Bajema I, Boletis JN, Cervera R, Doria A, Gordon C, Govoni M, et al. 2019 update of the EULAR recommendations for the management of systemic lupus erythematosus. Ann Rheum Dis. 2019;78:736–45.
Ren Z, Fan Y, Li A, Shen Q, Wu J, Ren L, Lu H, Ding S, Ren H, Liu C, et al. Alterations of the human gut microbiome in chronic kidney disease. Adv Sci. 2020;7:2001936.
Wei Y, Li Y, Yan L, Sun C, Miao Q, Wang Q, Xiao X, Lian M, Li B, Chen Y, et al. Alterations of gut microbiome in autoimmune hepatitis. Gut. 2020;69:569–77.
Chen C, Hemme C, Beleno J, Shi ZJ, Ning D, Qin Y, Tu Q, Jorgensen M, He Z, Wu L, Zhou J. Oral microbiota of periodontal health and disease and their changes after nonsurgical periodontal therapy. Isme j. 2018;12:1210–24.
Gergely P Jr, Grossman C, Niland B, Puskas F, Neupane H, Allam F, Banki K, Phillips PE, Perl A. Mitochondrial hyperpolarization and ATP depletion in patients with systemic lupus erythematosus. Arthritis Rheum. 2002;46:175–90.
Chodari L, Dilsiz Aytemir M, Vahedi P, Alipour M, Vahed SZ, Khatibi SMH, Ahmadian E, Ardalan M, Eftekhari A. Targeting mitochondrial biogenesis with polyphenol compounds. Oxid Med Cell Longev. 2021;2021:4946711.
Perl A, Hanczko R, Lai ZW, Oaks Z, Kelly R, Borsuk R, Asara JM, Phillips PE. Comprehensive metabolome analyses reveal N-acetylcysteine-responsive accumulation of kynurenine in systemic lupus erythematosus: implications for activation of the mechanistic target of rapamycin. Metabolomics. 2015;11:1157–74.
Shen P, Lin W, Deng X, Ba X, Han L, Chen Z, Qin K, Huang Y, Tu S. Potential implications of quercetin in autoimmune diseases. Front Immunol. 2021;12:689044.
Eftekhari A, Hasanzadeh A, Khalilov R, Hosainzadegan H, Ahmadian E, Eghbal MA. Hepatoprotective role of berberine against paraquat-induced liver toxicity in rat. Environ Sci Pollut Res Int. 2020;27:4969–75.
We thank all the generous volunteer subjects who enrolled in the study. We thank many clinical doctors from Department of Rheumatology and Immunology, the First Affiliated Hospital of Zhengzhou University, who were involved in this study. We thank Chao Liu from Shanghai Mobio Biomedical Technology.
This study was sponsored by grants from the National Key Research and Development Program of China (2018YFC2000500 and 2022YFC2303100), Research Project of Jinan Microecological Biomedicine Shandong Laboratory (JNL-2022001A and JNL-2022015B), and Henan Province Medical Science and Technology Program Provincial and Ministerial Joint Projects (SBGJ202003024).
Ethics approval and consent to participate
Ethical approval for this study was granted by the Ethics Committee of the First Affiliated Hospital of Zhengzhou University (2019-KY-200). Written informed consent of each participant was collected.
Consent for publication
All authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Guo, J., Cui, G., Huang, W. et al. Alterations in the human oral microbiota in systemic lupus erythematosus. J Transl Med 21, 95 (2023). https://doi.org/10.1186/s12967-023-03892-3
- Systemic lupus erythematosus
- Oral microbiota
- Diagnostic biomarker
- Disease activity