Next Article in Journal
The Chloroplast Genome of Wild Saposhnikovia divaricata: Genomic Features, Comparative Analysis, and Phylogenetic Relationships
Next Article in Special Issue
Microbial Richness of Marine Biofilms Revealed by Sequencing Full-Length 16S rRNA Genes
Previous Article in Journal
Tissue-Specific Variations in Transcription Factors Elucidate Complex Immune System Regulation
Previous Article in Special Issue
Short-Term Responses of Soil Microbial Communities to Changes in Air Temperature, Soil Moisture and UV Radiation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 5
10.3390/genes13050930
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gut Microbiota as a Potential Predictive Biomarker in Relapsing-Remitting Multiple Sclerosis

by
Vicente Navarro-López
1,2,3,*,
María Ángeles Méndez-Miralles
1,2,4,*,
Rosa Vela-Yebra
4,
Ana Fríes-Ramos
5,
Pedro Sánchez-Pellicer
2,
Beatriz Ruzafa-Costas
2,
Eva Núñez-Delegido
2,
Humberto Gómez-Gómez
2,
Sara Chumillas-Lidón
2,
Jose A. Picó-Monllor
2,6 and
Laura Navarro-Moratalla
2
1
Ph.D. Program in Health Sciences, Campus de los Jerónimos 135, UCAM-Universidad Católica San Antonio de Murcia, 30107 Murcia, Spain
2
MiBioPath Research Group, Department of Clinical Medicine, Campus de los Jerónimos 135, UCAM-Universidad Católica San Antonio de Murcia, 30107 Murcia, Spain
3
Infectious Disease Unit, University Hospital Vinalopó, Carrer Tonico Sansano Mora 14, 03293 Elche, Spain
4
Department of Neurology, University Hospital of Torrevieja, Carretera CV95, s/n, 03186 Alicante, Spain
5
Department of Neurology, University Hospital of Vinalopó, Carrer Tonico Sansano Mora 14, 03293 Elche, Spain
6
Department of Pharmacology, Pediatrics and Organic Chemistry, Faculty of Pharmacy, Universidad Miguel Hernández de Elche, 03202 Elche, Spain
*
Authors to whom correspondence should be addressed.
Genes 2022, 13(5), 930; https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050930
Submission received: 28 April 2022 / Revised: 17 May 2022 / Accepted: 20 May 2022 / Published: 23 May 2022
(This article belongs to the Special Issue When Genes Meet Microbial Ecology and Evolution)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background: The influence of the microbiome on neurological diseases has been studied for years. Recent findings have shown a different composition of gut microbiota detected in patients with multiple sclerosis (MS). The role of this dysbiosis is still unknown. Objective: We analyzed the gut microbiota of 15 patients with active relapsing-remitting multiple sclerosis (RRMS), comparing with diet-matched healthy controls. Method: To determine the composition of the gut microbiota, we performed high-throughput sequencing of the 16S ribosomal RNA gene. The specific amplified sequences were in the V3 and V4 regions of the 16S ribosomal RNA gene. Results: The gut microbiota of RRMS patients differed from healthy controls in the levels of the Lachnospiraceae, Ezakiella, Ruminococcaceae, Hungatella, Roseburia, Clostridium, Shuttleworthia, Poephyromonas, and Bilophila genera. All these genera were included in a logistic regression analysis to determine the sensitivity and the specificity of the test. Finally, the ROC (receiver operating characteristic) and AUC with a 95% CI were calculated and best-matched for Ezakiella (AUC of 75.0 and CI from 60.6 to 89.4) and Bilophila (AUC of 70.2 and CI from 50.1 to 90.4). Conclusions: There is a dysbiosis in the gut microbiota of RRMS patients. An analysis of the components of the microbiota suggests the role of some genera as a predictive factor of RRMS prognosis and diagnosis.

1. Introduction

Multiple sclerosis (MS) is a rare chronic autoimmune inflammatory and degenerative disease characterized by a multifocal demyelination of the central nervous system (CNS) [1]. Traditionally, incidence rates have been higher in areas away from the equator and in developed countries [2,3]. This geographical distribution has been related to differences in sun exposure and vitamin D levels, but it could also be related to the hygiene hypothesis that correlates excess sterilization, cleaning, vaccination, and the widespread use of antibiotics to the complication of proper development of immune responses, which could play an important role in the development of MS [4].
The disease typically evolves as relapses that consist of acute episodes of demyelination of the white matter, which manifest as an altered nerve conduction, secondary axonal damage, oligodendrocyte loss, and subsequent permanent damage of the involved CNS area [5]. This alteration is immunomodulated by myelin-autoreactive CD4+ T cells [6], although a recent report suggested that B cells could also be involved in the pathogenesis of this neurological disease [7]. The etiology of MS remains unknown, and it is postulated that it may developed in genetically predisposed individuals due to an alteration of the immune system triggered by an unknown environmental factor not yet identified. Some studies have evaluated the influence of the Epstein–Barr virus, smoking, low levels of vitamin D, or alcohol consumption [8]. Among them, a new one emerges as a potential factor in the etiology of MS: gut microbiota. An increasingly growing number of reports published over the last years have suggested the potential role of microbiota in MS [9].
MS is the most frequent demyelinating disease, which usually develops in young adults between the second and the third decade of life, being more frequent in females than in males [1]. The incidence is growing in western countries, but also the difference between both sexes (proportion of female/male of 2.3–3.5:1). There are some theories to try to explain these differences, but they probably are based on the difference in sexual hormones (progesterone and some estrogens are protectives factors in pregnancy), more susceptibility to vitamin D deficiency, other environmental risk factors, and more genetic predisposition with an overrepresentation of MS susceptibility after the random inactivation of one X chromosome in women [10].
To prevent disease activity and disability progression, new disease-modifying therapies (DMTs) have emerged based on controlling the altered immune response in MS patients. The first therapies that were effective, beta-interferon and glatiramer acetate, had immunomodulatory characteristics, just as teriflunomide did later (an immunosuppressive drug that interfered in pyrimidine synthesis), all with modest efficacy. Afterwards, some DMTs have appeared with different therapeutic immune targets. Dimethyl fumarate has anti-inflammatory and cytoprotective properties with moderate efficacy. Other DMTs work by blocking the access of activated lymphocytes to the CNS with two strategies: restricting their migration across the blood–brain barrier by blocking adhesion α4β1 integrin (natalizumab), or avoiding their egress from lymph nodes by inhibition of the sphingosine-1-phosphate receptor (fingolimod and siponimod). On the other hand, some DMTs (ocrelizumab and ofatumumab) are monoclonal antibodies that bind the CD20 antigen located on the surface of mature B cells and produce their depletion, causing chronic immunosuppression [11]. Lastly, immune reconstitution therapies have appeared which cause a more aggressive state of immunosuppression followed by reconstitution with a new profile of the immune system, including alemtuzumab, a monoclonal antibody binding the CD52 antigen expressed on the surface of T and B lymphocytes, and cladribine, a synthetic chlorinated deoxyadenosine analog that interferes with DNA synthesis and repair, causing the sustained reduction of circulating T and B lymphocytes [12].
Deciding which DMT should be provided and when for each patient is relevant to control disease progression and to improve the patient’s quality of life [13]. To date, there are scarce resources and tools in clinical practice that can be used for MS diagnosis and to predict an adequate response to treatment [14]. Most physicians rely on magnetic resonance imaging (MRI) scans to monitor the brain damage of patients under treatment. However, this approach requires monitoring brain damage for several months before determining treatment failure, a time where neuronal deterioration has already occurred. Additionally, other solutions exist, such as cerebrospinal fluid (CSF) analysis or blood tests, that are suboptimal due to their invasiveness or low specificity in predicting MS treatment response. Taking into consideration the importance of early personalized treatment to stop progression, gut microbiota may be an innovative approach in the identification of novel biomarkers used to predict MS diagnosis and the patient’s response to DMT. Gut microbiota has a population of around 100 trillion microorganisms mainly composed by bacteria, but it also contains archaea, eukaryotes, and viruses [15]. It begins to colonize the gastrointestinal tract at the moment of birth, and the initial composition depends on the mode of delivery, since the starting population comes from the maternal urogenital tract or the skin in vaginal or caesarean section delivery, respectively [16]. The composition of the microbiota is defined after our first years of life but could be changed by different environmental factors [17]. An increasing number of evidence is starting to correlate gut microbiota with several diseases, such as obesity [18] and inflammatory bowel disease [19], among others. In addition, it has been shown that gut microbiota plays a crucial role for the host metabolism and also in the development, maturation, and appropriate function of the immune system, preventing colonization by pathologic microorganisms [20] and helping to maintain the intestinal immune barrier by stimulating innate [21] and adaptive immune responses [22]. Moreover, the fact that MS is related with some gastrointestinal symptoms such as constipation, fecal incontinence, and intestinal permeability suggests a possible connection between gut microbiota and demyelinating diseases [23]. Recently, specific functions have started to be associated with different microorganisms, such as the potential role of Bacteroides fragilis to induce Foxp3+ Treg cells. These cells can ameliorate inflammation conditions in ulcerative colitis or experimental autoimmune encephalomyelitis (EAE), a murine model of MS [24]. In fact, several studies of EAE have shown a possible correlation between gut microbiota and MS. These studies have described how the use of broad-spectrum antibiotics is related with a slower progression of EAE and a decrease in proinflammatory cytokines and mesenteric Th17 cells [25], as well as how gut microbiota can influence the permeability of the blood–brain barrier [26]. These data have suggested that gut microbiota can modify a host’s immune response in EAE studies, influencing its origin, prevention, and possible treatment. Due to the increasing number of scientific studies that have identified gut dysbiosis in the stool samples of MS patients, we decided to investigate the possible role of the gut microbiota composition as a predictive factor of the disease in a prospective 24-month follow-up study with a cohort of patients with active relapsing-remitting multiple sclerosis (RRMS).

2. Material and Method

2.1. Study Population

The patients evaluated for inclusion in this study had previously been diagnosed with RRMS according to the 2017 McDonald criteria from clinical and radiological findings [27]. All the patients included were diagnosed with active RRMS, defined in our protocol as the presence of at least one relapse (defined as a new symptom or worsening of a previous one with objective findings reflecting focal or multifocal involvement of the CNS, developing acutely or subacutely, lasting greater than 24 h without a coexisting fever or infection, and followed by a full or incomplete recovery), at least one T1-gadolinium-enhancing (GdE) lesion, or new or unequivocally enlarging T2 lesions documented by magnetic resonance imaging (MRI) [28] in the two years prior to their inclusion in the study. The control group was selected from the healthy family members or partners who lived with the patients to homogenize both groups, especially matched by diet. All the participants provided sociodemographic and other clinical data, which included age, sex, body mass index (BMI), and the usual diet they followed during the last year. The exclusion criteria were the following: use of immunosuppressant drugs such as systemic corticosteroids, methotrexate, cyclosporine, or anti-tumor necrosis factor alpha (TNFα) prescribed for any comorbidity in the previous 3 months; systemic antibiotics or demonstration of a bacterial infection in the previous 2 weeks; and concomitant diagnosis of cirrhosis or intestinal bowel disease.
Furthermore, the MS patient group continued to be evaluated periodically for a 24-month follow-up period by an expert neurologist in MS. Moreover, brain MRI (and spinal cord MRI in some cases) was performed at least every 12 months to detect the presence of new lesions in the follow-up period. Expert neuroradiologists manually counted the number of new or unequivocally enlarging T2 lesions and T1-GdE lesions detected by MRI. Stool samples were obtained from all the participants and were immediately frozen and stored at −80 °C until they were processed for high-throughput sequencing. The Ethics Committees at the University Hospital of Vinalopó in Elche and the University Hospital of Torrevieja, Spain, approved the study in accordance with the Helsinki Declaration protocol, and all the subjects gave their written, informed consent before being included in the study.

2.2. High-Throughput Sequencing of Stool Samples of Patients with MS

DNA from the stool samples was isolated with the aid of a MagnaPure Compact System (Roche Life Science, Mannheim, Germany) to avoid bias in the DNA purification toward the misrepresentation of Gram-positive bacteria, following the Yuan et al., method [29]. For high-throughput sequencing, the hypervariable region of V3-V4 of the bacterial 16s ribosomal RNA gene was amplified using key-tagged eubacterial primers [30] and sequenced with a MiSeq Illumina Platform (Illumina, San Diego, CA, USA). The Illumina recommendations for library preparation and sequencing for metagenomics studies were used. The bacterial compositions of all the patients were compared with a group of healthy diet-matched individuals selected from the healthy partners living with them.

2.3. Bioinformatics

The resulting sequences were split considering the barcode introduced during the PCR reaction, while the R1 and R2 reads were overlapped using PEAR program version 0.9.122 (https://cme.h-its.org/exelixis/web/software/pear/doc.html, accessed on 4 September 2019), providing a single FASTQ file for each sample. Quality control of the sequences was performed using different steps. First, quality filtering (with a minimum threshold of Q20) was performed using the fastx tool kit version 0.013 (http://hannonlab.cshl.edu/fastx_toolkit/download.html, accessed on 4 September 2019). Finally, primer (16s rRNA primers) trimming and length selection (reads over 300nts) was performed with cutadapt version 1.4.123 (https://cutadapt.readthedocs.io/en/stable/, accessed on 4 September 2019). These FASTQ files were converted to FASTA files, and UCHIME program version 7.0.1001 (https://www.drive5.com/usearch/manual/uchime_algo.html, accessed on 4 September 2019) was used in order to remove chimeras that could arise during the amplification and sequencing step. These clean FASTA files were compared with BLAST24 against the NCBI 16s rRNA database using blast version 2.2.29+ (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 4 September 2019). The resulting XML files were processed using a python script from the Statistical Department of the Biotech Bioithas (Parc Científic de la Universitat d’Alacant, Alicante, Spain) in order to annotate each sequence at different phylogenetic levels (phylum, family, genera, and species).

2.4. Statistical Analysis

The descriptive quantitative variables were expressed as means ± SD (standard deviation) and as frequencies for qualitative variables. Comparisons between the groups were performed with a U Mann–Whitney test for quantitative data, while qualitative variables were analyzed using chi-squared and Fisher tests. An odds ratio (OR) with a 95% confidence interval (CI) was taken as a measure of effect size. A two-tailed p value of less than 0.05 was considered to indicate statistical significance. Finally, a logistic regression analysis model and the area under the receiver operating characteristic curve (AUCROC) using the variable case as dependent and the variable number of sequences of each genera as independent was performed. The predicted values were obtained, and the area under the ROC curve was estimated, as well as the cut-off point that maximized sensitivity and specificity. All the statistical analyses were performed using IBM SPSS statistics version 22 (SPSS Inc. Chicago, IL, USA) and R version 3.2.3 (https://cran.r-project.org/index.html, accessed on 14 October 2019). In the case of microbiota analysis, alpha diversity was conducted using a specaccum program in the vegan package, as implemented for R version 3.2.3.

3. Results

From May to June 2016, fifteen patients with a diagnosis of RRMS were enrolled in the study. Thirteen patients (86.6%) were female, with a mean age of 38.15 ± 8.08 years old. All the patients were under disease-modifying therapies (DMTs) when they were included in the study. However, 3/15 patients (20%) had stopped their DMTs weeks or months before they collected the stool samples because they were in a switching period until a new selected treatment started due to uncontrolled disease activity. Beta-interferon (6/15; 40.0%) and fingolimod (4/15; 26.7%) were the most frequent treatments administered. All the patients met the 2017 McDonald criteria of RRMS [27] at the time of inclusion in the study, myelitis being the most prevalent type of initial clinical presentation (6/15; 40.0%). The baseline characteristics of the RRMS patients are described in Table 1.
The compositions of the gut microbiota in all the patients were analyzed by high-throughput sequencing. This allowed us to describe “relapsing-remitting multiple sclerosis microbiota” based on the median values of the bacteria detected in all the patients (Supplementary Table S1). The analysis of the levels of species showed 112 different bacteria in the gut microbiota of the RRMS patient group. The most frequent genera detected among the analyzed stool samples were Bacteroides, Faecalibacterium, and Ruminococcus. Other common detected genera with mean values higher than 1% of the total stool samples were Gemmiger, Roseburia, Alistipes, Eubacterium, Oscillibacter, Lachnoclostridium, Barnesiella, Parabacteroides, Ruminiclostridium, Blautia, Phascolarctobacterium, Erysipelatoclostridium, and Bifidobacterium.
The bacterial compositions of all the RRMS patients were compared with a control group composed of healthy partners living with the patients to homogenize both groups by age and diet. The Shannon biodiversity index was calculated and then compared (mean, (interquartile range—IQR)) between the RRMS patients (3.05, (2.86–3.31)) and the healthy controls (2.43, (2.19–2.57)), showing significantly higher variability in the control group (p < 0.01) (Figure 1).
To search for possible genera that could allow the clustering of patient data, we performed a principal component analysis of each patient, comparing bacteria composition at the taxonomical genera level. No significant differences were found at the microbial structure level between both groups (Figure 2).
In addition, a logistic regression analysis was performed at the genera level of the microbiota composition and detected statistically significant differences between these two population groups with a p value < 0.05 in levels of the genera of Lachnospiraceae, Ezakiella, Ruminococcaceae, Hungatella, Roseburia, Clostridium, Shuttleworthia, Porphyromonas, and Bilophila (Table 2). The sensitivity and the specificity of the test, selecting the best cut-off point for the number of the sequences for each germ in the ROC curves (receiver operating characteristic) and analysis of the area under the curve (AUC) with a 95% CI, were calculated for all the above described genera and best-matched for Ezakiella (AUC 75.0; CI 60.6 to 89.4) and Bilophila (AUC 70.2; CI 50.1 to 90.4) (Figure 3 and Figure 4).
The Ezakiella genus was detected in eight (57.1%) cases and in one (7.1%) control. The median sequence numbers between the cases and controls were three and zero, respectively (p = 0.004). Regarding the Bilophila genus, sequences were found in 13 (92.9%) cases and in 11 (78.6%) controls, the median sequence numbers in cases and controls being 96 and 24, respectively (p = 0.036) (Supplementary Figure S1).
The different components of diet were registered for each patient and control included in the study. No difference was found in the comparative analysis of the main components of the diet in both study groups (Figure 5).

4. Discussion

During the last years, several studies have evaluated gut microbial composition in different cohorts of patients with MS. The results have suggested a different microbiota composition in MS patients compared to healthy populations [31,32,33,34,35]. In our study, this analysis showed an important dysbiosis in the group of patients with relapsing-remitting multiple sclerosis (RRMS) compared to healthy controls.
The most important contribution of our work was the study of the gut microbiota composition as a predictive factor of the disease course in patients with active RRMS. We found a statistically significant positive correlation between the number of sequences in the RRMS group and the levels of the Ezakiella and Bilophila genera. In the ROC curve, we found a cut-off point that established a sensitivity and specificity relation that better discriminated between cases of RRMS and healthy people.
Horton et al. [31], found an unspecified member of Coriobacteriales denominated SV_520, which was associated with an increased risk of experiencing disease activity (both clinical relapses and new lesions on MRI) in pediatric MS patients. On the other hand, their study showed some bacterial genera as protective factors, such as Butyricicoccus desmolans, Odoribacter splanchnicus, unidentified species in the Lachnospiraceae NK4A136 group, and Ruminococcaceae. Another study by Tremlett et al. [32] linked the absence of Fusobacteria with a shorter time of relapse, also in pediatric MS cases.
Among the 15 patients described here, we found decreased levels of Prevotella, as described previously in other studies with MS patients [33]. We also found low levels of Streptococcus genera, data previously reported in several studies that showed low butyrate production in the gut microbiomes of patients with inflammatory diseases [34]. Also, this data may be of interest considering the recent publication of Tankou et al. [35]. These authors published an interventional study evaluating the biological effect of the administration of an oral probiotic mix containing Streptococcus, among others, after the completion of 2 months of treatment. Although the number of MS patients included seems to be small, the investigators demonstrated an increase in the relative abundance of species in the gut microbiota and the induction of an anti-inflammatory peripheral immune response with a decrease in proinflammatory monocytes, as well as in activated monocyte and dendritic cell markers, present both in MS patients and healthy controls after the administration of the probiotic mix.
Cantarel et al., found that 15 patients with MS had decreased basal levels of the Bacteroidaceae family and the Faecalibacterium genera and increased levels of the Ruminococcus genera compared with a group of healthy subjects. After supplementation with vitamin D for 90 days, there were significant increased levels of the Faecalibacterium genera and Enterobacteriaceae family and decreased levels of the Ruminococcus genera when they were compared with healthy controls [36].
A larger study compared the gut microbiota of 60 MS patients with 43 healthy controls [37] and evidenced significant higher levels in archaea, especially Methanobrevibacter smithii and in the Akkermansia genera, among the group of MS patients. These bacteria have been associated with a potential role in inflammation. As in the Cantarel et al., study, the levels of phylum bacteroidetes were reduced in our MS patients, specifically a decrease in the butyrate producer, Butyricimonas. Due to the anti-inflammatory properties of butyrate, a reduction in the levels of these species of MS patients may contribute to disease, as occurs in other diseases where microbiota and intestinal permeability seem to be involved [38].
In another study, Tremlett et al. examined the gut microbiota of 18 MS pediatric patients [39] and found increased levels of Clostridium, Bilophila, Escherichia, and Shigella. Some strains of these bacteria have been associated with infection and inflammation [40,41]. Tremlett et al. also described a decrease in the levels of Eubacterium rectale and Corynebacterium. All these observational studies have concluded that MS patients have less biodiversity in the gut microbiota than healthy people and show significant variability in the levels of different bacteria genera, suggesting that may play a role in the inflammatory processes involved in the pathogenesis of the disease.
Finally, we must take some considerations with these results. Significant associations were found with the cross-reactivity of human GDP-L-fucose synthase peptides and immunodominant myelin basic protein peptides by CNS-infiltrating CD4+ T cells implicated in the myelin degradation process typical in MS patients [42]. These autoreactive CD4+T cells can also be activated by homologous GDP-L-fucose synthase from the gut microbiota, especially bacterial species related with MS pathogenesis in previous research [43]. These results coupled with the analysis here communicated and the possible cross-recognition of the GDP-L-fucose synthase enzyme from gut microbiota suggest a possible role of some components of this microbiota as a trigger of pathogenic autoimmune responses in MS. There is scientific evidence indicating that only a few gut bacteria activate this peptide [44]. All this information supports our hypothesis on the role of Bilophila and Ezakiella levels with disease activity in MS patients.
Another known factor that can affect and modify microbiota is diet. In fact, diet and age are the main factors that contribute to the composition of the microbiota. Although experimental information suggests that there are a variety of dietary interventions that could benefit patients with MS (vitamin D supplementation, low salt intake, etc.) [45], to date, its benefit has not been well-established in either relapsing or progressive MS. Data about the diet and its mean components were collected in both the patient and control groups in our study. Although the control group was made up mostly of the patients’ partners, having chosen nonconsanguineous controls who shared the same diet (no significant differences in the diets of both groups) brought us reliable data obtained through this control group (mainly male). Despite the fact that the influence of gender on the microbiota is important, the most determinant factor of microbiota composition is the diet [46].
So far, to our knowledge, there is no previous work that has studied the microbiota of patients with RRMS by comparing data with a control group matched by diet to avoid the bias of food intake on the analysis of gut microbiota. The results here communicated suggested a relation between higher levels of Ezakiella and Bilophila genera and RRMS when compared with data for healthy people. The positive correlation here described among these bacterial species of the gut microbiota and the relapses experienced by these patients despite adequate treatment may be relevant and could serve as a predictive biomarker of active RRMS and also a prognosis biomarker in RRMS patients. The main limitations of these results were the low number of analyzed cases and the greater percentage of male patients in the control group. These potential biases should be controlled in future clinical trials that include a greater number of patients and, among these, a greater presence of males to corroborate the data here presented.

5. Conclusions

Our study described the gut microbial composition of a case series of active relapsing-remitting MS patients and a control group. Diet data was registered in both groups, and the different components of diet and other related variables such as BMI were similar at baseline. The results suggested that high levels of Ezakiella and Bilophila genera could be a risk factor in patients with relapsing-remitting MS. The knowledge here communicated holds conceivable promise for developing novel prognosis biomarkers and more efficient therapeutic routes, including the modulation of gut microbiota in RRMS patients.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes13050930/s1, Figure S1: Analysis of cases and controls. Genera Ezakiella and Bilophila. Table S1: Microbiota composition at different taxonomic levels.

Author Contributions

Conceptualization, M.Á.M.-M., R.V.-Y. and V.N.-L.; methodology M.Á.M.-M.; R.V.-Y., A.F.-R. and V.N.-L.; software, H.G.-G., S.C.-L. and L.N.-M.; validation, V.N.-L. and H.G.-G.; formal analysis, H.G.-G. and L.N.-M.; investigation, P.S.-P.; resources, V.N.-L. and L.N.-M.; data curation P.S.-P., J.A.P.-M. and A.F.-R.; writing—original draft preparation, M.Á.M.-M. and V.N.-L.; writing—review and editing, E.N.-D., B.R.-C., R.V.-Y., A.F.-R., H.G.-G., J.A.P.-M., P.S.-P. and L.N.-M.; visualization, E.N.-D., B.R.-C., R.V.-Y., A.F.-R., S.C.-L., H.G.-G. and L.N.-M.; supervision, P.S.-P. and V.N.-L.; project administration, L.N.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research has received funding from European Union´s Horizon 2020 research and innovation programme under grant agreement No 867463.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committees for Clinical Research of the University Hospital of Vinalopó and the University Hospital of Torrevieja on 9 November 2017.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Calvo-Barreiro, L.; Eixarch, H.; Montalban, X.; Espejo, C. Combined therapies to treat complex diseases: The role of the gut microbiota in multiple sclerosis. Autoimmun. Rev. 2018, 17, 165–174. [Google Scholar] [CrossRef] [PubMed]
  2. Sawcer, S.; Franklin, R.J.M.; Ban, M. Multiple sclerosis genetics. Lancet Neurol. 2014, 13, 700–709. [Google Scholar] [CrossRef]
  3. Simpson, S.; Blizzard, L.; Otahal, P.; Van der Mei, I.; Taylor, B. Latitude is significantly associated with the prevalence of multiple sclerosis: A meta-analysis. J. Neurol. Neurosurg. Psychiatry 2011, 82, 1132–1141. [Google Scholar] [CrossRef] [PubMed]
  4. Fleming, J.; Fabry, Z. The hygiene hypothesis and multiple sclerosis. Ann. Neurol. 2007, 61, 85–89. [Google Scholar] [CrossRef] [PubMed]
  5. Mendibe, M.M.; Boyero, S.; Rodrigo, R.M.; Zarranz, J.J. Esclerosis múltiple y encefalitis autoinmunes. Leucodistrofias y otras enfermedades desmielinizantes. In Neurología, 5th ed.; Zarranz, J.J., Ed.; Elsevier: Madrid, Spain, 2013; pp. 451–475. [Google Scholar]
  6. Weiner, H.L. The challenge of multiple sclerosis: How do we cure a chronic heterogeneous disease? Ann. Neurol. 2009, 65, 239–248. [Google Scholar] [CrossRef] [PubMed]
  7. Krumbholz, M.; Derfuss, T.; Hohlfeld, R. B cells and antibodies in multiple sclerosis pathogenesis and therapy. Nat. Rev. Neurol. 2012, 8, 613–623. [Google Scholar] [CrossRef] [PubMed]
  8. Yadav, S.K.; Mindur, J.E.; Ito, K. Advances in the immunopathogenesis of multiple sclerosis. Curr. Opin. Neurol. 2015, 28, 206–219. [Google Scholar] [CrossRef]
  9. Chen, J.; Chia, N.; Kalari, K.R.; Yao, J.Z.; Novotna, M.; Paz, M.M.; Luckey, D.H.; Marietta, E.V.; Jeraldo, P.R.; Chen, X.; et al. Multiple sclerosis patients have a distinct gut microbiota compared to healthy controls. Sci. Rep. 2016, 6, 28484. [Google Scholar] [CrossRef] [Green Version]
  10. Hanne, F.H.; Gold, R.; Tintoré, M. Sex and gender issues in multiple sclerosis. Ther. Adv. Neurol. Disord. 2013, 6, 237–248. [Google Scholar]
  11. Hauser, S.L.; Cree, B.A.C. Treatment of multiple sclerosis: A review. Am. J. Med. 2020, 133, 1380–1390. [Google Scholar] [CrossRef]
  12. Lünemann, J.D.; Ruck, T.; Muraro, P.A.; Bar-Or, A.; Wiendl, H. Immune reconstitution therapies: Concepts for durable remission in multiple sclerosis. Nat. Rev. Neurol. 2019, 16, 56–62. [Google Scholar] [CrossRef] [PubMed]
  13. Río, J.; Nos, C.; Tintoré, M.; Borrás, C.; Galán, I.; Comabella, M.; Montalban, X. Assessment of different treatment failure criteria in a cohort of relapsing-remitting multiple sclerosis patients treated with interferon beta: Implications for clinical trials. Ann. Neurol. 2002, 52, 400–406. [Google Scholar] [CrossRef] [PubMed]
  14. Ziemssen, T.; Akgün, K.; Brück, W. Molecular biomarkers in multiple sclerosis. J. Neuroinflamm. 2019, 16, 272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Dave, M.; Higgins, P.D.; Middha, S.; Rioux, K.P. The human gut microbiome: Current knowledge, challenges, and future directions. Transl. Res. 2012, 160, 246–257. [Google Scholar] [CrossRef] [PubMed]
  16. Torrazza, R.M.; Neu, J. The developing intestinal microbiome and its relationship to health and disease in the neonate. J. Perinatol. 2011, 31, S29–S34. [Google Scholar] [CrossRef] [PubMed]
  17. Arumugam, M.; Raes, J.; Pelletier, E.; Le Paslier, D.; Yamada, T.; Mende, D.R.; Fernandes, G.R.; Tap, J.; Bruls, T.; Batto, J.M.; et al. Enterotypes of the human gut microbiome. Nature 2011, 473, 174–180. [Google Scholar] [CrossRef] [PubMed]
  18. Turnbaugh, P.J.; Hamady, M.; Yatsunenko, T.; Cantarel, B.L.; Duncan, A.; Ley, R.E.; Sogin, M.L.; Jones, W.J.; Roe, B.A.; Affourtit, J.P.; et al. A core gut microbiome in obese and lean twins. Nature 2009, 457, 480–484. [Google Scholar] [CrossRef] [Green Version]
  19. Ott, S.J.; Musfeldt, M.; Wenderoth, D.F.; Hampe, J.; Brant, O.; Fölsch, U.R.; Timmis, K.N.; Schreiber, S. Reduction in diversity of the colonic mucosa associated bacterial microflora in patients with active inflammatory bowel disease. Gut 2004, 53, 685–693. [Google Scholar] [CrossRef] [Green Version]
  20. Stecher, B.; Chaffron, S.; Käppeli, R.; Hapfelmeier, S.; Freedrich, S.; Weber, T.C.; Kirundi, J.; Suar, M.; McCoy, K.D.; von Mering, C.; et al. Like will to like: Abundances of closely related species can predict susceptibility to intestinal colonization by pathogenic and commensal bacteria. PLoS Pathog. 2010, 6, e1000711. [Google Scholar] [CrossRef]
  21. Clarke, T.B.; Davis, K.M.; Lysenko, E.S.; Zhou, A.Y.; Yu, Y.; Weiser, J.N. Recognition of peptidoglycan from the microbiota by Nod1 enhances systemic innate immunity. Nat. Med. 2010, 16, 228–231. [Google Scholar] [CrossRef] [Green Version]
  22. Berer, K.; Krishnamoorthy, G. Commensal gut flora and brain autoimmunity: A love or hate affair? Acta Neuropathol. 2012, 123, 639–651. [Google Scholar] [CrossRef] [PubMed]
  23. Gupta, G.; Gelfand, J.M.; Lewis, J.D. Increased risk for demyelinating diseases in patients with inflammatory bowel disease. Gastroenterology 2005, 129, 819–826. [Google Scholar] [CrossRef] [PubMed]
  24. Ochoa-Reparaz, J.; Mielcarz, D.W.; Ditrio, L.E.; Burroughs, A.R.; Begum-Haque, S.; Dasgupta, S.; Kasper, D.L.; Kasper, L.H. Central Nervous System Demyelinating Disease Protection by the Human Commensal Bacteroides fragilis Depends on Polysaccharide A Expression. J. Immunol. 2010, 185, 4101–4108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Yokote, H.; Miyake, S.; Croxford, J.L.; Oki, S.; Mizusawa, H.; Yamamura, T. NKT cell-dependent amelioration of a mouse model of multiple sclerosis by altering gut flora. Am. J. Pathol. 2008, 173, 1714–1723. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Braniste, V.; Al-Asmakh, M.; Kowal, C.; Anuar, F.; Abbaspour, A.; Tóth, M.; Korecka, A.; Bakocevic, N.; Ng, L.G.; Kundu, P.; et al. The gut microbiota influences blood-brain barrier permeability in mice. Sci. Transl. Med. 2014, 158, 263ra158. [Google Scholar] [CrossRef] [Green Version]
  27. Thompson, A.J.; Banwell, B.L.; Barkhof, F.; Carroll, W.M.; Coetzee, T.; Comi, G.; Correale, J.; Fazekas, F.; Filippi, M.; Freedman, M.S.; et al. Diagnosis of multiple sclerosis: 2017 revisions of the McDonald criteria. Lancet Neurol. 2018, 17, 162–173. [Google Scholar] [CrossRef]
  28. McGinley, M.P.; Goldschmidt, C.H.; Rae-Grant, A.D. Diagnosis and treatment of multiple sclerosis: A review. JAMA 2021, 325, 765–779. [Google Scholar] [CrossRef]
  29. Yuan, S.; Cohen, D.B.; Ravel, J.; Abdo, Z.; Forney, L.J. Evaluation of methods for the extraction and purification of DNA from the human microbiome. PLoS ONE 2012, 7, e33865. [Google Scholar] [CrossRef] [Green Version]
  30. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glöckner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41, e1. [Google Scholar] [CrossRef]
  31. Horton, M.K.; McCauley, K.; Fadrosh, D.; Fujimura, K.; Graves, J.; Ness, J.; Wheeler, Y.; Gorman, M.P.; Benson, L.A.; Weinstock-Guttman, B.; et al. Gut microbiome is associated with multiple sclerosis activity in children. Ann. Clin. Transl. Neurol. 2021, 8, 1867–1883. [Google Scholar] [CrossRef]
  32. Tremlett, H.; Fadrosh, D.W.; Faruqi, A.A.; Hart, J.; Roalstad, S.; Graves, J.; Lynch, S.; Waubant, E.; Aaen, G.; Belman, A.; et al. Gut microbiota composition and relapse risk in pediatric MS: A pilot study. J. Neurol. Sci. 2016, 363, 153–157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Shahi, S.K.; Freedman, S.N.; Mangalam, A.K. Gut microbiome in multiple sclerosis: The players involved and the roles they play. Gut Microbes 2017, 8, 607–615. [Google Scholar] [CrossRef]
  34. Forbes, J.D.; Van Domselaar, G.; Bernstein, C.N. The gut microbiota in immune-mediated inflammatory diseases. Front. Microbiol. 2016, 7, 1081. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Tankou, S.K.; Regev, K.; Healy, B.C.; Tjon, E.; Laghi, L.; Cox, L.M.; Kivisäkk, L.; Pierre, I.V.; Hrishikesh, L.; Gandhi, R.; et al. A probiotic modulates the microbiome and immunity in multiple sclerosis. Ann. Neurol. 2018, 83, 1147–1161. [Google Scholar] [CrossRef] [PubMed]
  36. Cantarel, B.L.; Waubant, E.; Chehoud, C.; Kuczynski, J.; DeSantis, T.Z.; Warrington, J.; Venkatesan, A.; Fraser, C.M.; Mowry, E.M. Gut microbiota in multiple sclerosis: Possible influence of immunomodulators. J. Investig. Med. 2015, 63, 729–734. [Google Scholar] [CrossRef]
  37. Jangi, S.; Gandhi, R.; Cox, L.M.; Li, N.; von Glehn, F.; Yan, R.; Patel, B.; Mazzola, M.A.; Liu, S.; Glanz, B.L.; et al. Alterations of the human gut microbiome in multiple sclerosis. Nat. Commun. 2016, 7, 12015. [Google Scholar] [CrossRef] [Green Version]
  38. Geirnaert, A.; Calatayud, M.; Grootaert, C.; Laukens, D.; Devriese, S.; Smagghe, G.; De Vos, M.; Boon, B.; Van de Wiele, T. Butyrate-producing bacteria supplemented in vitro to Crohn’s disease patient microbiota increased butyrate production and enhanced intestinal epithelial barrier integrity. Sci. Rep. 2017, 7, 11450. [Google Scholar] [CrossRef] [Green Version]
  39. Tremlett, H.; Fadrosh, D.W.; Faruqi, A.A.; Zhu, F.; Hart, J.; Roalstad, S.; Graves, J.; Lynch, S.; Waubant, E.; US Network of Pediatric MS Centers. Gut microbiota in early pediatric multiple sclerosis: A case-control study. Eur. J. Neurol. 2016, 23, 1308–1321. [Google Scholar] [CrossRef] [Green Version]
  40. Issa, M.; Ananthakrishnan, A.N.; Binion, D.G. Clostridium difficile and inflammatory bowel disease. Inflamm. Bowel. Dis. 2008, 14, 1432–1442. [Google Scholar] [CrossRef]
  41. Sansonetti, P.J. Rupture, Invasion and Inflammatory Destruction of the Intestinal Barrier by Shigella: The Yin and Yang of Innate Immunity. Can. J. Infect. Dis. Med. Microbiol. 2006, 17, 117–119. [Google Scholar] [CrossRef] [Green Version]
  42. Planas, R.; Santos, R.; Tomas-Ojer, P.; Cruciani, C.; Lutterotti, A.; Faigle, W.; Schaeren-Wiemers, N.; Espejo, C.; Eixarch, H.; Pinilla, C.; et al. GDP-l-fucose synthase is a CD4T+ T cell-specific autoantigen in DRB3*02:02 patients with multiple sclerosis. Sci. Transl. Med. 2018, 10, eaat4301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Mangalam, A.K.; Yadav, M.; Yadav, R. The emerging world of microbiome in autoimmune disorders: Opportunities and challenges. Indian J. Rheumatol. 2021, 16, 57. [Google Scholar] [CrossRef] [PubMed]
  44. Sabatino, J.J.; Zamvil, S.S. T cells take aim at a ubiquitous autoantigen in multiple sclerosis. Sci. Transl. Med. 2018, 10, eaau8826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Van den Hoogen, W.J.; Laman, J.D.; Hart, B.A. Modulation of Multiple Sclerosis and Its Animal Model Experimental Autoimmune Encephalomyelitis by Food and Gut Microbiota. Front. Immunol. 2017, 8, 1081. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Cuevas-Sierra, A.; Romo-Hualde, A.; Aranaz, P.; Goni, L.; Cuervo, M.; Martínez, J.A.; Milagro, F.I.; Riezu-Boj, J.I. Diet-and sex-related changes of gut microbiota composition and functional profiles after 4 months of weight loss intervention. Eur. J. Nutr. 2021, 60, 3279–3301. [Google Scholar] [CrossRef]
Genes 13 00930 g001 550
Figure 1. Boxplot based on Shannon diversity index: a comparison between patients with RRMS and healthy controls.
Figure 1. Boxplot based on Shannon diversity index: a comparison between patients with RRMS and healthy controls.
Genes 13 00930 g001
Genes 13 00930 g002 550
Figure 2. Principal component analysis of the microbiota composition comparing healthy controls and RRMS patients. Comparison of group samples was based on the variability of the bacteria composition at the genera taxonomical level for each RRMS patient and healthy controls.
Figure 2. Principal component analysis of the microbiota composition comparing healthy controls and RRMS patients. Comparison of group samples was based on the variability of the bacteria composition at the genera taxonomical level for each RRMS patient and healthy controls.
Genes 13 00930 g002
Genes 13 00930 g003 550
Figure 3. ROC curve (receiver operating characteristic), AUC (area under the curve), and the selection of the best cut-off point with sensitivity and specificity in Ezakiella genera.
Figure 3. ROC curve (receiver operating characteristic), AUC (area under the curve), and the selection of the best cut-off point with sensitivity and specificity in Ezakiella genera.
Genes 13 00930 g003
Genes 13 00930 g004 550
Figure 4. ROC curve (receiver operating characteristic), AUC (area under the curve), and the selection of the best cut-off point with sensitivity and specificity in Bilophila genus.
Figure 4. ROC curve (receiver operating characteristic), AUC (area under the curve), and the selection of the best cut-off point with sensitivity and specificity in Bilophila genus.
Genes 13 00930 g004
Genes 13 00930 g005 550
Figure 5. Dietary registry of RRMS patients and healthy controls. Percentage of the main dietary components recorded in the group of patients with active RRMS and healthy controls included in the study.
Figure 5. Dietary registry of RRMS patients and healthy controls. Percentage of the main dietary components recorded in the group of patients with active RRMS and healthy controls included in the study.
Genes 13 00930 g005
Table
Table 1. Demographics of RRMS patients and initial clinical presentation and treatment at baseline.
Table 1. Demographics of RRMS patients and initial clinical presentation and treatment at baseline.
PatientGenderAge (Years)Weight (kg)Height (m)BMI (kg/m2)DMTInitial Clinical PresentationEDSSRelapsesMRI Lesions
1Female38901.7230.42beta-interferonBrainstem2.500
2Female18571.7119.49beta-interferon Myelitis111
3Female38571.5422.77natalizumabMyelitis2.500
4Female35591.5125.88fingolimodMyelitis100
5Female351071.7833.77dimethyl fumarateMyelitis, brainstem3.510
6Female3475,51.6029.49beta-interferon Optic neuritis336
7Female311001.8030.86fingolimodHemiparesis3.525
8Male3385.51.7029.58natalizumabMyelitis6.500
9Female45591.6521.67fingolimodMyelitis2.528
10Female5392.71.5538.58beta-interferon Brainstem3.502
11Female37631.6922.05beta-interferon Brainstem1.500
12Male44921.7330.77teriflunomideHemiparesis3. 523
13Female29541.6021.09dimethyl fumarateHemiparesis103
14Female51631.6622.86beta-interferon Brainstem, optic neuritis2.514
15Female37591.5125.88glatiramer acetateBrainstem200
RRMS: relapsing-remitting multiple sclerosis; BMI: body mass index; DMT: disease-modifying therapy; EDSS: expanded disability status scale; MRI: magnetic resonance imaging.
Table
Table 2. Case/control ratio by genera with statistical significance difference.
Table 2. Case/control ratio by genera with statistical significance difference.
GeneraMultiplier 95% CIp-ValueAUC (95% CI)
Lachnospiraceae0.06 (0.01 to 0.35)0.00167.1 (46.7 to 87.6)
Ezakiella13.08 (2.16 to 79.24)0.00577.0 (62.0 to 92.0)
Ruminococcaceae12.65 (2.09 to 76.62)0.00668.1 (47.9 to 88.3)
Hungatella8.57 (1.41 to 51.89)0.0268.3 (51.3 to 85.4)
Roseburia0.12 (0.02 to 0.73)0.0268.1 (50.2 to 86.0)
Clostridium7.26 (1.20 to 43.99)0.0367.4 (47.5 to 87.2)
Shuttleworthia6.32 (1.04 to 38.29)0.0462.4 (42.5 to 82.3)
Ruminococcaceae6.23 (1.03 to 37.71)0.0462.4 (42.1 to 82.7)
Porphyromonas6.00 (0.99 to 36.33)0.0564.3 (46.0 to 82.5)
Bilophila5.93 (0.98 to 35.91)0.0573.0 (54.0 to 93.0)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Navarro-López, V.; Méndez-Miralles, M.Á.; Vela-Yebra, R.; Fríes-Ramos, A.; Sánchez-Pellicer, P.; Ruzafa-Costas, B.; Núñez-Delegido, E.; Gómez-Gómez, H.; Chumillas-Lidón, S.; Picó-Monllor, J.A.; et al. Gut Microbiota as a Potential Predictive Biomarker in Relapsing-Remitting Multiple Sclerosis. Genes 2022, 13, 930. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050930

AMA Style

Navarro-López V, Méndez-Miralles MÁ, Vela-Yebra R, Fríes-Ramos A, Sánchez-Pellicer P, Ruzafa-Costas B, Núñez-Delegido E, Gómez-Gómez H, Chumillas-Lidón S, Picó-Monllor JA, et al. Gut Microbiota as a Potential Predictive Biomarker in Relapsing-Remitting Multiple Sclerosis. Genes. 2022; 13(5):930. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050930

Chicago/Turabian Style

Navarro-López, Vicente, María Ángeles Méndez-Miralles, Rosa Vela-Yebra, Ana Fríes-Ramos, Pedro Sánchez-Pellicer, Beatriz Ruzafa-Costas, Eva Núñez-Delegido, Humberto Gómez-Gómez, Sara Chumillas-Lidón, Jose A. Picó-Monllor, and et al. 2022. "Gut Microbiota as a Potential Predictive Biomarker in Relapsing-Remitting Multiple Sclerosis" Genes 13, no. 5: 930. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13050930

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop