Data-driven multiple-level analysis of gut-microbiome-immune-joint interactions in rheumatoid arthritis
BMC Genomics volume 20, Article number: 124 (2019)
Rheumatoid arthritis (RA) is the most common autoimmune disease and affects about 1% of the population. The cause of RA remains largely unknown and could result from a complex interaction between genes and environment factors. Recent studies suggested that gut microbiota and their collective metabolic outputs exert profound effects on the host immune system and are implicated in RA. However, which and how gut microbial metabolites interact with host genetics in contributing to RA pathogenesis remains unknown. In this study, we present a data-driven study to understand how gut microbial metabolites contribute to RA at the genetic, functional and phenotypic levels.
We used publicly available disease genetics, chemical genetics, human metabolome, genetic signaling pathways, mouse genome-wide mutation phenotypes, and mouse phenotype ontology data. We identified RA-associated microbial metabolites and prioritized them based on their genetic, functional and phenotypic relevance to RA. We evaluated the prioritization methods using short-chain fatty acids (SCFAs), which were previously shown to be involved in RA etiology. We validate the algorithms by showing that SCFAs are highly associated with RA at genetic, functional and phenotypic levels: SCFAs ranked at top 3.52% based on shared genes with RA, top 5.69% based on shared genetic pathways, and top 16.94% based on shared phenotypes. Based on the genetic-level analysis, human gut microbial metabolites directly interact with many RA-associated genes (as many as 18.1% of all 166 RA genes). Based on the functional-level analysis, human gut microbial metabolites participate in many RA-associated genetic pathways (as many as 71.4% of 311 genetic pathways significantly enriched for RA), including immune system pathways. Based on the phenotypic-level analysis, gut microbial metabolites affect many RA-related phenotypes (as many as 51.3% of 978 phenotypes significantly enriched for RA), including many immune system phenotypes.
Our study demonstrates strong gut-microbiome-immune-joint interactions in RA, which converged on both genetic, functional and phenotypic levels.
Rheumatoid arthritis (RA) is one of the most common autoimmune diseases and affects over 1.3 million Americans and 1% of the worldwide population ((https://www.rheumatoidarthritis.org/ra/facts-and-statistics/), ). RA is complex, with genetic, epigenetic, and environmental factors contributing to disease susceptibility and progression . While significant progress has been made in understanding genetic, molecular, and cellular aspects of RA, relatively little is known about which environmental factors are important in RA susceptibility and how they interact with host genetics in the development of RA . Human and mouse model microbiome studies have shown that gut dysbiosis, an imbalance in the intestinal microbiota , is associated with RA [5,6,7,8,9,10,11]. Studies in mouse models have shown a requirement of gut microbiota for arthritis development [12,13,14].
Human gut microbiota contribute to human diseases and health via the cumulative effects of microbial metabolites [15,16,17]. It has become increasingly clear that the prodigious metabolite activities of gut microbiota strongly influence RA susceptibility and progression [6, 7, 18,19,20,21]. Short chain fatty acids (SCFA) are the primary end-products of fermentation of non-digestible dietary fiber by the gut microbiota. SCFAs have emerged as major mediators in linking nutrition, gut microbiota, and human health [22, 23]. Recent studies have shown that SCFAs play important roles in the suppression of inflammation in RA . Mice deficient for SFCA receptor showed exacerbated inflammation in modes of RA . Butyrate, one of the most abundant SCFAs, acts as an endogenous histone deacetylase (HDAC) inhibitor and has been shown to decrease inflammation in animal models of RA and other inflammatory diseases .
Although the link between microbial metabolism and RA has been recognized, the mechanisms underlying how microbial metabolites interact with human genetics in promoting or protecting against RA remain largely unknown. We previously demonstrated that data-driven computational approaches have potential in uncovering mechanistic links between microbial metabolites and human diseases (colorectal cancer and Alzheimer’s disease) [24,25,26]. Specific for RA, we previously developed a mechanism-based prediction system, mMetabolitePredict, for human metabolome biomarker discovery and applied it to identify and prioritize metabolomic biomarkers for RA . We found that among 259,170 prioritized chemicals/metabolites in human body, the subset of metabolites originated from human gut microbiota ranked highly . This finding motivated our current study (funded by Pfizer ASPIRE Rheumatology and Dermatology Research Award) to perform data-driven systematic analysis of which and how human gut microbial metabolites are involved in the immune-joint axis of the RA etiology at the genetic, functional and phenotypic levels. We evaluated the algorithms using SCFAs, which are known to have a role in the suppression of inflammation in RA [19, 20]. We evaluated whether SCFAs were ranked highly based on their genetic, functional and phenotypic relevance to RA. To the best of our knowledge, our study represents the first computational approach to comprehensively characterize the complex gut-microbiome-immune-joint interactions in RA. The unique informatics contribution is that we innovatively leveraged large amounts of publicly available data collected for other purpose and developed data-driven computational methods to understand gut-microbiome-gene-disease interactions. Our approaches are highly flexible and can be applied to any other diseases. The biomedical contribution of our study is that the identification of gut microbial metabolites and the understanding of their role in RA has potential in providing new insights into the basic mechanisms of disease etiology and enable new possibilities for disease diagnosis, prevention, and treatment.
Genetic connections: Microbial metabolites may be genetically involved in RA and interact with many RA-associated genes
For evaluation, we show that the genetics-based ranking algorithm ranked all three SCFAs consistently highly across three complementary disease genetics data resources (Table 1). SCFAs on average ranked in the top 3.52% among 127 gut microbial metabolites, with acetic acid ranked at top 1, butyric acid at top 2 and propionic acid at top 9. Our study shows that butyric acid regulates many RA-related genes, including IL10, IL2, IL6, and STAT4 (Table 2), suggesting the potential roles of SCFAs for their anti-inflammatory effects in protecting joint in RA.
Our studies show that 61 out of the 127 gut microbial metabolites directly interact with RA-associated genes. The top 10 microbial metabolites (ranked based on the number of shared genes with RA) and their shared RA gene are shown in Table 2. For example, acetic acid ranked at top 1 and shared 30 genes with RA (18.1% of all 166 RA genes). Many of the shared genes are immune-related, strongly suggesting the gut-microbiome-immune-joint interactions in RA.
Functional connections: Microbial metabolites may be functionally involved in RA and participate in many RA-associated genetic pathways
We identified genetic pathways significantly associated with RA and for each microbial metabolite. We then ranked metabolites based on the numbers of shared genetic pathways with RA. All three SCFAs ranked highly based on their pathway overlaps with RA (Table 3). A total of 311 pathways were significantly enriched for RA, among which butyric acid shared 222 pathways (71.4%), acetic acid shared 126 pathways (40.5%), and propionic acid shared 152 pathways (48.9%). Among 127 microbial metabolites, 116 metabolites shared at least one genetic pathway with RA. The top 20 ranked metabolites are shown in Table 4. The fact that the majority of RA-associated genetic pathways were regulated by SCFAs and other gut microbial metabolites indicates that human gut microbial metabolism is functionally involved in RA etiology.
We then ranked the shared pathways between RA and each metabolite by the balanced measure of their enrichment folds. For example, the pathway “IL27-mediated signaling events” was 10.31-fold enriched for RA and 6.48-fold enriched for butyric acid. The F1 combined enrichment fold of this pathways was 7.96. The top 10 genetic pathways regulated by both RA and butyric acid are shown in Table 5. As shown in the table, the majority of the top shared pathways are related to immune functions, strongly suggesting the gut-microbiome-immune interactions in RA at functional-level. We performed the same analysis for the other two SCFAs (acetic acid and propionic acid) and for methane (the top one ranked and non-SCFA metabolite). Table 6 shows the top ten ranked pathways shared between RA and acetic acid, propionic acid and methane. The results show that majority of the top shared pathways between RA and metabolites (SCFAs and non-SCFA methane) are related to immune functions, though the specific pathways for each metabolite vary. In summary, microbial metabolites may be involved in RA pathology through different immune pathways.
Phenotypic connections: Microbial metabolites may affect RA at the phenotypic level and affect many RA-related phenotypes
We examined the phenotypic connections between gut microbial metabolites and RA. As the evaluation, we showed that SCFAs were significantly associated with RA at the phenotypic level (Table 7). The top 20 metabolites ranked based on phenotypic overlaps with RA are shown in Table 8. For example, a total of 978 phenotypes were significantly enriched for RA-associated genes (166 genes from combined resources), among which butyric acid shared 502 phenotypes (51.3%) and propionic acid shared 335 phenotypes (34.3%) with RA. These results indicate that SCFAs and other microbial metabolites are phenotypically involved in RA.
Case study: Butyric acid is phenotypically involved in RA
Butyric acid is the most abundant metabolites of gut microbiota in the fermentation of dietary fiber. Our above analysis showed that butyric acid was highly associated with RA at both genetic, functional and phenotypic levels. We then examined how butyric acid was phenotypically involved in RA. Butyric acid shared 502 phenotypes with RA. We classified these shared phenotypes based on the Mouse Phenotype Ontology (MPO) classification schemes (3rd- and 4th-level classifications). These 502 phenotypes were classified into 52 3rd-level classes and 164 4th-level classes.
The top 10 3rd-level classes are shown in Fig. 1. The 3rd-level phenotype class “abnormal immune system physiology” ranked at top one. Among a total of 502 phenotypes shared between RA and butyric acid, 148 (23%) phenotypes belonged to this class. The top ten 4th-level classes of the shared phenotypes are shown in Fig. 2, among which seven phenotypes were directly related to immune functions, including “abnormal immune serum protein physiology”, “abnormal inflammatory response”, “abnormal cell-mediated immunity” and “abnormal adaptive immunity”.
Since the gene-phenotype annotations in MGI are largely mutational, our phenotypic-level analysis suggests potential causal phenotypic effects of butyric acid and other gut microbial metabolites on RA, though effects on various immune functions.
In this study, we performed data-driven analysis of the gut-microbiome-immune-joint interactions in RA. We showed that gut microbial metabolites were strongly involved in RA at both genetic, functional and phenotypic levels. Specifically, our study shows that gut microbial metabolites interact with many RA-associated genes, participate in RA-related immune pathways and affect RA-associated immunological phenotypes. As compared to our previous studies for colorectal cancer and Alzheimer’s diseases, the identified microbial metabolites as well as the genetic pathways involved are different in RA from ones in colorectal cancer and Alzheimer’s disease [24,25,26]. For example, SCFAs are highly related to RA through immune pathways. On the other hand, trimethylamine-n-oxide is highly related to colorectal cancer through many cancer-related pathways [24, 26]. To the best of our knowledge, our study represents the first computational approach to comprehensively characterize the complex gut-microbiome-immune- joint interactions in RA. While our study is pure ‘in silico’, it generated large amounts of data/knowledge/hypotheses that can facilitate other biomedical researchers to conduct hypothesis-driven functional studies of gut microbial metabolisms in RA. A few limitations inherent in the publicly available datasets warrant further discussion.
First, although our analysis suggests strong functional connections between gut microbial metabolites and RA, especially in our phenotypic analysis that used the causal/mutational gene-phenotype associations, our findings are largely associational. In order to translate the findings for RA diagnosis, prevention, and treatment, it is necessary to establish cause-effect relationships of the identified metabolite-gene-pathway-phenotype-RA associations and identify specific gut bacteria that produce the metabolites.
Second, HMDB is currently the most comprehensive human metabolome database, containing a total of 83,479 small molecule metabolites found in the human body. STITCH is currently the most comprehensive chemical genetics databases, containing genetic associations for 500,000 small molecules. However, HMDB contains only 172 metabolites originated from gut microbiota, among which 127 metabolites have associated genes in STITCH. The field of human microbiome research is fast moving, with an increasing number of microbial metabolites being identified and published in literature. We are currently developing text mining and natural language processing techniques to extract human gut microbial metabolites from published biomedical literature.
Third, our data-driven multi-level analysis is not specific to RA and can be applied to any other diseases. The only change to the algorithm is to replace RA-associated genes to genes associated with another disease. However, the lack of evaluation data (known microbial metabolites associated with diseases) in structured format prevented us from systematically evaluating how the algorithm perform in other diseases. Increasing number of published biomedical research articles have reported the associations among microbial metabolites, gut bacteria, and diseases. However, the knowledge is still buried in free-text documents with limited machine understandability. In order to systematic evaluate our algorithm in other diseases, we need either manually curate or develop natural language processing techniques to automatically extract disease-microbial metabolite associations from biomedical literature. We are actively pursuing the latter approach.
Fourth, studies have shown the importance of diet and associated changes in the gut microbiota in human diseases. Intestinal SFCAs are produced by gut microbiota digesting high fiber diet and are involved in human metabolism and health [22, 23]. Trimethylamine n-oxide (TMAO) is formed by gut microbiota in digesting red meat and high fat diet. Both human studies have shown that TMAO is mechanistically involved in cardiovascular diseases , renal disease  and colorectal cancer . However, the exact relationship among specific diet, bacteria and diseases remain largely unknown. In one of our recent studies, we developed network-based systems approach to investigate food-metabolite interactions in Alzheimer’s disease . The food metabolites can be produced by either human host or gut bacteria, however we currently lack the knowledge to differentiate these two. The identification and understanding how diet and food are associated with diseases by impacting gut microbiota will have great potential in treating and preventing human diseases, including RA.
The computational framework that we developed has built-in flexibility and capability for us to continuously incorporate new data as it becomes available in our future studies. We believe that our view of gut-microbiome-RA interactions will evolve as more data becomes available.
The cause of RA remains largely unknown and could result from a complex interaction between genes and environment factors. Recent studies suggested that gut microbiota and their collective metabolic outputs exert profound effects on the host immune system and are implicated in RA. However, which and how gut microbial metabolites interact with host genetics in contributing to RA pathogenesis remains unknown. In this study, we present a data-driven study to understand how gut microbial metabolites contribute to RA at the genetic, functional and phenotypic levels. Our in-silico data-driven study suggests strong gut-microbiome-immune-joint interactions in RA, which converge on both genetic, functional and phenotypic levels.
We used publicly available human disease genetics, human chemical genetics, human metabolome, genetic signaling pathways, mouse genome-wide mutation phenotypes, and mouse phenotype ontology to characterize the genetic, functional, and phenotypic connections between human gut microbial metabolites and RA.
RA genetics data
We used three data resources to obtain RA-associated genes: (1) we obtained 155 RA-associated genes form the Catalog of Published Genome-Wide Association Studies (GWAS catalog) (data accessed in June, 2017). The GWAS catalog is an exhaustive source of disease/trait-gene associations from published GWAS data and currently contains 34,790 disease/trait-gene pairs for 1655 common complex diseases/traits , 2) we obtained 16 RA-associated genes from the Online Mendelian Inheritance in Man database (OMIM) (data accessed in July, 2017). OMIM is the most comprehensive source of disease genetics for Mendelian disorders and currently includes 10,125 disease-gene pairs for 10,674 diseases/phenotypes ; and (3) we obtained 10 RA-associated genes from ClinVar (data accessed in July, 2017). ClinVar is a publicly available resource of reports of the relationships among human variations and phenotypes and currently contains 9873 disease-gene associations for 5240 diseases/phenotypes . We used these three complementary disease genetics resources to demonstrate the robustness of our findings.
The human metabolome database (HMDB)
HMDB contains detailed information about small molecule metabolites found in the human body . Currently, HMDB contains 83,479 metabolites. In this study, we focused on the 172 metabolites originated in human gut microbiota (data accessed in July, 2017).
Metabolite genetics data
We used STITCH (Search Tool for Interactions of Chemicals) database to obtain genes associated with gut microbial metabolites obtained from HMDB. STITCH contains data on the interactions between 500,000 small molecules and 9.6 million proteins from 2031 organisms . In this study, we used chemical-gene associations found in human body, which include 15,473,939 chemical-gene pairs for 473,602 chemicals, and 18,701 human genes (data accessed in July, 2017). Among the 172 microbial metabolites from HMDB, 127 were mapped.
to chemical names in STITCH. Genes associated with the mapped microbial metabolites were then obtained from STITCH. For example, we mapped butyric acid (in HMDB) to butyrate (in STITCH) and obtained a total of 793 butyrate-associated genes from STITCH.
Molecular pathway data
We used rich pathway information from the Molecular Signatures Database (MSigDB) to investigate how microbial metabolites were functionally related to RA. MSigDB is currently the most comprehensive resource of 17,779 annotated pathways and gene sets  (data accessed in July, 2017). For each microbial metabolite, we identified molecular pathways that were significantly enriched for both RA and the metabolite.
Genome-wide mutational phenotypes in mouse models and mouse phenotype ontology
The Mouse Genome Database (MGD) made available large amounts of phenotypic descriptions of systematic genetic knockouts in mouse models . Such large-scale systematic genetic knockouts are impossible to do in human. These strong causal gene-phenotype associations (318,709 gene-phenotype association annotations for 60,474 targeted mutant alleles and 12,104.
phenotypes) have been useful for screening functional effects of chemicals on disease phenotypes. We have recently developed computational algorithms to perform virtual phenotypic screening to prioritize drugs for diseases by matching mouse mutational phenotype profiles between drugs and diseases [39,40,41,42]. We validated our virtual screening drug candidates in experimental models of ovarian cancer [39, 40]. In this study, we used the same strategy to assess the phenotypic effects of gut microbial metabolites on RA-related phenotypes. For example, the microbial metabolite butyrate is associated with the gene IL17A and the knockout of IL17A in mouse models is associated with the phenotype “rheumatoid arthritis”. We then used the classification scheme of the Mouse Phenotype Ontology (MPO) at MGD to group identified phenotypes (e.g., “rheumatoid arthritis”, “abnormal cytokine level”) into classes (e.g., “abnormal immune system physiology”).
Analyze genetic connections between gut microbial metabolites and RA and prioritize metabolites based on their shared genes with RA
We hypothesize that if a metabolite interacts with many RA-associated genes, then this metabolite may be highly involved in RA. We prioritized microbial metabolites based on number of their shared genes with RA. We obtained RA-associated genes from three complementary disease genetics data resources (OMIM, the GWAS Catalog, and Clinvar): 15 RA-associated genes from OMIM, 155 genes from the GWAS Catalog, 10 genes from ClinVar, and 166 genes from these three resources combined. We obtained microbial metabolite-associated genes from the STITCH database. Metabolites were ranked based on the numbers of their shared genes with RA.
Animal studies showed that SCFAs had a role in the suppression of inflammation in RA [19, 20]. For evaluation, we examined whether SCFAs were ranked highly based on their genetic overlaps with RA. Both mean and median rankings of SCFAs among all metabolites were calculated. Significance was calculated by comparing mean rankings to random expectation, which is 50%. Rankings based on three disease genetics data resources were compared to demonstrate the robustness of our analysis and findings.
Analyze functional connections between gut microbial metabolites and RA and prioritize metabolites based on their shared genetic pathways with RA
We obtained RA-associated genes from the three disease genetics databases separately. Pathways associated with each gene were obtained from MSigDB. For each pathway, we assessed its probability of being associated with a given set of RA-associated genes (e.g., 15 RA genes from the OMIM database) as compared to its probability of being associated with the same number of randomly selected genes. The random process was repeated 1000 times and a t-test was used to assess the statistical significance. For instance, a total of 108 pathways were significantly associated with the 16 RA-associated genes from OMIM. The pathway “Cytokines and Inflammatory Response” had a significant 61-fold enrichment as compared to the random expectation. Similarly, we identified significantly enriched genetic pathways for each of the 127 microbial metabolites (with metabolite-associated genes obtained from the STITCH database). For example, a total of 748 pathways were significantly enriched for butyric acid. The pathway “Cytokines and Inflammatory Response” shows a significant 5-fold enrichment for butyric acid.
A shared genetic pathway may be associated with RA and a metabolite at different significance level. We developed a prioritization measure to identify shared pathways that rank highly for both RA and a metabolite. The intuition is that a shared pathway between RA and a metabolite ranks highly if and only if it ranks highly for both RA and the metabolite. The ranking of a shared pathway between RA and a metabolite s defined as: rank = 2*(ranking_ra * ranking_m)/(ranking_ra + 2 ranking_m), where ranking_ra is the enrichment fold of a pathway for RA; and ranking_m is the enrichment fold of the same pathway for the metabolite. For example, the pathway “cytokines and inflammatory response” showed a 61-fold enrichment for RA and a 5-fold enrichment for butyric acid. The combined ranking score of this shared pathway for both RA and butyric acid was 9.24. After identifying shared pathways, metabolites were then prioritized based on the numbers of shared pathways with RA.
The prioritization algorithm was evaluated using three known RA-associated SCFAs. Mean ranking, median rankings and the significance were calculated. Rankings based on three different disease genetics data resources were compared to demonstrate the robustness of the finding.
Analyze phenotypic connections between gut microbial metabolites and RA and prioritize metabolites based on their shared phenotypes with RA
We obtained RA-associated genes to their corresponding mouse gene homologs (e.g., IL17A = > Ctla) using human-mouse homolog mapping data from MGD . The mapped mouse genes were then linked to their corresponding mutational phenotypes in mouse models (e.g., IL17A = > rheumatoid arthritis, TNF = > abnormal inflammatory response) using gene-phenotype association annotations from MGD. For each mapped phenotype, we assessed its probability of being associated with RA-associated genes as compared to its probability associated with the same number of randomly selected genes. The random process was repeated 1000 times and a t-test was used to assess the statistical significance. As an example, the phenotype “abnormal T-helper 1 physiology” showed a significant 36-fold enrichment for RA as compared to random expectation. Similarly, we identified significantly enriched phenotypes for each gut microbial metabolite. For example, the phenotype “abnormal T-helper 1 physiology” shows a significant 1.7-fold enrichment for butyric acid. Phenotypes shared between RA and each metabolite were then prioritized as described above for prioritizing shared genetic pathways. After identifying shared phenotypes, metabolites were then prioritized based on the numbers of shared phenotypes with RA.
The prioritization algorithm was evaluated using three known RA-associated SCFAs. Mean ranking, median rankings and the significance were calculated. Rankings based on three different disease genetics data resources were compared to demonstrate the robustness of the finding.
Short-chain fatty acids
- GWAS catalog:
The Catalog of Published Genome-Wide Association Studies
The Online Mendelian Inheritance in Man database
The Human Metabolome Database
Search Tool for Interactions of Chemicals
The Molecular Signatures Database
The Mouse Genome Database
Firestein GS. Evolving concepts of rheumatoid arthritis. Nature. 2003;423(6937):356.
McInnes IB, Schett G. The pathogenesis of rheumatoid arthritis. N Engl J Med. 2011;365(23):2205–19.
Karlson EW, Deane K. Environmental and gene-environment interactions and risk of rheumatoid arthritis. Rheum Dis Clin N Am. 2012;38(2):405–26.
DeGruttola AK, Low D, Mizoguchi A, Mizoguchi E. Current understanding of dysbiosis in disease in human and animal models. Inflamm Bowel Dis. 2016;22(5):1137–50.
Taneja V. Arthritis susceptibility and the gut microbiome. FEBS Lett. 2014;588(22):4244–9.
Scher JU, Abramson SB. The microbiome and rheumatoid arthritis. Nat Rev Rheumatol. 2011;7(10):569.
Clemente JC, Manasson J, Scher JU. The role of the gut microbiome in systemic inflammatory disease. Bmj. 2018;360:j5145.
Scher JU, Sczesnak A, Longman RS, Segata N, Ubeda C, Bielski C, et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. elife. 2013;2.
Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, et al. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat Med. 2015;21(8):895.
Forbes JD, Chen CY, Knox NC, Marrie RA, El-Gabalawy H, de Kievit T, et al. A comparative study of the gut microbiota in immune-mediated inflammatory diseases—does a common dysbiosis exist? Microbiome. 2018;6(1):221.
Picchianti-Diamanti A, Panebianco C, Salemi S, Sorgi M, Di Rosa R, Tropea A, et al. Analysis of gut microbiota in rheumatoid arthritis patients: disease-related Dysbiosis and modifications induced by Etanercept. Int J Mol Sci. 2018;19(10):2938.
Block KE, Zheng Z, Dent AL, Kee BL, Huang H. Gut microbiota regulates K/BxN autoimmune arthritis through follicular helper T but not Th17 cells. J Immunol. 2016;196(4):1550–7.
Abdollahi-Roodsaz S, Joosten LA, Koenders MI, Devesa I, Roelofs MF, Radstake TR, et al. Stimulation of TLR2 and TLR4 differentially skews the balance of T cells in a mouse model of arthritis. J Clin Invest. 2008;118(1):205.
Wu HJ, Ivanov II, Darce J, Hattori K, Shima T, Umesaki Y, et al. Gut-residing segmented filamentous bacteria drive autoimmune arthritis via T helper 17 cells. Immunity. 2010;32(6):815–27.
Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, Pettersson S. Host-gut microbiota metabolic interactions. Science. 2012;336(6086):1262–7.
Tremaroli V, Backhed F. Functional interactions between the gut microbiota and host metabolism. Nature. 2012;489(7415):242–9.
Ray K. Gut microbiota: microbial metabolites as mimickers of human molecules. Nat Rev Gastroenterol Hepatol. 2017;14:630–1.
Maslowski KM, Mackay CR. Diet, gut microbiota and immune responses. Nat Immunol. 2011;12(1):5–9.
Maslowski KM, Vieira AT, Ng A, Kranich J, Sierro F, Yu D, Schilter HC, Rolph MS, Mackay F, Artis D, et al. Regulation of inflammatory responses by gut microbiota and chemoattractant receptor GPR43. Nature. 2009;461:1282–6.
Wang L, de Zoeten EF, Greene MI, Hancock WW. Immunomodulatory effects of deacetylase inhibitors: therapeutic targeting of FOXP3+ regulatory T cells. Nat. Rev Drug Discov. 2009;8:969–81.
Thorburn AN, Macia L, Mackay CR. Diet, metabolites, and “western-lifestyle” inflammatory diseases. Immunity. 2014;40(6):833–42.
Ríos-Covián D, Ruas-Madiedo P, Margolles A, Gueimonde M, De los Reyes-Gavilán CG, Salazar N. Intestinal short chain fatty acids and their link with diet and human health. Frontiers in microbiology. 2016;7:185.
Morrison DJ, Preston T. Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut Microbes. 2016;7(3):189–200.
Xu R, Wang Q, Li L. Genome-wide systems analysis reveals strong link between colorectal cancer and trimethylamine N-oxide (TMAO), a gut microbial metabolite of dietary meat and fat. BMC Genomics. 2015;16(Suppl 7):S4.
Xu R, Wang Q. Towards understanding brain-gut-microbiome connections in Alzheimer’s disease. BMC Syst Biol. 2016;10(63). https://doi.org/10.1186/s12918-016-0307-y.
Wang Q, Li L, Xu R. A systems biology approach to predict and characterize human gut microbial metabolites in colorectal cancer. Sci Rep. 2018;8(1):6225.
Wang Q, Xu R. MetabolitePredict: a de novo human metabolomics prediction system and its applications in rheumatoid arthritis. J Biomed Inform. 2017;71:222–8.
Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, DuGar B, et al. Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. 2011;472(7341):57.
Tang WW, Wang Z, Kennedy DJ, Wu Y, Buffa JA, Agatisa-Boyle B, et al. Gut microbiota-dependent trimethylamine N-oxide (TMAO) pathway contributes to both development of renal insufficiency and mortality risk in chronic kidney DiseaseNovelty and significance. Circ Res. 2015;116(3):448–55.
Bae, S, Ulrich, C. M, Neuhouser, M. L, Malysheva, O, Bailey, L. B, Xiao, L, ... & Miller, J. W. (2014). Plasma choline metabolites and colorectal cancer risk in the Women's Health Initiative observational study. Cancer research 2014 Dec 15;74(24):7442–7452.
Chen Y, Xu R. Context-sensitive network analysis identifies food metabolites associated with Alzheimer’s disease: an exploratory study. BMC Genomics. 2018. https://doi.org/10.1186/s12920-018-0459-2MGNM-D-10-00077.
MacArthur J, Bowler E, Cerezo M, Gil L, Hall P, Hastings E, et al. The new NHGRI-EBI catalog of published genome-wide association studies (GWAS catalog). Nucleic Acids Res. 2016;45(D1):D896–901.
Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM. Org: online Mendelian inheritance in man (OMIM), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2014;43(D1):D789–98.
Landrum MJ, Lee JM, Benson M, Brown G, Chao C, Chitipiralla S, et al. ClinVar: public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2015;44(D1):D862–8.
Wishart DS, Jewison T, Guo AC, Wilson M, Knox C, Liu Y. & Bouatra, S. (2013). HMDB 3.0-the human metabolome database in. Nucleic Acids Res. 2013;41(D1):D801–7.
Szklarczyk D, Santos A, Von Mering C, Jensen LJ, Bork P, Kuhn M. STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 2015;44(D1):D380–4.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739–40.
Eppig JT, Blake JA, Bult CJ, Kadin JA, Richardson JE, Mouse Genome Database Group. The mouse genome database (MGD): facilitating mouse as a model for human biology and disease. Nucleic Acids Res. 2014;43(D1):D726–36.
Nagaraj AB, Wang QQ, Joseph P, Zheng C, Chen Y, Kovalenko O, et al. Using a novel computational drug-repositioning approach (DrugPredict) to rapidly identify potent drug candidates for cancer treatment. Nature Oncogene. 2018;37(3):403.
Nagaraj AB, Joseph P, Kovalenko O, Wang Q, Xu R, DiFeo A. Evaluating class III antiarrhythmic agents as novel MYC targeting drugs in ovarian cancer. Gynecol Oncol. 2018;151(3):525–32.
Chen Y, Gao Z, Wang B, Xu R. Towards precision medicine-based therapies for glioblastoma: interrogating human disease genomics and mouse phenotypes. BMC Genomics. 2016;17(7):516.
Chen, Y., Cai, X., & Xu, R. (2015). Combining Human Disease Genetics and Mouse Model Phenotypes towards Drug Repositioning for Parkinson’s disease. In AMIA Annual Symposium Proceedings (Vol. 2015, p. 1851). American Medical Informatics Association.
Availability of data and material
Q.W. and R.X. have jointly conceived the study, designed the experiment, performed the experiment and wrote the manuscript. All authors have participated in study discussion and manuscript preparation. All of the authors have read and approved the final manuscript.
This work was supported by the Pfizer 2015 ASPIRE Rheumatology and Dermatology Research Award (WI206753, Xu). Additional support for Xu and wang are from Eunice Kennedy Shriver National Institute of Child Health & Human Development of the National Institutes of Health under the NIH Director’s New Innovator Award number DP2HD084068 (Xu), NIH National Institute of Aging (1 R01 AG057557–01, Xu), NIH National Institute of Aging (1 R01 AG061388–01, Xu), NIH National Institute of Aging (1 R56 AG062272–01, Xu), American Cancer Society Research Scholar Grant (RSG-16-049-01 - MPC, Xu), NIH Clinical and Translational Science Collaborative of Cleveland (1UL1TR002548–01, Konstan). Publication of this article was sponsored by the Pfizer 2015 ASPIRE Rheumatology and Dermatology Research Award (WI206753, Xu). The funding bodies had no role in in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The 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
Wang, Q., Xu, R. Data-driven multiple-level analysis of gut-microbiome-immune-joint interactions in rheumatoid arthritis. BMC Genomics 20, 124 (2019). https://doi.org/10.1186/s12864-019-5510-y