Genetically predicted causal effects of gut microbiota on spinal pain: a two-sample Mendelian randomization analysis

Background Observational studies have hinted at a correlation between the gut microbiota and spinal pain (SP). However, the impact of the gut microbiota on SP remains inconclusive. Methods In this study, we employed a two-sample Mendelian randomization (MR) analysis to explore the causal relationship between the gut microbiota and SP, encompassing neck pain (NP), thoracic spine pain (TSP), low back pain (LBP), and back pain (BP). The compiled gut microbiota data originated from a genome-wide association study (GWAS) conducted by the MiBioGen consortium (n = 18,340). Summary data for NP were sourced from the UK Biobank, TSP from the FinnGen Biobank, and LBP from both the UK Biobank and FinnGen Biobank. Summary data for BP were obtained from the UK Biobank. The primary analytical approach for assessing causal relationships was the Inverse Variance Weighted (IVW) method, supplemented by various sensitivity analyses to ensure result robustness. Results The IVW analysis unveiled 37 bacterial genera with a potential causal relationship to SP. After Benjamini-Hochberg corrected test, four bacterial genera emerged with a strong causal relationship to SP. Specifically, Oxalobacter (OR: 1.143, 95% CI 1.061–1.232, P = 0.0004) and Tyzzerella 3 (OR: 1.145, 95% CI 1.059–1.238, P = 0.0007) were identified as risk factors for LBP, while Ruminococcaceae UCG011 (OR: 0.859, 95% CI 0.791–0.932, P = 0.0003) was marked as a protective factor for LBP, and Olsenella (OR: 0.893, 95% CI 0.839–0.951, P = 0.0004) was recognized as a protective factor for low back pain or/and sciatica. No significant heterogeneity or horizontal pleiotropy was observed through alternative testing methods. Conclusion This study establishes a causal relationship between the gut microbiota and SP, shedding light on the “gut-spine” axis. These findings offer novel perspectives for understanding the etiology of SP and provide a theoretical foundation for potential interventions targeting the gut microbiota to prevent and treat SP.


Introduction
Spinal pain (SP) manifests as a prevalent symptom in various spine-related disorders, encompassing neck pain (NP) and back pain (BP), further categorized into thoracic spine pain (TSP) and low back pain (LBP).The prevalence in the general population is noteworthy, with NP affecting approximately 44%, TSP 15%, and LBP 56% (Manchikanti et al., 2004).The impact of SP on daily life and occupational functions surpasses that of other musculoskeletal disorders due to the spine's role as the central axis and core support of the human skeleton (Haldeman et al., 2012).The 2015 Global Burden of Disease study underscored SP as a paramount contributor to disability globally, especially affecting individuals aged 25-64 (GBD 2015 Disease andInjury Incidence andPrevalence Collaborators, 2016).Consequently, healthcare costs for SP treatment soared, reaching an estimated $134.5 billion in the United States in 2020 (de Luca et al., 2023).The pressing need for effective prevention and treatment strategies for SP is a critical concern in the realm of global public health.
Recent studies propose a potential relationship between the gut microbiota and the development of SP.Patients with NP exhibited higher levels of Bacteroidetes and Proteobacteria and significantly lower levels of Firmicutes compared to healthy controls (Qin, 2022).Similarly, individuals with BP showed a distinct gut microbiota composition, characterized by elevated concentrations of Adlercreutzia, Roseburia, and Uncl.Christensenellaceae (Dekker Nitert et al., 2020).These findings underscore a compelling relationship between gut microbiota and SP.While degenerative spinal diseases traditionally bear the blame for SP, a notable shift toward a younger demographic experiencing SP suggests alternative causative factors beyond structural spinal degeneration (Lei et al., 2015).Consequently, gut microbiota dysbiosis emerges as a significant potential contributor to SP.
In light of these observations, MR offers a promising avenue for exploring causal relationships.MR, first proposed by Smith and Ebrahim (2003), utilizes genetic variation as an instrumental variable, enabling the elimination of confounding bias and reverse causation (Emdin et al., 2017;Soremekun et al., 2022).In this study, we collected data on gut microbiota and SP from the genomewide association study (GWAS) database, employing MR analysis to discern the causal relationship between them.The study aims to provide innovative insights for the prevention and treatment of SP by unraveling the intricate interplay between gut microbiota dynamics and SP.

Study design
Mendelian Randomization, a pivotal method for evaluating causal relationships between exposure and outcome, stands as the cornerstone of our investigative framework.In quest of discerning whether gut microbiota assumes a protective or facilitative role in SP development, our study strategically designates gut microbiota as the exposure variable, with four SP types-NP, TSP, LBP, and BP-comprising the outcome variables.According to Davies et al. (2018), the implementation of a two-sample MR protocol necessitated strict adherence to three fundamental assumptions.First, the instrumental variables, genetic variation, should exhibit a robust relationship with gut microbiota.Second, genetic variation must remain unlinked to confounding factors, excluding those pertaining to gut microbiota and the four SP types.Third, the effects of genetic variation on the four SP types must be channeled through a singular pathway-namely, gut microbiota-with the exclusion of any alternative pathways.In Figure 1, the specific flowchart detailing the MR analysis is presented.

Exposure data
The summary statistics for gut microbiota utilized in this investigation originated from a GWAS conducted by the MiBioGen consortium in 2021 (Kurilshikov et al., 2021).The study encompassed a substantial cohort of 18,340 individuals, with over 70% representing European ancestry (n = 13,266).Within this cohort, participants were stratified into 24 distinct groups based on factors such as background, sex, and age.The composition of their fecal microbiota was meticulously Three assumptions and flowchart for Mendelian randomization study.
examined through the application of 16S rRNA gene sequencing.This comprehensive analysis led to the identification of 211 taxa, spanning across 131 genera, 35 families, 20 orders, 16 classes, and 9 phyla.Recognizing the significance of genuslevel precision in characterizing the microbiome, we opted to exclusively retain data at the genus level for subsequent phases of our study and detailed analysis.This strategic decision ensures a more granular and accurate exploration of each potentially pathogenic microbiome.

Outcome data
All GWAS summary statistics pertaining to SP were meticulously sourced from the EU Open GWAS platform. 1This comprehensive repository comprises data packages for the four primary types of SP: NP, TSP, LBP, and BP.After a thorough review of pertinent datasets, eight were judiciously selected for our analysis.Specifically, two GWAS summary statistics for NP, considering disease durations of 1 month and 3+ months, were incorporated into our study, both sourced from the UK Biobank (ukb-e).The GWAS summary statistics for TSP were obtained from the FinnGen Biobank.Similarly, the GWAS summary statistics for LBP were procured from the FinnGen Biobank.Recognizing that patients with LBP often exhibit radiating pain in the lower limbs (Koes et al., 2006), we expanded our dataset to include two GWAS summary statistics: one for low back pain with radicular pain from the UK Biobank (ukb-e) and another for low back pain with sciatica from the FinnGen Biobank.For the BP category, two GWAS summary statistics for BP, comparing disease durations of 1 month versus 3+ months, were selected from the UK Biobank (ukb-e).For a 1 https://gwas.mrcieu.ac.uk/ detailed breakdown of outcome summary statistics, please refer to Supplementary Table 1.

Instrumental variable (IV)
In order to ensure the precision and robustness of the study outcomes, we incorporated single nucleotide polymorphisms (SNPs) demonstrating robust relationships with the exposure factors.While the conventional threshold for controlling the false positive rate is typically set at P < 5e−08, this study opted for a more inclusive threshold of P < 1e−05 (Xie et al., 2023).This adjustment was made to augment the pool of screened SNPs, thereby enhancing result credibility, given the limited number of SNPs obtained under the more stringent condition.Subsequently, to acquire mutually independent instrumental variables, SNPs within a 10,000 kb clumping window were excluded to mitigate Linkage Disequilibrium (LD), specifically retaining SNPs with an R 2 < 0.001.Concurrently, F-statistics of the SNPs were computed, leading to the exclusion of weak instrumental variables characterized by F < 10 ( Burgess and Thompson, 2011).Finally, a meticulous harmonization process was applied to align the effector alleles of the exposure and outcome SNPs.This involved the exclusion of alleles that were incompatible (e.g., A/C paired with A/G) or exhibited intermediate allele frequencies.

Statistical analysis
The causal relationship between the gut microbiota and SP was primarily analyzed using the Inverse Variance Weighted (IVW) method.A potential causal relationship between the gut microbiota and SP was considered when the IVW result yielded a P < 0.05.Following multiple testing correction of the IVW results using the Benjamini-Hochberg False Discovery Rate (FDR) method, gut Hong et al. 10.3389/fmicb.2024.1357303microbiota with P < 0.05 were deemed to exhibit a significant causal relationship with SP.However, the validity of this conclusion relies on the absence of horizontal pleiotropy and heterogeneity (Burgess et al., 2016).Therefore, various additional methods were incorporated to scrutinize the potential presence or absence of horizontal pleiotropy and heterogeneity.
Initially, Cochrane's Q test was employed to assess heterogeneity in the IVs.A significance level of P < 0.05 denoted heterogeneity, necessitating the use of the random-effects IVW model; conversely, the fixed-effects IVW model was applied (Greco et al., 2015).In the IVW regression, where the intercept term is not considered, results become inaccurate in the presence of horizontal pleiotropy (Burgess et al., 2017).To address this, the MR-Egger regression method and the Weighted Median (WM) method were employed to correct for horizontal pleiotropy.MR-Egger assumes that the IVs of horizontal pleiotropy is greater than 50%, while WM assumes it to be less than 50%, both re-estimating results accordingly (Bowden et al., 2016).Both methods complemented the IVW method, correcting for bias introduced by horizontal pleiotropy and enhancing result credibility (Bowden et al., 2015;Slob and Burgess, 2020).MR-Egger and MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) were employed to test horizontal pleiotropy.An MR-Egger P-value of less than 0.05 indicated a nonzero intercept, signifying the presence of horizontal pleiotropy.In addition to testing horizontal pleiotropy, MR-PRESSO identified and corrected for outliers, thus improving the accuracy of the analysis results.Finally, a leave-one-out sensitivity analysis was conducted to assess the stability of the MR results and detect the presence of outliers in the SNPs.
All the aforementioned statistical analyses were performed using the R program (version 4.3.1),with the "VariantAnnotation, " "mrcieu/gwasglue, " "mrcieu/ieugwasr, " and "MRCIEU/TwoSampleMR" packages loaded and executed.The specific MR analysis codes are detailed in the Supplementary material.

SNP selection
From a pool of 119 bacterial genera, we meticulously screened and identified 1531 SNPs to serve as IVs in our analysis.Crucially, all of these SNPs exhibited F-statistics surpassing 10, with values ranging from 14.58 to 88.42.This range signifies a robust F-statistic profile, indicating the absence of bias arising from weak IVs.This stringent selection process ensures the reliability and strength of the IVs utilized in our study.

FIGURE 2
Forest plot of causality between gut microbiota and risk of neck pain.

Sensitivity analysis and Benjamini-Hochberg corrected test
The results of our MR analysis were subjected to a comprehensive set of sensitivity tests to ensure the robustness Frontiers in Microbiology 10 frontiersin.orgForest plot of causality between gut microbiota and risk of back pain.
and reliability of the findings.We employed the IVW test, MR-Egger's intercept test, and the MR-PRESSO global test (Supplementary Table 10).Importantly, the vast majority of IVW test, relying on Cochran's Q test, demonstrated a lack of heterogeneity (P > 0.05), indicating consistency across our analyses.Furthermore, neither MR-Egger's intercept test nor the MR-PRESSO global test detected evidence of horizontal pleiotropy (P > 0.05), reinforcing the reliability of our findings.Scatterplots illustrating the outcomes of the five MR analysis methods have been presented in Supplementary Figures 1-8.In addition, a leave-one-out analysis, as illustrated in Supplementary Figures 9-16, demonstrated that the exclusion of any single SNP did not significantly alter the overall results, confirming the stability and reliability of our findings.To control for false positives, we applied the Benjamini-Hochberg corrected test.Corrected for multiple testing, Oxalobacter (OR: 1.143, 95% CI 1.061-1.232,P = 0.0004) and Tyzzerella 3 (OR: 1.145, 95% CI 1.059-1.238,P = 0.0007) were identified as risk factors for LBP, while Ruminococcaceae UCG011 (OR: 0.859, 95% CI 0.791-0.932,P = 0.0003) was identified as a protective factor against LBP.Additionally, Ruminococcus gnavus group (OR: 0.998, 95% CI 0.996-0.999,P = 0.0004) was revealed as a protective factor specifically against lumbar and other intervertebral disk disorders with radiculopathy, and Olsenella (OR: 0.893, 95% CI 0.839-0.951,P = 0.0004) emerged as a protective factor against low back pain or/and sciatica.

Discussion
To the best of our knowledge, this study stands as the inaugural and thorough investigation employing publicly available GWAS databases to scrutinize the causal relationship between the gut microbiota and four prevalent types of SP-specifically, NP, TSP, LBP, and BP.Through a meticulous analysis of the data, we observed that all IVW results within the dataset for ' lumbar and other intervertebral disk disorders with radiculopathy ' were close to 1, with narrow 95% confidence intervals.This suggests the absence of a significant causal effect between the gut microbiota and this particular type of LBP.Consequently, we opted to exclude the data related to this specific type of LBP.Following a stringent selection process, we successfully identified a total of 37 bacterial genera exhibiting potential causal relationships with SP.Subsequently, through the application of the Benjamini-Hochberg corrected test, we corroborated the strong causal relationship of four bacterial genera with SP.The discoveries emerging from this study serve as an invaluable reference for subsequent research endeavors aimed at preventing and managing SP through precise interventions in the gut microbiota.
The study demonstrates that over 90% of cases involving SP lack a definitive diagnosis related to identifiable causes such as fractures, infections, tumors, etc (Haldeman et al., 2012).Consequently, clinical interventions face challenges in delivering precise and effective targeted treatments, often relying on symptomatic relief measures such as analgesics and soft tissue relaxation (Cheng, 2021).In recent years, growing attention has been directed toward understanding the intricate relationship between the gut microbiota and pain.The gut microbiota, with its array of metabolites, neurotransmitters, by-products, and other mediators, possesses the ability to modulate pain (Guo et al., 2019).In this context, we conducted a comprehensive exploration of the causal relationship between the gut microbiota and SP using Mendelian randomization analysis.Upon subjecting the MR results to the Benjamini-Hochberg corrected test, a strong causal relationship emerged between the gut microbiota and low back pain or/and lower limb pain.However, no significant causal relationship was discerned between the gut microbiota and NP, TSP, or BP.The outcomes of our study align with the observations of Arai et al. (2018), who reported a noteworthy positive relationship between the severity of constipation and the intensity of pain in individuals experiencing low back pain or/and lower limb pain.Notably, the severity of constipation is known to exert a significant influence on the composition of the gut microbiota (Zhu et al., 2014).We observed that Oxalobacter and Tyzzerella 3 act as risk factors for low back pain or/and lower limb pain, while Ruminococcaceae UCG011 and Olsenella exhibit protective effects.Rajasekaran et al. (2020) identified a higher abundance of Oxalobacter in degenerated intervertebral disks.Intervertebral disk degeneration is widely recognized as a crucial factor in the onset and progression of LBP (Meng et al., 2023).Hence, we posit that Oxalobacter might contribute to LBP by potentially inducing intervertebral disk degeneration.Currently, there is limited clinical or experimental research on Tyzzerella 3 concerning LBP, and direct evidence establishing Tyzzerella 3 as a risk factor is lacking.Some studies found that Tyzzerella can produce aromatic amines and is positively correlated with the occurrence of diabetes (Ruuskanen et al., 2022;Sugiyama et al., 2022).Therefore, Tyzzerella 3 might promote the development of LBP by influencing endocrine or endothelial functions.Regarding the protective bacterial genera identified in the development of LBP, relevant studies on changes in the abundance of Ruminococcaceae UCG011 and Olsenella in the feces of LBP patients or animal models are currently unavailable.Ruminococcaceae UCG011 belongs to the Ruminococcaceae family, and evidence indicates a significant reduction in Ruminococcaceae family bacteria in the feces of mice with lumbar disk herniation compared to the control group (Wang et al., 2021).Butyrate produced by Ruminococcaceae family bacteria can promote the synthesis of vitamin D in the body (Murdaca et al., 2021).Studies have found that vitamin D has a regulatory effect on the immune function of the body, as it can inhibit the differentiation of pro-inflammatory cells such as T helper (Th) 17 cells, thereby reducing the secretion of pro-inflammatory factors such as IL-17 (Murdaca et al., 2011;Yamamoto and Jørgensen, 2020).Therefore, we hypothesize that improving local chronic inflammation in the body may be a potential mechanism for the protective effect of Ruminococcaceae UCG011.Additionally, research indicates that Short-chain fatty acids (SCFAs) produced by Olsenella can mediate the immune balance between Th17 cells and regulatory T cells (Tregs) by activating Tregs' activity, thus balancing the levels of pro-inflammatory factors (derived from Th17 cells) and anti-inflammatory factors (derived from Treg cells) in the body (Kim, 2021).This could be a potential mechanism by which Olsenella resists the development of LBP.The involvement of the gut microbiota in the development of LBP is an exceedingly complex process.Therefore, the specific mechanisms by which the four identified positive gut microbiota affect low back pain or/and lower limb pain, require further experimental validation in subsequent studies.
In addition to the four bacterial genera strongly associated with SP previously mentioned, our study identified 33 other bacterial genera exhibiting potential causal relationship with SP.To ensure a comprehensive exploration of gut microbiota potentially linked to SP, we incorporated data from all gut microbiota groups with Benjamini-Hochberg corrected test p-values below 0.05.Among the eight bacterial genera potentially causally related to NP, the majority belonged to the Firmicutes phylum.Qin et al. (2019) reported symptom improvement in cervical spondylosis patients following manual therapy, accompanied by a significant increase in the proportions of Lactobacillus and Bifidobacterium in feces, both belonging to the Firmicutes phylum.
Unfortunately, limited clinical and foundational research exists on the gut microbiota's association with NP, and there is no evidence supporting the relationship between Lachnospiraceae UCG010, Escherichia-Shigella, Methanobrevibacter, Eubacterium nodatum group, Faecalibacterium, Subdoligranulum, Butyrivibrio, Ruminococcaceae UCG010, and NP.Potential protective factors for TSP include Butyrivibrio, Ruminococcaceae UCG011, and Eubacterium brachy group.However, similar to NP, research on the gut microbiota and TSP is scarce, lacking supporting evidence.Interestingly, through a comparative analysis of potential causally related bacterial genera in NP and TSP, Butyrivibrio was identified as a shared potential protective genus.Its production of butyrate, propionate, and other SCFAs inhibiting histone deacetylase activity may be one of the potential mechanisms underlying Butyrivibrio's analgesic effects (Descalzi et al., 2015;Tang et al., 2022).Research on the relationship between gut microbiota and low back pain or/and lower limb pain is more abundant.Yao et al. (2023) found a significant decrease in Lachnospiraceae NK4A136 group, Turicibacter, and Roseburia in the feces of rats in a LBP model compared to the control group.In various LBP animal models, Wang et al. (2021) also found lower abundance of Lachnospiraceae NK4A136 group in rats with lumbar disk herniation compared to the control group.A cross-sectional study revealed a significant reduction in skeletal muscle mass index in patients with LBP compared to normal individuals (Park et al., 2023).Sugimura et al. (2022) found a positive relationship between the abundance of Eisenbergiella and skeletal muscle mass, suggesting a reasonable hypothesis that Eisenbergiella may prevent LBP by increasing lumbar skeletal muscle mass.Consistent with the above findings, our study identified Lachnospiraceae NK4A136 group, Turicibacter, Roseburia, and Eisenbergiella as potential protective factors for LBP.In studies investigating the relationship between BP and gut microbiota, Dekker Nitert et al. (2020) found a higher abundance of the genus Uncl.Christensenellaceae in patients with BP.However, contrary to this conclusion, our results indicate that Christensenellaceae R.7 group is a potential protective factor for BP, a finding supported by Hollister et al. (2020).In summary, this study experimentally validated some of the bacterial genera with potential causal relationships to SP that were screened.This underscores the significance of exploring bacterial genera with potential causal relationships to SP, and the unvalidated genera warrant further investigation.
Whether it is NP, TSP, LBP, or BP, they represent distinct segments of SP.Fundamentally, these pains originate in the spinal column, suggesting a certain degree of interconnectedness among pains occurring in different spinal segments.To dissect this phenomenon, we conducted a comparative analysis of bacterial genera potentially influencing SP across different segments.Our investigation revealed that Ruminococcaceae UCG010 exhibits potential protective properties, common to both NP in the last month and BP in the last month.Similarly, Ruminococcaceae UCG011 emerges as a potential safeguard against TSP, LBP, and the amalgamation of low back pain or/and sciatica.Notably, Butyrivibrio surfaces as a potential protective agent against NP in the last month and TSP.Intriguingly, the Eubacterium nodatum group assumes a potential risk factor role in NP for 3+ months while concurrently acting as a potential protective factor for low back pain or/and sciatica.This nuanced outcome suggests that a singular bacterial genus might exert counteractive regulatory influences on pains manifesting in different spinal segments.In summation, the intricate modulation of pains across different spinal segments by the gut microbiota underscores a complex process.Shared or divergent mechanisms may underscore these phenomena, necessitating further exploration to unveil the underlying intricacies.Furthermore, to scrutinize potential disparities in causally related bacterial genera associated with SP across varying durations, we incorporated datasets reflecting distinct timeframes-1 month and 3 months.Regrettably, datasets for different durations of TSP and low back pain or/and lower limb pain, were unavailable.Consequently, our analysis encompassed datasets for 1 month and 3 months of NP and BP.MR outcomes reveal an absence of shared potential causally related bacterial genera between 1 month and 3 months of NP, a trend also mirrored in the context of BP.This suggests a dynamic nature, signifying that bacterial genera with potential causal relationship with SP may evolve during the course of its progression.Contrary to the findings of Dekker Nitert et al. ( 2020), who observed heightened abundance of Adlercreutzia, Roseburia, and Uncl.Christensenellaceae in patients with 1 month of BP, our MR results challenge this observation.Specifically, Adlercreutzia (P = 0.53) and Roseburia (P = 0.92) do not emerge as potential causally related bacterial genera for 1 month of BP.Instead, the Christensenellaceae R.7 group surfaces as a potential protective bacterial genus for this duration of BP.We posit that discrepancies in results may be attributed to the relatively modest sample size employed in the clinical study.In conclusion, it is imperative to delineate causally related bacterial genera associated with SP at distinct developmental stages.This endeavor holds paramount significance for tailoring precise interventions targeting SP across varying durations.
The robustness of this investigation stems from the application of MR to evaluate the causal relationship between gut microbiota and SP.By leveraging genetic variations as IVs, MR effectively mitigates external confounding factors, thereby bolstering internal validity owing to the inherent randomness and independence of these genetic determinants.Notably, subgroup analyses conducted for distinct durations of spinal pain (i.e., 1 month and 3 months) contribute valuable insights, offering implications for future interventions targeting specific gut microbiota to address SP across varying timeframes.However, our study contends with several limitations.Firstly, the exclusivity of data from European countries prompts the necessity for further experimental validation to ascertain the generalizability of MR results to populations in other regions, such as Asia and Africa.Secondly, the adherence to the conventional significance threshold for SNPs in MR studies (P < 5e−08) restricts the number of selected SNPs, potentially influencing subsequent sensitivity analyses and control for horizontal pleiotropy.Thirdly, the constraint of exposure data to the genus level hampers a more nuanced analysis of the causal relationship between higher taxonomic levels of gut microbiota and SP.Lastly, relying on public databases rather than original sources limits the scope of available datasets, resulting in the omission of subgroup analyses for TSP and low back pain or/and lower limb pain during SP subgroup analyses-a notable constraint in our study.

Conclusion
Through our two-sample Mendelian randomization (MR) analysis, we have successfully identified a comprehensive set of 37 bacterial genera potentially causally related to SP, encompassing NP, TSP, LBP, and BP.Employing the rigorous Benjamini-Hochberg corrected test, we have discerned four bacterial genera that exhibit robust causal relationships with SP.Notably, our analysis has unveiled a potential common regulatory mechanism, suggesting that certain bacterial genera may influence spinal pain across different segments.Furthermore, the dynamic nature of these potential causal relationships has been illuminated, indicating that, as the duration of SP progresses, the involved bacterial genera undergo changes.In conclusion, this study contributes valuable insights into the intricate interplay between the gut microbiota and SP.Such revelations not only serve as a significant reference for understanding these complex connections but also pave the way for future research endeavors aimed at preventing and treating SP through targeted interventions in the gut microbiota.

FIGURE 5
FIGURE 5 The author(s) declare financial support was received for the research, authorship, and/or publication of this article.This work was supported by the "Pioneer" and "Leading Goose" R&D Program of Zhejiang (No. 2022C03123), the General Program of the National Natural Science Foundation of China (81774442 and 82274672), the National Natural Science Foundation of China Youth Fund Project (82305426), Key Discipline Project of Highlevel TCM of National Administration of Traditional Chinese Medicine (GJXK2023-85), the Zhejiang Provincial Natural Science Foundation Project (Q23H270025), the Scientific research project of Zhejiang Provincial Department of Education (Y202351266), and the School-Level Scientific Research Project of Zhejiang Chinese Medicine University (2021RCZXZK03, 2022FSYYZQ13, and 2022GJYY045).

TABLE 2
MR results of causal effects between gut microbiome and thoracic spine pain.
TABLE 3 MR results of causal effects between gut microbiome and low back pain with/without lower limb pain.
FIGURE 4Forest plot of causality between gut microbiota and risk of low back pain with/without lower limb pain.MR, mendelian randomisation; SE, standard error; IVW, inverse variance weighted; OR, odds ratio; CI, confidence interval.
TABLE 4 MR results of causal effects between gut microbiome and back pain.