ORIGINAL RESEARCH article

Front. Nutr., 18 May 2023

Sec. Nutritional Epidemiology

Volume 10 - 2023 | https://doi.org/10.3389/fnut.2023.1103330

Increased selenium and decreased iron levels in relation to risk of coronary artery disease in patients with diabetes

  • 1. Department of Cardiology, The First Affiliated Hospital of Ningbo University, School of Medicine, Ningbo University, Ningbo, China

  • 2. Cardiovascular Disease Clinical Medical Research Center of Ningbo, Ningbo, China

  • 3. Department of Cardiology, Ningbo First Hospital, Ningbo, China

  • 4. Key Laboratory of Precision Medicine for Atherosclerotic Diseases of Zhejiang Province, Ningbo, China

Abstract

Background:

Observational studies have reported inconsistent associations between micronutrient levels and the risk of coronary artery disease (CAD) in diabetic patients. We aim to explore the causal association between genetically predicted concentrations of micronutrients (phosphorus, magnesium, selenium, iron, zinc, and copper) and CAD in patients with diabetes.

Methods:

Single nucleotide polymorphisms (SNPs) connected to serum micronutrient levels were extracted from the corresponding published genome-wide association studies (GWASs). Summary-level statistics for CAD in diabetic patients were obtained from a GWAS of 15,666 patients with diabetes. The primary analysis was carried out with the inverse variance weighted approach, and sensitivity analyses using other statistical methods were further employed to assess the robustness of the results.

Results:

Genetically predicted selenium level was causally associated with a higher risk of CAD in diabetic patients (odds ratio [OR]: 1.25; 95% confidence interval [CI]: 1.10–1.42; p = 5.01 × 10−4). While, genetically predicted iron concentrations in patients with diabetes were inversely associated with the risk of CAD (OR: 0.82; 95% CI: 0.75–0.90; p = 2.16 × 10−5). The association pattern kept robust in most sensitivity analyses. Nominally significant associations were observed for magnesium and copper with the risk of CAD in patients with diabetes. No consistent evidence was found for the causal associations between phosphorus and zinc levels, and the risk of CAD in patients with diabetes.

Conclusion:

We provide consistent evidence for the causal effect of increased selenium and decreased iron levels on CAD in patients with diabetes, highlighting the necessity of micronutrient monitoring and application in these patients.

Introduction

Coronary artery disease (CAD) remains the leading cause of death worldwide, especially in patients with diabetes (, ). Since CAD is responsible for more than 50% of diabetes-related mortality, it dictates the prognosis for diabetic patients (). Therefore, the 2019 European Society of Cardiology (ESC) guidelines have clarified the importance of preventing CAD in patients with diabetes ().

Growing evidence from observational studies and randomized controlled trials (RCTs) indicated that essential micronutrients may play a critical role in the development of CAD in people with diabetes, but the results were inconsistent (). For example, a meta-analysis including 40 prospective cohort studies with over 1 million individuals has shown that increasing dietary magnesium intake was associated with a reduced risk of diabetes and all-cause mortality, but not CAD or total cardiovascular diseases (CVDs) (). Observational studies found the negative or no association between selenium biomarkers and CAD (, ), however, the RCTs revealed that decreased heart disease mortality among individuals with diabetes was related to increased selenium concentration (). On the one hand, as observational studies based on reports of participants were subjected to confounding factors, which might be inaccurate leading to biased results (). On the other hand, due to the limits of the sample size, the evidence from RCTs may not be powerful enough to evaluate the causal effect of micronutrients on the risk of CAD in diabetic patients ().

Mendelian randomization (MR) approach can be applied to explore the potential causal association between exposure and disease by using genetic variants as instrumental variables (IVs) (). The constraints of observational studies are successfully resolved by the random assignment of genotype at conception and the non-influence of genetic variations by potential confounding variables (). In the current study, a two-sample MR analysis was conducted to investigate the causal associations between the risk of CAD in diabetic patients and circulating concentrations of six systematically selected micronutrients, including phosphorus, magnesium, selenium, iron, zinc, and copper.

Methods

Study design

A two-sample MR analysis was designed to estimate the causal relationship between genetically determined circulating micronutrient concentrations and the risk of CAD in diabetic patients (Figure 1). The following three core assumptions should be met by the single nucleotide polymorphisms (SNPs) chosen as IVs for circulating concentrations of micronutrients: (1) IVs should be closely related to the circulating concentrations of micronutrients, (2) IVs should be independent of any potential confounders, and (3) IVs should be associated with the risk of CAD in patients with diabetes only through the concentrations of micronutrients.

Figure 1

Genetic instrument selection

First, SNPs were obtained from recently published genome-wide association studies (GWASs) that independently affect these nutrient concentrations at the genome-wide significance level (p < 5 × 10−8). Then, the linkage disequilibrium tests were performed based on the European 1000 Genomes Project reference panel (r2 < 0.01). If two SNPs were in linkage disequilibrium, the one with smaller value of p would be kept. Considering palindromic SNPs, those with minor allele frequency larger than 0.42 were regarded as not inferable and removed. Specifically, SNPs linked to serum phosphorus levels were extracted from a large GWAS meta-analysis including 16,264 participants of European ancestry (). Six SNPs that achieved genome-wide significance in the joint analysis of the discovery (n = 15,366 participants) and replication (n = 8,463 participants) cohorts from European descent were utilized as genetic IVs for serum magnesium concentration (). The GWAS meta-analysis of log-transformed toenail selenium concentrations and standardized residuals of log-transformed blood selenium concentrations, which included up to 4,162 individuals in four United States studies, provided the genetic summary data for serum selenium levels (). The genetic association with serum iron levels was derived from the Genetics of Iron Status consortium, with up to 48,972 participants (). The GWAS meta-analysis employing two cohorts from Australia and the United Kingdom yielded the SNPs chosen as genetic IVs for zinc and copper concentrations ().

Data for outcome

The summary statistics for CAD in patients with diabetes were extracted from the recently published GWAS, including 15,666 patients of European ancestry with diabetes (3,968 CAD cases and 11,696 controls) from the United Kingdom Biobank (). The average age at diabetes diagnosis was 52.4 ± 12.2 for CAD cases (Male: 2,936; 74.0%) and it was 51.2 ± 12.6 for controls (Male: 7,037; 60.2%). The average age at visit was 62.7 ± 5.6 and 60.2 ± 7.0 for individuals with or without CAD, respectively.

All of the studies in our analyses have obtained relevant ethics review approvals, and all the participants included in the original studies provided written informed permission. All the data used in the current study had been publicly available.

Statistical analysis

The multiplicative random-effects inverse-variance weighted (IVW) method was employed as the primary analysis to evaluate the effect of genetically predicted micronutrient concentrations on the risk of CAD in diabetic individuals. Specifically, the causal estimate for each SNP was generated using the Wald estimator, and the corresponding standard error was calculated using the Delta method. Subsequently, the overall estimate was calculated by meta-analyzing all the estimates by the IVW method ().

To further validate the accuracy of the findings, the Maximum likelihood (), Weighted median (), MR-Egger regression (), and Mendelian Randomization Pleiotropy Residual Sum and Outlier (MR-PRESSO) methods were applied in follow-up sensitivity analyses (). For instance, the Maximum likelihood method could provide a greater empirical power of estimates as it assumed that the genetic association between risk factors and outcomes follows a bivariate normal distribution (). The Weighted median method could still produce reliable estimates even if ≤50% of the weight comes from the ineffective SNPs (). Intercept tests could be used in the MR-Egger regression to assess the potential horizontal pleiotropy (). The MR-PRESSO method was conducted to identify potential outliers and, after eliminating them, to provide relatively unbiased causal estimates (). In addition, scatter plots and leave-one-out analyses were performed to depict the associations of genetically determined micronutrient levels with CAD in patients with diabetes. However, sensitivity analyses and leave-one-out analyses could not be performed as the number of SNPs for zinc and copper was less than three. Cochran’s Q statistics and corresponding value of p were calculated to assess the degree of heterogeneity in the IVW analyses (). Considering the Bonferroni adjustment for multiple tests, a value of p of <0.008 (0.05/6 exposures) was deemed statistically significant. The value of ps between 0.008 and 0.05 were considered to indicate suggested associations. All the statistical analyses were conducted by R Software (version 4.1.1.; R Foundation for Statistical Computing, Vienna, Austria), the R package TwoSampleMR,1 and MR-PRESSO.2

Results

Two to seven SNPs genetically determining the serum phosphorus levels were identified as IVs for serum phosphorus, magnesium, selenium, iron, zinc, and copper levels, respectively (Table 1). In the MR analysis, all F-statistic values of the genetic tools were above the suggested threshold of 10 (Table 1).

Table 1

ExposureSNPChrPosEAOAEAFFMicronutrientsCAD in patients with diabetes
BetaSEp valueBetaSEp value
Phosphorusrs1697421121,823,292CT0.491000.0500.0051.14E−270.0330.0270.217
Phosphorusrs172657033122,048,644GA0.85360.0360.0064.32E−090.0420.0380.259
Phosphorusrs9469578633,706,479TC0.92430.0590.0091.11E−11−0.0210.0520.684
Phosphorusrs9475836136,133,659TC0.29490.0350.0053.45E−120.0360.0300.224
Phosphorusrs2970818124,606,168TA0.09350.0470.0084.38E−09−0.0290.0440.510
Magnesiumrs11144134977,499,796CT0.081210.0110.0018.20E−150.0620.0480.191
Magnesiumrs13146355477,412,140AG0.44250.0050.0016.30E−130.0330.0270.220
Magnesiumrs39255841130,760,335TC0.55360.0060.0015.20E−160.0070.0270.791
Magnesiumrs40720371155,162,067TC0.541000.0100.0012.00E−360.0120.0270.640
Magnesiumrs4483783169,100,899AG0.53160.0040.0011.30E−080.0190.0270.471
Seleniumrs921943579,020,653TC0.291190.2500.0201.90E−390.0690.0290.019
Seleniumrs567754579,120,593CT0.67670.1700.0208.40E−200.0560.0280.048
Seleniumrs3797535579,004,574TC0.10360.2100.0402.10E−15−0.0110.0490.816
Seleniumrs11951068579,008,491AG0.06310.2100.0401.90E−110.0760.0510.137
Seleniumrs705415578,996,137CT0.88390.2300.0404.60E−10−0.0080.0420.841
Seleniumrs65862822143,058,387CT0.85210.1200.0304.00E−09−0.0460.0350.194
Seleniumrs17899532143,062,826TC0.16170.1200.0303.40E−080.0410.0400.301
Ironrs1800562626,093,141AG0.076960.3280.0162.70E−97−0.0540.0490.268
Ironrs1799945626,091,179CG0.15450−0.1890.0101.10E−810.0690.0380.067
Ironrs8557912237,462,936AG0.55807−0.1810.0071.30E−1390.0280.0270.299
Ironrs81772403133,477,701TG0.6795−0.0660.0076.70E−200.0050.0280.860
Ironrs73858047100,235,970AC0.62950.0640.0071.40E−18−0.0030.0280.918
Zincrs21200191575,334,184TC0.79750.2870.0331.60E−180.0210.0340.537
Zincrs1532423886,268,313AG0.37470.1780.0266.40E−12−0.0060.0270.811
Copperrs117555013,691,528GA0.22380.1980.0325.00E−100.0030.0320.922
Copperrs27692641151,344,741GT0.16850.3130.0342.60E−200.0060.0340.866

Characteristics of the single-nucleotide polymorphisms associated with serum micronutrients levels and coronary artery disease in patients with diabetes.

SNP, single-nucleotide polymorphism; Chr, chromosome; Pos, position; EA, effect allele; OA, other allele; EAF, frequency of effect allele; F, F-statistics; SE, standard error; and CAD, coronary artery disease.

The primary findings of MR studies of genetically predicted circulation concentrations of micronutrients with the risk of CAD in individuals with diabetes were displayed in Figure 2. The random-effects IVW results indicated that genetically predisposition to one standard deviation increase in concentrations of serum copper, selenium, and magnesium was linked to 2% (odds ratio [OR], 1.02; 95% CI, 1.02–1.02 p = 3.75 × 10−49), 25% (OR, 1.25; 95% CI, 1.10–1.42; p = 5.01 × 10−4), and 41% (OR, 1.41; 95% CI, 1.14–1.73; p = 1.25 × 10−3) higher risk of CAD in diabetic patients, respectively (Figure 2). An 18% (OR, 0.82; 95% CI, 0.75–0.90, p = 2.16 × 10−5) reduced risk of CAD was observed in patients with diabetes when the genetically predicted serum iron content increase by one standard deviation (Figure 2). There was minimal proof that circulating concentrations of phosphorus and zinc were associated with the risk of CAD in patients with diabetes (Figure 2). The scatter plots also visually depicted the associations between micronutrients and CAD in diabetic patients (Supplementary Figures 1–6).

Figure 2

The association patterns of phosphorus, selenium, and iron based on sensitivity analyses were consistent with the IVW MR analyses, but not magnesium (Figure 2). In addition, stable correlations were found in the MR-PRESSO analysis of serum phosphorus (OR, 1.56; 95% CI, 0.84–2.90; p = 0.23), serum magnesium (OR, 1.34; 95% CI, 1.08–1.66; p = 0.05), serum selenium (OR, 1.25; 95% CI, 1.10–1.42; p = 0.01), and serum iron (OR, 0.82; 95% CI, 0.75–0.90; p = 0.01) with no outliers were revealed (Figure 2). Between the estimates of chosen SNPs, no evidence of heterogeneity for the relationships between micronutrients and CAD in diabetic patients was observed, and neither the MR-Egger intercept test nor the Cochrane’s Q test indicated any possible directional pleiotropy (all p > 0.05; Table 2). Leave-one-out analyses suggested that no single SNP significantly influenced the effect of serum micronutrient levels on CAD in diabetic patients (Supplementary Figures 7–10).

Table 2

MicronutrientsQ valuepQMR-Egger interceptpintercept
Phosphorus3.260.5160.1000.324
Magnesium1.480.8300.0160.691
Selenium4.780.5720.0440.454
Iron1.030.9050.0090.761
Zinc0.320.570NANA
Copper1.76E-040.989NANA

Heterogeneity and pleiotropy tests for the associations of micronutrients with coronary artery disease in patients with diabetes.

MR, Mendelian randomization; NA, not available.

Discussion

In this comprehensive MR analysis, genetic data from the largest published GWAS were leveraged to evaluate the relationship between genetic susceptibility to six micronutrients and the risk of CAD in diabetic patients. We provided consistent evidence that circulating selenium concentrations were genetically expected to be related with a higher risk of CAD, whereas iron concentrations were associated with a lower risk of CAD in patients with diabetes. The association pattern remained consistent when repeated in the majority of supplementary analyses. However, there was limited evidence to link the risk of CAD in diabetic patients with circulating levels of magnesium, phosphorus, zinc, and copper.

Selenium and CAD in patients with diabetes

According to the previous observational studies and RCTs, the association between selenium and CAD in patients with diabetes was inconsistent (, , ). A prospective study involving 3,897 diabetes in the Dongfeng-Tongji cohort suggested an inverse association between plasma levels of selenium and risk of cardiovascular diseases (CVDs) in patients with diabetes (). Selenium supplementation was not sufficient however, to reduce CAD mortality, according to the findings from a meta-analysis that included 16 RCTs (). Additionally, previous observational studies have reported no difference in circulating selenium concentrations between diabetic patients with and without CAD (). A positive association of selenium with diabetes was found in previous observational studies (), several RCTs (), and a MR study (). Numerous in vitro and animal investigations have revealed the mechanism for how selenium increases the risk of diabetes and CAD. As a member of the glutathione peroxidase (GPx) family, selenium serves as the center of redox (). Transgenic animal models have found increased GPx1 expression interferes with insulin signaling by removing hydrogen peroxide, leading to the development of insulin resistance, hyperglycemia, and obesity (, ). Selenoprotein P (SelP), a selenium-supply protein, is hypothesized to raise the risk of diabetes by promoting insulin resistance and dysregulating glucose metabolism (). In addition, Selk, a selenoprotein of the endoplasmic reticulum membrane, contribute to foam cell formation and atherogenesis by stabilizing expression of CD36 in macrophages during inflammation (, ). According to the results of our MR investigation, selenium may be associated in a directionally consistent manner with CAD in patients with diabetes. Given that diabetes is a known risk factor for CAD, this conclusion might seem intuitive. However, considering the majority of the individuals covered with this research were of European origin, the generalizability of our findings to other groups needs to be further investigated. The inverse association between levels of selenium and CVDs risk in Asian diabetic may due to the difference of dietary structure, lifestyle and genetic predisposition.

Iron and CAD in patients with diabetes

An inverse association between iron concentration and CAD in diabetic patients was observed in our MR analysis. A two-sample MR approach examining serum iron status for CAD risk in the general population revealed that serum iron concentration was linked to a lower chance of developing CAD (OR, 0.94; 95%CI, 0.88–1.00; p = 0.039), which is consistent with our findings in the diabetic population (). In addition, a recent two-sample MR study based on the data from United Kingdom Biobank discovered that high levels of iron status were protective against coronary atherosclerosis in the male population (). Furthermore, a meta-analysis of prospective studies involving 156,427 participants showed a negative association between serum iron and risk of coronary heart disease after excluding the study by Morrisson et al. (risk ratio [RR], 0.80; 95%CI, 0.73–0.87) (). Numerous observational studies have also demonstrated the protective effect of iron on CAD in diabetic individuals (, ). An inverse correlation between iron reserves and cardiovascular disease in patients with diabetes was reported by a cross-sectional and prospective observational study encompassing 38,671 people and 821 diabetes patients (OR, 0.81; 95%CI, 0.68–0.96; p = 0.018) (). Similarly, the results of an observational study including 424 consecutive men with type 2 diabetes mellitus showed high ferritin levels may reduce cardiovascular risk in men with diabetes (). Several plausible mechanisms have been hypothesized to elucidate the protective effect of high iron load on CAD. For instance, an animal study found that a high-iron diet attenuates atherosclerosis in mice lacking apolipoprotein E (). Similar, another recent animal study suggested that iron overload could diminish atherosclerosis in apolipoprotein E knockout mice by interfering with hepatic CD36 and fatty acid binding proteins-mediated fatty acid uptake and transport (). Furthermore, it has been proven that ferritin, a natural antioxidant, may reduce the risk of CAD in patients with diabetes by compensating for chronic systemic inflammation in diabetes (, ).

Other micronutrients and CAD in patients with diabetes

In the current MR study, we observed a nominally significant association between genetically predicted concentrations of magnesium and the risk of CAD in patients with diabetes, but the other four statistical models were not statistically significant. As a result, we preclude the presence of a stable causal association between serum magnesium concentration and outcome. We also observed a significant correlation between genetically predicted concentrations of copper and the risk of CAD in patients with diabetes in the main analysis; however, because there are only two genetic instruments for copper, we are unable to perform sensitivity analysis to assess the stability of the results. Meanwhile, the results of IVW (fixed effects) suggested no causal connection between copper and CAD in patients with diabetes. Thus, we are unable to tell whether there is a possible causal relationship between copper and the outcome. The results of the current MR study showed that little evidence approved the causal effects of genetically predicted concentrations of phosphorus and zinc on CAD risk in diabetic patients. There is a scarcity of observational epidemiology research on these micronutrient concentrations and the incidence of CAD in diabetic patients, and the results from the few available observational studies of the general population are inconclusive (, ). Thus, our results from the current MR study may imply that serum phosphorus and zinc levels should not be regarded as independent risk factors for CAD in patients with diabetes.

Strengths and limitations

The design of MR study, which avoids biases frequently seen in standard observational studies and provides the non-biased causal connection between exposure and outcome, is the main merit (). Besides, our MR study uses summary-level data from the large genetic consortium to date, which allows us to more accurately formulate our study hypothesis. Meanwhile, the statistical power in the current investigation is ensured by the estimated effects (F-statistics) of each instrumental variable exceeding the threshold. Moreover, sensitivity analyses based on multiple statistical models combined with leave-one-out analyses were employed to detect the stability of the main results, which offered additional reliable evidence.

It is crucial to acknowledge several potential limitations when interpreting our results. First, although MR-PRESSO analysis and MR-Egger intercept tests did not reveal any evidence of pleiotropy that might have influenced our results, potential horizontal pleiotropy cannot be completely excluded. Second, the current study was based on summary-level data and lacked subgroup-specific analyses, as there are no corresponding sex- or age-specific data sets in the consortium. Third, the majority of the individuals in our MR research were of European origin, which may restrict the generalizability of the primary findings to other groups. Therefore, the corresponding results should be cautious to make the conclusion.

Conclusion

The current study provides genetic evidence for the possible causal effects of increased selenium and decreased iron levels on the increased risk of CAD in patients with diabetes. Diet, supplements, or other methods to modify circulating selenium and iron concentrations may be effective strategies to prevent CAD in patients with diabetes.

Funding

This work was supported by grants from the Key Laboratory of Precision Medicine for Atherosclerotic Diseases of Zhejiang Province, China (Grant No. 2022E10026), National Natural Science Foundation of China (82200489), the Major Project of Science and Technology Innovation 2025 in Ningbo, China (Grant No. 2021Z134), the Key research and development project of Zhejiang Province, China (Grant No. 2021C03096), and Public Science and Technology Projects of Ningbo (202002N3175).

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. The patients/participants provided their written informed consent to participate in this study.

Author contributions

MT, NH, TH, and HC designed the study and wrote the analysis plan. NH and TH undertook analyses. MT and TH wrote the first draft of the manuscript with critical revisions from NH, JY, and HC. MT, TH, JY, HC, and NH interpreted the results in the study and gave final approval of the version to be published. All authors contributed to the article and approved the submitted version.

Acknowledgments

The authors thank all the investigators of CHARGE Consortium, GEFOS consortium, United Kingdom Biobank study, MAGIC, and GISC for providing the data publicly.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

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

References

Summary

Keywords

micronutrient, coronary artery disease, diabetes, causal association, selenium, iron

Citation

Tian M, Hu T, Ying J, Cui H and Huangfu N (2023) Increased selenium and decreased iron levels in relation to risk of coronary artery disease in patients with diabetes. Front. Nutr. 10:1103330. doi: 10.3389/fnut.2023.1103330

Received

29 November 2022

Accepted

03 May 2023

Published

18 May 2023

Volume

10 - 2023

Edited by

Xiaohua Liang, Children‘s Hospital of Chongqing Medical University, China

Reviewed by

Bojana B. Vidovic, University of Belgrade, Serbia; Jong Sung Kim, The University of Iowa, United States

Updates

Copyright

*Correspondence: Ning Huangfu,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics