The causal connections between gut microbiota and chronic pain and pain phenotypes through the gut-brain axis: A Mendelian randomization study and mediation analysis

Jia Chen, Xuemei Luo, Luyao Wang, Yiping Luo, Xia Deng, Yuan Tian, Qiao Lu, Genwen Sun, Michael Maes, Xu Zhang
doi: https://doi.org/10.1101/2024.09.28.24314553

Abstract

Recent studies have increasingly highlighted a potential link between gut microbiota (GM) and chronic pain. However, existing research has yielded mixed results, leaving many questions unanswered. We employed a two-sample Mendelian randomization design to investigate whether GM directly influences multisite chronic pain (MCP) and eight phenotypes. We then used a combination of meta-analysis, genetic colocalization analysis, and Bayesian weighted mendelian randomization to further strengthen the evidence for causal relationship. Additionally, we examined how GM might be associated with human brain structures. Finally, we used mediation analysis to explore whether the structure and function of specific brain regions play a role in the relationship between GM and MCP and the pain phenotypes.

We observed a significant causal association between 13 GM and chronic pain. Eight GM protected against pain phenotypes, whereas seven increased risk. The results indicate that GM could potentially affect the area, thickness, and volume of brain region cortices involved in sensory, emotional, and cognitive aspects of pain. Importantly, our findings indicated that the causal effect between the Genus Odoribacter and neck/shoulder pain could be mediated by the mean optical density in the left Fornix cres+Stria terminalis. In conclusion, different pain phenotypes were causally associated with different GM among European ancestry. Core GM may influence chronic pain by affecting complex brain networks. Notably, we constructed a basic gut-brain axis model for neck/shoulder pain, involving the Genus Odoribacter, left Fornix cres+Stria terminalis, and neck/shoulder pain.

Introduction

Chronic pain (CP) is a major public health issue with significant global attention. About 30% global population experiences CP [34]. While various causes exist, CP can be idiopathic (with unknown origin) or difficult to categorize pathophysiologically [82]. Notably, CP shows a genetic link, with pain phenotypes exhibiting heritability [49].

The gut microbiota (GM) is the human body’s most complex and densely populated microbial ecosystem [33,76]. Recent research highlights its emerging role in the pathophysiology and regulation of pain [62–63]. The gut-brain axis (GBA) is a complex network of communication channels between the gut and brain [66], disruptions in which might be linked to neurological diseases including pain [27, 66].

Preclinical studies have demonstrated a link between GM and visceral pain [25,53]. Probiotics can lessen visceral pain [80] and might even affect the development of brain regions involved in pain processing [45]. Animal studies further suggest that GM influences various chronic pain conditions, including pain caused by inflammation and nerve damage [40, 62]. Other findings hint at a possible connection to nociceptive pain, complex regional pain syndrome, and headaches [50]. This suggests that GM might directly or indirectly regulate neuronal excitability in the peripheral nervous system [61–62, 74].

Despite these findings, the causal effect of GM on pain remains unclear due to potential biases. Mendelian randomization (MR) has been employed to investigate the relationship between GM and specific pain conditions, including colorectal cancer [39], irritable bowel syndrome [42], intervertebral disc degeneration [22], osteoarthritis [38,94], ankylosing spondylitis [90], and back pain [29,78]. Not only are the conclusions of the above-mentioned studies inconsistent, but also pain with unclear pathophysiology is more common in clinical practice. This underscores the need to examine pain phenotypes.

While the microbiota’s role in pain is gaining traction, many questions remain unanswered. In chronic pain, key questions include whether the GM plays a role in chronic pain and different pain phenotypes? How can GBA models be constructed for chronic pain and different pain phenotypes? MR can analyze the potential causal relationship between GM and pain by using genetic variations as instrumental variables (IVs) [37]. Furthermore, mediation analysis [67] can help decompose the impact of GM on pain.

In this study, we aimed to (1) identify potential positive GM affecting multisite chronic pain (MCP) and eight pain phenotypes, (2) analyze the effect of potential positive GM on brain structure, and (3) develop a basic GBA model for MCP and pain phenotypes to provide practical suggestions for clinical practice.

Materials and methods

1.1 Study design

This study investigated the potential cause-and-effect link between GM and chronic pain using a two-sample MR approach. We utilized data on four relevant GM features and assessed their association with both MCP and eight specific pain phenotypes. We further evaluated the relationship between positive exposure and outcomes using meta-analysis and genetic colocalization analysis. Subsequently, we examined the relationship between potential positive GM and brain imaging-derived phenotypes, and brain structure. Finally, we explored whether the structure and function of specific brain regions impacted the relationship between positive GM and MCP, as well as the eight pain phenotypes, through mediation analysis.

1. Data sources and instruments

To minimize bias caused by sample overlap, we obtained exposure (GM) and outcome (MCP and pain phenotype) data from independent sources. Table 1 summarizes the data sources and sample sizes. Ethical approval for the original studies was publicly confirmed, and this analysis did not require separate ethical approval.

View this table:
Table 1 Description of GWAS summary statistics used for each phenotype.

2.1 MCP

We obtained data on MCP from the UK Biobank Consortium, a large collection of genetic and health information. This database includes genome-wide association studies (GWAS) for a large population of European ancestry (387,649 individuals) [31]. MCP referred to chronic pain lasting 3 months or more in the head, face, neck/shoulder, back, stomach/abdomen, hip, and knee.

2.2 Eight pain phenotypes

Data on eight individual pain conditions were collected from UK Biobank participants (ranging from 151,922 to 226,683 participants depending on the specific pain type) [49]. A validated questionnaire was utilized to assess participants’ experiences with pain in the past month that impacted their daily activities, including pain in the head, face, neck/shoulder pain, back, stomach/abdomen, hips, knees, and all over the body. We excluded participants who chose “prefer not to say” and focused on pain in seven specific locations (head, neck/shoulder, back, abdomen, face, knee, and hip). Individuals who reported pain in a specific location were classified as cases, regardless of other pain selections. Participants who reported “none of the above” were classified as controls. We also included “pain all over the body” as a potential indicator of chronic widespread pain, which is a distinct condition from localized chronic pain.

2.3 GM

Data 1: MiBioGen consortium study (14,306 Europeans): This study used 16S rRNA gene sequencing and analyzed 211 microbial taxa (9 phyla, 16 classes, 20 orders, 35 families, 131 genera) from European participants[36].

Data 2: German single-country GWAS study (8,956 Germans): This study analyzed gut bacteria using 16S rRNA gene sequencing[65]. We examined 233 different bacterial groups (4 phyla, 8 classes, 6 orders, 10 families, 29 genera; 65, 62, and 49 at 97% OTU, 99% OTU, and ASV levels) to understand their abundance. Additionally, we analyzed the presence or absence of 198 bacterial features (2 classes, 1 order, 2 families, 17 genera; 65, 62, and 49 on 97% OTU, 99% OTU, and ASV levels), based on taxonomic annotation or ASV/OTU clustering.

Data 3: Dutch Microbiome Project study (7,738 individuals): This study used shotgun metagenomic sequencing[44] and analyzed 207 microbial taxa (5 phyla, 10 classes, 13 orders, 26 families, 48 genera, and 105 species) and 205 functional pathways.

Data 4: FINRISK 2002 study (5,959 individuals): This study used shotgun metagenomic sequencing [60] and analyzed 471 distinct Genome Taxonomy Database (GTDB) taxa (11 phyla, 19 classes, 24 orders, 62 families, 146 genera, and 209 species).

2.4 Brain

2.4.1 Brain imaging-derived phenotypes

Summary statistics for brain imaging-derived (BID) phenotypes were obtained from a large GWAS conducted by Smith et al. [70] using data from the UK Biobank, which included 39,691 European individuals. This dataset comprises nearly 4,000 BID phenotypes encompassing brain structure, function, connectivity, and microstructure.

2.4.2 Human brain cortical data

We leveraged brain magnetic resonance imaging (MRI) data from a large-scale analysis (genome-wide association meta-analysis) involving over 51,000 individuals of European descent across 60 separate studies. This analysis was conducted by a research group within the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium. The study focused on two key aspects of the brain’s outer layer (cortex): the total surface area (SA) and average thickness (TH). Additionally, they examined 34 specific regions within the cortex known to have distinct functions, as defined by the Desikan-Killiany atlas [24]. For further details, please refer to the original ENIGMA study. In our analysis, we specifically used data on the SA and TH of these brain regions, accounting for overall brain size and genetic background.

2.4.3 Brain-wide volumes

Researchers developed a technique called genomic principal component analysis (PCA), which allowed them to analyze brain scans and identify general underlying patterns in brain structures (morphometry) across the entire brain; the analysis takes into account the potential underlying genetic influences [21]. This technique was applied to genetic data associated with 83 distinct brain volume measurements in over 36,000 UK Biobank participants. Unlike traditional methods that focus on genetic similarity between individuals, Genomic PCA identifies key genetic components that capture the overall variation in brain volume across different brain regions.

3. Genetic variants selection criteria

We identified genetic variations (single nucleotide polymorphisms, SNPs) that strongly influence GM using strict criteria, including a very low P-value (<1×10-5) [78], a 10000-kb window around a leading SNP, and minimal linkage disequilibrium (LD, r2 < 0.001, kb=10000) with other nearby SNPs to minimize confounding effects. Additionally, the chosen SNPs had a minor allele frequency (MAF) > 0.01 to ensure adequate statistical power. After this initial filtering, we further ensured that the variations were reflective of the effects of exposure on outcome; any SNPs with F-statistics < 10 or failure to harmonize, including those harboring incompatible alleles or palindromic characteristics with intermediate allele frequencies, were excluded. Correct inference direction from exposure to outcome was identified using the Steiger filtering approach, to lessen the influence of reverse association. Finally, we used the NHGRI-EBI Catalog (https://www.ebi.ac.uk/gwas/) to retrieve these SNPs and excluded SNPs related to confounding factors and outcomes (education, living habits, smoking, diet) [97]. Figure 1 summarizes this selection process.

Figure 1:
Figure 1:

Flow chart of study

Figure 2a
Figure 2a

The Forest plots of Meta analysis of Genus Adlercreutzia on MCP

Figure 2b
Figure 2b

The Forest plots of Meta analysis of Phylum Verrucomicrobia on MCP

Figure 2c
Figure 2c

The Forest plots of Meta analysis of Phylum Verrucomicrobia on Headache

Figure 2d
Figure 2d

The Forest plots of Meta analysis of Genus Alistipes on Facialpain

Figure 2e
Figure 2e

The Forest plots of Meta analysis of Genus Bifidobacterium on Facialpain

Figure 2f
Figure 2f

The Forest plots of Meta analysis of Class Alphaproteobacteria on Neckshoulderpain

Figure 2g
Figure 2g

The Forest plots of Meta analysis of Genus Odoribacter on Neckshoulderpain

Figure 2h
Figure 2h

The Forest plots of Meta analysis of Class Gammaproteobacteria on Abdominalpain

Figure 2i
Figure 2i

The Forest plots of Meta analysis of Genus Alistipes on Abdominalpain

Figure 2j
Figure 2j

The Forest plots of Meta analysis of Specie Eubacterium_eligens on Abdominalpain

Figure 2k
Figure 2k

The Forest plots of Meta analysis of Family Prevotellaceae on Backpain

Figure 2l
Figure 2l

The Forest plots of Meta analysis of Order Bifidobacteriales on Backpain

Figure 2m
Figure 2m

The Forest plots of Meta analysis of Genus Subdoligranulum on Hippain

Figure 2n
Figure 2n

The Forest plots of Meta analysis of Genus Dialister on Kneepain

Figure 2o
Figure 2o

The Forest plots of Meta analysis of Genus Desulfovibrio on pain all over body

Figure 2p
Figure 2p

The Forest plots of Meta analysis of Genus Ruminococcus2 on pain all over body

Figure 3a
Figure 3a

The scatter plots of MR: Genus Adlercreutzia on MCP

Figure 3b
Figure 3b

The scatter plots of MR: Phylum Verrucomicrobia on MCP

Figure 3c
Figure 3c

The scatter plots of MR: Phylum Verrucomicrobia on Headache

Figure 3d
Figure 3d

The scatter plots of MR: Genus Alistipes on Facialpain

Figure 3e
Figure 3e

The scatter plots of MR: Genus Bifidobacterium on Facialpain

Figure 3f
Figure 3f

The scatter plots of MR: Class Alphaproteobacteria on Neckshoulderpain

Figure 3g
Figure 3g

The scatter plots of MR: Genus Odoribacter on Neckshoulderpain

Figure 3h
Figure 3h

The scatter plots of MR: Class Gammaproteobacteria on Abdominalpain

Figure 3i
Figure 3i

The scatter plots of MR: Genus Alistipes on Abdominalpain

Figure 3j
Figure 3j

The scatter plots of MR: Specie Eubacterium_eligens on Abdominalpain

Figure 3k
Figure 3k

The scatter plots of MR: Family Prevotellaceae on Backpain

Figure 3l
Figure 3l

The scatter plots of MR: Order Bifidobacteriales on Backpain

Figure 3m
Figure 3m

The scatter plots of MR: Genus Subdoligranulum on Hippain

Figure 3n
Figure 3n

The scatter plots of MR: Genus Dialister on Kneepain

Figure 3o
Figure 3o

The scatter plots of MR: Genus Desulfovibrio on pain all over body

Figure 3p
Figure 3p

The scatter plots of MR: Genus Ruminococcus2 on pain all over body

Figure 4a
Figure 4a

The locuscomparer plot of coloc analysis: Genus Adlercreutzia on MCP

Figure 4b
Figure 4b

The locuscomparer plot of coloc analysis: Phylum Verrucomicrobia on MCP

Figure 4c
Figure 4c

The locuscomparer plot of coloc analysis: Phylum Verrucomicrobia on Headache

Figure 4d
Figure 4d

The locuscomparer plot of coloc analysis: Genus Alistipes on Facialpain

Figure 4e
Figure 4e

The locuscomparer plot of coloc analysis: Genus Bifidobacterium on Facialpain

Figure 4f
Figure 4f

The locuscomparer plot of coloc analysis: Genus Odoribacter on Neckshoulderpain

Figure 4g
Figure 4g

The locuscomparer plot of coloc analysis: Genus Alistipes on Abdominalpain

Figure 4h
Figure 4h

The locuscomparer plot of coloc analysis: Specie Eubacterium_eligens on Abdominalpain

Figure 4i
Figure 4i

The locuscomparer plot of coloc analysis: Order Bifidobacteriales on Backpain

Figure 4j
Figure 4j

The locuscomparer plot of coloc analysis: Genus Subdoligranulum on Hippain

Figure 4k
Figure 4k

The locuscomparer plot of coloc analysis: Genus Dialister on Kneepain

Figure 4l
Figure 4l

The locuscomparer plot of coloc analysis: Genus Desulfovibrio on pain all over body

Figure 4m
Figure 4m

The locuscomparer plot of coloc analysis: Genus Ruminococcus2 on pain all over body

4. Statistical analyses

We used R version 4.2.0 for all statistical analyses. We employed various statistical methodologies to evaluate the causal influence of GM on chronic pain phenotypes. Commonly used approaches included the multiplicative random effects–inverse-variance weighted (MRE-IVW) [56], weighted mode (WM) [26], weighted median (WME) [8], and MR-Egger regression [7]. Additionally, we incorporated more recent methods like contamination mixture (ConMix) [14], MR-robust adjusted profile score (MR-RAPS)[96], debiased inverse-variance weighted (dIVW) method [94], and cML-MA [93] to enhance the robustness of our findings. False discovery rate (FDR) correction was applied to strengthen the statistical power of the results. Given the potential bias in IVW estimation when pleiotropic IVs were utilized [13], several sensitivity analyses were conducted to address pleiotropy in the causal estimates, including Cochran’s Q statistic for heterogeneity (P>0.05) [9], MR-Egger regression for horizontal pleiotropy [7]and its intercept term (P>0.05), MR-PRESSO [84] for further evaluation of horizontal pleiotropy, leave-one-out analysis to evaluate the influence of individual SNPs on the observed associations, and power calculations using a non-centrality parameter-based approach [11] and publicly available mRnd web tool (http://cnsgenomics.com/shiny/mRnd/). Any results with evidence of confounding by pleiotropy were excluded (Online Methods). We also limited our analysis to pairs of variables supported by at least three genetic variants.

We used meta-analysis to identify the positive GM, which were then subjected to mediation analysis. To verify whether the screened positive GM and pain phenotype share the same causal variation in a given region (50 kb above and below each SNP), we employed genetic colocalization analysis (coloc). Furthermore, Bayesian weighted MR (BWMR) was implemented to strengthen the evidence for a causal relationship between the identified positive GM and pain phenotypes [95]. We also examined the impact of these potential positive GM on specific brain regions within the cerebral cortex.

To assess the potential mediating role of the intermediary effect (IE) on the relationship between the exposure (a) and the outcome (c’), we conducted a mediation analysis using an interactive mediation framework and Sobel tests (http://quantpsy.org/sobel/sobel.htm) (Figure 1(b-2). The IVW method was employed to estimate the causal effects, including IE (a*b), direct effect (DE; c’ - (a*b)), and total effect (TE; c’). Statistical significance was determined by a two-sided P-value threshold of 0.05.

We employed FDR correction to account for multiple comparisons when evaluating statistical significance. Associations with an FDR of <0.05 were considered statistically significant. However, we expanded the p-value of the meta-analysis to 0.1 when no p-value was less than 0.05. We acknowledge the exploratory nature of this study and report both raw P and FDR values.

5 Results

5.1 Preliminary results of MR: GM on MCP and eight pain phenotypes

5.1.1 GM on MCP

We observed that 4 taxa (2 orders and 2 families) were related to MCP in Data 1, 8 taxa (2 genera, 2 OTU97, 2 OTU99, and 2 TestASV) were related to MCP in Data 2, 11 taxa (1 phylum, 1 class, 1 order, 2 families, 2 genera, and 4 species) were related to MCP in Data 3, and 18 taxa (1 class, 2 orders, 1 family, 7 genera, and 7 unclassified) were related to MCP in Data 4. (see Supplementary 1)

5.1.2 GM on Headache

We observed that 7 taxa (1 phylum, 1 class, 1 order, 2 families, and 2 genera) were related to headache in Data 1, 9 taxa (1 class, 1 order, 1 family, 1 OTU97, 1 OTU99, and 4 TestASV) were related to headache in Data 2, 7 taxa (1 order, 1 family, 2 genera, and 3 species) were related to headache in Data 3, and 17 taxa (1 order, 2 families, 6 genera, 2 species, and 6 unclassified) were related to headache in Data 4.(see Supplementary 2)

5.1.3 GM on Facial Pain

We observed that 6 taxa (1 order, 2 families, and 3 genera) were related to facial pain in Data 1, 9 taxa (1 genus, 4 OTU97, 3 OTU99, and 1 TestASV) were related to facial pain in Data 2, 5 taxa (1 family, 1 genus, and 3 species) were related to facial pain in Data 3, and 22 taxa (2 phyla, 1 order, 2 families, 9 genera, 2 species, and 7 unclassified) were related to facial pain in Data 4. (see Supplementary 3)

5.1.4 GM on Neck/Shoulder Pain

We observed that 4 taxa (1 phylum, 2 families, and 1 genus) were related to neck/shoulder pain in Data 1, 5 taxa (1 class, 1 family, 1 OTU97, 1 OTU99, and 1 TestASV) were related to neck/shoulder pain in Data 2, 5 taxa (1 class, 1 order, 1 family, 1 genus, and 1 species) were related to neck/shoulder pain in Data 3, and 20 taxa (1 family, 10 genera, 2 species, and 7 unclassified) were related to neck/shoulder pain in Data 4.(see Supplementary 4)

5.1.5 GM on Abdominal Pain

We observed that 8 taxa (1 class, 2 orders, 2 families, and 3 genera) were related to abdominal pain in Data 1, 8 taxa (1 class, 1 genus, 3 OTU97, 2 OTU99, and 1 TestASV) were related to abdominal pain in Data 2, 3 taxa (3 species) were related to abdominal pain in Data 3, and 19 taxa (2 orders, 1 family, 9 genera, 3 species, and 4 unclassified) were related to abdominal pain in Data 4. (see Supplementary 5)

5.1.6 GM on Back Pain

We observed that 5 taxa (1 phylum, 1 class, 1 order, and 2 genera) were related to back pain in Data 1, 6 taxa (1 family, 1 genus, 2 OTU97, and 2 OTU99) were related to back pain in Data 2, 3 taxa (1 order, 1 family, and 1 species) were related to back pain in Data 3, and 16 taxa (2 classes, 2 families, 7 genera, 1 species, and 4 unclassified) were related to back pain in Data 4. (see Supplementary 6)

5.1.7 GM on Hip Pain

We observed that 7 taxa (1 family and 6 genera) were related to hip pain in Data 1, 8 taxa (3 OTU97, 3 OTU99, and 2 TestASV) were related to hip pain in Data 2, 8 taxa (1 family, 4 genera, and 3 species) were related to hip pain in Data 3, and 17 taxa (2 orders, 2 families, 8 genera, and 5 unclassified) were related to hip pain in Data 4. (see Supplementary 7)

5.1.8 GM on Knee Pain

We observed that 11 taxa (2 phyla, 2 classes, 3 orders, and 4 genera) were related to knee pain in Data 1, 3 taxa (3 OTU99) were related to knee pain in Data 2, 5 taxa (1 family and 4 genera) were related to knee pain in Data 3, and 21 taxa (2 orders, 5 families, 9 genera, and 5 unclassified) were related to knee pain in Data 4. (see Supplementary 8)

5.1.9 GM on Pain All Over the Body

We observed that 5 taxa (5 genera) were related to pain all over the body in Data 1, 11 taxa (1 class, 1 order, 1 family, 2 OTU97, 5 OTU99, and 1 TestASV) were related to pain all over the body in Data 2, 10 taxa (1 class, 2 orders, 3 families, 2 genera, and 2 species) were related to pain all over the body in Data 3, and 21 taxa (11 genera, 3 species, and 7 unclassified) were related to pain all over the body in Data 4. (see Supplementary 9) After adjusting the p-values by FDR, Genus Odoribacter was still significantly associated with neck/shoulder pain (PIVW = 1.203×10^-4, PFDR = 0.048).

5.2 Results of meta-analysis: GM on MCP and eight pain phenotypes

In the second step, we conducted a meta-analysis based on the preliminary results from Section 5.1 in this exploratory study (see Table 2, Figure 2 and 3).

View this table:
Table 2: Mendelian randomization results of positive gut microbiota abundance based on Meta analysis on MCP and 8 pain phenotypes
View this table:
Table 3: genetic colocalisation analysis results of positive gut microbiota abundance based on Meta analysis on MCP and 8 pain phenotypes

MCP: The abundance of Genus Adlercreutzia (PIVW = 3.71E-02, PMeta = 2.24E-02) was positively associated with MCP, while that of Phylum Verrucomicrobia (PIVW = 1.90E-02, PMeta = 1.12E-02) was negatively associated with MCP.

Headache: The abundance of Phylum Verrucomicrobia (PIVW = 2.02E-02, PMeta = 2.32E-02) was negatively associated with headache.

Facial Pain: The abundances of Genus Alistipes (PIVW = 1.17E-02, PMeta = 4.68E-02) and Genus Bifidobacterium (PIVW = 1.43E-02, PMeta = 4.48E-02) were positively associated with facial pain.

Neck/Shoulder Pain: The abundances of Class Alphaproteobacteria (PIVW = 4.92E-02, PMeta = 1.17E-02) and Genus Odoribacter (PIVW = 1.20E-04, PMeta = 5.54E-05) were negatively associated with neck/shoulder pain.

Abdominal Pain: The abundances of Class Gammaproteobacteria (PIVW = 3.10E-02, PMeta = 1.95E-02) and Species Eubacterium eligens (PIVW = 3.72E-02, PMeta = 1.18E-02) were positively associated with abdominal pain, while that of Genus Alistipes (PIVW = 3.46E-02, PMeta = 3.71E-02) was negatively associated with abdominal pain.

Back Pain: The abundance of Order Bifidobacteriales (PIVW = 3.99E-02, PMeta = 4.83E-02) was positively associated with back pain.

Pain All Over the Body: The abundance of Genus Ruminococcus2 (PIVW = 1.46E-03, PMeta = 2.81E-03) was positively associated with pain all over the body, while Genus Desulfovibrio (PIVW = 7.01E-03, PMeta = 1.49E-02) was negatively associated with pain all over the body.

In this exploratory study, we expanded the P-value threshold for the meta-analysis to 0.1. The following associations were observed:

Back Pain: The abundance of Family Prevotellaceae (PIVW = 1.15E-02, PMeta = 5.20E-02) was negatively associated with back pain.

Hip Pain: The abundance of Genus Subdoligranulum (PIVW = 4.09E-02, PMeta = 8.41E-02) was negatively associated with hip pain.

Knee Pain: The abundance of Genus Dialister (PIVW = 3.94E-02, PMeta = 7.69E-02) was negatively associated with knee pain.

5.3 Results of BWMR and coloc analysis: GM on MCP and eight pain phenotypes

To further strengthen the causal inference between the identified GM features and chronic pain phenotypes, we employed both BWMR and colocalization analyses. 13 positive GMs were still causally related to MCP and eight pain phenotypes in BWMR, except for Class Alphaproteobacteria on neck/shoulder pain (PBWMR = 5.73E-02) (see Table 2).

The PP.H4.abf values of the 13 positive GMs ranged from 0.001691934 to 0.055804256. No data were found in the target region for the remaining three positive GMs. Colocalization analysis findings further strengthened the causal relationship between positive GM and MCP and eight pain phenotypes (See Table 2 and 3, Figure 3 and 4).

5.4 Results of MR: GM on BID phenotypes and human cerebral cortex

5.4.1 GM on BID Phenotypes (see Supplementary 10)

5.4.1.1 Data 1

Genus Desulfovibrio: Its abundance was positively associated with the rostral middle frontal area in the right hemisphere (PIVW = 3.79E-02) and negatively associated with rfMRI connectivity (ICA100 edge 1064) (PIVW = 3.47E-02), for which detailed information was not found.

Genus Ruminococcus2: Its abundance was positively associated with the putamen volume in the right hemisphere (PIVW = 2.94E-02), insula volume in the left (PIVW = 3.19E-02) and right hemispheres (PIVW = 4.72E-02), and G+S-occipital-inf area in the right hemisphere (PIVW = 2.17E-02).

5.4.1.2 Data 2

Class Gammaproteobacteria: Its abundance was positively associated with the medial geniculate nucleus (MGN) volume in the right hemisphere (PIVW = 1.03E-02) and negatively associated with the S-intrapariet+P-trans volume in the left hemisphere (PIVW = 4.87E-02), as well as positively associated with connectivity ICA100 edge 819 (PIVW = 4.90E-02) and negatively associated with connectivity ICA100 edge 1063 (PIVW = 2.89E-02), for which detailed information was not found.

Family Prevotellaceae: Its abundance was positively associated with the weighted mean L1 in trace left posterior thalamic radiation (PIVW = 1.67E-02).

5.4.1.3 Data 3

Genus Odoribacter: Its abundance was positively associated with mean OD in the Fornix crescent+stria terminalis (left) (PIVW = 1.49E-04).

Genus Subdoligranulum: Its abundance was positively associated with the weighted mean L2 in trace right cingulate gyrus part cingulum (PIVW = 3.00E-02).

Genus Dialister: Its abundance was negatively associated with the mean TH of G-pariet-inf-Supramar in the left hemisphere (PIVW = 4.79E-02).

Order Bifidobacteriales: Its abundance was negatively associated with the hippocampal tail volume in the right hemisphere (PIVW = 3.11E-03).

5.4.1.4 Data 4

Genus Alistipes: Its abundance was positively associated with the weighted mean diffusion tensor mode (MO) in the left superior longitudinal fasciculus (PIVW = 3.01E-02) and mean OD in the Fornix crescent+stria terminalis (right) on fractional anisotropy (FA) skeleton (PIVW = 4.94E-02) and negatively associated with the mean TH of the posterior cingulate in the left hemisphere (PIVW = 1.49E-02), mean MO in the inferior cerebellar peduncle (left) on FA skeleton (PIVW = 4.72E-02), corticoamygdaloid transition volume in the right hemisphere (PIVW = 1.12E-03), and medial orbitofrontal cortex volume in the left hemisphere (PIVW = 4.02E-03).

Genus Bifidobacterium: Its abundance was positively associated with connectivity (ICA100 edge 405) (PIVW = 2.31E-02), for which detailed information was not found, and negatively associated with the mean TH of the parstriangularis (PIVW = 1.69E-03) and S-circular-insula-sup (PIVW = 4.09E-03) in the left hemisphere; mean TH of the G-insular-short in the right hemisphere (PIVW = 3.36E-03); mean OD in the body of the corpus callosum (PIVW = 3.06E-02) and Fornix crescent+stria terminalis (right) on the FA skeleton (PIVW = 4.47E-02); HATA volume (PIVW = 5.89E-03) and area of isthmus cingulate (PIVW = 2.66E-03) in the left hemisphere; G-front-middle volume (PIVW = 4.58E-02), S-circular-insula-inf volume (PIVW = 4.85E-02), and area of superior temporal (PIVW = 2.95E-02) in the right hemisphere; and rfMRI amplitudes (ICA100 node 27) (PIVW = 5.89E-03), for which detailed information was not found.

5.4.2 GM on the Human Cerebral Cortex (see Supplementary 11)

5.4.2.1 Data 1

Genus Desulfovibrio: Its abundance was positively associated with mean SA of the paracentral (PIVW = 4.57E-02) and mean TH of the frontal pole (PIVW = 2.60E-02) and negatively associated with mean SA of the superior parietal (PIVW = 3.48E-02).

Genus Ruminococcus2: Its abundance was positively associated with mean SA of the bankssts (PIVW = 3.83E-02) and fusiform (PIVW = 1.61E-02).

Phylum Verrucomicrobia: Its abundance was negatively associated with mean SA of the lateral orbitofrontal (PIVW = 3.04E-02).

5.4.2.2 Data 2

Class Gammaproteobacteria: Its abundance was positively associated with mean SA of the bankssts (PIVW = 4.80E-02), paracentral (PIVW = 4.47E-02), posterior cingulate (PIVW = 1.97E-03), and precentral (PIVW = 2.27E-02) and negatively associated with the mean TH of the isthmus cingulate (PIVW = 4.04E-02).

Family Prevotellaceae: Its abundance was negatively associated with mean TH of the inferior parietal (PIVW = 2.33E-03).

5.4.2.3 Data 3

Genus Odoribacter: Its abundance was positively associated with the mean TH of lateral occipital (PIVW = 3.22E-02) and negatively associated with the mean TH of superior frontal (PIVW = 3.46E-02).

Genus Subdoligranulum: Its abundance was negatively associated with the mean TH of the rostral middle frontal cortex (PIVW = 4.19E-02).

Genus Dialister: Its abundance was positively associated with the mean TH of middle temporal (PIVW = 4.16E-02) and negatively associated with the mean SA of inferior parietal (PIVW = 4.91E-02) and precentral (PIVW = 4.46E-02).

Species E. eligens: Its abundance was positively associated with the mean SA of posterior cingulate (PIVW = 3.02E-02).

Phylum Verrucomicrobia: Its abundance was positively associated with global full SA (PIVW = 4.83E-02) and the mean SA of caudal middle frontal (PIVW = 3.77E-02), insula (PIVW = 8.68E-03), superior temporal (PIVW = 3.40E-02), and TH of postcentral (PIVW = 2.87E-02) and rostral anterior cingulate (PIVW = 2.91E-02) and negatively associated with the mean SA of pericalcarine (PIVW = 3.32E-02).

Genus Adlercreutzia: Its abundance was positively associated with the mean SA of lateral occipital (PIVW = 3.39E-02).

5.4.2.4 Data 4

Genus Alistipes: Its abundance was positively associated with the mean SA of lateral occipital (PIVW = 2.61E-02) and negatively associated with the mean SA of fusiform (PIVW = 1.66E-03) and postcentral (PIVW = 4.62E-02).

Genus Bifidobacterium: Its abundance was positively associated with the mean TH of posterior cingulate (PIVW = 4.32E-02) and negatively associated with the mean SA of supramarginal (PIVW = 1.83E-03) and transverse temporal (PIVW = 7.67E-03).

5.4.2 GM on 83 Brain Regions (see Supplementary 12)

Class Gammaproteobacteria: Its abundance was negatively associated with the left superior parietal volume (PIVW = 3.07E-02).

Genus Bifidobacterium: Its abundance was negatively associated with the right caudal middle frontal volume (PIVW = 1.97E-02). No causal effects were observed for the other 13 genetic predispositions on 83 brain regions because their PIVWs were >0.05.

5.5 Results of mediation analysis

5.5.1 GM on MCP and eight pain phenotypes (See Table 4, Fig 5, Supplementary 13)

Finally, we conducted a mediation analysis to explore the primary mediating GBA model of MCP and eight pain phenotypes based on Sections 5.2 and 5.3. Our findings indicate that the causal effect between Genus Odoribacter and neck/shoulder pain could be mediated by mean OD in Fornix cres + Stria terminalis (left) (Table 4). For Genus Odoribacter - mean OD in Fornix cres+Stria terminalis (left) - neck/shoulder pain, TE = −0.020, IE = −0.003, DE = −0.017 (P = 0.027), IE_div_TE = 0.140 (0.013-0.267) (see Figure 5).

View this table:
Table 4 Three causal effects (a, b, c’) of the mediation analysis in IVW method
Figure 5
Figure 5

Mediator model: Genus Odoribacter to neckshoulder pain

6. Discussion

6.1. GM and pain phenotypes

This study is the first to analyze the causal association of GM with MCP and eight pain phenotypes by MR, incorporating human brain structures. Furthermore, we were able to establish a GBA model linking Genus Odoribacter to neck/shoulder pain. We identified 13 potential relevant GM. Importantly, different pain phenotypes were causally associated with specific GM. Seven positive and eight negative causal associations were found, among which Genus Odoribacter as a protective GM for neck/shoulder pain needs to be emphasized.

Previous MR studies did not identify GM overlapping with those within our study, except for Family Prevotellaceae. This discrepancy may stem from previous studies predominantly focusing on specific pain-related diseases as outcomes, while our study examined pain phenotypes. Notably, our study suggests that Family Prevotellaceae provides protection against back pain, whereas previous MR studies reported it increases the risk of low back pain [23,78]. This difference could be attributed to variations in study outcomes and methodologies; previous studies often utilized single GM databases, whereas we incorporated four different GM databases and conducted a meta-analysis.

Meanwhile, our results suggest Phylum Verrucomicrobia and Genus Dialister are protective factors for CP, while Genus Adlercreutzia is risk factor, consistent with previous observational studies [18,23,30,40–41,85]. Despite variations in associated pain sites, this reinforces the potential role of these GMs in pain pathogenesis. Conversely, previous observational studies involving Genus Desulfovibrio, Genus Subdoligranulum, and Genus Odoribacter have yielded inconsistent results [17,23,86–88]. Our study’s findings may be more compelling due to its ability to establish causation, which observational studies cannot confirm definitively as they often have small sample sizes.

Interestingly, despite its probiotic properties in inhibiting pain behaviors [55], Genus Bifidobacterium was found to increase the risk of facial pain in our study, aligning with previous inconsistent findings [23,54,59]. Similarly, Genus Adlercreutzia, recognized for its health benefits [10], has been associated with back pain [18] and identified as a risk factor for MCP in our study. A recent attempt to use Bifidobacterium in treating neuropathic pain was unsuccessful [28], suggesting that the effects of symbionts may be complicated through interactions with pathobionts that are associated with pain conditions. The role of GM in CP is complex, possibly involving differential effects within microbial networks [10]. Future research may benefit from constructing gut microbial networks based on key GMs to elucidate CP mechanisms.

6.2. The gut-brain axis

Pain is a dynamic interplay between sensory, emotional, and cognitive factors [48] processed within the pain matrix-a network of interconnected brain regions [15], especially CP. Emerging research suggests a potential influence of GM in modulating pain perception [45]. Since changes in brain structure can significantly affect its function [57], we investigated the potential connection between GM and brain regions involved in pain processing using brain imaging data. In general, except for Family Prevotellaceae, Genus Subdoligranulum, Order Bifidobacteriales, and Species E. eligens, the remaining positive GMs have been potentially causally linked to the brain regions primarily involved in sensory, emotional, and cognitive functions, which may help us partly understand how GM affect CP.

Firstly, we found that multiple positive GM could affect the area of pain-related cerebral cortex. In our study, the abundance of Phylum Verrucomicrobia, a protective factor against MCP and headache, was not only positively associated with the entire cortical area, but also positively correlated with the area of caudal middle frontal and insula. This finding extends the results of a recent study on CP utilizing the UK Biobank data [5]. Other studies suggest that the anterior insula plays a more important role in anticipating pain and its emotional features, while the posterior insula is more active during the physical sensations of pain [55, 58, 68]. The observed protective effect of Phylum Verrucomicrobia on these regions aligns with the possibility that these bacteria might help in maintaining normal brain function related to pain perception. This finding offers a potential target for future research on chronic pain management strategies.

Furthermore, our study revealed that Genus Desulfovibrio was not only positively correlated with the area of rostral middle frontal cortex which are crucial for emotional and cognitive processing, but also positively correlated with the paracentral area and negatively with the superiorparietal area belonging sensorimotor cortices. As a protective factor for chronic multisite pain, the interaction between these affected brain areas deserves further exploration.

Secondly, we found that multiple positive GM could affect the thickness and volume of pain-related cerebral cortex. CP patients often show reduced gray matter volume in several brain regions, including areas involved in processing pain, managing emotions, and controlling movement [3–4,6,20,32,35,46–47,51–52,64,70–71,73,91]. Certain GM may have either positive or negative effects on the volume of these brain regions. For example, Phylum Verrucomicrobia abundance was associated with thicker gray matter in the rostral-anterior-cingulate and postcentral regions, further supporting its protective role against CP. Additionally, multiple GM were found to significantly influence one brain region, such as the frontal lobe. This finding suggests that the observed decrease in gray matter volume in CP patients may result from the combined effects of multiple GM.

Thirdly, this study reveals that GM may influence chronic pain by affecting the integrity and direction of axonal fiber dispersion in white matter. CP patients often exhibit changes in the microstructure of white matter, like superior longitudinal fasciculus, corpus callosum, and cingulum bundle, which play crucial roles in pain severity, resilience and emotional response to pain [52,69,81,89]. Our results indicate that certain positive GM can positively (Genus Subdoligranulum-cingulum bundle (rh)) or negatively (Genus Bifidobacterium-corpus callosum) influence white matter in specific locations. Observational studies have also linked reduced white matter integrity in the splenium of the corpus callosum to chronic low back pain disability [12]. Our study found that Family Prevotellaceae abundance positively affected posterior-thalamic-radiation (lh),supporting its role as a protective factor against back pain.

Morphological changes in the limbic system are closely linked to persistent pain. Particularly, the Fornix cres+Stria terminalis region needs to be emphasized, which is positively affected by Genus Odoribacter and Genus Alistipes and negative affected by Genus Bifidobacterium. Notably, our findings suggest that the causal relationship between Genus Odoribacter abundance and neck/shoulder pain may be mediated by mean OD in the Fornix cres+Stria terminalis (left).

The fornix, crucial for normal cognitive function [80], has also been associated with pain abnormalities [43]. Similarly, the stria terminalis is integral to the limbic system, connecting the amygdala to the hypothalamus. Although the fornix and stria terminalis are not directly connected, they create an indirect anatomical pathway between the hippocampus and the amygdala. Through shared connections to structures like the hypothalamus and septal nuclei, these regions can influence each other, forming a corticolimbic system critical in the initiation, persistence, and amplification of chronic pain states [1,83], and influencing the transition from acute to chronic pain [91]. Continued activation of corticolimbic circuits results in functional and anatomical changes in the cortex, contributing to prolonged pain states [91].

Our study identified new avenues for investigating the GBA in chronic pain and pain phenotypes. The Fornix cres+Stria terminalis is highlighted as a new factor in the GBA and a novel key area for CP (Genus Odoribacter to neck/shoulder pain). Targeted therapies focusing on these areas may offer innovative approaches to pain management and treatment.

There are some contradictions in the results: Genus Alistipes appears to exert contrasting influences in abdominal and facial pain. Class Gammaproteobacteria abundance was positively correlated with the paracentral and precentral areas but acted as a risk factor for abdominal pain, and Genus Dialister abundance was negatively associated with the precentral area but served as a protective factor for knee pain. Various internal brain networks and regulatory systems collaboratively contribute to the experience and persistence of chronic pain [16,19,46]. Future research using multimodal fMRI to investigate how GM influences these brain network changes may provide insights into these findings.

Limitations and strengths

Our study has some advantages. First, we employed a comprehensive set of sensitivity analyses, including meta-analyses and colocalization analyses. These analyses provide strong evidence that our findings are robust. Second, we addressed potential biases by applying methods like MR-Egger and MR-PRESSO to assess pleiotropy (where a genetic variant influences multiple traits) and identify outliers. Those steps helped to ensure the reliability of our MR results. Finally, the analysis utilized powerful IVs (F-statistics > 10), minimizing the influence of weak instruments that could bias the MR results.

However, our study also has certain limitations. Firstly, the data originated from participants of European ancestry. Secondly, we did not conduct reverse MR analyses or MR analyses directly analyzing the relationship between pain phenotypes and brain structure, which will be explored in future studies. Thirdly, to ensure a sufficient number of IVs for GM, we employed a significance threshold of P<1×10-5, which is less stringent than the genome-wide significance level (P<5×10-8).

Conclusions

Among European ancestry, we identified 13 GM that were causally related to chronic pain; 8 GM decreased risk towards MCP and pain phenotypes, and 7 GM increased risk. Core GM may influence CP by affecting complex brain networks associated with emotional, cognitive, and motor functions. Specifically, we constructed a basic GBA model for neck/shoulder pain with Genus Odoribacter - Fornix cres+Stria terminalis (left) interactions playing a key role in neck/shoulder pain. Our study contributes to novel insights into the relationships between specific GM and pain phenotypes, supporting that GBA play a key role in the pathogenesis of CP. Consequently, the GBA is a new drug target to treat CP through neuroregulatory therapies and targeted treatments with protective symbionts.

Data Availability

All relevant data are within the manuscript and its Supplementary Material.

Acknowledgements

We wish to acknowledge the participants and investigators of the UK Biobank, FINRISK 2002 study, Enhancing NeuroImaging Genetics through Meta-Analysis Consortium, the Integrative Epidemiology Unit (IEU, https://gwas.mrcieu.ac.uk/), the GWAS Catalog (https://www.ebi.ac.uk/gwas/).

Footnotes

  • ↵# Co-first author

  • Funding: No funding was received for this research. Declaration of competing interest None declared.

  • Competing interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

  • Data availability: All relevant data are within the manuscript and its Supplementary Material.

Reference

  1. [1].↵
    Apkarian A.V., Baliki M.N., Geha P.Y.,2009. Towards a theory of chronic pain.Prog Neurobiol 87(2):81–97. doi: 10.1016/j.pneurobio.2008.09.018.
    OpenUrlCrossRefPubMedWeb of Science
  2. [2].
    Apkarian A.V., Bushnell M.C., Treede R.D., Zubieta J.K., 2005. Human brain mechanisms of pain perception and regulation in health and disease. Eur J Pain 9(4):463–84. doi: 10.1016/j.ejpain.2004.11.001.
    OpenUrlCrossRefPubMedWeb of Science
  3. [3].↵
    Apkarian A.V., Sosa Y., Sonty S., Levy R.M., Harden R.N., Parrish T.B., Gitelman D.R.,2004.Chronic Back Pain Is Associated with Decreased Prefrontal and Thalamic Gray Matter Density. J Neurosci 24(46):10410–5. doi: 10.1523/JNEUROSCI.2541-04.2004.
    OpenUrlAbstract/FREE Full Text
  4. [4].↵
    Bhatt R.R., Gupta A., Mayer E.A., Zeltzer L.K., 2020. Chronic pain in children: structural and resting-state functional brain imaging within a developmental perspective. Pediatr Res 88(6):840–849. doi: 10.1038/s41390-019-0689-9. Epub 2019 Dec 2.
    OpenUrlCrossRef
  5. [5].↵
    Bhatt R.R., Haddad E., Zhu A.H., Thompson P.M., Gupta A., Mayer E.A., Jahanshad N., 2024. Mapping Brain Structure Variability in Chronic Pain: The Role of Widespreadness and Pain Type and its Mediating Relationship with Suicide Attempt. Biol Psychiatry 95(5):473–481. doi: 10.1016/j.biopsych.2023.07.016.
    OpenUrlCrossRef
  6. [6].↵
    Borsook D., Erpelding N., Becerra L., 2013. Losses and gains: chronic pain and altered brain morphology. Expert Rev Neurother 13(11):1221–34. doi: 10.1586/14737175.2013.846218.
    OpenUrlCrossRefPubMed
  7. [7].↵
    Bowden J., Davey Smith G., Burgess S., 2015. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525.
    OpenUrlCrossRefPubMed
  8. [8].↵
    Bowden J., Davey Smith G., Haycock P.C., Burgess S., 2016. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314.
    OpenUrlCrossRefPubMed
  9. [9].↵
    Bowden J., Fabiola Del Greco M., Minelli G., Smith G.D., Sheehan N., Thompson J., 2017. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat. Med. 36, 1783–1802.
    OpenUrlCrossRefPubMed
  10. [10].↵
    Brandon-Mong G.J., Shaw G.T.W., Chen W.H., Chen C.C., Wang D., 2020. A network approach to investigating the key microbes and stability of gut microbial communities in a mouse neuropathic pain model.BMC Microbiol 20(1):295. doi: 10.1186/s12866-020-01981-7.
    OpenUrlCrossRef
  11. [11].↵
    Brion M.J., Shakhbazov K., Visscher P.M., 2013. Calculating statistical power in Mendelian randomization studies. Int. J. Epidemiol. 42, 1497–1501.
    OpenUrlCrossRefPubMedWeb of Science
  12. [12].↵
    Buckalew N., Haut M.W., Aizenstein H., Rosano C., Edelman K.D., Perera S., Marrow L., Tadic S., Venkatraman V., Weiner D., 2013.White matter hyperintensity burden and disability in older adults: is chronic pain a contributor? PM R 5(6):471–80; quiz 480. doi: 10.1016/j.pmrj.2013.03.004.
    OpenUrlCrossRefPubMed
  13. [13].↵
    Burgess S., Butterworth A., Thompson S.G., 2013. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37, 658–665.
    OpenUrlCrossRefPubMed
  14. [14].↵
    Burgess S., Foley C.N., Allara E., Staley J.R., Howson J.M.M., 2020. A robust and efficient method for Mendelian randomization with hundreds of genetic variants. Nat Commun 11(1):376. doi: 10.1038/s41467-019-14156-4.
    OpenUrlCrossRef
  15. [15].↵
    Bushnell M.C., Ceko M., Low L.A.,2013. Cognitive and emotional control of pain and its disruption in chronic pain. Nat Rev Neurosci 14(7):502–511. doi: 10.1038/nrn3516.
    OpenUrlCrossRefPubMed
  16. [16].↵
    Canavero S., Bonicalzi V., 2015. Pain myths and the genesis of central pain. Pain Med 16(2):240–248. doi: 10.1111/pme.12509.
    OpenUrlCrossRefPubMed
  17. [17].↵
    Chriswell M.E., Lefferts A.R., Clay M.R., Hsu A.R., Seifert J., Feser M.L., Rims C., Bloom M.S., Bemis E.A., Liu S., Maerz M.D., Frank D.N., Demoruelle M.K., Deane K.D., James E.A., Buckner J.H., Robinson W.H., Holers V.M., KuhnK.A., 2022. Clonal IgA and IgG autoantibodies from individuals at risk for rheumatoid arthritis identify an arthritogenic strain of Subdoligranulum. Sci Transl Med 14(668):eabn5166. doi: 10.1126/scitranslmed.abn5166.
    OpenUrlCrossRef
  18. [18].↵
    Dekker Nitert M., Mousa A., Barrett H. L., Naderpoor N., de Courten B., 2020. Altered gut microbiota composition is associated with Back pain in overweight and obese individuals. Front. Endocrinol. 11:605. doi: 10.3389/fendo.2020.00605.
    OpenUrlCrossRef
  19. [19].↵
    Fitzcharles M.A., Cohen S.P., Clauw D.J., Littlejohn G., Usui C., Häuser W., 2021. Nociplastic pain: towards an understanding of prevalent pain conditions. Lancet 397(10289):2098–2110. doi: 10.1016/S0140-6736(21)00392-5.
    OpenUrlCrossRefPubMed
  20. [20].↵
    Fritz J.H., Wittfeld K., Hegenscheid K., Schmidt C.O., Langner S., Lotze M.,2016.Chronic Back Pain Is Associated With Decreased Prefrontal and Anterior Insular Gray Matter: Results From a Population-Based Cohort Study. J Pain 17(1):111–8. doi: 10.1016/j.jpain.2015.10.003.
    OpenUrlCrossRef
  21. [21].↵
    Fürtjes A.E., Arathimos R., Coleman J.R.I., Cole J.H., Cox S.R., Deary I.J., de la Fuente J., Madole J.W., Tucker-Drob E.M., Ritchie S.J., 2023. General dimensions of human brain morphometry inferred from genome-wide association data.Hum Brain Mapp 44(8):3311–3323. doi: 10.1002/hbm.26283.
    OpenUrlCrossRef
  22. [22].↵
    Geng Z.M., Wang J., Chen G.D., Liu J.C., Lan J., Zhang Z.P., Miao J., 2023. Gut microbiota and intervertebral disc degeneration: a bidirectional two-sample Mendelian randomization study.J Orthop Surg Res 18(1):601. doi: 10.1186/s13018-023-04081-0.
    OpenUrlCrossRef
  23. [23].↵
    Goudman L., Demuyser T., Pilitsis J.G., Billot M., Roulaud M., Rigoard P., Moens M., 2024. Gut dysbiosis in patients with chronic pain: a systematic review and meta-analysis. Front Immunol 15:1342833. doi: 10.3389/fimmu.2024.1342833.
    OpenUrlCrossRef
  24. [24].↵
    Grasby K.L., Jahanshad N., Painter J.N., Colodro-Conde L., Bralten J., Hibar D.P., Lind P.A., et al., Genetics through Meta-Analysis Consortium - Genetics working group, 2020. The genetic architecture of the human cerebral cortex. Science 367(6484):eaay6690. doi: 10.1126/science.aay6690.
    OpenUrlAbstract/FREE Full Text
  25. [25].↵
    Grundy L., Erickson A., Brierley S.M.,2019. Visceral pain. Annu Rev Physiol 81: 261–284.doi: 10.1146/annurev-physiol-020518-114525.
    OpenUrlCrossRefPubMed
  26. [26].↵
    Hartwig F.P., Smith G.D., Bowden J., 2017. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int. J. Epidemiol. 46, 1985–1998.
    OpenUrlCrossRefPubMed
  27. [27].↵
    Heiss C.N., Olofsson L.E., 2019. The role of the gut microbiota in development, function and disorders of the central nervous system and the enteric nervous system. J Neuroendocrinol J Neuroendocrinol 31(5):e12684. doi: 10.1111/jne.12684.
    OpenUrlCrossRefPubMed
  28. [28].↵
    Huang J., Zhang C., Wang J., Guo Q., Zou W.,2019. Oral Lactobacillus reuteri LR06 or Bifidobacterium BL5b supplement do not produce analgesic effects on neuropathic and inflammatory pain in rats. Brain Behav 9(4):e01260. doi: 10.1002/brb3.1260.
    OpenUrlCrossRef
  29. [29].↵
    Hui J.N., Chen Y.J., Li C.E., Gou Y.F., Liu Y., Zhou R.X., Kang M.J., Liu C., Wang B.Y., Shi P.X., Cheng S.Q., Yang X.N., Pan C.Y., Jia Y.M., Cheng B.L., Liu H., Wen Y., Zhang F., 2023. Insight into the Causal Relationship between Gut Microbiota and Back Pain: A Two Sample Bidirectional Mendelian Randomization Study.Adv Genet (Hoboken) 4(4):2300192. doi: 10.1002/ggn2.202300192.
    OpenUrlCrossRef
  30. [30].↵
    Ivashkin V., Shifrin O., Maslennikov R., Poluektova E., Korolev A., Kudryavtseva A., Krasnov G., Benuni N., Barbara G., 2023. Eubiotic effect of rifaximin is associated with decreasing abdominal pain in symptomatic uncomplicated diverticular disease: results from an observational cohort study. BMC Gastroenterol 23(1):82. doi: 10.1186/s12876-023-02690-x.
    OpenUrlCrossRef
  31. [31].↵
    Johnston K.J.A., Ward J., Ray P.R., Adams M.J., McIntosh A.M., Smith B.H., Strawbridge R.J., Price T.J., Smith D.J., Nicholl B.I., Bailey M.E.S., 2021. Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank. PLoS Genet 17(4):e1009428. doi:10.1371/journal.pgen.1009428.
    OpenUrlCrossRef
  32. [32].↵
    Kang D., McAuley J.H., Kassem M.S., Gatt J.M., Gustin S.M.,2019. What does the grey matter decrease in the medial prefrontal cortex reflect in people with chronic pain? Eur J Pain 23(2):203–219. doi: 10.1002/ejp.1304.
    OpenUrlCrossRefPubMed
  33. [33].↵
    Karczewski J., Troost F.J., Konings I., Dekker J., Kleerebezem M., Brummer R.J.M., Wells J.M., 2010. Regulation of human epithelial tight junction proteins by Lactobacillus plantarum in vivo and protective effects on the epithelial barrier. Am J Physiol Gastrointest Liver Physiol 298(6):G851–9. doi: 10.1152/ajpgi.00327.2009.
    OpenUrlCrossRefPubMedWeb of Science
  34. [34].↵
    Kozak-Szkopek E., Broczek K., Slusarczyk P., Wieczorowska-Tobis K., Klich-Raczka A., Szybalska A., Mossakowska M., 2017. Prevalence of chronic pain in the elderly polish population - results of the polsenior study. Arch Med Sci 13(5):1197–1206. doi: 10.5114/aoms.2015.55270.
    OpenUrlCrossRef
  35. [35].↵
    Kuchinad A., Schweinhardt P., Seminowicz D.A., Wood P.B., Chizh B.A., Bushnell M.C.,2007. Accelerated Brain Gray Matter Loss in Fibromyalgia Patients: Premature Aging of the Brain? J Neurosci 27(15):4004–7. doi: 10.1523/JNEUROSCI.0098-07.2007.
    OpenUrlAbstract/FREE Full Text
  36. [36].↵
    Kurilshikov A., Medina-Gomez C., Bacigalupe R., et al., 2021. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat Genet 53(2): 156–165. doi:10.1038/s41588-020-00763-1.
    OpenUrlCrossRef
  37. [37].↵
    Lawlor D.A., Harbord R.M., Sterne J.A.C., Timpson N., Smith G.D., 2008. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology.Stat Med 27(8):1133–1163. doi: 10.1002/sim.3034.
    OpenUrlCrossRefPubMed
  38. [38].↵
    Lee Y.H., Song G.G., 2021. The Gut Microbiome and Osteoarthritis: A Two-Sample Mendelian Randomization Study.J Rheum Dis 28(2):94–100. doi: 10.4078/jrd.2021.28.2.94.
    OpenUrlCrossRef
  39. [39].↵
    Li H.F., Sheng D.H., Jin C.D., Zhao G.P., Zhang L., 2023. Identifying and ranking causal microbial biomarkers for colorectal cancer at different cancer subsites and stages: a Mendelian randomization study. Front Oncol 13:1224705. doi: 10.3389/fonc.2023.1224705.
    OpenUrlCrossRef
  40. [40].↵
    Li J.S., Su S.L., Xu Z., Zhao L.H., Fan R.Y., Guo J.M., Qian D.W., Duan J.A., 2022. Potential roles of gut microbiota and microbial metabolites in chronic inflammatory pain and the mechanisms of therapy drugs. Ther Adv Chronic Dis. 28:13:20406223221091177. doi: 10.1177/20406223221091177.
    OpenUrlCrossRef
  41. [41].↵
    Li Y., Zhang S.X., Yin X.F., Zhang M.X., Qiao J., Xin X.H., Chang M.J., Gao C., Li Y.F., Li X.F., 2021. The Gut Microbiota and Its Relevance to Peripheral Lymphocyte Subpopulations and Cytokines in Patients with Rheumatoid Arthritis. J Immunol Res 2021:6665563. doi: 10.1155/2021/6665563.
    OpenUrlCrossRef
  42. [42].↵
    Liu B., Ye D., Yang H., Song J., Sun X.H., He Z.X., Mao Y.Y., Hao G.F., 2023. Assessing the relationship between GM and irritable bowel syndrome: a two-sample Mendelian randomization analysis.BMC Gastroenterology 23(1):150. doi: 10.1186/s12876-023-02791-7.
    OpenUrlCrossRef
  43. [43].↵
    Liu J., Zhu J., Yuan F., Zhang X., Zhang Q., 2018.Abnormal brain white matter in patients with right trigeminal neuralgia: a diffusion tensor imaging study.J Headache Pain 19(1):46. doi: 10.1186/s10194-018-0871-1.
    OpenUrlCrossRef
  44. [44].↵
    Lopera-Maya E.A., Kurilshikov A., van der Graaf A., Hu S.X., Andreu-Sánchez S., Chen L.M., Vila A.V., Gacesa R., Sinha T., Collij V., Klaassen M.A.Y., Bolte L.A., Gois M.F.B., Neerincx P.B.T., Swertz M.A., LifeLines Cohort Study, Harmsen H.J.M., Wijmenga C., Fu J.Y., Weersma R.K., Zhernakova A., Sanna S., 2022. Effect of host genetics on the gut microbiome in 7,738 participants of the Dutch microbiome project. Nat Genet 54(2):143–151. doi: 10.1038/s41588-021-00992-y.
    OpenUrlCrossRef
  45. [45].↵
    Luczynski P., Tramullas M., Viola M., Shanahan F., Clarke G., O’Mahony S., Dinan T.G., Cryan J.F., 2017. Microbiota regulates visceral pain in the mouse. Elife 6:e25887. doi: 10.7554/eLife.25887.
    OpenUrlCrossRef
  46. [46].↵
    Martucci K.T., MacKey S.C., 2018. Neuroimaging of Pain: Human Evidence and Clinical Relevance of Central Nervous System Processes and Modulation. Anesthesiology 128(6):1241–1254. doi: 10.1097/ALN.0000000000002137.
    OpenUrlCrossRefPubMed
  47. [47].↵
    Mayer E.A., Gupta A., Kilpatrick L.A., Hong J.Y., 2015.Imaging brain mechanisms in chronic visceral pain. Pain 156 Suppl 1(0 1):S50–S63. doi: 10.1097/j.pain.0000000000000106.
    OpenUrlCrossRef
  48. [48].↵
    Melzack R., 2001. Pain and the neuromatrix in the brain. J Dent Educ 65(12):1378–1382.
    OpenUrlAbstract
  49. [49].↵
    Meng W., Adams M.J., Reel P., Rajendrakumar A., Huang Y., Deary I.J., Palmer C.N.A., McIntosh A.M., Smith B.H.,2020. Genetic correlations between pain phenotypes and depression and neuroticism. Eur J Hum Genet 28(3):358–366. doi: 10.1038/s41431-019-0530-2.
    OpenUrlCrossRefPubMed
  50. [50].↵
    Minerbi A., Shen S., 2022. Gut microbiome in anesthesiology and pain medicine. Anesthesiology 137(1):93–108. doi:10.1097/ALN.0000000000004204.
    OpenUrlCrossRef
  51. [51].↵
    Mutso A.A., Radzicki D., Baliki M.N., Huang L., Banisadr G., Centeno M.V., Radulovic J., Martina M., Miller R.J., Apkarian A.V.,2012. Abnormalities in hippocampal functioning with persistent pain.J Neurosci 32(17):5747–56. doi: 10.1523/JNEUROSCI.0587-12.2012.
    OpenUrlAbstract/FREE Full Text
  52. [52].↵
    Ng S.K., Urquhart D.M., Fitzgerald P.B., Cicuttini F.M., Hussain S.M., Fitzgibbon B.M.,2018.The Relationship Between Structural and Functional Brain Changes and Altered Emotion and Cognition in Chronic Low Back Pain Brain Changes A Systematic Review of MRI and fMRI Studies. Clin J Pain 34(3):237–261. doi: 10.1097/AJP.0000000000000534.
    OpenUrlCrossRefPubMed
  53. [53].↵
    O’Mahony S.M., Dinan T.G., Cryan J.F.,2017. The gut microbiota as a key regulator of visceral pain. Pain 158 Suppl 1:S19–S28. doi: 10.1097/j.pain.0000000000000779.
    OpenUrlCrossRef
  54. [54].↵
    Oh D.K., Na H.S., Jhun J.Y., Lee J.S., Um I.G., Lee S.Y., Park M.S., Cho M.L., Park S.H.,2023. Bifidobacterium longum BORI inhibits pain behavior and chondrocyte death, and attenuates osteoarthritis progression. PLoS One 18(6):e0286456. doi: 10.1371/journal.pone.0286456.
    OpenUrlCrossRef
  55. [55].↵
    Ostrowsky K., Magnin M., Ryvlin P., Isnard J., Guenot M., Mauguiere F., 2002.Representation of pain and somatic sensation in the human insula: A study of responses to direct electrical cortical stimulation. Cereb Cortex 12(4):376–85. doi: 10.1093/cercor/12.4.376.
    OpenUrlCrossRefPubMedWeb of Science
  56. [56].↵
    Palmer T.M., Sterne J.A.C., Harbord R.M., Lawlor D.A., Sheehan N.A., Meng S., Granell R., Smith G.D., Didelez V., 2011. Instrumental variable estimation of causal risk ratios and causal odds ratios in Mendelian randomization analyses. Am. J. Epidemiol. 173, 1392–1403.
    OpenUrlCrossRefPubMedWeb of Science
  57. [57].↵
    Pang J.C., Aquino K.M., Oldehinkel M., Robinson P.A., Fulcher B.D., Breakspear M., Fornito A., 2023.Geometric constraints on human brain function.Nature 618(7965):566–574. doi: 10.1038/s41586-023-06098-1.
    OpenUrlCrossRef
  58. [58].↵
    Ploghaus A., Tracey I., Gati J.S., Clare S., Menon R.S., Matthews P.M., Rawlins J.N., 1999. Dissociating pain from its anticipation in the human brain. Science 284(5422):1979–81. doi: 10.1126/science.284.5422.1979.
    OpenUrlAbstract/FREE Full Text
  59. [59].↵
    Pratt C., Campbell M.D., 2020. The Effect of Bifidobacterium on Reducing Symptomatic Abdominal Pain in Patients with Irritable Bowel Syndrome: A Systematic Review. Probiotics Antimicrob Proteins 12(3):834–839. doi: 10.1007/s12602-019-09609-7.
    OpenUrlCrossRef
  60. [60].↵
    Qin Y.W., Havulinna A.S., Liu Y., Jousilahti P., Ritchie S.C., Tokolyi A., Sanders J.G., Valsta L., Brożyńska M., Zhu Q.Y., Tripathi A., Vázquez-Baeza Y., Loomba R., Cheng S.S., Jain M., Niiranen T., Lahti L., Knight R., Salomaa V., Inouye M., Méric G., 2022.Combined effects of host genetics and diet on human gut microbiota and incident disease in a single population cohort. Nat Genet 54(2):134–142. doi: 10.1038/s41588-021-00991-z.
    OpenUrlCrossRef
  61. [61].↵
    Ramakrishna C., Corleto J., Ruegger P.M., Logan G.D., Peacock B.B., Mendonca S., Yamaki S., Adamson T., Ermel R., McKemy D., Borneman J., Cantin E.M., 2019. Dominant role of the gut microbiota in chemotherapy induced neuropathic pain. Sci Rep 9(1):20324. doi: 10.1038/s41598-019-56832-x.
    OpenUrlCrossRef
  62. [62].↵
    Ran Guo, Li-Hua Chen, Chungen Xing, Tong Liu, 2019. Pain regulation by gut microbiota: molecular mechanisms and therapeutic potential. Br J Anaesth 123 (5): 637–654. doi: 10.1016/j.bja.2019.07.026.
    OpenUrlCrossRef
  63. [63].↵
    Ren M., Lotfipour S., 2020. The role of the gut microbiome in opioid use, Behav Pharmacol 31(2&3):113–121. doi: 10.1097/FBP.0000000000000538.
    OpenUrlCrossRef
  64. [64].↵
    Rodriguez-Raecke R., Niemeier A., Ihle K., Ruether W., May A., 2009. Brain Gray Matter Decrease in Chronic Pain Is the Consequence and Not the Cause of Pain. J Neurosci 29(44):13746–50. doi: 10.1523/JNEUROSCI.3687-09.2009
    OpenUrlAbstract/FREE Full Text
  65. [65].↵
    Rühlemann DM.C., Hermes B.M., Bang C., Doms S., Moitinho-Silva L., Thingholm L.B., Frost F., Degenhardt F., Wittig M., Kässens DJ., Weiss F.U., Peters A., Neuhaus K., Völker U., Völzke H., Homuth G., Weiss S., Grallert H., Laudes M., Lieb W., Haller D., Lerch DM.M., Baines DJ.F., Franke A., 2021. Genome-wide association study in 8,956 German individuals identifies influence of ABO histo-blood groups on gut microbiome. Nat Genet 53(2):147–155. doi: 10.1038/s41588-020-00747-1.
    OpenUrlCrossRef
  66. [66].↵
    Russo R., Cristiano C., Avagliano C., Caro C.D., Rana G.L., Raso G.M., Canani R.B., Meli R., Calignano A., 2018. Gutebrain axis: role of lipids in the regulation of inflammation, pain and CNS diseases. Curr Med Chem 25(32):3930–3952. doi: 10.2174/0929867324666170216113756.
    OpenUrlCrossRef
  67. [67].↵
    Sanderson E., 2021. Multivariable Mendelian randomization and mediation. Cold Spring Harb Perspect Med 11(2):a038984. doi:10.1101/cshperspect.a038984.
    OpenUrlAbstract/FREE Full Text
  68. [68].↵
    Sawamoto N., Honda M., Okada T., Hanakawa T., Kanda M., Fukuyama H., Konishi J., Shibasaki H.,2000. Expectation of painenhances responses to nonpainful somatosensory stimulation in the anterior cingulate cortex and parietal operculum/posterior insula: An event-related functional magnetic resonance imaging study. J Neurosci 20(19):7438–45. doi: 10.1523/JNEUROSCI.20-19-07438.2000.
    OpenUrlAbstract/FREE Full Text
  69. [69].↵
    Schermuly I., Fellgiebel A., Wagner S., Yakushev I., Stoeter P., Schmitt R., Knickenberg R.J., Bleichner F., Beutel M.E., 2010. Association between cingulum bundle structure and cognitive performance: an observational study in major depression.Eur Psychiatry 25(6):355–60. doi: 10.1016/j.eurpsy.2010.05.001.
    OpenUrlCrossRefPubMedWeb of Science
  70. [70].↵
    Schmidt-Wilcke T., Ganssbauer S., Neuner T., Bogdahn U., May A.,2008. Subtle grey matter changes between migraine patients and healthy controls. Cephalalgia 28(1):1–4. doi: 10.1111/j.1468-2982.2007.01428.x.
    OpenUrlCrossRefPubMedWeb of Science
  71. [71].↵
    Seminowicz D.A., Labus J.S., Bueller J.A., Tillisch K., Naliboff B.D., Bushnell M.C., Mayer E.A.,2010. Regional Gray Matter Density Changes in Brains of Patients With Irritable Bowel Syndrome. Gastroenterology 139(1):48–57.e2. doi: 10.1053/j.gastro.2010.03.049.
    OpenUrlCrossRefPubMedWeb of Science
  72. [72].
    Seminowicz D.A., Moayedi M.,2017. The Dorsolateral Prefrontal Cortex in Acute and Chronic Pain. J Pain 18(9):1027–1035. doi: 10.1016/j.jpain.2017.03.008.
    OpenUrlCrossRefPubMed
  73. [73].↵
    Shiers S., Price T.J.,2020.Molecular, circuit, and anatomical changes in the prefrontal cortex in chronic pain. Pain. 161(8):1726–1729. doi: 10.1097/j.pain.000000000000189745.
    OpenUrlCrossRef
  74. [74].↵
    Silvestre M.P., Rodrigues A.M., Canhão H., Marques C., Teixeira D., Calhau C., Branco J., 2020. Cross-talk between diet-associated dysbiosis and hand osteoarthritis. Nutrients 12(11):3469. doi: 10.3390/nu12113469.
    OpenUrlCrossRef
  75. [75].
    Smith S.M., Douaud G., Chen W., Hanayik T., AlfaroAlmagro F., Sharp K., Elliott L.T., 2021. An expanded set of genome-wide association studies of brain imaging phenotypes in UK Biobank.Nat Neurosci 24(5):737–745. doi: 10.1038/s41593-021-00826-4.
    OpenUrlCrossRef
  76. [76].↵
    Sommer F., Bäckhed F., 2013. The gut microbiota--masters of host development and physiology. Nat Rev Microbiol 11(4):227–38. doi: 10.1038/nrmicro2974.
    OpenUrlCrossRefPubMed
  77. [80].
    Storey J.D.,2003. The positive false discovery rate: A Bayesian interpretation and the q-value. Ann. Stat. 31(6):2013–2035.DOI: 10.1214/aos/1074290335.
    OpenUrlCrossRefWeb of Science
  78. [77].↵
    Su M.C., Tang Y.D., Kong W.S., Zhang S.Y., Zhu T, 2023. Genetically supported causality between gut microbiota, gut metabolites and low back pain: a two-sample Mendelian randomization study. Front Microbiol 14:1157451. doi: 10.3389/fmicb.2023.1157451.
    OpenUrlCrossRef
  79. [78].
    Theodorou V., Ait Belgnaoui A., Agostini S., Eutamene H.,2014. Effect of commensals and probiotics on visceral sensitivity and pain in irritable bowel syndrome. Gut Microbe 5(3):430–436. doi: 10.4161/gmic.29796.
    OpenUrlCrossRefPubMed
  80. [79].↵
    Thomas A.G., Koumellis P., Dineen R.A., 2011.The Fornix in Health and Disease:An Imaging Review.Radiographics. 31(4):1107–21. doi: 10.1148/rg.314105729.
    OpenUrlCrossRefPubMedWeb of Science
  81. [80].↵
    Timmers I., Biggs E.E., Bruckert L., Tremblay-McGaw A.G., Zhang H., Borsook D., Simons L.E.,2024. Probing white matter microstructure in youth with chronic pain and its relation to catastrophizing using neurite orientation dispersion and density imaging. Pain 8. doi: 10.1097/j.pain.0000000000003269.
    OpenUrlCrossRef
  82. [81].↵
    Treede R.D, Rief W., Barke A., Aziz Q., Bennett M.I., Benoliel R., Cohen M., Evers S.,
  83. [82].↵
    Finnerup N.B., First M.B., Giamberardino M.A., Kaasa S., Kosek E., Lavand homme P., Nicholas M., Perrot S., Scholz J., Schug S., Smith B.H., Svensson P., Vlaeyen J.W.S., Wang Shuu-Jiun, 2015. A classification of chronic pain for ICD-11. Pain 156(6):1003–1007. doi: 10.1097/j.pain.0000000000000160.
    OpenUrlCrossRefPubMed
  84. [83].↵
    Vachon-Presseau E., Centeno M.V., Ren W., Berger S.E., Tetreault P., Ghantous M., Baria A., Farmer M., Baliki M.N., Schnitzer T.J., Apkarian A.V., 2016.The Emotional Brain as a Predictor and Amplifier of Chronic Pain.J Dent Res 95(6):605–12. doi: 10.1177/0022034516638027.
    OpenUrlCrossRefPubMed
  85. [84].↵
    Verbanck M., Chen C.Y., Neale B., D R., 2018. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 50 (5), 693–698.
    OpenUrlCrossRefPubMed
  86. [85].↵
    Wang S., Zang M., Yang X., Lv L., Chen L., Cui J., Liu Y., Xia Y., Zhou N., Yang Z., Li Y., Shi B., 2023. Gut microbiome in men with chronic prostatitis/chronic pelvic pain syndrome: profiling and its predictive significance. World J Urol 41(11):3019–3026. doi: 10.1007/s00345-023-04587-6.
    OpenUrlCrossRef
  87. [86].
    Wang X., Wu Y., Liu Y., Chen F., Chen S., Zhang F., Li S., Wang C., Gong Y., Huang R., Hu M., Ning Y., Zhao H., Guo X., 2023. Altered gut microbiome profile in patients with knee osteoarthritis. Front Microbiol 14:1153424. doi: 10.3389/fmicb.2023.1153424.
    OpenUrlCrossRef
  88. [87].↵
    Wu R., Wang D., Cheng L., Su R., Li B., Fan C., Gao C., Wang C., 2024. Impaired immune tolerance mediated by reduced Tfr cells in rheumatoid arthritis linked to gut microbiota dysbiosis and altered metabolites. Arthritis Res Ther 26(1):21. doi: 10.1186/s13075-023-03260-y.
    OpenUrlCrossRef
  89. [88].↵
    Wu Z., Liu H., Yan E., Zhang X., Wang Y., Huang C., He T., Miao L., Yang L., Jiang R., Qi C., Liu C., Wang D., Yang C.,2023. Desulfovibrio confers resilience to the comorbidity of pain and anxiety in a mouse model of chronic inflammatory pain. Psychopharmacology (Berl) 240(1):87–100. doi: 10.1007/s00213-022-06277-4.
    OpenUrlCrossRef
  90. [89].↵
    Xiao X., Ding M., Zhang Y.Q.,2021. Role of the Anterior Cingulate Cortex in Translational Pain Research.Neurosci Bull 37(3):405–422. doi: 10.1007/s12264-020-00615-2.
    OpenUrlCrossRefPubMed
  91. [90].↵
    Yang M.Y., Wan X.J., Zheng H.S., Xu K., Xie J.L., Yu H., Wang J.C., Xu P., 2023. No Evidence of a Genetic Causal Relationship between Ankylosing Spondylitis and Gut Microbiota: A Two-Sample Mendelian Randomization Study.Nutrients 15(4):1057. doi: 10.3390/nu15041057.
    OpenUrlCrossRef
  92. [91].
    Yang S., Chang M.C., 2019.Chronic Pain: Structural and Functional Changes in Brain Structures and Associated Negative Affective State.Int J Mol Sci, 20(13):3130. doi: 10.3390/ijms20133130.
    OpenUrlCrossRefPubMed
  93. [92].↵
    Ye T., Shao J., Kang H., 2021. Debiased Inverse-Variance Weighted Estimator in Two-Sample Summary-Data Mendelian Randomization. The Annals of Statistics, 49(4), 2079–2100. https://arxiv.org/abs/1911.09802.
    OpenUrl
  94. [93].↵
    Yin Q.Y., Zhu L., 2024. Does co-localization analysis reinforce the results of Mendelian randomization? Brain 147(1):e7–e8. doi: 10.1093/brain/awad295.
    OpenUrlCrossRef
  95. [94].↵
    Yu X.H., Yang Y.Q., Cao R.R., Bo L., Lei S.F., 2021. The causal role of gut microbiota in development of osteoarthritis.Osteoarthritis Cartilage 29(12):1741–1750. doi: 10.1016/j.joca.2021.08.003.
    OpenUrlCrossRef
  96. [95].↵
    Zhao J., Ming J., Hu X.H., Chen G., Liu J., Yang C., 2020. Bayesian weighted Mendelian randomization for causal inference based on summary statistics. Bioinformatics 36(5):1501–1508. doi:10.1093/bioinformatics/btz749.
    OpenUrlCrossRef
  97. [96].↵
    Zhao Q., Wang J., Hemani G., Bowden J., Small D.S., 2020. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann Statist 48(3):1742–1769.https://arxiv.org/abs/1801.09652.
    OpenUrl
  98. [97].
    Zhou R., Zhang L., Sun Y., Yan J., Jiang H., 2023. Causal Associations between Dietary Habits and Chronic Pain: A Two-Sample Mendelian Randomization Study. Nutrients 15(17):3709. doi: 10.3390/nu15173709.
    OpenUrlCrossRef
Back to top
Posted September 30, 2024.
Download PDF
Data/Code
Email
The causal connections between gut microbiota and chronic pain and pain phenotypes through the gut-brain axis: A Mendelian randomization study and mediation analysis
Jia Chen, Xuemei Luo, Luyao Wang, Yiping Luo, Xia Deng, Yuan Tian, Qiao Lu, Genwen Sun, Michael Maes, Xu Zhang
medRxiv 2024.09.28.24314553; doi: https://doi.org/10.1101/2024.09.28.24314553
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area