The Causal Effects of Blood Iron and Copper on Lipid Metabolism Disease: Evidence from Phenome-wide Mendelian Randomization Study
=================================================================================================================================

* Jingqi Zhou
* Chang Liu
* Michael Francis
* Yitang Sun
* Moon-Suhn Ryu
* Arthur Grider
* Kaixiong Ye

## Abstract

**Background** Blood levels of iron and copper, even within their normal ranges, have been associated with a wide range of clinical outcomes. Available epidemiological evidence on blood iron and copper association with potential clinical effects, such as lipid metabolism disorder, is inconsistent and scarce. This study aims to examine and disentangle the causal clinical effects of iron and copper.

**Methods** Genetic instruments for the blood level of iron and copper were curated from existing genome-wide association studies. Candidate clinical outcomes were identified based on a phenome-wide association study (PheWAS) between these genetic instruments and multiple phenotypes in 310,999 unrelated individuals of European ancestry from the UK Biobank. All signals passing stringent correction for multiple testing were followed by Mendelian randomization (MR) analyses, with replication in independent data sources where possible.

**Results** Genetically predicted higher blood levels of iron and copper are both associated with lower risks of iron deficiency anemia (OR = 0.75, 95% CI:0.67-0.85, *p* =1.90×10−06 for iron; OR= 0.88, 95% CI: 0.78-0.98, *p* =0.03 for cooper), lipid metabolism disorders and its two subcategories, hyperlipidemia (OR = 0.90, 95% CI:0.85-0.96, *p* =6.44×10−04 for iron; OR = 0.92, 95% CI:0.87-0.98, *p* = 5.51×10−03 for cooper) and hypercholesterolemia (OR = 0.90, 95% CI:0.84-0.95, *p* = 5.34×10−04 for iron; OR = 0.93, 95% CI:0.89–0.99, *p* =0.02 for cooper). Consistently, they are also associated with lower blood levels of total cholesterol and low-density lipoprotein cholesterol. Multiple sensitivity tests were applied to assess the pleiotropy and stability of the estimation, and consistent evidence across different approaches was obtained. Additionally, unique clinical effects of each blood mineral were also pinpointed, and sex-stratified MR analysis further revealed the clinical implication of iron and copper exhibits some degree of sex differences.

**Conclusions** Our comparative PheWAS-MR study of iron and copper comprehensively characterized their shared and unique clinical effects, highlighting their novel causal roles in hyperlipidemia and hypercholesterolemia. Given the modifiable and variable nature of blood mineral status, these findings warrant further investigation.

Keywords
*   iron
*   copper
*   phenome-wide association study
*   Mendelian randomization
*   lipid metabolism disorder

## Introduction

Iron (Fe) and copper (Cu) are two essential mineral nutrients for human health through its vital role in enzymatic reactions and cellular energy metabolism 1, 2. Physiological processes that maintain these minerals homeostasis are overwhelmed by dietary loss, inappropriate iron absorption, inflammation, infection, liver disease, and disordered erythropoiesis. Disturbance to the systemic iron or copper homeostasis, including deficiency, excess, and even fluctuations in the normal ranges, could have clinical implications3. Iron deficiency, the most widespread micronutrient deficiency worldwide, is well-established to cause anemia, while iron overload can lead to chronic liver disease, cirrhosis, and hepatocellular carcinoma 4, 5. Copper deficiency is responsible for the lethal Menkes syndrome, featured by developmental delay and intellectual disability, while copper overload results in Wilson’s disease, characterized by hepatic cirrhosis and neurological deficits 6. Furthermore, copper is essential for iron absorption, and copper deficiency impairs iron metabolism. The multiple interactions of copper and iron involving both the positive and negative aspects by which associated with various of disease 7. For instance, epidemiological studies have reported elevated blood levels of iron and copper are associated with a higher risk of type 2 diabetes anemia and osteoarthritis (OA) 8, 9. Higher iron and copper have been associated with a reduced risk of hyperlipidemia and an improved blood lipid profile 10, 11. However, conflicting associative patterns were also reported 12. Given the crucial role of minerals, studying the health risks of iron and copper has important clinical guiding significance. Due to the limited data available and the fact that most studies were observational, it is uncertain whether the associations are causal and not biased by reverse causality or residual confounding. Moreover, the causality of iron and copper effect on overall health has not been clearly established.

Mendelian randomization (MR), a complementary approach to epidemiological observations, utilizes genetic variants as instrumental variables to approximate the lifetime status of an exposure (e.g., the blood level of a mineral) and evaluates its causal effect on a clinical outcome. The random allocation of alleles at conception and the natural direction of causality from genetic variants to phenotypes protect MR estimates from confounding and reverse causality 13, 14. MR studies have been performed to evaluate the causality of specific mineral-outcome pairs, such as iron on stroke, coronary artery disease (CAD), Parkinson’s disease, gout, and inflammatory bowel disease 15, 16; copper on ischemic heart disease 17. Compared to these hypothesis-driven MR studies with an obvious bias towards cardiovascular diseases, phenome-wide association study coupled with MR (PheWAS-MR) enables an unbiased and hypothesis-free scan through a wide range of phenotypes (i.e., phenome) and prioritized candidate clinical outcomes for MR causal inferences. Most importantly, few existing studies simultaneously examine multiple blood minerals at a phenome-wide scale to disentangle their confounded clinical effects 18.

In this study, we systematically evaluated and compared the causal clinical effects of iron and copper. Genetic instruments for the blood level of iron and copper were curated from existing genome-wide association studies (GWAS) with the largest sample sizes. Candidate clinical outcomes were identified based on a phenome-wide association study between these genetic instruments and multiple disease outcomes from the UK Biobank. All signals passing stringent correction for multiple testing were followed by Mendelian randomization (MR) analyses, with replication in independent data sources where possible. Additionally, we examined the causal association between iron or copper levels and a subset of lipid profiles, which are essential biomarkers for lipid metabolism disease, in the secondary analysis with an independent dataset.

## Methods

The guideline of Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) (S1 Checklist) was followed in this study report. The UK Biobank project was approved by the North West Multi-Centre Research Ethics Committee and appropriate informed consent was obtained from participants. Data used in the project was accessed through an approved application to UK Biobank. Only a subset of unrelated White British individuals with high-quality genotype data (N = 310,999) were included in this study in order to minimize the complication of population structure. The study design of PheWAS-MR, replication, and sensitivity analysis was prespecified. A schematic of the overall study design is shown in Fig 1.

![Fig 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.17.20196956/F1.medium.gif)

[Fig 1:](http://medrxiv.org/content/early/2020/09/22/2020.09.17.20196956/F1)

Fig 1: 
A flow chart of the study design.

### Genetic instruments for blood iron and copper

Genetic instrumental variables for the blood levels of iron and copper were selected from previous GWAS (Table 1). Based on a meta-analysis of 48,972 subjects of European descent performed by the Genetics of Iron Status consortium 19, three single nucleotide polymorphisms (SNPs) were selected as genetic instruments for systemic iron status. All three of them are consistently associated with four iron status biomarkers, including serum iron, ferritin, transferrin, and transferrin saturation. Their effects on serum iron were used in MR analysis. They are very strong genetic instruments, with F-statistics ranging from 445 to 796 and collectively explaining approximately 3.8% of the variation in serum iron. They have been previously used in MR analysis of iron status 20. Genetic variants associated with erythrocyte copper levels were identified in an Australian cohort (N = 2603) 21. Two SNPs were selected for erythrocyte copper, which accounted for 4.60% of the phenotypic variance. For each blood mineral, a weighted polygenic risk score (GRS) was constructed for each UK Biobank participant by summing over selected SNPs weighted by their effect sizes.

View this table:
[Table 1:](http://medrxiv.org/content/early/2020/09/22/2020.09.17.20196956/T1)

Table 1: 
SNPs used as genetic instruments for blood minerals in the PheWAS and MR analyses.

### Study population

The UK Biobank is a population-based prospective cohort established to study genetic and environmental determinants of human diseases. More than 500,000 individuals aged 40-69 years were recruited between 2006 and 2010, all of whom gave written consent and underwent baseline measurements. Participants donated blood samples for genotyping and biomarker analysis. Moreover, UK Biobank participants consented access to their electronic health records. Only participants fulfilling the following criteria were included in our analysis: genetic ancestry of Caucasian, included in genetic principal component analysis, not outliers for heterogeneity and missing genotype rate, no sex chromosome aneuploidy, no high degree of genetic kinship (i.e., ten or more third-degree relatives identified), and self-reported sex matching genetic sex. Additionally, for remaining pairs of relatives (kinship coefficient > 0.0884), a minimum number of participants were removed so that all remaining participants are unrelated. Analysis was restricted to participants of Caucasian descent in order to maintain consistency with the European samples used to obtain genetic instruments for blood mineral status. A total of 310,999 unrelated individuals were included in our final analysis.

### Phenome wide association study

Three phenotypic datasets in the UK Biobank, including inpatient hospital records, cancer registry, and death registry, were included in this analysis. First, we used the International Classification of Diseases (ICD) versions 9 and 10 to identify cases in the hospital episode statistics, including both incident and prevalent cases, but not self-reported diagnoses. As ICD-9 / ICD-10 codes are not organized as independent phenotypes, we mapped them to the phecode system of distinct diseases or traits 22. The mapping process also automatically excludes patients with similar or potentially overlapping diseases from the corresponding control group. The online phecode map is accessible via the link [https://phewascatalog.org/phecodes_icd10](https://phewascatalog.org/phecodes_icd10). We restricted analysis to phecodes with sufficient numbers of cases to ensure more than 80% statistical power in subsequent MR analyses. The number of cases and controls in analyses of phenotypes for each blood mineral is shown in S1 Table.

For each mineral, case-control groups were generated for each phecode, and logistic regression was performed, separately, for each instrument SNP, adjusting for age, sex, genotyping array, assessment center, and the first 10 genetic principal components. Moreover, we explored the association between the weighted GRS of each blood mineral and the corresponding phecodes with logistic regression. Given that phenotypes investigated are not totally independent even in the phecode system, we applied the false discovery rate (FDR) method with a cutoff of 0.05 to correct for multiple testing.

### MR analyses

To assess whether blood mineral levels have causal effects on the candidate clinical outcomes identified by the PheWAS, we further conducted two-sample MR analyses, where the instrument-mineral and instrument-outcome associations were estimated in two different samples. In this study, the instrument-mineral effect was estimated in previous GWAS, while the instrument-outcome effect was estimated in the above-mentioned PheWAS. To estimate the causal effect of each mineral on a phenotype, for each genetic instrument, we calculated the ratio of its effect on the phenotype to its effect on the blood mineral. An overall effect estimate across multiple genetic instruments was achieved with three methods, the inverse variance-weighted (IVW) method, weighted median (WM) method, and MR-Egger 23, 24. Effect estimates, measured as OR of the outcome, were normalized to 1 standard deviation (SD) increment in each blood mineral. Additionally, sex-stratified IVW MR estimates were obtained with sex-specific instrument-outcome effects from genetically males and females, respectively. Mineral-outcome pairs exhibiting evidence of horizontal pleiotropy were excluded in this stratified analysis.

To validate some mineral-outcome relationships that we discovered based on UK Biobank, we performed replication analysis based on different cohorts, which were accessed through the MR-Base 25. Specifically, for the four blood lipids, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL), low-density lipoprotein cholesterol (LDL), and triglyceride, we extracted the instrument-outcome summary statistics from a previous GWAS by the Global Lipids Genetics Consortium (GLGC) 26.

To further evaluate if the effects of two minerals on a shared clinical outcome are independent of each other, we adopted a multivariable MR framework 27, in which the effects of multiple related risk factors can be estimated simultaneously. For instance, to test if the effects of iron and copper on a shared clinical outcome (e.g., hyperlipidemia) are independent, we included all the genetic instruments for iron and copper in the analysis. The effect sizes of instrument SNPs on serum iron and erythrocyte copper were extracted from previous GWAS 19, 21, while their effects on the shared clinical outcome were calculated in our curated UK Biobank dataset. We modelled the SNP effects on the shared outcome (βshared) against the SNP effects on iron (βFe), adjusted for the effects on copper (βCu), using a weighted linear regression model (βshared∼βFe + βCu), where the weights were defined by inverse standard errors of βshared.

### Sensitivity analyses

Several sensitivity analyses were performed to detect and correct for the presence of pleiotropy in the causal estimates. We began with the fixed-effect IVW MR approach, and calculated Cochran’s Q statistic for heterogeneity, with which a p-value ≤ 0.05 indicates the presence of pleiotropy across genetic instruments. Heterogeneity across SNP-outcome associations will lead to bias in the fixed-effect IVW estimate. If substantial heterogeneity exists, we proceeded to a random-effect IVW approach, which allows instruments to exhibit balanced horizontal pleiotropy.

Because the IVW estimate is essentially a weighted average of the Wald ratios obtained from each SNP, if any of the SNPs exhibit horizontal pleiotropy, then the causal effect estimate is liable to be biased. We further performed WM MR 24 and MR-Egger analyses 28. The WM analysis calculates the median of an empirical distribution of MR association estimates weighted by their precisions. It provides consistent estimates when more than half of the instruments are valid 24. The MR-Egger method provides an intercept test, with a non-zero intercept indicating the presence of unbalanced horizontal pleiotropy. It also provides an unbiased estimate of the causal effect, taking into account pleiotropy. When there was no heterogeneity across instruments, a fixed-effect model was used to estimate the causal effect. Otherwise a random-effect model was used 28.

Besides the above-mentioned sensitivity analyses that leverage information from multiple genetic instruments for one blood mineral, for each individual genetic instrument, we also evaluated their pleiotropic effects by searching the GWAS catalog and PhenoScanner V2 29, to identify any secondary phenotypes associated (*p* < 5 × 10−8) with them or their proxies (r2 > 0.8). All analyses were conducted using the MendelianRandomisation 30 and TwoSampleMR 25 packages, and the R programming language.

## Results

After quality controls on the full UK Biobank dataset, a subset of 310,999 unrelated individuals of European descent were included in the PheWAS analysis. Main demographic characteristics of the study population and basic statistics of all genetic instruments are summarized in S2 Table. All instrument SNPs satisfied the Hardy-Weinberg equilibrium test (*p* > 0.05). The association between the weighted GRS of each blood mineral and common confounding factors (i.e., age, sex, genotyping array, assessment center, and the genetic principal components) are presented in S3 Table. All these factors were controlled for in the PheWAS analysis.

PheWAS analysis was limited to phenotypes (i.e., phecodes) that had enough cases to ensure more than 80% statistical power in the subsequent MR analyses. Among them, 12 clinical outcomes, representing 6 disease categories, reached statistical significance at the 5% FDR threshold. Next, we performed MR analyses to examine the possible causal links between each of the significant outcomes identified in PheWAS. Mineral-outcome pairs exhibit consistent MR evidence for causal effects, without suggestions of heterogeneity and horizontal pleiotropy across genetic instruments based on Cochran’s Q statistic and MR-Egger’s intercept test were provided in test (S4 Tables).

### Shared and unique causal clinical effects of blood iron and copper

Iron and copper share some of causally associated clinical outcomes, including iron deficiency anemia, disorders of lipoid metabolism, hypercholesterolemia, and hyperlipidemia (Fig 2). IVW MR analyses provided significant evidence for causal effects of both higher blood iron and copper on reduced risks of iron deficiency anemia (OR per SD increase in serum iron = 0.75, 95%CI: 0.67-0.85, *p* =1.81×10−6; OR per SD increase in erythrocyte copper = 0.88, 95% CI: 0.79-0.99, *p* = 0.032) and disorders of lipoid metabolism (OR per SD increase in serum iron = 0.90, 95% CI: 0.85-0.96, *p* =6.61×10−4; OR per SD increase in erythrocyte copper = 0.92, 95% CI: 0.87-0.98, *p* = 4.94×10−3;), and its two subtypes, hyperlipidemia (OR = 0.90, CI: 0.85-0.96, *p* = 6.44×10−4; OR = 0.92, CI: 0.87-0.98, *p* = 5.51×10−3) and hypercholesterolemia (OR = 0.90, CI: 0.84-0.95, p = 5.34×10−4; OR = 0.93, CI: 0.88–0.99, *p* = 0.022). Specifically, unique causal clinical effects were identified for iron and copper. Per SD increase in serum iron is associated with enhanced risk of varicose veins (OR = 1.28, CI: 1.15-1.42, *p* = 4.34×10−6) and acquired foot deformities (OR = 1.21, CI: 1.09-1.35, *p* = 4.95×10−4), but with reduced risk of other anemia (OR = 0.72; CI: 0.65-0.79; *p* = 1.12×10−11). Per SD increment in erythrocyte copper is associated with increased risk of osteoarthrosis (OA) (OR = 1.07, CI: 1.02-1.13, *p* = 0.010).

![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.17.20196956/F2.medium.gif)

[Fig 2.](http://medrxiv.org/content/early/2020/09/22/2020.09.17.20196956/F2)

Fig 2. 
A forest plot showing significant mineral-outcome associations based on MR analysis. The causal estimates are from IVW MR and have no indications of pleiotropy. The odds ratios (ORs) with their 95% confidence intervals (CIs) are scaled to 1-SD increase in blood iron or copper level. Complete MR results are provided in S4 Tables.

### Causality of both Iron and Cooper on Lipid Metabolism Trait

We performed further investigation into the shared clinical effects of serum iron and erythrocyte copper on lipid metabolism disorders, by first evaluating if they are independent of each other and then by examining their causal effects on blood biomarkers of lipid metabolism, including TC, HDL, LDL, and triglyceride. First, we adopted a multivariable MR framework to estimate the direct effects of serum iron and erythrocyte copper on their shared clinical outcomes, by conditioning on the effects of the other mineral. We found that the causal effects of iron and copper are independent of each other and do not show much differences compared to those of univariate MR (S5 Table). Based on the shared but independent effects of iron and coper on lipid metabolism disorders, we further hypothesized that their genetically predicted levels are associated with blood biomarkers of lipid metabolism. We tested this hypothesis with the UK Biobank data with exclusion of participants under statin medication, and then in an independent replication cohort from the Global Lipids Genetics Consortium (GLGC) accessed through the MR-base database (Table 2). In both cohorts, genetically predicted higher levels of serum iron is associated with decreased levels of LDL (WM OR = 0.91, CI: 0.89-0.94, *p* = 5.27×10−11 in UKBB) and TC (WM OR = 0.91, CI: 0.88-0.94, *p* = 8.99×10−8 in UKBB), but elevated levels of triglycerides (WM OR = 1.05, CI: 1.02-1.07, *p* = 4.2×10−5 in UKBB). No significant effects were found for HDL. For erythrocyte copper, while no significant effects were found in UK Biobank, analysis in the GCLC cohort showed that per SD increase in copper is associated with reduced levels of LDL (OR = 0.95, CI: 0.92-0.99, *p* = 0.011 in GCLC) and TC (OR = 0.96, CI: 0.92-0.99, *p* =0.02 in GCLC). Overall, we found evidence for the shared effects of iron and copper on reducing blood levels of LDL and TC.

View this table:
[Table 2:](http://medrxiv.org/content/early/2020/09/22/2020.09.17.20196956/T2)

Table 2: 
MR analyses of iron and copper on blood lipids in UK Biobank and GLGC

### Interpretation of potential pleiotropy in PheWAS-MR results

To detect and correct for any possible pleiotropic effect of multiple instruments, we then conducted the different sets of sensitivity analysis. All causal clinical effects of iron and copper presented above in combined analysis, except LDL and TC, do not have evidence for heterogeneity or horizontal pleiotropy across genetic instruments, as indicated by Cochran’s Q statistic and MR-Egger’s intercept test (S4 Tables). For some relationships, causal estimates from another methods, WM MR and MR-Egger, are broadly consistent with those from IVW MR, although with wider confidence intervals as expected from their lower statistical powers 31. For instance, the causal estimates of serum iron on disorders of lipoid metabolism were consistent in the IVW MR (OR = 0.90, CI: 0.85-0.96, *p* = 6.61×10−4) and WM MR (OR= 0.91, 95% CI:0.85-0.97, *p* = 6.25×10−3), and there was no evidence of pleiotropy (*p* ≥ 0.05 for all related tests).

Besides sensitivity analyses that leverage information from multiple genetic instruments for one blood mineral, for each individual instrument SNP, we also evaluated their pleiotropic effects by searching the GWAS catalog and PhenoScanner V2 29. The HFE rs1800562 variant is associated with low-density lipoprotein cholesterol 26, and this SNP were identified as pleiotropic with potential confounders. Similar magnitudes of associations were observed when these variants were excluded from the MR analyses. We additionally performed a sensitivity analysis of sex-stratified IVW MR to identify sex-specific causal clinical effects of these minerals. We identified mineral-outcome pairs exhibiting sex-specific causal relationships (Fig 3, Supplementary Tables 6). That is, their MR estimates are significant only in one sex group, but not in the other or the sex-combined sample. Iron is associated with higher risks of myeloproliferative diseases (OR =2.20, CI: 1.37-3.53, *p* =1.14×10−3) in females and chronic liver disease (OR =1.47, CI: 1.07-2.03, *p* =0.017) in males. Specifically, copper increases the risk of diabetes mellitus (OR = 1.20, CI: 1.02-1.40, *p* = 0.027) and its subcategory, T2D (OR = 1.19, CI: 1.00-1.40, *p* = 0.044) in female. But serum iron and copper have no sex-specific effects on lipid metabolism disease.

![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/22/2020.09.17.20196956/F3.medium.gif)

[Fig 3.](http://medrxiv.org/content/early/2020/09/22/2020.09.17.20196956/F3)

Fig 3. 
A forest plot showing significant mineral-outcome associations based on sex-stratified MR analysis. Results in male, female, and the combined analyses are shown. The causal estimates are from IVW MR and have no indications of pleiotropy. The odds ratios (ORs) with their 95% confidence intervals (CIs) are scaled to 1-SD increase in blood mineral level. Complete MR results are provided in S6 Tables.

## Discussion

Our study adopted a phenome-wide approach to systematically evaluate and compare the clinical effects of blood iron and copper minerals. This is the first PheWAS-MR study for copper, and the first to perform systematic comparison across iron and copper. Our findings highlighted shared clinical effects, specifically the protective effects of iron and copper on lipid metabolism disorders and iron deficiency anemia. Some causal mineral-outcome relationships identified in this study have been well known and supported by existing mechanistic studies, while others are novel, awaiting further confirmation and future mechanistic exploration.

Most notably, we found that genetically predicted higher blood levels of iron and copper are both protective against lipid metabolism disorders and its two subcategories, hyperlipidemia and hypercholesterolemia. Consistently, they are also associated with lower blood TC and LDL. This is the first MR study to establish the causal protective effect of copper on lipid metabolism disorders. The interplays among lipid, iron, and copper metabolism have been recognized, but they are complicated and far away from being fully elucidated. Copper deficiency have long been linked to increased risks of hyperlipidemia and cardiovascular diseases in both humans and animal models 32, while copper supplementation in patients of hyperlipidemia was shown to improve the blood lipid profile 33. Epidemiological findings, however, are conflicting, with some suggesting negative associations of serum copper with both TC and LDL 34, and others showing positive associations 35. The relationship between iron status and blood lipids is similarly convoluted. On the one hand, some epidemiological studies found that serum iron is lower in hyperlipidemic patients 36, and that serum ferritin is negatively associated with LDL 11. Very likely related, elevated iron status is negatively associated with coronary artery disease in both observational and MR studies 37, 38. Moreover, studies on rats showed that iron deficiency upregulates lipogenic genes but downregulates apolipoprotein H and genes involved in the mitochondrial beta-oxidation, resulting in increased circulating lipids 39. On the other hand, anemic and/or iron-deficient patients and animal models tend to have lower blood TC and LDL 40, while serum ferritin has often been found to be positively associated with an unhealthy lipid profile (i.e., elevated TC, LDL, and triglyceride, but decreased HDL) 41. To gain mechanistic insights into how iron regulates lipid metabolism, we performed bioinformatic search for iron responsive elements in genes involved in cholesterol metabolism and fatty acid degradation 42. Multiple potential iron-responsive genes were identified, including LIPH and LDLRAP1 (S1-S2 Figs). Our study highlights iron and copper as risk factors of dyslipidemia and calls for future studies into their physiological mechanism.

Well-known mineral-outcome relationships with established physiological mechanisms provide support for the power and validity of our study. Iron is an integral component of the oxygen-carrying hemoglobin and is required for the process of erythropoiesis 4. In keeping with this, we found that genetically predicted higher level of serum iron is protective against anemia, while in women it increases the risks of myeloproliferative disease, polycythemia vera, and secondary polycythemia. Myeloproliferative diseases are bone marrow and blood disorders featured by abnormal hematopoiesis, while its most common subtype, polycythemia vera is characterized by erythrocytosis (i.e., excessive red blood cell production) 43. Very interestingly, iron deficiency without anemia is present in virtually all patients of polycythemia vera 44. Our findings suggest that iron deficiency is the result, rather than the cause, of erythrocytosis and reaffirm the current practice of not using iron supplement to treat patients of polycythemia vera. At the other extreme of systemic status, iron overload can cause tissue damage and its excessive accumulation in liver leads to cirrhosis and hepatocellular carcinoma 3. Consistently, our study revealed that elevated serum iron in men is positively associated with risks of chronic liver disease and cirrhosis, and one of its subtypes, other chronic nonalcoholic liver disease. For copper, its shared protective effect with iron on iron deficiency anemia also has relatively well-established mechanistic basis. First, the efflux of iron from enterocyte, hepatocyte, and macrophage into the blood circulation requires the action of two copper-dependent ferroxidases (i.e., hephaestin and ceruloplasmin), which oxidize ferrous iron into ferric iron so that it could be loaded onto transferrin. Besides enabling iron transport, copper is also required for hemoglobin biosynthesis, probably by assisting iron import into or utilization in mitochondria 45.

Our study revealed the novel causal roles of elevated blood iron in increased risk of varicose vein. Varicose vein, a common venous disease of the lower extremity, is characterized by incompetent valves, reflux, and venous wall dilation. Its etiological process involves the hydrostatic-pressure-induced activation of matrix degrading enzymes and inflammatory cascade 46. Iron overload and its causal HFE genetic variations have been associated with varicose vein 47.

Mechanistically, iron overload induces oxidative stress and hyperactivation of matrix degrading enzymes 48. Other significant mineral-outcome relationships are also worth attention. For iron, our finding of its status increases the risk of glossitis (i.e., tongue inflammation) has been reported in a previous MR study 49. However, it is fairly well established that patients of atrophic glossitis frequently suffer from deficiencies of nutrients, including iron, and corresponding nutrient supplementation is able to resolve oral symptoms 50. How to reconcile these findings should be a subject of future studies. Our findings open many new venues of research to elucidate the roles of blood iron and copper in these clinical conditions.

Our study has strengths and limitations. This is the first comparative PheWAS-MR of iron and copper, based on a large prospective cohort, to simultaneously evaluate their shared and unique causal clinical effects. Our approach is unbiased and hypothesis-free. The MR approach reduces the bias from confounding and reserve causality. We utilized strong genetic instruments and applied multiple sensitivity analyses to detect the presence of pleiotropy. Only MR estimates with no indications of pleiotropy were reported. For selected outcomes (i.e., blood lipids), we replicated the results with an independent dataset. The sex-stratified analysis enabled the identification of sex-specific relationships. Our results confirmed some previous MR findings while made novel discoveries, most of which are supported by existing epidemiological and mechanistic studies. Our study also has multiple limitations. Despite a large sample size (N = 310,999) of the overall cohort, the case numbers of specific outcomes are still small, limiting our statistical power for rare diseases and those with modest effects. Moreover, since our genetic instruments only proxy blood mineral levels in the normal ranges, some clinical effects of extreme levels, such as deficiency and overload, are likely missed. Another related issue is that MR estimation assumes a linear relationship between the exposure and the outcome, and as a result, non-linear effects might have been missed. Given the well-known threshold effects of mineral deficiency and overload, future studies with a non-linear model will likely reveal more clinical effects. Genetic instruments approximate the average effects over the life course, while the physiological relevance of a blood mineral could vary by life stages and is not captured by our study. Finally, UK Biobank is a middle- and old-aged cohort and we restricted analysis to those of European descent, future studies in cohorts of younger ages or other ethnic backgrounds are needed to confirm our findings and to search for more clinical consequences.

## Conclusions

In conclusion, our comparative PheWAS-MR study of iron and copper comprehensively characterized their shared and unique clinical effects. These known and novel mineral-outcome relationships provide profound insights into the heath impacts of normally varying blood minerals and into the etiologies of some clinical conditions, specifically hyperlipidemia and hypercholesterolemia. Our findings also emphasize the possibility and importance of managing blood minerals, probably through dietary adjustment, in the prevention and management of these medical conditions.

## Data Availability

Genetic and phenotypic data from UK Biobank could be accessed following an application and approval process. Summary statistics generated in this project are provided in the supplementary files.

## Funding

AG is supported by United States Department of Agriculture, National Institute of Food and Agriculture Hatch Funds. MSR is supported by the Allen Foundation Inc., and the United States Department of Agriculture (MIN-18-118). KY is supported by the University of Georgia Research Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Conflict of interest

The authors have declared that no competing interests exist.

## Data availability statement

Genetic and phenotypic data from UK Biobank could be accessed following an application and approval process. Summary statistics generated in this project are provided in the supplementary files.

## Author Contributions

Conceptualization: Jingqi Zhou, Kaixiong Ye.

Data Curation: Jingqi Zhou, Chang Liu, Michael Francis.

Formal Analysis: Jingqi Zhou, Chang Liu.

Funding Acquisition: Moon-Suhn Ryu, Arthur Grider, Kaixiong Ye.

Investigation: Jingqi Zhou, Chang Liu, Arthur Grider, Kaixiong Ye.

Methodology: Jingqi Zhou, Kaixiong Ye.

Project Administration: Jingqi Zhou, Kaixiong Ye.

Resources: Moon-Suhn Ryu, Arthur Grider, Kaixiong Ye.

Software: Jingqi Zhou, Chang Liu, Michael Francis, Yitang Sun.

Supervision: Kaixiong Ye.

Validation: Michael Francis, Yitang Sun, Kaixiong Ye.

Visualization: Jingqi Zhou, Chang Liu.

Writing – Original Draft Preparation: Jingqi Zhou, Chang Liu.

Writing – Review & Editing: Jingqi Zhou, Chang Liu, Michael Francis, Yitang Sun, Moon-Suhn Ryu, Arthur Grider, Kaixiong Ye.

## Supporting Information

S1 Checklist: STROBE checklist. STROBE, Strengthening the Reporting of Observational Studies in Epidemiology.

S1 Fig: Genes in the KEGG pathway of cholesterol metabolism predicted to carry iron response elements. Predictions with high, medium, and low confidence are shown in yellow, blue, and orange, respectively.

S2 Fig: Genes in the KEGG pathway of fatty acid degradation predicted to carry iron response elements. Predictions with high, medium, and low confidence are shown in yellow, blue, and orange, respectively.

S1 Table: The number of phenotypes and cases considered in each disease category.

S2 Table: Descriptive characteristics of the UK Biobank participants and the SNPs included in PheWAS analyses.

S3 Table: The association between the weighted GRS of each blood mineral and common confounding factors.

S4 Table: Lists of significant mineral-outcome relationships without indications of pleiotropy. S5 Table: Multivarate IVW MR results for iron and copper.

S6 Table: Lists of significant mineral-outcome relationships without indications of pleiotropy based on sex-stratified MR analysis.

## Acknowledgments

This study was conducted using the UK Biobank Resource under project 48818. The authors would like to thank all UK Biobank participants. We also want to thank other members of the Ye lab for stimulating discussions.

## Abbreviations

CAD
:   coronary artery disease
FDR
:   false discovery rate
GLGC
:   Global Lipids Genetics Consortium
ICD
:   International Classification of Diseases
IVW
:   inverse variance-weighted
GRS
:   polygenic risk score
GWAS
:   genome-wide association study
HDL
:   high-density lipoprotein cholesterol
LDL
:   low-density lipoprotein cholesterol
MR
:   mendelian randomization
OR
:   odds ratio
OA
:   Osteoarthrosis
PheWAS
:   phenome-wide association study
SNP
:   single nucleotide polymorphism
SD
:   standard deviation
T2D
:   type 2 diabetes
TC
:   total cholesterol
WM
:   weighted median

*   Received September 17, 2020.
*   Revision received September 17, 2020.
*   Accepted September 22, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## Reference

1.  1.Hunsaker EW and Franz KJ. Emerging Opportunities To Manipulate Metal Trafficking for Therapeutic Benefit. 2019:13528–13545.
    
    

2.  2.Waldron KJ, Rutherford JC, Ford D and Robinson NJ. Metalloproteins and metal sensing. Nature. 2009;460:823–830.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature08300&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19675642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268938300026&link_type=ISI) 

3.  3.Ferreira CR and Gahl WA. Disorders of metal metabolism. Transl Sci Rare Dis. 2017;2:101–139.
    
    

4.  4.Camaschella C, Pagani A, Nai A and Silvestri L. The mutual control of iron and erythropoiesis. Int J Lab Hematol. 2016;38 Suppl 1:20–6.
    
    

5.  5.Pietrangelo A. Mechanisms of iron hepatotoxicity. J Hepatol. 2016;65:226–227.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jhep.2016.01.037&link_type=DOI) 

6.  6.Bandmann O, Weiss KH and Kaler SG. Wilson’s disease and other neurological copper disorders. Lancet Neurol. 2015;14:103–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(14)70190-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25496901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

7.  7.Soma S, Latimer AJ, Chun H, Vicary AC, Timbalia SA, Boulet A, et al. Elesclomol restores mitochondrial function in genetic models of copper deficiency. Proc Natl Acad Sci U S A. 2018;115:8161–8166.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE1LzMyLzgxNjEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8yMi8yMDIwLjA5LjE3LjIwMTk2OTU2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

8.  8.Bao W, Rong Y, Rong S and Liu L. Dietary iron intake, body iron stores, and the risk of type 2 diabetes: a systematic review and meta-analysis. BMC Med. 2012;10:119.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1741-7015-10-119&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23046549&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

9.  9.Kennish L, Attur M, Oh C, Krasnokutsky S, Samuels J, Greenberg JD, et al. Age-dependent ferritin elevations and HFE C282Y mutation as risk factors for symptomatic knee osteoarthritis in males: a longitudinal cohort study. BMC Musculoskelet Disord. 2014;15:8.
    
    

10. 10.Lima SC, Arrais RF, Sales CH, Almeida MG, de Sena KC, Oliveira VT, et al. Assessment of copper and lipid profile in obese children and adolescents. Biol Trace Elem Res. 2006;114:19–29.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17205984&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

11. 11.Yu L, Yan J, Zhang Q, Lin H, Zhu L, Liu Q, et al. Association between Serum Ferritin and Blood Lipids: Influence of Diabetes and hs-CRP Levels. J Diabetes Res. 2020;2020:4138696.
    
    

12. 12.Li J, Bao W, Zhang T, Zhou Y, Yang H, Jia H, et al. Independent relationship between serum ferritin levels and dyslipidemia in Chinese adults: A population study. PLoS One. 2017;12:e0190310.
    
    

13. 13.Davies NM, Holmes MV and Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE2OiIzNjIvanVsMTJfMi9rNjAxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMjIvMjAyMC4wOS4xNy4yMDE5Njk1Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

14. 14.Davey Smith G and Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23:R89–98.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddu328&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25064373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000349825700013&link_type=ISI) 

15. 15.Pichler I, Del Greco MF, Gogele M, Lill CM, Bertram L, Do CB, et al. Serum iron levels and the risk of Parkinson disease: a Mendelian randomization study. PLoS Med. 2013;10:e1001462.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001462&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23750121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

16. 16.Yuan S and Larsson S. Causal associations of iron status with gout and rheumatoid arthritis, but not with inflammatory bowel disease. Clin Nutr. 2020.
    
    

17. 17.Kodali HP, Pavilonis BT and Schooling CM. Effects of copper and zinc on ischemic heart disease and myocardial infarction: a Mendelian randomization study. Am J Clin Nutr. 2018;108:237–242.
    
    

18. 18.Cheng WW, Zhu Q and Zhang HY. Mineral Nutrition and the Risk of Chronic Diseases: A Mendelian Randomization Study. Nutrients. 2019;11.
    
    

19. 19.Benyamin B, Esko T, Ried JS, Radhakrishnan A, Vermeulen SH, Traglia M, et al. Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis. Nat Commun. 2014;5:4926.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms5926&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25352340&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

20. 20.Gill D, Monori G, Tzoulaki I and Dehghan A. Iron Status and Risk of Stroke. Stroke. 2018;49:2815–2821.
    
    

21. 21.Evans DM, Zhu G, Dy V, Heath AC, Madden PA, Kemp JP, et al. Genome-wide association study identifies loci affecting blood copper, selenium and zinc. Hum Mol Genet. 2013;22:3998–4006.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddt239&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23720494&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325163900015&link_type=ISI) 

22. 22.Wu P, Gifford A, Meng X, Li X, Campbell H, Varley T, et al. Mapping ICD-10 and ICD-10-CM Codes to Phecodes: Workflow Development and Initial Evaluation. JMIR Med Inform. 2019;7:e14325.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2196/14325&link_type=DOI) 

23. 23.Burgess S, Butterworth A and Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37:658–65.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21758&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24114802&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

24. 24.Bowden J, Davey Smith G, Haycock PC and Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40:304–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21965&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27061298&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

25. 25.Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7.
    
    

26. 26.Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, et al. Discovery and refinement of loci associated with lipid levels. Nature Genetics. 2013;45:1274–1283.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2797&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24097068&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

27. 27.Sanderson E, Davey Smith G, Windmeijer F and Bowden J. An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol. 2019;48:713–727.
    
    

28. 28.Burgess S and Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32:377–389.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.1007/s10654-017-0255-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28527048&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

29. 29.Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. 2019;35:4851–4853.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

30. 30.Yavorska OO and Burgess S. MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int J Epidemiol. 2017;46:1734–1739.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyx034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28398548&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

31. 31.Verbanck M, Chen CY, Neale B and Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50:693–698.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0099-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29686387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

32. 32.DiNicolantonio JJ, Mangan D and O’Keefe JH. Copper deficiency may be a leading cause of ischaemic heart disease. Open Heart. 2018;5:e000784.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6Nzoib3BlbmhydCI7czo1OiJyZXNpZCI7czoxMToiNS8yL2UwMDA3ODQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8yMi8yMDIwLjA5LjE3LjIwMTk2OTU2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

33. 33.Alarcon-Corredor OM, Guerrero Y, Ramirez de Fernandez M, D’Jesus I, Burguera M, Burguera JL, et al. Effect of copper supplementation on lipid profile of Venezuelan hyperlipemic patients. Arch Latinoam Nutr. 2004;54:413–8.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15969266&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

34. 34.Bo S, Durazzo M, Gambino R, Berutti C, Milanesio N, Caropreso A, et al. Associations of dietary and serum copper with inflammation, oxidative stress, and metabolic variables in adults. J Nutr. 2008;138:305–10.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjk6IjEzOC8yLzMwNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzIyLzIwMjAuMDkuMTcuMjAxOTY5NTYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

35. 35.Zang X, Huang H, Zhuang Z, Chen R, Xie Z, Xu C, et al. The association between serum copper concentrations and cardiovascular disease risk factors in children and adolescents in NHANES. Environ Sci Pollut Res Int. 2018;25:16951–16958.
    
    

36. 36.Zhu Y, He B, Xiao Y and Chen Y. Iron metabolism and its association with dyslipidemia risk in children and adolescents: a cross-sectional study. Lipids Health Dis. 2019;18:50.
    
    

37. 37. Das De S, Krishna S and Jethwa A. Iron status and its association with coronary heart disease: systematic review and meta-analysis of prospective studies. Atherosclerosis. 2015;238:296–303.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25544180&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

38. 38.Gill D, Brewer CF, Monori G, Tregouet DA, Franceschini N, Giambartolomei C, et al. Effects of Genetically Determined Iron Status on Risk of Venous Thromboembolism and Carotid Atherosclerotic Disease: A Mendelian Randomization Study. J Am Heart Assoc. 2019;8:e012994.
    
    

39. 39.Davis MR, Rendina E, Peterson SK, Lucas EA, Smith BJ and Clarke SL. Enhanced expression of lipogenic genes may contribute to hyperglycemia and alterations in plasma lipids in response to dietary iron deficiency. Genes Nutr. 2012;7:415–25.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12263-011-0278-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22228222&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

40. 40.Özdemir A, Sevinç C, Selamet U and Türkmen F. The Relationship Between Iron Deficiency Anemia and Lipid Metabolism in Premenopausal Women. The American Journal of the Medical Sciences. 2007;334:331–333.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/MAJ.0b013e318145b107&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18004086&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

41. 41.Wrede CE, Buettner R, Bollheimer LC, Scholmerich J, Palitzsch KD and Hellerbrand C. Association between serum ferritin and the insulin resistance syndrome in a representative population. Eur J Endocrinol. 2006;154:333–40.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZWplIjtzOjU6InJlc2lkIjtzOjk6IjE1NC8yLzMzMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzIyLzIwMjAuMDkuMTcuMjAxOTY5NTYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

42. 42.Campillos M, Cases I, Hentze MW and Sanchez M. SIREs: searching for iron-responsive elements. Nucleic Acids Res. 2010;38:W360–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkq371&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20460462&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

43. 43.Spivak JL. Myeloproliferative Neoplasms. N Engl J Med. 2017;377:895–6.
    
    

44. 44.Ginzburg YZ, Feola M, Zimran E, Varkonyi J, Ganz T and Hoffman R. Dysregulated iron metabolism in polycythemia vera: etiology and consequences. Leukemia. 2018;32:2105–2116.
    
    

45. 45.Doguer C, Ha JH and Collins JF. Intersection of Iron and Copper Metabolism in the Mammalian Intestine and Liver. Compr Physiol. 2018;8:1433–1461.
    
    

46. 46.Kucukguven A and Khalil RA. Matrix metalloproteinases as potential targets in the venous dilation associated with varicose veins. Curr Drug Targets. 2013;14:287–324.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23316963&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 

47. 47.Fukaya E, Flores AM, Lindholm D, Gustafsson S, Zanetti D, Ingelsson E, et al. Clinical and Genetic Determinants of Varicose Veins. Circulation. 2018;138:2869–2880.
    
    

48. 48.Zamboni P, Scapoli G, Lanzara V, Izzo M, Fortini P, Legnaro R, et al. Serum iron and matrix metalloproteinase-9 variations in limbs affected by chronic venous disease and venous leg ulcers. Dermatol Surg. 2005;31:644-9; discussion 649.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00042728-200506000-00005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15996413&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F22%2F2020.09.17.20196956.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229646400005&link_type=ISI) 

49. 49.Gill D, Benyamin B, Moore LSP, Monori G, Zhou A, Koskeridis F, et al. Associations of genetically determined iron status across the phenome: A mendelian randomization study. PLoS Med. 2019;16:e1002833.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002833&link_type=DOI) 

50. 50.Chiang CP, Chang JY, Wang YP, Wu YH, Wu YC and Sun A. Atrophic glossitis: Etiology, serum autoantibodies, anemia, hematinic deficiencies, hyperhomocysteinemia, and management. J Formos Med Assoc. 2020;119:774–780.