Article Text


Gut Microbiota in Multiple Sclerosis
  1. Brandi L. Cantarel, PhD*†,
  2. Emmanuelle Waubant, MD, PhD,
  3. Christel Chehoud§,
  4. Justin Kuczynski, PhD§,
  5. Todd Z. DeSantis, MSc§,
  6. Janet Warrington, PhD§,
  7. Arun Venkatesan, MD, PhD,
  8. Claire M. Fraser, PhD†¶,
  9. Ellen M. Mowry, MD, MCR
  1. From the *Baylor Institute for Immunology Research, Baylor Health, Dallas, TX; †Institute for Genome Sciences, University of Maryland School of Medicine, Baltimore, MD; ‡Department of Neurology, UCSF School of Medicine, San Francisco; and §Second Genome Inc, South San Francisco, CA; and ∥Department of Neurology, Johns Hopkins University School of Medicine; and ¶Department of Medicine, University of Maryland School of Medicine, Baltimore, MD.
  1. Received July 28, 2014, and in revised form December 28, 2014.
  2. Accepted for publication February 6, 2015.
  3. The study was funded by NIH K23 NS067055 (to E.M.M.) and an Irene Perstein Award (to E.M.M.).
  4. E.M.M. is currently receiving free glatiramer acetate in support of the conduct of a clinical trial, and some of the authors (C.C., J.K., T.Z.D., and J.W.) are employees of Second Genome, Inc.
  5. Reprints: Ellen M. Mowry, MD, MCR, The Johns Hopkins Hospital, 600 N Wolfe St, Baltimore, MD 21287. E-mail: emowry1{at}
  6. Supplemental digital content is available for this article. Direct URL citation appears in the printed text and is provided in the HTML and PDF versions of this article on the journal’s Web site (

Possible Influence of Immunomodulators


Objectives Differences in gut bacteria have been described in several autoimmune disorders. In this exploratory pilot study, we compared gut bacteria in patients with multiple sclerosis and healthy controls and evaluated the influence of glatiramer acetate and vitamin D treatment on the microbiota.

Methods Subjects were otherwise healthy white women with or without relapsing-remitting multiple sclerosis who were vitamin D insufficient. Patients with multiple sclerosis were untreated or were receiving glatiramer acetate. Subjects collected stool at baseline and after 90 days of vitamin D3 (5000 IU/d) supplementation. The abundance of operational taxonomic units was evaluated by hybridization of 16S rRNA to a DNA microarray.

Results While there was overlap of gut bacterial communities, the abundance of some operational taxonomic units, including Faecalibacterium, was lower in patients with multiple sclerosis. Glatiramer acetate–treated patients with multiple sclerosis showed differences in community composition compared with untreated subjects, including Bacteroidaceae, Faecalibacterium, Ruminococcus, Lactobacillaceae, Clostridium, and other Clostridiales. Compared with the other groups, untreated patients with multiple sclerosis had an increase in the Akkermansia, Faecalibacterium, and Coprococcus genera after vitamin D supplementation.

Conclusions While overall bacterial communities were similar, specific operational taxonomic units differed between healthy controls and patients with multiple sclerosis. Glatiramer acetate and vitamin D supplementation were associated with differences or changes in the microbiota. This study was exploratory, and larger studies are needed to confirm these preliminary results.

Key Words
  • multiple sclerosis
  • gut microbiome
  • vitamin D
  • glatiramer acetate
  • autoimmunity

Statistics from

Key Words

Multiple sclerosis (MS), an autoimmune neurologic disease caused by a combination of genetic and environmental exposures, has increased in incidence in the past several decades,1suggesting changing environmental risk factors.

The gut microbiota is symbiotic with the human and is crucial to normal immune function.2Differences in gut bacteria exist in those with autoimmune disorders (eg, inflammatory bowel disease)2compared with healthy individuals, and antibiotic and probiotic regimens are being assessed as potential therapies. Altering gut bacteria influences the development and severity of experimental autoimmune encephalomyelitis, a mouse model of MS.3A more recent study has shown that unaltered commensal bacteria can trigger, after exposure to myelin oligodendrocyte glycoprotein, a spontaneous form of experimental autoimmune encephalomyelitis.4No study has assessed gut bacteria in patients with MS. We sought to evaluate the gut microbiota in healthy controls (HCs) and those with MS and to explore differences associated with MS treatment and supplementation with vitamin D, both of which are used commonly in MS clinical practice.


We capitalized on a vitamin D3 supplementation study (NCT01667796) to conduct this substudy. All subjects provided written informed consent, and the study was approved by the UCSF Institutional Review Board and thus is in compliance with the Declaration of Helsinki. Because the parent study was a small pilot study, strict inclusion and exclusion criteria were used to reduce heterogeneity that might otherwise have impacted the results. Subjects were white, non-Hispanic women aged 18 to 60 years with screening body mass index of 18 to 30 kg/m2and screening 25-hydroxyvitamin D level of 30 ng/mL or less. Subjects had to be willing to stop taking additional multivitamins or cod liver oil and to avoid tanning beds for the duration of the study. Subjects could not be pregnant or nursing and had to be willing to avoid pregnancy during the study. Those with known gastrointestinal disorders were excluded from the overall study because such diseases could impact the absorption of vitamin D; for similar reasons, subjects could not be on a no-fat diet. Those with known renal or liver disease, history of nephrolithiasis, history of hypercalcemia or hypercalciuria, screening serum calcium greater than 10 mg/dL, anemia (hemoglobin <11 g/dL) hyperthyroidism, infection with Mycobacterium species, sarcoidosis, or other serious chronic illness (including cancer, cardiac, HIV) were excluded. Furthermore, those who were taking thiazide diuretics, digoxin, diltiazem, verapamil, cimetidine, heparin, low-molecular-weight heparin, or medication associated with malabsorption were excluded. Finally, subjects were excluded if in the month prior to screening they had smoked cigarettes, used illicit drugs, or had taken nontopical steroids. Patients with MS had relapsing-remitting disease by McDonald criteria with Expanded Disability Status Scale scores of 3.0 or less and without major heat sensitivity, criteria that were chosen to minimize differences in sun exposure behaviors. Also, patients were taking either no MS medication or glatiramer acetate (GA) in order to minimize heterogeneity associated with medications, although this criterion was broadened after the current substudy was completed because of slow enrollment.

We asked subjects enrolled by March 2011 to participate in this substudy. 25-Hydroxyvitamin D levels (referred to herein as “vitamin D levels”) were obtained at the screening visit and, in all but 1 participant, after the study was completed. Participants were asked to collect their first morning stool at baseline and after 90 days of oral vitamin D3 5000 IU/d supplementation. Samples were shipped overnight on ice packs to the processing facility, where they were immediately stored at −80 °C. The UltraClean Fecal DNA Isolation Kit (MoBio, Carlsbad, CA) was used to batch-isolate total DNA, which was then amplified by polymerase chain reaction with bacterial 16S rRNA gene degenerate forward primer: 27 F.1 5′-AGRGTTTGATCMTGGCTCAG-3′ and a nondegenerate reverse primer: 1492R.jgi 5′-GGTTACCTTGTTACGACTT-3′. Amplified products were concentrated with a solid-phase reversible immobilization method and quantified with an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). PhyloChip Control Mix was added, followed by 35 cycles of bacterial 16S rRNA gene polymerase chain reaction amplification. The products were fragmented, labeled with biotin, and hybridized to the PhyloChip Array (version G3; Affymetrix).5The array contains representative sequences from operational taxonomic units (OTUs), analogous to a bacterial strain, which were clustered at 0.5% divergence using the Greengenes rRNA database.6A GeneArray scanner (Affymetrix) was used to wash, stain, and scan the PhyloChip arrays. Standard Affymetrix software was used to capture the scans. Hybridization values, the fluorescence intensity, for each taxon were calculated as a trimmed average; maximum and minimum values were removed prior to averaging. After maximum and minimum values were removed; the mean fluorescence intensity was log2 transformed before multiplying by 1000; doubling of the fluorescence intensity is thus indicated when the score changes by 1000.

Statistical Analyses

Based on the abundances of the OTUs determined by array hybridization, gut communities were compared by calculation of dissimilarities with weighted Unifrac. Dendrograms were generated by clustering the samples from the dissimilarity matrix hierarchically using the average-neighbor method. Two-dimensional ordination plotting by principal coordinates analysis was used to visualize intersample relationships by using the dissimilarity values to position the points relative to each other. The Prediction Analysis for Microarrays, using a nearest shrunken centroid method, was used to generate lists of significant taxa whose abundance characterizes each class.7The Adonis test was used to determine the significance of microbial differences associated with variables of interest.

To assess changes in OTU abundances associated with vitamin D supplementation, the within-person change in abundance was calculated for each OTU by OTUpost − OTUpre. Distributions of OTU changes in each group were compared using a Wilcoxon rank sum test with a Bonferroni correction, which adjusts for multiple comparisons.


We collected fecal samples from 15 subjects. A flow diagram depicting the sample sources is presented in Figure 1.The median (range) age of MS patients was 42 years (30–48 years) and of HCs was 38 years (29–51 years) (P = 0.95 [Wilcoxon rank sum test]). The mean 25-hydroxyvitamin D levels (measured in 1 batch by liquid chromatography–mass spectrometry at Heartland Assays [Ames, IA]) were similar between the 2 groups at baseline (25.9 ± 4.4 ng/mL in the MS subjects vs 23.2 ± 5.7 ng/mL in the HCs [P = 0.31, t test]) and after 90 days of supplementation (55.6 ± 17.0 ng/mL in the MS subjects vs 59.8 ± 11.7 ng/mL in the HCs [P = 0.58, t test]). When a Wilcoxon rank sum test was used, the differences also did not appear to be significant. Among MS patients (n = 7; 5 GA-treated), 1 submitted the initial stool sample late (excluded from pre–/post–vitamin D analyses), whereas 2 gave specimens only after supplementation. Thus, 4 MS subjects (2 GA-treated) had samples available before and after supplementation.


These samples were not available (n = 2) or not used (n = 1; sample given late) for the pre– versus post–vitamin D supplementation analysis. *One GA-treated subject (who did not contribute a pre–vitamin D supplementation sample) took antibiotics for an infection prior to the study (azithromycin for sinus infection), and **1 baseline sample was inadvertently stored at −20 °C for several days before being moved to −80 °C. There was no apparent impact on results of either affected sample (ie, they did not appear to be outliers).

Overall Comparison of MS Subjects Versus Healthy Controls

There was no apparent whole-community compositional shift in MS (Fig. 2; P = 0.15 when all samples were included; P = 0.74 when including only pre–vitamin D supplementation samples; P = 0.21 when including only post–vitamin D supplementation samples). Prior to supplementation, there were significantly fewer OTUs classified as Bacteroidaceae (37 OTUs) and Faecalibacterium (57 OTUs) and a higher number of Ruminococcus OTUs (38 OTUs) in MS compared with HCs (Supplementary Table 1, Using Prediction Analysis for Microarrays to identify a small subset of taxa with contrasting hybridization scores between the groups, OTUs within RikenellaceaeII, Lachnospiraceae, Porphyromonadaceae, Bacteroides, and Oscillibacter (Supplementary Table 2, collectively may predict whether samples prior to vitamin D supplementation came from MS patients or HCs.


PCoA by phenotype and supplementation. Principal component analysis is a transformation of Weighted Unifrac distance, a pairwise distance between samples based the calculation of the shared branches a the phylogenetic tree of the representative rRNA genes from 19,757 OTUs present in at least 1 sample, weighted by abundance. MS, affected; pre/post, refers to supplementation with vitamin D. Axis 1: 33% of variation explained. Axis 2: 21% of variation explained. All samples were used to generate the plots; for plot clarity, the excluded sample (late sample from 1 MS subject) was removed from this figure (see Supplementary Figure 1, for unaltered plot).

Overall Differences Between Treated Versus Untreated MS Subjects

Whole-community differences were observed when comparing GA-treated and untreated MS subjects using an Adonis test (P = 0.007) when all samples were considered (P = 0.66 when comparing only pre–vitamin D supplementation samples; P = 0.048 when comparing only post–vitamin D supplementation samples) (Fig. 3). Operational taxonomic units showing statistically significant different abundance between these groups (Supplementary Table 3, include Bacteroidaceae (500 OTUs), Faecalibacterium (186 OTUs), Ruminococcus (104 OTUs), Lactobacillaceae (77 OTUs), Clostridium (49 OTUs), and other Clostridiales (201 OTUs).


Differential clustering based on GA treatment. Principal component analysis is a transformation of Weighted Unifrac distance, a pairwise distance between samples based on the calculation of the shared branches a the phylogenetic tree of the representative rRNA genes from 19,757 OTUs present in at least 1 sample, weighted by abundance. MS, affected where patients are either on GA or untreated. Axis 1: 33% of variation explained. Axis 2: 21% of variation explained. All samples were used to generate the plots; for plot clarity, the excluded samples (HCs) were removed from this figure (see Supplementary Figure 2, for unaltered plot).

Effects of Vitamin D Supplementation

For the overall group, vitamin D supplementation was associated with a trend for a whole-community shift (P = 0.062), although the within-group P values were attenuated (P = 0.58 for the MS patients; P = 0.20 for HCs). Compared with baseline, in samples after supplementation, the abundance of Faecalibacterium (81 OTUs) and Enterobacteriaceae (244 OTUs) increased, whereas Ruminococcus decreased (53 OTUs) after vitamin D supplementation (Supplementary Table 4, At the individual level, when restricting analyses to only the OTUs found in both presupplementation and postsupplementation samples, most OTUs increased slightly in the HC samples following vitamin D supplementation (Fig. 4). In contrast, whereas the MS patients who were treated with GA showed few changes in OTU abundances, with the majority of OTUs remaining consistent in their abundance (median change near 0), untreated MS subjects showed a slight decrease of abundance for the majority of OTUs by Wilcoxon rank sum (P < 0.001) (Fig. 4).


Distribution of change of overall microbial abundance after vitamin D supplementation. Change in relative microbial OTU abundance in fecal samples between pre– and post–vitamin D supplementation samples in each subject group using only OTUs present in both samples. *Significant change in the distribution between HC and untreated.

While overall abundances did not markedly change in ways that distinguished the subject groups, a few specific OTUs distinguished each MS subgroup from the other 2 groups. Compared with HCs and GA-treated MS subjects, untreated MS subjects had an increase in the Akkermansia, Faecalibacterium, and Coprococcus genera after vitamin D supplementation (Fig. 5A). Compared with healthy participants and untreated MS subjects, those treated with GA had increases in Janthinobacterium and decreases in Eubacterium and Ruminococcus after high-dose vitamin D supplementation (Fig. 5B).


Distribution of statistically significant changes of OTU abundance in taxa. Change in relative abundance of taxa between the pre– and post–vitamin D supplement time points where statistically significant changes between (A) untreated MS and HCs or treated MS and (B) treated MS and untreated MS or HCs. *Significant change in the distribution between HC and untreated.


In this first exploratory study of gut bacteria in MS, the most important finding is that the gut bacterial communities differed in treated versus untreated MS patients. While we utilized all stool samples for the overall analysis because of the small sample size, which may lead to intersubject variability accounting for apparent differences, the difference remained statistically significant when considering only samples from MS subjects given after vitamin D supplementation. That the differences were not statistically significant before vitamin D supplementation probably is related to the small number of available samples from that time point, although we cannot exclude a true lack of difference between GA treated and untreated at that time. That MS patients receiving GA appear to have differences in gut bacteria compared with untreated patients is intriguing with respect to possible mechanisms of action of GA, although it clearly requires further investigation in a much larger cohort in order to establish the finding is replicable and to ascertain if the association is one of cause and effect. It is certainly plausible that GA may influence gut microbial populations. Glatiramer acetate has been shown to modulate the immune response in the gut in animal models of inflammatory bowel disease (which is known to be associated with gut microbial alterations), reducing the lymphocytic response to colonic extract in mesenteric lymph nodes.8The investigators of this study speculated that GA may locally displace microbial or self-antigen from antigen-presenting cells, thus reducing T-cell activation. Glatiramer acetate has also been shown to improve a mouse model of inflammatory bowel disease by inducing CD8+ T regulatory cells.9Regardless of the putative mechanism or direction of the relationship, the finding that the gut microbial profile may differ in treated versus untreated MS patients has immediate implications for design of studies of the gut microbiota in MS; investigators should carefully consider whether patients are treated or not and may even consider that different MS therapies could differentially impact gut microflora, in order to avoid heterogeneity and confounding.

In addition to noting the whole-community differences in GA-treated versus untreated MS patients, we distinguished that among MS patients, vitamin D supplementation altered the gut microbiota differently, depending on GA treatment status. Of great interest is that compared with healthy and GA-treated MS subjects, there was increased abundance of Faecalibacterium after vitamin D supplementation in MS subjects not on GA. Faecalibacterium is associated with a reduced inflammatory state due to its role in butyrate production,10and strategies to boost the abundance of Faecalibacterium have been hypothesized to be potentially therapeutic for patients with ulcerative colitis.11Several Coprococcus strains are also butyrate producers and may thus be anti-inflammatory.12 Akkermansia, which also increased in the untreated MS patients after vitamin D supplementation, is a mucin-degrading bacteria thought to be involved in immune tolerance toward commensal gut microbes, as evidenced by changes in gene expression of multiple pathways involving the immune response after germ-free mice were colonized by Akkermansia in 1 study.13In a mouse model of inflammatory bowel disease, Akkermansia-derived extracellular vesicles protected against the disease.14Little is known about the roles of Eubacteria (genera), Ruminococcus, and Janthinobacterium, all of which changed in abundance after vitamin D supplementation in treated MS subjects. The effect of vitamin D on the gut microbiota has been investigated in some mouse models of inflammatory bowel disease. Because bacteria do not possess vitamin D receptors, it is thought that any effects of vitamin D therapy are mediated by effects on the host that lead to secondary effects on gut bacteria.15In a mouse model of inflammatory bowel disease, knockout of Cyp27B1, the enzyme that converts 25-hydroxyvitamin D to the active 1, 25-dihydroxyvitamin D, is associated with elevated levels of proinflammatory Helicobacteraceae, which can induce inflammatory bowel disease symptoms.16The authors of the study suggest that disruption of vitamin D activity (through knockout models) increases the inflammatory milieu in the gut, leading to dysbiosis in which disease-inducing pathogens are able to outcompete commensal bacteria. While data regarding the effect of vitamin D in humans are sparse, our study suggests that further investigations of how vitamin D supplementation impacts gut microbial populations in MS patients are warranted.

This study also identified specific OTUs in bacterial genera that exhibited abundance differences in MS subjects compared with HCs, some of which were also identified in subjects with inflammatory bowel disorders.11,17Particularly of note the anti-inflammatory Faecalibacterium appeared to be less abundant in MS patients overall. This finding is consistent with the observation that Faecalibacterium is less abundant in patients with newly diagnosed, untreated Crohn disease than in HCs.18Patients with type 1 diabetes have been shown to have reduced butyrate-producing and mucin-degrading bacteria compared with controls such that if confirmed, the results herein may imply similarities in the mechanisms by which gut bacteria influence the development of autoimmune disease.19On the other hand, we were unable to define an overall MS-specific bacterial, community type. This finding could be viewed as consistent with studies of other autoimmune disease and supports the concept that instead of a major shift, an “unhealthy” community is relative to one’s “healthy” community of bacteria. However, it is still possible that the small sample size may have precluded detection of a specific overall MS community type.

The study has several limitations. The small sample size may preclude detecting other important differences, but the use of a conservative Bonferroni correction in this study limits concerns that significant results are due to chance. That only 1 sample was available per subject before and after supplementation may have introduced outliers; future studies should acquire more than 1 sample to represent each time point of interest. Whether the results would have differed if the patients had not had vitamin D insufficiency at baseline is unclear. Self-reported race and ethnicity were used, rather than genetic ancestry markers, to attempt to minimize heterogeneity in skin pigmentation that would have potentially confounded the results of the parent study. Larger studies may want to consider using such markers in the screening or analysis process to account for differences in host genetics that may influence a person’s microbial colonization but may also need to consider self-reported cultural groupings so as to account for differences in behaviors that may influence the microbiota.20Finally, the study would be strengthened by the use of metagenomic approaches. However, the study also has several strengths including the relative homogeneity of MS cases and the exclusion of conditions that might affect the gut microbiota.

Despite the study’s small size, our results have important implications. First, that certain OTUs may differ between healthy and MS subjects, and in particular that differences in some of these OTUs have also been reported in inflammatory bowel diseases or are associated with immune function, lends support for further evaluations of the gut microbiota in MS. Second, that MS therapies and vitamin D supplementation (widely used in current MS practice) may influence gut bacterial populations suggests that studies of the gut microbiota in MS should take these immunomodulatory therapies into account in their design. This study provides only preliminary data about the gut microbiota in MS, but unraveling the relationships of therapies to microbial populations in the gut and elsewhere in the body will be important to understanding the independent effects of these factors in MS.


View Abstract

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.