Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Causal effects of circulating cytokine concentrations on risk of Alzheimer’s disease: A bidirectional two-sample Mendelian randomization study

View ORCID ProfilePanagiota Pagoni, Laura D Howe, George Davey Smith, Yoav Ben-Shlomo, Evie Stergiakouli, Emma L Anderson
doi: https://doi.org/10.1101/2020.11.18.20232629
Panagiota Pagoni
1MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK
2Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Panagiota Pagoni
  • For correspondence: panagiota.pagoni@bristol.ac.uk
Laura D Howe
1MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK
2Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
George Davey Smith
1MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK
2Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yoav Ben-Shlomo
2Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Evie Stergiakouli
1MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK
2Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Emma L Anderson
1MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK
2Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

ABSTRACT

Background There is considerable interest in the role of neuroinflammation in the pathogenesis of Alzheimer’s disease. Evidence from observational studies suggests an association between cytokine concentrations and Alzheimer’s disease. However, establishing a causal role of cytokine concentrations on risk of Alzheimer’s disease is challenging due to bias from reverse causation and residual confounding.

Methods We used two-sample MR to explore causal effects of circulating cytokine concentrations on Alzheimer’s disease and vice versa, employing genetic variants associated with cytokine concentrations (N=8,293) and Alzheimer’s disease (71,880 cases / 383,378 controls) from the largest non-overlapping genome-wide association studies (GWAS) of European ancestry.

Results There was weak evidence to suggest that 1 standard deviation (SD) increase in levels of CTACK (CCL27) (OR= 1.09 95%CI: 1.01 to 1.19, p=0.03) increased risk of Alzheimer’s disease. There was also weak evidence of a causal effect of 1 SD increase in levels of MIP-1b (CCL4) (OR=1.04 95%CI: 0.99 to 1.09, p=0.08), Eotaxin (OR=1.08 95%CI: 0.99 to 1.17, p =0.10), GROa (CXCL1) (OR=1.04 95%CI: 0.99 to 1.10, p=0.15), MIG (CXCL9) (OR=1.17 95%CI: 0.97 to 1.41, p=0.10), IL-8 (Wald Ratio: OR=1.21 95%CI: 0.97 to 1.51, p=0.09) and IL-2 (Wald Ratio: OR=1.21 95%CI: 0.94 to 1.56, p=0.14) on greater risk of Alzheimer’s disease. There was little evidence of a causal effect of genetic liability to Alzheimer’s disease on circulating cytokine concentrations.

Conclusions Our study provides some evidence supporting a causal role of cytokines in the pathogenesis of Alzheimer’s disease. However, more studies are needed to elucidate the specific mechanistic pathways via which cytokines alter the risk of Alzheimer’s disease.

INTRODUCTION

It has been estimated that 47 million people are affected by Alzheimer’s disease and other forms of dementia worldwide and this number is likely to increase, mostly due to ageing of the population (1). Thus, dementia is a major public health challenge and priority.

Over the last decade attention has been drawn to the interplay between the central nervous system (CNS) and immune responses (i.e. neuroinflammation) (2) in the pathogenesis of Alzheimer’s disease, mostly due to evidence stemming from observational studies suggesting that inflammatory disorders (e.g. rheumatoid arthritis) (3) and chronic inflammation (e.g. periodontitis) (4) are associated with a higher risk of Alzheimer’s disease. The neuroinflammation hypothesis suggests that in response to production and deposition of Aβ in the brain, the CNS activates microglia to protect the cells and overall brain function. As part of this defense mechanism, secondary inflammatory mediators such as cytokines (5, 6), lipid metabolites and free radicals are generated to rehabilitate the homeostasis and ensure a healthy neural function (7). However, over-activation of microglia may occur and lead to an exaggerated release of pro- and anti-inflammatory mediators (8), resulting in deterioration or even initiation of neurological diseases.

Many observational studies have examined the association between cytokine concentrations and risk of Alzheimer’s disease. A recent meta-analysis of 170 studies reported elevated peripheral levels of C-reactive protein (CRP), interleukin-6 (IL-6), interleukin-1 beta (IL-1β), soluble tumour necrosis factor receptor 1 and 2 (sTNFR-1 & sTNFR-2), interleukin-10 (IL-10), monocyte chemoattractant protein-1 (MCP-1) and transforming growth factor-1 (TGF-1) in individuals with an Alzheimer’s disease diagnosis compared to healthy controls (9). However, deciphering the role of inflammatory markers in the pathogenesis of Alzheimer’s disease is challenging because of the potential bias due to reverse causation, which refers to the possibility of Alzheimer’s disease being the cause rather than the consequence of inflammation. Confounding is another key source of bias inherent in existing observational studies, as Alzheimer’s disease patients tend to be older, and thus are more at risk of conditions such as obesity and hypertension, which may increase systemic inflammation.

A promising method that could help overcome the limitations of observational studies is Mendelian randomization (MR). Mendelian randomization uses genetic variants as proxies for an exposure and allows investigation of the causal effect of the exposure on the outcome of interest (10, 11). Moreover, two-sample MR studies, where the effect estimates of the genetic variants of exposure and outcome are extracted from different genome-wide association studies (GWAS) (12) can be useful, as they allow the estimation of causal effects without requiring the exposure and outcome to be measured in the same participants. Deciphering the causal role of neuroinflammation in the pathogenesis of Alzheimer’s disease could provide valuable information towards possible therapeutic targets. However, to date, very few studies have used MR to examine this question (13-16), and existing MR studies have included only a small, select group of previously implicated cytokines such as TNF-a, interleukins and CRP. Considering the plethora of observational evidence indicating a role of cytokines in the pathogenesis of Alzheimer’s disease and the limitations of observational studies to infer causality, we aimed to examine the bidirectional causal effects of cytokines on risk of Alzheimer’s disease in a MR framework.

METHODS

GWAS Summary Data

We used the most comprehensive and largest publicly available GWAS meta-analysis on concentrations of 41 circulating cytokines including up to 8,293 individuals from three independent population cohorts (The Cardiovascular Risk in Young Finns Study (YFS), FINRISK 1997 and FINRISK 2002)(17). Cytokines GWAS were adjusted for sex, age and body mass index (BMI).

For Alzheimer’s disease the most recently conducted GWAS was used, which consists of three phases. In phase one, 24,087 clinically diagnosed cases from the International Genomics of Alzheimer’s Project (IGAP), the Alzheimer’s Disease Sequencing Project (ADSP) and the Alzheimer’s disease working group of the Psychiatric Genomics Consortium (PGZ-ALZ) and 55,058 matched controls were included. In phase two, cases consist of 47,793 proxy-cases as defined in the UK Biobank and 328,320 proxy-controls (18). Participants were considered proxy-cases if they had positively responded to the question ‘Has your mother or father ever suffered from Alzheimer’s disease/Dementia’. Finally, phase three is a meta-analysis of all individuals in phase one & two and therefore consists of 71,880 cases and 383,378 controls of European ancestry. In our main analysis, we used data from phase one. As these summary data were not publicly available, we used summary estimates from Korologou-Linden et al. (19), which corresponds to phase one of the Alzheimer’s GWAS. Effect estimates were adjusted for sex, age, genotyping array and assessment centre. This study, assessed the effect of genetic liability to Alzheimer’s disease on circulating cytokine concentrations, as captured by common genetic variants. When examining binary exposures in MR settings, causal inferences are valid for the continuous liability underlying the binary exposure (20). Thus, when we test for the causal effect of Alzheimer’s disease, we are essentially examining the effect of genetic liability for this exposure which can be present in an individual even when they do not have a diagnosis.

Instrument selection

For each exposure of interest, approximately independent genome-wide significant single nucleotide polymorphisms (SNPs) were identified (r2<0.01 within a 10,000 kb window, p<5×10−08). Eight cytokines had no SNPs available at the genome wide significance threshold, thus we relaxed the significance threshold to p<5×10−07 for these. When investigating the causal effects of cytokine concentrations on Alzheimer’s disease, we extracted the SD-scaled effect sizes and standard errors for cytokine SNPs from the publicly available cytokine GWAS, and corresponding log odds and standard errors from the Alzheimer’s disease GWAS. When investigating the reverse direction, we extracted log odds and standard errors for Alzheimer’s disease SNPs from the Alzheimer’s GWAS and corresponding SD-scaled effect sizes and standard errors from the cytokines GWAS.

Before estimating total causal effects, we harmonised the alleles of our datasets and further information about the procedure followed can be found in Supplementary Note 1.

Mendelian randomization analysis

Univariable MR was employed to estimate the total causal effect of each circulating cytokine concentration on Alzheimer’s disease, and vice versa. MR relies on three assumptions that the genetic variants should satisfy to be considered as valid instruments and therefore yield unbiased causal effect estimates. Genetic variants i) must be strongly associated with the exposure of interest, ii) independently of any confounders of the exposure – outcome association and iii) they are associated with the outcome only via the exposure (i.e. no horizontal pleiotropy) (21).

When a single variant was available as a proxy for the exposure of interest, the Wald ratio estimator was employed to quantify the causal effect. When multiple variants were available, the Inverse-Variance-Weighted (IVW) method was used to estimate the total causal effect, which is equivalent to fitting a weighted linear regression of the gene-outcome associations on the gene-exposure associations, with the intercept term constrained to zero. Therefore, IVW estimates assume that all genetic variants are valid instruments with no pleiotropic effects (22, 23).

Sensitivity analyses

We performed a series of sensitivity analyses to test the validity of the core assumptions which MR relies on (Table 1). A more comprehensive description of these sensitivity analyses can be found in Supplementary Note 2 and Supplementary Note 3.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1.

Description of sensitivity analyses conducted, and the theoretical context of each analysis.

Analyses including proxy cases of Alzheimer’s disease in UKB

The main analysis was conducted using summary data from phase one of the Alzheimer’s disease GWAS, which included only clinically diagnosed cases. To increase the statistical power of our analysis, we re-ran our analysis using summary data from phase three of Alzheimer’s GWAS, which includes proxy cases as defined in UKB (18). The decision to use phase three as a sensitivity analysis rather than as main analysis, was based on the following arguments: i. UK Biobank ‘cases’ included in this meta-analysis were not themselves diagnosed with Alzheimer’s disease (rather, it refers to the participants’ parent), ii. Cases are not doctor diagnosed but based on self-report and iii. The question asks about broad dementia, not specifically Alzheimer’s disease.

RESULTS

Selection of instruments & instrument strength

Out of the 41 cytokines we aimed to explore, 26 of them had at least one genetic variant available at the genome-wide significant threshold p<5×10−08 for use in our MR analyses. For an additional 8 cytokines, genetic variants were available at a more liberal threshold of p<5×10−07. For Alzheimer’s disease 29 genetic variants were available at the genome-wide significant threshold p<5×10−08, 26 of which were available in the cytokines GWAS. Additional information about the genetic instruments used in our analyses can be found in Table S1-S2. Notably, all selected instruments for cytokine concentrations and Alzheimer’s disease demonstrated an F-statistic larger than 10, indicating that weak instruments bias is unlikely to bias our results. Information on the number of instruments identified for each cytokine concentration, the threshold used for selecting instruments and the cumulative F-statistic per cytokine can be found in Table 2. Lastly, four genetic variants were identified as genetic instruments for more than one cytokine and Figure 1 illustrates the overlap. All these variants were included in our main analyses and their influence on our results was further explored in leave – one out analyses.

Figure 1.
  • Download figure
  • Open in new tab
Figure 1.

Graphical presentation of the number of genetic variants overlapping between cytokines

View this table:
  • View inline
  • View popup
Table 2.

Descriptive table of the number of participants for each circulating cytokine concentration, number of instruments identified for each cytokine, the threshold used for selecting instruments and the cumulative F-statistic

Causal effects of circulating cytokine concentrations on risk of Alzheimer’s disease

Overall, there was little evidence to support a causal effect of greater levels of circulating cytokines on risk of Alzheimer’s disease (Figure 2A-2C). Exception was CTACK (CCL27), where we observed weak evidence of a causal effect on Alzheimer’s disease (IVW: OR per 1 standard deviation (SD) increase =1.09 95%CI: 1.01 to 1.19, p=0.03). We also observed weak evidence of a causal effect of 1 SD increase in concentrations of MIP-1b (CCL4) (IVW: OR=1.04 95%CI: 0.99 to 1.09, p=0.08) and Eotaxin (IVW: OR=1.08 95%CI: 0.99 to 1.17, p =0.10) on risk of Alzheimer’s disease. Additionally, weak evidence was observed for an adverse effect of 1 SD increase in levels of GROa (CXCL1) (IVW: OR=1.04 95%CI: 0.99 to 1.10, p=0.15), MIG (CXCL9) (Wald Ratio: OR=1.17 95%CI: 0.97 to 1.41, p=0.10), IL-8 (Wald Ratio: OR=1.21 95%CI: 0.97 to 1.51, p=0.09) and IL-2 (Wald Ratio: OR=1.21 95%CI: 0.94 to 1.56, p=0.14).

Figure 2A.
  • Download figure
  • Open in new tab
Figure 2A.

Total causal effect of genetically predicted cytokine concentrations on the risk of Alzheimer’s disease, as estimated by Wald Ratio, IVW, MR-Egger and Weighted median estimators.

Figure 2B.
  • Download figure
  • Open in new tab
Figure 2B.

Total causal effect of genetically predicted cytokine concentrations on the risk of Alzheimer’s disease, as estimated by Wald Ratio, IVW, MR-Egger and Weighted median estimators.

Figure 2C.
  • Download figure
  • Open in new tab
Figure 2C.

Total causal effect of genetically predicted cytokine concentrations on the risk of Alzheimer’s disease, as estimated by Wald Ratio, IVW, MR-Egger and Weighted median estimators.

Sensitivity analyses

When more than three genetic instruments were available and MR-Egger and Weighted median could be estimated, we observed comparable results to the IVW. More specifically, for the cytokines we observed evidence for a causal effect of cytokines on Alzheimer’s disease, MR-Egger and Weighted median estimators yielded similar effect estimates to the IVW estimator (CTACK: MR-Egger slope: OR=0.99 95%CI: 0.83 to 1.17, p=0.88; Weighted median: OR=1.10 95%CI: 1.01 to 1.20, p =0.03; MIP-1b: MR-Egger slope: OR=1.02 95%CI: 0.95 to 1.11, p =0.58; Weighted median: OR=1.01 95%CI: 0.95 to 1.07, p =0.82; Eotaxin: MR-Egger slope: OR=1.06 95%CI: 0.86 to 1.31, p =0.58; Weighted median: OR=1.11 95%CI: 1.00 to 1.24, p =0.05).

There was no indication of heterogeneity as estimated by Cochran’s Q statistic and after iteratively removing genetic variants results remained virtually the same in the leave-one-out analysis (see Supplementary material Figure S1-S20). Steiger filtering suggested that none of the genetic variants used as instruments for cytokines explained more variation in the outcome, thus the larger sample size of the Alzheimer’s disease GWAS was unlikely to have affected our results and the direction of effects was correctly estimated.

When we restricted our analyses to variants located in closest proximity to the encoding gene of each cytokine, we obtained at least one cis variant for 11 cytokines that were initially included in our main analysis. Overall, results remained virtually the same, but with wider confidence intervals, as a smaller number of instruments was included (Table S3). For greater levels of the TRAIL cytokine we observed weak evidence of a detrimental effect (Wald Ratio: OR=1.23 95%CI: 0.97 to 1.55, p=0.07) on Alzheimer’s disease.

Analyses including proxy cases of Alzheimer’s disease in UKB

Overall, there was limited evidence to suggest a causal effect of circulating cytokine concentrations on risk of Alzheimer’s disease (Supplementary Figures S21A-S21C). As the sample size increased significantly with the addition of proxy-cases, confidence intervals were narrower compared to analyses including only cases of Alzheimer’s disease. For IL-2, a detrimental causal effect on Alzheimer’s disease risk was observed (Wald Ratio: OR per 1 SD increase =1.03 95%CI: 1.00 to 1.06, p=0.04), although the effect estimate was much smaller in magnitude compared to analyses including only diagnosed cased of Alzheimer’s disease.

Causal effects of genetic liability to Alzheimer’s disease on circulating cytokine concentrations

There was limited evidence to suggest a causal effect of genetic liability to Alzheimer’s disease on higher levels of circulating cytokine concentrations (Supplementary Figures 22A-22E). However, weak evidence was observed for a causal effect of genetic liability to Alzheimer’s disease on increased levels of circulating SCGFb (IVW: β= -0.06 95%CI: -0.12 to 0.01, p=0.08).

Sensitivity analyses

Overall, MR-Egger and Weighted median yield similar results to IVW estimators. Additionally, there was no indication of heterogeneity as estimated by Cochran’s Q statistic and leave-one-out analysis did not indicate any genetic variants as influential.

Analyses including proxy cases of Alzheimer’s disease in UKB

When we included the proxy-cases of Alzheimer’s disease in our analysis we observed limited evidence to suggest a causal effect of genetic liability to Alzheimer’s disease on circulating cytokine concentrations (Supplementary Figures S22A-S22E). However, the causal effect of genetic liability to Alzheimer’s disease on SCGFb concentrations was confirmed (IVW: β= -0.74 95%CI: -1.16 to -0.33, p<0.01). Notably, the estimated causal effect was greater in magnitude.

DISCUSSION

Within a two-sample MR framework, we investigated the effect of 34 circulating cytokine concentrations on Alzheimer’s disease risk, and vice versa. We observed some evidence for a detrimental effect of greater levels of CTACK (CCL27), MIP-1b (CCL4), Eotaxin, GROa (CXCL1), MIG (CXCL9), IL-8 and IL-2 on the risk of Alzheimer’s disease. In the reverse direction, we observed some evidence to suggest a causal effect of genetic liability to Alzheimer’s disease only on levels of SCGFb.

Very few studies have examined the causal effects of circulating cytokine concentrations of Alzheimer’s disease risk. Previous two-sample MR studies have investigated the causal role of a few cytokines (TNF-a, IL-18, IL-1ra, IL-6) on risk of Alzheimer’s disease (13-16) and reported limited evidence to support a causal role. In our study we used a larger GWAS of Alzheimer’s disease, thus increasing statistical power to observe causal effects. Despite the increased statistical power, our results were in line with these findings suggesting there is little evidence of a causal effect of circulating cytokines on risk of Alzheimer’s disease. Notably, none of the previous two-sample MR studies have investigated the causal effect of genetic liability of Alzheimer’s disease on increased levels of cytokines.

Chemokines are cytokines which regulate immune cell migration and are thought to be mediators of the peripheral monocytes into the inflamed CNS (30), thus are hypothesized to be involved in the pathogenesis of Alzheimer’s disease (31). CTACK (CCL4), MIG (CXCL9), GROa (CXCL1), MIP-1b (CCL4), Eotaxin (CCL11) and IL-8 belong to the chemokine family and evidence suggest that they could potentially play a role in Alzheimer’s disease. More specifically, CTACK (CCL27) is a chemokine involved in the CNS as it is expressed in the cerebral cortex and limbic regions which are mainly affected in Alzheimer’s disease (32). Previous studies have observed higher levels of circulating CTACK (CCL27) in Alzheimer’s disease patients compared to healthy controls (33, 34). However, future research is required to better characterise the role of CTACK (CCL27) in Alzheimer’s disease aetiology.

MIG (CXCL9) is another chemokine that is considered to play a role in the interplay between neurons and glial cells, and binds onto the CXCR3 receptor which has been previously reported to be involved in the pathogenesis of various CNS conditions (e.g. multiple sclerosis, glioma, bipolar disorder) (35, 36). Regarding Alzheimer’s disease, a cross-sectional study reported substantially higher levels of circulating MIG (CXCL9) in patients with Alzheimer’s disease compared to non-cognitively impaired and mildly-cognitively impaired participants (34). Additionally, a recent case-control study demonstrated evidence for an association between higher levels of MIG and Alzheimer’s disease in a Mexican population (33).

A case-control study found that GROa (CXCL1) is overexpressed in the brains of 23 Alzheimer’s disease patients, with no prior diagnosis of immunological diseases, hypertension, cardiac disease or diabetes, compared to age-matched controls (37). This result is supported by animal studies, where CXCL1 was found to drive the hypermethylation of Tau in the primary cortical neurons of mice (38, 39). This is supporting to our findings and together, it suggests a plausible causal role of CXCL1 in the pathogenesis of Alzheimer’s disease.

Xia et al. reported that receptors of MIP-1b (CCL4) were present on microglia and subpopulation of reactive astrocytes and neurons in brains of patients with Alzheimer’s disease compared to controls, thus they could potentially play a role in the progression of Alzheimer’s disease through glial-glial and glial-neuronal interactions (40). Moreover, higher levels of MIP-1b has been associated with cognitive decline in patients with Alzheimer’s disease (41). Additionally, several studies have reported that expression of CCL4 is greater in brains of HIV-infected patients with dementia compared to HIV-infected patients without dementia, which indicates that CCL4 possibly regulates an inflammatory process that indirectly affects neurons (42-44).

Higher levels of Eotaxin were identified in the cerebrospinal fluid and serum of Alzheimer’s disease cases compared to healthy controls (41, 45, 46). In contrast, in two cohort studies, higher levels of Eotaxin in the plasma were not associated with Alzheimer’s disease progression (47, 48). However, all the studies are of small sample size and thus underpowered to identify associations.

Evidence regarding the role of IL-8 in the Alzheimer’s disease pathogenesis is contradicting. A case-control study observed elevated levels of IL-8 in the CSF of Alzheimer’s disease patients compared to non-demented controls (41), when another case-control study reported significantly lower levels in both serum and CSF of Alzheimer’s disease patients (49). Additionally, IL-8 plasma and CSF levels have been found to be associated with higher levels of p-tau and with higher levels of CSF Aβ1-42, which are hallmarks of Alzheimer’s disease (50). However, evidence stems from studies with limited number of participants and further research is required to elucidate the role of these chemokines in the pathogenesis and progression of Alzheimer’s disease.

Strengths & limitations

The main strength of the present study is the use of two-sample MR which could aid in overcoming the drawbacks of traditional observational epidemiology (i.e. reverse causation, residual confounding), thus allowing the estimation of the effect size of the association between circulating cytokines and Alzheimer’s disease. Additionally, the plethora of cytokines which we could identify instruments for allowed us to explore causal associations between cytokine concentrations and Alzheimer’s disease. Following an agnostic approach could be beneficial in this setting as limited observational evidence exists for most of the studied cytokines. We also used the largest GWAS study of Alzheimer’s disease, which included 24,087 clinically diagnosed cases and thus, our study had statistical power to identify causal effects. The power to detect an odds ratio of 1.1 per 1 – SD increase in circulating cytokine concentration, with a = 0.05 and a coefficient of determination (R2) of 3% to 6%, was between 60% and 80%.

There are a few limitations of this study. As few of the genetic variants used in our MR analyses were used as instruments for more than one cytokine, we cannot exclude the possibility that our results were biased due to horizontal pleiotropic effects. We addressed this issue by estimating the total causal effects with alternative MR methods (MR-Egger and weighted median). Even though results remained virtually the same we could not completely exclude the possibility of pleiotropic effects as these methods are not reliable when a limited number of genetic variants is available, which is the case for many of the cytokines investigated.

Moreover, cytokines are complex in their activity as they act pleiotropically (i.e. single cytokine acts on several different cell types), are redundant in their activity (i.e. same process might be activated by multiple cytokines) and can act either synergistically or antagonistically (51). This complex activity of cytokines raises two main issues in our study. The first issue is the possibility that the observed causal effects are either confounded or mediated by one or more cytokines, and therefore not representing the direct causal effect of each cytokine (i.e. the causal effect after taking into account the causal effect of all other cytokines that might mediate or confound the observed total causal effect). A future direction in deciphering the direct causal effects of cytokines on risk of Alzheimer’s disease would be the employment of multivariable Mendelian randomization, which allows the estimation of direct causal effects (52, 53). MVMR could be used either by including cytokines which belong in the same structural family (i.e. chemokines, interleukins, haemopoietins, chemokines, interferons, TNF family, colony stimulating factors) in a single model or based on their anti/pro-inflammatory activity. However, we were not able to implement MVMR mostly due to the presence of correlation between the estimated genetic effects of the inflammatory markers, as they were obtained from the same participants. Additionally, the limited number of genetic instruments available for each cytokine would have resulted in weak instrument bias and would have limited the statistical power to identify direct effects in the MVMR framework (54). The second issue is that in our analyses we assume no gene-gene or gene-environment interactions and thus modelling the interplay between the multiple inflammatory markers examined is not feasible.

Lastly, observed causal effects might be the consequence of collider bias. Collider bias occurs when conditioning on a variable that is affected by both exposure and outcome and could lead to spurious associations (55, 56). In our study, the GWAS used to extract instruments for cytokine concentrations was adjusted for BMI, which is affected by both cytokines and Alzheimer’s disease, thus it could be considered as a collider. However, this is unlikely to affect our results as individuals on the cytokines GWAS were not selected based on their BMI measurement.

CONCLUSION

In a two-sample MR framework we observed some evidence to support a causal role of cytokines in the pathogenesis of Alzheimer’s disease. More studies are needed to elucidate the specific mechanistic pathways that underlie this process. Better understanding of these processes could potentially lead to novel therapeutic targets for affected individuals.

Data Availability

N/A

Footnotes

  • ↵* joint last authors

  • Funding:This work was supported by a grant from the BRACE Alzheimer’s Disease charity (BR16/028). PP, ELA, and ES work in a unit that receives funding from the University of Bristol and the UK Medical Research Council (MC_UU_00011/1, MC_UU_00011/3). ELA is funded by an MRC Skills Development Award from the UK Medical Research Council (MR/P014437/1). LDH is funded by a Career Development Award from the UK Medical Research Council (MR/M020894/1). This publication is the work of the authors, and ELA, will serve as a guarantor for the contents of this paper.

  • Competing interests: None declared

REGERENCES

  1. 1.↵
    Alzheimer’s A. 2015 Alzheimer’s disease facts and figures. Alzheimers Dement. 2015;11(3):332–84.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Heneka MT, Carson MJ, El Khoury J, Landreth GE, Brosseron F, Feinstein DL, et al. Neuroinflammation in Alzheimer’s disease. Lancet Neurol. 2015;14(4):388–405.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Chou RC, Kane M, Ghimire S, Gautam S, Gui J. Treatment for Rheumatoid Arthritis and Risk of Alzheimer’s Disease: A Nested Case-Control Analysis. CNS Drugs. 2016;30(11):1111–20.
    OpenUrlCrossRefPubMed
  4. 4.↵
    Leira Y, Dominguez C, Seoane J, Seoane-Romero J, Pias-Peleteiro JM, Takkouche B, et al. Is Periodontal Disease Associated with Alzheimer’s Disease? A Systematic Review with Meta-Analysis. Neuroepidemiology. 2017;48(1-2):21–31.
    OpenUrl
  5. 5.↵
    Mrak RE, Griffinbc WS. The role of activated astrocytes and of the neurotrophic cytokine S100B in the pathogenesis of Alzheimer’s disease. Neurobiol Aging. 2001;22(6):915–22.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    Tuppo EE, Arias HR. The role of inflammation in Alzheimer’s disease. Int J Biochem Cell Biol. 2005;37(2):289–305.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Solito E, Sastre M. Microglia function in Alzheimer’s disease. Front Pharmacol. 2012;3:14.
    OpenUrlCrossRefPubMed
  8. 8.↵
    Tejera D, Heneka MT. Microglia in Alzheimer’s disease: the good, the bad and the ugly. Curr Alzheimer Res. 2016;13(4):370–80.
    OpenUrlCrossRef
  9. 9.↵
    Shen XN, Niu LD, Wang YJ, Cao XP, Liu Q, Tan L, et al. Inflammatory markers in Alzheimer’s disease and mild cognitive impairment: a meta-analysis and systematic review of 170 studies. J Neurol Neurosurg Psychiatry. 2019;90(5):590–8.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23(R1):R89–98.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    Davey Smith G, Ebrahim S. What can mendelian randomisation tell us about modifiable behavioural and environmental exposures? BMJ. 2005;330(7499):1076–9.
    OpenUrlFREE Full Text
  12. 12.↵
    Lawlor DA. Commentary: Two-sample Mendelian randomization: opportunities and challenges. Int J Epidemiol. 2016;45(3):908–15.
    OpenUrlCrossRefPubMed
  13. 13.↵
    Andrews SJ, Goate A. Mendelian randomization indicates that TNF is not causally associated with Alzheimer’s disease. 2019:673749.
  14. 14.
    Tsui A, Davis D. Systemic inflammation and causal risk for Alzheimer’s dementia: Possibilities and limitations of a Mendelian randomization approach. Aging Med (Milton). 2018;1(3):249–53.
    OpenUrl
  15. 15.
    Larsson SC, Traylor M, Malik R, Dichgans M, Burgess S, Markus HS. Modifiable pathways in Alzheimer’s disease: Mendelian randomisation analysis. BMJ (Clinical research ed). 2017;359:j5375.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    Prins BP, Abbasi A, Wong A, Vaez A, Nolte I, Franceschini N, et al. Investigating the Causal Relationship of C-Reactive Protein with 32 Complex Somatic and Psychiatric Outcomes: A Large-Scale Cross-Consortium Mendelian Randomization Study. PLoS Med. 2016;13(6):e1001976.
    OpenUrlPubMed
  17. 17.↵
    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(1):40–50.
    OpenUrl
  18. 18.↵
    Jansen IE, Savage JE, Watanabe K, Bryois J, Williams DM, Steinberg S, et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet. 2019;51(3):404–13.
    OpenUrlPubMed
  19. 19.↵
    Korologou-Linden R, Anderson EL, Howe LD, Millard LAC, Ben-Shlomo Y, Williams DM, et al. The causes and consequences of Alzheimer’s disease: phenome-wide evidence from Mendelian randomization. medRxiv. 2020:2019.12.18.19013847.
  20. 20.↵
    Burgess S, Labrecque JA. Mendelian randomization with a binary exposure variable: interpretation and presentation of causal estimates. Eur J Epidemiol. 2018;33(10):947–52.
    OpenUrlCrossRefPubMed
  21. 21.↵
    Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27(8):1133–63.
    OpenUrlCrossRefPubMed
  22. 22.↵
    Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65.
    OpenUrlCrossRefPubMed
  23. 23.↵
    Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, Consortium E-I. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52.
    OpenUrlCrossRefPubMed
  24. 24.
    Burgess S, Thompson SG, Collaboration CCG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.
    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(2):512–25.
    OpenUrlCrossRefPubMed
  26. 26.
    Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40(4):304–14.
    OpenUrlCrossRefPubMed
  27. 27.
    Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.
    OpenUrlCrossRefPubMed
  28. 28.
    Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13(11):e1007081.
    OpenUrlCrossRef
  29. 29.
    Schmidt AF, Finan C, Gordillo-Marañón M, Asselbergs FW, Freitag DF, Patel RS, et al. Genetic drug target validation using Mendelian randomization. 2019:781039.
  30. 30.↵
    Zuena AR, Casolini P, Lattanzi R, Maftei D. Chemokines in Alzheimer’s Disease: New Insights Into Prokineticins, Chemokine-Like Proteins. Front Pharmacol. 2019;10:622.
    OpenUrl
  31. 31.↵
    Guedes JR, Lao T, Cardoso AL, El Khoury J. Roles of Microglial and Monocyte Chemokines and Their Receptors in Regulating Alzheimer’s Disease-Associated Amyloid-beta and Tau Pathologies. Front Neurol. 2018;9:549.
    OpenUrlCrossRefPubMed
  32. 32.↵
    Gunsolly C, Nicholson JD, Listwak SJ, Ledee D, Zelenka P, Verthelyi D, et al. Expression and regulation in the brain of the chemokine CCL27 gene locus. J Neuroimmunol. 2010;225(1-2):82–90.
    OpenUrlPubMed
  33. 33.↵
    Gongora-Rivera F, Gonzalez-Aquines A, Ortiz-Jimenez X, de la Garza CM, Salinas-Carmona M. Chemokine profile in Alzheimer’s disease: Results from a Mexican population. J Clin Neurosci. 2020;73:159–61.
    OpenUrl
  34. 34.↵
    Lee KS, Chung JH, Lee KH, Shin MJ, Oh BH, Hong CH. Bioplex analysis of plasma cytokines in Alzheimer’s disease and mild cognitive impairment. Immunol Lett. 2008;121(2):105–9.
    OpenUrl
  35. 35.↵
    Koper OM, Kaminska J, Sawicki K, Kemona H. CXCL9, CXCL10, CXCL11, and their receptor (CXCR3) in neuroinflammation and neurodegeneration. Adv Clin Exp Med. 2018;27(6):849–56.
    OpenUrl
  36. 36.↵
    Zhou YQ, Liu DQ, Chen SP, Sun J, Zhou XR, Xing C, et al. The Role of CXCR3 in Neurological Diseases. Curr Neuropharmacol. 2019;17(2):142–50.
    OpenUrl
  37. 37.↵
    Zhang K, Tian L, Liu L, Feng Y, Dong YB, Li B, et al. CXCL1 contributes to beta-amyloid-induced transendothelial migration of monocytes in Alzheimer’s disease. PLoS One. 2013;8(8):e72744.
    OpenUrlCrossRefPubMed
  38. 38.↵
    Xia M, Hyman BT. GROα/KC, a chemokine receptor CXCR2 ligand, can be a potent trigger for neuronal ERK1/2 and PI-3 kinase pathways and for tau hyperphosphorylation—a role in Alzheimer’s disease? Journal of Neuroimmunology. 2002;122(1):55–64.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    Shang Y, Tian L, Chen T, Liu X, Zhang J, Liu D, et al. CXCL1 promotes the proliferation of neural stem cells by stimulating the generation of reactive oxygen species in APP/PS1 mice. Biochemical and Biophysical Research Communications. 2019;515(1):201–6.
    OpenUrl
  40. 40.↵
    Xia MQ, Qin SX, Wu LJ, Mackay CR, Hyman BT. Immunohistochemical study of the beta-chemokine receptors CCR3 and CCR5 and their ligands in normal and Alzheimer’s disease brains. Am J Pathol. 1998;153(1):31–7.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    Taipa R, das Neves SP, Sousa AL, Fernandes J, Pinto C, Correia AP, et al. Proinflammatory and anti-inflammatory cytokines in the CSF of patients with Alzheimer’s disease and their correlation with cognitive decline. Neurobiol Aging. 2019;76:125–32.
    OpenUrl
  42. 42.↵
    Minagar A, Shapshak P, Fujimura R, Ownby R, Heyes M, Eisdorfer C. The role of macrophage/microglia and astrocytes in the pathogenesis of three neurologic disorders: HIV-associated dementia, Alzheimer disease, and multiple sclerosis. J Neurol Sci. 2002;202(1-2):13–23.
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.
    Schmidtmayerova H, Nottet HS, Nuovo G, Raabe T, Flanagan CR, Dubrovsky L, et al. Human immunodeficiency virus type 1 infection alters chemokine beta peptide expression in human monocytes: implications for recruitment of leukocytes into brain and lymph nodes. Proc Natl Acad Sci U S A. 1996;93(2):700–4.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    Sui Y, Li S, Pinson D, Adany I, Li Z, Villinger F, et al. Simian human immunodeficiency virus-associated pneumonia correlates with increased expression of MCP-1, CXCL10, and viral RNA in the lungs of rhesus macaques. Am J Pathol. 2005;166(2):355–65.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Choi C, Jeong JH, Jang JS, Choi K, Lee J, Kwon J, et al. Multiplex analysis of cytokines in the serum and cerebrospinal fluid of patients with Alzheimer’s disease by color-coded bead technology. J Clin Neurol. 2008;4(2):84–8.
    OpenUrlCrossRefPubMed
  46. 46.↵
    Soares HD, Potter WZ, Pickering E, Kuhn M, Immermann FW, Shera DM, et al. Plasma biomarkers associated with the apolipoprotein E genotype and Alzheimer disease. Arch Neurol. 2012;69(10):1310–7.
    OpenUrlCrossRefPubMed
  47. 47.↵
    Westin K, Buchhave P, Nielsen H, Minthon L, Janciauskiene S, Hansson O. CCL2 is associated with a faster rate of cognitive decline during early stages of Alzheimer’s disease. PLoS One. 2012;7(1):e30525.
    OpenUrlCrossRefPubMed
  48. 48.↵
    Leung R, Proitsi P, Simmons A, Lunnon K, Guntert A, Kronenberg D, et al. Inflammatory proteins in plasma are associated with severity of Alzheimer’s disease. PLoS One. 2013;8(6):e64971.
    OpenUrlCrossRefPubMed
  49. 49.↵
    Hesse R, Wahler A, Gummert P, Kirschmer S, Otto M, Tumani H, et al. Decreased IL-8 levels in CSF and serum of AD patients and negative correlation of MMSE and IL-1beta. BMC Neurol. 2016;16(1):185.
    OpenUrl
  50. 50.↵
    Bettcher BM, Johnson SC, Fitch R, Casaletto KB, Heffernan KS, Asthana S, et al. Cerebrospinal Fluid and Plasma Levels of Inflammation Differentially Relate to CNS Markers of Alzheimer’s Disease Pathology and Neuronal Damage. J Alzheimers Dis. 2018;62(1):385–97.
    OpenUrl
  51. 51.↵
    Zhang JM, An J. Cytokines, inflammation, and pain. Int Anesthesiol Clin. 2007;45(2):27–37.
    OpenUrlCrossRefPubMed
  52. 52.↵
    Burgess S, Dudbridge F, Thompson SG. Re: “Multivariable Mendelian Randomization: The Use of Pleiotropic Genetic Variants to Estimate Causal Effects”. American Journal of Epidemiology. 2015;181(4):290–1.
    OpenUrlCrossRefPubMed
  53. 53.↵
    Sanderson E, Davey Smith G, Windmeijer F, Bowden J. An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol. 2019;48(3):713–27.
    OpenUrl
  54. 54.↵
    Sanderson E, Spiller W, Bowden J. Testing and Correcting for Weak and Pleiotropic Instruments in Two-Sample Multivariable Mendelian Randomisation. bioRxiv. 2020:2020.04.02.021980.
  55. 55.↵
    Gkatzionis A, Burgess S. Contextualizing selection bias in Mendelian randomization: how bad is it likely to be? Int J Epidemiol. 2019;48(3):691–701.
    OpenUrl
  56. 56.↵
    Paternoster L, Tilling K, Davey Smith G. Genetic epidemiology and Mendelian randomization for informing disease therapeutics: Conceptual and methodological challenges. PLoS Genet. 2017;13(10):e1006944.
    OpenUrlCrossRefPubMed
Back to top
PreviousNext
Posted November 20, 2020.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Causal effects of circulating cytokine concentrations on risk of Alzheimer’s disease: A bidirectional two-sample Mendelian randomization study
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Causal effects of circulating cytokine concentrations on risk of Alzheimer’s disease: A bidirectional two-sample Mendelian randomization study
Panagiota Pagoni, Laura D Howe, George Davey Smith, Yoav Ben-Shlomo, Evie Stergiakouli, Emma L Anderson
medRxiv 2020.11.18.20232629; doi: https://doi.org/10.1101/2020.11.18.20232629
Reddit logo Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Causal effects of circulating cytokine concentrations on risk of Alzheimer’s disease: A bidirectional two-sample Mendelian randomization study
Panagiota Pagoni, Laura D Howe, George Davey Smith, Yoav Ben-Shlomo, Evie Stergiakouli, Emma L Anderson
medRxiv 2020.11.18.20232629; doi: https://doi.org/10.1101/2020.11.18.20232629

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Genetic and Genomic Medicine
Subject Areas
All Articles
  • Addiction Medicine (227)
  • Allergy and Immunology (504)
  • Anesthesia (110)
  • Cardiovascular Medicine (1238)
  • Dentistry and Oral Medicine (206)
  • Dermatology (147)
  • Emergency Medicine (282)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (531)
  • Epidemiology (10020)
  • Forensic Medicine (5)
  • Gastroenterology (499)
  • Genetic and Genomic Medicine (2452)
  • Geriatric Medicine (236)
  • Health Economics (479)
  • Health Informatics (1642)
  • Health Policy (752)
  • Health Systems and Quality Improvement (636)
  • Hematology (248)
  • HIV/AIDS (533)
  • Infectious Diseases (except HIV/AIDS) (11864)
  • Intensive Care and Critical Care Medicine (626)
  • Medical Education (252)
  • Medical Ethics (74)
  • Nephrology (268)
  • Neurology (2280)
  • Nursing (139)
  • Nutrition (352)
  • Obstetrics and Gynecology (454)
  • Occupational and Environmental Health (536)
  • Oncology (1245)
  • Ophthalmology (377)
  • Orthopedics (134)
  • Otolaryngology (226)
  • Pain Medicine (157)
  • Palliative Medicine (50)
  • Pathology (324)
  • Pediatrics (730)
  • Pharmacology and Therapeutics (312)
  • Primary Care Research (282)
  • Psychiatry and Clinical Psychology (2280)
  • Public and Global Health (4832)
  • Radiology and Imaging (837)
  • Rehabilitation Medicine and Physical Therapy (491)
  • Respiratory Medicine (651)
  • Rheumatology (285)
  • Sexual and Reproductive Health (238)
  • Sports Medicine (227)
  • Surgery (267)
  • Toxicology (44)
  • Transplantation (125)
  • Urology (99)