Abstract
Little is known regarding differences in the gut microbiomes of rheumatoid arthritis (RA) patients and healthy cohorts in China. This study aimed to identify differences in the fecal microbiomes of 66 Chinese patients with RA and 60 healthy Chinese controls. The V3-V4 variable regions of bacterial 16S rRNA genes were sequenced with the Illumina system to define the bacterial composition. The alpha-diversity index of the microbiome of the RA patients was significantly lower than that of the control group. The bacterial genera Bacteroides (p = 0.02202) and Escherichia-Shigella (p = 0.03137) were more abundant in RA patients. In contrast, Lactobacillus (p = 0.000014), Alloprevotella (p = 0.0000008615), Enterobacter (p = 0.000005759), and Odoribacter (p = 0.0000166) were less abundant in the RA group than in the control group. Spearman correlation analysis of blood physiological measures of RA showed that bacterial genera such as Dorea and Ruminococcus were positively correlated with RF-IgA and anti-CCP antibodies. Furthermore, Alloprevotella and Parabacteroides were positively correlated with the erythrocyte sedimentation rate, and Prevotella-2 and Alloprevotella were positively correlated with C-reactive protein, both biomarkers of inflammation. These findings suggest that the gut microbiota may contribute to RA development via interactions with the host immune system.
Introduction
Rheumatoid arthritis (RA) is a common, chronic, autoimmune, and inflammatory disease affecting up to 1% of adults worldwide. It is characterized by progressive disability, systemic complications, early death, and socioeconomic costs (Firestein, ; Cooles and Isaacs, ). The pathogenesis of RA is complex, and the cause remains obscure.
Theories of the relationship between microbiome dysbiosis and RA date from the toxemic factor hypothesis originally proposed at the turn of the twentieth century. This hypothesis posited that the overgrowth of gram-negative bacteria in the intestines leads to an increase in toxic metabolites that enter the circulation and ultimately promote joint inflammation (Brusca et al., ). It is controversial whether RA is initiated by pathogenic bacteria (Tung et al., ; Scher and Abramson, ). Many bacteria have been proposed to be related to the morbidity of RA, including Mycoplasma fermentans (Sato et al., ), Escherichia coli (Newkirk et al., ), and Proteus mirabilis (Ebringer et al., , ). Mycbacterium avium subspecies paratuberculosis (MAP), and mycobacterial antigens (Sharp et al., ), Some metabolites of these bacteria have been linked to RA. For instance, antibodies to Proteus urease (IRRET) and Proteus hemolysin (ESRRAL) might act as autoantibodies to hyaline cartilage (LRREI) and HLA-DR1/DR4 (EQRRAA), respectively (Ebringer and Rashid, ). Intriguingly, epidemiological and pathogenetic relationships between periodontitis and RA have emerged, mainly associated with Porphyromonas gingivalis, which has been associated with periodontitis, autoantibody formation, and joint inflammation (Mikuls et al., ; Hitchon et al., ; Scher and Abramson, ). In addition to certain infectious agents, intestinal commensal microbes may also be related to RA.
Trillions of bacteria inhabit the human body. Through various mechanisms, they interact with the physical nervous system, endocrine system, and immune system. Colonization with autoantigen ortholog-producing commensal bacteria may initiate and sustain chronic autoimmunity in genetically predisposed populations (Greiling et al., ). Notably, commensal bacteria appear to be important in suppressing the inflammatory response and promoting immunological tolerance. The microbiota can regulate the differentiation of lymphocyte subsets; moreover, the microbiota, or specific microbes, can affect immune functions. Bacteroides fragilis (Mazmanian et al., ; Round and Mazmanian, ) and Faecalibacterium prausnitzii (Sokol et al., ), which exist in the human intestine as members of the commensal microbiota, induce CD4+ T cells that secrete interleukin-10 (IL-10). Moreover, polysaccharide A of B. fragilis promotes immunological tolerance by activating Foxp3+ regulatory T cells and suppresses T helper 17 (Th17) cell responses via a Toll-like receptor-2-dependent pathway (Round et al., ). The colonization of the small intestines of mice with segmented filamentous bacteria (SFB) induces the appearance of CD4+ T helper cells, promotes the production of Th17 cells, and induces autoimmune diseases, such as arthritis, in susceptible mice (Ivanov et al., ). SFB colonization of the intestinal epithelium induces the expression of serum amyloid A (SAA) protein, which in turn stimulates lamina propria dendritic cells (DCs) to produce IL-6 and IL-23, thus inducing the differentiation of Th17 cells (Ivanov et al., ). Notably, SAA has been shown to play a role in chronic inflammation in patients with RA (Chambers et al., ), indicating that SFB may contribute to the interplay among the microbiome, Th17 cells, and autoimmune function. In addition, Lactobacillus is a probiotic genus of bacteria that reside in the intestine as commensal microbes and exert immunomodulatory functions.
Perturbation of the gut microbiome contributes to dysbiosis and leads to diseases, including immune dysfunction. The pathogenesis of autoimmune diseases is associated with the absence of tolerance checkpoints that normally prevent self-antigens from stimulating the relentless growth of self-reactive B and T lymphocytes (Goodnow, ). CD4+CD25+ regulatory T cells are crucial in regulating self-reactive T cells and preventing autoimmune disease (Kojima et al., ; Tung et al., ). The Single Nucleotide Polymorphisms (SNPs) in the negative regulators Protein Tyrosine Phosphatase Non-receptor type 2 and 22 (PTPN2/22) may lead a dysregulated immune response, and triggers continuing apoptosis in chronic inflammation in RA (Sharp et al., ). Recently, it has been reported that the single nucleotide polymorphisms (SNPs) of STAT4, PTPN2, PSORS1C1, and TRAF3IP2RA genes are associated with the clinical efficacy of tumor necrosis factor (TNF) inhibitors in the treatment of rheumatoid arthritis (RA) patients (Yang et al., ). Meanwhile, the activation of self-reactive CD4+ T cells is necessary, although not sufficient, for inducing autoimmune diseases. HLA-II molecules, especially in activating CD4+ T cells, are critical in triggering human immune response (Huang et al., ). Therefore, the regulation of T cell subsets is crucial in maintaining immune balance. The disturbed gut microbiome may also promote the expression of IL-23. Spore-forming species, especially Clostridium clusters IV and XIVa, enhance the expression of transforming growth factor-beta (TGF-β), indoleamine 2,3-dioxygenase, and Foxp3 and have been implicated in the induction or accumulation of regulatory T cells (Tregs) in the colonic mucosa (Atarashi et al., ; Chiba and Seno, ). TGF-β can increases CD39 expression on Tregs via the activation of TGFBRII/TGFBRI, SMAD2, and the transcription factor CREB, which is activated in a p38-dependent manner and induces CD39 expression by promoting ENTPD1 gene transcription (Peres et al., ). Regulatory T (Treg) cells which play a key role in the modulation of immune responses have an impaired function in RA. Foxp3 is a master regulator of Treg cells which its expression is under the tight control of epigenetic mechanisms (Parisa et al., ).
Many studies in the global community have focused on the connection between the gut microbiome and RA in recent years. These have revealed the depletion of Bifidobacteria, the Bacteroides-Porphyromonas-Prevotella group, the Bacteroides fragilis subgroup, and the Eubacterium rectale-Clostridium coccoides group in the microbiome of RA patients compared with that in patients with non-inflammatory fibromyalgia (Vaahtovuo et al., ). In addition, Lactobacillus communities are more abundant in early RA (Liu et al., ). However, there is limited knowledge of the differences in the gut microbiomes of RA patients and healthy cohorts in China. In this study, we recruited 126 individuals for analysis of the microbial composition in RA patients in order to identify a feasible diagnostic strategy for RA based on the fecal microbiome.
Materials and Methods
Ethics Statement
This study was approved by the Institutional Review Board (IRB) at Shanghai GuangHua Hospital of Integrated Traditional Chinese and Western Medicine. All participants gave written informed consent in line with the Declaration of Helsinki.
Recruitment of Subjects and Sample Collection
Sixty-six patients with RA were recruited from Guanghua Hospital in Shanghai, China. Sixty healthy individuals also agreed to participate. RA patients fulfilled the RA criteria of the American College of Rheumatology. Fecal samples were obtained from all study subjects. All samples were transported in liquid nitrogen and stored at −80°C until the extraction of bacterial DNA.
DNA Extraction and Amplification of V3-V4 Regions
Total DNA was extracted from fecal samples using the QIAamp DNA Stool Mini Kit (Qiagen, Duesseldorf, Germany) following the manufacturer's instructions (QIA amp DNA Stool Handbook 04/2010). The 16S rRNA high-throughput sequencing was performed by Majorbio (Shanghai, China) with the MiSeq instrument (Illumina Inc., San Diego, CA, USA). Variable regions V3-V4 of the 16S rRNA genes of bacteria were amplified with barcode-indexed primers 338F and 806R. The forward primer sequence was 5′-ACTCCTACGGGAGGCAGCA-3′, and the reverse primer sequence was 5′-GGACTACHVGGGTWTCTAAT-3′.
MiSeq Sequencing
Amplicons were purified by 2% gel extraction (AxyPrep DNA Gel Extraction Kit, Axygen Biosciences, Union City, CA, USA) according to the manufacturer's instructions. The concentrations of polymerase chain reaction (PCR) products were measured using QuantiFluorTM-ST (Promega, Madison, WI, USA). Products were normalized at equimolar concentrations and subjected to paired-end sequencing (2 × 250) on the Illumina MiSeq platform (Illumina Inc.) according to standard protocols.
OTU Clustering and Annotation
After sequencing, the raw paired-end reads with overlapping nucleotides were assembled by pandaseq. Subsequently, the reads were quality-filtered. OTUs were determined at a 97% similarity cut-off value, and clustering analysis was then conducted using Usearch (version 7.0, http://drive5.com/uparse/). The Silva database (Release128, http://www.arb-silva.de) was used to obtain specific taxonomic information corresponding to each OTU.
Fecal Microbiome Analysis
The α-diversity was measured by species richness from the rarefied OTU table, and β-diversity was estimated by computational analysis using unweighted UniFrac scores and was visualized with principal coordinate analysis (Tang et al., ). Hierarchical cluster analysis was performed with single-linkage, complete-linkage, and average-linkage clustering methods. The LEfSe was used to elucidate differences in bacterial taxa. An LDA score ≥ 2 indicated that the bacterial group was an important contributor to RA. LEfSe analysis was performed online in the Galaxy workflow framework (http://huttenhower.sph.harvard.edu/galaxy/root?tool_id=lefse_upload).
Quantitative Analysis of 16S rRNA Genes
Quantitative PCR (qPCR) was validated by constructing an artificial mixture of 16S rRNA genes according to published approaches (Becker et al., ; Bacchetti De Gregoris et al., ). The forward primer (5′-CATGTGGTTTAATTCGATGAT-3′) and the reverse primer (5′-AGCTGACGACAACCATGCAG-3′) were designed to detect Bacteroides.
Statistical Analyses
All statistical analyses were performed using the R package (v. 3.2.2). Student's t-tests were performed to assess differences in the α-diversity of the two groups. For the comparison of continuous variables, the Mann-Whitney U-test was conducted for two groups. Multiple test correction was performed using the Benjamini and Hochberg false discovery rate. Spearman's rank test was performed for correlation analysis.
Results
Sample Size and Sequencing Depth
A total of 126 samples from 66 RA patients and 60 healthy controls were analyzed in this study. The age and sex of all subjects are described in Table 1, and the sociodemographic details and pathological status of the RA patients are provided in Table 2. Detailed demographic and clinical characteristics of the study cohorts with RA are shown in Supplementary Table 1. Given a fact that RA occurs more frequently in women, female patients are more than male patients in this study. From these, 7,074,888 high-quality sequences with a median read length of 436 bp were acquired. Pan and core species analysis revealed increases in the number of species with increased sample sizes (Figure 1).
Table 1
| Group | Gender | Age | Number of subjects |
|---|---|---|---|
| RA | Female | 21–29 | 1 |
| (N = 66) | (N = 51) | 30–39 | 8 |
| 40–49 | 9 | ||
| 50–59 | 14 | ||
| 60–69 | 16 | ||
| 70–79 | 2 | ||
| >=80 | 1 | ||
| Male | 50–59 | 6 | |
| (N = 15) | 60~69 | 9 | |
| Control | Female | 20–29 | 5 |
| (N = 60) | (N = 31) | 30–39 | 2 |
| 40–49 | 11 | ||
| 50–59 | 10 | ||
| 60–69 | 3 | ||
| Male | 20–29 | 9 | |
| (N = 29) | 30–39 | 4 | |
| 40–49 | 6 | ||
| 50–59 | 6 | ||
| >=60 | 4 |
Age and sex characteristics of study subjects.
Table 2
| Variables | Values |
|---|---|
| Male n (%) | 15 (22.7%) |
| Age (years) | 54.95 ± 1.4 |
| Disease duration (years) | 32 ± 0.6 (male patients) |
| 33.8 ± 0.2 (female patients) | |
| Inflammatory markers | CR > 133 μmol/L (male patients) |
| CR > 106 μmol/L (female patients) | |
| ESR > 15 mm/h (male patients) | |
| ESR > 20 mm/h (female patients) | |
| CRP > 10 mg/L | |
| Serology | RF > 20 |
| Anti-CCP > 25 | |
| Joint involvement | 15.9 ± 14.9 (male patients) |
| 24.5 ± 17.4 (female patients) |
Sociodemographic characteristics and pathological status of the RA group.
Values are shown as mean ± SD.
SD, standard deviation; CR, Creatinine; ESR, Erythrocyte sedimentation rate; CRP, C-reactive protein; RF, Rheumatoid factor; CCP, Cyclic citrullinated peptide.
Figure 1
α-Diversity of the Gut Microbiome
The α-diversity, determined using the Chao1 and Sobs indexes, differed significantly between the RA group and the healthy control group (Figure 2). The gut microbiome was markedly more diverse in composition in the control group than in the RA group. The data suggest a relationship between a decline in the α-diversity of the gut microbiome and the development of RA.
Figure 2
β-Diversity of the Gut Microbiome
Principal coordinate analysis based on operational taxonomic unit (OTU) abundances provided an overview of the gut microbiome and reflected the β-diversities of the different groups. The β-diversity was clearly higher in the control group than in the RA group (Figure 3). Among the control group, the bacterial communities in some fecal samples were similar to those of patients with RA, while the microbial composition of other individuals in the control group was completely distinct from that of the RA group. Samples from both males and females were included in this analysis and were found to differ in β-diversity.
Figure 3
Microbiome Alterations
Bacteria in the phyla Bacteroidetes, Firmicutes, Proteobacteria, and Actinobacteria are the most enriched in the human intestine. The present data revealed the increased abundance of Bacteroidetes in RA patients, with decreased abundances of the other three bacterial phyla in the gut microbiome of RA patients (Figure 4A).
Figure 4
Of the total of 179 bacterial families, 126 were significantly different between the RA and control groups (p < 0.05; data not shown). Of the total of 438 bacterial genera, 221 differed markedly between the two groups (Supplementary Table 2). Significance tests were performed for the 20 most abundant bacterial genera that differed in abundance between the two groups, and the results are depicted in Figure 4B. RA group exhibited reductions in the genera Lactobacillus, Alloprevotella, Enterobacter, Clostridium sensu stricto-3, Ruminococcaceae UCG-014, Odoribacter, Rikenellaceae RC9, Enterococcus, Klebsiella, Desulfovibrio, Citrobacter, Akkermansia, Helicobacter, Rikenella, Staphylococcus, Coprococcus 1, Coriobacteriaceae UCG-002, and Rhodococcus. In contrast, nine bacterial genera, including Bacteroides, Escherichia-Shigella, Parasutterella, Flavonifractor, the Eubacterium xylanophium group, Tyzzerella, Sellimonas, and Oscillospira, were significantly increased in abundance in RA patients compared with that in healthy individuals (Supplementary Figure 1). The abundance of Bacteroides was confirmed by qPCR (Figure 4C). Circos plots were used to visualize the relationship between sample groups and specific bacterial families and genera (Supplementary Figure 2). Samples in the RA and control groups displayed similar bacterial members, although the abundances of many of these members differed substantially. Lactobacillaceae were evident in both groups, with the abundance being markedly higher in the control group.
Bacterial Taxa Related to RA
To identify the specific bacterial taxa associated with RA, the compositions of the gut microbiomes of RA patients and healthy individuals were compared using the LDA effect size (LEfSe) method. The resulting cladogram presents the structure of the gut microbiome and the predominant bacteria in the control and RA groups. In the RA group, several bacterial genera, including Bacteroides, Escherichia Shigella, Tyzzerella, Parasutterella, Flavonifractor, Eubacterium saphenum group, and Sellimonas, could be considered microbial markers, as shown in Supplementary Figure 3.
Functional Predictions Based on 16S rRNA Gene Sequencing
Inferred microbiome functions based on 16S rRNA gene sequences were also explored using the PICRUSt algorithm (Langille et al., ). Sixty Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, involving the biosynthesis of secondary metabolites, xenobiotics biodegradation and metabolism, metabolism of terpenoids and polyketides, immune system, and endocrine system, that differed in abundance between the RA group and the control group at q ≤ 0.05 are shown in the heatmap in Supplementary Figure 4. The implicated pathways did not obviously differ among the different bacterial phyla. Most pathways in the heatmap were enriched in the healthy controls but depleted in the RA group. Notably, the genus Bacteroides was depleted in most of these pathways. Interestingly, the pathway related to steroid hormone biosynthesis, which includes sex hormones, was enriched in Bacteroides but depleted in other bacterial genera.
Correlation Between Microbiota and Physiological Measures of RA
Spearman correlation analysis of rheumatoid factors (RFs), anti-CCP antibodies, inflammatory biomarkers, and the most abundant 50 bacterial genera among the intestinal microbiota was performed. As shown in Figure 5, positive correlations (correlation coefficient > 0) and negative correlations (correlation coefficient < 0) were observed between various bacteria and rheumatoid factors, anti-CCP antibodies, and inflammatory biomarkers.
Figure 5
Discussion
The differences in the gut microbiome between RA patients and healthy individuals have been demonstrated systematically since the beginning of this century. Nevertheless, details of the microbiome of RA patients in the Chinese population remained unclear. In the current study, we recruited a large cohort of 126 individuals to study the characteristics of the gut microbiome composition in RA patients in China.
The gut microbiome of these RA patients showed a decline in species richness as well as a moderate reduction in β-diversity, consistent with previous reports (Chen et al., ; Van de Wiele et al., ). The compositional characteristics of the gut microbiome in the RA group were like the overall characteristics of most samples in the control group but were distinct from some samples in the control group. This collection of control individuals contained both males and females with no clear defining characteristic, indicating that the composition of the gut microbiome differed markedly among those in the healthy cohort but not among patients with RA. The results indicate that many individuals with specific compositional characteristics of the gut microbiome are more susceptible to RA, which may aid in the diagnosis or determination of the susceptibility of individuals to RA via detection of the gut microbiome.
At the genus level, Lactobacillus was depleted in the RA group compared to its abundance in the control group. Lactobacilli are typical probiotic microorganisms found in the normal bacterial flora. They exhibit immunoregulatory functions in the host and play important roles in the maintenance of homeostasis in the intestinal microbiota by secreting immunomodulatory agents (Lorca et al., ; Zvanych et al., ). However, some species of Lactobacillus, such as L. bifidus, can cause joint swelling in mice (Abdollahi-Roodsaz et al., ). Several studies have reported a greater abundance of Lactobacillus in RA patients than in healthy controls in the gut microbiome and the oral and salivary microbiome (Liu et al., ; Van de Wiele et al., ). Conversely, in this study, the abundance of Lactobacillus was significantly reduced in RA patients compared to that in the control group. Many studies, including several randomized, double-blinded, placebo-controlled trials, have demonstrated that supplementation with L.GG, L. casei, or Lactobacillus-enriched vegetarian food can moderately improve health and attenuate inflammation in RA patients (Nenonen et al., ; Hatakka et al., ; de los Angeles Pineda et al., ; Alipour et al., ).
At the phylum level, RA patients displayed an increase in Bacteroidetes and decreases in Firmicutes, Proteobacteria, and Actinobacteria compared to levels in healthy individuals. Differences in the abundances of the genus Bacteroides between the RA group and control group accounted for the majority of differences within the phylum Bacteroidetes. Bacteroides species are reported to mediate chronic colitis, gastritis, and arthritis in HLA-B27 transgenic rats (Rath et al., ) and to cause osteomyelitis in RA (Crook and Gray, ). Studies in germ-free mice have demonstrated that the Bacteroides integrase antigen contributes to the recruitment and proliferation of low-avidity CD8+ T cells, which may be analogous to thymic CD4+ Tregs, and the response to chronic antigenic exposure in gut-associated lymphoid tissues (Hebbandi Nanjundappa et al., ). The enrichment of Bacteroidetes in RA patients in this study indicates that this condition may contribute to disease progression. Additionally, the expansion of Bacteroidetes may be attributed to a compensatory mechanism for downregulating autoimmune reactions.
The microbiota, especially certain microbes, can influence human physiology and pathology. Functional analysis of the KEGG pathways related to 16S rRNA gene sequences and differences in the microbiomes of RA patients and healthy controls provided evidence that some microbe-derived functions, including the biosynthesis of some secondary metabolites, degradation of xenobiotics, metabolism, and the endocrine system, may be depleted in patients with RA. However, these pathways were also depleted in Bacteroides, a bacterial genus enriched in the RA group. These findings indicate that the composition of the gut microbiota may contribute to the initiation or development of RA through altered microbial function. Bacteroides might play a crucial role in this process.
Moreover, there were significant relationships between physiological measures of RA and the microbiota, as shown by Spearman correlation analysis. For example, Sutterella, Megamonas, and Faecalibacterium were negatively correlated with RF levels. Other bacteria, such as Dorea and Alloprevotella, were positively correlated with the RF levels or anti-CCP antibodies. This indicates that specific bacteria may contribute to symptoms of RA. However, it is difficult to determine the exact influence of the bacterial population on RA, as isolation and in vitro culture of most microbial strains is challenging.
The development of RA is associated with gender. The morbidity of RA is higher in women than in men (van Vollenhoven, ). Previous studies demonstrated that sex hormones can induce autoimmune diseases through an effect on the function of the immune system (Ahmed et al., ). Many microbiota functions predicted using bacterial 16S rRNA gene sequencing involved the endocrine system and steroid hormones, which may play a role in RA. An experiment in mice showed that interferon-gamma is the predominant inflammatory cytokine after the activation of T cells in female mice, indicating that the Th1 response is more intense in female mice, while the Th2 response is dominant in male mice (Bebo et al., ). Female hormones can regulate the activity of T cells, B cells, natural killer cells, macrophages, endotheliocytes, and matrix cells. Moreover, T cells are the main target cells of female hormones. Thus, levels of female hormones can result in the dysbiosis of Tregs. However, women are not more susceptible to some other autoimmune diseases, including ankylosing spondylitis. A study of a non-obese mouse model of type 1 diabetes revealed that the sex-related microbiome can drive hormone-dependent regulation of autoimmunity (Calin et al., ). In the present study, LEfSe analysis of the gut microbiome compositions of males and females showed that Rikenellaceae, Porphyromonadaceae, Coriobacteriaceae, and certain genera were significantly more abundant in females, while Pasteurellaceae, Butyricicoccus, Clostridiaceae 1, Clostridium sensu stricto 1, Alisnonella, and other genera were more abundant in males according to linear discriminant analysis (LDA) scores (Supplementary Figure 5). The members of the Clostridiaceae enriched in males are reported to stimulate the differentiation of colonic regulatory T cells (Hasegawa and Inohara, ). However, we did not conclusively demonstrate that differences in the composition of the gut microbiome between males and females influence the development of RA in this study, which was limited by the sample size to some extent. Except for gender, given the complexity of clinical samples collection, we cannot rule out the influence of age on experimental results, although the results in this study is solid. We performed statistics on the ages of the RA group and the control group and found that the former was significantly higher than the latter (data not shown). Despite this, such age groups are closer to the real situation in clinical cases. Therefore, we believe that age and gender are the co-influence factors of the bacterial population of RA patients, but our results still have considerable guiding significance.
Studies of the gut microbiome in Chinese patients with RA can contribute to the recognition of the microbiome composition in patients with autoimmune diseases. Such altered bacterial collections can act as an adjuvant criterion for clinical diagnosis. Further research using different enterotypes and genders is needed, since the normal gut microbiome may exhibit marked differences in the enrichment of some bacteria, such as Bacteroides spp., Prevotella spp., and Ruminococcus spp. (Jeffery et al., ).
The contribution of the fungiome and virome to human health or disease has been highlighted recently. A study based on real-time PCR detected cytomegalovirus, Epstein–Barr virus, and human herpesvirus 6 in the sera and peripheral blood mononuclear cells of patients with RA (Alvarez-Lafuente et al., ). Intestinal fungiome and virome analyses are theoretically included in analyses of the microbiome, although the actual results may be limited by technological complexities. Further studies to reveal the detailed bacteriome, virome, and fungiome characteristics in RA patients with different enterotypes are therefore warranted.
Conclusions
In this study, we showed that the diversity and composition of the microbiome of RA patients differed from those in healthy control subjects in China. An increase in Bacteroides and Escherichia-Shigella, combined with decreases in Lactobacillus, Alloprevotella, Enterobacter, and Odoribacter appear to be characteristic of RA in patients from Shanghai, China. These findings differ from those of other studies performed outside of China. Spearman correlation analysis showed that Dorea and Ruminococcus were positively correlated with RF-IgA and anti-CCP antibodies. Similarly, Alloprevotella was positively correlated with rheumatoid factors (RF-IgM, RF-IgA, and RF-IgG) and inflammation biomarkers (CRP, ESR). The roles of these bacteria in the etiology and pathogenesis of RA require further exploration, as this analysis was limited by the small sample size. Future studies should also evaluate the significance of these microbial changes for the diagnosis, prognosis, and drug response of RA.
Statements
Data availability statement
The datasets generated for this study can be found in NCI SRA https://www.ncbi.nlm.nih.gov/sra/SRP223483.
Ethics statement
This study was approved by the Institutional Review Board (IRB) at Shanghai GuangHua Hospital of Integrated Traditional Chinese and Western Medicine. All participants gave written informed consent in line with the Declaration of Helsinki.
Author contributions
LX, Y-WT, ZP, HL, and XG designed the study. YS, QC, and PL conducted analyses and wrote the manuscript. YS, QC, LY, QL, and CL edited the manuscript. RX, DH, WJ, YB, YS, and KD provided sociodemographic and pathological data.
Funding
This work was supported by the Shanghai Changning District Health System Characteristic Specialized Construction Project Program (No. 20162003), Jieping Wu Medical Foundation (No. 320.6750.15232), and Science and Technology Commission Plan (No. CNKW2016Z04).
Acknowledgments
We thank Dr. Wenjuan Sun and Xiaodan Li (Shanghai Major Biological Engineering Technology & Service Co., Ltd., Shanghai, China) for providing excellent technical assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2019.00369/full#supplementary-material
References
1
Abdollahi-RoodsazS.JoostenL. A.KoendersM. I.DevesaI.RoelofsM. F.RadstakeT. R.et al. (2008). Stimulation of TLR2 and TLR4 differentially skews the balance of T cells in a mouse model of arthritis. J. Clin. Invest. 118, 205–216. 10.1172/JCI32639
2
AhmedS. A.PenhaleW. J.TalalN. (1985). Sex hormones, immune responses, and autoimmune diseases: mechanisms of sex hormone action. Am. J. Pathol. 121, 531–551.
3
AlipourB.Homayouni-RadA.Vaghef-MehrabanyE.SharifS. K.Vaghef-MehrabanyL.Asghari-JafarabadiM.et al. (2014). Effects of Lactobacillus casei supplementation on disease activity and inflammatory cytokines in rheumatoid arthritis patients: a randomized double-blind clinical trial. Int. J. Rheum. Dis. 17, 519–527. 10.1111/1756-185X.12333
4
Alvarez-LafuenteR.Fernandez-GutierrezB.de MiguelS.JoverJ. A.RollinR.LozaE.et al. (2005). Potential relationship between herpes viruses and rheumatoid arthritis: analysis with quantitative real time polymerase chain reaction. Ann. Rheum. Dis. 64, 1357–1359. 10.1136/ard.2004.033514
5
AtarashiK.TanoueT.ShimaT.ImaokaA.KuwaharaT.MomoseY.et al. (2011). Induction of colonic regulatory T cells by indigenous Clostridium species. Science331, 337–341. 10.1126/science.1198469
6
Bacchetti De GregorisT.AldredN.ClareA. S.BurgessJ. G. (2011). Improvement of phylum- and class-specific primers for real-time PCR quantification of bacterial taxa. J. Microbiol. Methods86, 351–356. 10.1016/j.mimet.2011.06.010
7
BeboB. F.SchusterJ. C.VandenbarkA. A.OffnerH. (1999). Androgens alter the cytokine profile and reduce encephalitogenicity of myelin-reactive T cells. J. Immunol. 162, 35–40.
8
BeckerS.BogerP.OehlmannR.ErnstA. (2000). PCR bias in ecological analysis: a case study for quantitative Taq nuclease assays in analyses of microbial communities. Appl. Environ. Microbiol. 66, 4945–4953. 10.1128/AEM.66.11.4945-4953.2000
9
BruscaS. B.AbramsonS. B.ScherJ. U. (2014). Microbiome and mucosal inflammation as extra-articular triggers for rheumatoid arthritis and autoimmunity. Curr. Opin. Rheumatol.26, 101–107. 10.1097/BOR.0000000000000008
10
CalinA.BrophyS.BlakeD. (2000). Impact of sex on inheritance of ankylosing spondylitis. Lancet355:1098. 10.1016/S0140-6736(05)72218-2
11
ChambersR.MacFarlaneD.WhicherJ.DieppeP. (1983). Serum amyloid-A protein concentration in rheumatoid arthritis and its role in monitoring disease activity. Ann. Rheum. Dis. 42, 665–667. 10.1136/ard.42.6.665
12
ChenJ.WrightK.DavisJ. M.JeraldoP.MariettaE. V.MurrayJ.et al. (2016). An expansion of rare lineage intestinal microbes characterizes rheumatoid arthritis. Genome Med. 8:43. 10.1186/s13073-016-0299-7
13
ChibaT.SenoH. (2011). Indigenous clostridium species regulate systemic immune responses by induction of colonic regulatory T cells. Gastroenterology141, 1114–1116. 10.1053/j.gastro.2011.07.013
14
CoolesF. A.IsaacsJ. D. (2011). Pathophysiology of rheumatoid arthritis. Curr. Opin. Rheumatol. 23, 233–240. 10.1097/BOR.0b013e32834518a3
15
CrookP. R.GrayJ. (1982). Bacteroides causing osteomyelitis in rheumatoid arthritis. Ann. Rheum. Dis. 41, 645–646. 10.1136/ard.41.6.645
16
de los Angeles PinedaM.ThompsonS. F.SummersK.de LeonF.PopeJ.ReidG. (2011). A randomized, double-blinded, placebo-controlled pilot study of probiotics in active rheumatoid arthritis. Med. Sci. Monit. 17:CR347. 10.12659/MSM.881808
17
EbringerA.CorbettM.MacafeeY.BaronP.PtaszynskaT.WilsonC.et al. (1985). Antibodies to proteus in rheumatoid arthritis. Lancet326, 305–307. 10.1016/S0140-6736(85)90352-6
18
EbringerA.RashidT. (2014). Rheumatoid arthritis is caused by a Proteus urinary tract infection. APMIS122, 363–368. 10.1111/apm.12154
19
EbringerA.RashidT.WilsonC. (2010). Rheumatoid arthritis, Proteus, anti-CCP antibodies and Karl Popper. Autoimmun. Rev. 9, 216–223. 10.1016/j.autrev.2009.10.006
20
FiresteinG. S. (2003). Evolving concepts of rheumatoid arthritis. Nature423, 356–361. 10.1038/nature01661
21
GoodnowC. C. (2007). Multistep pathogenesis of autoimmune disease. Cell130, 25–35. 10.1016/j.cell.2007.06.033
22
GreilingT. M.DehnerC.ChenX.HughesK.IniguezA. J.BoccittoM.et al. (2018). Commensal orthologs of the human autoantigen Ro60 as triggers of autoimmunity in lupus. Sci. Transl. Med. 10:eaan2306. 10.1126/scitranslmed.aan2306
23
HasegawaM.InoharaN. (2014). Regulation of the gut microbiota by the mucosal immune system in mice. Int. Immunol. 26, 481–487. 10.1093/intimm/dxu049
24
HatakkaK.MartioJ.KorpelaM.HerranenM.PoussaT.LaasanenT.et al. (2003). Effects of probiotic therapy on the activity and activation of mild rheumatoid arthritis–a pilot study. Scand. J. Rheumatol. 32, 211–215. 10.1080/03009740310003695
25
Hebbandi NanjundappaR.RonchiF.WangJ.Clemente-CasaresX.YamanouchiJ.Sokke UmeshappaC.et al. (2017). A gut microbial mimic that hijacks diabetogenic autoreactivity to suppress colitis. Cell171, 655–667.e17. 10.1016/j.cell.2017.09.022
26
HitchonC. A.ChandadF.FerucciE. D.WillemzeA.Ioan-FacsinayA.van der WoudeD.et al. (2010). Antibodies to Porphyromonas gingivalis are associated with anticitrullinated protein antibodies in patients with rheumatoid arthritis and their relatives. J. Rheumatol. 37, 1105–1112. 10.3899/jrheum.091323
27
HuangZ.NiuQ.YangB.ZhangJ.YangM.XuH.et al. (2018). Genetic polymorphism of rs9277535 in HLA-DP associated with rheumatoid arthritis and anti-CCP production in a Chinese population. Clin Rheumatol. 37, 1799–1805. 10.1007/s10067-018-4030-5
28
IvanovI. I.AtarashiK.ManelN.BrodieE. L.ShimaT.KaraozU.et al. (2009). Induction of intestinal Th17 cells by segmented filamentous bacteria. Cell139, 485–498. 10.1016/j.cell.2009.09.033
29
JefferyI. B.ClaessonM. J.O'TooleP. W.ShanahanF. (2012). Categorization of the gut microbiota: enterotypes or gradients?Nat. Rev. Microbiol. 10, 591–592. 10.1038/nrmicro2859
30
KojimaA.Tanaka-KojimaY.SakakuraT.NishizukaY. (1976). Prevention of postthymectomy autoimmune thyroiditis in mice. Lab. Invest. 34, 601–605.
31
LangilleM. G. I.ZaneveldJ.CaporasoJ. G.McdonaldD.KnightsD.ReyesJ.et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. 10.1038/nbt.2676
32
LiuX.ZouQ.ZengB.FangY.WeiH. (2013). Analysis of fecal Lactobacillus community structure in patients with early rheumatoid arthritis. Curr. Microbiol. 67, 170–176. 10.1007/s00284-013-0338-1
33
LorcaG. L.WadstromT.ValdezG. F.LjunghA. (2001). Lactobacillus acidophilus autolysins inhibit Helicobacter pylori in vitro. Curr. Microbiol. 42, 39–44. 10.1007/s002840010175
34
MazmanianS. K.LiuC. H.TzianabosA. O.KasperD. L. (2005). An immunomodulatory molecule of symbiotic bacteria directs maturation of the host immune system. Cell122, 107–118. 10.1016/j.cell.2005.05.007
35
MikulsT. R.PayneJ. B.ReinhardtR. A.ThieleG. M.MaziarzE.CannellaA. C.et al. (2009). Antibody responses to Porphyromonas gingivalis (P. gingivalis) in subjects with rheumatoid arthritis and periodontitis. Int. Immunopharmacol. 9, 38–42. 10.1016/j.intimp.2008.09.008
36
NenonenM.HelveT.RaumaA.HänninenO. (1998). Uncooked, lactobacilli-rich, vegan food and rheumatoid arthritis. Br. J. Rheumatol. 37, 274–281. 10.1093/rheumatology/37.3.274
37
NewkirkM. M.ZbarA.BaronM.MangesA. R. (2010). Distinct bacterial colonization patterns of Escherichia coli subtypes associate with rheumatoid factor status in early inflammatory arthritis. Rheumatology49, 1311–1316. 10.1093/rheumatology/keq088
38
ParisaZ.YariK.MostafaeiS.IranshahiN.AssarS.FekriA.et al. (2018). Analysis of helios gene expression and foxp3 tsdr methylation in the newly diagnosed rheumatoid arthritis patients. Immunol. Investig.10, 1–11. 10.1080/08820139.2018.1480029
39
PeresR. S.DonateP. B.TalbotJ.CecilioN. T.LoboP. R.MachadoC.aioC.et al. (2018). Tgf-β signalling defect is linked to low cd39 expression on regulatory t cells and methotrexate resistance in rheumatoid arthritis. J. Autoimmun.90, 49–58. 10.1016/j.jaut.2018.01.004
40
RathH. C.HerfarthH. H.IkedaJ. S.GrentherW. B.HammT. E.Jr.BalishE.et al. (1996). Normal luminal bacteria, especially Bacteroides species, mediate chronic colitis, gastritis, and arthritis in HLA-B27/human beta2 microglobulin transgenic rats. J. Clin. Invest. 98, 945–953. 10.1172/JCI118878
41
RoundJ. L.LeeS. M.LiJ.TranG.JabriB.ChatilaT. A.et al. (2011). The Toll-like receptor 2 pathway establishes colonization by a commensal of the human microbiota. Science332, 974–977. 10.1126/science.1206095
42
RoundJ. L.MazmanianS. K. (2010). Inducible Foxp3+ regulatory T-cell development by a commensal bacterium of the intestinal microbiota. Proc. Natl. Acad. Sci. U.S.A. 107, 12204–12209. 10.1073/pnas.0909122107
43
SatoN.OizumiT.KinbaraM.SatoT.FunayamaH.SatoS.et al. (2010). Promotion of arthritis and allergy in mice by aminoglycoglycerophospholipid, a membrane antigen specific to Mycoplasma fermentans. FEMS Immunol. Med. Microbiol. 59, 33–41. 10.1111/j.1574-695X.2010.00657.x
44
ScherJ. U.AbramsonS. B. (2011). The microbiome and rheumatoid arthritis. Nat. Rev. Rheumatol. 7, 569–578. 10.1038/nrrheum.2011.121
45
SharpR. C.BegS. A.NaserS. A. (2018). Polymorphisms in protein tyrosine phosphatase non-receptor type 2 and 22 (PTPN2/22) are linked to hyper-proliferative T-cells and susceptibility to mycobacteria in rheumatoid arthritis. Front Cell Infect Microbiol.8:11. 10.3389/fcimb.2018.00011
46
SokolH.PigneurB.WatterlotL.LakhdariO.BermudezhumaranL. G.GratadouxJ.et al. (2008). Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc. Natl. Acad. Sci. U.S.A. 105, 16731–16736. 10.1073/pnas.0804812105
47
TangR.WeiY.LiY.ChenW.ChenH.WangQ.et al. (2018). Gut microbial profile is altered in primary biliary cholangitis and partially restored after UDCA therapy. Gut67, 534–541. 10.1136/gutjnl-2016-313332
48
TungK.SmithS.MatznerP.KasaiK.OliverJ.FeuchterF.et al. (1987). Murine autoimmune oophoritis, epididymoorchitis, and gastritis induced by day 3 thymectomy. Autoantibodies. Am. J. Pathol. 126:303.
49
VaahtovuoJ.MunukkaE.KorkeamakiM.LuukkainenR.ToivanenP. (2008). Fecal microbiota in early rheumatoid arthritis. J. Rheumatol. 35, 1500–1505.
50
Van de WieleT.Van PraetJ. T.MarzoratiM.DrennanM. B.ElewautD. (2016). How the microbiota shapes rheumatic diseases. Nat. Rev. Rheumatol. 12, 398–411. 10.1038/nrrheum.2016.85
51
van VollenhovenR. F. (2009). Sex differences in rheumatoid arthritis: more than meets the eye. BMC Med. 7:12. 10.1186/1741-7015-7-12
52
WuJ.PetersB. A.DominianniC.ZhangY.PeiZ.YangL.et al. (2016). Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 10, 2435–2446. 10.1038/ismej.2016.37
53
YangM. J.HouY. L.YangX. L.WangC. X.ZhiL. X.YouC. G. (2018). Development and application of a PCR-HRM molecular diagnostic method of SNPs linked with TNF inhibitor efficacy. Diagnosis6, 277–286. 10.1515/dx-2018-0062
54
ZvanychR.LukendaN.KimJ. J.LiX.PetrofE. O.KhanW. I.et al. (2014). Small molecule immunomodulins from cultures of the human microbiome member Lactobacillus plantarum. J. Antibiot. 67, 85–88. 10.1038/ja.2013.126
Summary
Keywords
rheumatoid arthritis, gut microbiome, 16S rRNA gene sequencing, inflammation, biomarker
Citation
Sun Y, Chen Q, Lin P, Xu R, He D, Ji W, Bian Y, Shen Y, Li Q, Liu C, Dong K, Tang Y-W, Pei Z, Yang L, Lu H, Guo X and Xiao L (2019) Characteristics of Gut Microbiota in Patients With Rheumatoid Arthritis in Shanghai, China. Front. Cell. Infect. Microbiol. 9:369. doi: 10.3389/fcimb.2019.00369
Received
06 May 2019
Accepted
09 October 2019
Published
23 October 2019
Volume
9 - 2019
Edited by
Saleh A. Naser, University of Central Florida, United States
Reviewed by
Robert Cody Sharp, University of Florida Health, United States; Mohtashem Samsam, University of Central Florida College of Medicine, United States; Susu M. Zughaier, Qatar University, Qatar
Updates
Copyright
© 2019 Sun, Chen, Lin, Xu, He, Ji, Bian, Shen, Li, Liu, Dong, Tang, Pei, Yang, Lu, Guo and Xiao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiaokui Guo xkguo@shsmu.edu.cnLianbo Xiao 13701888178@163.com
This article was submitted to Microbiome in Health and Disease, a section of the journal Frontiers in Cellular and Infection Microbiology
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.