- Department of Neurology, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
Background: Previous observational studies have shown that circulating selenium levels are inversely associated with ischemic stroke (IS). Our aims were to evaluate the causal links between selenium levels and IS, and its subtypes by Mendelian randomization (MR) analysis.
Methods: We used the two-sample Mendelian randomization (MR) method to determine whether the circulating selenium levels are causally associated with the risk of stroke. We extracted the genetic variants (SNPs) associated with blood and toenail selenium levels from a large genome-wide association study (GWAS) meta-analysis. Inverse variance-weighted (IVW) method was used as the determinant of the causal effects of exposures on outcomes.
Results: A total of 4 SNPs (rs921943, rs6859667, rs6586282, and rs1789953) significantly associated with selenium levels were obtained. The results indicated no causal effects of selenium levels on ischemic stroke by MR analysis (OR = 0.968, 95% CI 0.914–1.026, p = 0.269). Meanwhile, there was no evidence of a causal link between circulating selenium levels and subtypes of IS.
Conclusion: The MR study indicated no evidence to support the causal links between genetically predicted selenium levels and IS. Our results also did not support the use of selenium supplementation for IS prevention at the genetic level.
Introduction
Ischemic stroke (IS) is one of the leading causes of death worldwide and a major cause of serious long-term disability (Campbell et al., 2019). Although IS mortality has been declining globally over the past 2 decades, the number of IS incidents, IS survivors, IS-related deaths, and overall disability-adjusted life years (DALY) lost remains significant and increases year by year (Krishnamurthi et al., 2013). Therefore, early identification of the subjects with a high risk of developing or relapsing IS is of great importance. In addition, the benefit of effective medication for IS (i.e., alteplase) is time-dependent, which limits the wide application of alteplase practice (Phipps and Cronin, 2020). The major challenge of developing new anti-stroke drugs is the presence of the blood–brain barrier and blood circulation gaps, as well as the complexity of signal transduction processes and inflammatory response (Amani et al., 2017; Saxton and Sabatini, 2017). Moreover, fast metabolization clearance from blood circulation and poor transport across the blood–brain barrier hinder the efficacy of most central venous system medications (Amani et al., 2017; Amani et al., 2019). All in all, further investigation of risk factors of IS and targeted therapy strategies is warranted.
The major modifiable risk factors of IS include hypertension, diabetes mellitus, hyperlipidemia, and smoking (Go et al., 2014; Feigin et al., 2016). In addition, some trace elements, particularly essential trace elements, have been reported to be associated with IS (Zecca et al., 2004; Scheiber et al., 2014). Selenium is one of the essential trace elements involved in human physiological processes, metabolism, antioxidant defense, immune regulation, and so on (Burk et al., 2014). The main functions of selenoproteins, the main functional form of selenium, in the neural cells are modulation of neurogenesis, regulation of Ca2+ channels, and maintenance of the redox balance (Cardoso et al., 2015). Reported in vitro studies show that selenium protects mitochondrial functional performance, stimulates mitochondrial biogenesis, and reduces infarct volume after focal cerebral ischemia, through an autophagy-dependent mechanism (Mehta et al., 2012).
Evidence from observational studies indicated that circulating selenium levels were inversely correlated with certain cardiovascular outcomes with a possible U-shaped association, and beneficial effects against IS were found in IS patients as well (Flores-Mateo et al., 2006; Stranges et al., 2010; Rees et al., 2013). However, results from clinical trials were controversial. Specifically, reports of the Selenium and Vitamin E Cancer Prevention Trial (SELECT) and Nutritional Prevention of Cancer Trial (NPC) found no beneficial effects on the incidence and mortality of coronary heart disease and stroke (Stranges et al., 2006; Lippman et al., 2009). In addition, results from a population-based survey revealed that blood selenium concentration might be inversely associated with the prevalence of stroke, and the relationship was non-linear (Hu et al., 2019). However, due to selection bias and reverse causation, the association between selenium levels and the risk of IS may be overestimated. In addition, whether selenium had different impacts on IS subtypes remains unclear. Mendelian randomization (MR), which uses genetic variants as instrumental variables, is a powerful method for inferring causal links between exposures and outcomes. MR analysis uses genetic variants associated with the selenium levels, as the random allocation in randomized controlled trials, to determine the causal effect of the selenium levels on IS, and vice versa (Davies et al., 2018). Since the genes are randomly allocated at conception, genetically predicted selenium levels are not associated with any potential confounders. In addition, random allocation at birth can also avoid the bias caused by reverse causation, as other factors, like disease status cannot affect the genes (Davies et al., 2018). MR analysis was established by three main assumptions (Emdin et al., 2017). First, instrumental variables were significantly associated with the exposure. Next, no links between instrumental variables and confounders were identified. Last, the impact of instrumental variables on outcome was only via exposure (Figure 1). Therefore, MR analysis could overcome the limitations of observational studies and provide insights into the association between selenium and IS. And our aims were to evaluate the causal links between selenium levels and IS and their subtypes by MR analysis.
FIGURE 1. Main assumptions of the Mendelian randomization study of selenium levels and ischemic stroke. IS, ischemic stroke; LVS, large-vessel atherosclerosis stroke; CES, cardio-embolic stroke; SVS, small-vessel occlusion stroke.
Materials and Methods
Data Sources
The genetic variants associated with selenium levels were obtained from a large genome-wide association study (GWAS) meta-analysis of blood selenium (n = 5,477) and toenail selenium (n = 4,162) levels in people of European ancestry (Evans et al., 2013; Cornelis et al., 2015). The genetic variants associated with IS were obtained from a large GWAS by the MEGASTROKE consortium with 34,217 cases and 406,111 controls (Malik et al., 2018). Based on the Trial of ORG 10172 in Acute Stroke Treatment (TOAST) classification, all IS cases could be further divided into large-vessel atherosclerosis stroke (LVS, n = 4,373), cardio-embolic stroke (CES, n = 7,193), and small-vessel occlusion stroke (SVS, n = 5,386) (Adams et al., 1993; Malik et al., 2018). To perform bidirectional MR analysis, the GWAS of the blood selenium level was used as the outcome dataset (Evans et al., 2013).
Sample overlap was calculated in percentages by dividing the number of participants in the GWAS of selenium levels by the number of participants in the respective cohorts in the GWAS of IS and its subtypes (Evans et al., 2013; Cornelis et al., 2015; Malik et al., 2018). An acceptable level of population overlaps between selenium and IS and its subtypes GWAS datasets was 0.22–0.63%.
Selection Criteria of Genetic Variants
We selected genetic variants associated with selenium levels, IS of all causes, LAS, CES, and SVS at genome-wide significance (p < 5 × 10−8) as instrumental variables. Then linkage disequilibrium was tested among the preliminarily selected single-nucleotide polymorphisms (SNPs), and those with r2 > 0.01 in the 1000 Genome Project of Europeans were excluded. The proportion of variance (R2) in the selenium levels explained by the selected genetic variants was calculated using the following formula: R2 = 2 × β2 × (1-EAF) × EAF, where β represents the estimated effect of the genetic variant and EAF represents the effect allele frequency (Palmer et al., 2012). In addition, F-statistic was calculated using the following formula: F = R2 × (N-k-1)/k (1-R2), where R2 represents the proportion of variance explained by the genetic variants, N represents the sample size, and k represents the number of included SNPs (Palmer et al., 2012). The SNPs with an F-statistic <10 were considered weak instruments and were excluded from the MR analysis (Burgess et al., 2011).
Then, the corresponding genetic variants were obtained from the dataset of outcomes (IS or selenium). If selenium-associated SNPs were not available in the outcome datasets, then a proxy SNP in linkage disequilibrium (r2 > 0.9) was searched online (https://ldlink.nci.nih.gov/) as replacement and used in the further analysis.
All genetic variants were searched in the PhenoScanner V2 database to assess whether those variants were significantly associated with the risk factors for IS and its subtypes (Kamat et al., 2019).
Statistical Analysis
All analyses were conducted by R software (version 4.0.3) with R packages TwoSampleMR, MRPRESSO, and MendelianRandomization (Yavorska and Burgess, 2017; Hemani et al., 2018; Verbanck et al., 2018). The estimated effect for blood and toenail selenium levels was presented as Z-score units per effect allele (Evans et al., 2013; Cornelis et al., 2015). Therefore, the Z-score was converted to β and standard error values by the formulas described previously (Kho et al., 2019). The inverse variance-weighted (IVW) method was used as the determinants of the causal effects of exposures on outcomes (Hemani et al., 2018). We also performed MR-Egger, simple median, weighted median, simple mode, weighted mode, robust adjusted profile score (RAPS), Bayesian weighted Mendelian randomization (BWMR), Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO), and Mendelian randomization least absolute shrinkage and selection operator (MR-LASSO) methods (Bowden et al., 2015; Bowden et al., 2016; Hartwig et al., 2017; Verbanck et al., 2018; Zhao et al., 2020). Sensitivity tests including the heterogeneity test (Cochrane’s Q test), pleiotropy test (MR-Egger intercept test), and leave-one-out test were performed (Bowden et al., 2015). Bonferroni correction (corrected p = 0.05/X/Y, where X represents the number of exposures and Y represents the number of outcomes) was used for multiple comparisons.
Power Calculation for Bidirectional Mendelian Randomization Analyses
Statistical power for the bidirectional MR analyses was calculated by mRnd (Brion et al., 2013). The minimum effect estimates of selenium levels required to achieve a power of 80% based on the sample size of the outcome datasets and the R2 by the IVs were calculated and is given in Supplementary Table S1.
Results
The Causal Effects of Selenium Levels on Ischemic Stroke
A total of 4 SNPs (rs921943, rs6859667, rs6586282, and rs1789953) significantly associated with selenium levels were obtained (Table 1). The 4 SNPs explained 5.9% of the variance in the selenium levels, and the corresponding F-statistic was about 151.8. Then, we used PhenoScanner V2 to find whether horizontal pleiotropy existed in the 4 SNPs (Kamat et al., 2019). We found that rs6586282 was significantly associated with plasma homocysteine levels, and rs921943 was associated with height. In MR analysis, the IVW method indicated no causal effects of selenium levels on IS of all causes (OR = 0.968, 95% CI 0.914–1.026, p = 0.269), LVS (OR = 1.015, 95% CI 0.881–1.170, p = 0.835), CES (OR = 1.031, 95% CI 0.922–1.154, p = 0.591), and SVS (OR = 0.984, 95% CI 0.861–1.124, p = 0.811) (Supplementary Table S2 and Figure 2). Heterogeneity tests indicated no heterogeneities of the genetic variants for IS of all causes (p = 0.626), LVS (p = 0.472), CES (p = 0.259), and SVS (p = 0.293) (Supplementary Table S3), and pleiotropy tests indicated no pleiotropy of the genetic variants for IS of all causes (p = 0.896), LVS (p = 0.874), CES (p = 0.669), and SVS (p = 0.802) (Supplementary Table S3). Leave-one-out analysis indicated that the results were still powerful and stable even if they excluded any single SNP (Supplementary Figure S1). Likewise, excluding the effect of rs6586282 did not significantly change the results of MR analysis (Supplementary Figure S1). Altogether, our results indicated no causal effects of selenium levels on IS and its subtypes by MR analysis.
FIGURE 2. Mendelian randomization analysis of the causal effects of selenium levels on ischemic stroke. A total of 4 SNPs significantly associated with selenium levels were obtained. MR, Mendelian randomization; IS, ischemic stroke; SNP, single-nucleotide polymorphism; OR, odds ratio; CI, confidential interval; IVW, inverse variance-weighted; RAPS, robust adjusted profile score; BWMR, Bayesian weighted Mendelian randomization; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier; MR-LASSO, Mendelian randomization least absolute shrinkage and selection operator; LVS, large-vessel atherosclerosis stroke; CES, cardio-embolic stroke; SVS, small-vessel occlusion stroke.
The Causal Effects of Ischemic Stroke on Blood Selenium
To further explore the association between the blood selenium level and IS and its subtypes, we further performed bidirectional MR analysis to estimate the causal effects of IS and its subtypes on blood selenium level. Overall, 9, 4, and 4 SNPs significantly associated with IS of all causes, LVS, and CES were obtained, respectively (Supplementary Table S2). No SNPs significantly associated with SVS were identified. After testing for linkage disequilibrium, 7, 2, and 3 SNPs significantly associated with IS of all causes, LVS, and CES remained, respectively (Table 2; Supplementary Tables S2–S4). By using the IVW method, our results indicated no causal effects of IS of all causes (OR = 0.920, 95% CI 0.622–1.360, p = 0.674), LVS (OR = 1.105, 95% CI 0.620–1.976, p = 0.732), and CES (OR = 0.962, 95% CI 0.787–1.176, p = 0.706) on the blood selenium level (Supplementary Tables S2–S4 and Figure 3). Sensitivity analysis indicated heterogeneities in the analysis of LVS (p = 0.027) and blood selenium level (Table 3). No heterogeneities were identified in the analysis of IS of all cause (p = 0.352) or CES (p = 0.692) (Table 3). The pleiotropy test indicated no pleiotropy (IS of all causes: p = 0.404; CES: p = 0.672) among the genetic variants (Table 3). Leave-one-out analysis indicated that the results of our analysis were powerful (Supplementary Figure S2). Altogether, our results indicated no causal effects of IS and its subtypes on the blood selenium level by MR analysis.
FIGURE 3. Mendelian randomization analysis of the causal effects of ischemic stroke on blood selenium levels. A total of 6, 2, and 3 SNPs significantly associated with IS of all causes, LVS, and CES were obtained in the reverse Mendelian randomization analysis. MR, Mendelian randomization; IS, ischemic stroke; SNP, single-nucleotide polymorphism; OR, odds ratio; CI, confidential interval; IVW, inverse variance-weighted; RAPS, robust adjusted profile score; BWMR, Bayesian weighted Mendelian randomization; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier; MR-LASSO, Mendelian randomization least absolute shrinkage and selection operator; LVS, large-vessel atherosclerosis stroke; CES, cardio-embolic stroke; SVS, small-vessel occlusion stroke.
Discussion
By bidirectional MR analysis based on the summarized data of the GWAS, we found that neither selenium levels were causally associated with IS and its subtypes nor IS and its subtypes were causally associated with selenium levels. The results of our analysis were robust with multiple statistical methods, such as heterogeneity test, pleiotropy test, and leave-one-out analysis.
To our knowledge, the present study is the first study to investigate the causal links between selenium levels and IS and its subtypes by using the bidirectional MR method. Previously, the association between selenium levels and IS was controversial and not well investigated. Prior studies have revealed the potential protective role of selenium in cardiovascular disease. In a case–control study with more than 1,000 Chinese subjects, lower concentrations of selenium were associated with a higher risk of IS (Wen et al., 2019). The inverse association between selenium levels and prevalence of IS was also observed in American subjects (Hu et al., 2019). Nevertheless, Wu et al. (2021) revealed no association between baseline serum selenium levels and stroke in a cohort study (Wei et al., 2004). In a meta-analysis including 12 observational studies, circulating selenium levels were inversely associated with the risk of stroke (Ding and Zhang, 2021). However, in a subgroup analysis, the negative association of selenium levels and stroke was confirmed in the retrospective study group, but not in the prospective study group (Ding and Zhang, 2021). Therefore, the association between selenium levels and IS was controversial and not well investigated. Studies which demonstrated the association between selenium levels and IS with different etiologies were rare. Mironczuk et al. (2021) reported a higher copper-to-selenium ratio in CES patients but a relatively low copper-to-selenium ratio in SVS patients.
The association between selenium levels and stroke is complicated. Selenium is an essential trace element of the human body and shows antioxidant activity by scavenging free radicals (Fang et al., 2002). In the rodent IS model, pretreatment of selenium had significant protective effects on the activity of catalase, superoxide dismutase, and glutathione peroxidase (Ansari et al., 2004). In addition, selenium pretreatment significantly improved hypoxia/ischemia-induced neuron death and reduced infarction volume by alleviating oxidative stress and maintaining mitochondrial function (Mehta et al., 2012). However, the beneficial effect of selenium could be attenuated or even eliminated because of the increasing inflammation and oxidative stress caused by stroke (Ding and Zhang, 2021). Moreover, excess blood selenium concentration (130–150 μg/L) might be associated with minimal mortality (Rayman, 2012).
Gender differences could be a reason for the null finding. Hu et al. (2021) reported a negative association between selenium levels and the first stroke in males but not in females. Different sources (plasma, whole blood, diet, and environment) of selenium used in different studies could be another reason for the null finding and the discrepancy between the present and previous studies (Hu et al., 2017; Merrill et al., 2017; Hu et al., 2019; Wen et al., 2019; Xiao et al., 2019; Hu et al., 2021). Then, regarding the effect of IS on selenium levels, lower selenium levels were observed among acute IS patients in a retrospective study (Angelova et al., 2008). But our analysis provided no evidence of causal effects of IS on selenium levels. Wu et al. (2021) reported genetically predicted selenium levels were negatively causally associated with total cholesterol and low-density lipoprotein cholesterol, which were risk factors for IS (Diener and Hankey, 2020). Furthermore, selenium was reported to be positively correlated with systemic arterial function (Chan et al., 2012). Because previous studies reported non-linear association (including J-shaped and U-shaped) between selenium levels and stroke, the links between selenium levels and IS are rather complicated and still need further investigation (Bleys et al., 2008; Hu et al., 2017; Hu et al., 2019; Hu et al., 2021).
Given the antioxidant activity of selenium and selenoproteins, selenium supplementation was proposed as a potential strategy for the prevention of multiple disorders, like IS, osteoarthritis, rheumatoid arthritis, hypothyroidism, and prostate cancer (Sanmartin et al., 2011). Regarding stroke, selenium supplementation directly into the brain induced the expression of antioxidant glutathione peroxidase 4, which further inhibited the ferroptosis of neurons in a brain hemorrhage model (Alim et al., 2019). In a clinical trial of 29,584 Chinese people, the group receiving selenium supplements for a period of 5 years had a reduction in stroke mortality (9%), but no statistical significance was identified (Mark et al., 1998). Through a secondary analysis of the Nutritional Prevention of Cancer Trial, Stranges et al demonstrated no beneficial effect of selenium supplementation on stroke or cardiovascular disease incidence (Stranges et al., 2006). By bidirectional MR analysis, our results did not support the effectiveness of selenium supplementation in the prevention of IS and its subtypes at the genetic level. Given the impact of selenium levels on blood lipids and arterial function (Chan et al., 2012; Wu et al., 2021), the efficacy of selenium supplementation in subjects with hyperlipidemia or atherosclerotic lesions needed further investigation.
There were some limitations to our study. First, only subjects with European ancestry were included in the MR analysis. The prevalence and incidence of IS vary with ethnicity and so do the proportions of the subtypes of IS (Kim and Kim, 2014). Studies of Western populations indicated CES was the most common subtype of IS, while studies in Asian countries reported a higher prevalence of LVS than CES (Kolominsky-Rabas et al., 2001; Tsai et al., 2013). And the ethnicity differences among the SNPs associated with selenium levels also exist (Supplementary Table S5). Therefore, the results of this study needed further validation in Asian or African people. Second, despite including the genetic variants significantly associated with selenium levels from the largest GWAS of selenium levels, only 4 SNPs were finally included in MR analysis. While the 4 SNPs explained approximately 5.9% of the variance of selenium levels and the F-statistic of each SNP was more than 10. Therefore, more genetic variants associated with selenium levels, both blood and toenail selenium levels, need to be identified in the future. Third, pleiotropy, which is inevitable in MR analysis, may overestimate the effect of the exposure on the outcome. To eliminate the impact of pleiotropy as much as possible, we sought to identify potential pleiotropic SNPs before the MR analysis. By PhenoScanner, we found one SNP significantly associated with homocysteine. In addition, we performed a pleiotropy test by MR-Egger intercept, and no pleiotropy was found in the present study. Fourth, regarding outcome datasets of selenium levels, only blood selenium levels were used in the MR analysis. So, the causal effects of IS and its subtypes on toenail selenium levels are still unclear. Last, although our analysis suggested no effectiveness of selenium supplementation for patients with IS at the genetic level, large randomized controlled trials are needed to investigate the efficacy and safety of selenium supplementation for IS patients.
Conclusion
In conclusion, our bidirectional MR study provides no evidence to support the causal links between genetically predicted selenium levels and IS. Our results also did not support the use of selenium supplementation for IS prevention at the genetic level. Clinical trials with high quality and large sample size are warranted to further elucidate the underlying association between selenium levels and IS and the clinical benefit of selenium supplementation for the prevention of IS.
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
YX and HF conceived and designed the study. WL and LZ collected the data. HF, WL and LP analyzed the data. YG, LZ and RZ interpreted the results. JY, BS and YX supervised the study. HF and WL wrote the manuscript. All authors approved the submitted version.
Funding
The present study was funded by the Education Department of Henan Province (No. 172102310086) and the Education Department of Henan Province (No. 17A320060).
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.
Acknowledgments
The authors thank Cornelis et al. and the MEGASTROKE consortium for providing summarized data of GWAS (Evans et al., 2013; Cornelis et al., 2015; Malik et al., 2018).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.782691/full#supplementary-material
Supplementary Figure S1 | Sensitivity analysis of the causal effects of selenium levels on ischemic stroke. (A) Leave-one-out analysis of the causal effects of selenium levels on ischemic stroke and its subtypes. (B) Scatter plot analysis of the causal effects of selenium levels on ischemic stroke and its subtypes. (C) Forest plot analysis of the causal effects of selenium levels on ischemic stroke and its subtypes. (D) Funnel plot analysis of the causal effects of selenium levels on ischemic stroke and its subtypes. IS, ischemic stroke; LVS, large-vessel atherosclerosis stroke; CES, cardio-embolic stroke; SVS, small-vessel occlusion stroke; MR, Mendelian randomization; SNP, single-nucleotide polymorphism.
Supplementary Figure S2 | Sensitivity analysis of the causal effects of ischemic stroke on blood selenium levels. (A) Leave-one-out analysis of the causal effects of ischemic stroke and its subtypes on blood selenium levels. (B) Scatter plot analysis of the causal effects of ischemic stroke and its subtypes on blood selenium levels. (C) Forest plot analysis of the causal effects of ischemic stroke on blood selenium levels. (D) Funnel plot analysis of the causal effects of ischemic stroke on blood selenium levels. IS, ischemic stroke; LVS, large-vessel atherosclerosis stroke; CES, cardio-embolic stroke; SVS, small-vessel occlusion stroke; MR, Mendelian randomization; SNP, single-nucleotide polymorphism.
References
Adams, H. P., Bendixen, B. H., Kappelle, L. J., Biller, J., Love, B. B., Gordon, D. L., et al. (1993). Classification of Subtype of Acute Ischemic Stroke. Definitions for Use in a Multicenter Clinical Trial. TOAST. Trial of Org 10172 in Acute Stroke Treatment. Stroke 24, 35–41. doi:10.1161/01.str.24.1.35
Alim, I., Caulfield, J. T., Chen, Y., Swarup, V., Geschwind, D. H., Ivanova, E., et al. (2019). Selenium Drives a Transcriptional Adaptive Program to Block Ferroptosis and Treat Stroke. Cell 177, 1262–1279. doi:10.1016/j.cell.2019.03.032
Amani, H., Habibey, R., Hajmiresmail, S. J., Latifi, S., Pazoki-Toroudi, H., and Akhavan, O. (2017). Antioxidant Nanomaterials in Advanced Diagnoses and Treatments of Ischemia Reperfusion Injuries. J. Mater. Chem. B 5, 9452–9476. doi:10.1039/c7tb01689a
Amani, H., Habibey, R., Shokri, F., Hajmiresmail, S. J., Akhavan, O., Mashaghi, A., et al. (2019). Selenium Nanoparticles for Targeted Stroke Therapy through Modulation of Inflammatory and Metabolic Signaling. Sci. Rep. 9, 6044. doi:10.1038/s41598-019-42633-9
Angelova, E. A., Atanassova, P. A., Chalakova, N. T., and Dimitrov, B. D. (2008). Associations between Serum Selenium and Total Plasma Homocysteine during the Acute Phase of Ischaemic Stroke. Eur. Neurol. 60, 298–303. doi:10.1159/000157884
Ansari, M. A., Ahmad, A. S., Ahmad, M., Salim, S., Yousuf, S., Ishrat, T., et al. (2004). Selenium Protects Cerebral Ischemia in Rat Brain Mitochondria. Bter 101, 73–86. doi:10.1385/BTER:101:1:73
Bleys, J., Navas-Acien, A., and Guallar, E. (2008). Serum Selenium Levels and All-Cause, Cancer, and Cardiovascular Mortality Among US Adults. Arch. Intern. Med. 168, 404–410. doi:10.1001/archinternmed.2007.74
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
Brion, M.-J. A., Shakhbazov, K., and Visscher, P. M. (2013). Calculating Statistical Power in Mendelian Randomization Studies. Int. J. Epidemiol. 42, 1497–1501. doi:10.1093/ije/dyt179
Burgess, S., Thompson, S. G., and Crp Chd Genetics Collaboration, (2011). Avoiding Bias from Weak Instruments in Mendelian Randomization Studies. Int. J. Epidemiol. 40, 755–764. doi:10.1093/ije/dyr036
Burk, R. F., Hill, K. E., Motley, A. K., Winfrey, V. P., Kurokawa, S., Mitchell, S. L., et al. (2014). Selenoprotein P and Apolipoprotein E Receptor‐2 Interact at the Blood‐brain Barrier and Also within the Brain to Maintain an Essential Selenium Pool that Protects against Neurodegeneration. FASEB j. 28, 3579–3588. doi:10.1096/fj.14-252874
Campbell, B. C. V., De Silva, D. A., Macleod, M. R., Coutts, S. B., Schwamm, L. H., Davis, S. M., et al. (2019). Ischaemic Stroke. Nat. Rev. Dis. Primers 5, 70. doi:10.1038/s41572-019-0118-8
Cardoso, B. R., Roberts, B. R., Bush, A. I., and Hare, D. J. (2015). Selenium, Selenoproteins and Neurodegenerative Diseases. Metallomics 7, 1213–1228. doi:10.1039/c5mt00075k
Chan, Y.-H., Siu, C.-W., Yiu, K.-H., Chan, H.-T., Li, S.-W., Tam, S., et al. (2012). Adverse Systemic Arterial Function in Patients with Selenium Deficiency. J. Nutr. Health Aging 16, 85–88. doi:10.1007/s12603-011-0086-5
Cornelis, M. C., Fornage, M., Foy, M., Xun, P., Gladyshev, V. N., Morris, S., et al. (2015). Genome-wide Association Study of Selenium Concentrations. Hum. Mol. Genet. 24, 1469–1477. doi:10.1093/hmg/ddu546
Davies, N. M., Holmes, M. V., and Davey Smith, G. (2018). Reading Mendelian Randomisation Studies: a Guide, Glossary, and Checklist for Clinicians. BMJ 362, k601. doi:10.1136/bmj.k601
Diener, H.-C., and Hankey, G. J. (2020). Primary and Secondary Prevention of Ischemic Stroke and Cerebral Hemorrhage. J. Am. Coll. Cardiol. 75, 1804–1818. doi:10.1016/j.jacc.2019.12.072
Ding, J., and Zhang, Y. (2021). Relationship between the Circulating Selenium Level and Stroke: A Meta-Analysis of Observational Studies. J. Am. Coll. Nutr., 1–9. doi:10.1080/07315724.2021.1902880
Emdin, C. A., Khera, A. V., and Kathiresan, S. (2017). Mendelian Randomization. Jama 318, 1925–1926. doi:10.1001/jama.2017.17219
Evans, D. M., Zhu, G., Dy, V., Heath, A. C., Madden, P. A. F., Kemp, J. P., et al. (2013). Genome-wide Association Study Identifies Loci Affecting Blood Copper, Selenium and Zinc. Hum. Mol. Genet. 22, 3998–4006. doi:10.1093/hmg/ddt239
Fang, Y.-Z., Yang, S., and Wu, G. (2002). Free Radicals, Antioxidants, and Nutrition. Nutrition 18, 872–879. doi:10.1016/s0899-9007(02)00916-4
Feigin, V. L., Roth, G. A., Naghavi, M., Parmar, P., Krishnamurthi, R., Chugh, S., et al. (2016). Global burden of Stroke and Risk Factors in 188 Countries, during 1990-2013: a Systematic Analysis for the Global Burden of Disease Study 2013. Lancet Neurol. 15, 913–924. doi:10.1016/S1474-4422(16)30073-4
Flores-Mateo, G., Navas-Acien, A., Pastor-Barriuso, R., and Guallar, E. (2006). Selenium and Coronary Heart Disease: a Meta-Analysis. Am. J. Clin. Nutr. 84, 762–773. doi:10.1093/ajcn/84.4.762
Go, A. S., Mozaffarian, D., Roger, V. L., Benjamin, E. J., Berry, J. D., Blaha, M. J., et al. (2014). Heart Disease and Stroke Statistics-2014 Update. Circulation 129, e28. doi:10.1161/01.cir.0000441139.02102.80
Hartwig, F. P., Davey Smith, G., and Bowden, J. (2017). Robust Inference in Summary Data Mendelian Randomization via the Zero Modal Pleiotropy assumption. Int. J. Epidemiol. 46, 1985–1998. doi:10.1093/ije/dyx102
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, e34408. doi:10.7554/eLife.34408
Hu, H., Bi, C., Lin, T., Liu, L., Song, Y., Wang, B., et al. (2021). Sex Difference in the Association between Plasma Selenium and First Stroke: a Community-Based Nested Case-Control Study. Biol. Sex. Differ. 12, 39. doi:10.1186/s13293-021-00383-2
Hu, X. F., Sharin, T., and Chan, H. M. (2017). Dietary and Blood Selenium Are Inversely Associated with the Prevalence of Stroke Among Inuit in Canada. J. Trace Elem. Med. Biol. 44, 322–330. doi:10.1016/j.jtemb.2017.09.007
Hu, X. F., Stranges, S., and Chan, L. H. M. (2019). Circulating Selenium Concentration Is Inversely Associated with the Prevalence of Stroke: Results from the Canadian Health Measures Survey and the National Health and Nutrition Examination Survey. Jaha 8, e012290. doi:10.1161/JAHA.119.012290
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
Kho, P. F., Glubb, D. M., Thompson, D. J., Spurdle, A. B., and O'Mara, T. A. (2019). Assessing the Role of Selenium in Endometrial Cancer Risk: A Mendelian Randomization Study. Front. Oncol. 9, 182. doi:10.3389/fonc.2019.00182
Kim, B. J., and Kim, J. S. (2014). Ischemic Stroke Subtype Classification: an Asian Viewpoint. J. Stroke 16, 8–17. doi:10.5853/jos.2014.16.1.8
Kolominsky-Rabas, P. L., Weber, M., Gefeller, O., Neundoerfer, B., and Heuschmann, P. U. (2001). Epidemiology of Ischemic Stroke Subtypes According to TOAST Criteria. Stroke 32, 2735–2740. doi:10.1161/hs1201.100209
Krishnamurthi, R. V., Feigin, V. L., Forouzanfar, M. H., Mensah, G. A., Connor, M., Bennett, D. A., et al. (2013). Global and Regional burden of First-Ever Ischaemic and Haemorrhagic Stroke during 1990-2010: Findings from the Global Burden of Disease Study 2010. Lancet Glob. Health 1, e259–e281. doi:10.1016/S2214-109X(13)70089-5
Lippman, S. M., Klein, E. A., Goodman, P. J., Lucia, M. S., Thompson, I. M., Ford, L. G., et al. (2009). Effect of Selenium and Vitamin E on Risk of Prostate Cancer and Other Cancers. JAMA 301, 39–51. doi:10.1001/jama.2008.864
Malik, R., Chauhan, G., Chauhan, G., Traylor, M., Sargurupremraj, M., Okada, Y., 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
Mark, S. D., Wang, W., Mark, J. F., Fraumeni, J. F., Li, J.-Y., Taylor, P. R., et al. (1998). Do nutritional Supplements Lower the Risk of Stroke or Hypertension? Epidemiology 9, 9–15. doi:10.1097/00001648-199801000-00005
Mehta, S. L., Kumari, S., Mendelev, N., and Li, P. A. (2012). Selenium Preserves Mitochondrial Function, Stimulates Mitochondrial Biogenesis, and Reduces Infarct Volume after Focal Cerebral Ischemia. BMC Neurosci. 13, 79. doi:10.1186/1471-2202-13-79
Merrill, P. D., Ampah, S. B., He, K., Rembert, N. J., Brockman, J., Kleindorfer, D., et al. (2017). Association between Trace Elements in the Environment and Stroke Risk: The Reasons for Geographic and Racial Differences in Stroke (REGARDS) Study. J. Trace Elem. Med. Biol. 42, 45–49. doi:10.1016/j.jtemb.2017.04.003
Mirończuk, A., Kapica-Topczewska, K., Socha, K., Soroczyńska, J., Jamiołkowski, J., Kułakowska, A., et al. (2021). Selenium, Copper, Zinc Concentrations and Cu/Zn, Cu/Se Molar Ratios in the Serum of Patients with Acute Ischemic Stroke in Northeastern Poland-A New Insight into Stroke Pathophysiology. Nutrients 13, 2139. doi:10.3390/nu13072139
Palmer, T. M., Lawlor, D. A., Harbord, R. M., Sheehan, N. A., Tobias, J. H., Timpson, N. J., et al. (2012). Using Multiple Genetic Variants as Instrumental Variables for Modifiable Risk Factors. Stat. Methods Med. Res. 21, 223–242. doi:10.1177/0962280210394459
Phipps, M. S., and Cronin, C. A. (2020). Management of Acute Ischemic Stroke. BMJ 368, l6983. doi:10.1136/bmj.l6983
Rayman, M. P. (2012). Selenium and Human Health. The Lancet 379, 1256–1268. doi:10.1016/S0140-6736(11)61452-9
Rees, K., Hartley, L., Day, C., Flowers, N., Clarke, A., and Stranges, S. (2013). Selenium Supplementation for the Primary Prevention of Cardiovascular Disease. Cochrane Database Syst. Rev. 2013, CD009671. doi:10.1002/14651858.CD009671.pub2
Sanmartin, C., Plano, D., Font, M., and Palop, J. A. (2011). Selenium and Clinical Trials: New Therapeutic Evidence for Multiple Diseases. Cmc 18, 4635–4650. doi:10.2174/092986711797379249
Saxton, R. A., and Sabatini, D. M. (2017). mTOR Signaling in Growth, Metabolism, and Disease. Cell 168, 960–976. doi:10.1016/j.cell.2017.02.004
Scheiber, I. F., Mercer, J. F. B., and Dringen, R. (2014). Metabolism and Functions of Copper in Brain. Prog. Neurobiol. 116, 33–57. doi:10.1016/j.pneurobio.2014.01.002
Stranges, S., Marshall, J. R., Trevisan, M., Natarajan, R., Donahue, R. P., Combs, G. F., et al. (2006). Effects of Selenium Supplementation on Cardiovascular Disease Incidence and Mortality: Secondary Analyses in a Randomized Clinical Trial. Am. J. Epidemiol. 163, 694–699. doi:10.1093/aje/kwj097
Stranges, S., Navas-Acien, A., Rayman, M. P., and Guallar, E. (2010). Selenium Status and Cardiometabolic Health: State of the Evidence. Nutr. Metab. Cardiovasc. Dis. 20, 754–760. doi:10.1016/j.numecd.2010.10.001
Tsai, C.-F., Thomas, B., and Sudlow, C. L. M. (2013). Epidemiology of Stroke and its Subtypes in Chinese vs white Populations: a Systematic Review. Neurology 81, 264–272. doi:10.1212/WNL.0b013e31829bfde3
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
Wei, W.-Q., Abnet, C. C., Qiao, Y.-L., Dawsey, S. M., Dong, Z.-W., Sun, X.-D., et al. (2004). Prospective Study of Serum Selenium Concentrations and Esophageal and Gastric Cardia Cancer, Heart Disease, Stroke, and Total Death. Am. J. Clin. Nutr. 79, 80–85. doi:10.1093/ajcn/79.1.80
Wen, Y., Huang, S., Zhang, Y., Zhang, H., Zhou, L., Li, D., et al. (2019). Associations of Multiple Plasma Metals with the Risk of Ischemic Stroke: A Case-Control Study. Environ. Int. 125, 125–134. doi:10.1016/j.envint.2018.12.037
Wu, Q., Sun, X., Chen, Q., Zhang, X., and Zhu, Y. (2021). Genetically Predicted Selenium Is Negatively Associated with Serum TC, LDL-C and Positively Associated with HbA1C Levels. J. Trace Elem. Med. Biol. 67, 126785. doi:10.1016/j.jtemb.2021.126785
Xiao, Y., Yuan, Y., Liu, Y., Yu, Y., Jia, N., Zhou, L., et al. (2019). Circulating Multiple Metals and Incident Stroke in Chinese Adults. Stroke 50, 1661–1668. doi:10.1161/STROKEAHA.119.025060
Yavorska, O. O., and Burgess, S. (2017). MendelianRandomization: an R Package for Performing Mendelian Randomization Analyses Using Summarized Data. Int. J. Epidemiol. 46, 1734–1739. doi:10.1093/ije/dyx034
Zecca, L., Youdim, M. B. H., Riederer, P., Connor, J. R., and Crichton, R. R. (2004). Iron, Brain Ageing and Neurodegenerative Disorders. Nat. Rev. Neurosci. 5, 863–873. doi:10.1038/nrn1537
Keywords: selenium, stroke, trace element, cause, Mendelian randomization (MR)
Citation: Fang H, Liu W, Zhang L, Pei L, Gao Y, Zhao L, Zhang R, Yang J, Song B and Xu Y (2022) A Bidirectional Mendelian Randomization Study of Selenium Levels and Ischemic Stroke. Front. Genet. 13:782691. doi: 10.3389/fgene.2022.782691
Received: 24 September 2021; Accepted: 07 March 2022;
Published: 13 April 2022.
Edited by:
Liangcai Zhang, Janssen Research and Development, United StatesReviewed by:
Kei Hang Katie Chan, City University of Hong Kong, Hong Kong SAR, ChinaCheryl D. Cropp, Samford University, United States
Copyright © 2022 Fang, Liu, Zhang, Pei, Gao, Zhao, Zhang, Yang, Song and Xu. 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: Yuming Xu, eHV5dW1pbmdAenp1LmVkdS5jbg==
†These authors have contributed equally to this work