A multi-omics bidirectional mendelian randomization study and meta-analysis on the causal relationship between gut microbiota, inflammatory proteins, and fibromyalgia ====================================================================================================================================================================== * Mengqi Niu * Jing Li * Victoria Sarafian * Michael Maes ## Abstract **Background** Fibromyalgia (FM) is a chronic disorder characterized by widespread pain and immune dysregulation. Emerging evidence suggests that gut microbiota and inflammatory proteins may contribute to the development of FM. **Objective** The aim of this study was to investigate the causal relationships between gut microbiota, inflammatory proteins (cytokines/chemokines), and FM using bidirectional Mendelian randomization (MR) and meta-analysis approaches. **Methods** MR analyses were conducted using genetic data from European populations, employing methods such as MR-IVW, MR-Egger, and MR-weighted median. Reverse MR was also performed, with FM treated as the exposure. A meta-analysis was conducted to consolidate the findings. **Results** *Ruminococcus gauvreauii* was identified as a risk factor for FM, while *Enterorhabdus*, *Parabacteroides*, *Butyricicoccus*, and *Prevotella 9* were found to be protective. Five inflammatory proteins—C-X-C motif chemokine 5 (CXCL5), S100-A12, Leukemia inhibitory factor receptor (LIFR), Monocyte chemoattractant protein 2 (MCP-2/CCL8), and Tumor necrosis factor (TNF-α)—exhibited protective associations, while Natural killer cell receptor 2B4 (NKCR-2B4/CD244) and Interleukin-12 subunit beta (IL-12β) were associated with an increased risk of FM. **Conclusion** This study highlights the role of gut microbiota and inflammatory proteins (cytokines/chemokines) in the pathogenesis of FM. Through Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, the findings suggest their involvement in immune regulation, inflammatory responses, and viral pathways. These findings provide new insights into potential therapeutic targets for modulating gut health and immune responses, opening new avenues for future research and clinical interventions. Keywords * Fibromyalgia * Gut Microbiota * Inflammatory Proteins * Mendelian Randomization Study * Meta Analysis * Multi-omics ## 1. Introduction Fibromyalgia (FM) is a condition characterized by widespread chronic musculoskeletal pain, stiffness in specific areas, and a heightened sensitivity to pressure, known as tender points in soft tissue (1). Despite pain being its primary and prominent feature, FM is characterized by a complex array of symptoms including fatigue, sleep disturbances, and functional symptoms (medical symptoms that cannot be explained by structural or pathological definitions) (2). Globally, FM is estimated to affect 2.7% of the population, with higher prevalence rates observed in women, individuals over 50 years old, those with lower socioeconomic status, and obese individuals (3). However, due to the lack of consensus on the pathogenic mechanisms underpinning FM, there is currently no universally effective treatment (2). A recent meta-analysis shows that FM is characterized by indicants of peripheral immune activation with changes in immunoregulatory cytokines and pro-inflammatory cytokines (4). FM patients exhibit enhanced circulating inflammatory cytokines and chemokines released by peripheral blood mononuclear cells (PBMCs), indicating activation of both innate and adaptive immunity (5–7). Current research revealed elevated concentrations of M1 macrophage cytokines, such as IL-6, IL-8, IL-1β, and TNF-α, in the serum of FM patients (4, 5, 8–10). Immune cells, including mast cells, monocytes, and neutrophils, which mediate the inflammatory process, may also contribute to the immune pathogenesis of FM (11). In patients with major depression, acute COVID-19 infection, Long-COVID, schizophrenia and end-stage renal disease, FM symptoms are associated with immune activation as indicated by increased levels of pro-inflammatory mediators (12–16). Interestingly, plasma levels of soluble gp130 (sgp130), the natural inhibitor of IL-6 transsignalling, is significantly decreased in FM (17), which together with increased IL-6 levels suggest that the immune activation in FM is mainly driven by increased IL-6 signaling (4). Additionally, in animal studies, resident macrophages in muscle tissue have been shown to play a role in the development of chronic widespread muscle pain (18). Increasing evidence suggests that neurogenic inflammatory processes occurring in peripheral tissues, the spinal cord, and the brain contribute to the pathophysiology of FM (5, 19, 20). All of these translate into the clinical symptoms reported by FM patients, such as swelling and sensory dullness, which also impact central symptoms including cognitive changes and fatigue. Furthermore, physiological mechanisms related to stress and emotions are considered upstream driving factors of neurogenic inflammation in FM (7, 21). In 2007, it was published that FM symptoms of chronic fatigue syndrome are associated with increased translocation of Gram-negative gut commensal bacteria, indicating leaky gut (22). In schizophrenia patients, FM symptoms are strongly associated with breakdown of the paracellular pathways (leaky gut) and increased bacterial translocation (13). Other findings suggested that gut dysbiosis may play a role in FM. Thus, studies by Minerbi et al. (23) and Minerbi and Fitzcharles et al. (24) reported associations between gut microbiota and FM. In addition, a systematic review showed that intestinal microbiota is involved in the pathophysiology of FM (25). However, until now it has remained unclear whether there are causal relationships between immune activation, the gut-immune axis, and FM. Recently developed Mendelian randomization (MR) designs have been used to study whether specific exposures have causal effects on specific outcomes (26). As an analytical method, MR designs utilize lineage genetic variations (single nucleotide polymorphisms, SNPs) as instrumental variables for the target exposure (27). Based on the random assortment during meiosis and the fixed allocation of genetic variations at conception, MR studies are less susceptible to confounding factors and reverse causation biases. Additionally, the use of summary statistics from genome-wide association studies (GWAS) in two-sample MR designs significantly enhances the statistical power of causal inference (28). Therefore, in this study, we use summary statistics from GWAS to conduct MR research. The a priori hypothesis is to discover clear multi-omics causal relationships between FM and immune activation as assessed with cytokines/chemokines and changes in the gut microbiota indicating an inflammatory milieu or leaky gut. ## 2. Materials and Methods ### 2.1. Study Design In our analysis, we used publicly available data from shared repositories, eliminating the need for additional ethical approval for this study. The MR analysis must meet three assumptions: (1) SNPs are associated with exposure; (2) SNPs are independent of confounders in the exposure-outcome relationship (independence assumption); (3) SNPs influence the outcome solely through exposure (29). **Figure 1** illustrates the workflow of this study.  [Figure 1](http://medrxiv.org/content/early/2024/09/14/2024.09.13.24313599/F1) Figure 1 The flowchart of bidirectional mendelian randomization study and meta-analysis explore the causal relationship between gut microbiota, inflammatory proteins, and fibromyalgia. PQTL: protein quantitative trait loci ### 2.2. Data Sources This study utilized three types of exposure data and two types of outcome data. To minimize bias due to population heterogeneity, only aggregated data from European populations were used. SNPs related to the human gut microbiome composition were obtained from the latest data from the MiBioGen international consortium ([https://mibiogen.gcc.rug.nl/](https://mibiogen.gcc.rug.nl/)). This GWAS included 211 microbial taxa classified at six levels, from phylum to species (30). We sourced data on 91 inflammatory protein types from the study by Zhao et al (31). Summary statistics for FM outcomes were sourced from the FinnGen consortium (released in 2023, [https://r10.finngen.fi/](https://r10.finngen.fi/)) and the study by Moscati et al (32). ### 2.3. Genetic Instrument Selection (1) Based on previous studies, a significance threshold (p < 1×10[[) was applied for gut microbiota (33), while a stricter p-value (p < 5×10[[) was used for inflammatory proteins (34) to ensure a moderate number of SNPs for all traits. (2) To ensure independence of the instruments used for exposure, we excluded SNPs in linkage disequilibrium (LD) (r² < 0.001, clumping window = 10,000 kb) (35). (3) Palindromic SNPs with intermediate minor allele frequencies (MAF) were excluded (36). (4) To determine the existence of weak instrument bias, an F-statistic was computed to measure the power of the instrumental variable. An F-statistic above 10 indicated a low likelihood of weak instrument bias (37). ### 2.4. Two-Sample Mendelian Randomization This study primarily relied on random effects inverse variance weighting (IVW) to assess the causal effect of gut microbiota, and inflammatory proteins on FM. The IVW estimates were derived by aggregating the Wald ratios of all genetic variants. IVW is based on the assumption that all SNPs do not exhibit horizontal pleiotropy, thus providing the most accurate causal effect estimate. To ensure robustness, we also employed MR-Egger and weighted median (WM) methods for supplementary analysis (28, 38). Sensitivity analyses are essential for examining horizontal pleiotropy and heterogeneity, which can significantly impact MR estimates. In this study, we used three methods to detect and correct these issues: Cochran’s Q test, MR-Egger intercept test, and MR-PRESSO (Pleiotropy RESidual Sum and Outlier) analysis. A p-value < 0.05 from Cochran’s Q test indicated heterogeneity in the results (39). The MR-Egger intercept was used to detect directional pleiotropy and bias due to invalid instruments (38). Finally, MR-PRESSO was used to identify and correct for potential pleiotropic outliers (40). To assess the robustness of the results, we employed the leave-one-out method, sequentially excluding each SNP to determine if any single SNP drove the estimates (38). We stringently screened the gut microbiota and cytokines/chemokines with potential causal relationships to FM based on several criteria: (1) Significant p-value (p < 0.05) from the primary MR IVW analysis; (2) Consistent direction and magnitude across three MR methods; (3) No significant heterogeneity or horizontal pleiotropy in MR results. ### 2.5. Meta-analysis and genetic directionality assessment To comprehensively estimate the robustness of gut microbiota and inflammatory proteins identified by the aforementioned criteria, we used the FinnGen consortium’s summary statistics for FM as the preliminary analysis. We replicated the IVW analysis in another FM cohort (study by Moscati et al.). The meta-analysis of the results from the two MR analyses identified gut microbiota and immune proteins causally related to FM. Additionally, to rule out bias from reverse causality and further evaluate the robustness of the causal relationship between gut microbiota and immune proteins, and FM, we conducted reverse MR analysis within the preliminary analysis cohort (FinnGen consortium) for significant results (P < 0.05) from the meta-analysis. This provided final evidence for causal relationships. Statistical analyses were performed using R software (version 4.3.2) with the packages MendelianRandomization, MRPRESSO, TwoSampleMR, Meta, among others. For binary outcomes, results were reported as odds ratios (OR) with corresponding 95% confidence intervals (CI) to quantify the strength of causal relationships. For continuous outcomes, results were reported as β coefficients, standard errors (SE), and corresponding 95% CI. ### 2.6. Protein-protein interaction (PPI) annotation and enrichment To better explore the interactions between proteins and the pathogenesis, we performed PPI annotation and enrichment analysis using STRING ([https://cn.string-db.org](https://cn.string-db.org)) and Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases. ## 3. Results ### 3.1. Causal relationship between gut microbiota and fibromyalgia After excluding unnamed bacteria, we identified 196 gut microbiotas. In the preliminary analysis cohort, six gut microbiota showed statistically significant results (**Table S1**), including family *Enterobacteriaceae* (OR = 0.68, 95% CI: 0.47-0.98, p = 0.041), genus *Butyricicoccus* (OR = 0.68, 95% CI: 0.49-0.95, p = 0.023), genus *Coprococcus1* (OR = 0.69, 95% CI: 0.52-0.93, p = 0.014), genus *Eggerthella* (OR = 1.31, 95% CI: 1.08-1.59, p = 0.006), genus *Ruminococcaceae* UCG005 (OR = 1.37, 95% CI: 1.04-1.81, p = 0.025), and order Enterobacteriales (OR = 0.68, 95% CI: 0.47-0.98, p = 0.041). Significant gut microbiota also passed sensitivity tests including Cochran’s Q test (p > 0.05), MR-Egger intercept test (p > 0.05), MR-PRESSO analysis (p > 0.05), and leave-one-out analysis (**Table S2, Figure S1**). We replicated the MR analysis using another GWAS dataset for FM. As expected, similar trends were observed for the candidate gut microbiota in another GWAS dataset (study by Moscati et al.), but results were not significant due to large sample size differences. To strengthen the credibility of the estimates, we performed a meta-analysis of MR results from the preliminary and replication cohorts, identifying six gut microbiota influencing FM (**Figure 2**). Specifically, higher levels of genus *Enterorhabdus* (OR = 0.78, 95% CI: 0.61-1.00, p = 0.048), genus *Parabacteroides* (OR = 0.61, 95% CI: 0.41-0.91, p = 0.016), genus *Butyricicoccus* (OR = 0.73, 95% CI: 0.55-0.96, p = 0.027), genus *Candidatus Soleaferrea* (OR = 0.83, 95% CI: 0.70-1.00, p = 0.046), and genus *Prevotella 9* (OR = 0.83, 95% CI: 0.69-0.99, p = 0.035) reduced the risk of FM. However, higher levels of genus *Ruminococcus gauvreauii* group (OR = 1.56, 95% CI: 1.19-2.04, p = 0.001) increased susceptibility to FM. In the reverse MR analysis (**Table 1**), we found a reverse causal relationship between FM and the genus *Candidatus Soleaferrea* (p < 0.05), indicating that genus *Candidatus Soleaferrea* did not pass the reverse causality validation.  [Figure 2](http://medrxiv.org/content/early/2024/09/14/2024.09.13.24313599/F2) Figure 2 Meta results of gut microbiota and fibromyalgia mendelian randomization analysis. View this table: [Table 1.](http://medrxiv.org/content/early/2024/09/14/2024.09.13.24313599/T1) Table 1. Inverse Mendelian randomization (MR) results of fibromyalgia and significant gut microbiota after meta-analysis. ### 3.2. Causal relationship between 91 inflammatory proteins and fibromyalgia In the preliminary analysis cohort, six inflammatory proteins showed statistically significant results (**Table S3**), including Natural killer cell receptor 2B4 (NKCR-2B4/CD244; OR = 1.18, 95% CI: 1.02-1.38, p = 0.031), CUB domain-containing protein 1 (OR = 0.78, 95% CI: 0.61-0.99, p = 0.037), S100-A12 (OR = 0.78, 95% CI: 0.61-0.99, p = 0.037), Interleukin-12 subunit beta (IL-12β; OR = 1.11, 95% CI: 1.02-1.21, p = 0.016), Monocyte chemoattractant protein 2 (MCP-2,/CCL8; OR = 0.90, 95% CI: 0.84-0.97, p = 0.008), and Tumor necrosis factor (TNF-ɑ; OR = 0.77, 95% CI: 0.61-0.96, p = 0.020). Significant immune proteins also passed sensitivity tests including Cochran’s Q test (p > 0.05), MR-Egger intercept test (p > 0.05), MR-PRESSO analysis (p > 0.05), and leave-one-out analysis (**Table S4, Figure S2**). MR analysis in the replication cohort showed similar trends for candidate immune proteins, but results were not significant due to large sample size differences. We performed a meta-analysis of MR results from the preliminary and replication cohorts, ultimately identifying seven immune proteins influencing FM (**Figure 3**). Specifically, higher levels of C-X-C motif chemokine 5 (CXCL5; OR = 0.89, 95% CI: 0.80-1.00, p = 0.046), S100-A12 (OR = 0.76, 95% CI: 0.62-0.93, p = 0.007), Leukemia inhibitory factor receptor (LIFR; OR = 0.85, 95% CI: 0.73-0.98, p = 0.025), MCP-2/CCL8 (OR = 0.93, 95% CI: 0.87-0.99, p = 0.025), and TNF-α (OR = 0.79, 95% CI: 0.65-0.95, p = 0.013) reduced the risk of FM.  [Figure 3](http://medrxiv.org/content/early/2024/09/14/2024.09.13.24313599/F3) Figure 3 Meta results of inflammatory proteins and fibromyalgia mendelian randomization analysis. However, higher levels of NKCR-2B4/CD244 (OR = 1.20, 95% CI: 1.04-1.38, p = 0.012) and IL-12β (OR = 1.10, 95% CI: 1.02-1.18, p = 0.012) increased susceptibility to FM. All significant results in the meta-analysis passed the reverse causality validation in subsequent analyses (**Table 2**). View this table: [Table 2.](http://medrxiv.org/content/early/2024/09/14/2024.09.13.24313599/T2) Table 2. Inverse Mendelian randomization (MR) results of fibromyalgia and significant inflammatory protein after meta-analysis. For the genes CXCL5, S100-A12, LIFR, CCL8, and TNF, the functional enrichment analyses using GO biological processes and KEGG pathways are shown in **Table S5**. According to results after False Discovery Rate (FDR), the GO analysis shows significant enrichment mainly in the “neutrophil chemotaxis.” KEGG pathway analysis indicates notable enrichment in “viral protein interaction with cytokine and cytokine receptor”. For the genes CD244 and IL-12β, the GO biological process and KEGG functional enrichment analyses are presented in **Table S6**. Post-FDR analysis reveals significant enrichment primarily in “positive regulation of interferon-gamma production” for GO biological processes, while KEGG analysis highlights significant enrichment in the “JAK-STAT signaling pathway.” ## 4. Discussion ### 4.1. Gut microbiota and fibromyalgia Increasing evidence suggests that symbiotic gut bacteria play a significant role in FM (22–25). In our study, we identified *Enterorhabdus, Parabacteroides, Butyricicoccus,* and *Prevotella 9* as protective factors against FM, while the *R. gauvreauii* group was identified as a risk factor. Research on *Enterorhabdus* is currently limited. A study investigating the antidepressant effects of Ziyan Green Tea on mildly stressed mice through fecal metabolomics revealed a significant negative correlation between *Enterorhabdus* and the pro-inflammatory cytokine IL-6 (41). This suggests that *Enterorhabdus* may help inhibit inflammatory responses, thereby acting as a protective factor. *Parabacteroides* has been shown to produce short-chain fatty acids (SCFAs), such as butyrate, which play critical roles in immune regulation (25, 42, 43). *Parabacteroides* may also help repair the intestinal barrier, reduce gut inflammation, and improve metabolic abnormalities like insulin resistance (44), highlighting its potential as a protective factor in FM. *Butyricicoccus* is one of the main producers of butyrate, a key immunomodulatory SCFA that inhibits dendritic cell differentiation and promotes the development of anti-inflammatory regulatory T cells, underscoring the potential role of *Butyricicoccus* in FM (23, 25, 45). Research on *Prevotella 9* in FM is scarce, but studies have suggested that *Prevotella* may act as a potential probiotic (46, 47). For example, in multiple sclerosis, treatment with *Prevotella histicola* increased the frequency and number of CD4FoxP3 regulatory T cells in the periphery and gut, while decreasing pro-inflammatory IFN-γ and IL17-producing CD4 T cells in the CNS, suggesting that *Prevotella histicola* suppresses disease progression by enhancing anti-inflammatory immune responses and inhibiting pro-inflammatory ones (47). Similarly, *Prevotella* may have protective properties through its anti-inflammatory activities. *Ruminococcus*, on the other hand, may activate inflammatory signaling pathways through lipid dysregulation (especially PEA and glycerophosphocholine) and NLRP3 inflammasome activation (45). Importantly, *R. gauvreauii* may induce a systemic immune response by activating IL-6 and TNF-α signaling (48, 49). Moreover, *R. gauvreauii* may be associated with an increased risk for gastro-duodenal ulcers (50). Studies on *Ruminococcus* in depression have shown conflicting results. Zhao et al. found higher levels of *Ruminococcus* in major depressive disorder (MDD) patients, positively correlating with IL-1β, suggesting a pro-inflammatory role that may be important to the pathophysiology of FM (51). Conversely, other studies have shown higher relative abundance of *Ruminococcus* in healthy controls compared to MDD patients (52). In 2007 (Maes et al., 2007), it was reported that the FM symptoms of chronic fatigue syndrome are associated with increased translocation of Gram-negative gut commensal bacteria and gut barrier dysfunction (22). Dysbiosis, particularly the imbalance between pathobionts (which promote inflammation and damage to epithelial cells and tight junctions) and symbionts (health-promoting microbiota, which engage in anti-inflammatory activities, support gut homeostasis and tight junctions, and produce SCFAs and vitamins), is a potential cause of gut barrier dysfunction and bacterial translocation, playing a crucial role in FM (53–56). Consequently, our results may suggest that FM is associated with a specific enterotype, namely with *R. gauvreauii* as a risk factor and *Enterorhabdus, Parabacteroides, Butyricicoccus,* and *Prevotella 9* as protective factors. These results deserve replication in new cohorts of patients with FM. ### 4.2. Cytokines/chemokines in fibromyalgia In this study, five inflammatory proteins—CXCL5, S100-A12, LIFR, MCP-2/CCL8, and TNF-α—were identified as protective factors against FM. Conversely, NKCR-2B4/CD244 and IL-12β were identified as risk factors. Functional enrichment analysis of CXCL5, S100-A12, LIFR, CCL8, and TNF genes revealed significant enrichment in the GO biological process “Cellular response to chemical stimulus,” with other genes also enriched in “Neutrophil chemotaxis,” “Leukocyte migration,” and “Cytokine-mediated signaling pathway”. KEGG pathway analysis indicated that CXCL5, CCL8, LIFR, and TNF genes were enriched in pathways such as “Viral protein interaction with cytokine and cytokine receptor,” “Cytokine-cytokine receptor interaction,” and “Chemokine signaling pathway”. For CD244 and IL-12β, GO biological process analysis showed enrichment in “Positive regulation of interferon-gamma production”; CD244 was also enriched in “Positive regulation of natural killer cell activation” and “Positive regulation of natural killer cell mediated immunity”. KEGG analysis revealed that CD244 was enriched in pathways like “JAK-STAT signaling pathway,” “Inflammatory bowel disease,” and “Th1 and Th2 cell differentiation”. CXCL5, a CXCR2 ligand chemokine, is involved in recruiting leukocytes to tissues, contributing to neuroinflammation, blood-brain barrier breakdown, and neutrophil infiltration into the brain (57). Previous studies have shown that regulating the IL-6 pathway can alleviate pain-like behaviors in FM by reducing pro-inflammatory cytokines such as TNF-α and CXCL5. These observations may seem contradictory to our findings but reflect the complex interplay between immune activation and imbalance in FM, where both immune-inflammatory response system (IRS) and compensatory immune-regulatory response system (CIRS) are activated, leading to imbalances in regulatory and pro-inflammatory cytokines (4). S100A12, a member of the S100 calcium-binding protein family, is mainly secreted by neutrophils, macrophages, and smooth muscle cells. Under inflammatory conditions, S100A12 binds to receptor for advanced glycation end-products (RAGE), signaling through JNK, ERK, and p38 to activate transcription factors AP-1 and nuclear factor (NF)-κB, leading to the release of pro-inflammatory cytokines (58). Elevated S100A12 levels have been observed in the saliva of FM patients (59). Previous studies also identified S100A12 as a diagnostic biomarker for rheumatoid arthritis and depression (60). In a 2002 study by Maes et al, it was found that patients with depressive and anxiety symptoms during the early postpartum period exhibited reduced serum LIFR levels, which may contribute to immune activation (61). Another study indicated that long-term treatment with atypical antipsychotics could increase the anti-inflammatory capacity of serum in schizophrenia patients by elevating serum LIFR levels (62). MCP-2/CCL8, a member of the CC chemokine subfamily, has been reported as an agonist for C-C chemokine receptor types 2 and 4. It plays a crucial role in regulating leukocyte chemotaxis and inflammatory diseases (57). A study by Bäckryd et al. found high levels of the chemokine CCL8 in the cerebrospinal fluid of patients with peripheral neuropathic pain compared to healthy controls (63). In a mouse study, MCP-2/CCL8 levels were increased in oxidative stress-induced allergic inflammation, and anti-inflammatory treatment was shown to reduce MCP-2/CCL8 levels (64). Another study revealed that MCP-2/CCL8 levels were higher in MDD patients with suicidal ideation (SI) compared to those without SI and healthy controls (65). MCP-2/CCL8 may contribute to the neurobiology of suicidality in depressed elderly individuals by enhancing neuroinflammatory responses, oxidative stress, apoptotic pathways, impairing neurogenesis, and disrupting serotonergic neurotransmission (65). TNF-α is one of the most critical cytokines released by macrophages and microglia and is active in neuropathic pain pathways. It modulates pain signaling by binding to TNFR1 (TNF receptor 1) and plays a pivotal role in inflammation at pain sites (66). Research on TNF levels in FM patients has produced mixed results, with some studies reporting lower TNF levels (8, 67, 68) and others finding elevated TNF levels (69, 70). CD244, also known as natural killer cell receptor 2B4, is a transmembrane protein belonging to the immunoglobulin superfamily and plays a vital role in immune regulation. KEGG analysis revealed that CD244 is primarily enriched in the JAK-STAT signaling pathway. Previous studies have shown that this pathway plays a significant role in FM pathogenesis, leading to chemokine overexpression and glial cell activation (66). IL-12β is a pro-inflammatory cytokine and a component of IL-12, involved in activating Th1 immune responses. GO biological process analysis revealed that CD244, along with IL-12β, enhances inflammatory responses through the positive regulation of interferon-gamma production, which may be a risk factor for FM. Previous studies have found elevated IL-12 in patients with FM syndrome (66, 71). Recently, it has been shown that IL-12 also induces inflammatory hypernociception (72). A study demonstrated mechanical hypernociception in rats injected with IL-12, mediated by endothelin-1 and endothelin-B receptor interaction. This finding is significant as pain at the IL-12 injection site has been reported in human treatment (73). IL-12 has also been identified as an essential pro-inflammatory factor in other diseases, such as COVID-19, major depression, psoriasis, and systemic lupus erythematosus (74–78). To summarize, the immunological characteristics observed in the present investigation of FM cannot be solely attributed to a simple “inflammation.” The development of FM may be assigned to an intricate imbalance between the IRS and the CIRS. A recent meta-analysis has substantiated that FM is distinguished by alterations and disparities in immunoregulatory and pro-inflammatory cytokines, resulting in peripheral immunological activation (4). Based on the analysis of enrichment and annotation, the present study indicates that FM may be associated with elevated production of IFN-γ and activation of natural killer cells, as well as reduced chemotaxis and diminished defense against viruses. Indeed, CXCL5, CCL8, and TNF participate in the hsa04061 KEGG pathway (KEGG PATHWAY: hsa04061 (genome.jp)), specifically referred to as “viral protein interaction with cytokine and cytokine receptor”. DNA viruses could undermine the host immune system network in this pathway. A recent study by Maes et al. (2024) demonstrated a significant association between the replication of human herpes type 6 (HHV-6) and the manifestation of FM symptoms in individuals with Long COVID. In addition, reactivation of HHV-6 can potentially impact the integrity of the gut tight junctions, resulting in immune-inflammatory reactions (79). Another report indicated that individuals suffering from FM may exhibit HHV-6 infection (80). Consequently, the replication of HHV-6 may play a role in the immune-inflammatory mechanisms associated with FM by activation the immune-inflammatory response system (81), impairing chemotaxis and inducing antiviral responses via elevated production of IFN-γ (82). Consequently, the elevation of IFN-γ production has the potential to elicit peripheral pain responses via interactions between neurons and microglia within the spinal cord (82, 83). The immune response in FM may be exacerbated by the heightened pro-inflammatory effects of *R. gauvreauii*, alongside a diminished presence of gut bacteria that typically provide a protective barrier against inflammation and the phenomenon of leaky gut. This includes notable genera such as *Enterorhabdus, Parabacteroides, Butyricicoccus,* and *Prevotella 9*. ### 4.3. Limitations When interpreting the results of this study, certain limitations must be acknowledged. This study was conducted on a European population, which may introduce ethnic variability in the findings. Further validation is required through cellular, animal, and clinical studies. ## 5. Conclusion This study presents a comprehensive bidirectional MR analysis and meta-analysis investigating the causal relationships between gut microbiota, inflammatory proteins, and FM. Our findings identify *R. gauvreauii* as a potential risk factor for FM, while *Enterorhabdus*, *Parabacteroides*, *Butyricicoccus*, and *Prevotella 9* appear to act as protective factors. Additionally, five inflammatory proteins—CXCL5, S100-A12, LIFR, MCP-2/CCL8, and TNF-α —were associated with protective roles, whereas NKCR-2B4/CD244 and IL-12β were identified as risk factors. The study provides new insights into the complex interactions between gut microbiota, immune responses, and FM pathophysiology. The dysregulation of immune-inflammatory and compensatory immune-regulatory systems appears central to FM development, with potential links to viral pathways such as HHV-6 reactivation. In conclusion, this study underscores the importance of gut microbiota and immune modulation in FM, offering promising avenues for future therapeutic interventions. ### Ethical statement Our analyses were based on publicly available data that had been approved by relevant review boards. Therefore, no additional ethical approval is required. ### Funding There was no specific funding for this specific study. ### Conflict of interest The authors have no conflicts of interest with any commercial or other association in connection with the submitted article. ## Supporting information Supplementary Figure 1-2 [[supplements/313599_file05.docx]](pending:yes) Supplementary Table 1-6 [[supplements/313599_file06.xlsx]](pending:yes) ## Author’s contributions All the contributing authors have participated in the preparation of the manuscript. ## Data Access Statement The datasets used and/or analyzed during the current study are presented in the manuscript. Summary statistics for GWAS are publicly available. ## Acknowledgment The research community is thanked for making summary statistics from genome-wide association studies publicly available. * Received September 13, 2024. * Revision received September 13, 2024. * Accepted September 14, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Wolfe F, Smythe HA, Yunus MB, Bennett RM, Bombardier C, Goldenberg DL, et al. The American College of Rheumatology 1990 Criteria for the Classification of Fibromyalgia. Report of the Multicenter Criteria Committee. Arthritis and rheumatism. 1990;33(2):160–72. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.1780330203&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2306288&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990CP51800002&link_type=ISI) 2. 2.Sarzi-Puttini P, Giorgi V, Marotto D, Atzeni F. Fibromyalgia: an update on clinical characteristics, aetiopathogenesis and treatment. Nature reviews Rheumatology. 2020;16(11):645–60. 3. 3.Favretti M, Iannuccelli C, Di Franco M. Pain Biomarkers in Fibromyalgia Syndrome: Current Understanding and Future Directions. International journal of molecular sciences. 2023;24(13). 4. 4.Andrés-Rodríguez L, Borràs X, Feliu-Soler A, Pérez-Aranda A, Angarita-Osorio N, Moreno-Peral P, et al. Peripheral immune aberrations in fibromyalgia: A systematic review, meta-analysis and meta-regression. Brain, behavior, and immunity. 2020;87:881–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2019.12.02&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31887417&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 5. 5.Mendieta D, De la Cruz-Aguilera DL, Barrera-Villalpando MI, Becerril-Villanueva E, Arreola R, Hernández-Ferreira E, et al. IL-8 and IL-6 primarily mediate the inflammatory response in fibromyalgia patients. Journal of neuroimmunology. 2016;290:22–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneuroim.2015.11.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26711564&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 6. 6.Bote ME, García JJ, Hinchado MD, Ortega E. An exploratory study of the effect of regular aquatic exercise on the function of neutrophils from women with fibromyalgia: role of IL-8 and noradrenaline. Brain, behavior, and immunity. 2014;39:107–12. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2013.11.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24257330&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 7. 7.Siracusa R, Paola RD, Cuzzocrea S, Impellizzeri D. Fibromyalgia: Pathogenesis, Mechanisms, Diagnosis and Treatment Options Update. International journal of molecular sciences. 2021;22(8). 8. 8.Kosek E, Altawil R, Kadetoff D, Finn A, Westman M, Le Maître E, et al. Evidence of different mediators of central inflammation in dysfunctional and inflammatory pain--interleukin-8 in fibromyalgia and interleukin-1 β in rheumatoid arthritis. Journal of neuroimmunology. 2015;280:49–55. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneuroim.2015.02.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25773155&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 9. 9.Bote ME, García JJ, Hinchado MD, Ortega E. Inflammatory/stress feedback dysregulation in women with fibromyalgia. Neuroimmunomodulation. 2012;19(6):343–51. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000341664&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22986514&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 10. 10.Imamura M, Targino RA, Hsing WT, Imamura S, Azevedo RS, Boas LS, et al. Concentration of cytokines in patients with osteoarthritis of the knee and fibromyalgia. Clinical interventions in aging. 2014;9:939–44. 11. 11.Mastrangelo F, Frydas I, Ronconi G, Kritas SK, Tettamanti L, Caraffa A, et al. Low-grade chronic inflammation mediated by mast cells in fibromyalgia: role of IL-37. Journal of biological regulators and homeostatic agents. 2018;32(2):195–8. 12. 12.Maes M, Almulla AF, Zhou B, Algon AAA, Sodsai P. In major dysmood disorder, physiosomatic, chronic fatigue and fibromyalgia symptoms are driven by immune activation and increased immune-associated neurotoxicity. Scientific reports. 2024;14(1):7344. 13. 13.Maes M, Andrés-Rodríguez L, Vojdani A, Sirivichayakul S, Barbosa DS, Kanchanatawan B. In Schizophrenia, Chronic Fatigue Syndrome- and Fibromyalgia-Like Symptoms are Driven by Breakdown of the Paracellular Pathway with Increased Zonulin and Immune Activation-Associated Neurotoxicity. CNS & neurological disorders drug targets. 2023;22(2):215–25. 14. 14.Asad HN, Al-Hakeim HK, Moustafa SR, Maes M. A Causal-Pathway Phenotype of Chronic Fatigue Syndrome due to Hemodialysis in Patients with End-Stage Renal Disease. CNS & neurological disorders drug targets. 2023;22(2):191–206. 15. 15.Al-Hakeim HK, Khairi Abed A, Rouf Moustafa S, Almulla AF, Maes M. Tryptophan catabolites, inflammation, and insulin resistance as determinants of chronic fatigue syndrome and affective symptoms in long COVID. Frontiers in molecular neuroscience. 2023;16:1194769. 16. 16.Al-Jassas HK, Al-Hakeim HK, Maes M. Intersections between pneumonia, lowered oxygen saturation percentage and immune activation mediate depression, anxiety, and chronic fatigue syndrome-like symptoms due to COVID-19: A nomothetic network approach. Journal of affective disorders. 2022;297:233–45. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jad.2021.10.039&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34699853&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 17. 17.Maes M, Libbrecht I, Van Hunsel F, Lin AH, De Clerck L, Stevens W, et al. The immune-inflammatory pathophysiology of fibromyalgia: increased serum soluble gp130, the common signal transducer protein of various neurotrophic cytokines. Psychoneuroendocrinology. 1999;24(4):371–83. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0306-4530(98)00087-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10341365&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000079762900001&link_type=ISI) 18. 18.Gregory NS, Brito RG, Fusaro M, Sluka KA. ASIC3 Is Required for Development of Fatigue-Induced Hyperalgesia. Molecular neurobiology. 2016;53(2):1020–30. 19. 19.Slade GD, Conrad MS, Diatchenko L, Rashid NU, Zhong S, Smith S, et al. Cytokine biomarkers and chronic pain: association of genes, transcription, and circulating proteins with temporomandibular disorders and widespread palpation tenderness. Pain. 2011;152(12):2802–12. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pain.2011.09.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22000099&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296916600020&link_type=ISI) 20. 20.Sturgill J, McGee E, Menzies V. Unique cytokine signature in the plasma of patients with fibromyalgia. Journal of immunology research. 2014;2014:938576. 21. 21.Littlejohn G, Guymer E. Neurogenic inflammation in fibromyalgia. Seminars in immunopathology. 2018;40(3):291–300. 22. 22.Maes M, Mihaylova I, Leunis JC. Increased serum IgA and IgM against LPS of enterobacteria in chronic fatigue syndrome (CFS): indication for the involvement of gram-negative enterobacteria in the etiology of CFS and for the presence of an increased gut-intestinal permeability. Journal of affective disorders. 2007;99(1-3):237–40. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jad.2006.08.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17007934&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245154300029&link_type=ISI) 23. 23.Minerbi A, Gonzalez E, Brereton NJB, Anjarkouchian A, Dewar K, Fitzcharles MA, et al. Altered microbiome composition in individuals with fibromyalgia. Pain. 2019;160(11):2589–602. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 24. 24.Minerbi A, Fitzcharles MA. Gut microbiome: pertinence in fibromyalgia. Clinical and experimental rheumatology. 2020;38 Suppl 123(1):99–104. 25. 25.Palma-Ordóñez JF, Moreno-Fernández AM, Ramírez-Tejero JA, Durán-González E, Martínez-Lara A, Cotán D. Implication of intestinal microbiota in the etiopathogenesis of fibromyalgia: A systematic review. International journal of rheumatic diseases. 2024;27(1):e15021. 26. 26.Smith GD, Ebrahim S. ’Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? International journal of epidemiology. 2003;32(1):1–22. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyg070&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12689998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182341300001&link_type=ISI) 27. 27.Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genetic epidemiology. 2013;37(7):658–65. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21758&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24114802&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 28. 28.Pierce BL, Burgess S. Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. American journal of epidemiology. 2013;178(7):1177–84. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwt084&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23863760&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325151700023&link_type=ISI) 29. 29.Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Statistics in medicine. 2008;27(8):1133–63. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.3034&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17886233&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 30. 30.Kurilshikov A, Medina-Gomez C, Bacigalupe R, Radjabzadeh D, Wang J, Demirkan A, et al. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nature genetics. 2021;53(2):156–65. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-00763-1&link_type=DOI) 31. 31.Zhao JH, Stacey D, Eriksson N, Macdonald-Dunlop E, Hedman Å K, Kalnapenkis A, et al. Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets. Nat Immunol. 2023;24(9):1540–51. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41590-023-01588-w&link_type=DOI) 32. 32.Moscati A, Faucon AB, Arnaiz-Yépez C, Lönn SL, Sundquist J, Sundquist K, et al. Life is pain: Fibromyalgia as a nexus of multiple liability distributions. American journal of medical genetics Part B, Neuropsychiatric genetics : the official publication of the International Society of Psychiatric Genetics. 2023;192(7-8):171–82. 33. 33.Li P, Wang H, Guo L, Gou X, Chen G, Lin D, et al. Association between gut microbiota and preeclampsia-eclampsia: a two-sample Mendelian randomization study. BMC medicine. 2022;20(1):443. 34. 34.Lin X, Gao W, Huang C, Wu M, She X. Causal relationship between inflammatory proteins and glioblastoma: a two-sample bi[directional mendelian randomization study. Frontiers in genetics. 2024;15:1391921. 35. 35.Sun W, Wu X, Yang H, Yuan S, Chen J, Fang Y, et al. Identifying causal associations between women’s reproductive traits and risk of schizophrenia: a multivariate validated two-sample Mendelian randomization analysis. BMC psychiatry. 2024;24(1):161. 36. 36.Li YS, Xia YG, Liu YL, Jiang WR, Qiu HN, Wu F, et al. Metabolic-dysfunction associated steatotic liver disease-related diseases, cognition and dementia: A two-sample mendelian randomization study. PloS one. 2024;19(2):e0297883. 37. 37.Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ (Clinical research ed). 2018;362:k601. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE2OiIzNjIvanVsMTJfMi9rNjAxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTQvMjAyNC4wOS4xMy4yNDMxMzU5OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 38. 38.Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. European journal of epidemiology. 2017;32(5):377–89. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.1007/s10654-017-0255-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28527048&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 39. 39.Cohen JF, Chalumeau M, Cohen R, Korevaar DA, Khoshnood B, Bossuyt PM. Cochran’s Q test was useful to assess heterogeneity in likelihood ratios in studies of diagnostic accuracy. Journal of clinical epidemiology. 2015;68(3):299–306. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jclinepi.2014.09.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25441698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 40. 40.Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nature genetics. 2018;50(5):693–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0099-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29686387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 41. 41.Jia W, Tang Q, Zou Y, Yang Y, Wu W, Xu W. Investigating the antidepressant effect of Ziyan green tea on chronic unpredictable mild stress mice through fecal metabolomics. Frontiers in microbiology. 2023;14:1256142. 42. 42.Gudi R, Suber J, Brown R, Johnson BM, Vasu C. Pretreatment with Yeast-Derived Complex Dietary Polysaccharides Suppresses Gut Inflammation, Alters the Microbiota Composition, and Increases Immune Regulatory Short-Chain Fatty Acid Production in C57BL/6 Mice. The Journal of Nutrition. 2020;150(5):1291–302. 43. 43.Li M, Yue H, Wang Y, Guo C, Du Z, Jin C, et al. Intestinal microbes derived butyrate is related to the immunomodulatory activities of Dendrobium officinale polysaccharide. International Journal of Biological Macromolecules. 2020;149:717–23. 44. 44.Liu D, Zhang S, Li S, Zhang Q, Cai Y, Li P, et al. Indoleacrylic acid produced by Parabacteroides distasonis alleviates type 2 diabetes via activation of AhR to repair intestinal barrier. BMC biology. 2023;21(1):90. 45. 45.Jia Z, Liu X, Liao W. Unraveling the association between gut microbiota and chemotherapy efficacy: a two-sample Mendelian randomization study. Microbiology spectrum. 2024;12(8):e0394823. 46. 46.Shahi SK, Jensen SN, Murra AC, Tang N, Guo H, Gibson-Corley KN, et al. Human Commensal Prevotella histicola Ameliorates Disease as Effectively as Interferon-Beta in the Experimental Autoimmune Encephalomyelitis. Frontiers in immunology. 2020;11:578648. 47. 47.Shahi SK, Freedman SN, Murra AC, Zarei K, Sompallae R, Gibson-Corley KN, et al. Prevotella histicola, A Human Gut Commensal, Is as Potent as COPAXONE® in an Animal Model of Multiple Sclerosis. Frontiers in immunology. 2019;10:462. 48. 48.Xu T, Ge Y, Du H, Li Q, Xu X, Yi H, et al. Berberis kansuensis extract alleviates type 2 diabetes in rats by regulating gut microbiota composition. Journal of ethnopharmacology. 2021;273:113995. 49. 49.Crost EH, Coletto E, Bell A, Juge N. Ruminococcus gnavus: friend or foe for human health. FEMS microbiology reviews. 2023;47(2). 50. 50.Wu K, Luo Q, Liu Y, Li A, Xia D, Sun X. Causal relationship between gut microbiota and gastrointestinal diseases: a mendelian randomization study. Journal of translational medicine. 2024;22(1):92. 51. 51.Zhao H, Jin K, Jiang C, Pan F, Wu J, Luan H, et al. A pilot exploration of multi-omics research of gut microbiome in major depressive disorders. Translational psychiatry. 2022;12(1):8. 52. 52.Maes M, Vasupanrajit A, Jirakran K, Klomkliew P, Chanchaem P, Tunvirachaisakul C, et al. Exploration of the Gut Microbiome in Thai Patients with Major Depressive Disorder Shows a Specific Bacterial Profile with Depletion of the Ruminococcus Genus as a Putative Biomarker. Cells. 2023;12(9). 53. 53.Simeonova D, Ivanovska M, Murdjeva M, Carvalho AF, Maes M. Recognizing the Leaky Gut as a Trans-diagnostic Target for Neuroimmune Disorders Using Clinical Chemistry and Molecular Immunology Assays. Current topics in medicinal chemistry. 2018;18(19):1641–55. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2174/1568026618666181115100610&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30430944&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 54. 54.Rudzki L, Maes M. The Microbiota-Gut-Immune-Glia (MGIG) Axis in Major Depression. Molecular neurobiology. 2020;57(10):4269–95. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12035-020-01961-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32700250&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 55. 55.Slyepchenko A, Maes M, Jacka FN, Köhler CA, Barichello T, McIntyre RS, et al. Gut Microbiota, Bacterial Translocation, and Interactions with Diet: Pathophysiological Links between Major Depressive Disorder and Non-Communicable Medical Comorbidities. Psychotherapy and psychosomatics. 2017;86(1):31–46. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000448957&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27884012&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 56. 56.Maes M, Vasupanrajit A, Jirakran K, Klomkliew P, Chanchaem P, Tunvirachaisakul C, et al. Adverse childhood experiences and reoccurrence of illness impact the gut microbiome, which affects suicidal behaviours and the phenome of major depression: towards enterotypic phenotypes. Acta neuropsychiatrica. 2023;35(6):328–45. 57. 57.Wang LY, Tu YF, Lin YC, Huang CC. CXCL5 signaling is a shared pathway of neuroinflammation and blood-brain barrier injury contributing to white matter injury in the immature brain. Journal of neuroinflammation. 2016;13:6. 58. 58.Vega-Ramírez MT, Becerril-Villanueva E, Maldonado-García JL, Pavón L, Pérez-Sánchez G. S100 proteins: a new frontier in fibromyalgia research. Molecular brain. 2024;17(1):29. 59. 59.Bazzichi L, Ciregia F, Giusti L, Baldini C, Giannaccini G, Giacomelli C, et al. Detection of potential markers of primary fibromyalgia syndrome in human saliva. Proteomics Clinical applications. 2009;3(11):1296–304. 60. 60.Wang Z, Meng Z, Chen C. Screening of potential biomarkers in peripheral blood of patients with depression based on weighted gene co-expression network analysis and machine learning algorithms. Frontiers in psychiatry. 2022;13:1009911. 61. 61.Maes M, Verkerk R, Bonaccorso S, Ombelet W, Bosmans E, Scharpé S. Depressive and anxiety symptoms in the early puerperium are related to increased degradation of tryptophan into kynurenine, a phenomenon which is related to immune activation. Life sciences. 2002;71(16):1837–48. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0024-3205(02)01853-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12175700&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 62. 62.Maes M, Bocchio Chiavetto L, Bignotti S, Battisa Tura GJ, Pioli R, Boin F, et al. Increased serum interleukin-8 and interleukin-10 in schizophrenic patients resistant to treatment with neuroleptics and the stimulatory effects of clozapine on serum leukemia inhibitory factor receptor. Schizophrenia research. 2002;54(3):281–91. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0920-9964(00)00094-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11950553&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 63. 63.Bäckryd E, Lind AL, Thulin M, Larsson A, Gerdle B, Gordh T. High levels of cerebrospinal fluid chemokines point to the presence of neuroinflammation in peripheral neuropathic pain: a cross-sectional study of 2 cohorts of patients compared with healthy controls. Pain. 2017;158(12):2487–95. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 64. 64. Costa Carvalho JL, de Brito AA, de Oliveira AP, de Castro Faria Neto HC, Pereira TM, de Carvalho RA, et al. The chemokines secretion and the oxidative stress are targets of low-level laser therapy in allergic lung inflammation. Journal of biophotonics. 2016;9(11-12):1208–21. 65. 65.Liu MN, Tsai SJ, Yeh HL, Wu CC, Lin CP. MCP-2/CCL8 Level Associated With Suicidal Ideation in Elderly Men With Major Depression. Archives of suicide research : official journal of the International Academy for Suicide Research. 2020;24(3):467–76. 66. 66.Ellergezen P, Alp A, Çavun S, Çelebi M, Macunluoğlu AC. Pregabalin inhibits proinflammatory cytokine release in patients with fibromyalgia syndrome. Archives of rheumatology. 2023;38(2):307–14. 67. 67.Hernandez ME, Becerril E, Perez M, Leff P, Anton B, Estrada S, et al. Proinflammatory cytokine levels in fibromyalgia patients are independent of body mass index. BMC research notes. 2010;3(1):156. 68. 68.Blanco I, Janciauskiene S, Nita I, Fernández-Bustillo E, Cárcaba V, Gallo C, et al. Low plasma levels of monocyte chemoattractant protein-1 (MCP-1), tumor necrosis factor-alpha (TNFalpha), and vascular endothelial growth factor (VEGF) in patients with alpha1-antitrypsin deficiency-related fibromyalgia. Clinical rheumatology. 2010;29(2):189–97. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19924498&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 69. 69.Tsilioni I, Russell IJ, Stewart JM, Gleason RM, Theoharides TC. Neuropeptides CRH, SP, HK-1, and Inflammatory Cytokines IL-6 and TNF Are Increased in Serum of Patients with Fibromyalgia Syndrome, Implicating Mast Cells. The Journal of pharmacology and experimental therapeutics. 2016;356(3):664–72. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoianBldCI7czo1OiJyZXNpZCI7czo5OiIzNTYvMy82NjQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8xNC8yMDI0LjA5LjEzLjI0MzEzNTk5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 70. 70.Bazzichi L, Rossi A, Massimetti G, Giannaccini G, Giuliano T, De Feo F, et al. Cytokine patterns in fibromyalgia and their correlation with clinical manifestations. Clinical and experimental rheumatology. 2007;25(2):225–30. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17543146&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246882200214&link_type=ISI) 71. 71.Feng J, Zhang Z, Wu X, Mao A, Chang F, Deng X, et al. Discovery of potential new gene variants and inflammatory cytokine associations with fibromyalgia syndrome by whole exome sequencing. PloS one. 2013;8(6):e65033. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0065033&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23762283&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 72. 72.Verri WA, Jr., Cunha TM, Parada CA, Wei XQ, Ferreira SH, Liew FY, et al. IL-15 mediates immune inflammatory hypernociception by triggering a sequential release of IFN-gamma, endothelin, and prostaglandin. Proceedings of the National Academy of Sciences of the United States of America. 2006;103(25):9721–5. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAzLzI1Lzk3MjEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8xNC8yMDI0LjA5LjEzLjI0MzEzNTk5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 73. 73.Verri WA, Jr., Cunha TM, Parada CA, Poole S, Liew FY, Ferreira SH, et al. Antigen-induced inflammatory mechanical hypernociception in mice is mediated by IL-18. Brain, behavior, and immunity. 2007;21(5):535–43. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2006.11.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17194566&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247652200004&link_type=ISI) 74. 74.Sagulkoo P, Plaimas K, Suratanee A, Colado Simão AN, Vissoci Reiche EM, Maes M. Immunopathogenesis and Immunogenetic Variants in COVID-19. Current pharmaceutical design. 2022;28(22):1780–97. 75. 75.Brunoni AR, Supasitthumrong T, Teixeira AL, Vieira EL, Gattaz WF, Benseñor IM, et al. Differences in the immune-inflammatory profiles of unipolar and bipolar depression. Journal of affective disorders. 2020;262:8–15. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jad.2019.10.037&link_type=DOI) 76. 76.Cataldi C, Mari NL, Lozovoy MAB, Martins LMM, Reiche EMV, Maes M, et al. Proinflammatory and anti-inflammatory cytokine profiles in psoriasis: use as laboratory biomarkers and disease predictors. Inflammation research : official journal of the European Histamine Research Society [et al]. 2019;68(7):557–67. 77. 77.Maes M, Ruckoanich P, Chang YS, Mahanonda N, Berk M. Multiple aberrations in shared inflammatory and oxidative & nitrosative stress (IO&NS) pathways explain the co-association of depression and cardiovascular disorder (CVD), and the increased risk for CVD and due mortality in depressed patients. Progress in neuro-psychopharmacology & biological psychiatry. 2011;35(3):769–83. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pnpbp.2010.06.008&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20561554&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F14%2F2024.09.13.24313599.atom) 78. 78.Guimarães PM, Scavuzzi BM, Stadtlober NP, Franchi Santos L, Lozovoy MAB, Iriyoda TMV, et al. Cytokines in systemic lupus erythematosus: far beyond Th1/Th2 dualism lupus: cytokine profiles. Immunology and cell biology. 2017;95(9):824–31. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ICB.2017.53&link_type=DOI) 79. 79.Maes M, Almulla AF, Tang X, Stoyanova K, Vojdani A. From human herpes virus-6 reactivation to autoimmune reactivity against tight junctions and neuronal antigens, to inflammation, depression, and chronic fatigue syndrome due to Long COVID. Journal of medical virology. 2024;96(8):e29864. 80. 80.Krumina A, Chapenko S, Kenina V, Mihailova M, Logina I, Rasa S, et al. The role of HHV-6 and HHV-7 infections in the development of fibromyalgia. Journal of neurovirology. 2019;25(2):194–207. 81. 81.Vojdani A, Vojdani E, Saidara E, Maes M. Persistent SARS-CoV-2 Infection, EBV, HHV-6 and Other Factors May Contribute to Inflammation and Autoimmunity in Long COVID. Viruses. 2023;15(2). 82. 82.Liu CC, Gao YJ, Luo H, Berta T, Xu ZZ, Ji RR, et al. Interferon alpha inhibits spinal cord synaptic and nociceptive transmission via neuronal-glial interactions. Scientific reports. 2016;6:34356. 83. 83.Tsuda M, Masuda T, Kitano J, Shimoyama H, Tozaki-Saitoh H, Inoue K. IFN-gamma receptor signaling mediates spinal microglia activation driving neuropathic pain. Proceedings of the National Academy of Sciences of the United States of America. 2009;106(19):8032–7. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA2LzE5LzgwMzIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8xNC8yMDI0LjA5LjEzLjI0MzEzNTk5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)