- Research
- Open access
- Published:
Metabolic disorder and intestinal microflora dysbiosis in chronic inflammatory demyelinating polyradiculoneuropathy
Cell & Bioscience volume 13, Article number: 6 (2023)
Abstract
Objective
Chronic inflammatory demyelinating polyradiculoneuropathy (CIDP) is a rare acquired immune-mediated neuropathy. Although microbial infection is potentially a contributing factor, a causative link between CIDP and microbial infection remains unclear. There is also no definitive biomarker for CIDP diagnostics and therapies. The present study aimed to characterize the serum metabolic profile and gut microbiome structure in CIDP.
Methods
Targeted metabolomics profiling of serum, using liquid chromatography-mass spectrometry, and metagenomics sequencing of stool samples from a cohort of CIDP and non-CIDP subjects were performed to evaluate serum metabolic profiles and gut microbiome structure in CIDP subjects relative to healthy controls.
Results
Metabolome data revealed that the bile acids profile was perturbed in CIDP with bile acids and arachidonic acid enriched significantly in CIDP versus non-CIDP controls. Metagenome data revealed that opportunistic pathogens, such as Klebsiella pneumonia and Megamonas funiformis, and genes involved in bacterial infection were notably more abundant in CIDP subjects, while gut microbes related to biotransformation of secondary bile acids were abnormal in CIDP versus non-CIDP subjects. Correlation analysis revealed that changes in secondary bile acids were associated with altered gut microbes, including Bacteroides ovatus, Bacteroides caccae, and Ruminococcus gnavus.
Conclusion
Bile acids and arachidonic acid metabolism were disturbed in CIDP subjects and might be affected by the dysbiosis of gut microbial flora. These findings suggest that the combination of bile acids and arachidonic acid could be used as a CIDP biomarker and that modulation of gut microbiota might impact the clinical course of CIDP.
Introduction
Chronic inflammatory demyelinating polyradiculoneuropathy (CIDP) is a rare, acquired, and immune-mediated demyelinating neuropathy causing limb weakness and sensory deficits [1]. The cause of CIDP remains unknown, and diagnosis is challenging. The disease course is steadily progressive for more than 8 weeks and is characterized by strong heterogeneity in terms of clinical presentation, prognosis, and treatment response [2]. CIDP diagnosis is usually based on a progressive or relapsing course over 2 months, electrophysiological evidence of peripheral demyelination, as well as response to immune-modulating therapies [3,4,5]. Auto-antibodies against ganglioside proteins [6] and autoreactive T-cell responses against myelin antigens [7, 8] are reported to be involved in the immunopathogenesis of CIDP, indicating different autoimmune targets are likely to be relevant. However, none of the immunopathological findings are specific for CIDP, meaning that the results of the diagnostic tests should be carefully interpreted as misdiagnosis commonly occurs. In addition, monitoring disease activity in order to guide treatment management is also a difficult clinical problem, as nerve conduction studies and electromyography do not adequately reflect functional disability, especially in severe disease courses [9, 10]. Hence, there is an urgent need for providing valid biomarkers to diagnose CIDP.
Another type of immune-mediated demyelinating neuropathy, Guillain Barre syndrome (GBS), also named acute inflammatory demyelinating polyradiculoneuropathy, was reported to be related to microbial infection [11, 12]. The post-infectious nature of GBS is today well-recognized in the majority of presenting cases, and approximately 70% of GBS cases are preceded by an infectious respiratory or gastrointestinal illness [13]. Molecular mimicry between microbial and axolemmal surface molecules after Campylobacter jejuni infection is known to be the pathophysiological basis of neural involvement leading to GBS [14, 15]. In contrast to GBS, data on preceding infections for CIDP are not only lacking but also heterogeneous. Earlier studies indicated that up to 30% of CIDP cases are associated with preceding infection [16,17,18,19], while more recent studies suggested 10% to 20% [20, 21]. However, the causative link between CIDP and microbial infection remains unproven.
To address the above knowledge gap, an unbiased approach was taken to identify serum metabolites and gut microbial flora that were significantly changed in CIDP subjects relative to non-CIDP controls. To that end, untargeted metabolomics profiling of serum, using liquid chromatography-mass spectrometry (LC–MS), was performed to discover potential biomarkers for CIDP, and metagenomics analysis of gut microbial species and genes in stool samples was performed to discover the association between CIDP and gut microbial flora. Additionally, Pearson’s correlation analysis was performed to explore the association between the bile acids (BAs) profile and the microbiotic abundance in the gut. Our findings provide an improved understanding of perturbations of the metabolome-microbiome interface in CIDP, including identification of potential biomarkers to monitor ongoing CIDP activity and modulation of gut microbiota to prevent CIDP.
Methods
Study population and sample handling
A total of 64 subjects, including 31 CIDP patients and 33 age- and sex-matched healthy controls (Additional file 1: Table S1), were enrolled in this study from December 2020 to February 2022 at the Qilu Hospital of Shandong University. Diagnosis of CIDP was based on established guidelines for diagnosis and treatment of CIDP [3, 22]. Serum and stool samples were collected from subjects, immediately immersed in liquid nitrogen, and then frozen at − 80 °C. All samples were kept frozen until further analysis.
Metabolite extraction
Metabolites were extracted from serum essentially according to previously reported methods [23, 24]. In short, samples were extracted by adding precooled methanol and acetonitrile (2:1, v/v), and internal standards mixes 1 and 2 were added for quality control of sample preparation. After purified, the metabolites were prepared for LC-MS analysis.
LC–MS/MS analysis
The untargeted LC-MS/MS used a Waters 2D UPLC (Waters, USA) tandem Q Exactive high-resolution mass spectrometer (Thermo Fisher Scientific, USA) for separation and detection of metabolites. Chromatographic separation was performed on a Waters ACQUITY UPLC BEH C18 column (1.7 μm, 2.1 mm × 100 mm), with a column temperature of 45 °C. The mobile phase consisted of 0.1% formic acid and acetonitrile in the positive mode and 10 mM ammonium formate and acetonitrile in the negative mode. The mass spectrometric settings for positive or negative ionization modes were as follows: spray voltage, 3.8/ − 3.2 kV; sheath gas flow rate, 40 arbitrary units (arb); aux gas flow rate, 10 arb; aux gas heater temperature, 350 °C; capillary temperature, 320 °C. The full scan range was 70–1050 m/z with a resolution of 70,000, and the automatic gain control (AGC) target for MS acquisitions was set to 3e6 with a maximum ion injection time of 100 ms.
Principal component analysis (PCA) analysis of metabolites
A PCA model [25] was established for the comparative analysis of CIDP group and non-CIDP group to observe the distribution and separation trends between samples. The data were subjected to log2 conversion before the PCA model was established, with scaling by the Pareto scaling method.
Partial least squares-discriminant analysis (PLS-DA) analysis of metabolites
A PLS-DA [26] model between the comparative CIDP group and non-CIDP group was established after log2-log conversion of the metabolite data, using the Par scaling method and sevenfold cross validation. The PLS-DA model [26] was subjected to 200 response permutation tests. Variable importance in projection (VIP) was used to measure the influence intensity and interpretation ability of each metabolite, and variables with a VIP greater than 1 were considered to have a significant effect on the classification of sample categories.
DNA extraction and metagenomics sequencing and de novo assembly
Metagenomic extraction and metagenomics sequencing were performed at the BGI Co., Ltd. (Shenzhen, China). Stool metagenomic DNA extractions were performed using the MolPure® Stool DNA Kit (YESEN, Shanghai, China). Metagenomic libraries were sequenced on the DNBSEQ platform. High-quality short reads of each DNA sample were assembled by MEGAHIT [27].
Gene prediction and functional annotation
MetaGeneMark [28] was used for ab initio prediction of metagenomic genes. Functional annotations were made using BacMet [29] (20180311), KEGG [30] (v101) databases, COG [31] (20201125), and SwissProt [32] (release-2021_04).
Non-metric multidimensional scale (NMDS) analysis and analysis of similarities (ANOSIM)
NMDS was conducted as previously described [33], and ANOSIM was conducted by the “vegan” packages of R (v3.3.1) based on Bray-Curtis distance [34]. Species diversity analysis was conducted as previously described [35].
Statistical analyses
All statistical analyses were performed using GraphPad Prism 9.0 software. A two-tailed t-test was performed to compare the CIDP and non-CIDP groups. Correlations between the gut microbial abundance and BAs profile were estimated using Pearson’s correlation analysis.
Results
Broad metabolic shifts in CIDP Patients
To characterize the metabolomic profile in CIDP, each serum sample was analyzed by LC–MS/MS in non-targeted mode, and MS/MS was performed using sensitive, high-resolution MS to collect data from both positive and negative ionization modes to improve metabolite coverage. A total of 4739 different metabolites were detected in positive ion mode, among which 1736 metabolites were identified, and 1504 different metabolites were detected in negative ion mode, with 744 metabolites identified. To determine whether there were significant differences in metabolites between the CIDP and non-CIDP groups, a PCA model was established to observe the distribution and separation trend of the two groups of samples (Fig. 1A), and PLS-DA was also performed (Fig. 1B). Both the PCA and PLS-DA models showed that the major patterns of serum metabolites were largely separated for CIDP patients versus non-CIDP controls.
Further statistical analysis revealed that 597 metabolites identified in positive ion mode increased in CIDP compared with non-CIDP, and 481 metabolites decreased; additionally, 208 metabolites increased in CIDP, and 182 metabolites decreased in negative ion mode, indicating broad metabolic shifts in CIDP versus non-CIDP controls, as also demonstrated by the volcano plots (Fig. 1C).
Enrichments in bile acids and arachidonic acid in CIDP
Quantification and comparison of serum metabolites identified some with significant changes in CIDP patients (Additional file 1: Table S2 and S3). Analysis of the enriched serum metabolites revealed that BAs changed significantly in the CIDP group versus non-CIDP controls (Table 1). For primary BAs, the concentration of cholic acid (CA) was markedly increased in the CIDP group by 9.61-fold and 36.46-fold, respectively, in the negative and positive ion modes. Concentrations of taurocholic acid (TCA) and taurochenodeoxycholic acid 3-sulfate (TCDCS) increased by 38.30-fold and about eightfold, and only glycocholic acid (GCA) decreased in the CIDP group. Among secondary BAs, deoxycholic acid (DCA) and taurodeoxycholic acid (TDCA) were increased in the CIDP group, while tauroursodeoxycholic acid (TUDCA), glycodeoxycholic acid (GDCA), lithocholic acid taurine conjugate (TLCA), and taurolithocholic acid 3-sulfate (TLCS) were lower in the CIDP group than in the non-CIDP group, indicating that secondary BA metabolic pathways that associated with by microorganisms were abnormal in CIDP versus non-CIDP subjects.
We also analyzed the ratio of conjugated BAs to free BAs in the serum (Fig. 2A, B), and results showed that this ratio was lower in the CIDP group than in the non-CIDP group, indicating that host BA homeostasis was disrupted in CIDP. Since DCA is the microbial metabolic product of CA, and lithocholic acid (LCA) is the microbial metabolic product of chenodeoxycholic acid (CDCA), we calculated the relative CA-to-DCA ratio and relative CDCA-to-LCA ratio in the CIDP and non-CIDP groups to indirectly evaluate the effects of the gut microbiota. Results showed notable increases in the CA/DCA and CDCA/LCA ratios in the serum of the CIDP group (Fig. 2C, D). Our results indicate that CIDP markedly inhibited intestinal secondary BA formation.
Arachidonic acid (AA) and the mediator in AA metabolism (thromboxane b2, TXB2) are reported to be involved in different neurological disorders and cardiovascular diseases [36,37,38]. In this study, metabolomics data showed that AA increased by about twofold and TXB2 increased by about 42-fold in CIDP versus non-CIDP subjects and that the precursors of AA (linoleic acid and 8z, 11z, 14z-eicosatrienoic acid) also increased (Table 1). As some CIDP subjects developed cardiovascular disease (Additional file 1: Table S1), we further evaluated whether the increased AA and TXB2 were associated with cardiovascular disease. However, statistical analysis revealed no significant differences in the relative abundance of AA (p = 0.8995) and TXB2 (p = 0.4702) in CIDP patients with or without cardiovascular disease, whereas the relative abundance of AA (p < 0.0001 and p < 0.0001) and TXB2 (p = 0.0007 and p = 0.0003) differed significantly between the non-CIDP group and CIDP group regardless of the presence of cardiovascular disease (Fig. 3), indicating an association between CIDP and significantly elevated levels of AA and TXB2.
Metabolic pathways are perturbed in CIDP subjects
Metabolic pathway enrichment analysis was performed based on the KEGG database, and differentially abundant metabolites and their pathways identified in negative and positive ion modes were evaluated (Additional file 1: Table S2 and Table S3). Results showed that 20 metabolic pathways presented significant differences in CIDP compared with non-CIDP in the negative ion mode, and 14 metabolic pathways were perturbed in the positive ion mode (Fig. 4). The metabolites that differed in abundance between the CIDP and non-CIDP groups were mainly involved in primary BA biosynthesis, bile secretion, linoleic acid metabolism, and taurine and hypotaurine metabolism, with BA and AA appearing most frequently in the differentially enriched metabolic pathways (Fig. 4), indicating that BA and AA might be two key metabolites for CIDP.
Gut microbiome composition changes in CIDP
In this study, 28 qualified stool samples from 31 CIDP subjects and 33 qualified stool samples from non-CIDP subjects were obtained and tested. To determine whether the sample size of stool samples was reasonable, a rarefaction curve was drawn with sample (randomly sampled) size and taxonomic abundance of the samples. The rarefaction curve of CIDP and non-CIDP group tended to be steady (Additional file 1: Fig S1), indicating that the sampling quantities for the CIDP and non-CIDP groups were reasonable and sufficient. In addition, both ANOSIM and NMDS statistical analyses showed that the inter-group difference between CIDP and non-CIDP group was greater than the intra-group difference at the family level, genus level, and species level (Additional file 1: Fig S2), indicating that grouping was significant.
Microbiotic abundance at the family level showed that Enterobacteriaceae, Selenomonadaceae and Acidaminococcaceae was significantly higher in group CIDP than in group non-CIDP, while Bacteroidaceae, Veillonellaceae, Tannerellaceae and Rikenellaceae were decreased in group CIDP (Fig. 5A). At the genus level, the most enriched genera in group CIDP were Klebsiella (Enterobacteriaceae), Megamonas (Selenomonadaceae), and Phascolarctobacterium (Acidaminococcaceae), while Bacteroides (Bacteroidaceae) and Phocaeicola (Bacteroidaceae) had reduced abundance (Fig. 5B). At the species level, the most enriched species in group CIDP were Klebsiella pneumonia (Enterobacteriaceae), Escherichia coli (Enterobacteriaceae), Megamonas funiformis (Selenomonadaceae), and Phascolarctobacterium faecium (Acidaminococcaceae), while Phocaeicola dorei (Bacteroidaceae), Phocaeicola vulgatus (Bacteroidaceae), Bacteroides uniformis (Bacteroidaceae), and Bacteroides sp. A1C1 (Bacteroidaceae) decreased (Fig. 5C). Overall, the microbiotic abundance data revealed gut microflora dysbiosis in the CIDP group, since potential pathogens such as Klebsiella were enriched, while indigenous microorganisms such as Bacteroides decreased.
Further analysis based on the microbe abundance in individual samples showed that almost all CIDP samples (except CIDP 32) had 1–3 types of above enriched microorganisms with levels higher than 1% at family, genus, and species levels, including Klebsiella pneumoniae, Escherichia coli, Megamonas funiformis, and Phascolarctobacterium faecium at levels > 1% in around 13 of the samples (Additional file 1: Fig S3). KEGG annotation of the metagenome also revealed that genes classified as associated with bacterial infection and disease were more abundant in CIDP subjects (Additional file 1: Fig S4). As Klebsiella pneumonia, Escherichia coli, and Megamonas funiformis are well-known opportunistic pathogens, these data indicate that CIDP may be related to intestinal microbial infection.
Enriched genes involved in intestinal microbial infection in CIDP subjects
To further investigate whether there is the intestinal microbial pathogenic factor in CIDP, genes in metagenomes were annotated and analyzed. Consisting with the enrichment of opportunistic pathogens Klebsiella pneumonia, Escherichia coli, and Megamonas funiformis in CIDP, eighteen genes involved in bacterial infection, nine genes involved in bacterial invasion, three genes involved in quorum sensing were enriched in CIDP (Table 2), 70 genes involved in virulence factors were also enriched in CIDP (Additional file 1: Table S4). KEGG analysis revealed that pathways associated with infection such as bacterial invasion of epithelial cells, pathogenic bacteria infection, quorum sensing and bacterial secretion system, were enriched in CIDP groups (Additional file 1: Table S5). Previous study showed that type II secretion system [39], type III secretion system [40], type IV secretion system [41], type VI secretion system [42], and type VII secretion system [43, 44], played important roles during bacterial infection. In this study, other 259 genes involved in type II secretion system, type III secretion system, type IV secretion system and type VI secretion system were also enriched in CIDP (Additional file 1: Table S6). The enrichment of so many virulence factor-encoding genes and infection-related genes in CIDP further indicates that CIDP may be related to intestinal microbial infection.
Association between gut microbiotic structure and BA profiles
Secondary BAs are produced by gut microbes, such as Bacteroides, Clostridium, Lactobacillus, and Ruminococcus, via transformation of host-derived primary BAs [45]. In this study, the secondary BA profiles of CIDP subjects and non-CIDP subjects had notable differences (Table 1), indicating that the secondary BA metabolic pathways of microorganisms were abnormal in CIDP. Metagenome data revealed that gut microbes associated with the transformation of secondary BAs, including several genus of Bacteroides, Parabacteroides and Ruminococcus, were decreased in CIDP subjects versus non-CIDP subjects (Fig. 6A), which may result in abnormal transformation of secondary BAs. To further elucidate any correlation between these bacterial species and secondary BA transformation, Pearson’s correlation analyses were performed. Results revealed that the decreased TUDCA level in the serum of CIDP subjects significantly and positively correlated with Bacteroides sp. PHL_2737, Bacteroides sp. M10, and Bacteroides ovatus; the decreased GDCA level significantly and positively correlated with Bacteroides caccae and Bacteroides coprosuis; the decreased TLCA level significantly and positively correlated with Bacteroides caccae; and the increased TCDCS significantly and positively correlated with Ruminococcus gnavus (Fig. 6B, C). These results suggested that changes in the secondary BA profile might be associated with an altered composition of gut microbiota.
Discussion
Although Campylobacter jejuni infection is known to be the pathophysiological basis of neural involvement leading to GBS [14, 15], surprisingly, this pathogen was not detected in any of the stool samples from CIDP subjects in our study. However, our study revealed that gut microbiome structures differed significantly between the CIDP group and non-CIDP group at the family, genus, and species levels. Klebsiella pneumonia (Enterobacteriaceae), Escherichia coli (Enterobacteriaceae), Megamonas funiformis (Selenomonadaceae), and Phascolarctobacterium faecium (Acidaminococcaceae) were the most enriched species in group CIDP, and Wilcoxon rank-sum test analysis indicated that Klebsiella pneumoniae was the most significantly differing species between the CIDP and non-CIDP groups in terms of abundance. Many virulence factor-encoding genes and infection-related genes were enriched in CIDP subjects (Table 2, Additional file 1: Tables S5-S6). Besides, KEGG analysis of metagenome data also revealed that genes belonging to pathways associated with bacterial invasion of epithelial cells and bacterial infection were enriched in CIDP groups (Additional file 1: Table S5 and Figure S4). Klebsiella pneumonia and Escherichia coli are well-known opportunistic pathogens, and Megamonas funiformis may be associated with inflammatory bowel disease [46], colorectal cancer [47], multiple system atrophy [48], and myasthenia gravis [49]. Phascolarctobacterium faecium is potentially associated with the metabolic state and mood of the host [50]. These data support a link between the gut microbial infection and CIDP, although the function of these opportunistic pathogens in CIDP remains unclear.
Currently, there are no definitive biomarkers for CIDP. In this study, the serum BA profile was notably altered in CIDP subjects versus healthy subjects. BA metabolic disorders are associated with type 2 diabetes, type 1 diabetes mellitus, dyslipidemia, non-alcoholic fatty liver disease, and some neurodegenerative disorders [51,52,53,54,55]. Moreover, previous studies suggest that altered BA profiles are found with cirrhosis, schizophrenia, and Alzheimer's disease and that BAs may be candidate diagnostic biomarkers for these diseases [56,57,58]. In our study, host-derived primary BAs (CA, TCA, and TCDCS) were increased by over eightfold in the serum of CIDP subjects, whereas microbial-derived secondary BAs (GDCA, TUDCA, TLCA, and TLCS) were decreased. Notably, BA administration, particularly that of TUDCA and ursodeoxycholic acid (UDCA), contributed to neurologic symptom improvements in animal models of Alzheimer’s, Parkinson’s, and Huntington disease [59,60,61]. Here, TUDCA was also noticeably decreased in CIDP subjects; however, the relationship between TUDCA and CIDP remains unclear. Additionally, other differences were found in the serum BA profile of the CIDP group, as the conjugated serum BAs to free BAs ratio was decreased, while TCA/DCA and TCDCA/LCA ratios were increased. Interestingly, dehydrocholic acid (DHCA), which was considered as synthetic BA for clinical use previously, increased by about 47-fold in CIDP versus non-CIDP (Table 1), while CIDP patients were not treated with DHCA drug. These data suggested that the levels of some specific serum BAs may change with the development of CIDP and further suggested that dysmetabolism of BAs may be used as a surrogate biomarker to distinguish CIDP patients from the healthy population. However, their predictive value still needs further investigation and validation in larger prospective cohort studies.
Multiple studies have focused on the effect of intestinal microorganisms on secondary BA conversion as host-derived primary BAs are deconjugated via microbial bile salt hydrolases (BSH) and then transformed into secondary BAs through microbial hydroxysteroid dehydrogenases, which are widely found in the gut microbiota, including Bacteroides, Clostridium, Bifidobacterium, Eubacterium, Lactobacillus, Peptostreptococcus, Listeria, and Enterococcus [62,63,64,65]. However, the effects of gut microbiota in AA biosynthesis are unclear, although many environmental microorganisms have been reported to produce AA, including fungi of the genera Mortierella [66,67,68], Diasporangium [69], and Pythium [70], and bacteria such as Aureispira maritime [71]. In this study, AA metabolism was found abnormal in CIDP subjects, and the fecal metagenome data revealed that the number of eukaryotic microbes in the host gut was slightly increased in the CIDP group (0.035%) versus non-CIDP group (0.021%), although the differences were not statistically significant (Additional file 1: Fig S5). Additionally, KEGG level 3 analyses revealed increased gene abundance associated with steroid metabolism, unsaturated fatty acid biosynthesis, and linolenic acid and glycerolipid metabolism in the gut microbes of CIDP subjects (Additional file 1: Table S5), suggesting that intestinal microbes in CIDP patients might enhance the degradation of food into short chain fatty acids, including linolenic acid. Furthermore, metabolome data showed that the levels of serum linoleic acid and 8z, 11z, 14z-eicosatrienoic acid (an intermediate of AA synthesized from linoleic acid) were increased in the CIDP group (Table 1). In summary, we speculate that the dysbiosis of gut microbiota in CIDP subjects may result in the abnormal AA content in host serum as intestinal microbes could ferment food into short chain fatty acids, including linoleic acid, to provide raw materials for intestinal microbes or the host to synthesize AA (Fig. 7). However, how eukaryotic microorganisms in the gut affect AA synthesis in CIDP subjects remains to be further elucidated.
Microbes in healthy human guts belong to a small number of families, mainly dominated by Bacteroidetes and Firmicutes [72, 73]. Here, metagenomics data revealed that the gut microbial population changed significantly in CIDP subjects. For example, the abundance of Bacteroides (Bacteroidaceae, Bacteroidetes) was much lower in the stool samples of the CIDP subjects than non-CIDP subjects, while pathogenic bacteria such as Klebsiella pneumonia (Enterobacteriaceae) were enriched in the stool samples of the CIDP subjects. KEGG annotation of metagenomic data also revealed that genes involved in bacterial infection and disease were enriched in CIDP subjects. These findings indicate that the normal microbiota-host interactions are disturbed with CIDP.
The past few years have brought exciting new insights in the field of probiotic therapies. Guo et al. reviewed the potentials and challenges of using Clostridium species as probiotics to impact human heath through affecting microbe-host interactions [74]. Foley et al. investigated the role of BSHs encoded by Lactobacillus gasseri and Lactobacillus acidophilus, which were used as probiotics, and revealed that BA type and BSH substrate preferences affected the growth of these two species [75]. These reports contribute to the understanding of probiotic bacterial survival mechanisms in various BA-rich niches and provide information for future use of probiotics as therapeutic tools for manipulating gut microbiota. Given that our findings suggest abnormal microbiota-host interactions are present in CIDP patients, we propose using probiotics to improve these interactions in such patients through the gut-BA-host axis. Indeed, further development of probiotics efficient in secondary BA conversion could provide an augmenting therapy for CIDP.
Conclusion
This study revealed broad metabolic shifts in CIDP subjects versus control subjects, with perturbations in the serum BA profile and significantly increased AA in CIDP. Additionally, the gut microbiome composition was significantly changed in CIDP subjects, and levels of gut microbes related to secondary BA biotransformation were abnormal in CIDP. Our findings suggest that CIDP may be related to alterations in the gut microbiome, including potential infection with pathogens, and that the dysbiosis of gut microbial flora impacts BA and AA metabolism in CIDP subjects. Our findings also suggest that a combination of BAs and AA could be used as a surrogate biomarker of CIDP and that modulation of gut microbiota could impact the CIDP clinical course.
Availability of data and materials
Metagenome data were deposited to NCBI with SRA data: PRJNA878638. Metabolomics data have been deposited to the EMBL-EBI MetaboLights database with the identifier MTBLS5828.
References
Broers MC, Bunschoten C, Nieboer D, Lingsma HF, Jacobs BC. Incidence and prevalence of chronic inflammatory demyelinating polyradiculoneuropathy: a systematic review and meta-analysis. Neuroepidemiology. 2019;52(3–4):161–72.
Latov N. Diagnosis and treatment of chronic acquired demyelinating polyneuropathies. Nat Rev Neurol. 2014;10(8):435–46.
Lehmann HC, Burke D, Kuwabara S. Chronic inflammatory demyelinating polyneuropathy: update on diagnosis, immunopathogenesis and treatment. J Neurol Neurosurg Psychiatry. 2019;90(9):981–7.
Mathey EK, Park SB, Hughes RA, Pollard JD, Armati PJ, Barnett MH, et al. Chronic inflammatory demyelinating polyradiculoneuropathy: from pathology to phenotype. J Neurol Neurosurg Psychiatry. 2015;86(9):973–85.
Van den Bergh PY, Hadden RD, Bouche P, Cornblath DR, Hahn A, Illa I, et al. European Federation of Neurological Societies/Peripheral Nerve Society guideline on management of chronic inflammatory demyelinating polyradiculoneuropathy: report of a joint task force of the European Federation of Neurological Societies and the Peripheral Nerve Society - first revision. Eur J Neurol. 2010;17(3):356–63.
Klehmet J, Marschenz S, Ruprecht K, Wunderlich B, Buttner T, Hiemann R, et al. Analysis of anti-ganglioside antibodies by a line immunoassay in patients with chronic-inflammatory demyelinating polyneuropathies (CIDP). Clin Chem Lab Med. 2018;56(6):919–26.
Diederich JM, Staudt M, Meisel C, Hahn K, Meinl E, Meisel A, et al. Neurofascin and compact myelin antigen-specific t cell response pattern in chronic inflammatory demyelinating polyneuropathy subtypes. Front Neurol. 2018;9:171.
Staudt M, Diederich JM, Meisel C, Meisel A, Klehmet J. Differences in peripheral myelin antigen-specific T cell responses and T memory subsets in atypical versus typical CIDP. BMC Neurol. 2017;17(1):81.
Gruter T, Motte J, Fisse AL, Bulut Y, Kose N, Athanasopoulos D, et al. Pathological spontaneous activity as a prognostic marker in chronic inflammatory demyelinating polyneuropathy. Eur J Neurol. 2020;27(12):2595–603.
Kerasnoudis A, Pitarokoili K, Behrendt V, Gold R, Yoon MS. Correlation of nerve ultrasound, electrophysiological and clinical findings in chronic inflammatory demyelinating polyneuropathy. J Neuroimaging. 2015;25(2):207–16.
Hughes RA, Cornblath DR. Guillain-Barre syndrome. Lancet. 2005;366(9497):1653–66.
Sharma A, Lal V, Modi M, Vaishnavi C, Prabhakar S. Campylobacter jejuni infection in Guillain-Barre syndrome: a prospective case control study in a tertiary care hospital. Neurol India. 2011;59(5):717–21.
Jacobs BC, Rothbarth PH, van der Meche FG, Herbrink P, Schmitz PI, de Klerk MA, et al. The spectrum of antecedent infections in Guillain-Barre syndrome: a case-control study. Neurology. 1998;51(4):1110–5.
Willison HJ, Jacobs BC, van Doorn PA. Guillain-Barre syndrome. Lancet. 2016;388(10045):717–27.
Shahrizaila N, Lehmann HC, Kuwabara S. Guillain-Barre syndrome. Lancet. 2021;397(10280):1214–28.
Gorson KC, Allam G, Ropper AH. Chronic inflammatory demyelinating polyneuropathy: clinical features and response to treatment in 67 consecutive patients with and without a monoclonal gammopathy. Neurology. 1997;48(2):321–8.
McCombe PA, Pollard JD, McLeod JG. Chronic inflammatory demyelinating polyradiculoneuropathy. A clinical and electrophysiological study of 92 cases. Brain. 1987;110(Pt 6):1617–30.
Simmons Z, Albers JW, Bromberg MB, Feldman EL. Presentation and initial clinical course in patients with chronic inflammatory demyelinating polyradiculoneuropathy: comparison of patients without and with monoclonal gammopathy. Neurology. 1993;43(11):2202–9.
Simmons Z, Wald JJ, Albers JW. Chronic inflammatory demyelinating polyradiculoneuropathy in children: I. Presentation, electrodiagnostic studies, and initial clinical course, with comparison to adults. Muscle Nerve. 1997;20(8):1008–15.
Kuitwaard K, Bos-Eyssen ME, Blomkwist-Markens PH, van Doorn PA. Recurrences, vaccinations and long-term symptoms in GBS and CIDP. J Peripher Nerv Syst. 2009;14(4):310–5.
Rajabally YA, Peric S, Bozovic I, Loo LK, Kalac A, Palibrk A, et al. Antecedent infections and vaccinations in chronic inflammatory demyelinating polyneuropathy: a European collaborative study. Muscle Nerve. 2021;64(6):657–61.
Berger AR, Bradley WG, Brannagan TH, Busis NA, Cros DP, Dalakas MC, et al. Guidelines for the diagnosis and treatment of chronic inflammatory demyelinating polyneuropathy. J Peripher Nerv Syst. 2003;8(4):282–4.
Dunn WB, Broadhurst D, Begley P, Zelena E, Francis-McIntyre S, Anderson N, et al. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat Protoc. 2011;6(7):1060–83.
Sarafian MH, Gaudin M, Lewis MR, Martin FP, Holmes E, Nicholson JK, et al. Objective set of criteria for optimization of sample preparation procedures for ultra-high throughput untargeted blood plasma lipid profiling by ultra performance liquid chromatography-mass spectrometry. Anal Chem. 2014;86(12):5766–74.
Ben Salem K, Ben AA. Principal component analysis (PCA). Tunis Med. 2021;99(4):383–9.
Gromski PS, Muhamadali H, Ellis DI, Xu Y, Correa E, Turner ML, et al. A tutorial review: Metabolomics and partial least squares-discriminant analysis–a marriage of convenience or a shotgun wedding. Anal Chim Acta. 2015;879:10–23.
Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31(10):1674–6. https://doi.org/10.1093/bioinformatics/btv033.
Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010;38(12): e132.
Pal C, Bengtsson-Palme J, Rensing C, Kristiansson E, Larsson DG. BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 2014;42:D737-43.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Galperin MY, Makarova KS, Wolf YI, Koonin EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 2015;43:D261-9.
Poux S, Arighi CN, Magrane M, Bateman A, Wei CH, Lu Z, et al. On expert curation and scalability: UniProtKB/Swiss-Prot as a case study. Bioinformatics. 2017;33(21):3454–60.
Woods CT, Robertson S, Sinclair WH, Collier NF. Non-metric multidimensional performance indicator scaling reveals seasonal and team dissimilarity within the National Rugby League. J Sci Med Sport. 2018;21(4):410–5.
Wang ST, Meng XZ, Zhang JH, Dai YF, Shen Y, Xu XY, et al. 16S rRNA sequencing analysis of the correlation between the intestinal microbiota and body-mass of grass carp (Ctenopharyngodon idella). Comp Biochem Physiol Part D Genomics Proteomics. 2020;35: 100699.
Walters KE, Martiny JBH. Alpha-, beta-, and gamma-diversity of bacteria varies across habitats. PLoS ONE. 2020;15(9): e0233872.
Hoxha M. Arachidonic acid mediators and their role in neurological disease. CNS Neurol Disord Drug Targets. 2022;21(2):106–7.
Hoxha M, Rovati GE, Cavanillas AB. The leukotriene receptor antagonist montelukast and its possible role in the cardiovascular field. Eur J Clin Pharmacol. 2017;73(7):799–809.
Szczuko M, Koziol I, Kotlega D, Brodowski J, Drozd A. The role of thromboxane in the course and treatment of ischemic stroke: review. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms222111644.
von Tils D, Bladel I, Schmidt MA, Heusipp G. Type II secretion in Yersinia-a secretion system for pathogenicity and environmental fitness. Front Cell Infect Microbiol. 2012;2:160.
Hauser AR. The type III secretion system of Pseudomonas aeruginosa: infection by injection. Nat Rev Microbiol. 2009;7(9):654–65.
Skoog EC, Martin ME, Barrozo RM, Hansen LM, Cai LP, Lee SJ, et al. Maintenance of type IV secretion function during Helicobacter pylori infection in Mice. mBio. 2020. https://doi.org/10.1128/mBio.03147-20.
Loeven NA, Perault AI, Cotter PA, Hodges CA, Schwartzman JD, Hampton TH, et al. The Burkholderia cenocepacia Type VI secretion system effector TecA Is a virulence factor in mouse models of lung infection. mBio. 2021;12(5):e0209821.
Bowman L, Palmer T. The type VII secretion system of Staphylococcus. Annu Rev Microbiol. 2021;75:471–94.
Kim YS, Yang CS, Nguyen LT, Kim JK, Jin HS, Choe JH, et al. Mycobacterium abscessus ESX-3 plays an important role in host inflammatory and pathological responses during infection. Microbes Infect. 2017;19(1):5–17.
Winston JA, Theriot CM. Diversification of host bile acids by members of the gut microbiota. Gut Microbes. 2020;11(2):158–71.
Sakon H, Nagai F, Morotomi M, Tanaka R. Sutterella parvirubra sp. nov., and Megamonas funiformis sp. nov., isolated from human faeces. Int J Syst Evol Microbiol. 2008;58(4):970–5.
Ren X, Xu J, Zhang Y, Chen G, Zhang Y, Huang Q, et al. Bacterial alterations in post-cholecystectomy patients are associated with colorectal cancer. Front Oncol. 2020;10:1418.
Wan L, Zhou X, Wang C, Chen Z, Peng H, Hou X, et al. Alterations of the gut microbiota in multiple system atrophy patients. Front Neurosci. 2019;13:1102.
Liu P, Jiang Y, Gu S, Xue Y, Yang H, Li Y, et al. Metagenome-wide association study of gut microbiome revealed potential microbial marker set for diagnosis of pediatric myasthenia gravis. BMC Med. 2021;19(1):159.
Wu F, Guo X, Zhang J, Zhang M, Ou Z, Peng Y. Phascolarctobacterium faecium abundant colonization in human gastrointestinal tract. Exp Ther Med. 2017;14(4):3122–6.
Chavez-Talavera O, Tailleux A, Lefebvre P, Staels B. Bile acid control of metabolism and inflammation in obesity, type 2 diabetes Dyslipidemia, and Nonalcoholic fatty liver disease. Gastroenterology. 2017;152(7):1679–94.
Woodhams L, Al-Salami H. The roles of bile acids and applications of microencapsulation technology in treating type 1 diabetes mellitus. Ther Deliv. 2017;8(6):401–9.
Sipe LM, Chaib M, Pingili AK, Pierre JF, Makowski L. Microbiome, bile acids, and obesity: how microbially modified metabolites shape anti-tumor immunity. Immunol Rev. 2020;295(1):220–39.
MahmoudianDehkordi S, Arnold M, Nho K, Ahmad S, Jia W, Xie G, et al. Altered bile acid profile associates with cognitive impairment in Alzheimer’s disease-An emerging role for gut microbiome. Alzheimers Dement. 2019;15(1):76–92.
Zangerolamo L, Vettorazzi JF, Rosa LRO, Carneiro EM, Barbosa HCL. The bile acid TUDCA and neurodegenerative disorders: an overview. Life Sci. 2021;272: 119252.
Han X, Wang J, Gu H, Guo H, Cai Y, Liao X, et al. Predictive value of serum bile acids as metabolite biomarkers for liver cirrhosis: a systematic review and meta-analysis. Metabolomics. 2022;18(7):43.
Koike S, Miyaji Y, Sano H, Aikawa N, Kai M, Kasahara S, et al. Simultaneous determination of five bile acids as potential biomarkers for Alzheimer’s disease in mouse brain and plasma. Anal Sci. 2021;37(8):1165–70.
Tao Y, Zheng F, Cui D, Huang F, Wu X. A combination of three plasma bile acids as a putative biomarker for schizophrenia. Acta Neuropsychiatr. 2021;33(1):51–4.
Cuevas E, Burks S, Raymick J, Robinson B, Gomez-Crisostomo NP, Escudero-Lourdes C, et al. Tauroursodeoxycholic acid (TUDCA) is neuroprotective in a chronic mouse model of Parkinson’s disease. Nutr Neurosci. 2022;25(7):1374–91.
Keene CD, Rodrigues CM, Eich T, Chhabra MS, Steer CJ, Low WC. Tauroursodeoxycholic acid, a bile acid, is neuroprotective in a transgenic animal model of Huntington’s disease. Proc Natl Acad Sci U S A. 2002;99(16):10671–6.
Lo AC, Callaerts-Vegh Z, Nunes AF, Rodrigues CM, D’Hooge R. Tauroursodeoxycholic acid (TUDCA) supplementation prevents cognitive impairment and amyloid deposition in APP/PS1 mice. Neurobiol Dis. 2013;50:21–9.
Jones BV, Begley M, Hill C, Gahan CG, Marchesi JR. Functional and comparative metagenomic analysis of bile salt hydrolase activity in the human gut microbiome. Proc Natl Acad Sci U S A. 2008;105(36):13580–5.
Ridlon JM, Kang DJ, Hylemon PB. Bile salt biotransformations by human intestinal bacteria. J Lipid Res. 2006;47(2):241–59.
Ridlon JM, Alves JM, Hylemon PB, Bajaj JS. Cirrhosis, bile acids and gut microbiota: unraveling a complex relationship. Gut Microbes. 2013;4(5):382–7.
Kriaa A, Bourgin M, Potiron A, Mkaouar H, Jablaoui A, Gerard P, et al. Microbial impact on cholesterol and bile acid metabolism: current status and future prospects. J Lipid Res. 2019;60(2):323–32.
Dediukhina EG, Chistiakova TI, Vainshtein MB. Biosynthesis of arachidonic acid by micromycetes (review). Prikl Biokhim Mikrobiol. 2011;47(2):125–34.
Sakuradani E, Ando A, Ogawa J, Shimizu S. Improved production of various polyunsaturated fatty acids through filamentous fungus Mortierella alpina breeding. Appl Microbiol Biotechnol. 2009;84(1):1–10.
Sakuradani E, Shimizu S. Single cell oil production by Mortierella alpina. J Biotechnol. 2009;144(1):31–6.
Zhao M, Dai C-c, Guan X-y, Tao J. Genome shuffling amplifies the carbon source spectrum and improves arachidonic acid production in Diasporangium sp. Enzyme Microbial Technol. 2009;45(6–7):419–25.
Cmwthgr DR. Fungal production of eicosapentaenoic and arachidonic acids from industrial waste streams and crude soybean oil. Bioresour Technol. 1999;67:101–10.
Saelao S, Kanjana-Opas A, Kaewsuwan S. Optimization of biomass and arachidonic acid production by Aureispira maritima using response surface methodology. J Am Oil Chem Soc. 2010;88(5):619–29.
Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, et al. Diversity of the human intestinal microbial flora. Science. 2005;308(5728):1635–8.
Ley RE, Peterson DA, Gordon JI. Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell. 2006;124(4):837–48.
Guo P, Zhang K, Ma X, He P. Clostridium species as probiotics: potentials and challenges. J Anim Sci Biotechnol. 2020;11:24.
Foley MH, O’Flaherty S, Allen G, Rivera AJ, Stewart AK, Barrangou R, et al. Lactobacillus bile salt hydrolase substrate specificity governs bacterial fitness and host colonization. Proc Natl Acad Sci U S A. 2021. https://doi.org/10.1073/pnas.2017709118.
Acknowledgements
We would like to thank all of the participants for donating samples.
Funding
This work was supported by the Natural Science Foundation of Shandong Province (No. ZR2020MC011), the Key Research and Development Project of Shandong Province (No. 2021ZDSYS27), the Innovation Project of the Shandong Academy of Medical Sciences (No. 2022), and the Academic Promotion Programme of Shandong First Medical University (No. 2019LJ001).
Author information
Authors and Affiliations
Contributions
JF: Formal analysis, Validation, Writing—original draft. JS YC: Data curation, Investigation. CY: Methodology. QW and JH: Resource, Writing—review & editing. GC: Project administration, Supervision, Writing—review & editing. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was reviewed and approved by the Ethics Committee of Qilu Hospital of Shandong University (ref. 2020047). Each participant signed an informed consent form.
Consent for publication
We have obtained consent to publish from the participants to report patient data.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
The cohort of CIDP patients and healthy controls enrolled in this study. Table S2. Enriched metabolites and their pathways identified in the negative ion mode. Table S3. Enriched metabolites and their pathways identified in the positive ion mode. Table S4. Enriched genes involved in virulence factors in CIDP subjects. Table S5. KEGG 3 pathways with significant changes in gene abundance. Table S6. Enriched genes involved in secretion system in CIDP subjects. Figure S1. Rarefaction curve Boxplot. Abscissa represents for sample size while ordinate represents for number of species in sample. Diversity is limited when sample size is small, which is not reliable to represent for the entire microbiota structure. When rarefaction curve tends to be steady, it indicates that the sampling quantity is sufficient. Figure S2. NMDS and Anosim analysis. Anosim analysis at the family-level (A), genus-level (C) and the species-level (E). Between represents the distance between CIDP and non-CIDP groups, and the remaining boxes represent the distance within the corresponding group. NMDS analysis at the family-level (B), genus-level (D) and the species-level (F). Scales on X-axis and Y-axis on NMDS graph are the projection axes of samples in 2D. Figure S3. Microbial classifications and abundance in individual samples shown at family-level (A), genus-level (B) and the species-level (C). Figure S4. Box plot of enriched pathways in CIDP of KEGG secondary classification. The gene number can be viewed on the X-axis, the secondary classification can be viewed on Y-axis. Figure S5. Comparisons of the relative Eukaryota-Bacteria ratio in the stool of non-CIDP and CIDP group. Data are shown as mean ± SEM.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Fu, J., Shan, J., Cui, Y. et al. Metabolic disorder and intestinal microflora dysbiosis in chronic inflammatory demyelinating polyradiculoneuropathy. Cell Biosci 13, 6 (2023). https://doi.org/10.1186/s13578-023-00956-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13578-023-00956-1