<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="1.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Immunol.</journal-id>
<journal-title-group>
<journal-title>Frontiers in Immunology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Immunol.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">1664-3224</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fimmu.2026.1747643</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Research</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Dysbiosis of intestinal microbiota in patients with neuromyelitis optica spectrum disorders</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Du</surname><given-names>Qin</given-names></name>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname><given-names>Xiaofei</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/1901077/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Shi</surname><given-names>Ziyan</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/2308927/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Chen</surname><given-names>Hongxi</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/1674130/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhang</surname><given-names>Ying</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/3305099/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname><given-names>Rui</given-names></name>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Mou</surname><given-names>Zichao</given-names></name>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Kong</surname><given-names>Lingyao</given-names></name>
<uri xlink:href="https://loop.frontiersin.org/people/1993188/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Zhou</surname><given-names>Hongyu</given-names></name>
<xref ref-type="corresp" rid="c001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/776977/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Funding acquisition" vocab-term-identifier="https://credit.niso.org/contributor-roles/funding-acquisition/">Funding acquisition</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
</contrib-group>
<aff id="aff1"><institution>Department of Neurology, West China Hospital, Sichuan University</institution>, <city>Chengdu</city>, <state>Sichuan</state>,&#xa0;<country country="cn">China</country></aff>
<author-notes>
<corresp id="c001"><label>*</label>Correspondence: Hongyu Zhou, <email xlink:href="mailto:zhouhy@scu.edu.cn">zhouhy@scu.edu.cn</email></corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-02-27">
<day>27</day>
<month>02</month>
<year>2026</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2026</year>
</pub-date>
<volume>17</volume>
<elocation-id>1747643</elocation-id>
<history>
<date date-type="received">
<day>16</day>
<month>11</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>04</day>
<month>02</month>
<year>2026</year>
</date>
<date date-type="rev-recd">
<day>28</day>
<month>01</month>
<year>2026</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2026 Du, Wang, Shi, Chen, Zhang, Wang, Mou, Kong and Zhou.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>Du, Wang, Shi, Chen, Zhang, Wang, Mou, Kong and Zhou</copyright-holder>
<license>
<ali:license_ref start_date="2026-02-27">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</license-p>
</license>
</permissions>
<abstract>
<sec>
<title>Objective</title>
<p>This study aimed to explore the specific microbial signatures and metabolomic profiles of fecal microbiota in patients with neuromyelitis optica spectrum disorders (NMOSD) and assess the effects of immunosuppressants on their gut microbiota using a longitudinal cohort study.</p>
</sec>
<sec>
<title>Methods</title>
<p>We enrolled 21 treatment-na&#xef;ve NMOSD patients and 21 matched healthy controls (HCs). Fecal microbial composition and metabolomic profiles were compared between groups using 16S rRNA gene sequencing and ultra-high-performance liquid chromatography-mass spectrometry. Subsequently, fecal samples from NMOSD patients were collected and reassessed after immunosuppressant treatment.</p>
</sec>
<sec>
<title>Results</title>
<p>The gut microbial composition and metabolomic profiles of NMOSD patients were distinct from those of HCs. The &#x3b1;-diversity metrics were significantly higher in NMOSD patients than in HCs (P &lt;0.001). Microbiome alterations in NMOSD patients were characterized by increased abundances of <italic>Streptococcus</italic> and <italic>Ruminococcus</italic>, and decreased abundances of <italic>Faecalibacterium</italic>, <italic>Ralstonia</italic>, and <italic>Pseudomonas</italic> at the genus level (all with linear discriminant analysis scores &gt; 4 and P &lt; 0.001). Additionally, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States analysis identified 19 differentially abundant metabolites and 44 altered metabolic pathways in NMOSD patients compared to HCs. Immunosuppressive treatment for over six months may reduce these differences, shifting the gut microbiota composition and metabolite profiles of NMOSD patients closer to those of HCs.</p>
</sec>
<sec>
<title>Interpretation</title>
<p>Our study revealed significant gut microbiome dysbiosis and metabolic abnormalities in patients with NMOSD, which were markedly alleviated after six months of immunosuppressive treatment. These preliminary findings suggest the gut microbiota biomarkers could serve as potential therapeutic targets in the future.</p>
</sec>
</abstract>
<kwd-group>
<kwd>dysbiosis</kwd>
<kwd>gut microbiota</kwd>
<kwd>immunosuppressant</kwd>
<kwd>metabolomic profiles</kwd>
<kwd>neuromyelitis optica spectrum disorders</kwd>
</kwd-group>
<funding-group>
<funding-statement>The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the National Natural Science Foundation of China (grant numbers 82201494 and 82471376).</funding-statement>
</funding-group>
<counts>
<fig-count count="9"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="37"/>
<page-count count="11"/>
<word-count count="4775"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Multiple Sclerosis and Neuroimmunology</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>Neuromyelitis optica spectrum disorder (NMOSD) is a severe inflammatory demyelinating disease of the central nervous system, characterized by recurrent attacks of transverse myelitis and optic neuritis (<xref ref-type="bibr" rid="B1">1</xref>). Although the etiology of NMOSD remains elusive, aquaporin-4 (AQP4) has been identified as a disease-specific serum autoantibody marker (<xref ref-type="bibr" rid="B2">2</xref>). Growing evidence indicates that environmental factors play a significant role in the pathogenesis and progression of inflammatory and autoimmune diseases.</p>
<p>The gut microbiota, often called the second brain, may influence brain activity through the gut-microbiota-brain axis under both physiological and pathological conditions. The gut microbiome contributes substantially to these disorders by affecting the immune system and metabolic pathways (<xref ref-type="bibr" rid="B3">3</xref>). Dysbiosis of gastrointestinal microbiota affects the differentiation of proinflammatory T-cells and may enhance organ-specific autoimmunity (<xref ref-type="bibr" rid="B4">4</xref>). Extensive homology between gut microbiota and AQP4 protein demonstrates molecular mimicry in the pathogenesis of NMOSD (<xref ref-type="bibr" rid="B5">5</xref>).Compared with healthy controls (HCs), antibody responses against <italic>Helicobacter pylori</italic> and gastrointestinal antigens were more frequently observed in NMOSD patients, suggesting that alterations in the gastrointestinal environment may contribute to the initiation or progression of NMOSD (<xref ref-type="bibr" rid="B6">6</xref>).Our previous research found that NMOSD patients had increased abundances of the pathogenic genera <italic>Streptococcus</italic> and <italic>Flavonifractor</italic> compared with HCs (<xref ref-type="bibr" rid="B7">7</xref>).</p>
<p>Previous studies have primarily examined group differences at a cross-sectional level, which may not fully account for treatment-related effects. To address this limitation, we conducted a case-control study comparing gut microbiome and metabolomic profiles between treatment-na&#xef;ve NMOSD patients and HCs. Additionally, we performed longitudinal sampling (before and six months after immunosuppressant therapy) to evaluate how these medications alter the gut microbiota in NMOSD patients.</p>
</sec>
<sec id="s2">
<title>Methods</title>
<sec id="s2_1">
<title>Study design and participants</title>
<p>A total of 42 subjects, comprising 21 treatment-na&#xef;ve NMOSD patients (no treatment, NT group) and 21 HCs, were consecutively recruited in our study from West China Hospital, Sichuan University between January 2021 and June 2023. Patients were defined as treatment-na&#xef;ve under two criteria: 1) being immunotherapy-na&#xef;ve at initial diagnosis; or 2) having discontinued all immunomodulatory therapies after a continuous use of less than six months and maintaining a treatment-free washout period for at least six months. All patients met the 2015 international diagnostic criteria for NMOSD and were seropositive for AQP4-IgG as determined by cell-based assays (<xref ref-type="bibr" rid="B1">1</xref>, <xref ref-type="bibr" rid="B8">8</xref>). Importantly, none of the patients had received any immunotherapy at enrollment. HCs were matched for age, gender, body mass index (BMI), dietary habits, and geographical location. The detailed inclusion and exclusion criteria are presented in <xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>. We excluded subjects with acute infections, gastrointestinal disorders, malignancy, other autoimmune diseases, or those who had taken probiotics, antibiotics, or corticosteroids within one month prior to enrollment.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Flowchart of participant enrollment. The study included patients with NMOSD and HCs. NMOSD patients underwent assessments both before and after immunosuppressive therapy. HC, healthy control; NMOSD, neuromyelitis optica spectrum disorder.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g001.tif">
<alt-text content-type="machine-generated">Flowchart showing study participant selection: 21 NMOSD patients' stool samples collected before immunotherapy, with exclusions for immunosuppressant status, yielding 5 for post-immunotherapy; 21 healthy controls added; three final groups included.</alt-text>
</graphic></fig>
<p>Immunosuppressant in our study was mycophenolate mofetil (MMF) (20-30mg/kg/day). Demographic information and clinical characteristics, including age, gender, BMI, and other autoantibodies were collected at enrollment. The Medical Ethics Committee of West China Hospital, Sichuan University approved the study, and all participants provided informed consent prior to their inclusion in this study.</p>
</sec>
<sec id="s2_2">
<title>Fecal sample collection and 16S rRNA gene sequencing</title>
<p>All collected fecal samples were immediately cryopreserved at -80&#xb0;C until processing. Following thawing, genomic DNA was extracted and amplified by PCR using primers (338F: 5&#x2019;-ACTCCTACGGGAGGCAGCAG-3&#x2019;; 806R: 5&#x2019;-GGACTA CHVGGGTWTCTAAT-3&#x2019;) targeting the V3-V4 hypervariable regions of the bacterial 16S rRNA gene, as previously described (<xref ref-type="bibr" rid="B9">9</xref>). The PCR amplicons were purified using an AxyPrep DNA Gel Extraction Kit (Axygen, USA) and prepared for sequencing with the Illumina TruSeq DNA PCR-Free Library Preparation Kit to construct 16S rRNA gene libraries.</p>
</sec>
<sec id="s2_3">
<title>Data processing and bioinformatics analysis</title>
<p>Sequencing was performed on Illumina MiSeq platforms, and the resulting data were analyzed using Quantitative Insights Into Microbial Ecology (QIIME, v2.0) with default parameters for Illumina processing (<xref ref-type="bibr" rid="B10">10</xref>). Raw sequences underwent quality filtering prior to assembly, removing: (1) reads with primer mismatches, (2) reads overlaps containing &gt; 5% mismatches, and (3) reads shorter than 100 bp. Qualified paired-end reads were merged using FLASH (v1.2.11) (<xref ref-type="bibr" rid="B11">11</xref>). followed by chimera removal with UCHIME. High-quality sequences were clustered into operational taxonomic units (OTUs) at a &#x2265; 97% similarity threshold using VSEARCH (v2.15.0) (<xref ref-type="bibr" rid="B12">12</xref>), with representative sequences selected from each cluster.</p>
<p>Microbial &#x3b1;-diversity was evaluated using Observed species and Chao1 indices (richness), and Shannon and Simpson indices (diversity) (<xref ref-type="bibr" rid="B13">13</xref>). &#x3b2;-diversity was evaluated to examine microbial community heterogeneity between NMOSD patients and HCs through principal coordinates analysis (PCoA) based on Bray-Curtis distances and permutational multivariate analysis of variance (PERMANOVA; adonis2 function, 999 permutations).</p>
<p>Linear discriminant analysis effect size (LEfSe) was performed to identify differentially abundant taxa, applying a non-parametric Kruskal-Wallis test (&#x3b1;=0.05) followed by Linear discriminant analysis (LDA) with an effect size threshold of 4.0 (<xref ref-type="bibr" rid="B14">14</xref>). Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) was employed to predict metagenomic functional content, with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis performed to identify metabolic pathways associated with taxonomic composition differences (<xref ref-type="bibr" rid="B15">15</xref>).</p>
</sec>
<sec id="s2_4">
<title>Sample collection and metabolomics profiling analysis</title>
<p>Fresh fecal samples (around 400 mg) were collected upon enrollment, suspended in preservation buffer (Longseegen Stool Storage Kit, Longsee Biomedical Corp., China) following manufacturer&#x2019;s protocols, and cryopreserved at -80&#xb0;C. After homogenization, 200 &#xb5;L of fecal sample was vacuum-dried and resuspended in 800 &#xb5;L methanol. All samples were vortexed for 30 s, sonicated for 10 min, and incubated at &#x2212;20 &#xb0;C for 2 h to precipitate proteins. After centrifugation (13,000 rpm, 4&#xb0;C, 15 min), the supernatants were collected, vacuum-dried, and reconstituted in 200 &#xb5;L methanol/water (1:1, v/v).</p>
<p>Metabolomics analysis was performed by ultra-high-performance liquid chromatography and mass spectrometry (UHPLC-MS). The datasets were analyzed through pattern recognition methods using MetaboAnalyst 3.0. Univariate analysis (t-test) identified statistically significant features between groups, while multivariate partial-least-squares discrimination analysis (PLS-DA) assessed inter-subject variability and highlighted key discriminatory metabolites. Metabolites with PLS-DA-derived VIP &gt; 1.0 were prioritized as key contributors. Pathway analysis of these differential metabolites was performed using the KEGG database (<xref ref-type="bibr" rid="B16">16</xref>).</p>
</sec>
<sec id="s2_5">
<title>Statistical analysis</title>
<p>Quantitative data are described as the median (range) or mean &#xb1; SD (standard deviation). The Mann-Whitney U test or Student&#x2019;s t-test were used to compare variables between groups. The false discovery rate (FDR) was calculated to correct for multiple comparisons and assess screened difference variables (<xref ref-type="bibr" rid="B17">17</xref>). All analyses were performed using GraphPad Prism (v8.0) (GraphPad Software, Inc., San Diego, CA, USA) or SPSS Statistics software (v27.0) (SPSS Inc., Chicago, IL, USA). A p-value &lt; 0.05 was considered statistically significant.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<p>Our analysis included twenty-one treatment-na&#xef;ve NMOSD patients (NT group) and twenty-one matched healthy controls (HCs). Fecal samples were collected from all participants at enrollment. Of these 21 NMOSD patients, 5 subsequently underwent re-examination of fecal samples six months after immunosuppressant treatment (IST group). <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref> summarizes the detailed demographic information and clinical characteristics of the participants. The clinical manifestations and severity of each NMOSD patient are shown in <xref ref-type="supplementary-material" rid="SF1"><bold>Supplementary Table S1</bold></xref>.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Demographic and clinical characteristics of patients with NMOSD and healthy controls.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Characteristics</th>
<th valign="middle" align="center">NMOSD</th>
<th valign="middle" align="center">HCs</th>
<th valign="middle" align="center">P values</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">Stool samples, n</td>
<td valign="middle" align="left">21</td>
<td valign="middle" align="left">21</td>
<td valign="middle" align="left"/>
</tr>
<tr>
<td valign="middle" align="left">Female, n (%)</td>
<td valign="middle" align="left">19 (90.5)</td>
<td valign="middle" align="left">19 (90.5)</td>
<td valign="middle" align="left">1.00</td>
</tr>
<tr>
<td valign="middle" align="left">Age, years, mean &#xb1; SD</td>
<td valign="middle" align="left">45.1&#xa0;&#xb1;&#xa0;13.7</td>
<td valign="middle" align="left">47.1&#xa0;&#xb1;&#xa0;14.2</td>
<td valign="middle" align="left">0.27</td>
</tr>
<tr>
<td valign="middle" align="left">BMI, kg/m2, mean &#xb1; SD</td>
<td valign="middle" align="left">22.5&#xa0;&#xb1;&#xa0;3.6</td>
<td valign="middle" align="left">22.97&#xa0;&#xb1;&#xa0;2.57</td>
<td valign="middle" align="left">0.57</td>
</tr>
<tr>
<td valign="middle" align="left">Age at onset, years, mean &#xb1; SD</td>
<td valign="middle" align="left">37.3&#xa0;&#xb1;&#xa0;13.8</td>
<td valign="middle" align="center">&#x2013;</td>
<td valign="middle" align="center">&#x2013;</td>
</tr>
<tr>
<td valign="middle" align="left">Disease durations, years, mean &#xb1; SD</td>
<td valign="middle" align="left">8.5&#xa0;&#xb1;&#xa0;6.7</td>
<td valign="middle" align="center">&#x2013;</td>
<td valign="middle" align="center">&#x2013;</td>
</tr>
<tr>
<td valign="middle" align="left">Serum AQP4-IgG, n (%)</td>
<td valign="middle" align="left">21 (100)</td>
<td valign="middle" align="center">&#x2013;</td>
<td valign="middle" align="center">&#x2013;</td>
</tr>
<tr>
<td valign="middle" align="left">EDSS scores, median, range</td>
<td valign="middle" align="left">3.0 (0-7.0)</td>
<td valign="middle" align="center">&#x2013;</td>
<td valign="middle" align="center">&#x2013;</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>AQP4, aquaporin 4; BMI, body mass index; EDSS, expanded disability status scale; HCs, healthy controls; NMOSD, neuromyelitis optica spectrum disorder; SD, standard deviation.</p></fn>
</table-wrap-foot>
</table-wrap>
<sec id="s3_1">
<title>Alpha diversity</title>
<p>The &#x3b1;-diversity was compared to assess differences in microbial community structure between treatment-na&#xef;ve NMOSD patients and HCs. All &#x3b1;-diversity metrics (Chao1, observed OTUs, Shannon index, and Simpson index) were significantly higher in NMOSD patients than in HCs (P &lt; 0.001), indicating greater gut microbial richness and diversity in treatment-na&#xef;ve NMOSD patients. Furthermore, MMF-treated NMOSD patients showed reduced &#x3b1;-diversity compared to both treatment-na&#xef;ve patients and HCs, as demonstrated by the Chao1, observed species, Shannon index, and Simpson index (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2</bold></xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Alpha diversity of gut microbiota in NMOSD patients and HCs. Treatment-na&#xef;ve NMOSD patients demonstrated significantly higher gut microbial richness and diversity compared to HCs. Following IST, microbial diversity was significantly reduced in NMOSD patients relative to both pretreatment levels and HCs. HC, healthy control; IST, immunosuppressive therapy; NMOSD, neuromyelitis optica spectrum disorder; NT, no treatment.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g002.tif">
<alt-text content-type="machine-generated">Set of four box plots displays diversity indices for three groups: healthy controls (HC), NMOSD non-treated (NT), and NMOSD immunosuppressive therapy (IST). Top left shows Shannon index, top right Simpson index, bottom left Chao1, and bottom right observed OTUs. NMOSD NT group exhibits highest diversity in all indexes.</alt-text>
</graphic></fig>
</sec>
<sec id="s3_2">
<title>Beta diversity</title>
<p>To compare fecal microbial community structures between groups, &#x3b2;-diversity was assessed by calculating pairwise Bray-Curtis distances. PCoA based on the resulting distance matrices revealed distinct clustering between treatment-na&#xef;ve NMOSD patients and HCs, with clear separation along the primary axes (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3</bold></xref>). PERMANOVA confirmed significant compositional differences between groups (P &lt; 0.001, R&#xb2; = 0.13). LDA further identified phylogenetic clustering patterns. Specific genera (<italic>Streptococcus</italic> and <italic>Ruminococcus</italic>) showed significantly higher abundance in NMOSD patients compared to HCs, while <italic>Faecalibacterium</italic>, <italic>Pseudomonas</italic>, and <italic>Ralstonia</italic> were more abundant in HCs (LDA score &gt; 4, P &lt; 0.05; <xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>). Collectively, intestinal microbiota profiling demonstrated gut dysbiosis in NMOSD patients, characterized by increased &#x3b1;-diversity and altered community composition. These results establish significant differences between the gut microbiota of NMOSD patients and HCs.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>&#x3b2;-Diversity analysis of gut microbiota between treatment-na&#xef;ve NMOSD patients and HCs. Principal coordinate analysis (PCoA) based on unweighted UniFrac distances revealed distinct clustering patterns between NMOSD patients and HCs. Each point represents an individual sample (red circles = HCs; blue circles = NMOSD), with inter-point distances reflecting microbial community dissimilarity. HC, healthy control; NMOSD, neuromyelitis optica spectrum disorder.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g003.tif">
<alt-text content-type="machine-generated">Scatter plot showing two distinct clusters: blue dots representing NMOSD at lower left and red dots representing HC at upper right. Axes are labeled Component 1 (12.50 percent) and Component 2 (8.70 percent).</alt-text>
</graphic></fig>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Identification of differentially abundant gut microbes using LEfSe analysis. <bold>(A)</bold> Cladogram showing phylogenetic distribution of gut microbiota in treatment-na&#xef;ve NMOSD patients versus HCs. <bold>(B)</bold> LDA scores (log10 &gt; 4, p &lt; 0.05) identifying significant bacterial differences between treatment-na&#xef;ve NMOSD patients and HCs. <bold>(C)</bold> Cladogram showing phylogenetic distribution of gut microbiota in immunosuppressant-treated NMOSD patients versus HCs. <bold>(D)</bold> LDA scores (log10 &gt; 2, p &lt; 0.05) identifying significant bacterial differences between immunosuppressant-treated NMOSD patients and HCs. HC, healthy control; IST, immunosuppressive therapy; LDA, linear discriminant analysis; LEFSe, linear discriminant analysis effect size; NMOSD, neuromyelitis optica spectrum disorder; NT, no treatment.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g004.tif">
<alt-text content-type="machine-generated">Four-panel scientific figure depicts gut microbiota comparisons between NMOSD subgroups and healthy controls using cladograms and bar plots. Panel A shows a circular cladogram with colored branches highlighting taxa differences between NMOSD NT (red) and healthy controls (green), with an adjacent legend. Panel B presents a horizontal bar chart displaying linear discriminant analysis (LDA) scores for bacterial taxa enriched in each group. Panel C features another cladogram comparing NMOSD IST (green) and healthy controls (red). Panel D has a horizontal bar plot of LDA scores for two genera, Enterobacter and Klebsiella, distinguishing NMOSD IST and controls.</alt-text>
</graphic></fig>
<p>A comparative analysis of IST groups revealed only two significant differences in taxonomic biomarkers at the genus level between immunosuppressant-treated NMOSD patients and HCs. Specifically, <italic>Enterobacter</italic> was enriched in the NMOSD IST group (LDA score = 2.70, P &lt; 0.05), whereas <italic>Klebsiella</italic> showed high abundance in HCs (LDA score = 2.66, P &lt; 0.05) (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4B</bold></xref>).</p>
</sec>
<sec id="s3_3">
<title>Community analysis</title>
<p>Community-level analysis of gut microbiota revealed that Bacteroidetes, Firmicutes, Actinobacteria, and Proteobacteria were the predominant phyla in both treatment-na&#xef;ve NMOSD patients and healthy controls (HCs), collectively comprising &gt; 90% of the intestinal microbiota (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>). At the genus level, <italic>Bacteroides</italic>, <italic>Faecherlibacterium</italic>, <italic>Prevotella</italic>, <italic>Roseburia</italic>, and <italic>Blautia</italic> represented the five most abundant taxa (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Gut microbiota community structure comparison between NMOSD patients and HCs. <bold>(A)</bold> Phylum-level and <bold>(B)</bold> genus-level composition of gut microbial communities in treatment-na&#xef;ve NMOSD patients versus HCs. HC, healthy control; NMOSD, neuromyelitis optica spectrum disorder.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g005.tif">
<alt-text content-type="machine-generated">Two stacked bar charts display microbial diversity comparisons between NMOSD and healthy controls (HC). Chart A shows phyla-level relative frequencies, while chart B details genus-level relative abundances using distinct colors. Legends identify each taxon.</alt-text>
</graphic></fig>
</sec>
<sec id="s3_4">
<title>Overall metabolomics analysis</title>
<p>The metabolomes of 21 fecal samples from NMOSD patients and HCs were characterized and compared. Features with VIP scores &gt; 1.0 in multivariate analysis and p value &lt; 0.05 in univariate analysis were identified as the most important metabolites and were visualized in a heatmap (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6A</bold></xref>) and volcano plot (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6B</bold></xref>). A total of 19 key metabolites involved in amino acid, lipid, purine, and vitamin metabolism were altered in the feces of NMOSD patients compared with HCs (<xref ref-type="table" rid="T2"><bold>Table&#xa0;2</bold></xref>). The levels of D-Proline, L-Tyrosine, L-Methionine, L-Glutamate, N&amp;omega-Acetylhistamine, Val-Ala, L-Leucyl-L-Alanine, Val-Leu, Thr-Ile, Leu-Ile, Phe-Ile, Adenosine, Adenine, Hypoxanthine, Linoelaidic Acid, Palmitoyl Ethanolamide, 5(S)-HETE, and vitamin B1 were elevated in NMOSD fecal samples. In contrast, DL-phenylalanine levels were decreased in NMOSD (<xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7A</bold></xref>).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Metabolic profiles in treatment-na&#xef;ve NMOSD patients versus HCs. <bold>(A)</bold> Heatmap and <bold>(B)</bold> volcano plot displaying fecal metabolite profiles. Rows represent individual metabolites, columns represent study participants. Color gradients indicate relative metabolite abundance (red: increased; blue: decreased). HC, healthy control; NMOSD, neuromyelitis optica spectrum disorder.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g006.tif">
<alt-text content-type="machine-generated">Panel A shows a heatmap of gene expression data with hierarchical clustering; samples are grouped by two categories (NMOSD and HC) and values are color-coded from red (high) to blue (low). Panel B displays a volcano plot with each point representing a gene, showing log2 fold change on the x-axis and negative log10 p-value on the y-axis; significant upregulated genes are indicated in red, downregulated in green, non-significant in gray.</alt-text>
</graphic></fig>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Fecal identified differential metabolites between NMOSD patients and healthy controls.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Metabolite</th>
<th valign="middle" align="left">P value</th>
<th valign="middle" align="left">VIP</th>
<th valign="middle" align="left">FC</th>
<th valign="middle" align="left">Pathways</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">D-Proline</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.615</td>
<td valign="middle" align="left">0.027</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">L-Tyrosine</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.734</td>
<td valign="middle" align="left">5.068</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">L-Methionine</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.775</td>
<td valign="middle" align="left">3.099</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">L-Glutamate</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.459</td>
<td valign="middle" align="left">0.985</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">DL-Phenylalanine</td>
<td valign="middle" align="left">0.019</td>
<td valign="middle" align="left">1.351</td>
<td valign="middle" align="left">-0.487</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Val Ala</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.425</td>
<td valign="middle" align="left">3.998</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">L-Leucyl-L-Alanine</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.815</td>
<td valign="middle" align="left">5.301</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Val Leu</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.668</td>
<td valign="middle" align="left">0.462</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Thr Ile</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.377</td>
<td valign="middle" align="left">4.322</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Leu Ile</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.577</td>
<td valign="middle" align="left">4.455</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Phe Ile</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.639</td>
<td valign="middle" align="left">4.598</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Adenosine</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.694</td>
<td valign="middle" align="left">3.634</td>
<td valign="middle" align="left">Purine metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Adenine</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.378</td>
<td valign="middle" align="left">1.714</td>
<td valign="middle" align="left">Purine metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Hypoxanthine</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.440</td>
<td valign="middle" align="left">3.691</td>
<td valign="middle" align="left">Purine metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Linoelaidic Acid</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">1.651</td>
<td valign="middle" align="left">1.260</td>
<td valign="middle" align="left">Fatty acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Palmitoyl Ethanolamide</td>
<td valign="middle" align="left">&lt;0.001</td>
<td valign="middle" align="left">2.124</td>
<td valign="middle" align="left">3.666</td>
<td valign="middle" align="left">Fatty acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">5(S)-HETE</td>
<td valign="middle" align="left">0.004</td>
<td valign="middle" align="left">1.468</td>
<td valign="middle" align="left">2.322</td>
<td valign="middle" align="left">Fatty acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">N&amp;omega-Acetylhistamine</td>
<td valign="middle" align="left">0.017</td>
<td valign="middle" align="left">1.005</td>
<td valign="middle" align="left">0.819</td>
<td valign="middle" align="left">Other</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>FC, fold change; NMOSD, neuromyelitis optica spectrum disorder<bold>;</bold> VIP, variable importance in the projection.</p></fn>
</table-wrap-foot>
</table-wrap>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>PLS-DA VIP analysis of differential metabolites in treatment-na&#xef;ve NMOSD patients versus HCs. <bold>(A)</bold> VIP scores of significant metabolites (red: increased; green: decreased levels). <bold>(B)</bold> Enriched metabolic pathways showing pathway impact versus -log10(p-value) (top 10 pathways highlighted). HC, healthy control; NMOSD, neuromyelitis optica spectrum disorder; PLS-DA, partial least squares discriminant analysis; VIP, variable importance in projection.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g007.tif">
<alt-text content-type="machine-generated">Panel A shows a dot plot of VIP scores for features labeled by ID, with heatmap tiles indicating higher (red) or lower (green) abundance in NMOSD versus healthy controls; Panel B presents a KEGG enrichment bubble plot highlighting various pathways, where bubble size reflects feature count and color gradient from green to red indicates p.adjust values for statistical significance.</alt-text>
</graphic></fig>
</sec>
<sec id="s3_5">
<title>Metabolic pathway analysis</title>
<p>To identify biologically meaningful patterns based on the metabolomics data in feces, pathway profiling was conducted through the KEGG metabolic library using Metaboanalyst 3.0. Taken together, 44 gut microbe-related KEGG pathways were identified as differentially enriched between NMOSD group and HC group, among which 10 remained significant after FDR correction (FDR-adjusted P &lt; 0.05, <xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7B</bold></xref>). Elevated L-glutamate and L-tyrosine were enriched in alcoholism, central carbon metabolism in cancer, protein digestion and absorption, cocaine addiction, aminoacyl-tRNA biosynthesis, and amphetamine addiction. Meanwhile, L-tyrosine was also associated with Parkinson&#x2019;s disease and thiamine metabolism. Increased L-methionine was involved in central carbon metabolism in cancer, protein digestion and absorption, and aminoacyl-tRNA biosynthesis. Elevated adenosine was linked to alcoholism, Parkinson&#x2019;s disease, purine metabolism, and the cAMP signaling pathway. Higher hypoxanthine and thiamine levels were associated with purine metabolism and thiamine metabolism, respectively. Additionally, elevated (S)-lactate was enriched in central carbon metabolism in cancer and the cAMP signaling pathway, while increased adenine was linked to purine metabolism.</p>
</sec>
<sec id="s3_6">
<title>The impact of immunosuppressants on the gut microbiota of NMOSD patients</title>
<p>After six months of immunosuppressive therapy, fecal samples from treatment-na&#xef;ve NMOSD patients were re-collected and analyzed. Metabolomics analysis revealed that only two key metabolites were altered in the feces of NMOSD patients receiving immunosuppressants (IST group) compared with HCs (<xref ref-type="table" rid="T3"><bold>Table&#xa0;3</bold></xref>), as shown in a heatmap (<xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8A</bold></xref>) and volcano plot (<xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8B</bold></xref>). The levels of D-proline and oleic acid ethyl ester were elevated in the IST group (<xref ref-type="fig" rid="f9"><bold>Figure&#xa0;9A</bold></xref>). In total, 19 gut microbe-related KEGG pathways were differentially enriched between the IST and HC groups after FDR correction (P &lt; 0.05; <xref ref-type="fig" rid="f9"><bold>Figure&#xa0;9B</bold></xref>), including tyrosine metabolism, citrate cycle (TCA cycle), oxidative phosphorylation, alanine/aspartate/glutamate metabolism, pyruvate metabolism, butanoate metabolism, carbon fixation in prokaryotes, sulfur metabolism, cAMP signaling pathway, GABAergic synapse, glucagon signaling pathway, central carbon metabolism in cancer, propanoate metabolism, nicotinate/nicotinamide metabolism, two-component system, glyoxylate/dicarboxylate metabolism, phenylalanine metabolism, arginine/proline metabolism, and chlorocyclohexane/chlorobenzene degradation. D-proline was enriched in arginine/proline metabolism, while the remaining pathways were associated with increased succinate levels.</p>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Fecal identified differential metabolites between NMOSD IST group patients and healthy controls.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Metabolite</th>
<th valign="middle" align="left">P value</th>
<th valign="middle" align="left">VIP</th>
<th valign="middle" align="left">FC</th>
<th valign="middle" align="left">Pathways</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">D-Proline</td>
<td valign="middle" align="left">0.043</td>
<td valign="middle" align="left">2.360</td>
<td valign="middle" align="left">5.233</td>
<td valign="middle" align="left">Amino acid metabolism</td>
</tr>
<tr>
<td valign="middle" align="left">Oleic Acid ethyl ester</td>
<td valign="middle" align="left">0.049</td>
<td valign="middle" align="left">2.085</td>
<td valign="middle" align="left">4.619</td>
<td valign="middle" align="left">Fatty acid metabolism</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>FC, fold change; IST, immunosuppressive therapy, NMOSD, neuromyelitis optica spectrum disorder; VIP, variable importance in the projection.</p></fn>
</table-wrap-foot>
</table-wrap>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Metabolite profiles in immunosuppressant-treated NMOSD patients versus HCs. <bold>(A)</bold> Heatmap and <bold>(B)</bold> volcano plot displaying fecal metabolite profiles. Rows correspond to individual metabolites; columns represent study participants. Color gradients indicate relative abundance (red: increased; blue: decreased). HC, healthy control; IST, immunosuppressive therapy; NMOSD, neuromyelitis optica spectrum disorder.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g008.tif">
<alt-text content-type="machine-generated">Panel A displays a heatmap with hierarchical clustering illustrating differences in metabolite expression between two groups, colored in red and blue. Panel B shows a volcano plot with log2 fold change, highlighting upregulated (red), downregulated (green), and non-significant (grey) features.</alt-text>
</graphic></fig>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>PLS-DA VIP analysis of differential metabolites in immunosuppressant-treated NMOSD patients versus HCs. <bold>(A)</bold> VIP scores of significant metabolites (red: increased; green: decreased levels). <bold>(B)</bold> Enriched metabolic pathways showing pathway impact versus -log10(p-value) (top 2 pathways highlighted). HC, healthy control; NMOSD, neuromyelitis optica spectrum disorder; PLS-DA, partial least squares discriminant analysis; VIP, variable importance in projection.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-17-1747643-g009.tif">
<alt-text content-type="machine-generated">Panel A shows a dot plot ranking features by VIP scores with a heatmap indicating relative abundance for healthy controls and NMOSD IST groups. Panel B displays a KEGG enrichment dot plot with pathways on the y-axis, enrichment fold on the x-axis, dot color representing adjusted p-value, and dot size indicating count.</alt-text>
</graphic></fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>In the present study, we identified a striking gut microbiota imbalance in treatment-na&#xef;ve NMOSD patients compared to HCs, along with dysregulated microbiota-associated metabolic pathways. Community analysis and &#x3b2;-diversity revealed significant differences in microbial composition between the two groups. The &#x3b1;-diversity of fecal microbiota in NMOSD patients was significantly higher than that in HCs, suggesting a potential link between gut microbiota and NMOSD pathogenesis. Our findings align with prior Chinese cohort studies (<xref ref-type="bibr" rid="B18">18</xref>, <xref ref-type="bibr" rid="B19">19</xref>) reporting altered fecal microbiome &#x3b1;-diversity in NMOSD patients versus HCs. Although higher microbiota diversity is generally beneficial, elevated diversity and richness have also been observed in diseases like depression (<xref ref-type="bibr" rid="B20">20</xref>) and autism (<xref ref-type="bibr" rid="B21">21</xref>).</p>
<p>In our study, an overgrowth of opportunistic pathogens (e.g., <italic>Streptococcus</italic> and <italic>Ruminococcus</italic>) characterizes gut dysbiosis in NMOSD patients. Prior studies in Chinese NMOSD cohorts also reported associations with <italic>Streptococcus</italic> (<xref ref-type="bibr" rid="B7">7</xref>, <xref ref-type="bibr" rid="B18">18</xref>, <xref ref-type="bibr" rid="B19">19</xref>). Notably, Gong et&#xa0;al. found <italic>Streptococcus</italic> abundance positively correlated with disease severity (<xref ref-type="bibr" rid="B18">18</xref>). The immunomodulatory mechanism of <italic>Streptococcus</italic> remains unclear; however, its overrepresentation correlates with reduced short-chain fatty acid (SCFA) levels in NMOSD patients, which may enhance CD4+ T cell activity and pro-inflammatory responses (<xref ref-type="bibr" rid="B18">18</xref>, <xref ref-type="bibr" rid="B22">22</xref>). In addition, increased <italic>Streptococcus</italic> abundance may promote inflammation by elevating pro-inflammatory cytokines (TNF-&#x3b1;, IL-6, and IFN-&#x3b3;), which can also induce Th1/Th17 cell differentiation in humans (<xref ref-type="bibr" rid="B20">20</xref>, <xref ref-type="bibr" rid="B23">23</xref>). Furthermore, <italic>Streptococcus</italic> species such as <italic>S. pyogenes</italic> were shown to reduce Treg frequency in tumor microenvironments and impair Treg function via APC-derived IL-12 (<xref ref-type="bibr" rid="B24">24</xref>). These findings suggest <italic>Streptococcus</italic> may play a pivotal role in NMOSD pathogenesis. This contrasts with prior reports implicating <italic>Clostridium perfringens</italic> through molecular mimicry (<xref ref-type="bibr" rid="B25">25</xref>). We propose that genetic and environmental differences between Caucasian and Asian populations may underlie these discrepant findings.</p>
<p>We observed increased <italic>Ruminococcus</italic> abundance in the gut microbiota of NMOSD patients, though its mechanistic role requires further investigation. Xie et&#xa0;al. similarly reported elevated <italic>Ruminococcaceae</italic> levels during NMOSD acute phases (<xref ref-type="bibr" rid="B19">19</xref>). <italic>Ruminococcus</italic> is an anaerobic bacterial genus with context-dependent roles in disease. While some species promote health through SCFA production, others may exacerbate disease when overabundant (<xref ref-type="bibr" rid="B26">26</xref>). This genus is enriched in multiple disorders, including Crohn&#x2019;s disease, inflammatory bowel disease (IBD), spondylarthritis, and asthma (<xref ref-type="bibr" rid="B27">27</xref>). Notably, IBD patients exhibit transient <italic>Ruminococcus</italic> blooms correlating with disease activity (<xref ref-type="bibr" rid="B28">28</xref>), and its abundance associates with pro-inflammatory cytokines (e.g., IL-6, TNF-&#x3b1;) and LPS levels (<xref ref-type="bibr" rid="B29">29</xref>). Together with prior studies, our findings suggest that <italic>Streptococcus</italic> and <italic>Ruminococcus</italic> overgrowth may contribute to NMOSD pathogenesis via inflammatory mechanisms.</p>
<p>Furthermore, our data revealed depletion of beneficial commensal microbes (e.g., <italic>Faecalibacterium</italic>) in NMOSD patients. As a dominant bacterial genus in the healthy human gut, <italic>Faecalibacterium</italic> suppresses inflammation by producing SCFAs and IL-10. It also enhances gut barrier integrity, inhibiting pathogen invasion (<xref ref-type="bibr" rid="B30">30</xref>). The roles of <italic>Pseudomonas</italic> and <italic>Ralstonia</italic> in NMOSD remain unclear and warrant investigation. These findings suggest that gut microbiota dysbiosis, characterized by diminished commensal bacteria and increased opportunistic pathogens, may contribute to NMOSD pathogenesis through dysregulation of inflammatory mediators. However, further studies are needed to establish whether this dysbiosis initiates NMOSD or arises secondary to the disease.</p>
<p>To date, little is known about the potential link between gut microbiome metabolic pathways and NMOSD pathogenesis. Functional analyses revealed that gut microbiota alterations may contribute to NMOSD development through associated metabolic pathways. In this study, metabolomic profiling of fecal samples identified 19 significantly altered metabolites in treatment-na&#xef;ve NMOSD patients. Furthermore, PICRUSt analysis uncovered 10 dysregulated microbiome-associated metabolic pathways related to NMOSD. Notably, elevated fecal levels of glucogenic amino acids (including D-proline, L-glutamate, and L-methionine) and glucogenic/ketogenic amino acids (particularly L-tyrosine) were observed in NMOSD patients. These findings suggest potential disruptions in glucose and energy metabolism, as these amino acids may serve as alternative energy substrates. Glutamate metabolism has been previously reported to modulate both resting and activated T-cell function (<xref ref-type="bibr" rid="B31">31</xref>). Proline, a key component of mucus glycoproteins in colonic epithelial cells (<xref ref-type="bibr" rid="B32">32</xref>), showed elevated levels in stool samples, potentially indicating alterations in mucin synthesis and intestinal barrier function. These findings suggest a potential link between amino acid metabolism and gut homeostasis. We hypothesize that impaired amino acid absorption may lead to their accumulation in the intestinal lumen. Most notably, the significantly enriched aminoacyl-tRNA biosynthesis pathway, involving L-methionine, L-tyrosine and D-proline, indicates dysregulated amino acid turnover and protein biosynthesis in NMOSD patients.</p>
<p>Our analysis revealed skewed vitamin B metabolism in fecal samples from NMOSD patients. Costliow et&#xa0;al. suggested that under conditions of thiamine limitation, the biosynthesis of thiamine is crucial for the growth and competitive fitness of <italic>B. thetaiotaomicron</italic> (<xref ref-type="bibr" rid="B33">33</xref>), a <italic>Bacteroides</italic> species (phylum Bacteroidetes) in the human gut microbiota. We hypothesize that elevated fecal thiamine levels in NMOSD patients may sustain microbial survival by altering gut vitamin B metabolism. However, the immunomodulatory role of vitamin B remains poorly understood. Investigating vitamin B-mediated immune regulation could provide insights into NMOSD pathogenesis and reveal novel therapeutic targets.</p>
<p>Metabolites associated with purine metabolism were dysregulated in NMOSD patients compared to HCs, including hypoxanthine, adenine, and adenosine. Adenosine is a key immunosuppressive mediator secreted by Treg cells, which suppresses self-reactive immune responses, promotes transplant tolerance, and mitigates autoimmune disorders (<xref ref-type="bibr" rid="B34">34</xref>&#x2013;<xref ref-type="bibr" rid="B36">36</xref>). Additionally, lipid metabolism, which regulates critical cellular processes such as proliferation, differentiation, inflammation, and apoptosis, was significantly altered in NMOSD patients. Key perturbed metabolites in this pathway included linoelaidic acid, palmitoyl ethanolamide, and 5(S)-HETE. These findings suggest a potential link between lipid carbon-chain metabolism and NMOSD pathology. Thus, our study provides mechanistic insights into how gut microbiota-derived metabolic changes may contribute to NMOSD pathogenesis.</p>
<p>Notably, metabolism can be influenced by drug intake (<xref ref-type="bibr" rid="B37">37</xref>). Intriguingly, in the present study, not only did the gut microbiome composition in treatment-na&#xef;ve NMOSD patients closely resemble that of HCs after immunosuppressive therapy, but principal component analysis also revealed significant differences in overall metabolic profiles between NMOSD patients and HCs. Following therapy, the gut microbial metabolites in these patients became nearly identical to those in HCs. This suggests that metabolic alterations may arise not only from the disease itself but also from therapeutic interventions. A limitation of this study is the small cohort of treatment-na&#xef;ve NMOSD patients, which limits the statistical power to precisely evaluate medication-induced changes in fecal metabolites. Future studies should enroll a larger treatment-na&#xef;ve cohort to better delineate treatment-specific metabolic shifts.</p>
<p>Nevertheless, our study has several limitations. First, as a single-center case-control study, the sample size was limited, particularly in the longitudinal immunosuppressive treatment group. This precluded robust subgroup analyses to assess how disease status or clinical presentation specifically affects microbiota abundance. Second, although early longitudinal data from five patients suggested clinical benefits from immunosuppressants, the scarcity of serial samples prevented a comparative analysis of microbial or metabolomic dynamics between treatment responders and non-responders. Third, while cross-sectional analyses indicate that therapy is associated with a microbiota shift towards a state resembling healthy controls, the most direct evidence would come from paired pre- and post-treatment samples. The limited availability of such longitudinal specimens restricted our ability to perform this definitive within-patient analysis. Therefore, our findings necessitate further validation in larger, multi-center prospective studies with comprehensive longitudinal sampling. In addition, dietary habits may influence microbial community dynamics; thus, standardized dietary questionnaires need to be incorporated into gut microbiota research. Finally, the follow-up period of this longitudinal cohort study was limited to 6 months, which may not fully capture the long-term effects of immunosuppressive therapy on gut microbiota in NMOSD patients.</p>
</sec>
<sec id="s5" sec-type="conclusions">
<title>Conclusions</title>
<p>In conclusion, our research reveals that intestinal dysbiosis and aberrant metabolic pathways exist in NMOSD patients. Immunotherapy may shift the gut microbiota of NMOSD patients toward that of HCs. These findings could not only help elucidate the underlying mechanism by which microbiome dysbiosis contributes to NMOSD pathogenesis but also provide insights into potential individualized treatments for the disease. Further studies with larger sample sizes and extended follow-up durations are warranted to confirm our findings.</p>
</sec>
</body>
<back>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p></sec>
<sec id="s7" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>The studies involving humans were approved by The Medical Ethics Committee of West China Hospital, Sichuan University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants&#x2019; legal guardians/next of kin.</p></sec>
<sec id="s8" sec-type="author-contributions">
<title>Author contributions</title>
<p>QD: Writing &#x2013; original draft. XW: Writing &#x2013; review &amp; editing. ZS: Writing &#x2013; review &amp; editing. HC: Writing &#x2013; review &amp; editing. YZ: Writing &#x2013; review &amp; editing. RW: Writing &#x2013; review &amp; editing. ZM: Writing &#x2013; review &amp; editing. LK: Writing &#x2013; review &amp; editing. HZ: Funding acquisition, Writing &#x2013; review &amp; editing.</p></sec>
<sec id="s10" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p></sec>
<sec id="s11" sec-type="ai-statement">
<title>Generative AI statement</title>
<p>The author(s) declared that generative AI was not used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p></sec>
<sec id="s12" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p></sec>
<sec id="s13" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fimmu.2026.1747643/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fimmu.2026.1747643/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table1.docx" id="SF1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"><label>Supplementary Table&#xa0;1</label>
<caption>
<p>Clinical manifestations and severity of each patient with NMOSD. AQP4, aquaporin 4; APS, Area postrema syndrome; BS, brainstem syndrome; EDSS, expanded disability status scale; NMOSD, neuromyelitis optica spectrum disorder; ON, optic neuritis; TM, transverse myelitis.</p>
</caption></supplementary-material></sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Wingerchuk</surname> <given-names>DM</given-names></name>
<name><surname>Banwell</surname> <given-names>B</given-names></name>
<name><surname>Bennett</surname> <given-names>JL</given-names></name>
<name><surname>Cabre</surname> <given-names>P</given-names></name>
<name><surname>Carroll</surname> <given-names>W</given-names></name>
<name><surname>Chitnis</surname> <given-names>T</given-names></name>
<etal/>
</person-group>. 
<article-title>International consensus diagnostic criteria for neuromyelitis optica spectrum disorders</article-title>. <source>Neurology</source>. (<year>2015</year>) <volume>85</volume>:<page-range>177&#x2013;89</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1212/WNL.0000000000001729</pub-id>, PMID: <pub-id pub-id-type="pmid">26092914</pub-id>
</mixed-citation>
</ref>
<ref id="B2">
<label>2</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Hubbard</surname> <given-names>JA</given-names></name>
<name><surname>Hsu</surname> <given-names>MS</given-names></name>
<name><surname>Seldin</surname> <given-names>MM</given-names></name>
<name><surname>Binder</surname> <given-names>DK</given-names></name>
</person-group>. 
<article-title>Expression of the astrocyte water channel aquaporin-4 in the mouse brain</article-title>. <source>ASN Neuro</source>. (<year>2015</year>) <volume>7</volume>:<elocation-id>1759091415605486</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1177/1759091415605486</pub-id>, PMID: <pub-id pub-id-type="pmid">26489685</pub-id>
</mixed-citation>
</ref>
<ref id="B3">
<label>3</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Richards</surname> <given-names>JL</given-names></name>
<name><surname>Yap</surname> <given-names>YA</given-names></name>
<name><surname>McLeod</surname> <given-names>KH</given-names></name>
<name><surname>Mackay</surname> <given-names>CR</given-names></name>
<name><surname>Mari&#xf1;o</surname> <given-names>E</given-names></name>
</person-group>. 
<article-title>Dietary metabolites and the gut microbiota: an alternative approach to control inflammatory and autoimmune diseases</article-title>. <source>Clin Trans Immunol</source>. (<year>2016</year>) <volume>5</volume>:<fpage>e82</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/cti.2016.29</pub-id>, PMID: <pub-id pub-id-type="pmid">27350881</pub-id>
</mixed-citation>
</ref>
<ref id="B4">
<label>4</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Ivanov</surname> <given-names>II</given-names></name>
<name><surname>Atarashi</surname> <given-names>K</given-names></name>
<name><surname>Manel</surname> <given-names>N</given-names></name>
<name><surname>Brodie</surname> <given-names>EL</given-names></name>
<name><surname>Shima</surname> <given-names>T</given-names></name>
<name><surname>Karaoz</surname> <given-names>U</given-names></name>
<etal/>
</person-group>. 
<article-title>Induction of intestinal Th17 cells by segmented filamentous bacteria</article-title>. <source>Cell</source>. (<year>2009</year>) <volume>139</volume>:<page-range>485&#x2013;98</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cell.2009.09.033</pub-id>, PMID: <pub-id pub-id-type="pmid">19836068</pub-id>
</mixed-citation>
</ref>
<ref id="B5">
<label>5</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Varrin-Doyer</surname> <given-names>M</given-names></name>
<name><surname>Spencer</surname> <given-names>CM</given-names></name>
<name><surname>Schulze-Topphoff</surname> <given-names>U</given-names></name>
<name><surname>Nelson</surname> <given-names>PA</given-names></name>
<name><surname>Stroud</surname> <given-names>RM</given-names></name>
<name><surname>Cree</surname> <given-names>BA</given-names></name>
<etal/>
</person-group>. 
<article-title>Aquaporin 4-specific T cells in neuromyelitis optica exhibit a Th17 bias and recognize Clostridium ABC transporter</article-title>. <source>Ann Neurol</source>. (<year>2012</year>) <volume>72</volume>:<fpage>53</fpage>&#x2013;<lpage>64</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/ana.23651</pub-id>, PMID: <pub-id pub-id-type="pmid">22807325</pub-id>
</mixed-citation>
</ref>
<ref id="B6">
<label>6</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Long</surname> <given-names>Y</given-names></name>
<name><surname>Gao</surname> <given-names>C</given-names></name>
<name><surname>Qiu</surname> <given-names>W</given-names></name>
<name><surname>Hu</surname> <given-names>X</given-names></name>
<name><surname>Shu</surname> <given-names>Y</given-names></name>
<name><surname>Peng</surname> <given-names>F</given-names></name>
<etal/>
</person-group>. 
<article-title>Helicobacter pylori infection in Neuromyelitis Optica and Multiple Sclerosis</article-title>. <source>Neuroimmunomodulation</source>. (<year>2013</year>) <volume>20</volume>:<page-range>107&#x2013;12</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1159/000345838</pub-id>, PMID: <pub-id pub-id-type="pmid">23295676</pub-id>
</mixed-citation>
</ref>
<ref id="B7">
<label>7</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Shi</surname> <given-names>Z</given-names></name>
<name><surname>Qiu</surname> <given-names>Y</given-names></name>
<name><surname>Wang</surname> <given-names>J</given-names></name>
<name><surname>Fang</surname> <given-names>Y</given-names></name>
<name><surname>Zhang</surname> <given-names>Y</given-names></name>
<name><surname>Chen</surname> <given-names>H</given-names></name>
<etal/>
</person-group>. 
<article-title>Dysbiosis of gut microbiota in patients with neuromyelitis optica spectrum disorders: A cross sectional study</article-title>. <source>J Neuroimmunol</source>. (<year>2020</year>) <volume>339</volume>:<fpage>577126</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jneuroim.2019.577126</pub-id>, PMID: <pub-id pub-id-type="pmid">31841737</pub-id>
</mixed-citation>
</ref>
<ref id="B8">
<label>8</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Jarius</surname> <given-names>S</given-names></name>
<name><surname>Probst</surname> <given-names>C</given-names></name>
<name><surname>Borowski</surname> <given-names>K</given-names></name>
<name><surname>Franciotta</surname> <given-names>D</given-names></name>
<name><surname>Wildemann</surname> <given-names>B</given-names></name>
<name><surname>Stoecker</surname> <given-names>W</given-names></name>
<etal/>
</person-group>. 
<article-title>Standardized method for the detection of antibodies to aquaporin-4 based on a highly sensitive immunofluorescence assay employing recombinant target antigen</article-title>. <source>J Neurol Sci</source>. (<year>2010</year>) <volume>291</volume>:<page-range>52&#x2013;6</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jns.2010.01.002</pub-id>, PMID: <pub-id pub-id-type="pmid">20117794</pub-id>
</mixed-citation>
</ref>
<ref id="B9">
<label>9</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Fadrosh</surname> <given-names>DW</given-names></name>
<name><surname>Ma</surname> <given-names>B</given-names></name>
<name><surname>Gajer</surname> <given-names>P</given-names></name>
<name><surname>Sengamalay</surname> <given-names>N</given-names></name>
<name><surname>Ott</surname> <given-names>S</given-names></name>
<name><surname>Brotman</surname> <given-names>RM</given-names></name>
<etal/>
</person-group>. 
<article-title>An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform</article-title>. <source>Microbiome</source>. (<year>2014</year>) <volume>2</volume>:<fpage>6</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/2049-2618-2-6</pub-id>, PMID: <pub-id pub-id-type="pmid">24558975</pub-id>
</mixed-citation>
</ref>
<ref id="B10">
<label>10</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Caporaso</surname> <given-names>JG</given-names></name>
<name><surname>Kuczynski</surname> <given-names>J</given-names></name>
<name><surname>Stombaugh</surname> <given-names>J</given-names></name>
<name><surname>Bittinger</surname> <given-names>K</given-names></name>
<name><surname>Bushman</surname> <given-names>FD</given-names></name>
<name><surname>Costello</surname> <given-names>EK</given-names></name>
<etal/>
</person-group>. 
<article-title>QIIME allows analysis of high-throughput community sequencing data</article-title>. <source>Nat Methods</source>. (<year>2010</year>) <volume>7</volume>:<page-range>335&#x2013;6</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nmeth.f.303</pub-id>, PMID: <pub-id pub-id-type="pmid">20383131</pub-id>
</mixed-citation>
</ref>
<ref id="B11">
<label>11</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Bolger</surname> <given-names>AM</given-names></name>
<name><surname>Lohse</surname> <given-names>M</given-names></name>
<name><surname>Usadel</surname> <given-names>B</given-names></name>
</person-group>. 
<article-title>Trimmomatic: a flexible trimmer for Illumina sequence data</article-title>. <source>Bioinf (Oxford England)</source>. (<year>2014</year>) <volume>30</volume>:<page-range>2114&#x2013;20</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btu170</pub-id>, PMID: <pub-id pub-id-type="pmid">24695404</pub-id>
</mixed-citation>
</ref>
<ref id="B12">
<label>12</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Dong</surname> <given-names>LN</given-names></name>
<name><surname>Wang</surname> <given-names>JP</given-names></name>
<name><surname>Liu</surname> <given-names>P</given-names></name>
<name><surname>Yang</surname> <given-names>YF</given-names></name>
<name><surname>Feng</surname> <given-names>J</given-names></name>
<name><surname>Han</surname> <given-names>Y</given-names></name>
</person-group>. 
<article-title>Faecal and mucosal microbiota in patients with functional gastrointestinal disorders: Correlation with toll-like receptor 2/toll-like receptor 4 expression</article-title>. <source>World J Gastroenterol</source>. (<year>2017</year>) <volume>23</volume>:<page-range>6665&#x2013;73</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.3748/wjg.v23.i36.6665</pub-id>, PMID: <pub-id pub-id-type="pmid">29085211</pub-id>
</mixed-citation>
</ref>
<ref id="B13">
<label>13</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Price</surname> <given-names>MN</given-names></name>
<name><surname>Dehal</surname> <given-names>PS</given-names></name>
<name><surname>Arkin</surname> <given-names>AP</given-names></name>
</person-group>. 
<article-title>FastTree 2&#x2013;approximately maximum-likelihood trees for large alignments</article-title>. <source>PloS One</source>. (<year>2010</year>) <volume>5</volume>:<fpage>e9490</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0009490</pub-id>, PMID: <pub-id pub-id-type="pmid">20224823</pub-id>
</mixed-citation>
</ref>
<ref id="B14">
<label>14</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Segata</surname> <given-names>N</given-names></name>
<name><surname>Izard</surname> <given-names>J</given-names></name>
<name><surname>Waldron</surname> <given-names>L</given-names></name>
<name><surname>Gevers</surname> <given-names>D</given-names></name>
<name><surname>Miropolsky</surname> <given-names>L</given-names></name>
<name><surname>Garrett</surname> <given-names>WS</given-names></name>
<etal/>
</person-group>. 
<article-title>Metagenomic biomarker discovery and explanation</article-title>. <source>Genome Biol</source>. (<year>2011</year>) <volume>12</volume>:<fpage>R60</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/gb-2011-12-6-r60</pub-id>, PMID: <pub-id pub-id-type="pmid">21702898</pub-id>
</mixed-citation>
</ref>
<ref id="B15">
<label>15</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Langille</surname> <given-names>MG</given-names></name>
<name><surname>Zaneveld</surname> <given-names>J</given-names></name>
<name><surname>Caporaso</surname> <given-names>JG</given-names></name>
<name><surname>McDonald</surname> <given-names>D</given-names></name>
<name><surname>Knights</surname> <given-names>D</given-names></name>
<name><surname>Reyes</surname> <given-names>JA</given-names></name>
<etal/>
</person-group>. 
<article-title>Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences</article-title>. <source>Nat Biotechnol</source>. (<year>2013</year>) <volume>31</volume>:<page-range>814&#x2013;21</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nbt.2676</pub-id>, PMID: <pub-id pub-id-type="pmid">23975157</pub-id>
</mixed-citation>
</ref>
<ref id="B16">
<label>16</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Tankou</surname> <given-names>SK</given-names></name>
<name><surname>Regev</surname> <given-names>K</given-names></name>
<name><surname>Healy</surname> <given-names>BC</given-names></name>
<name><surname>Tjon</surname> <given-names>E</given-names></name>
<name><surname>Laghi</surname> <given-names>L</given-names></name>
<name><surname>Cox</surname> <given-names>LM</given-names></name>
<etal/>
</person-group>. 
<article-title>A probiotic modulates the microbiome and immunity in multiple sclerosis</article-title>. <source>Ann Neurol</source>. (<year>2018</year>) <volume>83</volume>:<page-range>1147&#x2013;61</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/ana.25244</pub-id>, PMID: <pub-id pub-id-type="pmid">29679417</pub-id>
</mixed-citation>
</ref>
<ref id="B17">
<label>17</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Li</surname> <given-names>Y</given-names></name>
<name><surname>Wang</surname> <given-names>HF</given-names></name>
<name><surname>Li</surname> <given-names>X</given-names></name>
<name><surname>Li</surname> <given-names>HX</given-names></name>
<name><surname>Zhang</surname> <given-names>Q</given-names></name>
<name><surname>Zhou</surname> <given-names>HW</given-names></name>
<etal/>
</person-group>. 
<article-title>Disordered intestinal microbes are associated with the activity of Systemic Lupus Erythematosus</article-title>. <source>Clin Sci (London England: 1979)</source>. (<year>2019</year>) <volume>133</volume>:<page-range>821&#x2013;38</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1042/CS20180841</pub-id>, PMID: <pub-id pub-id-type="pmid">30872359</pub-id>
</mixed-citation>
</ref>
<ref id="B18">
<label>18</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Gong</surname> <given-names>J</given-names></name>
<name><surname>Qiu</surname> <given-names>W</given-names></name>
<name><surname>Zeng</surname> <given-names>Q</given-names></name>
<name><surname>Liu</surname> <given-names>X</given-names></name>
<name><surname>Sun</surname> <given-names>X</given-names></name>
<name><surname>Li</surname> <given-names>H</given-names></name>
<etal/>
</person-group>. 
<article-title>Lack of short-chain fatty acids and overgrowth of opportunistic pathogens define dysbiosis of neuromyelitis optica spectrum disorders: A Chinese pilot study</article-title>. <source>Mult Scler</source>. (<year>2019</year>) <volume>25</volume>:<page-range>1316&#x2013;25</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1177/1352458518790396</pub-id>, PMID: <pub-id pub-id-type="pmid">30113252</pub-id>
</mixed-citation>
</ref>
<ref id="B19">
<label>19</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Xie</surname> <given-names>Q</given-names></name>
<name><surname>Sun</surname> <given-names>J</given-names></name>
<name><surname>Sun</surname> <given-names>M</given-names></name>
<name><surname>Wang</surname> <given-names>Q</given-names></name>
<name><surname>Wang</surname> <given-names>M</given-names></name>
</person-group>. 
<article-title>Perturbed microbial ecology in neuromyelitis optica spectrum disorder: Evidence from the gut microbiome and fecal metabolome</article-title>. <source>Mult Scler Relat Disord</source>. (<year>2024</year>) <volume>92</volume>:<fpage>105936</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.msard.2024.105936</pub-id>, PMID: <pub-id pub-id-type="pmid">39418776</pub-id>
</mixed-citation>
</ref>
<ref id="B20">
<label>20</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Jiang</surname> <given-names>H</given-names></name>
<name><surname>Ling</surname> <given-names>Z</given-names></name>
<name><surname>Zhang</surname> <given-names>Y</given-names></name>
<name><surname>Mao</surname> <given-names>H</given-names></name>
<name><surname>Ma</surname> <given-names>Z</given-names></name>
<name><surname>Yin</surname> <given-names>Y</given-names></name>
<etal/>
</person-group>. 
<article-title>Altered fecal microbiota composition in patients with major depressive disorder</article-title>. <source>Brain Behav Immun</source>. (<year>2015</year>) <volume>48</volume>:<page-range>186&#x2013;94</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.bbi.2015.03.016</pub-id>, PMID: <pub-id pub-id-type="pmid">25882912</pub-id>
</mixed-citation>
</ref>
<ref id="B21">
<label>21</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Finegold</surname> <given-names>SM</given-names></name>
<name><surname>Dowd</surname> <given-names>SE</given-names></name>
<name><surname>Gontcharova</surname> <given-names>V</given-names></name>
<name><surname>Liu</surname> <given-names>C</given-names></name>
<name><surname>Henley</surname> <given-names>KE</given-names></name>
<name><surname>Wolcott</surname> <given-names>RD</given-names></name>
<etal/>
</person-group>. 
<article-title>Pyrosequencing study of fecal microflora of autistic and control children</article-title>. <source>Anaerobe</source>. (<year>2010</year>) <volume>16</volume>:<page-range>444&#x2013;53</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.anaerobe.2010.06.008</pub-id>, PMID: <pub-id pub-id-type="pmid">20603222</pub-id>
</mixed-citation>
</ref>
<ref id="B22">
<label>22</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Su</surname> <given-names>J</given-names></name>
<name><surname>Wang</surname> <given-names>Y</given-names></name>
<name><surname>Yan</surname> <given-names>M</given-names></name>
<name><surname>He</surname> <given-names>Z</given-names></name>
<name><surname>Zhou</surname> <given-names>Y</given-names></name>
<name><surname>Xu</surname> <given-names>J</given-names></name>
<etal/>
</person-group>. 
<article-title>The beneficial effects of Polygonatum sibiricum Red. superfine powder on metabolic hypertensive rats via gut-derived LPS/TLR4 pathway inhibition</article-title>. <source>Phytomed: Int J Phytother Phytopharmacol</source>. (<year>2022</year>) <volume>106</volume>:<fpage>154404</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.phymed.2022.154404</pub-id>, PMID: <pub-id pub-id-type="pmid">36075182</pub-id>
</mixed-citation>
</ref>
<ref id="B23">
<label>23</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Assa</surname> <given-names>A</given-names></name>
<name><surname>Butcher</surname> <given-names>J</given-names></name>
<name><surname>Li</surname> <given-names>J</given-names></name>
<name><surname>Elkadri</surname> <given-names>A</given-names></name>
<name><surname>Sherman</surname> <given-names>PM</given-names></name>
<name><surname>Muise</surname> <given-names>AM</given-names></name>
<etal/>
</person-group>. 
<article-title>Mucosa-associated ileal microbiota in new-onset pediatric crohn&#x2019;s disease</article-title>. <source>Inflammatory Bowel Dis</source>. (<year>2016</year>) <volume>22</volume>:<page-range>1533&#x2013;9</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1097/MIB.0000000000000776</pub-id>, PMID: <pub-id pub-id-type="pmid">27271491</pub-id>
</mixed-citation>
</ref>
<ref id="B24">
<label>24</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Hirayama</surname> <given-names>M</given-names></name>
<name><surname>Nishikawa</surname> <given-names>H</given-names></name>
<name><surname>Nagata</surname> <given-names>Y</given-names></name>
<name><surname>Tsuji</surname> <given-names>T</given-names></name>
<name><surname>Kato</surname> <given-names>T</given-names></name>
<name><surname>Kageyama</surname> <given-names>S</given-names></name>
<etal/>
</person-group>. 
<article-title>Overcoming regulatory T-cell suppression by a lyophilized preparation of Streptococcus pyogenes</article-title>. <source>Eur J Immunol</source>. (<year>2013</year>) <volume>43</volume>:<fpage>989</fpage>&#x2013;<lpage>1000</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/eji.201242800</pub-id>, PMID: <pub-id pub-id-type="pmid">23436617</pub-id>
</mixed-citation>
</ref>
<ref id="B25">
<label>25</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Zamvil</surname> <given-names>SS</given-names></name>
<name><surname>Spencer</surname> <given-names>CM</given-names></name>
<name><surname>Baranzini</surname> <given-names>SE</given-names></name>
<name><surname>Cree</surname> <given-names>BAC</given-names></name>
</person-group>. 
<article-title>The gut microbiome in neuromyelitis optica</article-title>. <source>Neurotherapeut: J Am Soc Exp Neurother</source>. (<year>2018</year>) <volume>15</volume>:<fpage>92</fpage>&#x2013;<lpage>101</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s13311-017-0594-z</pub-id>, PMID: <pub-id pub-id-type="pmid">29280091</pub-id>
</mixed-citation>
</ref>
<ref id="B26">
<label>26</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Crost</surname> <given-names>EH</given-names></name>
<name><surname>Coletto</surname> <given-names>E</given-names></name>
<name><surname>Bell</surname> <given-names>A</given-names></name>
<name><surname>Juge</surname> <given-names>N</given-names></name>
</person-group>. 
<article-title>Ruminococcus gnavus: friend or foe for human health</article-title>. <source>FEMS Microbiol Rev</source>. (<year>2023</year>) <volume>47</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/femsre/fuad014</pub-id>, PMID: <pub-id pub-id-type="pmid">37015876</pub-id>
</mixed-citation>
</ref>
<ref id="B27">
<label>27</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Henke</surname> <given-names>MT</given-names></name>
<name><surname>Kenny</surname> <given-names>DJ</given-names></name>
<name><surname>Cassilly</surname> <given-names>CD</given-names></name>
<name><surname>Vlamakis</surname> <given-names>H</given-names></name>
<name><surname>Xavier</surname> <given-names>RJ</given-names></name>
<name><surname>Clardy</surname> <given-names>J</given-names></name>
</person-group>. 
<article-title>Ruminococcus gnavus, a member of the human gut microbiome associated with Crohn&#x2019;s disease, produces an inflammatory polysaccharide</article-title>. <source>Proc Natl Acad Sci United States America</source>. (<year>2019</year>) <volume>116</volume>:<page-range>12672&#x2013;7</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1904099116</pub-id>, PMID: <pub-id pub-id-type="pmid">31182571</pub-id>
</mixed-citation>
</ref>
<ref id="B28">
<label>28</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Hall</surname> <given-names>AB</given-names></name>
<name><surname>Yassour</surname> <given-names>M</given-names></name>
<name><surname>Sauk</surname> <given-names>J</given-names></name>
<name><surname>Garner</surname> <given-names>A</given-names></name>
<name><surname>Jiang</surname> <given-names>X</given-names></name>
<name><surname>Arthur</surname> <given-names>T</given-names></name>
<etal/>
</person-group>. 
<article-title>A novel Ruminococcus gnavus clade enriched in inflammatory bowel disease patients</article-title>. <source>Genome Med</source>. (<year>2017</year>) <volume>9</volume>:<fpage>103</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s13073-017-0490-5</pub-id>, PMID: <pub-id pub-id-type="pmid">29183332</pub-id>
</mixed-citation>
</ref>
<ref id="B29">
<label>29</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Li</surname> <given-names>N</given-names></name>
<name><surname>Wang</surname> <given-names>H</given-names></name>
<name><surname>Pei</surname> <given-names>H</given-names></name>
<name><surname>Wu</surname> <given-names>Y</given-names></name>
<name><surname>Li</surname> <given-names>L</given-names></name>
<name><surname>Ren</surname> <given-names>Y</given-names></name>
<etal/>
</person-group>. 
<article-title>Genus_Ruminococcus and order_Burkholderiales affect osteoporosis by regulating the microbiota-gut-bone axis</article-title>. <source>Front Microbiol</source>. (<year>2024</year>) <volume>15</volume>:<fpage>1373013</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmicb.2024.1373013</pub-id>, PMID: <pub-id pub-id-type="pmid">38835486</pub-id>
</mixed-citation>
</ref>
<ref id="B30">
<label>30</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Miquel</surname> <given-names>S</given-names></name>
<name><surname>Mart&#xed;n</surname> <given-names>R</given-names></name>
<name><surname>Rossi</surname> <given-names>O</given-names></name>
<name><surname>Berm&#xfa;dez-Humar&#xe1;n</surname> <given-names>LG</given-names></name>
<name><surname>Chatel</surname> <given-names>JM</given-names></name>
<name><surname>Sokol</surname> <given-names>H</given-names></name>
<etal/>
</person-group>. 
<article-title>Faecalibacterium prausnitzii and human intestinal health</article-title>. <source>Curr Opin Microbiol</source>. (<year>2013</year>) <volume>16</volume>:<page-range>255&#x2013;61</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.mib.2013.06.003</pub-id>, PMID: <pub-id pub-id-type="pmid">23831042</pub-id>
</mixed-citation>
</ref>
<ref id="B31">
<label>31</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Yang</surname> <given-names>Z</given-names></name>
<name><surname>Matteson</surname> <given-names>EL</given-names></name>
<name><surname>Goronzy</surname> <given-names>JJ</given-names></name>
<name><surname>Weyand</surname> <given-names>CM</given-names></name>
</person-group>. 
<article-title>T-cell metabolism in autoimmune disease</article-title>. <source>Arthritis Res Ther</source>. (<year>2015</year>) <volume>17</volume>:<fpage>29</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s13075-015-0542-4</pub-id>, PMID: <pub-id pub-id-type="pmid">25890351</pub-id>
</mixed-citation>
</ref>
<ref id="B32">
<label>32</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Monle&#xf3;n</surname> <given-names>D</given-names></name>
<name><surname>Morales</surname> <given-names>JM</given-names></name>
<name><surname>Barrasa</surname> <given-names>A</given-names></name>
<name><surname>L&#xf3;pez</surname> <given-names>JA</given-names></name>
<name><surname>V&#xe1;zquez</surname> <given-names>C</given-names></name>
<name><surname>Celda</surname> <given-names>B</given-names></name>
</person-group>. 
<article-title>Metabolite profiling of fecal water extracts from human colorectal cancer</article-title>. <source>NMR Biomed</source>. (<year>2009</year>) <volume>22</volume>:<page-range>342&#x2013;8</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/nbm.1345</pub-id>, PMID: <pub-id pub-id-type="pmid">19006102</pub-id>
</mixed-citation>
</ref>
<ref id="B33">
<label>33</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Costliow</surname> <given-names>ZA</given-names></name>
<name><surname>Degnan</surname> <given-names>PH</given-names></name>
</person-group>. 
<article-title>Thiamine acquisition strategies impact metabolism and competition in the gut microbe bacteroides thetaiotaomicron</article-title>. <source>mSystems</source>. (<year>2017</year>) <volume>2</volume>:<elocation-id>e00116-17</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/mSystems.00116-17</pub-id>, PMID: <pub-id pub-id-type="pmid">28951891</pub-id>
</mixed-citation>
</ref>
<ref id="B34">
<label>34</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Fletcher</surname> <given-names>JM</given-names></name>
<name><surname>Lonergan</surname> <given-names>R</given-names></name>
<name><surname>Costelloe</surname> <given-names>L</given-names></name>
<name><surname>Kinsella</surname> <given-names>K</given-names></name>
<name><surname>Moran</surname> <given-names>B</given-names></name>
<name><surname>O'Farrelly</surname> <given-names>C</given-names></name>
<etal/>
</person-group>. 
<article-title>CD39+Foxp3+ regulatory T Cells suppress pathogenic Th17 cells and are impaired in multiple sclerosis</article-title>. <source>J Immunol (Baltimore Md: 1950)</source>. (<year>2009</year>) <volume>183</volume>:<page-range>7602&#x2013;10</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.4049/jimmunol.0901881</pub-id>, PMID: <pub-id pub-id-type="pmid">19917691</pub-id>
</mixed-citation>
</ref>
<ref id="B35">
<label>35</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Mandapathil</surname> <given-names>M</given-names></name>
<name><surname>Hilldorfer</surname> <given-names>B</given-names></name>
<name><surname>Szczepanski</surname> <given-names>MJ</given-names></name>
<name><surname>Czystowska</surname> <given-names>M</given-names></name>
<name><surname>Szajnik</surname> <given-names>M</given-names></name>
<name><surname>Ren</surname> <given-names>J</given-names></name>
<etal/>
</person-group>. 
<article-title>Generation and accumulation of immunosuppressive adenosine by human CD4+CD25highFOXP3+ regulatory T cells</article-title>. <source>J Biol Chem</source>. (<year>2010</year>) <volume>285</volume>:<page-range>7176&#x2013;86</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1074/jbc.M109.047423</pub-id>, PMID: <pub-id pub-id-type="pmid">19858205</pub-id>
</mixed-citation>
</ref>
<ref id="B36">
<label>36</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Sitkovsky</surname> <given-names>M</given-names></name>
<name><surname>Lukashev</surname> <given-names>D</given-names></name>
<name><surname>Deaglio</surname> <given-names>S</given-names></name>
<name><surname>Dwyer</surname> <given-names>K</given-names></name>
<name><surname>Robson</surname> <given-names>SC</given-names></name>
<name><surname>Ohta</surname> <given-names>A</given-names></name>
</person-group>. 
<article-title>Adenosine A2A receptor antagonists: blockade of adenosinergic effects and T regulatory cells</article-title>. <source>Br J Pharmacol</source>. (<year>2008</year>) <volume>153 Suppl 1</volume>:<page-range>S457&#x2013;464</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/bjp.2008.23</pub-id>, PMID: <pub-id pub-id-type="pmid">18311159</pub-id>
</mixed-citation>
</ref>
<ref id="B37">
<label>37</label>
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Johnson</surname> <given-names>CH</given-names></name>
<name><surname>Patterson</surname> <given-names>AD</given-names></name>
<name><surname>Idle</surname> <given-names>JR</given-names></name>
<name><surname>Gonzalez</surname> <given-names>FJ</given-names></name>
</person-group>. 
<article-title>Xenobiotic metabolomics: major impact on the metabolome</article-title>. <source>Annu Rev Pharmacol Toxicol</source>. (<year>2012</year>) <volume>52</volume>:<fpage>37</fpage>&#x2013;<lpage>56</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev-pharmtox-010611-134748</pub-id>, PMID: <pub-id pub-id-type="pmid">21819238</pub-id>
</mixed-citation>
</ref>
</ref-list>
<fn-group>
<fn id="n1" fn-type="custom" custom-type="edited-by">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2351410">Mahsa Ghajarzadeh</ext-link>, Johns Hopkins University, United States</p></fn>
<fn id="n2" fn-type="custom" custom-type="reviewed-by">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2803778">Mohammadali Nahayati</ext-link>, University of Cincinnati, United States</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1447620">Shaoying Tan</ext-link>, Hong Kong Polytechnic University, Hong Kong SAR, China</p></fn>
</fn-group>
</back>
</article>