Skip to main content

Serum iron status and the risk of breast cancer in the European population: a two-sample Mendelian randomisation study

Abstract

Background

Previous observational studies have provided conflicting results on the association between serum iron status and the risk of breast cancer. Considering the relevance of this relationship to breast cancer prevention, its elucidation is warranted.

Object

We used a two-sample Mendelian randomisation (MR) study to explore the causal relationship between serum iron status and the risk of breast cancer.

Method

To select single nucleotide polymorphisms (SNPs) that could be used as instrumental variables for iron status, we used the Genetics of Iron Status consortium, which includes 11 discovery and 8 replication cohorts, encompassing 48,972 individuals of European descent. Moreover, we used the OncoArray network to select SNPs that could be considered instrumental variables for the outcome of interest (breast cancer); this dataset included 122,977 individuals of European descent with breast cancer and 105,974 peers without breast cancer. Both conservative (SNPs associated with overall iron status markers) and liberal (SNPs associated with the levels of at least one iron status marker) approaches were used as part of the MR analysis. For the former, we used an inverse-variance weighted (IVW) method, whereas for the latter, we used the IVW, MR-Egger regression, weighted median and simple mode methods.

Results

When the conservative approach was used, iron status showed no significant association with the risk of breast cancer or any of its subtypes. However, when the liberal approach was used, transferrin levels were found to be positively associated with the risk of ER-negative breast cancer based on the simple mode method (OR for MR, 1.225; 95% CI, 1.064, 1.410; P = 0.030). Nevertheless, the levels of the other iron status markers showed no association with the risk of breast cancer or its subtypes (P > 0.05).

Conclusion

In our MR study, the liberal approach suggested that changes in the concentration of transferrin could increase the risk of ER-negative breast cancer, although the levels of other iron status markers had no effect on the risk of breast cancer or its subtypes. This should be verified in future studies.

Introduction

The morbidity of breast cancer is increasing, affecting the quality of life of patients and their families and increasing their economic burden. Oxidative stress, a known cause of breast cancer, may have a role in this process [1]. Iron is a necessary micronutrient for the human body [2], and it plays an important role in various physiological processes, such as electron transfer, oxygen transport, immune function, DNA synthesis and energy production [3, 4]. Moreover, iron catalyses the generation of reactive oxygen species, and it could thus increase both oxidative stress and oncogene activation. Therefore, levels of iron in the body may affect the development of breast cancer [2, 5, 6].

Previous epidemiological studies have reported that higher levels of iron may be associated with a modestly increased risk of breast cancer [7, 8]. In contrast, studies by Quintana Pacheco et al. and Cade et al. have found that high levels of iron show an inverse relationship with breast cancer risk [9, 10]. Further, study by Kabat et al. has shown that iron status is not associated with the risk of breast cancer [11]. The link between iron status and breast cancer risk is therefore under much debate.

Previous studies have reported the presence of confounding factors, such as post-natal living environment, behaviour and habits, social status and environmental factors, that affect the association between iron status and breast cancer risk [12]. Indeed, such confounding factors could affect the causal inferences obtained in traditional epidemiological studies. Moreover, it is possible that breast cancer and iron levels have a reverse causal relationship, which could further affect the conclusions of traditional observational epidemiological studies.

Mendelian randomisation (MR) — a process in which genetic variations closely related to the exposure variable are used as instrumental variables — can help overcome the limitations of traditional epidemiological studies and allow us to make causal inferences regarding the effect of a particular exposure on an outcome [13]. Because alleles follow random distribution during gametogenesis, fertilised eggs have random genetic variants; thus, the genetic variations associated with the outcome or exposure are not affected by confounding factors or reverse causality [12].

To our knowledge, there have been no MR-based studies examining the relationship between iron status and breast cancer risk, even though the delineation of this relationship could have significant utility in breast cancer prevention and treatment. To this end, in the present study, we used MR to investigate whether iron status is related to the incidence of breast cancer using publicly available data from genome-wide association studies (GWASs).

Methods

We applied two-sample MR to summary data from the respective GWASs (Fig. 1). The original studies had been performed after obtaining informed consent from participants and had also received ethical approval. We used single nucleotide polymorphisms (SNPs) that showed strong relationships with total serum iron status as instrumental variables to explore the effect of iron status on breast cancer risk. By adjusting for the relationship of the instrumental variables (SNPs) with iron status and breast cancer, we estimated the effect of systemic iron status on breast cancer risk. To investigate whether iron status acts as a potential bias or mediator in the relationship between breast cancer risk and other breast cancer risk factors, we analysed the relationship between the SNPs associated with iron status and breast cancer risk factors.

Fig. 1
figure 1

Related databases and analysis methods in MR analysis. The publicly available summary data of SNP phenotypes were obtained from the largest meta-GWAS databases. The effect of iron status on asthma was estimated using a conservative approach (IVs: only SNPs connected with the concentrations of ferritin, serum iron and transferrin saturation were increased, concentrations of transferrin was decreased and systemic iron status would be increased) and a liberal approach (IVs: one of the SNPs was associated with breast cancer). In conservative instruments set, the inverse-variance-weighted (IVW) method was used, and in liberal instruments set, the IVW, MR Egger regression, weighted median and simple mode methods were used. MR, Mendelian randomization; IVs, instrumental variables; SNP, single nucleotide polymorphism; MR Egger, Mendelian randomization–Egger regression method; BMI, body mass index

Data sources: exposure

We obtained summary data from the largest meta-GWAS containing information on iron status, the Genetics of Iron Status (GIS) consortium, and performed MR analysis. This database included data on serum iron levels, transferrin levels, ferritin levels, and transferrin saturation. Data from 11 discovery and 8 replication cohorts encompassing 48,972 individuals of European ancestry were collected [14]. Benyamin et al. had performed genome-wide analyses within each cohort based on a uniform analysis plan after adjustment for factors such as principal component scores and age. Furthermore, the thresholds of the population stratification and quality control for each cohort were imputed into the score > 0.5, Hardy–Weinberg Equilibrium (HWE) ≥ 10-6, minor allele frequency (MAF) > 0.01 [14].

Data sources: outcome

Publicly available summary data on breast cancer were obtained from the largest meta-GWAS containing data on breast cancer risk — the OncoArray network — which contains data regarding five cancers: breast, ovarian, prostate, lung, and colorectal cancer [15]. The breast cancer database contained information on 122,977 individuals of European ancestry with breast cancer (cases) and 105,974 peers (controls). Of the case cohort, 69,501 individuals had ER-positive and 21,468 individuals had ER-negative breast cancer.

The thresholds of the population stratification and quality control for the GWAS analysis for the cohort were imputed into the score > 0.3, HWE ≥ 10−12, and MAF > 0.01. This information has been previously reported [16, 17].

Selection of instrumental variables

Instrumental variables that showed a strong association with iron status (P < 5 × 10−8) were selected from the GIS consortium (Tables 1, 2, 3, 4) [14]. All the SNPs could be found in the OncoArray network and showed linkage equilibrium (all pairwise r2 ≤ 0.01). None of the SNPs that were chosen as instrumental variables were correlated with breast cancer (P > 0.05). The corresponding statistical indicators (beta and SE) were obtained from the GIS consortium and OncoArray network databases. SNPs for iron status from the GIS consortium corresponded to the European population, with a sample size of 48,972 cases. SNPs for breast cancer from the OncoArray network database corresponded to the same European ethnic group, with a sample size of 122,977.

Table 1 Effect of iron concentration GWAS identified variants on breast cancer
Table 2 Effect of transferrin concentration GWAS identified variants on breast cancer
Table 3 Effect of ferritin concentration GWAS identified variants on breast cancer
Table 4 Effect of transferrin saturation concentration GWAS identified variants on breast cancer

To select instrumental variables, we used two analysis methods: conservative and liberal variable analyses. For conservative variable analyses, 3 SNPs (rs855791, rs1800562, rs1799945) strongly associated with the increasing concentrations of ferritin, serum iron and transferrin saturation and decreasing concentrations of transferrin (P < 5 × 10−8) were selected. Increased concentrations of ferritin and serum iron; increased transferrin saturation and decreased concentrations of transferrin would cause an improvement in systemic iron status [18]. Therefore, the genetic instrumental variables would have a coincident relationship with iron status via these four markers.

For liberal variable analyses, the SNPs strongly affiliated with at least one of the iron markers (5 SNPs for serum iron, 9 SNPs for transferrin, 5 SNPs for ferritin and 5 SNPs for transferrin saturation) in the GWAS (P < 5 × 10−8) were selected. Previous study has shown that there is a causal relationship between body mass index (BMI) and breast cancer [19]. Therefore, in order to rule out the possible confounding effect on the relationship between iron status and breast cancer, SNPs associated with BMI were retrieved from the GWAS database, GIANT, and we verified whether any of these were selected as instrumental variables in our study. In addition, we ensured that none of the instrumental variables were associated with the risk of breast cancer (all P > 0.05).

Validation of selected instrumental variables

To be validated as instrumental variables, the selected SNPs had to meet three important criteria [20]. First, they had to be strongly associated with the exposure (iron status). Second, they had to not be associated with any confounders of the relationship between the exposure (iron status) and outcome (breast cancer). Finally, they should have influenced the outcome (breast cancer) only via the exposure (iron status) (Fig. 2).

Fig. 2
figure 2

Diagram of the hypothesis of instrumental variables in Mendelian randomization study. IV, instrumental variable; BMI, body mass index

To minimise any possible weak instrumental variable bias, we ensured that the F statistic was above 10, indicating that the study had sufficient strength. We ensured that the F statistic for the instrumental variables was above 10, and this was used to impose restrictions on the possible bias [21]. To limit the possibility of bias related to population stratification, we ensured that both exposure and outcome cohorts included individuals of European descent. We performed three analyses to resolve issues associated with pleiotropy. First, we assessed the SNPs that were known asthma risk factors associated with breast cancer. Second, to analyse MR estimates, we used two approaches — the conservative approach (primary analysis) and the liberal approach (secondary analysis). Third, unknown directional pleiotropy was assessed using the MR-Egger test.

MR analyses

Two-sample MR was performed for testing the causal relationship between iron status and breast cancer. Moreover, breast cancer was subdivided as ER-positive and ER-negative breast cancer. Conservative and liberal methods were used in the MR analyses. As part of the conservative analysis, we used the inverse-variance-weighted (IVW) [21] method to conduct MR analyses. As part of the liberal analysis, we used the IVW, MR-Egger regression [22], weighted median, and simple mode methods [23] to estimate the effect of iron status on breast cancer risk. All of the data were selected from the OncoArray network and GIS consortium which were publicly available GWAS data (Fig. 1). R package ‘TwoSampleMR’ was used for MR analysis. Beta, SE and P values were obtained for the MR analysis. Beta values were transformed to odds ratios (ORs) using the following formula: OR = exp (beta). The 95% confidence interval (CI) was computed as follows: CI = exp (beta ± 1.96 × SE). Due to the limited consistency in publicly available data, the relationship between instrumental variables and other potential confounders, such as exercise and drinking, was difficult to assess. Hence, we used the GWAS catalogue database (https://www.ebi.ac.uk/gwas) to search for other phenotypes (including exercise, drinking, smoking, C-reactive protein levels, white blood cell counts) related to the selected instrumental variable SNPs and manually removed these SNPs from the MR analysis to rule out possible pleiotropic effects.

Sensitivity analyses

For the sensitivity analyses, we used the IVW and MR-Egger tests to evaluate heterogeneity and displayed the results using forest plots showing the value for each SNP and Cochran’s Q statistics [24, 25]. In addition, by deleting one SNP at a time and recomputing estimates for the overall instrument variable pool, a leave-one-out analysis was performed to identify SNPs that had a greater or non-proportional effect. To ensure that the MR analysis results were more robust, we also performed an MR-Egger statistical sensitivity analysis, which limited the pleiotropic effects of the instrumental variables. In MR-Egger regression, the intercept, as an indicator of the average pleiotropic deviation, is allowed to be freely estimated [26]. In the conservative approach, the MR-Egger method regression was not performed to test pleiotropic effects, because only 3 SNPs were used and one SNP was removed for LD, making the data insufficient for analysis [27]. For the same reason, the pleiotropic effects of iron and transferrin saturation on overall breast cancer risk could not be tested.

All above analyses were performed by R, version 3.6.1.

Results

Instrumental variables for iron status and breast cancer risk

Tables 1, 2, 3 and 4 show the connections of iron status with SNPs that were used as instrumental variables in the liberal and conservative analyses. For the liberal analyses, Tables 1, 2, 3 and 4 show the characteristics of the genetic variation related to iron concentration (3 SNPs for overall breast cancer and 5 SNPs for ER-positive breast cancer and ER-negative breast cancer), transferrin concentration (9 SNPs for ER-positive breast cancer and 8 SNPs for overall breast cancer and ER-negative breast cancer), ferritin (log) concentration (4 SNPs for ER-positive breast cancer and 5 SNPs for overall breast cancer and ER-negative breast cancer) and transferrin saturation concentration (3 SNPs for overall breast cancer and 5 SNPs for ER-positive breast cancer and ER-negative breast cancer), respectively. For the conservative analyses, we used 3 SNPs association with iron levels, transferrin levels, ferritin levels and transferrin saturation (rs1800562, rs1799945, and rs855791). F statistics for all the instrumental variables ranged from 40 (rs651007, ABO gene) to 3346 (rs8177240, TF gene), showing that all the SNPs were strong instrumental variables (Tables 1, 2, 3, 4).

The genetic instrument and breast cancer risk

None of the individual SNPs were associated with BMI, which is a confounding factor for breast cancer. In short, none of the selected variables associated with breast cancer risk in the liberal or conservative analyses were associated with increasing BMI (all P > 0.05) (Table 5).

Table 5 P value for the relationship between BMI and all genetic variations

Effect of iron status on breast cancer

Figure 3 shows the results of MR in examining the association of genetically predicted iron status with the risk of breast cancer. The ORs for the risk of breast cancer and its subtypes per SD increase in the levels of each iron marker are displayed. In the conservative analysis, none of the four iron status was found to be associated with overall breast cancer, ER-positive breast cancer and ER-negative breast cancer risk (all P > 0.05). In the liberal analysis, we found a positive correlation between transferrin levels and ER-negative breast cancer risk based on the simple mode (OR, 1.225; 95% CI, 1.064, 1.410; P = 0.030). However, the status of other iron markers had no association with the risk of breast cancer or its subtypes (P > 0.05).

Fig. 3
figure 3

MR estimates the association of genetically predicted iron status and the risk of breast cancer. The results showed ORs for breast cancer and their subtype per SD increase in every iron status biomarker. ad for iron, transferrin, ferritin, and transferrin saturation. MR, Mendelian randomization; MR-Egger, Mendelian randomization–Egger regression method; OR: odds ratio; 95% CI, 95% confidence interval; ER, oestrogen receptor

Sensitivity analyses

For MR estimates, we used the liberal instruments method, and the heterogeneity for overall breast cancer, ER-positive breast cancer, and ER-negative breast cancer did not show statistical significance (all P > 0.05). Using the IVW method, we found no evidence of heterogeneity for the associations between the status of the four iron markers (iron, transferrin, ferritin and transferrin saturation) and breast cancer risk (for overall breast cancer: Q 0.02, 7.59, 1.12 and 0.08; for ER-positive breast cancer: Q 4.99, 6.21, 2.07 and 3.16; for ER-negative breast cancer: Q 3.65, 7.68, 2.81 and 4.66, respectively; all P > 0.05). Moreover, when using the MR-Egger method in the liberal analysis, we did not identify aggregated directional pleiotropy between the levels of the four iron markers and breast cancer risk (ER-positive breast cancer: intercept −0.005, 0.003, 0.008 and 0.012; ER-negative breast cancer: intercept 0.010, 0.005, −0.001 and 0.017; overall breast cancer: intercept 0.0005 [transferrin] and −0.0004 [ferritin]; all P > 0.05) (Supplementary Figs. 1, 2 and 3).

Importantly, the MR estimates did not radically change when the leave-one-out analyses were performed, although the estimated direction showed differences (Supplementary Figs. 4, 5, 6 and 7). Previous GWASs have reported that rs174577 is also related with low-density lipoprotein cholesterol (LDL-D), high-density lipoprotein cholesterol (HDL-C), triglyceride (TG) and total cholesterol (TCHO) levels; rs4921915 is also related with TG and TCHO, and rs1800562 and rs651007 are related with LDL-C and TCHO [28]. Nevertheless, removing the 4 SNPs (rs174577, rs4921915, rs1800562 and rs651007) did not change the trends in the obtained results (Supplementary Figs. 4, 5, 6 and 7).

Discussion

Previous traditional observational studies have shown inconsistent results regarding the relationship between iron status and breast cancer [7,8,9,10,11], likely owing to differences in race, ethnicity and sample sizes. Moreover, reverse causality and residual confusion may also be present in these studies, and unconsidered confounding factors or unknown risk factors may influence the observed correlations between iron status and breast cancer. To overcome these limitations, in this study, we conducted two-sample MR to estimate the causal relationship between iron status and breast cancer risk using summary statistics data from the largest meta-GWASs of the European population. Our results showed that serum transferrin levels were positively associated with the risk of ER-positive breast cancer, but that other iron statuses had no association with the risk of breast cancer or its subtypes. To the best of our knowledge, the present study is the first to use MR to investigate the association between iron status and breast cancer.

In the causal estimation of the relationship between serum transferrin and ER-negative breast cancer risk, we used the simple mode approach in addition to the three traditional methods. A positive correlation [OR, 1.225 (1.084, 1.366); P = 0.030] was observed between transferrin levels and ER-negative breast cancer risk. It has been speculated that a high concentration of transferrin and transferrin receptors could increase the transport efficiency of iron, increasing intracellular iron concentration. This increase in iron concentration could lead to lipid peroxidation, gene mutations, DNA strand breakage and oncogene activation, thus leading to an increased risk of breast cancer [29, 30].

In addition, the simple mode approach is a new method wherein the mode-based estimate can be used to obtain a single causal estimation from multiple genetic instruments. It provides robustness to horizontal pleiotropy in a different manner than the IVW, MR-Egger and weighted median methods. The simple mode provides better detection capability than the MR-Egger estimate, although its detection capability is inferior to that of the IVW and weighted median methods [17]. In addition, the causal estimates of the simple mode for the relationship between transferrin and ER-negative breast cancer were marginally significant; hence, this relationship needs to be studied further.

As there is some risk of pleotropic effects in MR estimates [31], we searched for the possibility of this secondary effect using SNPs from the PubMed database. Our online search identified 4 SNPs — rs1800562 in HFE, rs174577 in FADS2, rs651007 in ABO and rs4921915 in NAT2 — related to LDL-C, TCHO, and/or TG levels, which have been reported to influence breast cancer risk [31, 32]. Nevertheless, the removal of these SNPs from our analysis did not cause any substantial changes in the MR estimation results, indicating that our results are unlikely to be biassed by blood lipid levels. In order to further test the robustness of our findings with regard to potential pleiotropic effects, we increased the number of SNPs available for analysis by relaxing the selection criteria for instrumental variables. In the sensitivity analysis, we found no significant differences via both the conservative and liberal approaches. The slight difference in estimation and confidence interval width between the different MR analysis methods may be accidental or a result of differences in measurement error, instead of reflecting actual differences. In addition, in the pleiotropic test, the MR-Egger method detected no bias. The public GWAS data on both exposure and outcome came from European cohorts, minimising the population bias. Further, our calculation results for the leave-one-out MR analysis were similar to the main MR estimates. Taken together, it appeared that there was no serious bias in the overall analysis and conclusion of our study.

Breast cancer is a heterogeneous disease with different histopathology and molecular subtypes, each with their own clinical prognosis and risk factors [33]. Previous studies have suggested that obesity may increase the risk of breast cancer [34]. Furthermore, BMI is related to iron concentrations and breast cancer risk [35, 36]. Data on BMI-related genetics were obtained from the GIANT consortium for 339,224 individuals of European descent [37]. None of the 11 SNPs we selected as instrumental variables were significantly associated with BMI (Table 5).

The change in iron status and the risk of breast cancer may be caused by common exposure factors. For example, inflammation affects iron status, increasing serum ferritin concentrations and reducing serum iron concentration [28, 38], and it could also increase oxidative tumour stress. Hence, inflammation may lead to an increase in iron status and the risk of breast cancer. However, due to the limited literature on this topic, further research is needed to fully explore these relationships.

There are some advantages to our study. We evaluated summary data from the largest meta-GWASs available, the GIS consortium and the OncoArray network. All the data extracted were from individuals of European descent, reducing any ancestry-related bias. In addition, we used two analysis methods to select instrumental variables — conservative and liberal analyses — which guarantees the robustness of our causal estimation.

However, there are some limitations to our present study. First, due to the limitations of the publicly available GWAS databases, it was difficult to perform hierarchical analysis according to factors such as age and sex in the combined exposures and outcome databases. Second, we used liberal instruments, which although provide more power to the study, also make the study particularly vulnerable to effects of pleiotropy. In this study, although we tried to reduce the effects of pleiotropy, bias due to the unknown biological functions of the SNPs concerned with iron status may be inevitable. Moreover, under ideal condition, for MR analysis, a sample size as large as possible should be used to make the results more reliable. Although the data for iron status and breast cancer were obtained from the GWAS databases with the largest sample size in the world, the sample size may still not have been fully ideal. The cohorts used in our study are all from European races, which contributing to reducing the bias brought by races. However, it is unknown whether the results are suitable for other races. Therefore, more researches should be conducted.

Conclusion

Our MR study indicates that changes in the serum transferrin concentration could increase the risk of ER-negative breast cancer, whereas the status of the other three iron statuses had no association with breast cancer. As the liberal instrument was relatively weak, these findings need to be verified in further studies.

Availability of data and materials

OncoArray network https://epi.grants.cancer.gov/gameon/

GWAS catalogue database https://www.ebi.ac.uk/gwas

GIANT consortium http://portals.broadinstitute.org/collaboration/giant/index.php/Main_Page

Abbreviations

MR:

Mendelian randomisation

SNPs:

Single nucleotide polymorphisms

IVW:

Inverse-variance weighted

GWAS:

Genome-wide association study

GIS:

Genetics of Iron Status

HWE:

Hardy–Weinberg Equilibrium Test

MAF:

Minor allele frequency

LDL-C:

LDL cholesterol

HDL-C:

HDL cholesterol

TG:

Triglyceride

TCHO:

Total cholesterol

MBE:

Mode-based estimate

References

  1. Gurer-Orhan H, Ince E, Konyar D, Saso L, Suzen S. The role of oxidative stress modulators in breast cancer. Curr Med Chem. 2018;25(33):4084–101. https://doi.org/10.2174/0929867324666170711114336.

    Article  CAS  PubMed  Google Scholar 

  2. Toyokuni S. Iron-induced carcinogenesis: the role of redox regulation. Free Radic Biol Med. 1996;20(4):553–66. https://doi.org/10.1016/0891-5849(95)02111-6.

    Article  CAS  PubMed  Google Scholar 

  3. Ponka P, Beaumont C, Richardson DR. Function and regulation of transferrin and ferritin. Semin Hematol. 1998;35(1):35–54.

    CAS  PubMed  Google Scholar 

  4. Rouault TA, Tong WH. Iron-sulfur cluster biogenesis and human disease. Trends Genet. 2008;24(8):398–407. https://doi.org/10.1016/j.tig.2008.05.008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Mccord JM. Iron, free radicals, and oxidative injury. Semin Hematol. 1998;35(1):5–12.

    CAS  PubMed  Google Scholar 

  6. Huang X. Iron overload and its association with cancer risk in humans: evidence for iron as a carcinogenic metal. Mutation Research-Fundamental and Molecular Mechanisms of Mutagenesis. 2003;533(1-2):153–71. https://doi.org/10.1016/j.mrfmmm.2003.08.023.

    Article  CAS  PubMed  Google Scholar 

  7. Cui Y, Vogt S, Olson N, Glass AG, Rohan TE. Levels of zinc, selenium, calcium, and iron in benign breast tissue and risk of subsequent breast cancer. Cancer Epidemiol Biomarkers Prev. 2007;16(8):1682–5. https://doi.org/10.1158/1055-9965.EPI-07-0187.

    Article  CAS  PubMed  Google Scholar 

  8. Gaur A, Collins H, Wulaningsih W, Holmberg L, Garmo H, Hammar N, et al. Iron metabolism and risk of cancer in the Swedish AMORIS study. Cancer Causes Control. 2013;24(7):1393–402. https://doi.org/10.1007/s10552-013-0219-8.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Quintana Pacheco DA, Sookthai D, Graf ME, Schübel R, Johnson T, Katzke VA, et al. Iron status in relation to cancer risk and mortality: findings from a population-based prospective study. Int J Cancer. 2018;143(3):561–9. https://doi.org/10.1002/ijc.31384.

    Article  CAS  PubMed  Google Scholar 

  10. Cade J, Thomas E, Vail A. Case-control study of breast cancer in south east England: nutritional factors. J Epidemiol Community Health. 1998;52(2):105–10. https://doi.org/10.1136/jech.52.2.105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Kabat GC, Miller AB, Jain M, Rohan TE. Dietary iron and heme iron intake and risk of breast cancer: a prospective cohort study. Cancer Epidemiol Biomarkers Prev. 2007;16(6):1306–8. https://doi.org/10.1158/1055-9965.EPI-07-0086.

    Article  CAS  PubMed  Google Scholar 

  12. Ebrahim S, Davey SG. Mendelian randomization: can genetic epidemiology help redress the failures of observational epidemiology? Hum Genet. 2008;123(1):15–33. https://doi.org/10.1007/s00439-007-0448-6.

    Article  PubMed  Google Scholar 

  13. Smith GD, Ebrahim S. “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22. https://doi.org/10.1093/ije/dyg070.

    Article  PubMed  Google Scholar 

  14. Benyamin B, Esko T, Ried JS, et al. Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis. Nat Commun. 2014;5:4926. https://doi.org/10.1038/ncomms5926.

    Article  CAS  PubMed  Google Scholar 

  15. Amos CI, Dennis J, Wang Z, Byun J, Schumacher FR, Gayther SA, et al. The OncoArray consortium: a network for understanding the genetic architecture of common cancers. Cancer Epidemiol Biomarkers Prev. 2017;26(1):126–35. https://doi.org/10.1158/1055-9965.EPI-16-0106.

    Article  PubMed  Google Scholar 

  16. Michailidou K, Lindstrom S, Dennis J, et al. Association analysis identifies 65 new breast cancer risk loci. Nature. 2017;551(7678):92-+. https://doi.org/10.1038/nature24284.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Milne RL, Kuchenbaecker KB, Michailidou K, et al. Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer. Nat Genet. 2017;49(12):1767–78. https://doi.org/10.1038/ng.3785.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Wish JB. Assessing iron status: beyond serum ferritin and transferrin saturation. Clin J Am Soc Nephrol. 2006;1(Suppl 1):S4–8. https://doi.org/10.2215/CJN.01490506.

    Article  CAS  PubMed  Google Scholar 

  19. Qian F, Wang S, Mitchell J, McGuffog L, Barrowdale D, Leslie G, et al. Height and body mass index as modifiers of breast cancer risk in BRCA1/2 mutation carriers: a Mendelian randomization study. J Natl Cancer Inst. 2019;111(4):350–64. https://doi.org/10.1093/jnci/djy132.

    Article  CAS  PubMed  Google Scholar 

  20. Zheng J, Baird D, Borges MC, Bowden J, Hemani G, Haycock P, et al. Recent developments in Mendelian randomization studies. Curr Epidemiol Rep. 2017;4(4):330–45. https://doi.org/10.1007/s40471-017-0128-6.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, et al. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21(3):223–42. https://doi.org/10.1177/0962280210394459.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Burgess S, Scott RA, Timpson NJ, Smith GD, Thompson SG. Consortium epic-interact. using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52. https://doi.org/10.1007/s10654-015-0011-z.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Hartwig FP, Davey SG, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98. https://doi.org/10.1093/ije/dyx102.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Del Greco MF, Minelli C, Sheehanc NA, Thompsonc JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Statistics in Medicine. 2015;34(21):2926–40. https://doi.org/10.1002/sim.6522.

    Article  Google Scholar 

  25. Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28(1):30–42. https://doi.org/10.1097/EDE.0000000000000559.

    Article  PubMed  Google Scholar 

  26. Bowden J, Smith GD, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25. https://doi.org/10.1093/ije/dyv080.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method (vol 32, pg 377, 2017). Eur J Epidemiol. 2017;32(5):391–2. https://doi.org/10.1007/s10654-017-0276-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45(11):1274–83. https://doi.org/10.1038/ng.2797.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Andrews NC, Schmidt PJ. Iron homeostasis. Annu Rev Physiol. 2007;69(1):69–85. https://doi.org/10.1146/annurev.physiol.69.031905.164337.

    Article  CAS  PubMed  Google Scholar 

  30. Kabat GC, Rohan TE. Does excess iron play a role in breast carcinogenesis? An unresolved hypothesis. Cancer Causes Control. 2007;18(10):1047–53. https://doi.org/10.1007/s10552-007-9058-9.

    Article  PubMed  Google Scholar 

  31. Holmes MV, Ala-Korpela M, Smith GD. Mendelian randomization in cardiometabolic disease: challenges in evaluating causality. Nat Rev Cardiol. 2017;14(10):577–90. https://doi.org/10.1038/nrcardio.2017.78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Beeghly-Fadiel A, Khankari NK, Delahanty RJ, Shu XO, Lu Y, Schmidt MK, et al. A Mendelian randomization analysis of circulating lipid traits and breast cancer risk. Int J Epidemiol. 2019;49(4):1117–31. https://doi.org/10.1093/ije/dyz242.

    Article  PubMed Central  Google Scholar 

  33. Weigelt B, Reis JS. Histological and molecular types of breast cancer: is there a unifying taxonomy? Nat Rev Clin Oncol. 2009;6(12):718–30. https://doi.org/10.1038/nrclinonc.2009.166.

    Article  CAS  PubMed  Google Scholar 

  34. Zimta AA, Tigu AB, Muntean M, Cenariu D, Slaby O, Berindan-Neagoe I. Molecular links between central obesity and breast cancer. Int J Mol Sci. 2019;20(21):5364. https://doi.org/10.3390/ijms20215364.

    Article  CAS  PubMed Central  Google Scholar 

  35. Zhao L, Zhang X, Shen Y, Fang X, Wang Y, Wang F. Obesity and iron deficiency: a quantitative meta-analysis. Obes Rev. 2015;16(12):1081–93. https://doi.org/10.1111/obr.12323.

    Article  CAS  PubMed  Google Scholar 

  36. Guo Y, Warren AS, Shu XO, et al. Genetically predicted body mass index and breast cancer risk: Mendelian randomization analyses of data from 145,000 women of European descent. PLoS Med. 2016;13(8):e1002105. https://doi.org/10.1371/journal.pmed.1002105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Locke AE, Kahali B, Berndt SI, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518(7538):197–206. https://doi.org/10.1038/nature14177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Butensky JE, Harmatz P, Lee M, et al. Altered iron metabolism in children with human immunodeficiency virus disease. Pediatr Hematol Oncol. 2009;26(2):69–84. https://doi.org/10.1080/08880010902754826.

    Article  Google Scholar 

Download references

Acknowledgements

We are extremely grateful to the GIS consortium, the OncoArray network, the GIANT consortium and all the participants in our study.

Funding

This study was funded by the Guangxi Degree and Postgraduate Education Reform Project (grant number JGY2020056), the Key Laboratory of Early Prevention and Treatment of Regional High-Incidence Tumours, Ministry of Education (grant number GKE2015-ZZ16).

Author information

Authors and Affiliations

Authors

Contributions

Yuang Mao and Yanling: designed the research and revised the manuscript; all authors: analysis and interpretation of data; Chenyang Hou and Qingzhi Hou: wrote the manuscript; Yuang Mao: responsible for data integrity and accuracy of data analysis and the authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yanling Hu or Yuang Mao.

Ethics declarations

Ethics approval and consent to participate

We use publicly available data and do not require ethical approval.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Fig. 1.

The SNP effects on iron status biomarkers and breast cancer for Scatterplot. iron (a), transferrin (b), ferritin (c), and transferrin saturation (d). MR Egger, Mendelian randomization–Egger regression method.

Additional file 2: Supplementary Fig. 2.

The SNP effects on iron status biomarkers and ER-positive breast cancer for Scatterplot. iron (a), transferrin (b), ferritin (c), and transferrin saturation (d). MR Egger, Mendelian randomization–Egger regression method, ER, estrogen receptor.

Additional file 3: Supplementary Fig. 3.

The SNP effects on iron status biomarkers and ER-negative breast cancer for Scatterplot. iron (a), transferrin (b), ferritin (c), and transferrin saturation (d). MR Egger, Mendelian randomization–Egger regression method; ER, estrogen receptor.

Additional file 4: Supplementary Fig. 4.

Odds ratio of breast cancer, ER-positive breast cancer, and ER-negative breast cancer risk per standard deviation increase in iron excluding 1 SNP at per time estimated by the inverse variance weighted and Wald ratio (only for breast cancer). SNP: single nucleotide polymorphism; OR: odds ratio; 95% CI, 95% confidence interval; ER, estrogen receptor.

Additional file 5: Supplementary Fig. 5.

Odds ratio of breast cancer, ER-positive breast cancer, and ER-negative breast cancer risk per standard deviation increase in transferrin excluding 1 SNP at per time estimated by the inverse variance weighted. SNP: single nucleotide polymorphism; OR: odds ratio; 95% CI, 95% confidence interval; ER, estrogen receptor.

Additional file 6: Supplementary Fig. 6.

Odds ratio of breast cancer, ER-positive breast cancer, and ER-negative breast cancer risk per standard deviation increase in ferritin excluding 1 SNP at per time estimated by the inverse variance weighted. SNP: single nucleotide polymorphism; OR: odds ratio; 95% CI, 95% confidence interval; ER, estrogen receptor.

Additional file 7: Supplementary Fig. 7.

Odds ratio of breast cancer, ER-positive breast cancer, and ER-negative breast cancer risk per standard deviation increase in Transferrin saturation excluding 1 SNP at per time estimated by the inverse variance weighted and Wald ratio (only for breast cancer). SNP: single nucleotide polymorphism; OR: odds ratio; 95% CI, 95% confidence interval; ER, estrogen receptor.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hou, C., Hou, Q., Xie, X. et al. Serum iron status and the risk of breast cancer in the European population: a two-sample Mendelian randomisation study. Genes Nutr 16, 9 (2021). https://doi.org/10.1186/s12263-021-00691-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12263-021-00691-7

Keywords