Alterations of the gut microbiome in Chinese patients with systemic lupus erythematosus
© The Author(s) 2016
Received: 26 August 2016
Accepted: 23 November 2016
Published: 8 December 2016
Systemic lupus erythematosus (SLE) in patients from Spain is associated with intestinal dysbiosis. This study explores whether the alteration of the gut microbiome in SLE patients from China is consistent with the intestinal dysbiosis of SLE patients from Spain.
The depletion of Firmicutes and the enrichment of Bacteroidetes in SLE patients from China were consistent with the SLE patients from Spain. Furthermore, we found that nine genera of gut microbiota were SLE-related microorganisms in Chinese subjects. Genera Rhodococcus, Eggerthella, Klebsiella, Prevotella, Eubacterium, Flavonifractor and Incertae sedis were significantly enriched, while genera Dialister and Pseudobutyrivibrio were significantly depleted in SLE patients. Receiver operating characteristic analysis indicated that the nine genera have the potential to distinguish SLE patients from healthy controls.
Comparing the dysbiosis of the gut microbiome among SLE patients from China or Spain, may indicate that the gut microbiome profiles of SLE patients are more influenced by disease than ethnicity.
KeywordsGut microbiome Systemic lupus erythematosus Chinese patients Illumina Miseq Genus level
The human gut harbors diverse and abundant microbes that benefit us by training the immune system , protecting against opportunistic pathogens , and harvesting nutrients and energy . Recently, increased evidence indicates that intestinal dysbiosis is associated with human diseases, such as cancer , liver cirrhosis , and autoimmune diseases [6, 7]. Systemic lupus erythematosus (SLE), a complex autoimmune disease characterized by its higher frequency in women, is associated with the gut microbiome. Compared to healthy controls, the SLE patients from Spain were characterized by increased Bacteroidetes levels and a lower Firmicutes/Bacteroidetes ratio . In addition, the species Lactobacillus was applied in management of SLE . Therefore, the gut microbiome can serve as new biomarkers or therapy for SLE.
It is well known that the human gut microbiome could be susceptible to the host’s genotype , age , sex  and diet . The association of the gut microbiome with diseases may be diverse due to the above factors. For example, the data of the association between the gut microbiome and obesity varied in different ethnicities [14–16]. The reason may be that pronounced differences in the gut microbiome existed among different ethnicities . China and Spain belong to two different continents with dissimilar genes and diets. Therefore, the alterations of the gut microbiome associated with SLE should be variable in SLE patients from China versus Spain.
In this study, Illumina Miseq sequencing was employed to investigate the alterations of the gut microbiome in SLE patients from China. The overall goal of the study was to ascertain whether the gut microbiome alteration of SLE in a Chinese population could be used as a biomarker and to evaluate differences compared to SLE biomarkers from Spain.
Demographic and clinical chemistry characteristics of human subjects
Age mean ± SD [min, max]
46.0 ± 1.8 [25, 61]
39.9 ± 4.3 [18, 62]
43.5 ± 2.4 [22, 68]
42.7 ± 1.9 [20, 56]
BMI (kg/m2) mean ± SD [min, max]
21.5 ± 0.6 [16.4, 28.8]
21.2 ± 1.2 [18.3, 27.7]
22.1 ± 1.0 [20.5, 27.7]
21.6 ± 0.2 [20.6, 25.8]
ESR (mm/h) ± SD [min, max]
14.8 ± 3.3 [3, 38]
14.6 ± 3.9 [8, 46]
6.9 ± 0.6 [4, 15]
7.2 ± 0.4 [3, 12]
SLEDAI ± SD [min, max]
7.5 ± 0.5 [3, 14]
6.7 ± 0.8 [4, 10]
Disease duration (years) ± SD [min, max]
7.9 ± 1.2 [1, 28]
5.0 ± 1.6 [1, 16]
Illumina Miseq sequencing of 16S rRNA gene-based amplicons and data processing
Total DNA was extracted from thawed fecal samples using the QIAamp® Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer protocols. The V3–V4 regions of the bacterial 16S rRNA gene sequences were amplified from the diluted DNA extracts with the primers 319F (5′-ACTCCTACGG GAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). PCR amplification was performed in a 30 μl mixture containing 0.5 μl of DMSO, 1.0 μl of forward primer (10 mM), 1.0 μl of reverse primer (10 mM), 5.0 μl of DNA sample, 7.5 μl of ddH2O and 15.0 μl of Phusion High-Fidelity PCR Master Mix with HF Buffer (NEB). The reactions were hot-started at 98 °C for 30 s, followed by 30 cycles of 98 °C for 15 s, 58 °C for 15 s, and 72 °C for 15 s, with a final extension step at 72 °C for 1 min. PCR products were purified using a QIAquick Gel Extraction kit (Qiagen, Valencia, CA, USA). The amplicon library was prepared using a TruSeq™ DNA sample preparation kit (Illumina Inc, San Diego, CA). The sequencing reaction was conducted using Illumina MiSeq platforms (Illumina Inc, San Diego, CA).
After sequencing, the data were analyzed on the Quantitative Insights Into Microbial Ecology (QIIME, www.qiime.org) platform using the default parameters . Before assembly, sequence reads were first filtered to remove low-quality or ambiguous reads, including reads lacking exact matching with the primer, reads containing ambiguous character (N), and reads with an average quality score <25. Only two reads with a sequence overlap longer than 20 bp were assembled. The assembled sequence reads with <400 bp or >500 bp were discarded. High-quality sequences were binned into 16S rRNA Operational Taxonomic Units (OTUs) and defined at ≥97% sequence homology. Chimera detection and removal were assessed using the GOLD reference databased and uchime . Taxonomic affiliation of each OTU was performed with QIIME against the SILVA database . Alpha diversity (PD_whole_tree, observed species, Simpson, Shannon, Singles, Doubles) and beta diversity (unweighted UniFrac) were performed in QIIME to 1000 reads, according to the OTU table rarefied and 20,000 reads were finally extracted from each sample for other analyses.
To test whether gut microbial species could be differentiated between SLE patients and healthy controls, a metric multidimensional scaling method based on projection known as principal coordinates analysis (PCoA) was used. Each sample was mapped based on the overall microbial composition and assessed for similarities.
Two datasets were created from all study subjects to investigate whether microbial biomarkers were available for analysis. The first dataset (train) was used to create the diagnostic algorithm. The above algorithm was then applied to the second dataset (test) for the purpose of evaluating how accurately microbial biomarkers were identified. All study subjects were randomly and independently divided into train or test datasets.
To visualize the potential of significantly altered genera as SLE biomarkers, receiver operating characteristic (ROC) curves were plotted in the statistical programming language R (V.3.1.2). ROC curves are an effective method of evaluating the quality or performance of diagnostic tests and are widely used in gut microbiomes to evaluate the performance of many microbial biomarkers. Area under the curves (AUCs) of ROC was calculated to evaluate the performance of the fitted logistic regression models. The AUCs were based on the predicted probability of SLE for each individual, using the multivariate logistic regression coefficient estimates together with the individual’s transformed relative abundances for each bacterial taxon included in the analysis.
Microbial phylum and genera concentrations between SLE patients and healthy controls were compared using a Mann–Whitney non-parametric test in SPSS software 16.0.
Structure and composition of the gut microbiome in SLE
SLE-associated microbiota changes
We next investigated whether the relative abundances of the gut microbiome differed between SLE patients and healthy controls at the genus level (Fig. 3). In the phylum Bacteroidetes, only the levels of the genera Prevotella were higher in SLE patients than in healthy controls. Five significantly altered levels of genera were found in phylum Firmicutes; genera Pseudobutyrivibrio and Dialister decreased in SLE patients, while Eubacterium, Flavonifractor, and Incertae sedis increased in SLE patients. Genera Klebsiella of the phylum Proteobacteria was found in lower amounts in healthy controls than in SLE patients. Genera Rhodococcus and Eggerthella of the phylum Actinobacteria were more abundant in SLE patients.
Evaluation of the performance of significantly changed genera
The clinical and immunological abnormalities of SLE among different ethnicities were known to widely vary because of diverse pathogenesis (genetic, environmental, hormonal, etc.,) . Therefore, the alterations of gut microbiome associated with SLE should display differences among separate ethnicities.
Firstly, a significant depletion of Firmicutes and enrichment of Bacteroidetes in SLE patients were found both in Spanish and Chinese subjects. Phylum Firmicutes and Bacteroidetes account for more than 90% of all phylogenetic species and are involved in host metabolism and immunity [1, 23]. Previous studies have described that the depletion of Firmicutes and the enrichment of Bacteroidetes in the human gut microbiota were associated with some disorders . Other studies also showed that an altered balance between Bacteroidetes and Firmicutes was correlated with clinical states . Therefore, the influence of SLE on human gut microbiome at the phylum level was not characterized by ethnicity. Therefore, the disorders of SLE are strongly related to the depletion of Firmicutes and the enrichment of Bacteroidetes in Spanish and Chinese patients.
Secondly, certain differences in gut microbiome were found among SLE patients from Spain and China. At the family level, SLE patients from Spain showed a depletion of Lachnospiraceae and Ruminococcaceae and an enrichment of Bacteroidaceae and Prevotellaceae . Our study indicated that only the family Prevotellaceae showed significant increases in SLE patients versus healthy controls from China (data not shown). However, this study employed Illumina Miseq sequencing to uncover some differences at the genus level between SLE patients and healthy controls from China. Significantly altered genera in SLE patients of China were mostly related to the intestinal disturbance of autoimmune disorders. Genus Rhodococcus infection was observed with Crohn’s disease . Genera Eggerthella, Klebsiella, Prevotella and Dialister were strongly positive for rheumatoid arthritis or ankylosing spondylitis [26–29]. Genus Pseudobutyrivibrio, reduced in SLE patients, displayed decreased levels in the gut microbiota profile of psoriatic arthritis . Genus Eubacterium was enriched in SLE patients of China, but decreased in multiple sclerosis patients . The disorders of the gut microbiome at the genera level in SLE patients from China were similar to the gut microbiome profiles of other autoimmune diseases. Thus, a combination of these genera could distinguish SLE patients from healthy controls.
As mentioned above, similarities at the phylum level and differences at family and genus levels of SLE-related gut microbiota existed between Spanish and Chinese subjects. The sensitivity of different taxon levels in evaluating microbial community may result in the above phenomenon. The gut microbiota appears to impose a greater influence on the detailed taxon of microorganisms. Therefore, the correlation of SLE with gut microbiota at the phylum level was more stable than that at the family or genera level. However, the differences of SLE-related gut microbiota at family and genera levels between Spanish and Chinese populations may be attributed to the differences of the host’s genetics and diet among the two ethnic groups.
Alterations of the gut microbiome in SLE patients from China and Spain were shown to be consistent at the phylum level with little differences at the family level. However, our study further revealed an alteration of the gut microbiome in SLE patients at the genus level. Significantly altered genera could thus have the ability to discriminate between SLE patients and healthy controls.
CW carried out study conception and design. ZH carried out acquisition of data, participated in the analysis and interpretation of data and drafted the manuscript. TS participated in the analysis and interpretation of data. HL, ZX participated in stool collection from the participants of this work. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and material
Please contact author for data requests.
Ethics approval and consent to participate
All procedures performed in this study involving human participants were in accordance with the ethical standards of the ethics Committee of Zhejiang Chinese Medical University. Informed consent to participate and publish was obtained from all individual participants included in the study.
This study was financially supported by National Natural Science Foundation of China (Nos. 81302936 and 81173249).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Round JL, Mazmanian SK. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol. 2009;9:313–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Kamada N, Chen GY, Inohara N, Núñez G. Control of pathogens and pathobionts by the gut microbiota. Nat Immunol. 2013;14:685–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Flint HJ, Scott KP, Louis P, Duncan SH. The role of the gut microbiota in nutrition and health. Nat Rev Gastroenterol Hepatol. 2012;9:577–89.View ArticlePubMedGoogle Scholar
- Zackular JP, Rogers MA, Ruffin MT, Schloss PD. The human gut microbiome as a screening tool for colorectal cancer. Cancer Prev Res. 2014;7:1112–21.View ArticleGoogle Scholar
- Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, Guo J, Le Chatelier E, Yao J, Wu L. Alterations of the human gut microbiome in liver cirrhosis. Nature. 2014;13:59–64.View ArticleGoogle Scholar
- Costello ME, Ciccia F, Willner D, Warrington N, Robinson PC, Gardiner B, Marshall M, Kenna TJ, Triolo G, Brown MA. Brief report: intestinal dysbiosis in ankylosing spondylitis. Arthritis Rheumatol. 2015;67:686–91.View ArticleGoogle Scholar
- Guo Z, Zhang J, Wang Z, Ang KY, Huang S, Hou Q, Su X, Qiao J, Zheng Y, Wang L. Intestinal microbiota distinguish gout patients from healthy humans. Sci Rep. 2016;6:20602.View ArticlePubMedPubMed CentralGoogle Scholar
- Hevia A, Milani C, López P, Cuervo A, Arboleya S, Duranti S, Turroni F, González S, Suárez A, Gueimonde M. Intestinal dysbiosis associated with systemic lupus erythematosus. Mbio. 2014;5:e01548-01514.View ArticleGoogle Scholar
- Mu Q, Zhang H, Luo XM. SLE: another autoimmune disorder influenced by microbes and diet? Front Immunol. 2015;6:608.View ArticlePubMedPubMed CentralGoogle Scholar
- Benson AK, Kelly SA, Legge R, Ma F, Low SJ, Kim J, Zhang M, Oh PL, Nehrenberg D, Hua K. Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. Proc Natl Acad Sci USA. 2010;107:18933–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP. Human gut microbiome viewed across age and geography. Nature. 2012;486:222–7.PubMedPubMed CentralGoogle Scholar
- Markle JG, Frank DN, Mortin-Toth S, Robertson CE, Feazel LM, Rolle-Kampczyk U, von Bergen M, McCoy KD, Macpherson AJ, Danska JS. Sex differences in the gut microbiome drive hormone-dependent regulation of autoimmunity. Science. 2013;339:1084–8.View ArticlePubMedGoogle Scholar
- Xu Z, Knight R. Dietary effects on human gut microbiome diversity. Br J Nutr. 2015;113:S1–5.View ArticlePubMedGoogle Scholar
- Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444:1022–3.View ArticlePubMedGoogle Scholar
- Patil DP, Dhotre DP, Chavan SG, Sultan A, Jain DS, Lanjekar VB, Gangawani J, Shah PS, Todkar JS, Shah S. Molecular analysis of gut microbiota in obesity among Indian individuals. J Biosci. 2012;37:647–57.View ArticlePubMedGoogle Scholar
- Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto J-M. Enterotypes of the human gut microbiome. Nature. 2011;473:174–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Tan EM, Cohen AS, Fries JF, Masi AT, Mcshane DJ, Rothfield NF, Schaller JG, Talal N, Winchester RJ. The 1982 revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheumatol. 1982;25:1271–7.View ArticleGoogle Scholar
- Petri M, Orbai AM, Alarcón GS, Gordon C, Merrill JT, Fortin PR, Bruce IN, Isenberg D, Wallace DJ, Nived O. Derivation and validation of the Systemic Lupus International Collaborating Clinics classification criteria for systemic lupus erythematosus. Arthritis Rheumatol. 2012;64:2677–86.View ArticleGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.View ArticlePubMedGoogle Scholar
- Li WG, Ye ZZ, Yin ZH, Zhang K. Clinical and immunological characteristics in 552 systemic lupus erythematosus patients in a southern province of China. Int J Rheum Dis. 2015. doi:10.1111/1756-185X.12480 [Epub ahead of print].Google Scholar
- Tremaroli V, Backhed F. Functional interactions between the gut microbiota and host metabolism. Nature. 2012;489:242–9.View ArticlePubMedGoogle Scholar
- Maranduba CMD, De Castro SBR, de Souza GT, Rossato C, da Guia FC, Valente MAS, Rettore JVP, Maranduba CP, de Souza CM, do Carmo AMR, et al. Intestinal microbiota as modulators of the immune system and neuroimmune system: impact on the host health and homeostasis. Clin Dev Immunol. 2015;2015:931574.Google Scholar
- Guerrero C, Tort J, Perez J, Andres M, Espejo E. Rhodococcus equi infection in a patient with Crohn’s disease treated with infliximab. J Infect. 2015;70:689–90.View ArticlePubMedGoogle Scholar
- Chen J, Wright K, Davis JM, Jeraldo P, Marietta EV, Murray J, Nelson H, Matteson EL, Taneja V. An expansion of rare lineage intestinal microbes characterizes rheumatoid arthritis. Genome Med. 2016;8:1.View ArticleGoogle Scholar
- Ebringer R, Cooke D, Cawdell DR, Cowling P, Ebringer A. Ankylosing spondylitis: klebsiella and HL-A B27. Rheumatol Rehabil. 1977;16:190–6.View ArticlePubMedGoogle Scholar
- Scher JU, Sczesnak A, Longman RS, Segata N, Ubeda C, Bielski C, Rostron T, Cerundolo V, Pamer EG, Abramson SB, et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. Elife. 2013;2:e01202.View ArticlePubMedPubMed CentralGoogle Scholar
- Tito RY, Cypers H, Joossens M, Varkas G, Van Praet L, Glorieus E, Van den Bosch F, De Vos M, Raes J. Dialister as microbial marker of disease activity in spondyloarthritis. Arthritis Rheumatol. 2016.
- Scher JU, Ubeda C, Artacho A, Attur M, Isaac S, Reddy SM, Marmon S, Neimann A, Brusca S, Patel T, et al. Decreased bacterial diversity characterizes the altered gut microbiota in patients with psoriatic arthritis, resembling dysbiosis in inflammatory bowel disease. Arthritis Rheumatol. 2015;67:128–39.View ArticlePubMedPubMed CentralGoogle Scholar
- Tremlett H, Fadrosh D, Lynch S, Hart J, Graves J, Lulu S, Aaen G, Belman A, Benson L, Casper C, et al. Gut microbiome in early pediatric multiple sclerosis: a case–control study. Mult Scler. 2014;20:338.View ArticleGoogle Scholar