ORIGINAL RESEARCH article

Front. Aging Neurosci., 01 March 2022

Sec. Parkinson’s Disease and Aging-related Movement Disorders

Volume 14 - 2022 | https://doi.org/10.3389/fnagi.2022.811059

Genetically Predicted Levels of Circulating Inflammatory Cytokines and the Risk and Age at Onset of Parkinson’s Disease: A Two-Sample Mendelian Randomization Study

    YZ

    Yating Zhao 1

    XZ

    Xiaoqian Zhang 1

    NG

    Na Guo 1

    DT

    Dandan Tian 1

    CZ

    Chenguang Zhang 1

    CM

    Changqing Mu 1

    CH

    Chen Han 1

    RZ

    Ruixia Zhu 1

    JZ

    Jian Zhang 2,3

    XL

    Xu Liu 1*

  • 1. Department of Neurology, First Affiliated Hospital of China Medical University, Shenyang, China

  • 2. Key Laboratory of Cell Biology, Ministry of Public Health, Department of Cell Biology, China Medical University, Shenyang, China

  • 3. Key Laboratory of Medical Cell Biology, Ministry of Education, China Medical University, Shenyang, China

Article metrics

View details

8

Citations

4,9k

Views

2,9k

Downloads

Abstract

Parkinson’s disease (PD) is widely considered to be a disabling neurodegenerative disorder, which has been ranked second worldwide just after Alzheimer’s disease. Until present, a wide range of studies has focused on the role of circulating inflammatory cytokines in the development of PD. However, the causal relationship between circulating inflammatory cytokines and the risk and age at the onset of PD has not been elucidated. Hence, to evaluate the effects of circulating inflammatory cytokines on the risk or age at the onset of PD more accurately, we conducted this two-sample Mendelian randomization (MR) study involving summary statistics from genome-wide association studies (GWASs). Totally, we included a GWAS for inflammatory cytokines (8,293 participants), a meta-analysis of GWASs for PD risk (482,730 participants), and a GWAS dataset for age at the onset of PD (17,996 patients with PD). A total of 149 and 131 polymorphisms for exploring relationships between 19 inflammatory cytokines and the risk and age at the onset of PD were obtained as instrumental variants. Then, we used a total of five MR methods, including inverse-variance weighted (IVW), Wald ratio, MR Egger regression, weighted median, and MR-pleiotropy residual sum and outlier (MR-PRESSO) methods. Finally, we found a causal association between circulating levels of macrophage inflammatory protein-1 beta (MIP1b) and PD risk in the IVW method (OR: 1.06; 95% CI: 1.02–1.10; P = 0.001). Meanwhile, other MR estimates by weighted median and MR-PRESSO methods yielded similar effect estimates. Besides, we identified a suggestive association of interleukin-16 (IL-16) levels with PD risk (OR: 1.08; 95% CI: 1.00–1.17; P = 0.037). For age at PD onset, there was no evidence supporting its correlation with inflammatory cytokines. Our findings implied that MIP1b and IL-16 may be novel biomarkers and promising therapeutic targets for PD development.

Introduction

Parkinson’s disease (PD) is considered to be an aging-related neurodegenerative disease, characterized by a broad spectrum of clinical features, including tremor, bradykinesia, rigidity, postural instability, and dysautonomia (Jankovic, 2008). During the past decades, the number of individuals suffering from PD increased from 2.5 million in 1990 to 6.1 million in 2016 globally, and it brought great distress and economic burden to society and families (GBD 2016 Parkinson’s Disease Collaborators, 2018). In the United States alone, the overall annual cost related to PD was estimated to be US$51.9 billion in 2017 and was projected to continually increase and eventually surpass US$79 billion by 2037 (Yang et al., 2020). Thus, huge efforts have been made for exploring possible biomarkers and clarifying the pathogenesis for PD development (Schapira and Jenner, 2011; Lotankar et al., 2017).

Presently, an increasing number of evidence supported the vital role of inflammation in the pathogenesis of PD (Marogianni et al., 2020; Kline et al., 2021). It has been reported that systematic inflammation may destroy the permeability of the blood-brain barrier (BBB), activate microglia, and trigger neuroinflammation, which ultimately led to the degeneration of dopaminergic neurons and PD occurrence (Villaran et al., 2010; Joshi and Singh, 2018). Thus, as an indispensable part of the pathophysiological process of inflammation, inflammatory cytokines are suggested to be involved in PD development. Previous observational studies have demonstrated that circulating levels of several inflammatory cytokines were higher in individuals with PD than in healthy controls. For example, Pochmann et al. (2018) indicated that higher levels of monocyte chemotactic protein-1 (MCP-1) were detected in patients with PD than those in the control group. In a meta-analysis involving 2,654 individuals, elevated concentrations of interleukin-6 (IL-6) and IL-1β were observed in patients with PD but not in healthy controls (Qin et al., 2016). Moreover, several genetic polymorphisms near genes encoding MCP-1 and IL-1β were also indicated to be related to an increased risk of PD (Wahner et al., 2007; Wang et al., 2019). Additionally, it has been reported that higher expression of tumor necrosis factor-α (TNF-α) was associated with earlier onset of PD (Lindenau et al., 2017). However, due to the confounding factors, reverse causation, and limited sample size in these observational studies, the relationships between inflammatory cytokines and the risk and age at the onset of PD may be misled.

Mendelian randomization (MR) is considered to be a robust method, which could overcome the limitations of observational studies mentioned earlier by using genetic variants as instrumental variables and large-scale data from genome-wide association studies (GWASs). Hence, in this study, we introduced the MR method to assess the causal relationship between inflammatory cytokines and the risk and age at the onset of PD.

Materials and Methods

Study Design

In Figure 1, the overall design of this two-sample MR analysis is displayed. The study was based on publicly available data from GWASs on circulating inflammatory cytokines and PD, with detailed information listed in Supplementary Table 1. Since we used summary statistics from published studies, no additional ethical approval was required.

FIGURE 1

FIGURE 1

Schematic representation of two-sample Mendelian randomization analyses for circulating levels of inflammatory cytokines and risk of Parkinson’s disease.

Data Sources and Instrument Selection

Concerning PD risk, we included summary data from a GWAS meta-analysis involving IPDGC-NeuroX, UK Biobank, SGPD and IPDGC study with 33,674 cases and 449,056 controls (contributing studies outlined in Supplementary Table 1; Nalls et al., 2019). As for age at the onset of PD, estimates were extracted from a GWAS dataset including 17 cohorts with 17,996 cases with PD (Blauwendraat et al., 2019).

Regarding the instrument variables of inflammatory cytokines, we obtained genetic variants from one recently published GWAS on circulating inflammatory cytokines involving 8,293 European participants (Ahola-Olli et al., 2017). First, we obtained 6,368 single-nucleotide polymorphisms (SNPs) according to a genome-wide significance threshold of P < 5 × 10–8 and the false discovery rate (FDR) less than 5%. Then, for each cytokine, we pruned all selected SNPs in linkage disequilibrium (LD) (r2 < 0.1 in the European 1,000 genomes reference panel), retaining 199 genetic polymorphisms with the lowest P-value as an independent instrument. Next, to avoid pleiotropic effects, 38 SNPs associated with concentrations of more than one cytokine were excluded, leaving 161 polymorphisms focusing on 19 inflammatory cytokines for use. Finally, since 30 SNPs related to the circulating levels of cytokines were not available in the PD risk or age of PD onset datasets, we selected other SNPs in high LD for replacement in further MR analysis (r2 varying between 0.8 and 1). In the end, a total of 149 SNPs and 131 SNPs for detecting relationships between inflammatory cytokines and the risk and age at the onset of PD were included in this study separately. The detailed information of 19 cytokines and their related SNPs used as instrument variables is displayed in Supplementary Table 2.

Statistical Analyses

First, for those cytokines with only one related SNP, the Wald ratio method was used to compute the MR estimates for the association between one-SD elevated circulating cytokine levels and PD risk or age of PD onset, respectively. When over one variant was included, the inverse-variance weighted (IVW) method, as principal analysis, was used to evaluate the potential causal relationships between circulating levels of inflammatory cytokines and the risk and age at the onset of PD. For further assessing the heterogeneity of the causal estimates across different SNPs, Cochran’s Q-test was applied in the IVW method. Subsequently, we evaluated the possibility that the overall MR estimate was driven by a single SNP using leave-one-out analysis (Noyce et al., 2019). Specifically, we reran the MR analysis by using all SNPs and excluding only one outlier SNP at each time. Furthermore, three supplementary analyses, including MR Egger regression, weighted median, and MR-pleiotropy residual sum and outlier (MR-PRESSO) approaches, were performed for investigating the presence of pleiotropy and further evaluating the causal relationships (Bowden et al., 2015, 2016; Verbanck et al., 2018). In addition, we calculated the F-statistics to quantify the strength of the selected instruments, all of which were above the threshold of the F-statistics (F > 10) typically recommended for MR analyses. As 19 exposures were involved, we set the statistically significant P-value threshold to 2.63 × 10–3 (0.05/19) after the Bonferroni correction. The suggestive association was identified as a P-value between the conventional significance level (0.05) and the Bonferroni-corrected significance level. All above statistical analyses were conducted by using the R-statistical software (version 4.0.3) with related R packages, including MR, two-sample MR, as well as MR-PRESSO.

Results

Circulating Inflammatory Cytokines and Parkinson’s Disease Risk

As displayed in Table 1 and Supplementary Figure 1, by using a significant threshold of 2.63 × 10–3 (Bonferroni correction for the correlation test of 19 cytokines), we found that genetically predicted one-SD increment in circulating MIP1b levels was associated with 6% higher risk of PD based on 74 SNPs in the IVW method (OR: 1.06; 95% CI: 1.02–1.10; P = 0.001). Furthermore, we did not observe any significant heterogeneity as measured by Cochran’s Q-test (I2 = 0.5%, P = 0.465). Subsequent leave-one-out analysis showed that no single SNP dominated the IVW point estimate (Supplementary Figure 2). Similarly, a significant association was found between circulating levels of MIP1b and PD risk in the MR-PRESSO method (OR: 1.06; 95% CI: 1.02–1.10; P = 0.002). Consistent with these results, suggestive evidence of an adverse effect of circulating MIP1b levels on PD risk was also observed by the weighted median method (OR: 1.08; 95% CI: 1.02–1.14; P = 0.014). Although the significant statistical association was not detected by MR Egger regression, the estimate was directionally consistent with other principal and supplementary analyses (Table 1 and Figure 2). Additionally, there was no evidence for potential pleiotropy according to the intercept assessed by MR Egger regression; meanwhile, no outlier SNP was detected using the MR-PRESSO method.

TABLE 1

CytokinesNo. of SNPsOR (95% CI)P for associationHeterogeneity test (I2, P)MR-Egger (intercept, P)P for MR-PRESSO global test
MIP1b
Inverse variance weighted741.06 (1.02–1.10)0.0010.5%, 0.465
MR egger741.06 (0.99–1.14)0.088–0.001, 0.908
Weighted median741.08 (1.02–1.14)0.014
MR-PRESSO (raw, 0 outliers)741.06 (1.02–1.10)0.0020.491
TRAIL
Inverse variance weighted250.98 (0.91–1.05)0.55641.7%, 0.016
MR egger250.93 (0.84–1.03)0.1770.019, 0.200
Weighted median250.99 (0.91–1.07)0.813
MR-PRESSO (raw, 0 outliers)250.98 (0.91–1.05)0.5620.016
IL18
Inverse variance weighted81.07 (0.97–1.19)0.18236.7%, 0.136
MR egger81.39 (1.07–1.79)0.012–0.061, 0.037
Weighted median81.10 (0.98–1.23)0.103
MR-PRESSO (raw, 0 outliers)81.07 (0.97–1.19)0.2240.144
MCP1
Inverse variance weighted70.97 (0.86–1.10)0.6570.0%, 0.776
MR egger71.13 (0.86–1.49)0.376–0.022, 0.233
Weighted median71.01 (0.86–1.18)0.923
MR-PRESSO (raw, 0 outliers)70.97 (0.89–1.06)0.5690.726
GROa
Inverse variance weighted60.98 (0.91–1.06)0.6340.0%, 0.455
MR egger61.01 (0.81–1.27)0.898–0.009, 0.755
Weighted median60.94 (0.85–1.04)0.224
MR-PRESSO (raw, 0 outliers)60.98 (0.91–1.06)0.6440.451
Eotaxin
Inverse variance weighted50.94 (0.82–1.09)0.4410.0%, 0.805
MR egger50.94 (0.49–1.82)0.8550.001, 0.988
Weighted median51.01 (0.85–1.21)0.894
MR-PRESSO (raw, 0 outliers)50.94 (0.86–1.04)0.2930.743
TNFb
Inverse variance weighted41.02 (0.95–1.10)0.5410.0%, 0.434
MR egger41.11 (0.79–1.56)0.532–0.088, 0.612
Weighted median41.03 (0.95–1.12)0.430
MR-PRESSO (raw, 0 outliers)41.02 (0.95–1.10)0.5680.644
CTACK
Inverse variance weighted41.03 (0.92–1.15)0.59321.8%, 0.280
MR egger41.05 (0.77–1.44)0.746–0.007, 0.877
Weighted median41.05 (0.93–1.17)0.444
MR-PRESSO (raw, 0 outliers)41.02 (0.92–1.15)0.6300.363
IL16
Inverse variance weighted31.08 (1.00–1.17)0.0370.0%, 0.705
MR egger31.06 (0.94–1.19)0.3390.010, 0.613
Weighted median31.07 (0.99–1.17)0.095
IL2ra
Inverse variance weighted31.03 (0.95–1.11)0.5380.0%, 0.574
MR egger30.91 (0.69–1.20)0.5060.062, 0.376
Weighted median31.03 (0.95–1.12)0.511
IP10
Inverse variance weighted20.92 (0.73–1.15)0.4480.0%, 0.934
IFNg1
Wald ratio10.88 (0.56–1.40)0.599
IL10
Wald ratio10.87 (0.59–1.28)0.474
IL12p70
Wald ratio11.13 (0.77–1.65)0.531
IL17
Wald ratio10.99 (0.65–1.51)0.979
MCSF
Wald ratio10.83 (0.65–1.06)0.127
MIF
Wald ratio10.98 (0.79–1.23)0.877
MIG
Wald ratio11.04 (0.82–1.32)0.754
RANTES
Wald ratio10.95 (0.62–1.44)0.803

MR analyses of genetically predicted levels of circulating inflammatory cytokines and risk of Parkinson’s disease.

SNP, single-nucleotide polymorphism; OR, odds ratio; CI, confidence interval; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier.

FIGURE 2

FIGURE 2

Scatterplot of genetic association with circulating levels of MIP1b against the genetic association with PD risk. Each black dot indicates an SNP, plotted by the estimate of SNP on the MIP1b levels and the estimate of SNP on PD risk with standard error bars. The slope of the line represents the causal relationship, and each method has a different line. PD, Parkinson’s disease; SNP, single-nucleotide polymorphism.

Regarding interleukin-16 (IL-16), we identified a suggestive association between circulating IL-16 levels and PD risk in IVW analysis. Specifically, for one SD increment of IL-16 levels, the OR of PD risk was 1.08 (95% CI: 1.00–1.17; P = 0.037). Further analysis showed a lack of evidence of heterogeneity among SNPs in the suggestive association of IL-16 with PD risk (I2 = 0, P = 0.705). Meanwhile, additional supplementary analysis displayed a similar causal trend with the estimation using the IVW method. In addition, no potential pleiotropy was detected using the MR Egger method (Table 1).

Apart from MIP1b and IL-16, the other 17 cytokines (TRAIL, MCP1, GROa, eotaxin, TNFb, CTACK, IL2ra, IP10, IFNg1, IL10, IL12p70, IL18, IL17, MCSF, MIF, MIG, and RANTES) were not shown to be associated with PD risk in the main IVW analysis and three supplementary analyses (Table 1). For each cytokine, no marked heterogeneity was found between related SNPs, except for TRAIL (I2 = 41.7%, P = 0.016). Meanwhile, the P-values for the intercepts from Egger regression did not demonstrate any pleiotropy, while IL-18 was the exception (P = 0.037). Further, the MR-PRESSO test showed that no outlier was found for all cytokines, except for TRAIL (P = 0.016).

Circulating Inflammatory Cytokines and Age at the Onset of Parkinson’s Disease

Except for identifying the correlation between inflammatory cytokines and PD risk mentioned earlier, we also detected the causal relationships between the circulating levels of these cytokines and the age of PD onset. Unfortunately, there was no evidence to support the causal relationships between these 19 cytokines and age at the onset of PD (Table 2). Moreover, no statistical evidence of instrument heterogeneity was detected using Cochran’s Q-test except for IP10 (I2 = 79.1%, P = 0.029). When applying MR Egger regression, evidence of potential pleiotropy was only observed in GROa (P = 0.029), while other inflammatory cytokines did not show any pleiotropy. Besides, no SNP was detected as an outlier through the MR-PRESSO test (Table 2).

TABLE 2

CytokinesNo. of SNPsBeta (95% CI)P for associationHeterogeneity test (I2, P)MR-egger (intercept, P)P for MR-PRESSO global test
MIP1b
Inverse variance weighted65–0.02 (–0.26–0.21)0.8540.0%, 0.516
MR egger650.00 (–0.49–0.48)0.986–0.004, 0.934
Weighted median65–0.04 (–0.40–0.32)0.840
MR-PRESSO (raw, 0 outliers)65–0.02 (–0.26–0.21)0.8530.534
TRAIL
Inverse variance weighted21–0.01 (–0.36–0.33)0.9380.0%, 0.939
MR Egger21–0.26 (–0.77–0.26)0.3240.095, 0.208
Weighted median210.04 (–0.42–0.50)0.861
MR-PRESSO (raw, 0 outliers)21–0.01 (–0.27–0.25)0.9190.950
IL18
Inverse variance weighted7–0.19 (–0.80–0.43)0.55320.5%, 0.273
MR Egger7–0.44 (–2.43–1.54)0.6610.061, 0.787
Weighted median7–0.29 (–0.98–0.40)0.411
MR-PRESSO (raw, 0 outliers)7–0.19 (–0.80–0.43)0.5750.328
MCP1
Inverse variance weighted7–0.39 (–1.19–0.41)0.3420.0%, 0.917
MR egger70.33 (–1.56–2.22)0.732–0.106, 0.409
Weighted median7–0.13 (–1.09–0.83)0.796
MR-PRESSO (raw, 0 outliers)7–0.39 (–0.86-0.08)0.1530.922
GROa
Inverse variance weighted6–0.52 (–1.14–0.11)0.10433.9%, 0.182
MR egger60.90 (–0.47–2.28)0.197–0.397, 0.029
Weighted median6–0.43 (–1.07–0.21)0.184
MR-PRESSO (raw, 0 outliers)6–0.52 (–1.14–0.11)0.1650.216
Eotaxin
Inverse variance weighted4–0.10 (–1.12–0.91)0.8440.0%, 0.966
MR egger40.57 (–4.24–5.38)0.816–0.087, 0.779
Weighted median40.00 (–1.20–1.21)0.995
MR-PRESSO (raw, 0 outliers)4–0.10 (–0.40–0.20)0.5560.966
CTACK
Inverse variance weighted40.12 (–0.54–0.78)0.71511.5%, 0.335
MR egger40.67 (–1.08–2.42)0.453–0.165, 0.500
Weighted median40.25 (–0.51–1.01)0.516
MR-PRESSO (raw, 0 outliers)40.12 (–0.54–0.78)0.7390.390
IL16
Inverse variance weighted30.41 (–0.11–0.93)0.1250.0%, 0.839
MR egger30.59 (–0.22–1.40)0.151–0.079, 0.558
Weighted median30.40 (–0.17–0.96)0.166
TNFb
Inverse variance weighted20.03 (–0.53–0.60)0.9050.0%, 0.902
IL2ra
Inverse variance weighted20.45 (–0.33–1.22)0.10750.2%, 0.156
IP10
Inverse variance weighted20.43 (–2.65–3.51)0.54779.1%, 0.029
IFNg1
Wald ratio1–0.79 (–4.24–2.67)0.654
IL10
Wald ratio1–0.06 (–2.57–2.46)0.965
IL12p70
Wald ratio10.44 (–2.12–2.99)0.738
IL17
Wald ratio1–1.18 (–3.80–1.45)0.379
MCSF
Wald ratio1–0.02 (–1.56–1.52)0.983
MIF
Wald ratio1–1.09 (–2.60–0.41)0.153
MIG
Wald ratio1–0.53 (–2.12–1.06)0.510
RANTES
Wald ratio10.87 (–1.56–3.30)0.484

MR analysis of genetically predicted levels of circulating inflammatory cytokines and age at onset of Parkinson’s disease.

SNP, single-nucleotide polymorphism; CI, confidence interval; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier.

Discussion

Parkinson’s disease is a progressive and disabling neurodegenerative disease that mainly affects individuals in their later years of life and its course may vary from 6.9 to 14.3 years (Macleod et al., 2014; Marras et al., 2018). It has been reported that in 2016 alone, PD caused 211,296 deaths and 3.2 million disability-adjusted life-years globally (GBD 2016 Parkinson’s Disease Collaborators, 2018). Despite the long course of the disease and its huge impact on life expectance for elders, the intervention for PD prevention and therapy is still deficient, and the biological mechanism underlying PD etiology is not yet well understood. Thus, we took advantage of a two-sample MR analysis to investigate the causal role of circulating inflammatory cytokines in the risk and age at the onset of PD. In this study, higher MIP1b levels showed a causal relationship with the increased risk of PD, whereas another cytokine, IL-16, displayed a suggestive association.

Macrophage inflammatory protein-1 beta, also known as CCL4, is a chemotactic protein with 69 amino acids produced by various cells, including natural killer cells, T cells, B cells, and neutrophils (Menten et al., 2002; Mukaida et al., 2020). A previous study by Zhu et al. (2014) demonstrated that in the transgenic mice model of Alzheimer’s disease, MIP1b mRNA was 18 times higher than that of wild-type mice, indicating that MIP1b might be involved in the development of Alzheimer’s disease. In addition, MIP1b was also reported to be related to other neurodegenerative diseases, including multiple system atrophy and amyotrophic lateral sclerosis (Compta et al., 2019; Martinez et al., 2020). This evidence implied that MIP1b might participate in the initiation and progression of a neurodegenerative disorder. As for PD, in the past few decades, limited observational studies provided evidence of a possible association between MIP1b and PD risk. According to a previous case-control study involving 50 participants by Calvani et al. (2020), circulating levels of MIP1b were 2-fold higher in patients with PD than in controls, suggesting that it could be considered a possible biomarker for PD occurrence. Additionally, according to Brockmann’s study of 142 patients with PD, higher serum levels of MIP1b may be associated with more severe non-motor symptoms of PD, including cognitive impairment, sleep behavior disorder, and orthostatic dysfunction (Brockmann et al., 2017). Although these studies indicated an association between circulating levels of MIP1b and the occurrence and progression of PD, these results from observational studies were ambiguous and questionable. The reasons were as follows: due to the small sample size and potential confounding, the results of observational studies may be biased. Besides, the relationship between circulating levels of MIP1b and PD risk obtained from observational studies could not be determined as a causal association or even may even be a reverse causality. Therefore, we conducted this MR analysis and reached a reliable result that one-SD increment in the circulating levels of MIP1b may increase the risk of PD by 6%, indicating that lowering the levels of MIP1b may be a promising therapeutic strategy for PD.

From the perspective of biological mechanism, MIP1b may participate in the pathophysiology process of PD through the following pathways. On the one hand, MIP1b could induce the activation of macrophages by binding to its specific receptor, CCR5 (Cheung et al., 2009; Saika et al., 2012). Then, the activated macrophage might release high levels of pro-inflammatory cytokines, such as IL-1β, TNF-α, IL-6, and MIP-1α (Saika et al., 2012; Shapouri-Moghaddam et al., 2018). These cytokines could further disrupt the permeability of brain microvascular endothelial cells and the integrity of BBB (Al-Obaidi and Desa, 2018; Morris et al., 2018). On the other hand, MIP1b not only can immobilize glycosaminoglycans on the apical surface of microvascular endothelial cells but can also bind to the basal surfaces of endothelial cells as well as the subendothelial matrix. Both of the connections contributed to the establishment and maintenance of MIP1b gradients to help various circulating inflammatory cells pass through the BBB and enter the central nervous system (CNS) (Shukaliak and Dorovini-Zis, 2000). Subsequently, in response to these inflammatory mediators, microglia were activated and functionally polarized toward the pro-inflammatory M1 phenotype, which further induced the secretion of inflammatory cytokines, such as TNF-α, MCP-1, and IL-1β (Qin et al., 2007; Jo et al., 2017; Pajares et al., 2020). These chronically elevated inflammatory cytokines and activated inflammatory cells may, in turn, lead to chronic self-perpetuating neuroinflammation and eventually end up in the development of PD (Collins et al., 2012).

Regarding IL-16, it is a protein with 631 amino acids and is generally accepted to be a chemotactic cytokine for T cells (Cruikshank et al., 2000; Skundric et al., 2015). In our study, we found a suggestive association that one-SD increment in the circulating levels of IL-16 was associated with an 8% increase in PD risk. This could be explained by the following reasons: Similar to MIP1b, elevated IL-16 levels might also promote macrophages to release higher levels of pro-inflammatory cytokines and contribute to BBB breakdown (Persidsky et al., 2006; Huang et al., 2019). Besides, IL-16 could induce the activation of protein kinase C (PKC) and its translocation from cytosol to the membrane (Parada et al., 1996). Later, the PKC might phosphorylate the tight junction proteins of BBB leading to BBB damage (Li et al., 2018). Moreover, IL-16 could potentiate the inflammatory response by stimulating T cells to produce more pro-inflammatory cytokines (Skundric et al., 2015). Altogether, the production of pro-inflammatory cytokines, activation of immune cells, and continuous chronic inflammatory response could lead to PD occurrence.

There were several advantages involved in our study. This was the first MR analysis for clarifying the causality of inflammatory cytokines levels with the risk and age at the onset of PD. Another merit was the two-sample MR method that provided robust causal associations by minimizing confounding factors and avoiding reverse causality. Besides, we used data from publicly available GWASs, which included a large number of patients with PD and controls, thus affording a great power to explore the causal relationship. Finally, except for the IVW method, we also used three supplementary analyses, including MR Egger regression, weighted median, and MR-PRESSO approaches. Directional consistency of the results across multiple MR approaches reinforced the reliability of associations. Nevertheless, the results of this MR analysis should be noted in the context of the following limitations. First, since potential pleiotropy, outlier, and heterogeneity can be detected, the relationship between inflammatory cytokines, including TRAIL, IP10, IL-18, and GROa, and the risk or age at the onset of PD should be interpreted with caution. Second, we used summary data from GWASs of mostly European adults, which limited the generality of the observed causal associations to other populations with different genetic backgrounds. Third, cytokines are not something constitutively expressed in the body. During inflammation, specific inflammatory cytokines are induced to a high level but for only a short period. All of those dynamic changes cannot be addressed by MR analysis. Fourth, we focused on the causal effect of peripheral cytokines levels on the risk or age at the onset of PD and supported that systematic peripheral inflammation and chronic cytokine disturbance may contribute to the development of PD. However, a recent study by Wijeyekoon et al. (2020) suggested that some cytokines in the CSF did not correlate well with blood. Thus, our analysis cannot reveal the potential causal relationship between the levels of inflammatory cytokines in CSF and the risk or age at the onset of PD. Last but not the least, not all inflammatory cytokines were analyzed in our MR analysis due to the exclusion criteria and the limited number of cytokines analyzed in the previous GWAS.

Conclusion

In conclusion, this MR analysis showed a causal relationship between circulating levels of MIP1b and PD risk and a suggestive association between IL-16 and PD. Additionally, there was no evidence in support of causality between inflammatory cytokines and age at the onset of PD. Our findings brought new insights into the pathogenesis of PD and indicated that MIP1b and IL-16 could be regarded as novel biomarkers and potential therapeutic targets for PD development.

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.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author contributions

XL designed the research. YZ and XZ had full access to all the data in the study and took responsibility for the integrity of the data and the accuracy of the data analysis. YZ and NG wrote the manuscript and performed the data analysis. All authors contributed to the statistical analysis, critically reviewed the manuscript during the writing process, and approved the final version to be published.

Funding

This work was supported by the National Natural Science Foundation of China (grant nos. 81400950, 81801053, and 81872044) and the Natural Science Foundation of Liaoning Province (grant no. 2019-MS-365).

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.811059/full#supplementary-material

Supplementary Figure 1

Forest plot of the causal effect estimates of MIP1b-associated SNPs on PD risk. Each black dot indicates the causal estimate of MIP1b on PD using an SNP by the Wald ratio method with 95% confidence interval bars. The MR estimate using all SNPs with the IVW method is also shown (red dot with 95% confidence interval bars). PD, Parkinson’s disease; SNP, single-nucleotide polymorphism; MR, Mendelian randomization; IVW, inverse-variance weighted.

Supplementary Figure 2

Leave-one-out analysis of MIP1b-associated SNPs on PD risk. It was conducted to determine whether the causal effect of circulating MIP1b levels on PD risk was disproportionately influenced by a single SNP. The black dot and bars indicate the MR estimate and 95% confidence interval when the specific SNP is removed. The overall MR analysis, including all SNPs, is also presented for comparison (red dot with 95% confidence interval bars). SNP, single-nucleotide polymorphism; PD, Parkinson’s disease; MR, Mendelian randomization.

References

  • 1

    Ahola-OlliA. V.WurtzP.HavulinnaA. S.AaltoK.PitkanenN.LehtimakiT.et al (2017). Genome-wide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors.Am. J. Hum. Genet.1004050. 10.1016/j.ajhg.2016.11.007

  • 2

    Al-ObaidiM. M. J.DesaM. N. M. (2018). Mechanisms of blood brain barrier disruption by different types of bacteria, and bacterial-host interactions facilitate the bacterial pathogen invading the brain.Cell Mol. Neurobiol.3813491368. 10.1007/s10571-018-0609-2

  • 3

    BlauwendraatC.HeilbronK.VallergaC. L.Bandres-CigaS.von CoellnR.PihlstromL.et al (2019). Parkinson’s disease age at onset genome-wide association study: defining heritability, genetic loci, and alpha-synuclein mechanisms.Mov. Disord.34866875. 10.1002/mds.27659

  • 4

    BowdenJ.Davey SmithG.BurgessS. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression.Int. J. Epidemiol.44512525. 10.1093/ije/dyv080

  • 5

    BowdenJ.Davey SmithG.HaycockP. C.BurgessS. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator.Genet. Epidemiol.40304314. 10.1002/gepi.21965

  • 6

    BrockmannK.SchulteC.Schneiderhan-MarraN.ApelA.Pont-SunyerC.VilasD.et al (2017). Inflammatory profile discriminates clinical subtypes in LRRK2-associated Parkinson’s disease.Eur. J. Neurol.24427e426. 10.1111/ene.13223

  • 7

    CalvaniR.PiccaA.LandiG.MariniF.BiancolilloA.Coelho-JuniorH. J.et al (2020). A novel multi-marker discovery approach identifies new serum biomarkers for Parkinson’s disease in older people: an EXosomes in PArkiNson Disease (EXPAND) ancillary study.Geroscience4213231334. 10.1007/s11357-020-00192-2

  • 8

    CheungR.MalikM.RavynV.TomkowiczB.PtasznikA.CollmanR. G. (2009). An arrestin-dependent multi-kinase signaling complex mediates MIP-1beta/CCL4 signaling and chemotaxis of primary human macrophages.J. Leukoc. Biol.86833845. 10.1189/jlb.0908551

  • 9

    CollinsL. M.ToulouseA.ConnorT. J.NolanY. M. (2012). Contributions of central and systemic inflammation to the pathophysiology of Parkinson’s disease.Neuropharmacology6221542168. 10.1016/j.neuropharm.2012.01.028

  • 10

    ComptaY.DiasS. P.GiraldoD. M.Perez-SorianoA.MunozE.SauraJ.et al (2019). Cerebrospinal fluid cytokines in multiple system atrophy: a cross-sectional Catalan MSA registry study.Parkinsonism. Relat. Disord.65312. 10.1016/j.parkreldis.2019.05.040

  • 11

    CruikshankW. W.KornfeldH.CenterD. M. (2000). Interleukin-16.J. Leukoc. Biol.67757766.

  • 12

    GBD 2016 Parkinson’s Disease Collaborators (2018). Global, regional, and national burden of Parkinson’s disease, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016.Lancet Neurol.17939953. 10.1016/S1474-4422(18)30295-3

  • 13

    HuangY.DuK. L.GuoP. Y.ZhaoR. M.WangB.ZhaoX. L.et al (2019). IL-16 regulates macrophage polarization as a target gene of mir-145-3p.Mol. Immunol.10719. 10.1016/j.molimm.2018.12.027

  • 14

    JankovicJ. (2008). Parkinson’s disease: clinical features and diagnosis.J. Neurol. Neurosurg. Psychiatry79368376.

  • 15

    JoM.KimJ. H.SongG. J.SeoM.HwangE. M.SukK. (2017). Astrocytic Orosomucoid-2 modulates microglial activation and neuroinflammation.J. Neurosci.3728782894. 10.1523/JNEUROSCI.2534-16.2017

  • 16

    JoshiN.SinghS. (2018). Updates on immunity and inflammation in Parkinson disease pathology.J. Neurosci. Res.96379390. 10.1002/jnr.24185

  • 17

    KlineE. M.HouserM. C.HerrickM. K.SeiblerP.KleinC.WestA.et al (2021). Genetic and environmental factors in Parkinson’s disease converge on immune function and inflammation.Mov. Disord.362536. 10.1002/mds.28411

  • 18

    LiH.ZhuL.FengJ.HuX.LiC.ZhangB. (2018). Hydrogen sulfide decreases blood-brain barrier damage via regulating protein Kinase C and tight junction after cardiac arrest in rats.Cell Physiol. Biochem.479941006. 10.1159/000490166

  • 19

    LindenauJ. D.AltmannV.Schumacher-SchuhA. F.RiederC. R.HutzM. H. (2017). Tumor necrosis factor alpha polymorphisms are associated with Parkinson’s disease age at onset.Neurosci. Lett.658133136. 10.1016/j.neulet.2017.08.049

  • 20

    LotankarS.PrabhavalkarK. S.BhattL. K. (2017). Biomarkers for Parkinson’s disease: recent advancement.Neurosci. Bull.33585597. 10.1007/s12264-017-0183-5

  • 21

    MacleodA. D.TaylorK. S.CounsellC. E. (2014). Mortality in Parkinson’s disease: a systematic review and meta-analysis.Mov. Disord.2916151622.

  • 22

    MarogianniC.SokratousM.DardiotisE.HadjigeorgiouG. M.BogdanosD.XiromerisiouG. (2020). Neurodegeneration and inflammation-an interesting interplay in Parkinson’s disease.Int. J. Mol. Sci.21:8421. 10.3390/ijms21228421

  • 23

    MarrasC.BeckJ. C.BowerJ. H.RobertsE.RitzB.RossG. W.et al (2018). Prevalence of Parkinson’s disease across North America.NPJ Parkinsons Dis.4:21. 10.1038/s41531-018-0058-0

  • 24

    MartinezH. R.Escamilla-OcanasC. E.Camara-LemarroyC. R.Gonzalez-GarzaM. T.Moreno-CuevasJ.Garcia SarreonM. A. (2020). Increased cerebrospinal fluid levels of cytokines monocyte chemoattractant protein-1 (MCP-1) and macrophage inflammatory protein-1beta (MIP-1beta) in patients with amyotrophic lateral sclerosis.Neurologia35165169. 10.1016/j.nrl.2017.07.020

  • 25

    MentenP.WuytsA.Van DammeJ. (2002). Macrophage inflammatory protein-1.Cytokine Growth Factor Rev.13455481.

  • 26

    MorrisG.FernandesB. S.PuriB. K.WalkerA. J.CarvalhoA. F.BerkM. (2018). Leaky brain in neurological and psychiatric disorders: drivers and consequences.Aust. N. Z. J. Psychiatry52924948. 10.1177/0004867418796955

  • 27

    MukaidaN.SasakiS. I.BabaT. (2020). CCL4 Signaling in the tumor microenvironment.Adv. Exp. Med. Biol.12312332. 10.1007/978-3-030-36667-4_3

  • 28

    NallsM. A.BlauwendraatC.VallergaC. L.HeilbronK.Bandres-CigaS.ChangD.et al (2019). Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: a meta-analysis of genome-wide association studies.Lancet Neurol.1810911102. 10.1016/S1474-4422(19)30320-5

  • 29

    NoyceA. J.Bandres-CigaS.KimJ.HeilbronK.KiaD.HemaniG.et al (2019). The Parkinson’s disease mendelian randomization research portal.Mov. Disord.3418641872. 10.1002/mds.27873

  • 30

    PajaresM.RojoA. I.MandaG.BoscaL.CuadradoA. (2020). Inflammation in Parkinson’s disease: mechanisms and therapeutic implications.Cells9:1687. 10.3390/cells9071687

  • 31

    ParadaN. A.CruikshankW. W.DanisH. L.RyanT. C.CenterD. M. (1996). IL-16- and other CD4 ligand-induced migration is dependent upon protein kinase C.Cell Immunol.168100106. 10.1006/cimm.1996.0054

  • 32

    PersidskyY.RamirezS. H.HaorahJ.KanmogneG. D. (2006). Blood-brain barrier: structural components and function under physiologic and pathologic conditions.J. Neuroimmune Pharmacol.1223236. 10.1007/s11481-006-9025-3

  • 33

    PochmannD.PeccinP. K.da SilvaI. R. V.DornelesG. P.PeresA.NiqueS.et al (2018). Cytokine modulation in response to acute and chronic aquatic therapy intervention in Parkinson disease individuals: a pilot study.Neurosci. Lett.6743035. 10.1016/j.neulet.2018.03.021

  • 34

    QinL.WuX.BlockM. L.LiuY.BreeseG. R.HongJ. S.et al (2007). Systemic LPS causes chronic neuroinflammation and progressive neurodegeneration.Glia55453462. 10.1002/glia.20467

  • 35

    QinX. Y.ZhangS. P.CaoC.LohY. P.ChengY. (2016). Aberrations in peripheral inflammatory cytokine levels in parkinson disease: a systematic review and meta-analysis.JAMA Neurol.7313161324. 10.1001/jamaneurol.2016.2742

  • 36

    SaikaF.KiguchiN.KobayashiY.FukazawaY.KishiokaS. (2012). CC-chemokine ligand 4/macrophage inflammatory protein-1beta participates in the induction of neuropathic pain after peripheral nerve injury.Eur. J. Pain1612711280. 10.1002/j.1532-2149.2012.00146.x

  • 37

    SchapiraA. H.JennerP. (2011). Etiology and pathogenesis of Parkinson’s disease.Mov. Disord.2610491055.

  • 38

    Shapouri-MoghaddamA.MohammadianS.VaziniH.TaghadosiM.EsmaeiliS. A.MardaniF.et al (2018). Macrophage plasticity, polarization, and function in health and disease.J. Cell. Physiol.23364256440. 10.1002/jcp.26429

  • 39

    ShukaliakJ. A.Dorovini-ZisK. (2000). Expression of the beta-chemokines RANTES and MIP-1 beta by human brain microvessel endothelial cells in primary culture.J. Neuropathol. Exp. Neurol.59339352. 10.1093/jnen/59.5.339

  • 40

    SkundricD. S.CruikshankW. W.MontgomeryP. C.LisakR. P.TseH. Y. (2015). Emerging role of IL-16 in cytokine-mediated regulation of multiple sclerosis.Cytokine75234248. 10.1016/j.cyto.2015.01.005

  • 41

    VerbanckM.ChenC. Y.NealeB.DoR. (2018). Publisher Correction: Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases.Nat. Genet.50:1196. 10.1038/s41588-018-0164-2

  • 42

    VillaranR. F.Espinosa-OlivaA. M.SarmientoM.De PablosR. M.ArguellesS.Delgado-CortesM. J.et al (2010). Ulcerative colitis exacerbates lipopolysaccharide-induced damage to the nigral dopaminergic system: potential risk factor in Parkinson‘s disease.J. Neurochem.11416871700. 10.1111/j.1471-4159.2010.06879.x

  • 43

    WahnerA. D.SinsheimerJ. S.BronsteinJ. M.RitzB. (2007). Inflammatory cytokine gene polymorphisms and increased risk of Parkinson disease.Arch. Neurol.64836840. 10.1001/archneur.64.6.836

  • 44

    WangY.ZhouM.WangY.JiangD.DengX. (2019). Association of polymorphisms in the MCP-1 and CCR2 genes with the risk of Parkinson’s disease.J. Neural. Transm.12614651470. 10.1007/s00702-019-02072-2

  • 45

    WijeyekoonR. S.MooreS. F.FarrellK.BreenD. P.BarkerR. A.Williams-GrayC. H. (2020). Cerebrospinal fluid cytokines and neurodegeneration-associated proteins in Parkinson’s disease.Mov. Disord.3510621066. 10.1002/mds.28015

  • 46

    YangW.HamiltonJ. L.KopilC.BeckJ. C.TannerC. M.AlbinR. L.et al (2020). Current and projected future economic burden of Parkinson’s disease in the U.S.NPJ Parkinsons Dis.6:15. 10.1038/s41531-020-0117-1

  • 47

    ZhuM.AllardJ. S.ZhangY.PerezE.SpanglerE. L.BeckerK. G. (2014). Age-related brain expression and regulation of the chemokine CCL4/MIP-1beta in APP/PS1 double-transgenic mice.J. Neuropathol. Exp. Neurol.73362374. 10.1097/NEN.0000000000000060

Summary

Keywords

cytokines, inflammation, Mendelian randomization, macrophage inflammatory protein-1 beta, interleukin-16, Parkinson’s disease

Citation

Zhao Y, Zhang X, Guo N, Tian D, Zhang C, Mu C, Han C, Zhu R, Zhang J and Liu X (2022) Genetically Predicted Levels of Circulating Inflammatory Cytokines and the Risk and Age at Onset of Parkinson’s Disease: A Two-Sample Mendelian Randomization Study. Front. Aging Neurosci. 14:811059. doi: 10.3389/fnagi.2022.811059

Received

08 November 2021

Accepted

31 January 2022

Published

01 March 2022

Volume

14 - 2022

Edited by

Robert Petersen, Central Michigan University, United States

Reviewed by

Pingyi Xu, First Affiliated Hospital of Guangzhou Medical University, China; Zhaoqi Yan, Gladstone Institutes, United States

Updates

Copyright

*Correspondence: Xu Liu,

†These authors have contributed equally to this work and share first authorship

This article was submitted to Parkinson’s Disease and Aging-related Movement Disorders, a section of the journal Frontiers in Aging Neuroscience

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics