- 1Department of Neurology and Psychiatry, Longyou People’s Hospital Affiliated with Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Quzhou, Zhejiang, China
- 2Department of Urology, The Quzhou Affiliated Hospital of Wenzhou Medical University, Quzhou People’s Hospital, Quzhou, Zhejiang, China
Background: Evidence from observational studies indicates that inflammatory proteins play a vital role in Guillain-Barre Syndrome (GBS). Nevertheless, it is unclear how circulating inflammatory proteins are causally associated with GBS. Herein, we conducted a two-sample Mendelian randomization (MR) analysis to systematically explore the causal links of genetically determined systemic inflammatory proteins on GBS.
Methods: A total of 8,293 participants of European ancestry were included in a genome-wide association study of 41 inflammatory proteins as instrumental variables. Five MR approaches, encompassing inverse-variance weighted, weighted median, MR-Egger, simple model, and weighted model were employed to explore the causal links between inflammatory proteins and GBS. MR-Egger regression was utilized to explore the pleiotropy. Cochran’s Q statistic was implemented to quantify the heterogeneity. Furthermore, we performed single-cell RNA sequencing analysis and predicted potential drug targets through molecular docking technology.
Results: By applying MR analysis, four inflammatory proteins causally associated with GBS were identified, encompassing IFN-γ (OR:1.96, 95%CI: 1.02-3.78, PIVW=0.045), IL-7 (OR:1.86, 95%CI: 1.07-3.23, PIVW=0.029), SCGF-β (OR:1.56, 95%CI: 1.11-2.19, PIVW=0.011), and Eotaxin (OR:1.99, 95%CI: 1.01-3.90, PIVW=0.046). The sensitivity analysis revealed no evidence of pleiotropy or heterogeneity. Additionally, significant genes were found through single-cell RNA sequencing analysis and several anti-inflammatory or neuroprotective small molecular compounds were identified by utilizing molecular docking technology.
Conclusions: Our MR analysis suggested that IFN-γ, IL-7, SCGF-β, and Eotaxin were causally linked to the occurrence and development of GBS. These findings elucidated potential causal associations and highlighted the significance of these inflammatory proteins in the pathogenesis and prospective therapeutic targets for GBS.
1 Introduction
Guillain-Barre Syndrome (GBS), also known as acute inflammatory demyelinating polyneuropathy, is an autoimmune disorder affecting the peripheral nervous system. It is primarily triggered by an abnormal immune response to an infectious pathogen (1). The clinical manifestations encompass a progressive, symmetrical paralysis that ascends, accompanied by flaccid limbs and varying degrees of sensory impairments (2). The pathogenesis of GBS remains largely elusive, however, it has been established that the immune inflammatory response plays a crucial role in the onset and progression of this condition (3).
Currently, accumulating evidence suggests that systemic and local release of inflammatory cytokines as well as their involvement in immune-mediated demyelination and axonal damage of peripheral nerves, play a crucial role in the pathogenesis of GBS (4). A previous study has demonstrated that inflammatory proteins play a vital role in the occurrence and biological mechanism of GBS and its experimental autoimmune neuritis in animal models (5). Another previous study has demonstrated that, in GBS, inflammatory cytokines can attract effector cells to the peripheral nerve roots, thereby triggering the production of nitric oxide, oxygen free radicals, and other substances that contribute to the demyelination of nerve roots (6). Interestingly, experimental research has manifested that serum circulating levels of IL-6 and TNF-α are higher than normal in individuals with GBS, thereby predicting the disease severity and other clinical characteristics of GBS (3). In addition, another research conducted by A. R. Exley et al. revealed that raised TNF-α concentrations were found in the GBS patients, suggesting that TNF-α may be a vital factor in the pathogenesis of GBS (7). In contrast, another study performed by E. Doncel-Perez et al. suggested that the levels of IL-2, TNF- β, and TNF-α were significantly lower in GBS individuals when compared with control participants, indicating an association between early growth response gene-2 (EGR2) and these inflammatory regulators in the recovery of GBS (8). Nevertheless, whether these inflammatory proteins have causal associations with GBS remains to be elucidated. Thus, given the susceptibility of observational studies to underlying control bias and confounders (9), further investigation into the potential causal effects of genetically predicted circulating inflammatory proteins and GBS is still required.
Mendelian randomization (MR) analysis can provide an effective means to assess causal inferences in observational studies (10). Conducting randomized controlled trials (RCTs) to detect the influence of circulating inflammatory proteins on the risk of GBS is impractical resulting from ethical conditions. Furthermore, the causal links of circulating inflammatory proteins on GBS remain vague in observational research on account of potential confounding factors. Therefore, alternative methods that facilitate casual inference can offer prospective insights into whether particular inflammatory proteins represent underlying risky factors. Conceptually, MR shares similarities with RCTs as genetic variants are randomly assigned before birth, thereby minimizing confounders and reducing reverse causality bias (11). For example, a previous MR research demonstrated that genetically determined levels of specific inflammatory cytokines exhibited causal associations with the risk of ischemic stroke (12). Furthermore, another MR study revealed that targeted interventions of relevant inflammatory regulators could reduce the risk of proliferative diabetic retinopathy (13).
As of now, existing research evidence is restricted to the effect of a single inflammatory protein on GBS outcome, whether from observational studies or RCTs. To date, the causal links between circulating inflammatory proteins and GBS have yet to be established through MR study. Hence, in the present study, by performing a two-sample MR and single-cell RNA sequencing analysis, we would methodically estimate the potential causal associations between genetically determined circulating inflammatory proteins and GBS. By integrating the findings of observational studies, MR analyses, and single-cell RNA sequencing analysis, we are expected to elucidate the underlying causal roles of specific circulating inflammatory proteins in the pathogenesis of GBS. Moreover, understanding the intricate correlations between systemic inflammatory proteins and the occurrence and progression in depth may be attributed to the development of targeted therapeutic interventions for GBS in the future.
2 Materials and methods
2.1 Study design
The flowchart and design for this study are displayed in Figure 1. Our research is a two-sample Mendelian-randomization (TSMR) study aimed at analyzing the causal relationship between 41 circulating inflammatory proteins and GBS utilizing genome-wide association study (GWAS) summary statistics. Detailed characteristics and information can be found in Supplementary Table S1. To acquire reliable results and augment the robustness of MR analysis findings, three core assumptions demand to be met when conducting TSMR analysis. First, the association between genetic variants and systemic inflammatory proteins is highly significant. Second, the SNPs are independent of potential confounding factors. Third, the instrumental variables can only influence outcomes through exposure factors, and the presence of horizontal pleiotropy is disregarded (14). Before MR analysis, we conducted the relevance assumption test and independent assumption test. We selected genetic instrumental variables for each systemic inflammatory protein to detect the causal association between every inflammatory protein and GBS. The application of publicly attainable GWAS summary datasets obviated the requirement for ethical approval and patient informed consent. All of the MR analysis and sensitivity analysis programs were conducted with the R packages of “Two Sample MR”, “fdrtool”, and “MR-PRESSO” (version 4.2.3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g001.jpg)
Figure 1. Study design and flowchart of two-sample Mendelian randomization for systemic inflammatory proteins and Guillain-Barre Syndrome. pQTL, proteins quantitative trait loci; SNP, single nucleotide polymorphism.
2.2 Data sources
2.2.1 Systemic inflammatory proteins
These 41 systemic inflammatory proteins database was obtained from one of the most comprehensive inflammatory proteins studies by Ahola-Olli et al. (15), encompassing the Cardiovascular Risk in Young Finns Study, FINRISK1997, and FINRISK2002. There was no overlap among these cohorts. Finally, 8,293 Finnish participants from three independent population-based cohorts were included. The mean age of the individuals was 52, and the mean BMI was 27.2kg/m2. The total follow-up span was from 1980 to 2011. These distributions of 41 inflammatory proteins were normalized by utilizing converse transformation. The explicit information and definition of 41 circulating inflammatory proteins were demonstrated in Supplementary Table S2.
2.2.2 Guillain-Barre Syndrome
The summary datasets for GBS participants were obtained from the publicly attainable GWAS database, encompassing 213 cases and 215,718 controls of European descent from the Finnish independent cohort. The analysis was adjusted for both specific principal components and relevant covariates specific to each GWAS, while the overall analysis accounted for principal components. GBS individuals were diagnosed by experienced neurological physicians at each cooperating institution. The diagnosed data were investigated through gathering medical records. The detailed description of GBS GWAS datasets is exhibited in Figure 1 and Supplementary Table S1.
2.3 Selection standards for instrumental variables
We applied comprehensive GWAS summary statistics for 41 systemic inflammatory proteins to extract instrumental variables through a series of control steps. The criteria for eligibility included: (1) given the scarcity of SNPs reaching genome-wide significance, we widened the cutoff value to P<5e-06 to extract qualified instrumental variables, consistent with Z. Song. et al.’s study (16). (2) the influence of robust linkage disequilibrium (LD) between SNPs was eliminated by employing an LD cutoff for the extracted SNPs (r²<0.001, distance>10,000kb), assuring SNPs should be independently associated with the exposure (17). (3) intermediate allele frequency palindromic SNPs were removed due to allele frequencies lacking in the GWAS data of systemic inflammatory proteins. (4) to remove bias arising from weak instrumental variables, we calculated the F-statistics for each SNP to assess statistical strength. The F-statistic was used to assess the strength of each SNP as an instrumental variable, with a value exceeding 10 indicating a robust instrument (18). Furthermore, for each exposure, we harmonized instrumental variables to assure compatibility and consistency between diverse data sources and variables. The brief selection criteria of genetic instrumental variables in this study are shown in Figure 1.
2.4 Statistical analysis
To estimate the causal effects of systemic inflammatory proteins linked to the risk of GBS by uniting various SNPs, we conducted a two-sample Mendelian randomization analysis utilizing five analysis models. The primary approach with the highly statistical power is the inverse-variance weighted (IVW) model which considers all genetic variants to be valid instrumental variables (19). This approach was regarded as the dominant approach in assessing the potential causal relationships between genetically predicted systemic inflammatory proteins and GBS. Furthermore, additional MR analysis methods, encompassing MR-Egger, weighted median, MR pleiotropy residual sum and outlier (MR-PRESSO), simple mode, and weighted mode, were incorporated as auxiliary analysis methods to IVW to provide more robust estimates across a broader range of scenarios (20). The causal effect may be measured via the MR-Egger method if genetic alterations exhibit directional pleiotropy (21). We applied the additional analyses of MR approaches with diverse modeling assumptions and strengths (such as weighted median, simple mode, etc.) to augment the robustness and stability of the findings.
2.5 Sensitivity analysis
To validate and quantify the strength of each instrumental variable concerning the exposure of interest, various sensitivity analysis approaches were also required. Thus, diverse approaches encompassing MR-Egger regression, Cochran’s Q statistics, MR-PRESSO, and leave-one-out (LOO) analysis were applied to validate the robustness of the significant findings (PIVW<0.05). Firstly, the MR-Egger intercept test was implemented to detect the presence of horizontal pleiotropy, with statistical significance defined as P-values for the intercept being less than 0.05 (22). Secondly, the heterogeneity can be assessed by utilizing Cochran’s Q statistic. The absence of heterogeneity among SNPs was indicated by a P-value greater than 0.05 in Cochran’s Q statistics (23). Thirdly, the MR-PRESSO global test was applied to investigate potential outliers as plausible pleiotropic biases and eliminate the influence of pleiotropy by removing the specific SNPs that deviated from normality. Finally, we evaluated the overall robustness with a LOO approach to our results (24). Additionally, the MR Steiger directionality test assures that effective genetic variants can explain a greater proportion of the variance in exposure compared to the outcome. If the genetic tool fails to meet this criterion, it is determined that the genetic variant has a bidirectional effect (25). Eventually, to more accurately identify true signals rather than false positive results caused by chance, we employed FDR (False Discovery Rate) correction to reduce the false positive rate (26). The presence of inflammatory proteins with a P<0.05, although not meeting the FDR-adjusted level of significance, indicated a suggestively causal association with GBS. Therefore, underlying eligible candidate inflammatory proteins for participation in the risk of GBS were identified by the following standards: (1) P-value for the primary MR analysis was significant (PIVW<0.05). (2) The findings of the IVW method were consistent with other approaches in direction and amplitude (scatter plots). (3) The LOO analysis did not identify any significant influential points.
2.6 Confounding analysis
Although we evaluated the heterogeneity and horizontal pleiotropy of the Mendelian randomization (MR) results through multiple sensitivity analysis approaches to explore any SNPs that violated Assumptions 2 and 3 of MR (Figure 1), there may also be little residual confounding instrumental variables. Therefore, to fulfill the independent assumption of MR, SNPs associated with confounders, encompassing blood pressure, body mass index (BMI), and smoking, were excluded by PhenoScannerV2 online tool (http://www.phenoscanner.medschl.cam.ac.uk/) screening, which was an important mean identifying confounding factors. If any SNP was found to be associated with these confounding factors (P<5e-06), MR analysis would be re-conducted after excluding this SNP to ensure the robustness of the findings.
2.7 Mapping candidate SNPs to genes and protein-protein interaction network analysis
We used the SNP-nexus website platform (http://www.snp-nexus.org/v4/) to perform functional annotation of SNPs associated with eligible inflammatory proteins by mapping each queried variant to its closest gene (27). Additionally, we visualized and predicted molecular interaction and PPI network utilizing STRING and GeneMANIA database. The degree algorithm of STRING software was employed to rank the significant proteins in PPI networks. GeneMANIA database can provide a protein-protein interactive network on gene and protein pathways, co-localization, co-expression, and functional assays with pinpoint accuracy of prediction algorithm (28).
2.8 Single-cell RNA sequencing analysis
We gained the single-cell RNA-sequencing (scRNA-seq) data of neuroepithelium tissue from the Panglao DB database (https://panglaodb.se/), a user-centric single-cell sequencing dataset for the scientific community that focused on single-cell RNA sequencing tests from mice and humans (29). We utilized the “Sample” module to search the datasets. “Human” and “Tissue” were set as screening conditions. The “Seurat” R package was employed to analyze the scRNA-seq data (30). This database encompassed 1368 single-cell RNA-sequencing dataset samples. The publicly attainable datasets utilized in the present study had acquired the necessary ethical approvals.
2.9 Molecular docking technology and potential therapeutic drugs prediction
With the RCSD Protein Data Bank (http://www.pdb.org), we acquired the crystal structures of targets (31). PyMOL software was utilized to isolate the primary ligand from the target protein and eliminate extraneous water molecules, phosphates, and other inactive ligands associated with the target protein (32). The 3D structural files of ligands were obtained from TCMSP, and with Chem3D software, their energy was minimized. The protein and active components were stored in the pdb format, while the AutoDockTools 1.5.6 software was employed to convert the pdb format of both the active components and target proteins into the pdbqt format. The parameters for the active pocket were configured, and AutoDock was employed for docking. The active component was deemed to exhibit favorable target binding activity when the binding energy reached or fell below -5.0kJ/mol (33). Discovery Studio 2019 was utilized to visualize the docking results, thereby predicting the potential drug targets.
3 Results
3.1 Selection of the instrumental variables
We performed a two-sample MR analysis to estimate the causal links of genetically determined inflammatory proteins on the risk of GBS. Taking into consideration the limited genetic variance, as well as the restricted number of SNPs and low statistical powers, the MR analysis was performed by widening the cutoff to P-value<5e-06. The number of identified instrumental variables for 41 circulating inflammatory proteins ranged from 3 to 22 (Supplementary Table S2). Furthermore, the F-statistics of all instrumental variables surpassed ten, suggesting that these instrumental variables were sufficiently valid for the MR analysis and that the weak instrumental bias was improbable to occur.
3.2 MR estimates and sensitivity analysis findings of genetically predicted 41 inflammatory proteins on GBS
In this study, to detect the causal links of 41 inflammatory proteins on the risk of GBS, we utilized the IVW approach as the dominant method, complemented by MR-Egger, weighted mode, simple mode, and weighted median. To visually represent the strength of instrumental variables and make the compiled data more intuitive, we created heatmaps to present the estimates of the MR findings for five methods (Figure 2). The detailed MR estimates of diverse approaches are presented in Supplementary Table S3. Moreover, in order to interpret causality more rigorously, the significance threshold was adjusted according to the FDR correction. The PIVW-value of FDR correction was listed in Supplementary Table S4. In addition, due to the IVW approach being susceptible to weak instrumental bias, sensitivity analyses were implemented to assure the reliability and stability of the causality. The sensitivity analysis findings for evaluating the reliability of our MR estimates are demonstrated in Supplementary Table S3. Meanwhile, the findings of MR-PRESSO global test showed no outliers and horizontal pleiotropy (MR-PRESSO global test P>0.05, Supplementary Table S4), suggesting the robustness and credibility of the MR estimates.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g002.jpg)
Figure 2. The heatmaps of five Mendelian randomization analysis methods. Different color blocks represent different odds ratio values. OR, odds ratio.
3.3 Identification of specific inflammatory proteins linked to the risk of GBS
After combining complementary and sensitivity analyses, four eligible inflammatory proteins that met the rigorous screening standards were identified as the candidate inflammatory proteins linked to the risk of GBS. Specifically, as shown in Figure 3, IFN-γ (OR:1.96, 95%CI: 1.02-3.78, PIVW=0.045), IL-7 (OR:1.86, 95%CI: 1.07-3.23, PIVW=0.029), SCGF-β (OR:1.56, 95%CI: 1.11-2.19, PIVW=0.011), and Eotaxin (OR:1.99, 95%CI: 1.01-3.90, PIVW=0.046) were involved in the occurrence and development of GBS. These four cytokines suggested potential causal associations with an elevated risk of GBS. However, after the FDR correction test, they lost the causal associations with GBS. In short, our analysis identified four possible risky factors related to GBS. Furthermore, to augment the reliability and robustness of the findings and improve the interpretability of the study results, we performed reverse MR analysis between genetically predicted GBS and 41 inflammatory proteins. The MR estimates of GBS on systemic inflammatory proteins were displayed in Figure 4 and Supplementary Table S5, indicating that genetically determined GBS showed no causal association with systemic inflammatory proteins. Additionally, scatter plots across various tests of links of four genetically determined inflammatory proteins with the risk of GBS were shown in Figures 5A–D. The findings of LOO analysis similarly manifested that no single SNP had a substantial impact on the stability of analysis outcomes (Figures 6A–D). Additionally, Table 1 presented the sensitivity analysis findings, suggesting the reliability and robustness of our MR analysis results (all P-value>0.05). Moreover, as shown in the funnel plots (Supplementary Figure S1), the distribution of SNPs was symmetrical, indicating that estimates were not violated.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g003.jpg)
Figure 3. Forest plot of the MR estimates for the associations between 41 inflammatory proteins and Guillain-Barre Syndrome. The inverse variance weighted method is considered the main approach. CI, confidence interval; OR, odds ratio.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g004.jpg)
Figure 4. Forest plot of the reverse MR estimates. The inverse variance weighted method is considered the main approach. CI, confidence interval; OR, odds ratio.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g005.jpg)
Figure 5. Scatter plots results of MR analysis between four identified inflammatory proteins and Guillain-Barre Syndrome. (A) IFN-γ, (B) IL-7, (C) SCGF-β, (D) Eotaxin. SNP, single nucleotide polymorphism.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g006.jpg)
Figure 6. Leave-one-out sensitivity analysis results of MR analysis between four identified inflammatory proteins and Guillain-Barre Syndrome. (A) IFN-γ, (B) IL-7, (C) SCGF-β, (D) Eotaxin.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-t001.jpg)
Table 1. Horizontal pleiotropy and heterogeneity analysis for the four identified inflammatory proteins.
3.4 Confounding analysis and protein-protein interactive network
For these four candidate inflammatory proteins, we further manually detected these protein quantity trait loci (pQTL) for the second common trait (alcohol consumption, smoking, BMI, and blood pressure). Looking closely at the PhenoScannerV2 online platform (http://www.phenoscanner.medschl.cam.ac.uk/), we discovered that 41 pQTL were not related to any confounding factors. The corresponding gene information of instrumental variables unraveling these four inflammatory proteins linked to the risk of GBS were shown in Supplementary Table S6. The Manhattan plot exhibited the distribution of genetic locus associated with four identified inflammatory proteins and significant genes were labelled in Figure 7A. Briefly, Manhattan plots revealed that the most significant genes, encompassing POLR1C, ACKR2, VIPR1, SLC34A1, SHANK1, and CLEC11A, were primarily connected with the modulation and function of target genes associated with IFN-γ, IL-7, SCGF-β, and Eotaxin. Moreover, we applied the GeneMANIA online tool to establish and visualize a protein-protein interactive (PPI) network of overlapping hub genes (Figure 7B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g007.jpg)
Figure 7. Manhattan plot exhibited the traits of genetic instrumental variables associated with four candidate inflammatory proteins and significant genes were labelled (A). The most significant gene POLR1C expressed levels and distribution in neuroepithelium tissues through single-cell RNA sequencing analysis (C). Protein-protein interactive (PPI) network by using GeneMANIA platform (B, D). Chr, chromosome.
3.5 Single-cell RNA sequencing analysis results
Single-cell sequencing technology provided cell-specific genetic information and revealed the detailed function and role of target genes. We analyzed scRNA-seq data of the neuroepithelium sample SRA426370 from the Panglao DB database along with cell clustering results and cell type information. After subjecting the data to meticulous quality control procedures, we employed the t-SNE technique to visualize the high-dimensional scRNA-seq data. Data downloaded from Panglao DB demonstrated that the most significant gene POLR1C was expressed in neuroepithelium tissue to varying degrees. As demonstrated in Figure 7C, neuroepithelium cells were segmented into diverse cell clusters, and the gene POLR1C was designated as the overlay expression. Therefore, we might infer that POLR1C may played a vital role in the biological process of GBS, which aligned with the results of the Manhattan plot exhibition (Figure 7A).
3.6 Molecular docking results
Firstly, Figure 7D exhibited the protein-protein interactive network of four candidate target proteins. Secondly, based on the Protein Data Bank (PDB) database, we screened several anti-inflammatory or neuroprotective small molecular compounds, encompassing stigmasterol, saffronin, quercetin, kaempferol, and naringenin. These compounds were selected as the key active ingredients to conduct molecular docking with IL-7. The key components and hub protein were validated by molecular docking utilizing the binding energy that reached or fell below -5.0 kcal/mol as the criterion (34). Finally, the binding energy ranged from -6.4 to -5.1 kcal/mol, and the RMSD was less than 2, suggesting that the significant active ingredients of anti-inflammatory or neuroprotective small molecular compounds were well with the target proteins. This suggests that the findings of our study are reliable, as shown in Figure 8 and Supplementary Table S7.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1456663/fimmu-15-1456663-HTML/image_m/fimmu-15-1456663-g008.jpg)
Figure 8. Molecular docking. Note. Binding mode of proteins and ligands. (A) Binding mode of IL-7 with stigmasterol; (B) Binding mode of IL-7 with saffronin; (C) Binding mode of IL-7 with quercetin; (D) Binding mode of IL-7 with kaempferol; (E) Binding mode of IL-7 with naringenin.
4 Discussion
In this study, we elucidated the causal effects between 41 genetically determined circulating inflammatory proteins and the risk of GBS by using genetic variation as instrumental variables in a two-sample Mendelian randomization analysis. Through strict screening standards and extensive sensitivity analysis methods, four candidate inflammatory proteins (IFN-γ, IL-7, SCGF-β, and Eotaxin) were identified as the significant inflammatory cytokines causally associated with the occurrence and progression of GBS. Furthermore, we employed the Panglao DB database to conduct the single-cell RNA sequencing analysis and discovered significant genes, encompassing POLR1C, ACKR2, VIPR1, SLC34A1, SHANK1, and CLEC11A, were enriched in diverse cell clusters of neuroepithelium tissue. Meanwhile, for the candidate key proteins, we utilized molecular docking techniques to identify several anti-inflammatory or neuroprotective small molecular compounds, encompassing stigmasterol, saffronin, quercetin, kaempferol, and naringenin.
In the present study, our findings manifest that genetically predicted circulating level of IFN-γ has a causally positive association with the risk of GBS. By coincidence, a previous study conducted by H. L. Zhang et al. revealed that IFN-γ facilitated the progression of GBS, which aligned with our results (35). Interestingly, in the research, they also showed that the cytokine IFN-γ played a pivotal role in Th-1-mediated autoimmune disorders by directing the immune response towards a Th1 phenotype through promoting T cell differentiation and inhibiting the development of Th2 cells, as observed in autoimmune disorders like GBS (35). Additionally, they presented a comprehensive overview of the current understanding regarding the inflammatory and immunoregulatory role of IFN-γ in GBS and experimental autoimmune neuritis (EAN), indicating that targeting IFN-γ could be an underlying therapeutic strategy for GBS (35). However, another experimental study performed by R. Press et al. demonstrated that pretreatment levels of blood mononuclear cells (MNC) spontaneously secreting IL-10 were found to be higher in the acute phase of GBS compared to control patients with aseptic meningitis, other neurological disorders, diabetic neuropathy, and healthy individuals (36). The levels of IFN-γ-secreting blood MNC did not show any increase throughout the course of GBS. They concluded that the involvement of IFN-γ in GBS pathogenesis was found to be negligible, while the up-regulation of IL-10 during the early phase of GBS and its association with axonal damage suggested a potential role in the pathogenetic of GBS (36). To sum up, multiple studies have shown different or even opposite results on whether IFN-γ is involved in the biological process of GBS both in laboratory settings and in living organisms. In particular, our findings suggest that genetically determined IFN-γ (OR:1.96, 95%CI: 1.02-3.78, PIVW=0.045) is the most significant risk factor for GBS. Hence, given the ambiguity and inconsistency of these results, further studies are required to elucidate the accurate biological mechanism and validate whether IFN-γ could be treated as a promising therapeutic target of GBS in the future.
Eotaxin, also called CCL11, is a member of the eotaxin family, which is a chemotactic agent for eosinophils that plays a role in innate immunity (37). A previous study conducted by P. Lavandoski et al. demonstrated that CCL11 has the capability to induce cellular senescence and sustain a pro-inflammatory response in human embryonic fibroblasts through the promotion of oxidative stress-induced DNA damage. Furthermore, they reported that the signaling pathway of CCL11 might play a crucial role in establishing a pro-aging environment within the lungs of asthmatic patients (38). Meanwhile, another similar research performed by D. Wu et al. showed that CCL11 played a pivotal role in facilitating the migration of eosinophils to the lungs of individuals with asthma (39). Nevertheless, to date, the association between CCL11 and neurological disorders has been poorly reported. Previous studies performed by Y. Aso et al. confirmed that serum CCL11 was elevated in patients with diabetic sensorimotor polyneuropathy (DSPN) and was associated with peripheral nerve dysfunction, in particular sensory nerve conduction velocity, suggesting that serum CCL11 may be a potential biomarker for DSPN (40). In addition, another research reported by D. Nazarinia et al. manifested that eotaxins may exert both physiological and pathological functions within the central nerve system, while CCL11 can induce neuronal cytotoxic effects by stimulating reactive oxygen species (ROS) production in microglia cells (41). CCL11 plays a vital role in psychosomatic and neuroinflammation, and the plasma levels of CCL11 are elevated in conditions characterized by neuroinflammation and neurodegenerative diseases (41). Disappointingly, the epidemiological evidence supporting the links between CCL11 and GBS was few resulting from its reliance on a case-control study design and constrained by small sample sizes. In the current study, we found that genetically predicted circulating level of CCL11 was causally associated with an elevated risk of GBS (OR:1.99, 95%CI: 1.01-3.90, PIVW=0.046), suggesting that further exploration may be needed to validate the findings of observational studies and investigate the accurate biological mechanism for CCL11 on GBS in the future.
IL-7, also referred to as lymphopoietin, belongs to the IL-2 family cytokines (42). IL-7 is secreted by stromal cells in the bone marrow and thymus of primary lymphoid tissues, as well as epithelial cells in the liver and gut, and endothelial cells (43). IL-7 binds to the IL-7 receptor which is composed of a high-affinity α-submit known as CD127 and the common γ-chain, playing a crucial role in the biological process of T cells development (44). Recently, there has been increasing evidence that IL-7 is involved in a multitude of autoimmune disorders, such as type I diabetes, multiple sclerosis, rheumatoid arthritis, and systemic lupus erythematosus (45). Previous studies showed that targeting components in IL-7 signaling pathways may have potential significance for treating multiple autoimmune diseases (45). As far as we know, GBS is also classified as an autoimmune-mediated peripheral neuropathy. However, we found very little research on the relationship between IL-7 and GBS by reviewing the literature. Our study findings manifested that the circulating level of IL-7 was positively linked to the risk of GBS (OR:1.86, 95%CI: 1.07-3.23, PIVW=0.029). In the future, in order to elucidate the precise mechanisms between IL-7 and GBS, we are expected to conduct in-depth clinical and experimental studies in this field. Finally, the epidemiological evidence of observational studies for the correlation between SCGF-β and GBS was also few. In the present study, we observed that a high genetically determined serum level of SCGF-β was found to be causally linked to an elevated risk of GBS (OR:1.56, 95%CI: 1.11-2.19, PIVW=0.011). On the one hand, one possible explanation was that, in our study, the sample size of GBS individuals was limited, and potential confounding factors were inevitable. On the other hand, the variations of these findings may be attributed to the different choices of instrumental variables and GWAS summary data. Additionally, our research demonstrated that several significant genes, including POLR1C, ACKR2, and VIPR1, were enriched in the diverse cell clusters of neuroepithelium tissues through single-cell RNA sequencing analysis. Hence, it is expected to further validate these findings with bulk RNA analysis in the future.
Our study possessed several advantages. Firstly, the main strength of our study is that it is the first to utilize the MR analysis approach as well as single-cell RNA sequencing analysis to systematically evaluate the causal associations between 41 inflammatory proteins and GBS. This MR design can mitigate confounding biases often presented in observational studies and provide more robust evidence of causality between exposure and outcome. Secondly, due to unbalanced horizontal pleiotropy being a common bias that may distort the MR causal estimates, we extracted the instrumental variables of inflammatory proteins by using the most extensive and updated GWAS attainable datasets. Thirdly, we identified candidate inflammatory proteins causally linked to increased risk of GBS by multiple standards and rigorous control steps, encompassing MR-Egger intercept test, Cochran’s Q statistics, MR-PRESSO, and LOO sensitivity analysis. The sensitivity analysis revealed no horizontal pleiotropy and heterogeneity, suggesting the robustness and credibility of the MR estimates. Furthermore, we employed molecular docking technology to predict underlying candidate drug targets for the treatment of GBS from genetic insights.
Nevertheless, there still exist several limitations in the present study to be considered when interpreting our findings. Initially, the extraction of instrumental variables was conducted by utilizing a widened cutoff (p<5e-06), and this may cause results bias and false-positive variants. Similarly, a previous study has also utilized the identical cutoff when investigating the causal links between 41 circulating cytokines and the risk of ankylosing spondylitis (46). Secondly, the outcome database was only from the FinnGen data and the participants of GWAS only contained European descent. Due to the distinction in genetic backgrounds and living environments, the findings of our study should be cautiously extrapolated to other populations. Thus, additional research involving various ethnic populations is essential to validate our results. Thirdly, despite performing multiple sensitivity analyses to validate MR assumptions, we cannot entirely eliminate the impact of horizontal pleiotropy and heterogeneity. Fourthly, owing to the effects of different circulating inflammatory proteins on the body that may exhibit complicated interaction during the development and evolvement of disease, performing an advanced algorithm for MR study is expected to untangle the complex interaction network of inflammatory proteins. Additionally, our study is restricted to public databases and lacks real-world clinical research, suggesting that these potential links are required to be validated through conducting further experimental and clinical research.
5 Conclusions
In conclusion, we explored the associations between 41 inflammatory proteins and GBS utilizing MR analysis. Our study established suggestively causal links between several specific inflammatory proteins (IFN-γ, IL-7, SCGF-β, and Eotaxin) and the risk of GBS. In particular, we utilized molecular docking technology to predict underlying candidate drug targets. In the future, further validation is warranted to assess the potential of these inflammatory proteins as therapeutic targets or biomarkers for GBS.
Data availability statement
The RNA-seq data presented in this study are available on the link (https://www.ncbi.nlm.nih.gov/Traces/study/?query_key=1&WebEnv=MCID_66d5b9081ba3e568a3bf9dc5&o=acc_s%3Aa). The GWAS dataset for Guillain-Barre Syndrome is obtained from the Open GWAS Project (http://gwas.mrcieu.ac.uk/). Further inquiries can be directed to the corresponding author.
Author contributions
JL: Investigation, Methodology, Software, Validation, Visualization, Writing – original draft. RP: Conceptualization, Data curation, Funding acquisition, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. RP received grants from Science and Technology of Key Project of QuZhou (ZD2022023).
Acknowledgments
We appreciated all the genetics consortiums for making the GWAS summary data publicly available.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1456663/full#supplementary-material
References
1. van den Berg B, Walgaard C, Drenthen J, Fokke C, Jacobs BC, van Doorn PA. Guillain-Barre syndrome: pathogenesis, diagnosis, treatment and prognosis. Nat Rev Neurol. (2014) 10:469–82. doi: 10.1038/nrneurol.2014.121
2. Baars AE, Kuitwaard K, de Koning LC, Luijten LWG, Kok WM, Eftimov F, et al. SARS-CoV-2 vaccination safety in Guillain-Barre syndrome, chronic inflammatory demyelinating polyneuropathy, and multifocal motor neuropathy. Neurology. (2023) 100:e182–91. doi: 10.1212/WNL.0000000000201376
3. Huang P, Xu M, He XY. Correlations between microRNA-146a and immunoglobulin and inflammatory factors in Guillain-Barre syndrome. J Int Med Res. (2020) 48:300060520904842. doi: 10.1177/0300060520904842
4. Nyati KK, Prasad KN. Role of cytokines and Toll-like receptors in the immunopathogenesis of Guillain-Barre syndrome. Mediators Inflammation. (2014) 2014:758639. doi: 10.1155/2014/758639
5. Lu MO, Zhu J. The role of cytokines in Guillain-Barre syndrome. J Neurol. (2011) 258:533–48. doi: 10.1007/s00415-010-5836-5
6. Ethemoglu O, Calik M. Effect of serum inflammatory markers on the prognosis of adult and pediatric patients with Guillain-Barre syndrome. Neuropsychiatr Dis Treat. (2018) 14:1255–60. doi: 10.2147/NDT.S162896
7. Exley AR, Smith N, Winer JB. Tumour necrosis factor-alpha and other cytokines in Guillain-Barre syndrome. J Neurol Neurosurg Psychiatry. (1994) 57:1118–20. doi: 10.1136/jnnp.57.9.1118
8. Doncel-Perez E, Mateos-Hernandez L, Pareja E, Garcia-Forcada A, Villar M, Tobes R, et al. Expression of early growth response gene-2 and regulated cytokines correlates with recovery from Guillain-Barre syndrome. J Immunol. (2016) 196:1102–7. doi: 10.4049/jimmunol.1502100
9. Xu S, Li X, Zhang S, Qi C, Zhang Z, Ma R, et al. Oxidative stress gene expression, DNA methylation, and gut microbiota interaction trigger Crohn's disease: a multi-omics Mendelian randomization study. BMC Med. (2023) 21:179. doi: 10.1186/s12916-023-02878-8
10. Wang X, Cheng H, Feng M, Jiang B, Ren C, Chen Q, et al. Causality of genetically proxied immunophenotypes on cardiovascular diseases: a Mendelian randomization study. Front Immunol. (2024) 15:1344773. doi: 10.3389/fimmu.2024.1344773
11. Ji D, Chen WZ, Zhang L, Zhang ZH, Chen LJ. Gut microbiota, circulating cytokines and dementia: a Mendelian randomization study. J Neuroinflamm. (2024) 21:2. doi: 10.1186/s12974-023-02999-0
12. Li Y, Lu J, Wang J, Deng P, Meng C, Tang H. Inflammatory cytokines and risk of ischemic stroke: A Mendelian randomization study. Front Pharmacol. (2021) 12:779899. doi: 10.3389/fphar.2021.779899
13. Shi Q, Wang Q, Wang Z, Lu J, Wang R. Systemic inflammatory regulators and proliferative diabetic retinopathy: A bidirectional Mendelian randomization study. Front Immunol. (2023) 14:1088778. doi: 10.3389/fimmu.2023.1088778
14. Pan R, Xiao M, Wu Z, Liu J, Wan L. Associations of genetically predicted circulating levels of cytokines with telomere length: a Mendelian randomization study. Front Immunol. (2023) 14:1276257. doi: 10.3389/fimmu.2023.1276257
15. Ahola-Olli AV, Wurtz P, Havulinna AS, Aalto K, Pitkanen N, Lehtimaki T, et al. Genome-wide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors. Am J Hum Genet. (2017) 100:40–50. doi: 10.1016/j.ajhg.2016.11.007
16. Song Z, Li X, Xie J, Han F, Wang N, Hou Y, et al. Associations of inflammatory cytokines with inflammatory bowel disease: a Mendelian randomization study. Front Immunol. (2023) 14:1327879. doi: 10.3389/fimmu.2023.1327879
17. Song H, Lei N, Zeng L, Li X, Jiang C, Feng Q, et al. Mendelian randomization analysis identified tumor necrosis factor as being associated with severe COVID-19. Front Pharmacol. (2023) 14:1171404. doi: 10.3389/fphar.2023.1171404
18. Zheng H, Shi YZ, Liang JT, Lu LL, Chen M. Modifiable factors for migraine prophylaxis: A mendelian randomization analysis. Front Pharmacol. (2023) 14:1010996. doi: 10.3389/fphar.2023.1010996
19. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758
20. Mi J, Liu Z, Jiang L, Li M, Wu X, Zhao N, et al. Mendelian randomization in blood metabolites identifies triglycerides and fatty acids saturation level as associated traits linked to pancreatitis risk. Front Nutr. (2022) 9:1021942. doi: 10.3389/fnut.2022.1021942
21. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
22. Zhong H, Liu S, Zhu J, Xu TH, Yu H, Wu L. Elucidating the role of blood metabolites on pancreatic cancer risk using two-sample Mendelian randomization analysis. Int J Cancer. (2024) 154:852–62. doi: 10.1002/ijc.34771
23. Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. (2017) 28:30–42. doi: 10.1097/EDE.0000000000000559
24. Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. (2017) 32:377–89. doi: 10.1007/s10654-017-0255-x
25. Zhang X, Zhou J, Xie Z, Li X, Hu J, He H, et al. Exploring blood metabolites and thyroid disorders: a bidirectional mendelian randomization study. Front Endocrinol (Lausanne). (2023) 14:1270336. doi: 10.3389/fendo.2023.1270336
26. Yao Z, Guo F, Tan Y, Zhang Y, Geng Y, Yang G, et al. Causal relationship between inflammatory cytokines and autoimmune thyroid disease: a bidirectional two-sample Mendelian randomization analysis. Front Immunol. (2024) 15:1334772. doi: 10.3389/fimmu.2024.1334772
27. Oscanoa J, Sivapalan L, Gadaleta E, Dayem Ullah AZ, Lemoine NR, Chelala C. SNPnexus: a web server for functional annotation of human genome sequence variation (2020 update). Nucleic Acids Res. (2020) 48:W185–92. doi: 10.1093/nar/gkaa420
28. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. (2010) 38:W214–20. doi: 10.1093/nar/gkq537
29. Franzen O, Gan LM, Bjorkegren JLM. PanglaoDB: a web server for exploration of mouse and human single-cell RNA sequencing data. Database (Oxford). (2019) 2019. doi: 10.1093/database/baz046
30. Hao Y, Hao S, Andersen-Nissen E, Mauck 3WM, Zheng S, Butler A, et al. Integrated analysis of multimodal single-cell data. Cell. (2021) 184:3573–87.e29. doi: 10.1016/j.cell.2021.04.048
31. Deng P, Liang H, Wang S, Hao R, Han J, Sun X, et al. Combined metabolomics and network pharmacology to elucidate the mechanisms of Dracorhodin Perchlorate in treating diabetic foot ulcer rats. Front Pharmacol. (2022) 13:1038656. doi: 10.3389/fphar.2022.1038656
32. Warnecke A, Sandalova T, Achour A, Harris RA. PyTMs: a useful PyMOL plugin for modeling common post-translational modifications. BMC Bioinf. (2014) 15:370. doi: 10.1186/s12859-014-0370-6
33. Tao S, Li J, Wang H, Ding S, Han W, He R, et al. Anti-colon cancer effects of Dendrobium officinale Kimura & Migo revealed by network pharmacology integrated with molecular docking and metabolomics studies. Front Med (Lausanne). (2022) 9:879986. doi: 10.3389/fmed.2022.879986
34. Feng C, Zhao M, Jiang L, Hu Z, Fan X. Mechanism of modified Danggui Sini decoction for knee osteoarthritis based on network pharmacology and molecular docking. Evid Based Complement Alternat Med. (2021) 2021:6680637. doi: 10.1155/2021/6680637
35. Zhang HL, Wu L, Wu X, Zhu J. Can IFN-gamma be a therapeutic target in Guillain-Barre syndrome? Expert Opin Ther Targets. (2014) 18:355–63. doi: 10.1517/14728222.2014.882899
36. Press R, Deretzi G, Zou LP, Zhu J, Fredman P, Lycke J, et al. IL-10 and IFN-gamma in Guillain-Barre syndrome. Network members of the Swedish epidemiological study group. J Neuroimmunol. (2001) 112:129–38. doi: 10.1016/S0165-5728(00)00388-X
37. Fan Z, Sun X, Chen X, Liu H, Miao X, Guo Y, et al. C-C motif chemokine CCL11 is a novel regulator and a potential therapeutic target in non-alcoholic fatty liver disease. JHEP Rep. (2023) 5:100805. doi: 10.1016/j.jhepr.2023.100805
38. Lavandoski P, Pierdona V, Maurmann RM, Grun LK, Guma F, Barbe-Tuana FM. Eotaxin-1/CCL11 promotes cellular senescence in human-derived fibroblasts through pro-oxidant and pro-inflammatory pathways. Front Immunol. (2023) 14:1243537. doi: 10.3389/fimmu.2023.1243537
39. Wu D, Zhou J, Bi H, Li L, Gao W, Huang M, et al. CCL11 as a potential diagnostic marker for asthma? J Asthma. (2014) 51:847–54. doi: 10.3109/02770903.2014.917659
40. Aso Y, Iijima T, Jojima T, Saito M, Tanuma D, Kase M, et al. High serum levels of CCL11/Eotaxin-1 are associated with diabetic sensorimotor polyneuropathy and peripheral nerve function but not with cardiac autonomic neuropathy in people with type 2 diabetes. Postgrad Med. (2024) 136:318–24. doi: 10.1080/00325481.2024.2347196
41. Nazarinia D, Behzadifard M, Gholampour J, Karimi R, Gholampour M. Eotaxin-1 (CCL11) in neuroinflammatory disorders and possible role in COVID-19 neurologic complications. Acta Neurol Belg. (2022) 122:865–9. doi: 10.1007/s13760-022-01984-3
42. Ariel A, Hershkoviz R, Cahalon L, Williams DE, Akiyama SK, Yamada KM, et al. Induction of T cell adhesion to extracellular matrix or endothelial cell ligands by soluble or matrix-bound interleukin-7. Eur J Immunol. (1997) 27:2562–70. doi: 10.1002/eji.1830271015
43. Kroncke R, Loppnow H, Flad HD, Gerdes J. Human follicular dendritic cells and vascular cells produce interleukin-7: a potential role for interleukin-7 in the germinal center reaction. Eur J Immunol. (1996) 26:2541–4. doi: 10.1002/eji.1830261040
44. Meyer A, Parmar PJ, Shahrara S. Significance of IL-7 and IL-7R in RA and autoimmunity. Autoimmun Rev. (2022) 21:103120. doi: 10.1016/j.autrev.2022.103120
45. Wang XS, Li BZ, Hu LF, Wen PF, Zhang M, Pan HF, et al. Perspectives of the relationship between IL-7 and autoimmune diseases. Clin Rheumatol. (2013) 32:1703–9. doi: 10.1007/s10067-013-2360-x
Keywords: systemic inflammatory proteins, Guillain-Barre Syndrome, Mendelian randomization, single-cell RNA-seq, molecular docking
Citation: Liu J and Pan R (2024) Causal effects of systemic inflammatory proteins on Guillain-Barre Syndrome: insights from genome-wide Mendelian randomization, single-cell RNA sequencing analysis, and network pharmacology. Front. Immunol. 15:1456663. doi: 10.3389/fimmu.2024.1456663
Received: 02 July 2024; Accepted: 20 August 2024;
Published: 09 September 2024.
Edited by:
Francisco Lopez-Munoz, Camilo José Cela University, SpainReviewed by:
Jiachen Liu, Washington University in St. Louis, United StatesBikash Khadka, Nepal Mediciti Hospital, Nepal
Roshan Aryal, Institute of Medicine, Nepal
Copyright © 2024 Liu and Pan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Renbing Pan, perfectprb@163.com