- Open Access
Alterations in the gut microbiota and metabolite profiles in the context of neuropathic pain
Molecular Brain volume 14, Article number: 50 (2021)
The aim of this study was to explore the relationships among gut microbiota disturbances and serum and spinal cord metabolic disorders in neuropathic pain. 16S rDNA amplicon sequencing and serum and spinal cord metabolomics were used to identify alterations in the microbiota and metabolite profiles in the sham rats and the chronic constriction injury (CCI) model rats. Correlations between the abundances of gut microbiota components at the genus level, the levels of serum metabolites, and pain-related behavioural parameters were analysed. Ingenuity pathway analysis (IPA) was applied to analyse the interaction networks of the differentially expressed serum metabolites. First, we found that the composition of the gut microbiota was different between rats with CCI-induced neuropathic pain and sham controls. At the genus level, the abundances of Helicobacter, Phascolarctobacterium, Christensenella, Blautia, Streptococcus, Rothia and Lactobacillus were significantly increased, whereas the abundances of Ignatzschineria, Butyricimonas, Escherichia, AF12, and Corynebacterium were significantly decreased. Additionally, 72 significantly differentially expressed serum metabolites and 17 significantly differentially expressed spinal cord metabolites were identified between the CCI rats and the sham rats. Finally, correlation analysis showed that changes in the gut microbiota was significantly correlated with changes in serum metabolite levels, suggesting that dysbiosis of the gut microbiota is an important factor in modulating metabolic disturbances in the context of neuropathic pain. In conclusion, our research provides a novel perspective on the potential roles of the gut microbiota and related metabolites in neuropathic pain.
Neuropathic pain is defined as pain caused by a lesion or disease of the somatosensory nervous system and has a prevalence of 6.9–10% in the general population [1, 2]. Neuropathic pain may result from exposure to toxins or many different kinds of diseases, including metabolic diseases; neurodegenerative, vascular or autoimmune disorders; tumours; trauma; infection; and hereditary diseases . Neuropathic pain is manifested as spontaneous or evoked pain, hyperalgesia (an increased pain response to a noxious stimulus), and allodynia (a painful response elicited by a normally non-nociceptive stimulus), which affect the health and life quality of patients. However, the mechanisms underlying neuropathic pain remain unclear, and the pharmacologic agents recommended as first-line treatments exhibit less than satisfactory analgesic effects [4, 5]. Therefore, a comprehensive understanding of the pathogenesis of neuropathic pain is urgently needed to develop more effective therapeutic strategies.
The adult human gut is colonized by a large variety of commensal microorganisms, which are collectively called the gut microbiota . The gut microbiota plays an essential role in bidirectional communication between the gut and the central nervous system, including the brain and spinal cord, via immunological, hormonal and neuronal signals [7, 8]. Dysfunction of this bidirectional communication is involved in numerous diseases of the nervous system, such as Parkinson’s disease, Alzheimer's disease, stroke, multiple sclerosis and spinal cord injury [9,10,11,12,13]. Additionally, the gut microbiota has been reported to be a pivotal regulator that directly or indirectly mediates the development of neuropathic pain through a complex network of immune, metabolic, endocrine, and neural signalling pathways . To date, the specific mechanisms underlying gut microbiota-mediated pain progression remain largely unknown.
The gut microbiota regulates several metabolic and neurological signalling pathways in the host that can be associated with neuropathic pain . The spinal dorsal horn (SDH) is the primary centre for the processing and transmission of pain perception and plays an essential role in the initiation and maintenance of neuropathic pain. The chronic constriction injury (CCI) model, a classical model that was initially developed by Bennett and Xie, is the most widely applied model of peripheral nerve injury-induced neuropathic pain [15, 16]. In this study, we explored the relationship among gut microbiota disturbances and serum and spinal cord metabolic disorders in the CCI model to reveal the possible mechanism through which the gut microbiota mediates neuropathic pain 16S rDNA sequencing and emtaboloimcs.
Materials and methods
Sixteen 5- to 6-week-old male SD rats weighing 180–220 g (Hunan SJA Laboratory Animal Co., Ltd, SCXK (Hunan) 2016-0002) were housed in a specific pathogen-free (SPF) laboratory and randomly divided into the sham (control) group and CCI group (n = 8 in each group). As previously described in the literature, the rats were anaesthetized with pentobarbitone sodium and fixed to an operating table . Then, the skin of the lower 1/3 of the left thigh was incised to expose the sciatic nerve. The nerve was ligated by making four knots with 4–0 silk thread at 1-mm intervals until the surrounding muscles twitched briefly. In the sham rats, the left sciatic was exposed but not subjected to any additional manipulation.
Pain-related behavioural parameters, including the mechanical withdrawal threshold (MWT) and thermal withdrawal latency (TWL), were assessed the day before and on the 3rd, 7th, 11th, and 15th days after the operation.
To assess the MWT, the rats were placed in a transparent plexiglass box (22 cm × 12 cm × 22 cm) with a metal mesh floor for a 30-min adaptation period. An electronic von Frey anaesthesiometer (IITC Life Science Instruments, Woodland Hills, CA, USA) was used to stimulate the plantar surface of the left hind paw of each rat five times at five-minute intervals with increasing intensity. Lifting or licking of the paw indicated a positive response. The threshold on the display screen was recorded, and the average value was calculated as the MWT.
To assess the TWL, the rats were placed in a transparent plexiglass box (20 cm × 15 cm × 18 cm) with a glass plate on the bottom for a 30-min adaptation period. A radiant heat source (Model 390, IITC Life Science Instruments, Woodland Hills, CA, USA) was used to stimulate the plantar surface of the left hind paw of each rat three times at ten-minute intervals. When the rat lifted or licked its paw, the stop button was pressed, and the displayed time was recorded. The maximum irradiation time was 20 s to avoid burning. The average value was calculated as the TWL.
Sample collection and preparation
All samples were collected after completion of the behavioural tests on the 15th postoperative day. Fresh faecal samples were collected in dry and clean tubes via abdominal massage. Each sample was frozen in liquid nitrogen and stored at − 80 °C for 16S rDNA amplicon sequencing. Then, rats were anaesthetized with pentobarbitone sodium. Ipsilateral spinal cord samples were collected and washed in polybutylene succinate (PBS). Blood samples were obtained from the abdominal aorta and centrifuged at 3000×g for 10 min to extract serum. The spinal cord and serum samples were frozen in liquid nitrogen and stored at − 80 °C for UPLC-Q-TOF/MS analysis.
16S rDNA amplicon sequencing
Total genomic DNA was extracted from each sample using the cetyltrimethylammonium bromide (CTAB) method and purified . The DNA purity and concentration were measured using agarose gel electrophoresis. The V3-V4 region of the 16S rRNA gene was amplified by PCR using specific barcoded primers. PCR was performed in a 30 μL reaction, which included 15 μL Phusion® High-Fidelity PCR Master Mix (New England Biolabs, Ipswich, MA, USA), 0.2 μM forward and reverse primers, and 10 ng template DNA. The PCR products were analysed using 2% agarose gel electrophoresis, purified with an AxyPrepDNA Gel Extraction Kit (Axygen Bioscience, Union City, USA) and sequenced on an Illumina MiSeq/HiSeq 2500 platform to generate paired-end reads.
Operational taxonomic units (OTUs) were clustered with 97% similarity, and taxonomic information was annotated using the Ribosomal Database Project (RDP) classifier. Beta diversity analysis, including principal component analysis (PCA), principal coordinate analysis (PCoA) and weighted UniFrac distance, was performed with the Quantitative Insights into Microbial Ecology (QIIME) software package. Differences between groups were analysed using T-tests, linear discriminant analysis (LDA) effect size (LEFSe), and analysis of similarity (ANOSIM).
Serum and spinal cord metabolomics
Serum (100 μL) and spinal cord (60 mg) samples were collected, mixed with adequate amounts of precooled acetonitrile/methanol (1:1, v/v), and centrifuged for 20 min at 4 °C and 14,000×g to collect the supernatant. For LC–MS analysis, the samples were separated using ultra-high-performance liquid chromatography (UHPLC, 1290 Infinity LC, Agilent Technologies, Santa Clara, CA, USA). Electrospray ionization (ESI) was used for detection in both positive and negative ion modes. Mass spectrometry analysis and metabolite identification were performed using an Agilent 6550 iFunnel Q-TOF spectrometer (Agilent Technologies, Santa Clara, CA, USA) and a Triple TOF 6600 mass spectrometer (SCIEX, Framingham, MA, USA), respectively.
The raw data were converted to the mzXML format using Proteowizard (http://proteowizard.sourceforge.net/) and imported into XCMS software for further analysis, including retention time correction, peak alignment and picking. Following Pareto-scaling preprocessing, the data were subjected to multivariate data analysis, including PCA, partial least square discriminant analysis (PLS-DA), orthogonal PLS-DA (OPLS-DA), and univariate statistical analysis, including fold change (FC) analysis and T-tests. Mean of metabolite concentrations in each group was used to calculate FC values. Differentially expressed metabolites were selected according to the following screening criteria: variable importance in projection (VIP) > 1 and P < 0.05.
Ingenuity pathway analysis (IPA) (www.ingenuity.com) was applied to analyse the interaction networks of the differentially expressed serum metabolites. Briefly, metabolite identifiers and corresponding FC expression data were exported to IPA software, and interaction networks with a score > 2 were displayed.
Correlations between the abundances of gut microbiota components at the genus level, the levels of serum metabolites, and pain-related behavior were analysed. The correlation coefficients were calculated using the Spearman algorithm in R Version 3.4.2. A correlation network was constructed with Cytoscape Version 3.5.1.
The behavioural data are expressed as the mean ± standard deviation and were statistically analysed by T-test using SPSS 21.0 software. P values of less than 0.05 were considered statistically significant.
CCI-induced neuropathic pain results in behavioural changes
The pain-related behavioural parameters, including the MWT and TWL, of the sham and CCI groups were measured preoperatively and 3, 7, 11, and 15 days postoperatively. Compared with those of the sham group, the ipsilateral MWT and TWL of the CCI group were significantly reduced on the 3rd postoperative day and remained at a low level during the 15-day period of observation (P < 0.01). As shown in Fig. 1a and b, the rats in the CCI group displayed obvious hyperalgesia.
CCI significantly alters the gut microbiota composition
To study whether the gut microbiota plays a role in CCI-induced neuropathic pain, 16S rDNA amplicon sequencing of faecal samples from the sham group (n = 8) and the CCI group (n = 8) was performed. The species accumulation curves of samples from the two groups were used to analyse the richness and evenness of the microbiological composition (Fig. 2a). The species accumulation curves indicated that the sample size was sufficient (Fig. 2b). Beta diversity analysis showed significant differences in gut microbiota composition between faecal samples from the sham and CCI groups. The PCA and PCoA results showed a clear separation of gut microbiota between rats from the sham and CCI groups (Fig. 2c, d). ANOSIM and weighted UniFrac analysis showed that the community structure and composition of the gut microbiota in the CCI group were different from those in the sham group, indicating that CCI induced dysbiosis of the gut microbiota (Fig. 2e, f).
Based on the results of species annotation, the 10 species with the greatest abundances in each group at different taxonomic levels (phylum, class, order, family, genus and species) were selected. At the phylum level, the abundances of Firmicutes and Verrucomicrobia were markedly increased, while the abundances of Bacteroidetes and Proteobacteria were significantly decreased in the intestines of the CCI rats (Fig. 3a). The Firmicutes/Bacteroidetes (F/B) ratio is considered an indicator of health status and is thought to reflect the degree of dysbiosis of the gut microbiota . The F/B ratio in the CCI group was significantly higher than that in the sham group (Fig. 3b).
To identify the specific gut microbiota components associated with neuropathic pain, LEFSe analysis was used to differentiate the components of the gut microbiota of the CCI rats from those of the gut microbiota of the sham rats. LEFSe analysis revealed 39 distinguishing components at different taxon levels, including 13 species enriched in the CCI rats and 26 species enriched in the sham rats (LDA > 2, P < 0.05, Fig. 3c, d). At the phylum level, the 13 species enriched in the CCI group belonged to Firmicutes (n = 10), Actinobacteria (n = 2) and Proteobacteria (n = 1), while the 26 species enriched in the sham group belonged to Proteobacteria (n = 8), Bacteroidetes (n = 5), Cyanobacteria (n = 5), Actinobacteria (n = 2), Firmicutes (n = 1) and others (n = 5).
At the genus level, the abundances of Lactobacillus (LDA = 4.30, P = 0.036), Helicobacter (LDA = 3.79, P = 0.011), Blautia (LDA = 3.60, P = 0.012), Christensenella (LDA = 3.26, P = 0.035), Phascolarctobacterium (LDA = 2.89, P = 0.008), Streptococcus (LDA = 2.78, P = 0.005), and Rothia (LDA = 2.51, P = 0.045) were significantly increased in the gut microbiota of the CCI-induced rats. The results also showed obvious decreases in the abundances of Escherichia (LDA = 3.91, P = 0.009), Corynebacterium (LDA = 2.95, P = 0.011), Ignatzschineria (LDA = 2.85, P = 0.001), AF12 (LDA = 2.68, P = 0.008), and Butyricimonas (LDA = 2.53, P = 0.002).
To explore the pathological significance of gut microbiota dysbiosis induced by CCI, the correlations between the abundances of the differential gut microbiota and pain-related behavioural parameters, including the MWT and TWL, were analysed using Spearman correlation analysis. The abundances of Ignatzschineria, Butyricimonas, Escherichia and AF12 were positively correlated with the MWT and TWL, while the abundances of Streptococcus, Phascolarctobacterium, Helicobacter, Blautia, Rothia, Lactobacillus and Christensenella were negatively correlated with the MWT and TWL (Spearman r > 0.5 and P < 0.05) (Fig. 3e).
CCI alters metabolism in the serum and spinal cord
CCI-induced neuropathic pain alters metabolism in the serum.
To identify the metabolic disorders induced by CCI, untargeted metabolomics analysis was employed to evaluate the metabolic changes in serum samples between the sham group (n = 8) and the CCI group (n = 8). A total of 11,971 positive-mode features and 9460 negative-mode features were identified in the metabolic profiles of all samples. Multivariate statistics, including PCA, PLS-DA and OPLS-DA, were applied to analyse the data. An obvious separation between the sham rats and the CCI rats in the PCA and OPLS-DA score plot was shown, indicating that CCI induced severe metabolic dysfunction (Fig. 4a–d). Furthermore, 200 permutation tests showed that these patterns had good reliability (Fig. 4e, f).
From the OPLS-DA model, potential metabolites were identified based on specific screening conditions (VIP > 1 and P < 0.05). Seventy-two significantly differential metabolites, including 41 upregulated metabolites and 31 downregulated metabolites, were identified in the CCI rats compared with the sham rats (Table 1). A clustered heat map was used to visualize the trends of changes in significantly differentially expressed metabolites between the two groups (Fig. 5).
IPA software was applied to analyse the interactions between the differentially expressed metabolites (Table 2). The results showed that these differentially expressed serum metabolites were involved in several molecular and cellular functions and were specifically enriched in these two categories. The interaction network with the top score was free radical scavenging, lipid metabolism, and molecular transport, with 18 focus molecules being involved in several canonical pathways and biological functions, including neuropathic pain signalling in dorsal horn neurons, neuroinflammation signalling pathway, release of neurotransmitters, and depolarization (Fig. 6a). The interaction network with the second highest score was immunological disease, inflammatory disease, and the inflammatory response, with 15 focus molecules being involved in the neuroinflammatory signalling pathway and release of arachidonic acid (Fig. 6b).
CCI-induced neuropathic pain alters metabolism in the spinal cord.
Untargeted metabolomics analysis was further used to identify metabolic dysfunction in the spinal cord induced by CCI. A total of 11,978 positive-mode features and 9218 negative-mode features were identified in the metabolic profiles of spinal cord samples from the sham rats and the CCI rats. PCA, PLS-DA and OPLS-DA were used to assess the variability between the two groups. In the OPLS-DA model, 17 significantly differential metabolites with the criteria of VIP > 1 and P < 0.05, including 4 upregulated metabolites and 13 downregulated metabolites, were identified in the CCI rats compared with the sham rats (Table 1, Fig. 7a).
Six common metabolites were identified by intersecting the differentially expressed metabolites in serum and spinal samples from the sham and CCI rats (Fig. 7b). Among these differentially expressed metabolites, 5, including N6,N6,N6-trimethyl-l-lysine, 3-methylhistidine, allantoin, d-quinovose and d(–)-beta-hydroxy butyric acid (BHB), exhibited consistent expression trends (Fig. 7c). In addition, 2-hydroxybutanoic acid (2-HB) was found at lower concentrations in both serum and spinal cord samples from CCI rats than those from sham rats, although the VIP values were slightly less than 1 (Fig. 7c).
Correlation of differentially abundant gut microbiota components and the levels of serum metabolites
To further study the functional significance of the metabolic dysfunction induced by gut microbiota disturbances in CCI rats, Spearman correlation analysis was used to calculate the association between the levels of 12 differentially abundant gut bacteria at the genus level and the levels of 73 differentially expressed metabolites in serum samples. Spearman correlation coefficients were visualized using a cluster heat map, with blue/green and red representing negative and positive correlations, respectively (|r|> 0.5 and P < 0.05, Fig. 8a, b). For example, the level of N6,N6,N6-trimethyl-l-lysine was significantly positively correlated with the abundances of Ignatzschineria and AF12 but negatively correlated with the abundances of Blautia, Streptococcus, and Phascolarctobacterium. In addition, we performed correlation analysis of gut microbiota disturbances, the levels of serum metabolites, and pain-related behavioural parameters (|r|> 0.5 and P < 0.05, Fig. 8c). For instance, the D(−)-BHB level was positively correlated with the TWL and the abundances of Butyricimonas and Ignatzschineria.
In recent years, several studies have revealed that the gut microbiota is a crucial modulator of the pathogenesis of neuropathic pain. It was found that in a mouse model of chemotherapy-induced peripheral neuropathy (CIPN), the abundance of Akkermansia muciniphila was decreased, promoting barrier dysfunction and systemic inflammation and subsequently driving pain sensitivity . Ding et al. reported that changes in the gut microbiota attenuated the development of CCI-induced neuropathic pain by modulating proinflammatory and anti-inflammatory T cells . However, the detailed characteristics of the gut microbiota in the context of neuropathic pain remains largely unknown. In this study, we revealed changes in the gut microbiota composition in a classical CCI model and the potential mechanism of metabolic dysfunction mediated by the gut microbiota using 16S rDNA and metabolomics sequencing.
We demonstrated that the characteristics of the gut microbiota were changed in rats with CCI-induced neuropathic pain compared with sham controls. The F/B ratio is considered a representative indicator of health status and reflects the degree of gut microbiota dysbiosis. An increased F/B ratio was previously found to be associated with numerous nervous system diseases, such as stroke, cognitive impairments and sleep deprivation [21,22,23]. In our research, the F/B ratio was significantly increased in CCI rats, and this increase corresponded to the status of neuropathic pain. We confirmed that within the phylum Firmicutes, the abundances of Phascolarctobacterium, Christensenella, and Blautia, which belong to the order Clostridiales, were increased in the CCI group. These bacteria have been reported to be implicated in other neurological and psychiatric disorders [24,25,26]. In addition, several pathogenic bacteria, including Streptococcus and Helicobacter, were enriched in the CCI rats, and the abundances of these bacteria were significantly negative correlated with the MWT and TWL. We also found that Ignatzschineria and Butyricimonas were enriched in sham rats and that the abundances of these bacteria were positively correlated with the MWT and TWL. These results are essential for further studies investigating the key role of the gut microbiota in neuropathic pain progression.
The gut microbiota regulates a series of key metabolic functions of the body, and dysbiosis of the gut microbiota is an important factor in modulating host metabolic disturbances . The result of the present study confirmed significant differences in metabolite levels in the serum and spinal cord between sham and CCI rats, suggesting that neuropathic pain may involve extensive metabolic disturbances. In our study, 72 differentially expressed serum metabolites, primarily those that regulate lipid metabolism and the inflammatory response, were identified in CCI rats. Seventeen differentially expressed metabolites, mainly those that mediate amino acid metabolism and energy production, were identified in spinal cord samples from CCI rats. We primarily focused on metabolites with consistent expression trends in serum and spinal samples and their relationship with the gut microbiota.
The abundances of Ignatzschineria and Butyricimonas were decreased in CCI rats, and these changes had a significant positive effect on the MWT and TWL. Butyricimonas is a butyrate-producing genus, and butyrate is a short-chain fatty acid with anti-inflammatory and anti-oxidative stress activities . Recently, it was reported that a reduction in the abundance of Butyricimonas induces reductions in butyrate levels and promotes neuroinflammation in some nervous system diseases, such as multiple sclerosis, Parkinson’s disease, and depression [29,30,31]. We found that the abundance of Butyricimonas was positively correlated with the level of BHB. BHB is one of the most abundant ketone bodies and is primarily synthesized in the liver via the oxidation of fatty acids . BHB is not only a passive energy carrier but also a key signalling molecule that participates extensively in neuronal function, lipid metabolism, and gene expression [32, 33]. The ketogenic diet is increasingly being used to treat neurological disorders, including seizures, Alzheimer’s disease, Parkinson’s disease, ischaemic stroke, and chronic pain [34,35,36,37,38]. We confirmed that CCI brought about a synchronous reduction in BHB levels in both the serum and spinal cord. Recently, it was also reported that BHB significantly relieves mechanical and thermal allodynia and improves locomotor function in mice with spinal cord injury, possibly by suppressing histone deacetylation and NLRP3 inflammasome activation and protecting mitochondrial function . In addition, BHB can increase the concentration of gamma-aminobutyric acid (GABA), which is the main inhibitory neurotransmitter in the central nervous system and contributes to neuropathic pain [40, 41]. Therefore, we hypothesized that Butyricimonas might play a role in the development of neuropathic pain via ketone body metabolism.
Ignatzschineria species are aerobic, rod-shaped, non-spore-forming, gram-negative bacteria belonging to the phylum Proteobacteria . At present, the roles of Ignatzschineria in human diseases are less reported than those of other bacteria and are underreported to some degree since phenotypic identification is challenging [43, 44]. Correlation analysis revealed that the abundance of Ignatzschineria was positively correlated with the levels of BHB, 3-methylhistidine, 2-HB, and N6,N6,N6-trimethyl-l-lysine, differentially expressed metabolites in the serum and spinal cord with the same decreasing tendencies. 3-Methylhistidine is one of the major l-histidine derivatives formed by the post-translational methylation of l-histidine residues . In addition to that of 3-methylhistidine, the expression of l-histidine and l-anserine, which are involved in histamine metabolism, was significantly reduced in the spinal cords of CCI rats. l-Histidine is a histamine precursor and is decarboxylated by the enzyme histidine decarboxylase to generate histamine, which acts as an important neurotransmitter to regulate a series of biological functions of the nervous system. Histamine exerts disparate effects via distinct receptor subtypes (H1-4) with different pharmacological and signal transduction properties . For example, presynaptic and post synaptic activation of the H3 receptor can reduce neuronal excitability, inhibit neuronal inflammation and produce pain relief . Several studies have demonstrated the analgesic roles of histamine in the central nervous system in multiple rodent models of neuropathic pain [47,48,49]. 2-HB is an organic acid derived from α-hydroxybutyrate, which is an early marker of insulin resistance and impaired glucose regulation . In our study, we found that the level of 2-HB was significantly decreased in both the sera and spinal cords of CCI rats, which might represent an adaptive response to lipid oxidation or oxidative stress [51, 52]. N6,N6,N6-trimethyl-l-lysine is a methylated derivative of lysine that serves as the substrate for carnitine biosynthesis. Therefore, we hypothesized that Ignatzschineria is closely related to histamine metabolism, ketone body metabolism and carnitine biosynthesis in neuropathic pain progression.
In addition, several pathogenic bacteria were enriched in the CCI rats, and the abundances of these bacteria were significantly negatively correlated with the MWT and TWL. Helicobacter, a specific genus of Proteobacteria, is involved in the occurrence and development of some nervous system diseases via multiple mechanisms, including neurotoxicity, neuroinflammation, microelement deficiency and so on [53, 54]. Our study showed that the abundance of Helicobacter was markedly increased in CCI rats and significantly positively associated with the levels of l-tyrosine and dopamine. The levels of several important metabolites involved in tyrosine metabolism, including l-tyrosine, dopamine and 1,4-dihydroxybenzene, were significantly elevated in the sera of CCI rats. l-tyrosine is the precursor of dopamine, and the latter is a key neurotransmitter in pain transmission . Streptococcus is associated with inflammatory pain and some neurotransmitters, such as serotonin [56, 57]. In our research, the abundance of Streptococcus was found to be positively correlated with the level of anthranilic acid and negatively associated with the level of kynurenic acid. Anthranilic acid and kynurenic acid are important metabolic products of l-tryptophan.
This study provides a novel perspective on the potential roles of the gut microbiota and related metabolites in neuropathic pain. The function and specific mechanism of the gut microbiota in neuropathic pain are extremely sophisticated. However, the limitations of this study need to be considered for future studies. First, faecal microbiota transplantation could be used to further verify the roles of the specific bacteria in neuropathic pain. Second, our studies were carried out in animal models, and the results require confirmation in the patients. Despite these limitations, our study makes a number of contributions to the field.
Availability of data and materials
Data openly available in a public repository that issues datasets with DOIs.
Jensen TS, Baron R, Haanpaa M, Kalso E, Loeser JD, Rice AS, Treede RA. A new definition of neuropathic pain. Pain. 2011;152(10):2204–5. https://0-doi-org.brum.beds.ac.uk/10.1016/j.pain.2011.06.017.
van Hecke O, Austin SK, Khan RA, Smith BH, Torrance N. Neuropathic pain in the general population: a systematic review of epidemiological studies. Pain. 2014;155(4):654–62. https://0-doi-org.brum.beds.ac.uk/10.1016/j.pain.2013.11.013.
Scholz J, Finnerup NB, Attal N, Aziz Q, Baron R, Bennett MI, Benoliel R, Cohen M, Cruccu G, Davis KD, et al. The IASP classification of chronic pain for ICD-11: chronic neuropathic pain. Pain. 2019;160(1):53–9. https://0-doi-org.brum.beds.ac.uk/10.1097/j.pain.0000000000001365.
Kamerman PR, Wadley AL, Davis KD, Hietaharju A, Jain P, Kopf A, Meyer AC, Raja SN, Rice AC, Smith BH, et al. World Health Organization essential medicines lists: where are the drugs to treat neuropathic pain? Pain. 2015;156(5):793–7. https://0-doi-org.brum.beds.ac.uk/10.1097/01.j.pain.0000460356.94374.a1.
Finnerup NB, Attal N, Haroutounian S, McNicol E, Baron R, Dworkin RH, Gilron I, Haanpää M, Hansson P, Jensen TS, et al. Pharmacotherapy for neuropathic pain in adults: a systematic review and meta-analysis. Lancet Neurol. 2015;14(2):162–73. https://0-doi-org.brum.beds.ac.uk/10.1016/S1474-4422(14)70251-0.
Sender R, Fuchs S, Milo R. Revised estimates for the number of human and bacteria cells in the body. Plos Biol. 2016;14(8):e1002533. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pbio.1002533.
Catanzaro R, Anzalone M, Calabrese F, Milazzo M, Capuana M, Italia A, Occhipinti S, Marotta F. The gut microbiota and its correlations with the central nervous system disorders. Panminerva Med. 2015;57(3):127–43.
Silva YP, Bernardi A, Frozza RL. The role of short-chain fatty acids from gut microbiota in gut-brain communication. Front Endocrinol (Lausanne). 2020;11:25. https://0-doi-org.brum.beds.ac.uk/10.3389/fendo.2020.00025.
Zheng W, He R, Yan Z, Huang Y, Huang W, Cai Z, Su Y, Liu S, Deng Y, Wang Q, et al. Regulation of immune-driven pathogenesis in Parkinson’s disease by gut microbiota. Brain Behav Immun. 2020;87:890–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bbi.2020.01.009.
Bostanciklioglu M. The role of gut microbiota in pathogenesis of Alzheimer’s disease. J Appl Microbiol. 2019;127(4):954–67. https://0-doi-org.brum.beds.ac.uk/10.1111/jam.
Battaglini D, Pimentel-Coelho PM, Robba C, Dos SC, Cruz FF, Pelosi P, Rocco PRM. Gut microbiota in acute ischemic stroke: from pathophysiology to therapeutic implications. Front Neurol. 2020;11:598. https://0-doi-org.brum.beds.ac.uk/10.3389/fneur.2020.00598.
Saresella M, Marventano I, Barone M, La Rosa F, Piancone F, Mendozzi L, d’Arma A, Rossi V, Pugnetti L, Roda G, et al. Alterations in circulating fatty acid are associated with gut microbiota dysbiosis and inflammation in multiple sclerosis. Front Immunol. 2020;11:1390. https://0-doi-org.brum.beds.ac.uk/10.3389/fimmu.2020.01390.
Jogia T, Ruitenberg MJ. Traumatic spinal cord injury and the gut microbiota: current insights and future challenges. Front Immunol. 2020;11:704. https://0-doi-org.brum.beds.ac.uk/10.3389/fimmu.2020.00704.
Lin B, Wang Y, Zhang P, Yuan Y, Zhang Y, Chen G. Gut microbiota regulates neuropathic pain: potential mechanisms and therapeutic strategy. J Headache Pain. 2020;21(1):103. https://0-doi-org.brum.beds.ac.uk/10.1186/s10194-020-01170-x.
Voderova L, Sandera S, Drobil T, Plevova B, Vaculin S. Skin temperature contribution to the decrease in withdrawal latency following chronic constriction injury. Physiol Behav. 2020;227:113147. https://0-doi-org.brum.beds.ac.uk/10.1016/j.physbeh.2020.113147.
Bennett GJ, Xie YK. A peripheral mononeuropathy in rat that produces disorders of pain sensation like those seen in man. Pain. 1988;33(1):87–107. https://0-doi-org.brum.beds.ac.uk/10.1016/0304-3959(88)90209-6.
Gerstein AC, Fu MS, Mukaremera L, Li Z, Ormerod KL, Fraser JA, Berman J, Nielsen K, et al. Polyploid titan cells produce haploid and aneuploid progeny to promote stress adaptation. MBIO. 2015;6(5):1315–40. https://doi.org/10.1128/mBio.01340-15.
Feng J, Zhao F, Sun J, Lin B, Zhao L, Liu Y, Jin Y, Li S, Li A, Wei Y. Alterations in the gut microbiota and metabolite profiles of thyroid carcinoma patients. Int J Cancer. 2019;144(11):2728–45. https://0-doi-org.brum.beds.ac.uk/10.1002/ijc.32007.
Ramakrishna C, Corleto J, Ruegger PM, Logan GD, Peacock BB, Mendonca S, Yamaki S, Adamson T, Ermel R, McKemy D, et al. Dominant role of the gut microbiota in chemotherapy induced neuropathic pain. Sci Rep. 2019;9(1):20324. https://0-doi-org.brum.beds.ac.uk/10.1038/s41598-019-56832-x.
Ding W, You Z, Chen Q, Yang L, Doheny J, Zhou X, Li N, Wang S, Hu K, Chen L, et al. Gut microbiota influences neuropathic pain through modulating proinflammatory and anti-inflammatory T cells. Anesth Analg. 2020. https://0-doi-org.brum.beds.ac.uk/10.1213/ANE.0000000000005155.
Park MJ, Pilla R, Panta A, Pandey S, Sarawichitr B, Suchodolski J, Sohrabji F. Reproductive senescence and ischemic stroke remodel the gut microbiome and modulate the effects of estrogen treatment in female rats. Transl Stroke Res. 2020;11(4):812–30. https://0-doi-org.brum.beds.ac.uk/10.1007/s12975-019-00760-5.
Yang X, Yu D, Xue L, Li H, Du J. Probiotics modulate the microbiota-gut-brain axis and improve memory deficits in aged SAMP8 mice. Acta Pharm Sin B. 2020;10(3):475–87. https://0-doi-org.brum.beds.ac.uk/10.1016/j.apsb.2019.07.001.
Benedict C, Vogel H, Jonas W, Woting A, Blaut M, Schurmann A, Cedernaes J. Gut microbiota and glucometabolic alterations in response to recurrent partial sleep deprivation in normal-weight young individuals. Mol Metab. 2016;5(12):1175–86. https://0-doi-org.brum.beds.ac.uk/10.1016/j.molmet.2016.10.003.
Barandouzi ZA, Starkweather AR, Henderson WA, Gyamfi A, Cong XS. Altered composition of gut microbiota in depression: a systematic review. Front Psychiatry. 2020;11:541. https://0-doi-org.brum.beds.ac.uk/10.3389/fpsyt.2020.00541.
Iglesias-Vazquez L, Van Ginkel RG, Arija V, Canals J. Composition of gut microbiota in children with autism spectrum disorder: a systematic review and meta-analysis. Nutrients. 2020;12(3):792. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12030792.
Kuznetsova V, Tyakht A, Akhmadishina L, Odintsova V, Klimenko N, Kostryukova E, Vakhitova M, Grigoryeva T, Malanin S, Vladimirtsev V, et al. Gut microbiome signature of Viliuisk encephalomyelitis in Yakuts includes an increase in microbes linked to lean body mass and eating behaviour. Orphanet J Rare Dis. 2020;15(1):327. https://0-doi-org.brum.beds.ac.uk/10.1186/s13023-020-01612-4.
Pascale A, Marchesi N, Marelli C, Coppola A, Luzi L, Govoni S, Giustina A, Gazzaruso C. Microbiota and metabolic diseases. Endocrine. 2018;61(3):357–71. https://0-doi-org.brum.beds.ac.uk/10.1007/s12020-018-1605-5.
Fu X, Liu Z, Zhu C, Mou H, Kong Q. Nondigestible carbohydrates, butyrate, and butyrate-producing bacteria. Crit Rev Food Sci Nutr. 2019;59(sup1):S130–52. https://0-doi-org.brum.beds.ac.uk/10.1080/10408398.2018.1542587.
Jangi S, Gandhi R, Cox LM, Li N, von Glehn F, Yan R, Patel B, Mazzola MA, Liu S, Glanz B, et al. Alterations of the human gut microbiome in multiple sclerosis. Nat Commun. 2016;7:12015. https://0-doi-org.brum.beds.ac.uk/10.1038/ncomms12015.
Jang JH, Yeom MJ, Ahn S, Oh JY, Ji S, Kim TH, Park HJ. Acupuncture inhibits neuroinflammation and gut microbial dysbiosis in a mouse model of Parkinson’s disease. Brain Behav Immun. 2020;89:641–55. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bbi.2020.08.015.
Yang C, Qu Y, Fujita Y, Ren Q, Ma M, Dong C, Hashimoto K. Possible role of the gut microbiota-brain axis in the antidepressant effects of (R)-ketamine in a social defeat stress model. Transl Psychiatry. 2017;7(12):1294. https://0-doi-org.brum.beds.ac.uk/10.1038/s41398-017-0031-4.
Newman JC, Verdin E. Beta-hydroxybutyrate: a signaling metabolite. Annu Rev Nutr. 2017;37:51–76. https://0-doi-org.brum.beds.ac.uk/10.1146/annurev-nutr-071816-064916.
Han YM, Ramprasath T, Zou MH. beta-hydroxybutyrate and its metabolic effects on age-associated pathology. Exp Mol Med. 2020;52(4):548–55. https://0-doi-org.brum.beds.ac.uk/10.1038/s12276-020-0415-z.
Olson CA, Vuong HE, Yano JM, Liang QY, Nusbaum DJ, Hsiao EY. The gut microbiota mediates the anti-seizure effects of the ketogenic diet. Cell. 2018;173(7):1728–41. https://0-doi-org.brum.beds.ac.uk/10.1016/j.cell.2018.06.051.
Broom GM, Shaw IC, Rucklidge JJ. The ketogenic diet as a potential treatment and prevention strategy for Alzheimer’s disease. Nutrition. 2019;60:118–21. https://0-doi-org.brum.beds.ac.uk/10.1016/j.nut.2018.10.003.
Phillips M, Murtagh D, Gilbertson LJ, Asztely F, Lynch C. Low-fat versus ketogenic diet in Parkinson’s disease: a pilot randomized controlled trial. Mov Disord. 2018;33(8):1306–14. https://0-doi-org.brum.beds.ac.uk/10.1002/mds.27390.
Shaafi S, Sharifi-Bonab M, Ghaemian N, Mokhtarkhani M, Akbari H. Early motor-behavioral outcome of ischemic stroke with ketogenic diet preconditioning: interventional animal study. J Stroke Cerebrovasc Dis. 2019;28(4):1032–9. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jstrokecerebrovasdis.
Ruskin DN, Kawamura M, Masino SA. Reduced pain and inflammation in juvenile and adult rats fed a ketogenic diet. PLoS ONE. 2009;4(12):e8349. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0008349.
Qian J, Zhu W, Lu M, Ni B, Yang J. D-beta-hydroxybutyrate promotes functional recovery and relieves pain hypersensitivity in mice with spinal cord injury. Br J Pharmacol. 2017;174(13):1961–71. https://0-doi-org.brum.beds.ac.uk/10.1111/bph.13788.
Suzuki Y, Takahashi H, Fukuda M, Hino H, Kobayashi K, Tanaka J, Ishii E. Beta-hydroxybutyrate alters GABA-transaminase activity in cultured astrocytes. Brain Res. 2009;1268:17–23. https://0-doi-org.brum.beds.ac.uk/10.1016/j.brainres.2009.02.074.
Yudkoff M, Daikhin Y, Melo TM, Nissim I, Sonnewald U, Nissim I. The ketogenic diet and brain metabolism of amino acids: relationship to the anticonvulsant effect. Annu Rev Nutr. 2007;27:415–30. https://0-doi-org.brum.beds.ac.uk/10.1146/annurev.nutr.27.061406.093722.
Gupta AK, Dharne MS, Rangrez AY, Verma P, Ghate HV, Rohde M, Patole MS, Shouche YS. Ignatzschineria indica sp. nov. and Ignatzschineria ureiclastica sp. nov., isolated from adult flesh flies (Diptera: Sarcophagidae). Int J Syst Evol Microbiol. 2011;61(Pt 6):1360–9. https://0-doi-org.brum.beds.ac.uk/10.1099/ijs.0.018622-0.
Heddema E, Janssen F, van Westreenen H. A case of Ignatzschineria bacteraemia in an unconscious man from the Netherlands. JMM Case Rep. 2016;3(3):e5043. https://0-doi-org.brum.beds.ac.uk/10.1099/jmmcr.0.005043.
Tsang CC, Tang J, Fong J, Kinne J, Lee HH, Joseph M, Jose S, Schuster R, Tang Y, Sivakumar S, Chen J, et al. Ignatzschineria cameli sp. nov., isolated from necrotic foot tissue of dromedaries (Camelus dromedarius) and associated maggots (Wohlfahrtia species) in Dubai. Int J Syst Evol Microbiol. 2018;68(11):3627–34. https://doi.org/10.1099/ijsem.0.003046.
Holecek M. Histidine in health and disease: metabolism, physiological importance, and use as a supplement. Nutrients. 2020;12(3):848. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12030848.
Obara I, Telezhkin V, Alrashdi I, Chazot PL. Histamine, histamine receptors, and neuropathic pain relief. Br J Pharmacol. 2020;177(3):580–99. https://0-doi-org.brum.beds.ac.uk/10.1111/bph.14696.
Yu J, Lou GD, Yue JX, Tang YY, Hou WW, Shou WT, Ohtsu H, Zhang SH, Chen Z. Effects of histamine on spontaneous neuropathic pain induced by peripheral axotomy. Neurosci Bull. 2013;29(3):261–9. https://0-doi-org.brum.beds.ac.uk/10.1007/s12264-013-1316-0.
Huang L, Adachi N, Nagaro T, Liu K, Arai T. Histaminergic involvement in neuropathic pain produced by partial ligation of the sciatic nerve in rats. Reg Anesth Pain Med. 2007;32(2):124–9. https://0-doi-org.brum.beds.ac.uk/10.1016/j.rapm.2006.11.009.
Yu J, Tang YY, Wang RR, Lou GD, Hu TT, Hou WW, Yue JX, Ohtsu H, Shi LY, Zhang SH, et al. A critical time window for the analgesic effect of central histamine in the partial sciatic ligation model of neuropathic pain. J Neuroinflammation. 2016;13(1):163. https://0-doi-org.brum.beds.ac.uk/10.1186/s12974-016-0637-0.
Recber T, Orgul G, Aydin E, Tanacan A, Nemutlu E, Kir S, Beksac MS. Metabolic infrastructure of pregnant women with methylenetetrahydrofolate reductase polymorphisms: a metabolomic analysis. Biomed Chromatosr. 2020;34(8):e4842. https://0-doi-org.brum.beds.ac.uk/10.1002/bmc.4842.
Ciavardelli D, Piras F, Consalvo A, Rossi C, Zucchelli M, Di Ilio C, Frazzini V, Caltagirone C, Spalletta G, Sensi S. Medium-chain plasma acylcarnitines, ketone levels, cognition, and gray matter volumes in healthy elderly, mildly cognitively impaired, or Alzheimer’s disease subjects. Neurobiol Aging. 2016;43:1–12. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neurobiolaging.2016.03.005.
Xu Y, Han J, Dong J, Fan X, Cai Y, Li J, Wang T, Zhou J, Shang J. Metabolomics characterizes the effects and mechanisms of quercetin in nonalcoholic fatty liver disease development. Int J Mol Sci. 2019;20(5):1220. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20051220.
Budzynski J, Klopocka M. Brain-gut axis in the pathogenesis of Helicobacter pylori infection. World J Gastroenterol. 2014;20(18):5212–25. https://0-doi-org.brum.beds.ac.uk/10.3748/wjg.v20.i18.5212.
Doulberis M, Kotronis G, Gialamprinou D, Polyzos SA, Papaefthymiou A, Katsinelos P, Kountouras J. Alzheimer’s disease and gastrointestinal microbiota; impact of Helicobacter pylori infection involvement. Int J Neurosci. 2020. https://0-doi-org.brum.beds.ac.uk/10.1080/00207454.2020.1738432.
Sagheddu C, Aroni S, De Felice M, Lecca S, Luchicchi A, Melis M, Lisa Muntoni A, Romano R, Palazzo E, Guida F, et al. Enhanced serotonin and mesolimbic dopamine transmissions in a rat model of neuropathic pain. Neuropharmacology. 2015;97:383–93. https://0-doi-org.brum.beds.ac.uk/10.1016/j.neuropharm.2015.06.003.
Chakravarthy K, Faltus R, Robinson G, Sevilla R, Shin J, Zielstorff M, Byford A, Leccese E, Caniga M, Hseih SC, et al. Etanercept ameliorates inflammation and pain in a novel mono-arthritic multi-flare model of streptococcal cell wall induced arthritis. BMC Musculoskelet Disord. 2014;15:409. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2474-15-409.
Guo R, Chen LH, Xing C, Liu T. Pain regulation by gut microbiota: molecular mechanisms and therapeutic potential. Br J Anaesth. 2019;123(5):637–54. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bja.2019.07.026.
We thank Phd Xueqing Duan of Guizhou University of Traditional Chinese Medicine for help with our study.
The present study was supported in part by a grant from doctor startup funds of Guizhou University of traditional Chinese medicine (No (2019) 125).
Ethics approval and consent to participate
The study was reviewed and approved by Animal Use and Ethics Committee, Affiliated Hospital of Guangzhou University of Chinese Medicine.
Consent for publication
The authors declare no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, P., Wang, C., Ren, Yn. et al. Alterations in the gut microbiota and metabolite profiles in the context of neuropathic pain. Mol Brain 14, 50 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s13041-021-00765-y
- Neuropathic pain
- Gut microbiota
- Metabolite profiles
- CCI model