Yalan Li
Jun Lu
Jie Wang
Peizhi Deng
Changjiang Meng1,2
Haibo Tang- 1Center of Clinical Pharmacology, The Third Xiangya Hospital, Central South University, Changsha, China
- 2Department of Cardiology, The Third Xiangya Hospital, Central South University, Changsha, China
- 3School of Medicine, Hunan Normal University, Changsha, China
- 4Department of Metabolic and Bariatric Surgery, The Third Xiangya Hospital, Central South University, Changsha, China
Background: Observational studies have revealed the association between some inflammatory cytokines and the occurrence of ischemic stroke, but the causal relationships remain unclear.
Methods: We conducted a two-sample Mendelian randomization (MR) analysis to assess the causal effects of thirty inflammatory cytokines and the risk of ischemic stroke. For exposure data, we collected genetic variants associated with inflammatory cytokines as instrumental variables (IVs) from a genome-wide association study (GWAS) meta-analysis from Finland (sample size up to 8,293). For the outcome data, we collected summary data of ischemic stroke from a large-scale GWAS meta-analysis involved 17 studies (34,217 cases and 406,111 controls). We further performed a series of sensitivity analyses as validation of primary MR results.
Results: According to the primary MR estimations and further sensitivity analyses, we established one robust association after Bonferroni correction: the odds ratio (95% CI) per unit change in genetically increased IL-4 was 0.84 (0.89–0.95) for ischemic stroke. The chemokine MCP3 showed a nominally significant association with ischemic stroke risk (OR: 0.93, 95% CI: 0.88–0.99, unadjusted p < 0.05). There was no evidence of a causal effect of other inflammatory cytokines and the risk of ischemic stroke.
Conclusions: Our study suggested that genetically increased IL-4 levels showed a protective effect on the risk of ischemic stroke, which provides important new insights into the potential therapeutic target for preventing ischemic stroke.
Introduction
Stroke has become one of the leading causes of death and long-term disability in humans (D 2019 Stroke Collabora, 2021). Ischemic strokes account for approximately 70% of the incidence of stroke and have an extremely high morbidity and mortality rate. Globally, more than nine million people suffer from ischemic stroke each year, placing a vast medical burden on society (Phipps and Cronin, 2020; Powers, 2020), suggesting an urgent need to find new treatment strategies to control stroke. Ischemic stroke results from the disruption of blood supply to the brain due to various causes and the corresponding neurological deficits. In addition to the traditional risk factors such as hypertension, diabetes, smoking, alcohol consumption, and obesity, inflammatory factors play a crucial role in ischemic stroke, with both adverse and beneficial effects (Bonaventura et al., 2016).
Inflammation promotes infarct enlargement and is accountable for its resolution with implications for remodeling and repair (Lambertsen et al., 2019). This conflicting outcome may be due to genetic variation in molecules involved in inflammatory and metabolic pathways (Alfieri et al., 2020). Although experimental evidence suggests that targeting some of these inflammatory cytokines holds promise for treating ischemic stroke (Lambertsen et al., 2019), studies on the connection between circulating inflammatory cytokines and stroke risk remain scarce. Some observational studies on the role of inflammatory factors in the occurrence and recurrence of ischemic stroke have failed to find a causal relationship (Georgakis et al., 2019a).
Mendelian randomization (MR), as a widely used analytical method for causal inference, is a specific case of instrumental variable analysis. When MR assumptions are established, they could be used to identify and quantify causal relationships between exposures and outcomes of interest. This design is not susceptible to confounders and reverse causality bias (Smith and Ebrahim, 2003). In this study, we performed a two-sample MR analysis to test whether there was a causal relationship between inflammatory cytokines and the risk of ischemic stroke.
Methods
Study Design
The workflow of our study is shown in Figure 1. First, we extracted genetic variants as IVs (instrumental variables) for thirty inflammatory cytokines. Second, we collected the summary data including all SNPs from the largest genome-wide association studies (GWASs) for ischemic stroke; Third, we performed two-sample MR analyses with five MR methods [e.g., inverse-variance weighted (IVW), MR-Egger regression, weighted median, MR-Robust adjusted profile score (MR-RAPS), and MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO)]. Fourthly, we conducted a series of sensitivity analyses, including Cochran’s Q Test, Egger intercept test, and MR PRESSO global test to evaluate the heterogeneity and horizontal pleiotropy of MR results.
FIGURE 1. Diagram of Mendelian randomization framework in this study.
Selection of Genetic Predictors of Inflammatory Cytokines
We obtained the genetic predictors from the most comprehensive cytokine-related GWAS meta-analysis for three independent cohorts [the Cardiovascular Risk in Young Finns Study (YFS), FINRISK 1997, and FINRISK 2002], including up to 8,293 Finnish participants (Ahola-Olli et al., 2017). The associations between genetic variants and inflammatory cytokines were adjusted for age, sex, and body mass index. We selected SNPs associated with inflammatory cytokines as IVs using a genome-wide p-value threshold (5 × 10–6, IV assumption 1, Supplementary Figure S1). A stringent condition (linkage disequilibrium threshold of r2 < 0.001 and distance located 10,000 kb apart from each other) was set to ensure that the selection of IVs is conditionally independent of each other. F statistic represents the strength of the relationship between IVs and inflammatory cytokines. Generally, F > 10 may attenuate bias produced by weak IVs (Burgess and Thompson, 2011).
Selection of Ischemic Stroke
We collected summary data of ischemic stroke from a large-scale GWAS meta-analysis containing 17 studies that involved only European participants (34,217 cases and 406,111 controls) (Malik et al., 2018). The largest available GWAS summary statistics were extracted from the MR-Base database (∼8,000,000 genetic variants) (Hemani et al., 2018). The participants had an identical genetic background (all Europeans), and to our knowledge, there was no sample overlap between the exposure and outcome GWASs.
Statistical Analysis
Two-Sample Mendelian Randomization Analysis
As shown in Supplementary Figure S1, we estimated the causal association between inflammatory cytokines and ischemic stroke using a classic MR model: βcausal effect = βZY/βZX (βZX and βZY represent the regression coefficient of SNPs on inflammatory cytokines and ischemic stroke, respectively) (Smith and Ebrahim, 2003; Zheng et al., 2017). MR analysis also relies on three IV assumptions (Supplementary Figure S1) (D 2019 Stroke Collabora, 2021): must be truly associated with inflammatory cytokines (in this study, defined as the genetic association p < 5 × 10–6) (Phipps and Cronin, 2020); not associated with confounders of inflammatory cytokines and ischemic stroke [in this study, we conducted a phenome-wide association test to assess the relationships of IVs with potential confounders such as body mass index, blood pressure, hypertension, and plasma lipid levels using PhenoScanner V2 (Kamat et al., 2019)]; and (Powers, 2020) should only be associated with ischemic stroke through inflammatory cytokines.
To evaluate the causal effects of inflammatory cytokines on the risk of ischemic stroke by combining multiple SNPs, we conducted a two-sample Mendelian randomization (Burgess et al., 2015) analysis using five primary methods, including IVW (Johnson and Uk, 2012), MR-Egger regression (Bowden et al., 2015), weighted median (Bowden et al., 2016), MR-RAPS (Zhao et al., 2020), and MR-PRESSO (Verbanck et al., 2018). The IVW is a conventional method to combine each Wald ratio estimation for multiple SNPs, with the largest statistical power among all MR methods. The weighted median estimator could provide valid estimation even when including 50% of the invalid genetic instruments. With the criterion relaxed, the MR-Egger regression provides a valid estimation even when all genetic instruments are invalid (presence of horizontal pleiotropy across SNPs). Nevertheless, it still requires the InSIDE assumption (Instrument Strength Independent of Direct Effect assumption) to be satisfied (Bowden et al., 2015). Notably, the MR-Egger regression has less power and provides wider confidence intervals. Since we used a relatively higher significant threshold (p < 5 × 10–6) to select genetic variants, we further performed the MR-RAPS to obtain MR estimations using potentially weak instruments. The MR-PRESSO regresses the SNP-outcome estimates against the SNP-exposure estimates to identify outlier SNPs and outputs a corrected MR estimate.
Sensitivity MR Analyses
We examined the heterogeneity of the MR results using the Cochran’s Q-test (Burgess et al., 2017) and evaluated the horizontal pleiotropy by testing whether the intercept in MR-Egger regression and MR-PRESSO global test. Once heterogeneity or horizontal pleiotropy was noted, we recomputed MR estimates after removing the outlier SNPs identified by MR-PRESSO.
Based on the 5 MR methods mentioned above, we took the IVW results as the primary MR estimates and considered the consistency of the results across other MR methods. Here, we defined the evidence for a potential causal effect when the following criteria were met (D 2019 Stroke Collabora, 2021): MR results of IVW passed the multiple comparisons adjusted p-value < 0.0017 (0.05/30) after Bonferroni correction (Phipps and Cronin, 2020); other MR methods showed a similar magnitude and same direction with IVW; and (Powers, 2020) there was no evidence of heterogeneity and horizontal pleiotropy (e.g., P heterogeneity, P intercept, and p-value for MR-PRESSO global test >0.05).
MR analysis was performed in R (version 4.0.3) with R packages “TwoSampleMR” (12), “mr-raps” (Zhao et al., 2020), and “MR-PRESSO” (Verbanck et al., 2018). p values were two-sided, and the statistical significance was set at the adjusted p-value < 0.0017.
Results
Participant Characteristics and Genetic Instruments
The characteristics of the participants from YFS and FINRISK and meta-analysis of GWAS for ischemic stroke are shown in Table 1. We selected 4–20 SNPs, 3–17 SNPs, 3–17 SNPs as instruments for interleukins, chemokines, and other cytokines, respectively (Supplementary Table S1), with an F statistic ranging from 18.7–99.7, reflecting a strong instrument strength for inflammatory cytokines.
TABLE 1. Characteristics of inflammatory cytokines and ischemic stroke datasets.
Estimation of Causal Effects of Inflammatory Cytokines on Ischemic Stroke
Before we conducted MR analysis, we identified an invalid variant rs7088799, which was an instrument variant of IL-10 and IL-12p70, and also associated with body mass index, blood pressure, and hypertension (p < 5 × 10–6) using PhenoScanner V2 (Supplementary Table S2). There was no evidence that any other instrument variants violating the IV assumption 2. Table 2 and Figure 2 show the overall results of MR analysis for increased inflammatory cytokines on the risk of ischemic stroke. There was evidence to support that (p < 0.0017 for IVW after Bonferroni correction) genetically increased IL-4 was associated with a lower risk of ischemic stroke (OR: 0.89, 95% CI: 0.84–0.95). The results from other MR methods showed good consistency with IVW (p < 0.05 in weighted median and MR-RAPS, and MR-Egger showed a similar effect size). Similarly, genetically determined MCP3 showed a negative association with ischemic stroke risk (OR: 0.93, 95% CI: 0.88–0.99) at a nominally significant threshold (p < 0.05). However, there was little evidence to support causal relationships between other inflammatory cytokines and the risk of ischemic stroke.
TABLE 2. Two-sample MR estimations showing the effects of interleukins, chemokines, and other cytokines on the ischemic stroke.
FIGURE 2. The association between genetically determined inflammatory cytokines [(A) interleukin (B) chemokine (C) other cytokines] and the risk of ischemic stroke.
Validation of MR Results Using Sensitivity Analyses
We carried out a series of sensitivity analyses to evaluate the heterogeneity and potential horizontal pleiotropy (Table 2). Cochran’s Q-test showed clear evidence (Ph < 0.05) for the presence of heterogeneity for the MR results of IL-6, IL-16, and TNFα (Table 2). MR-PRESSO identified outliers and provided a corrected estimation. After removing the outlier SNPs, the recomputed MR estimations of IL-6, IL-16 were similar to the results mentioned above, and MR-PRESSO failed to provide any results for TNFα because it proxied by only three SNPs (Supplementary Table S3). The MR-Egger intercept test showed no horizontal pleiotropy (Pintercept > 0.05) for all inflammatory cytokines.
Discussion
In this study, we performed a two-sample MR analysis to test if genetic evidence supported a causal relationship of inflammatory cytokines with the risk of ischemic stroke. Understanding the effect of inflammatory factors in ischemic stroke will help us further insights into how inflammation contributes to the initiation and progression of ischemic stroke. We found that genetically predicted IL-4 levels were negatively associated with ischemic stroke, which adds to epidemiological evidence to the role of inflammatory factor-targeted drug therapy in the prevention and treatment of ischemic stroke.
To our knowledge, there are three studies evaluating the association between partial interleukins and the risk of ischemic stroke so far (Jenny et al., 2019; Lin et al., 2020; Yuan et al., 2020). Jenny et al. reported that IL-6 was linked to an increased incidence of ischemic stroke and no significant associations of IL-8 and IL-10 with ischemic stroke risk (Jenny et al., 2019). Nevertheless, there were unmeasured confounders or misclassification of covariates that could lead to bias (Jenny et al., 2019). In addition, it was shown that IL-6 might promote early clinical deterioration of ischemic stroke, and TNF-α did not play a role in early clinical deterioration (Vila et al., 2000). Lin et al. found that IL-1ra is positively linked to cardioembolic stroke while IL-6 was negatively linked to stroke and coronary artery disease (Yuan et al., 2020). However, another study reported that genetically elevated IL-1Ra, soluble IL-6 receptor (sIL-6R), and C-reactive protein (CRP) levels are not causally associated with ischemic stroke. Our study found no statistically significant effect of IL-6 and IL-1ra on ischemic stroke. This inconsistency in results may be due to differences in sample size and different proportions of stroke subtypes in various populations (Tso et al., 2007). Besides, our study confirmed that increased IL-4 was associated with a reduced risk of ischemic stroke. The possible explanation is that the anti-inflammatory factor IL-4 acts by binding to the IL-4 receptor, enhancing the IL-4 signaling pathway, reducing the incidence of ischemic stroke, and promoting recovery after ischemic stroke (Xiong et al., 2011; Ferreira et al., 2014; Liu et al., 2016; Lively et al., 2016; Chen et al., 2020). Bis et al. reported that IL-1b promotes the ischemic stroke risk (Bis et al., 2008). However, in our study, the association between IL-1b and ischemic stroke was attenuated (OR:1.08, 95%CI: 0.99–1.16, p = 0.07). The lack of significance was probably explained by insufficient statistical power because we noticed that the confidence interval of IL-1b was wider than most other cytokines.
Most current studies focus on the effect of chemokines on stroke rather than ischemic stroke. Early animal experiments suggested elevated MCP3 was associated with stroke, especially in aging mice (Townshend et al., 2015). Extensive research using animal models has shown a vital role for monocyte chemotactic protein-1 (MCP1) in atherogenesis and atheroprogression (Aiello et al., 1999; Inoue et al., 2002; Lin et al., 2014). By binding to the receptor CCR2, MCP1 attracts monocytes by being upregulated under chronic inflammatory conditions. There was an observational study involving 17,180 individuals to determine the relationship between circulating levels of MCP1 and incident stroke in the general population (Georgakis et al., 2019b). It was reported that higher circulating MCP1 increases the risk of stroke. Further MR analysis supported this finding in the MEGASTROKE dataset, however, it showed no statistically significant correlation between MCP1 and ischemic stroke (OR: 1.07, 95% CI: 0.97–1.18, p = 0.17) in the UK Biobank (Georgakis et al., 2019c). In our study, we found no evidence to support that genetically determined MCP1 was associated with the risk of ischemic stroke. This might be partially attributable to different standards for selecting IVs and proportions of stroke subtypes in different populations.
Our study systematically evaluated the association between inflammatory cytokines and ischemic stroke risk using summary-level data from a large-scale GWAS meta-analysis. There are also several limitations to this study. Firstly, we extracted IVs using a relatively higher threshold (p < 5 × 10–6), which weak IVs may bias. However, considering both F-statistics and MR-RAPS results, it appears unlikely that weak IVs could have influenced our findings. Secondly, the association of inflammatory cytokines with stroke risks may not be a linear relationship, and we could not obtain the individual-level data to perform a further non-linear MR analysis. Thirdly, the population in this study is restricted to individuals of European ancestry, which means that our findings may not be extendable to other individuals of ancestry.
Conclusion
In summary, our study supported evidence that genetically determined IL-4 levels are related to the reduced risk of ischemic stroke, suggesting that regulation and intervention of certain inflammatory factors might represent an effective strategy for the future treatment and prevention of ischemic stroke.
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.
Author Contributions
YL and JL: Formal analysis, statistical analysis, and writing the original draft. JW and PD: Formal analysis, visualization. CM: Data collection, data curation. HT: Methodology, writing review and editing, and supervision. All authors contributed to the article and approved the submitted version.
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/fphar.2021.779899/full#supplementary-material
Abbreviations
CCR2, C-C Motif Chemokine Receptor 2; CI, confidence interval; CRP, C-reactive protein; GWAS, genome-wide association study; IL, Interleukin; InSIDE, Instrument Strength Independent of Direct Effect; IVs, instrumental variables; IVW, inverse-variance weighted; MCP, monocyte chemotactic protein; MR, Mendelian randomization; MR-PRESSO, MR-Pleiotropy Residual Sum and Outlier; MR-RAPS, MR-Robust adjusted profile score; OR, odds ratio; sIL-6R, soluble IL-6 receptor; SNPs, Single-nucleotide polymorphisms; TNFα, tumor necrosis factor α; YFS, Young Finns Study.
References
Ahola-Olli, A. V., Würtz, P., Havulinna, A. S., Aalto, K., Pitkänen, N., Lehtimäki, T., et al. (2017). Genome-wide Association Study Identifies 27 Loci Influencing Concentrations of Circulating Cytokines and Growth Factors. Am. J. Hum. Genet. 100, 40–50. doi:10.1016/j.ajhg.2016.11.007
Aiello, R. J., Bourassa, P. A., Lindsey, S., Weng, W., Natoli, E., Rollins, B. J., et al. (1999). Monocyte Chemoattractant Protein-1 Accelerates Atherosclerosis in Apolipoprotein E-Deficient Mice. Arterioscler Thromb. Vasc. Biol. 19, 1518–1525. doi:10.1161/01.atv.19.6.1518
Alfieri, D. F., Lehmann, M. F., Flauzino, T., de Araújo, M. C. M., Pivoto, N., Tirolla, R. M., et al. (2020). Immune-Inflammatory, Metabolic, Oxidative, and Nitrosative Stress Biomarkers Predict Acute Ischemic Stroke and Short-Term Outcome. Neurotox Res. 38, 330–343. doi:10.1007/s12640-020-00221-0
Bis, J. C., Heckbert, S. R., Smith, N. L., Reiner, A. P., Rice, K., Lumley, T., et al. (2008). Variation in Inflammation-Related Genes and Risk of Incident Nonfatal Myocardial Infarction or Ischemic Stroke. Atherosclerosis 198, 166–173. doi:10.1016/j.atherosclerosis.2007.09.031
Bonaventura, A., Liberale, L., Vecchié, A., Casula, M., Carbone, F., Dallegri, F., et al. (2016). Update on Inflammatory Biomarkers and Treatments in Ischemic Stroke. Int. J. Mol. Sci. 17. doi:10.3390/ijms17121967
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian Randomization with Invalid Instruments: Effect Estimation and Bias Detection through Egger Regression. Int. J. Epidemiol. 44, 512–525. doi:10.1093/ije/dyv080
Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet. Epidemiol. 40, 304–314. doi:10.1002/gepi.21965
Burgess, S., Bowden, J., Fall, T., Ingelsson, E., and Thompson, S. G. (2017). Sensitivity Analyses for Robust Causal Inference from Mendelian Randomization Analyses with Multiple Genetic Variants. Epidemiology 28, 30–42. doi:10.1097/EDE.0000000000000559
Burgess, S., Scott, R. A., Timpson, N. J., Davey Smith, G., and Thompson, S. G. (2015). Using Published Data in Mendelian Randomization: a Blueprint for Efficient Identification of Causal Risk Factors. Eur. J. Epidemiol. 30, 543–552. doi:10.1007/s10654-015-0011-z
Burgess, S., and Thompson, S. G. (2011). Avoiding Bias from Weak Instruments in Mendelian Randomization Studies. Int. J. Epidemiol. 40, 755–764. doi:10.1093/ije/dyr036
Chen, X., Zhang, J., Song, Y., Yang, P., Yang, Y., Huang, Z., et al. (2020). Deficiency of Anti-inflammatory Cytokine IL-4 Leads to Neural Hyperexcitability and Aggravates Cerebral Ischemia-Reperfusion Injury. Acta Pharm. Sin B 10, 1634–1645. doi:10.1016/j.apsb.2020.05.002
Ferreira, R., Lively, S., and Schlichter, L. C. (2014). IL-4 Type 1 Receptor Signaling Up-Regulates KCNN4 Expression, and Increases the KCa3.1 Current and its Contribution to Migration of Alternative-Activated Microglia. Front Cel Neurosci 8, 183. doi:10.3389/fncel.2014.00183
GBD 2019 Stroke Collaborators (2021). Global, Regional, and National burden of Stroke and its Risk Factors, 1990-2019: a Systematic Analysis for the Global Burden of Disease Study 2019. Lancet Neurol. 20, 795–820. doi:10.1016/s1474-4422(21)00252-0
Georgakis, M. K., Gill, D., Rannikmäe, K., Traylor, M., Anderson, C. D., Lee, J. M., et al. (2019). Genetically Determined Levels of Circulating Cytokines and Risk of Stroke. Circulation 139, 256–268. doi:10.1161/CIRCULATIONAHA.118.035905
Georgakis, M. K., Malik, R., Björkbacka, H., Pana, T. A., Demissie, S., Ayers, C., et al. (2019). Circulating Monocyte Chemoattractant Protein-1 and Risk of Stroke: Meta-Analysis of Population-Based Studies Involving 17 180 Individuals. Circ. Res. 125, 773–782. doi:10.1161/CIRCRESAHA.119.315380
Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base Platform Supports Systematic Causal Inference across the Human Phenome. Elife 7, 7. doi:10.7554/eLife.34408
Inoue, S., Egashira, K., Ni, W., Kitamoto, S., Usui, M., Otani, K., et al. (2002). Anti-monocyte Chemoattractant Protein-1 Gene Therapy Limits Progression and Destabilization of Established Atherosclerosis in Apolipoprotein E-Knockout Mice. Circulation 106, 2700–2706. doi:10.1161/01.cir.0000038140.80105.ad
Jenny, N. S., Callas, P. W., Judd, S. E., McClure, L. A., Kissela, B., Zakai, N. A., et al. (2019). Inflammatory Cytokines and Ischemic Stroke Risk: The REGARDS Cohort. Neurology 92, e2375–e84. doi:10.1212/wnl.0000000000007416
Kamat, M. A., Blackshaw, J. A., Young, R., Surendran, P., Burgess, S., Danesh, J., et al. (2019). PhenoScanner V2: an Expanded Tool for Searching Human Genotype-Phenotype Associations. Bioinformatics 35, 4851–4853. doi:10.1093/bioinformatics/btz469
Lambertsen, K. L., Finsen, B., and Clausen, B. H. (2019). Post-stroke Inflammation-Target or Tool for Therapy? Acta Neuropathol. 137, 693–714. doi:10.1007/s00401-018-1930-z
Lin, J., Kakkar, V., and Lu, X. (2014). Impact of MCP-1 in Atherosclerosis. Curr. Pharm. Des. 20, 4580–4588. doi:10.2174/1381612820666140522115801
Lin, J., Wang, Y., Wang, Y., and Pan, Y. (2020). Inflammatory Biomarkers and Risk of Ischemic Stroke and Subtypes: A 2-sample Mendelian Randomization Study. Neurol. Res. 42, 118–125. doi:10.1080/01616412.2019.1710404
Liu, X., Liu, J., Zhao, S., Zhang, H., Cai, W., Cai, M., et al. (2016). Interleukin-4 Is Essential for Microglia/Macrophage M2 Polarization and Long-Term Recovery after Cerebral Ischemia. Stroke 47, 498–504. doi:10.1161/strokeaha.115.012079
Lively, S., Hutchings, S., and Schlichter, L. C. (2016). Molecular and Cellular Responses to Interleukin-4 Treatment in a Rat Model of Transient Ischemia. J. Neuropathol. Exp. Neurol. 75, 1058–1071. doi:10.1093/jnen/nlw081
Malik, R., Chauhan, G., Traylor, M., Sargurupremraj, M., Okada, Y., Mishra, A., et al. (2018). Multiancestry Genome-wide Association Study of 520,000 Subjects Identifies 32 Loci Associated with Stroke and Stroke Subtypes. Nat. Genet. 50, 524–537. doi:10.1038/s41588-018-0058-3
Phipps, M. S., and Cronin, C. A. (2020). Management of Acute Ischemic Stroke. BMJ 368, l6983. doi:10.1136/bmj.l6983
Powers, W. J. (2020). Acute Ischemic Stroke. N. Engl. J. Med. 383, 252–260. doi:10.1056/NEJMcp1917030
Smith, G. D., and Ebrahim, S. (2003). 'Mendelian Randomization': Can Genetic Epidemiology Contribute to Understanding Environmental Determinants of Disease. Int. J. Epidemiol. 32, 1–22. doi:10.1093/ije/dyg070
Townshend, C. A., Koellhoffer, E. C., Grenier, J., Chauhan, A., DiMauro, S., Venna, V., et al. (2015).Abstract W P98: Monocyte Chemotactic Protein-3 Elevated in Stoke and Aging, Stroke, 46. AWP98. doi:10.1161/str.46.suppl_1.wp98
Tso, A. R., Merino, J. G., and Warach, S. (2007). Interleukin-6 174G/C Polymorphism and Ischemic Stroke: a Systematic Review. Stroke 38, 3070–3075. doi:10.1161/strokeaha.107.492231
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018). Detection of Widespread Horizontal Pleiotropy in Causal Relationships Inferred from Mendelian Randomization between Complex Traits and Diseases. Nat. Genet. 50, 693–698. doi:10.1038/s41588-018-0099-7
Vila, N., Castillo, J., Dávalos, A., and Chamorro, A. (2000). Proinflammatory Cytokines and Early Neurological Worsening in Ischemic Stroke. Stroke 31, 2325–2329. doi:10.1161/01.str.31.10.2325
Xiong, X., Barreto, G. E., Xu, L., Ouyang, Y. B., Xie, X., and Giffard, R. G. (2011). Increased Brain Injury and Worsened Neurological Outcome in Interleukin-4 Knockout Mice after Transient Focal Cerebral Ischemia. Stroke 42, 2026–2032. doi:10.1161/strokeaha.110.593772
Yuan, S., Lin, A., He, Q. Q., Burgess, S., and Larsson, S. C. (2020). Circulating Interleukins in Relation to Coronary Artery Disease, Atrial Fibrillation and Ischemic Stroke and its Subtypes: A Two-Sample Mendelian Randomization Study. Int. J. Cardiol. 313, 99–104. doi:10.1016/j.ijcard.2020.03.053
Zhao, Q., Wang, J., Hemani, G., Bowden, J., and Small, D. S. (2020). Statistical Inference in Two-Sample Summary-Data Mendelian Randomization Using Robust Adjusted Profile Score. Ann. Statist. 48, 48. doi:10.1214/19-AOS1866
Keywords: Mendelian randomization, inflammatory cytokines, ischemic stroke, causal inference, instrumental variable (IV)
Citation: Li Y, Lu J, Wang J, Deng P, Meng C and Tang H (2022) Inflammatory Cytokines and Risk of Ischemic Stroke: A Mendelian Randomization Study. Front. Pharmacol. 12:779899. doi: 10.3389/fphar.2021.779899
Received: 20 September 2021; Accepted: 06 December 2021;
Published: 17 January 2022.
Edited by:
Cristiano Fava, University of Verona, ItalyReviewed by:
Minoli A. Perera, Northwestern University, United StatesAlexander A. Mongin, Albany Medical College, United States
Copyright © 2022 Li, Lu, Wang, Deng, Meng and Tang. 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: Haibo Tang, thb_9108@csu.edu.cn
†These authors have contributed equally to this work