Genetically Predicted Fibroblast Growth Factor 23 and Major Cardiovascular Diseases, Their Risk Factors, Kidney Function, and Longevity: A Two-Sample Mendelian Randomization Study

Introduction: Fibroblast growth factor 23 (FGF23), a potential biomarker for kidney function, is related to cardiovascular disease (CVD) and diabetes, although it is unclear whether the relation is causal. This study evaluated the associations of genetically predicted FGF23 with major CVDs, their risk factors, kidney function, and longevity using Mendelian randomization (MR).

Methods: This is a two-sample MR study using summary statistics from large genome-wide association studies. Primary outcomes included coronary artery disease (CAD), myocardial infarction, heart failure, and atrial fibrillation. Secondary outcomes included cardiovascular risk factors, kidney function, and longevity. We used four single-nucleotide polymorphisms (SNPs) predicting FGF23, excluding rs2769071 in the ABO gene, which likely violates the MR exclusion-restriction assumption. We used inverse-variance weighted (IVW) as the primary statistical method to assess associations of FGF23 with the outcomes. Sensitivity analyses included weighted median (WM) and MR-Egger. We repeated the analyses including all five SNPs. Last, we validated the positive findings from the main analyses in a smaller study, i.e., FinnGen.

Results: Using IVW, genetically predicted higher FGF23 was inversely associated with risk of CAD [odds ratio (OR): 0.69 per logtransformed FGF23 (pg/ml) increase, 95% confidence interval (CI): 0.52–0.91] and type 2 diabetes mellitus (T2DM) (OR: 0.70, 95% CI: 0.52–0.96), but not with the other outcomes. The WM and MR-Egger estimates were directionally consistent.

Conclusion: This study suggests that genetically predicted higher FGF23 may be protective against CAD and T2DM. Future studies should explore the underlying mechanisms related to the potential protective effect of FGF23. FGF23 was unlikely a cause of poorer renal function.

Introduction

Cardiovascular disease (CVD) is the leading cause of mortality globally (World Health Organization, 2017). According to the Global Burden of Disease Study, an estimated 17.8 million people, around one-third of worldwide deaths, were from CVD in 2017 (Jagannathan et al., 2019). However, causes of CVD remain incompletely understood (Ezzati et al., 2015). Given that kidney diseases are often linked to CVD (Anavekar et al., 2004), recent studies have started to explore the role of kidney-function-related biomarkers, such as fibroblast growth factor 23 (FGF23) in CVD. FGF23 is mainly secreted by osteoblasts and osteocytes and is responsible for phosphate homeostasis (Kocełak et al., 2012). Previous observational studies showed that higher FGF23 was associated with higher risk of major CVDs (Batra et al., 2016), such as hypertension (Fyfe-Johsnon et al., 2016), coronary artery disease (CAD) (Lutsey et al., 2014), myocardial infarction (MI) (Di Giuseppe et al., 2015), atrial fibrillation (AF) (Mathew et al., 2014), and heart failure (HF) (Lutsey et al., 2014). However, a meta-analysis of FGF23 with risk of CVD in 17 general population cohorts suggested no causal relation (Marthi et al., 2018). These discrepancies may indicate the possibility of confounding, reverse causation, and selection bias given that most of these studies were observational.

Mendelian randomization (MR) is a potentially more credible design compared with conventional observational studies given the use of genetic variants randomly allocated at conception and, hence, is more resistant to confounding (Davies et al., 2018; Burgess et al., 2021). Previous MR studies have mainly focused on the associations of genetically predicted FGF23 with bone-related phenotypes, which showed that FGF23 is inversely related to bone mineral density and osteoporosis (Wang et al., 2020; Yokomoto-Umakoshi et al., 2020). Other MR studies also assessed its relation with CAD, stroke, blood pressure, and lipids (Yokomoto-Umakoshi et al., 2020; Zheng et al., 2020), but these studies have not assessed the relation with HF, AF, or type 2 diabetes mellitus (T2DM). Furthermore, since some of these previous studies were conducted in genome-wide association studies (GWAS), which have differences in disease and control definition, analysis model, and study design of the included studies, or may have included potentially invalid instrument such as a highly pleiotropic single-nucleotide polymorphism (SNP) in the ABO gene (Folkersen et al., 2020; Wang et al., 2020; Yokomoto-Umakoshi et al., 2020; Zheng et al., 2020). In view of the limited number of FGF23 instruments, such discrepancies can influence the overall assessment of causality for MR studies, as per indicated in previous MR studies such as growth differentiation factor 15 where there are inconsistent findings between CARDIoGRAM and UK Biobank (Au Yeung et al., 2019). To comprehensively evaluate the associations of genetically predicted FGF23 with CVD and its risk factors, we conducted a two-sample MR study using summary statistics from GWAS (Lawlor, 2016). Given that FGF23 is a potential biomarker for kidney function (Robinson-Cohen et al., 2018), we also explored its associations with kidney function for completeness. Finally, as previous studies suggested that FGF23 may be related to longevity from animal studies (Bär et al., 2018), we also assessed its association with longevity.

Materials and Methods

This is an MR study which is based on three assumptions. First, the genetic instruments, i.e., SNPs should predict FGF23. Second, the SNPs should not be associated with potential confounders. Last, the SNPs should only affect the outcome through affecting FGF23 (Figure 1).

FIGURE 1

Figure 1. Directed acyclic graph illustrating Mendelian randomization assumptions. Three assumptions should be satisfied: (1) the genetic variants, i.e., SNPs, should predict FGF23; (2) the SNPs should not be associated with potential confounders; (3) the SNPs should only affect the outcomes through affecting FGF23. SNPs, single nucleotide polymorphisms; CAD, coronary artery disease; MI, myocardial infarction; HF, heart failure; AF, atrial fibrillation; SBP, systolic blood pressure; DBP, diastolic blood pressure; FG, fasting glucose; HbA_1c, glycated hemoglobin; T2DM, type 2 diabetes mellitus; BMI, body mass index; eGFRcrea, estimated glomerular filtration rate from creatinine; eGFRcys, estimated glomerular filtration rate from cystatin C; UACR, urinary albumin-to-creatinine ratio; CKD, chronic kidney disease.

Data Sources

Exposure

Genetic predictors of FGF23 were obtained from a meta-analysis of GWAS consisting of seven studies with a total of 16,624 individuals of European ancestry (Robinson-Cohen et al., 2018). The mean age of the study populations ranged from 36.4 to 78.0 years old, and 45.5% (N = 7,572) were male. Circulating FGF23 was detected by enzyme-linked immunosorbent assay (ELISA; Kainos Laboratories Inc., Tokyo, Japan) for intact fibroblast growth factor 23 (iFGF23) in 82.5% participants (N = 13,716, five cohorts) and by ELISA kit (Immutopics, San Clemente, CA, United States) for C-terminal fragment fibroblast growth factor 23 (cFGF23) in 17.5% participants (N = 2,908, two cohorts). The GWAS excluded participants with estimated glomerular filtration rate based on creatinine (eGFRcrea) <30 ml/min per 1.73 m² based on the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation. FGF23 was log-transformed, and the GWAS was adjusted for age, sex, and the top 10 principal components of ancestry in linear regression (Robinson-Cohen et al., 2018). We considered that SNPs strongly associated with FGF23 (p-value < 5 × 10^–8) and excluded SNPs in high linkage disequilibrium (LD) (r² < 0.001). We also excluded the SNP (rs2769071) in the ABO gene given that it is known to be highly pleiotropic and, hence, may be an invalid instrument (Supplementary Figure 1; Li and Schooling, 2020), giving four SNPs in the main analysis, and all five in a supplementary analysis.

Outcomes

We extracted summary statistics for the outcomes from the relevant GWAS and the IEU GWAS database (Hemani et al., 2018). The primary outcomes included CAD (Nikpay et al., 2015), MI (Nikpay et al., 2015), HF (Shah et al., 2020), and AF (Roselli et al., 2018). We also included cardiovascular risk factors as secondary outcomes, including blood pressure [systolic blood pressure (SBP), diastolic blood pressure (DBP) (Mitchell et al., 2019)], body mass index (BMI) (Yengo et al., 2018), glycaemic traits [fasting glucose (FG) (Lagou et al., 2021), glycated hemoglobin (HbA_1c) (Wheeler et al., 2017)], and T2DM (restricted to European UK Biobank participants) (Mahajan et al., 2018), but not lipids and stroke because they have been thoroughly investigated in previous studies (Yokomoto-Umakoshi et al., 2020; Zheng et al., 2020). Given that FGF23 is closely related to kidney function (Robinson-Cohen et al., 2018), we also included kidney function as secondary outcomes. These traits included GFRcrea (Wuttke et al., 2019), eGFR based on serum cystatin C (eGFRcys) (Gorski et al., 2017), urinary albumin-to-creatinine ratio (UACR) (Teumer et al., 2019), and chronic kidney disease (CKD) (Wuttke et al., 2019). Last, given that the previous animal studies suggest that reduced FGF23 may decrease lifespan (Bär et al., 2018), we also included parental attained age (i.e., combination of mother and father’s age if alive or at death) in the UK Biobank (Pilling et al., 2017) and longevity (surviving to the 90th percentile age) (Deelen et al., 2019) as secondary outcomes. A proxy SNP in high LD (r² ≥ 0.8) was identified if the target SNP was not available for the outcome. LD proxies were obtained using 1,000-genome European phase 3 data. Palindromic SNPs were retained with minor allele frequency <0.3. Details of the data sources are in Table 1.

TABLE 1

Table 1. Information of outcomes included in the study.

Statistical Analyses

We aligned the SNPs on the same allele for exposure and outcome data using allele letter and effect allele frequency (Hartwig et al., 2016). We assessed instrument strength with FGF23 using the F-statistic where an F-statistic > 10 indicates that a weak instrument bias is unlikely (Bowden et al., 2016; Davies et al., 2018). We calculated the variance (R²) of FGF23 explained by the instruments using an equation used in previous MR studies (Yarmolinsky et al., 2018; Au Yeung et al., 2021) and used this information to compute the overall F-statistic (Yarmolinsky et al., 2018). As our primary statistical method, we assessed the role of FGF23 using inverse-variance weighted (IVW) with multiplicative random effects. IVW assumes all SNPs are valid or have balanced pleiotropy (Bowden and Holmes, 2019). We assessed heterogeneity of the Wald ratios (SNP on outcome divided by on exposure) using Cochran’s Q-test where high heterogeneity may indicate the presence of invalid instruments. Details of the statistical analysis methods are in Table 2.

TABLE 2

Table 2. Details of statistical analysis methods used in this Mendelian randomization study.

Sensitivity Analyses

The weighted median (WM) gives valid estimates as long as 50% of the weight is derived from valid SNPs (Bowden et al., 2016). MR-Egger allows for all SNPs to be invalid, if the instrument strength independent of direct effect (InSIDE) assumption is met, i.e., pleiotropic effects of the SNPs are not associated with the strength of SNP on exposure. The MR-Egger intercept (p-value < 0.05) indicates the presence of horizontal pleiotropy, i.e., SNPs affect the outcomes through genetic pathways independent of the exposure (Bowden et al., 2015). Directionally consistent findings from different methods may strengthen our findings (Lawlor et al., 2016). We also repeated the analyses including the SNP in the ABO gene to assess how the inclusion of this SNP may have affected the results. Last, we used FinnGen (January 14, 2020 release) with a sample size of up to 96,499 to validate any positive findings identified from the main analyses (Supplementary Figure 1).

Power Calculation

To correct for multiple testing, a Bonferroni corrected p-value of 0.0125 (i.e., 0.05/4) was considered as statistical significance for the primary outcomes. At the significance level of 0.0125 and 1.5% variance of FGF23 explained by four SNPs, we calculated statistical power for each primary outcome (Freeman et al., 2013; Burgess, 2014; Supplementary Figure 2).

All analyses were performed using R version 4.0.2 with the R packages (TwoSampleMR, version 0.5.5) (Hemani et al., 2018).

Ethics Approval

This study only used publicly available data, so no ethical approval is needed.

Results

We included four SNPs (i.e., rs17216707, rs11741640, rs17479566, and rs9925837), which explained 1.5% of the variance of FGF23 in the main analysis (Supplementary Figure 1). The overall F-statistic of the four SNPs was 63, and F-statistic of each SNP was larger than 10, indicating that a weak instrument bias is unlikely. Details of the instruments can be found in Supplementary Table 1 and Table 1.

Using IVW, higher FGF23 was associated with lower CAD risk [odds ratio (OR): 0.69 per natural log transformed FGF23 increase, 95% confidence interval (CI): 0.52–0.91], with directionally consistent findings from WM and MR-Egger. Similar findings were observed for MI. However, FGF23 was not associated with AF and HF although with wide 95% CIs. There was no evidence for heterogeneity based on the Cochran’s Q-test or MR-Egger intercept (Figure 2).

FIGURE 2

Figure 2. The associations of genetically predicted FGF23 and major cardiovascular diseases using Mendelian randomization. No. of SNPs, number of single nucleotide polymorphisms; IVW, inverse-variance weighted; WM, weighted median; CVD, cardiovascular disease; CAD, coronary artery disease; MI, myocardial infarction; HF, heart failure; AF, atrial fibrillation. ^∗Bonferroni-corrected P-value < 0.0125.

Figures 3A,B show the association of FGF23 with the secondary outcomes, including cardiovascular risk factors, kidney function, and longevity. We found an inverse association of FGF23 with T2DM risk (IVW OR: 0.70 per natural log transformed FGF23 increase, 95% CI: 0.52–0.96), with directionally consistent findings from WM and MR-Egger. However, FGF23 was not associated with CKD or with glycemic traits, blood pressure, BMI, kidney function, or longevity. Supplementary Table 2 gives associations of the SNPs with all outcomes.

FIGURE 3

Figure 3. The associations of genetically predicted FGF23 and cardiovascular risk factors, kidney function, and longevity using Mendelian randomization. (A) Type 2 diabetes mellitus, chronic kidney disease, and longevity; (B) Glycaemic traits, blood pressure traits, BMI, kidney function and longevity. No. of SNPs, number of single nucleotide polymorphisms; IVW, inverse-variance weighted; WM, weighted median; T2DM, type 2 diabetes mellitus; CKD, chronic kidney disease; FG, fasting glucose; HbA_1c, glycated hemoglobin; SBP, systolic blood pressure; DBP, diastolic blood pressure; BMI, body mass index; eGFRcrea, estimated glomerular filtration rate from creatinine; eGFRcys, estimated glomerular filtration rate from cystatin C; UACR, urinary albumin-to-creatinine ratio. ^∗P-value < 0.05.

We also assessed the associations of FGF23 with CAD and T2DM using a separate study (FinnGen study) as a verification (Supplementary Table 3), where we found directionally consistent findings using IVW and WM (Supplementary Table 4). However, the MR-Egger estimate was in the opposite direction with wide 95% CI.

In sensitivity analysis where we included rs2769071, a variant in the ABO gene, associations with CAD, MI, and T2DM were attenuated, while FGF23 remained unrelated to the other outcomes. However, the Cochran’s Q-statistics indicated heterogeneity for many analyses, suggesting including rs2769071 may invalidate the analyses, in particular, for IVW (Supplementary Figures 3, 4A,B and Supplementary Tables 2, 5).

Discussions

In this MR study, which included a comprehensive range of cardiovascular outcomes, their risk factors, kidney function, and longevity, inconsistent with previous observational studies, we found that genetically predicted FGF23 was inversely associated with CAD, MI, and T2DM but not with other cardiovascular outcomes or risk factors. We also found no strong evidence that FGF23 is related to kidney function. Last, we did not find evidence for an effect of genetically predicted FGF23 on longevity. Our study adds by clarifying, using genetic evidence, associations of FGF23 with CVD, T2DM, and kidney function, unbiased by confounding, as well as explaining the discrepant results in previous studies.

Previous observational studies found that FGF23 is positively associated with risk of CVD (Lutsey et al., 2014; Mathew et al., 2014; Di Giuseppe et al., 2015). However, these findings were not consistently observed in studies in animals and humans (Takashi et al., 2017; Liu et al., 2018; Pastor-Arroyo et al., 2018). Inconsistent findings may be due to confounding, for example, by socioeconomic position and lifestyle factors. Using a design, which is more resistant to confounding, we found an inverse association of FGF23 with CAD and MI. These findings are consistent with animal studies where lack of FGF23 increased the risk of age-related diseases (Bär et al., 2018).

The mechanisms underlying the effects of FGF23 remain unclear. FGF23 may decrease calcium reabsorption, increase phosphate excretion, and reduce vitamin D production (Liu and Quarles, 2007). Given that calcium is increasingly recognised as a cause of CAD (Xu et al., 2017), this may be one of the possible pathways by which FGF23 reduces the risk of CAD. Another possible pathway is coagulation, where the FGF23 increasing allele G of rs11741640 is strongly associated with prolonged activated partial thromboplastin time (p-value: 4.96 × 10^–93) in Biobank Japan (Hemani et al., 2018; Kanai et al., 2018). Alternatively, the protective effect on CAD could be mediated via reduced risk of T2DM (Ahmad et al., 2015). Regarding the inverse association of FGF23 in T2DM, although FGF23 may reduce vitamin D (Gutiérrez, 2010), a recent MR study showed no strong evidence for an association of FGF23 in vitamin D levels (Wang et al., 2020). Taking into account the potential protective effect of vitamin D on T2DM risk (Yuan et al., 2019), any potential protective effect of FGF23 on T2DM is unlikely via vitamin D pathways. The inconsistent associations with glycemic traits, such as HbA_1c and FG, also warrant future investigations. Better understanding of the role of FGF23 and closely related membrane-bound protein (e.g., klotho) (Razzaque, 2009) may also shed light on the mechanisms of medications known to have pleiotropic effects, such as metformin, whose pharmacological target, AMPK (Luo et al., 2020b), is a regulator of FGF23 production (Glosse et al., 2018).

Fibroblast growth factor 23 was not clearly associated with kidney function or CKD (Figures 2, 3A,B and Supplementary Figures 3, 4A,B), although associations have been seen in some (Larsson et al., 2003; Fliser et al., 2007) but not all studies (Levin et al., 2014; Alderson et al., 2016). These inconsistent findings may be due to selection bias as some of these studies were conducted in CKD patients. Alternatively, FGF23 could be a consequence of CKD instead of a cause (Fauconnier et al., 2019). The lack of association of FGF23 with kidney-related outcomes partly supports this argument. Whether CKD impacts FGF23 (i.e., reverse causation) can be explored using a bi-directional MR design in the future, when suitable studies are available such as genetic summary statistics of FGF23 becomes available.

Findings differed with or without inclusion of the SNP from the ABO gene (Figures 2, 3A,B and Supplementary Figures 3, 4A,B), which may help explain discrepancies across different MR studies (Folkersen et al., 2020). Other reasons accounting for discrepancies may include the choice of the outcome GWAS. For example, the inverse association of FGF23 with CAD was only observed in this MR but not a previous study, which included UK Biobank (Yokomoto-Umakoshi et al., 2020). Given that only a few SNPs were available for FGF23, the magnitude and direction of genetic associations can be influenced more easily due to differences in disease definition used in the respective GWAS, the inclusion of UK Biobank, and different analytic models, which requires further investigations.

Despite using a study design less susceptible to confounding than typical observational studies, there were some limitations. First, MR has stringent assumptions. We chose instruments for FGF23 from the most recent GWAS; previous studies suggested that FGF23 is predictive of bone mineral density, which is a known effect of FGF23 (Wang et al., 2020). Given that genetic variants were randomly allocated at conception, these instruments are unlikely to be confounded. It is more difficult to assess violation of the exclusion restriction assumption (i.e., instruments affect the outcomes other than via affecting the exposure) given the limited statistical power of sensitivity analyses with a small number of SNPs. This is particularly problematic for MR-Egger and could have explained the differences with other sensitivity analyses in the presence of outliers. Nevertheless, given that ABO gene is highly pleiotropic, we removed this SNP in the main analysis to reduce the likelihood of violation of the exclusion restriction assumption. Heterogeneity was higher when this SNP was included (Supplementary Figures 3, 4A,B). However, we could not rule out the possibility of horizontal pleiotropy in other SNPs. For example, rs11741640 was associated with some hematological traits (e.g., hemoglobin), and our previous MR studies suggest that hemoglobin may play a role in venous thromboembolism but not CAD (Zhong et al., 2016; Luo et al., 2020a). As such, replication of our findings using larger GWAS of FGF23 with more genetic instruments, including rare variants, which may have larger effect sizes, is warranted. Second, MR studies are also subject to selection bias, which may explain the null findings concerning HF and AF, which usually occurred at older ages (Schooling et al., 2021b), although FGF23 was not associated with longevity making this explanation less likely. Third, some of the studies in the FGF23 GWAS were also included in the outcome GWAS, although the proportion overlap was generally very small given the small sample size of the FGF23, apart from the eGFRcys GWAS (45%) (Supplementary Table 6). This implies the presence of weak instrument bias, the MR estimates would be biased toward null, although the F-statistics indicated low evidence for weak instrument bias (Burgess et al., 2016). Fourth, it is also increasingly recognised that covariable adjustment, either in the form of analysis or study design, may bias the MR estimates, which is one limitation regarding the use of summary statistics (Hartwig et al., 2021; Schooling et al., in press). Furthermore, it is possible that the effect of FGF23 on cardiovascular outcomes was sex specific, as evident in its relation with bone mineral density (Manolagas et al., 2013; Kodrič et al., 2019). However, we were unable to conduct sex-specific analyses given the lack of comprehensive summary statistics. Last, we were unable to conduct a bi-directional MR study to assess whether poorer kidney function increased FGF23 due to lack of relevant genetic summary statistics.

In conclusion, this MR study suggests that genetically predicted higher FGF23 may be protective against CAD and T2DM. Future studies should explore the mechanisms underpinning such relations, which may help explain the discrepant results concerning the role of FGF23 in CVD, T2DM, and kidney diseases across different studies, and further elucidate the role of FGF23 in these diseases.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

This study only used publicly available summary statistics from relevant genome-wide association studies (GWAS) and UK Biobank and hence no ethics approval was required. Respective ethics approval have been obtained by the GWAS and the UK Biobank investigators.

Author Contributions

SLAY and YL designed the study. YL wrote the analysis plan and interpreted the results, with feedback from SL and SLAY. YL wrote the first draft of the manuscript with critical feedback and revisions from SL, SLAY, and CMS. SLAY was the guarantor of the work. All authors contributed to the article and approved the submitted version.

Conflict of Interest

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

Publisher’s Note

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

Acknowledgments

Summary data on coronary artery disease and myocardial infarction have been contributed by CARDIoGRAMplusC4D investigators and have been downloaded from http://www.cardiogramplusc4d.org/. Summary data on heart failure have been contributed by The Heart Failure Molecular Epidemiology for Therapeutic Targets (HERMES) investigators and have been downloaded from https://www.hermesconsortium.org/. Summary data on atrial fibrillation have been contributed by 2018 AF HRC GWAS investigators and have been downloaded from Cardiovascular Disease Knowledge Portal https://cvd.hugeamp.org/. Summary data on fasting glucose and HbA_1c have been contributed by MAGIC investigators and have been downloaded from http://www.magicinvestigators.org/. Summary data on type 2 diabetes mellitus have been contributed by DIAMANTE investigators and have been downloaded from http://diagram-consortium.org/. Summary data on systolic blood pressure and diastolic blood pressure have been contributed by MRC-IEU investigators in which quality control filtering was conducted by R. Mitchell, G. Hemani, T. Dudding, L. Corbin, S. Harrison, and L. Paternoster as described in the published protocol (doi: 10.5523/bris.1ovaau5sxunp2cv8rcy88688v) and have been fetched from https://gwas.mrcieu.ac.uk/. Summary data on body mass index have been contributed by GIANT investigators and have been downloaded from https://portals.broadinstitute.org/collaboration/giant/index.php/GIANT_consortium_data_files. Summary data on creatinine-based estimation of GFR, cystatin C-based estimation of GFR, urinary albumin-to-creatinine ratio and chronic kidney disease have been contributed by CKDGen investigators and have been downloaded from https://ckdgen.imbi.uni-freiburg.de/#. Summary data on parental attained age have been contributed by Pilling et al. (2017), and have been downloaded from https://www.ebi.ac.uk/gwas/publications/29227965. Summary data on longevity (surviving to the 90th percentile age) have been contributed by Deelen et al. (2019), and have been downloaded from https://www.ebi.ac.uk/gwas/publications/31413261. We want to acknowledge the participants and investigators of the FinnGen study; summary data on coronary heart disease and type 2 diabetes mellitus from FinnGen study have been downloaded from https://www.finngen.fi/en. We thank Li Shun for the helpful discussions concerning the pleiotropic effects of ABO gene.

Supplementary Material

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

Supplementary Figure 1 | Study design of this Mendelian randomization study of genetically predicted FGF23 and cardiovascular diseases, their risk factors, kidney function, and longevity.

Supplementary Figure 2 | Power estimates for this Mendelian randomization study of FGF23 in major cardiovascular outcomes.

Supplementary Figure 3 | The associations of genetically predicted FGF23 and major cardiovascular diseases using Mendelian randomization, including the rs2769071 in ABO gene.

Supplementary Figure 4 | The associations of genetically predicted FGF23 on cardiovascular risk factors, kidney function, and longevity using Mendelian randomization, including the rs2769071 in ABO gene.

Supplementary Table 1 | Genetic predictors of log-transformed FGF23.

Supplementary Table 2 | Genetic associations of major CVDs, their risk factors, kidney function, and longevity.

Supplementary Table 3 | Genetic associations of CVDs and T2DM in FinnGen study.

Supplementary Table 4 | The associations of genetically predicted FGF23 and CAD and T2DM in FinnGen study using Mendelian randomization.

Supplementary Table 5 | Genetic association of rs2769071 in ABO gene with major CVD, their risk factors, kidney function, and longevity.

Supplementary Table 6 | Participant overlap between the FGF23 genome wide association studies (GWAS) and the outcome GWAS.

References

Ahmad, O. S., Morris, J. A., Mujammami, M., Forgetta, V., Leong, A., Li, R., et al. (2015). A Mendelian randomization study of the effect of type-2 diabetes on coronary heart disease. Nat. Commun. 6:7060. doi: 10.1038/ncomms8060

PubMed Abstract | CrossRef Full Text | Google Scholar

Alderson, H. V., Ritchie, J. P., Middleton, R., Larsson, A., Larsson, T. E., and Kalra, P. A. (2016). FGF-23 and Osteoprotegerin but not Fetuin-A are associated with death and enhance risk prediction in non-dialysis chronic kidney disease stages 3–5. Nephrology 21, 566–573. doi: 10.1111/nep.12664

PubMed Abstract | CrossRef Full Text | Google Scholar

Anavekar, N. S., Mcmurray, J. J., Velazquez, E. J., Solomon, S. D., Kober, L., Rouleau, J.-L., et al. (2004). Relation between renal dysfunction and cardiovascular outcomes after myocardial infarction. N. Engl. J. Med. 351, 1285–1295. doi: 10.1056/nejmoa041365

PubMed Abstract | CrossRef Full Text | Google Scholar

Au Yeung, S. L., Borges, M. C., Lawlor, D. A., and Schooling, C. M. (2021). The impact of lung function on cardiovascular diseases and cardiovascular risk factors: a two sample bi-directional Mendelian randomization study. Thorax doi: 10.1136/thoraxjnl-2020-215600 Epub ahead of print

CrossRef Full Text | PubMed Abstract | Google Scholar

Au Yeung, S. L., Luo, S., and Schooling, C. M. (2019). The impact of GDF-15, a biomarker for metformin, on the risk of coronary artery disease, breast and colorectal cancer, and type 2 diabetes and metabolic traits: a Mendelian randomisation study. Diabetologia 62, 1638–1646. doi: 10.1007/s00125-019-4913-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Bär, L., Feger, M., Fajol, A., Klotz, L.-O., Zeng, S., Lang, F., et al. (2018). Insulin suppresses the production of fibroblast growth factor 23 (FGF23). Proc. Natl. Acad. Sci. U.S.A. 115, 5804–5809.

Google Scholar

Batra, J., Buttar, R. S., Kaur, P., Kreimerman, J., and Melamed, M. L. (2016). FGF-23 and cardiovascular disease: review of literature. Curr. Opin. Endocrinol. Diabetes Obesity 23:423. doi: 10.1097/med.0000000000000294

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525. doi: 10.1093/ije/dyv080

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314. doi: 10.1002/gepi.21965

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowden, J., Del Greco, M. F., Minelli, C., Davey Smith, G., Sheehan, N. A., and Thompson, J. R. (2016). Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I 2 statistic. Int. J. Epidemiol. 45, 1961–1974.

Google Scholar

Bowden, J., and Holmes, M. V. (2019). Meta-analysis and Mendelian randomization: a review. Res. Synthesis Methods 10, 486–496. doi: 10.1002/jrsm.1346

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S. (2014). Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int. J. Epidemiol. 43, 922–929. doi: 10.1093/ije/dyu005

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., Davies, N. M., and Thompson, S. G. (2016). Bias due to participant overlap in two-sample Mendelian randomization. Genet. Epidemiol. 40, 597–608. doi: 10.1002/gepi.21998

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgess, S., O’donnell, C. J., and Gill, D. (2021). Expressing results from a mendelian randomization analysis: separating results from inferences. JAMA Cardiol. 6, 7–8.

Google Scholar

Davies, N. M., Holmes, M. V., and Smith, G. D. (2018). Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ 362:k601. doi: 10.1136/bmj.k601

PubMed Abstract | CrossRef Full Text | Google Scholar

Deelen, J., Evans, D. S., Arking, D. E., Tesi, N., Nygaard, M., Liu, X., et al. (2019). A meta-analysis of genome-wide association studies identifies multiple longevity genes. Nat. Commun. 10:3669.

Google Scholar

Di Giuseppe, R., Kühn, T., Hirche, F., Buijsse, B., Dierkes, J., Fritsche, A., et al. (2015). Plasma fibroblast growth factor 23 and risk of cardiovascular disease: results from the EPIC-Germany case-cohort study. Eur. J. Epidemiol. 30, 131–141. doi: 10.1007/s10654-014-9982-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ezzati, M., Obermeyer, Z., Tzoulaki, I., Mayosi, B. M., Elliott, P., and Leon, D. A. (2015). Contributions of risk factors and medical care to cardiovascular mortality trends. Nat. Rev. Cardiol. 12:508. doi: 10.1038/nrcardio.2015.82

PubMed Abstract | CrossRef Full Text | Google Scholar

Fauconnier, C., Roy, T., Gillerot, G., Roy, C., Pouleur, A.-C., and Gruson, D. (2019). FGF23: clinical usefulness and analytical evolution. Clin. Biochem. 66, 1–12. doi: 10.1016/j.clinbiochem.2019.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Fliser, D., Kollerits, B., Neyer, U., Ankerst, D. P., Lhotta, K., Lingenhel, A., et al. (2007). Fibroblast growth factor 23 (FGF23) predicts progression of chronic kidney disease: the Mild to Moderate Kidney Disease (MMKD) Study. J. Am. Soc. Nephrol. 18, 2600–2608. doi: 10.1681/asn.2006080936

PubMed Abstract | CrossRef Full Text | Google Scholar

Folkersen, L., Gustafsson, S., Wang, Q., Hansen, D. H., Hedman, ÅK., Schork, A., et al. (2020). Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals. Nat. Metab. 2, 1135–1148.

Google Scholar

Freeman, G., Cowling, B. J., and Schooling, C. M. (2013). Power and sample size calculations for Mendelian randomization studies using one genetic instrument. Int. J. Epidemiol. 42, 1157–1163. doi: 10.1093/ije/dyt110

PubMed Abstract | CrossRef Full Text | Google Scholar

Fyfe-Johsnon, A. L., Alonso, A., Selvin, E., Bower, J. K., Pankow, J. S., Agarwal, S. K., et al. (2016). Serum fibroblast growth factor-23 and incident hypertension: the atherosclerosis risk in communities study. J. Hypertension 34:1266. doi: 10.1097/hjh.0000000000000936

PubMed Abstract | CrossRef Full Text | Google Scholar

Gan, W., Clarke, R. J., Mahajan, A., Kulohoma, B., Kitajima, H., Robertson, N. R., et al. (2017). Bone mineral density and risk of type 2 diabetes and coronary heart disease: a Mendelian randomization study. Wellcome Open Res. 2:68. doi: 10.12688/wellcomeopenres.12288.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Glosse, P., Feger, M., Mutig, K., Chen, H., Hirche, F., Hasan, A. A., et al. (2018). AMP-activated kinase is a regulator of fibroblast growth factor 23 production. Kidney Int. 94, 491–501. doi: 10.1016/j.kint.2018.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorski, M., Van Der Most, P. J., Teumer, A., Chu, A. Y., Li, M., Mijatovic, V., et al. (2017). 1000 Genomes-based meta-analysis identifies 10 novel loci for kidney function. Sci. Rep. 7:45040.

Google Scholar

Gutiérrez, O. M. (2010). Fibroblast growth factor 23 and disordered vitamin d metabolism in chronic kidney disease: updating the “Trade-off” Hypothesis. Clin. J. Am. Soc. Nephrol. 5, 1710–1716. doi: 10.2215/cjn.02640310

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartwig, F. P., Davies, N. M., Hemani, G., and Davey Smith, G. (2016). Two-sample Mendelian Randomization: Avoiding the Downsides of A Powerful, Widely Applicable but Potentially Fallible Technique. Oxford: Oxford University Press.

Google Scholar

Hartwig, F. P., Tilling, K., Davey Smith, G., Lawlor, D. A., and Borges, M. C. (2021). Bias in two-sample Mendelian randomization when using heritable covariable-adjusted summary associations. Int. J. Epidemiol. doi: 10.1093/ije/dyaa266 Epub ahead of print

CrossRef Full Text | PubMed Abstract | Google Scholar

Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base platform supports systematic causal inference across the human phenome. Elife 7:e34408.

Google Scholar

Jagannathan, R., Patel, S. A., Ali, M. K., and Narayan, K. M. V. (2019). Global updates on cardiovascular disease mortality trends and attribution of traditional risk factors. Curr. Diabetes Rep. 19:44.

Google Scholar

Kanai, M., Akiyama, M., Takahashi, A., Matoba, N., Momozawa, Y., Ikeda, M., et al. (2018). Genetic analysis of quantitative traits in the Japanese population links cell types to complex human diseases. Nat. Genet. 50, 390–400. doi: 10.1038/s41588-018-0047-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Kocełak, P., Olszanecka-Glinianowicz, M., and Chudek, J. (2012). Fibroblast growth factor 23–structure, function and role in kidney diseases. Adv. Clin. Exp. Med. 21, 391–401.

Google Scholar

Kodrič, K., Zupan, J., Kranjc, T., Komadina, R., Mlakar, V., Marc, J., et al. (2019). Sex-determining region Y (SRY) attributes to gender differences in RANKL expression and incidence of osteoporosis. Exp. Mol. Med. 51, 1–16. doi: 10.1038/s12276-019-0294-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Lagou, V., Mägi, R., Hottenga, J.-J., Grallert, H., Perry, J. R., Bouatia-Naji, N., et al. (2021). Sex-dimorphic genetic effects and novel loci for fasting glucose and insulin variability. Nat. Commun. 12:24.

Google Scholar

Larsson, T., Nisbeth, U., Ljunggren, Ö, Jüppner, H., and Jonsson, K. B. (2003). Circulating concentration of FGF-23 increases as renal function declines in patients with chronic kidney disease, but does not change in response to variation in phosphate intake in healthy volunteers. Kidney Int. 64, 2272–2279. doi: 10.1046/j.1523-1755.2003.00328.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawlor, D. A. (2016). Commentary: two-sample Mendelian randomization: opportunities and challenges. Int. J. Epidemiol. 45:908. doi: 10.1093/ije/dyw127

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawlor, D. A., Tilling, K., and Davey Smith, G. (2016). Triangulation in aetiological epidemiology. Int. J. Epidemiol. 45, 1866–1886.

Google Scholar

Levin, A., Rigatto, C., Barrett, B., Madore, F., Muirhead, N., Holmes, D., et al. (2014). Biomarkers of inflammation, fibrosis, cardiac stretch and injury predict death but not renal replacement therapy at 1 year in a Canadian chronic kidney disease cohort. Nephrol. Dialysis Trans. 29, 1037–1047.

Google Scholar

Li, S., and Schooling, C. (2020). A phenome-wide association study of ABO blood groups. BMC Med. 18:334. doi: 10.1186/s12916-020-01795-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, E. S., Thoonen, R., Petit, E., Yu, B., Buys, E. S., Scherrer-Crosbie, M., et al. (2018). Increased circulating FGF23 does not lead to cardiac hypertrophy in the male Hyp mouse model of XLH. Endocrinology 159, 2165–2172. doi: 10.1210/en.2018-00174

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., and Quarles, L. D. (2007). How fibroblast growth factor 23 works. J. Am. Soc. Nephrol. 18, 1637–1647. doi: 10.1681/asn.2007010068

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, S., Au Yeung, S. L., Zuber, V., Burgess, S., and Schooling, C. M. (2020a). Impact of Genetically Predicted Red Blood Cell Traits on Venous Thromboembolism: Multivariable Mendelian Randomization Study Using UK Biobank. J. Am. Heart Assoc. 9:e016771.

Google Scholar

Luo, S., Schooling, C. M., Wong, I. C. K., and Yeung, S. L. A. (2020b). Evaluating the impact of AMPK activation, a target of metformin, on risk of cardiovascular diseases and cancer in the UK Biobank: a Mendelian randomisation study. Diabetologia 63, 2349–2358. doi: 10.1007/s00125-020-05243-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Lutsey, P. L., Alonso, A., Selvin, E., Pankow, J. S., Michos, E. D., Agarwal, S. K., et al. (2014). Fibroblast growth factor-23 and incident coronary heart disease, heart failure, and cardiovascular mortality: the Atherosclerosis Risk in Communities study. J. Am. Heart Assoc. 3:e000936.

Google Scholar

Mahajan, A., Wessel, J., Willems, S. M., Zhao, W., Robertson, N. R., Chu, A. Y., et al. (2018). Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes. Nat. Genet. 50, 559–571.

Google Scholar

Manolagas, S. C., O’brien, C. A., and Almeida, M. (2013). The role of estrogen and androgen receptors in bone health and disease. Nat. Rev. Endocrinol. 9:699. doi: 10.1038/nrendo.2013.179

PubMed Abstract | CrossRef Full Text | Google Scholar

Marthi, A., Donovan, K., Haynes, R., Wheeler, D. C., Baigent, C., Rooney, C. M., et al. (2018). Fibroblast growth factor-23 and risks of cardiovascular and noncardiovascular diseases: a meta-analysis. J. Am. Soc. Nephrol. 29, 2015–2027. doi: 10.1681/asn.2017121334

PubMed Abstract | CrossRef Full Text | Google Scholar

Mathew, J. S., Sachs, M. C., Katz, R., Patton, K. K., Heckbert, S. R., Hoofnagle, A. N., et al. (2014). Fibroblast growth factor-23 and incident atrial fibrillation: the Multi-Ethnic Study of Atherosclerosis (MESA) and the Cardiovascular Health Study (CHS). Circulation 130, 298–307. doi: 10.1161/circulationaha.113.005499

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, R., Elsworth, B., Mitchell, R., Raistrick, C., Paternoster, L., and Hemani, G. (2019). MRC IEU UK Biobank GWAS Pipeline Version 2. Bristol: University of Bristol.

Google Scholar

Nikpay, M., Goel, A., Won, H.-H., Hall, L. M., Willenborg, C., Kanoni, S., et al. (2015). A comprehensive 1000 Genomes–based genome-wide association meta-analysis of coronary artery disease. Nat. Genet. 47:1121. doi: 10.1038/ng.3396

PubMed Abstract | CrossRef Full Text | Google Scholar

Pastor-Arroyo, E.-M., Gehring, N., Krudewig, C., Costantino, S., Bettoni, C., Knöpfel, T., et al. (2018). The elevation of circulating fibroblast growth factor 23 without kidney disease does not increase cardiovascular disease risk. Kidney Int. 94, 49–59. doi: 10.1016/j.kint.2018.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Pilling, L. C., Kuo, C.-L., Sicinski, K., Tamosauskaite, J., Kuchel, G. A., Harries, L. W., et al. (2017). Human longevity: 25 genetic loci associated in 389,166 UK biobank participants. Aging (Albany NY) 9:2504. doi: 10.18632/aging.101334

PubMed Abstract | CrossRef Full Text | Google Scholar

Razzaque, M. S. (2009). The FGF23-Klotho axis: endocrine regulation of phosphate homeostasis. Nat. Rev. Endocrinol. 5, 611–619. doi: 10.1038/nrendo.2009.196

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson-Cohen, C., Bartz, T. M., Lai, D., Ikizler, T. A., Peacock, M., Imel, E. A., et al. (2018). Genetic variants associated with circulating fibroblast growth factor 23. J. Am. Soc. Nephrol. 29, 2583–2592.

Google Scholar

Roselli, C., Chaffin, M. D., Weng, L.-C., Aeschbacher, S., Ahlberg, G., Albert, C. M., et al. (2018). Multi-ethnic genome-wide association study for atrial fibrillation. Nat. Genet. 50, 1225–1233.

Google Scholar

Schooling, C., Zhao, J., Au Yeung, S., and Kwok, M. (in press). Letter in response to “Bias in two-sample Mendelian randomization when using heritable covariable-adjusted summary associations” – “Interpreting Mendelian randomization studies pre-adjusted for the heritable covariable survival to recruitment. Int. J. Epidemiol.

Google Scholar

Schooling, C. M., Lopez, P. M., Yang, Z., Zhao, J. V., Au Yeung, S. L., and Huang, J. V. (2021b). Use of multivariable mendelian randomization to address biases due to competing risk before recruitment. Front. Genet. 11:852. doi: 10.3389/fgene.2020.610852

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, S., Henry, A., Roselli, C., Lin, H., Sveinbjörnsson, G., Fatemifar, G., et al. (2020). Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure. Nat. Commun. 11:163.

Google Scholar

Takashi, Y., Kinoshita, Y., Hori, M., Ito, N., Taguchi, M., and Fukumoto, S. (2017). Patients with FGF23-related hypophosphatemic rickets/osteomalacia do not present with left ventricular hypertrophy. Endocrine Res. 42, 132–137. doi: 10.1080/07435800.2016.1242604

PubMed Abstract | CrossRef Full Text | Google Scholar

Teumer, A., Li, Y., Ghasemi, S., Prins, B. P., Wuttke, M., Hermle, T., et al. (2019). Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria. Nat. Commun. 10:4130.

Google Scholar

Wang, Y., Wang, H., and Chen, P. (2020). higher fibroblast growth factor 23 levels are causally associated with lower bone mineral density of heel and femoral neck: evidence from two-sample mendelian randomization analysis. Front. Public Health 8:467. doi: 10.3389/fpubh.2020.00467

PubMed Abstract | CrossRef Full Text | Google Scholar

Wheeler, E., Leong, A., Liu, C.-T., Hivert, M.-F., Strawbridge, R. J., Podmore, C., et al. (2017). Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: a transethnic genome-wide meta-analysis. PLoS Med. 14:e1002383. doi: 10.1371/journal.pmed.1002383

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization (2017). Cardiovascular Diseases (CVDs) [Online]. Available online at: https://www.who.int/en/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds) (Accessed 05 April, 2020)

Google Scholar

Wuttke, M., Li, Y., Li, M., Sieber, K. B., Feitosa, M. F., Gorski, M., et al. (2019). A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat. Genet. 51:957.

Google Scholar

Xu, L., Lin, S. L., and Schooling, C. M. (2017). A Mendelian randomization study of the effect of calcium on coronary artery disease, myocardial infarction and their risk factors. Sci. Rep. 7:42691.

Google Scholar

Yarmolinsky, J., Bonilla, C., Haycock, P. C., Langdon, R. J., Lotta, L. A., Langenberg, C., et al. (2018). Circulating selenium and prostate cancer risk: a Mendelian randomization analysis. JNCI J. Natl. Cancer Institute 110, 1035–1038. doi: 10.1093/jnci/djy081

PubMed Abstract | CrossRef Full Text | Google Scholar

Yengo, L., Sidorenko, J., Kemper, K. E., Zheng, Z., Wood, A. R., Weedon, M. N., et al. (2018). Meta-analysis of genome-wide association studies for height and body mass index in∼ 700000 individuals of European ancestry. Hum. Mol. Genet. 27, 3641–3649. doi: 10.1093/hmg/ddy271

PubMed Abstract | CrossRef Full Text | Google Scholar

Yokomoto-Umakoshi, M., Umakoshi, H., Miyazawa, T., Ogata, M., Sakamoto, R., and Ogawa, Y. (2020). Investigating the causal effect of fibroblast growth factor 23 on osteoporosis and cardiometabolic disorders: a Mendelian randomization study. Bone 143:115777. doi: 10.1016/j.bone.2020.115777

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, S., Jiang, X., Michaëlsson, K., and Larsson, S. C. (2019). Genetic Prediction of Serum 25-Hydroxyvitamin D, calcium, and parathyroid hormone levels in relation to development of type 2 diabetes: a mendelian randomization study. Diabetes Care 42, 2197–2203. doi: 10.2337/dc19-1247

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, K., Lin, L., Cui, P., Liu, T., Chen, L., Yang, C., et al. (2020). Association of fibroblast growth factor 23 with ischemic stroke and its subtypes: a mendelian randomization study. Front. Genet. 11:608517. doi: 10.3389/fgene.2020.608517

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, Y., Lin, S. L., and Schooling, C. M. (2016). The effect of hematocrit and hemoglobin on the risk of ischemic heart disease: a Mendelian randomization study. Prev. Med. 91, 351–355. doi: 10.1016/j.ypmed.2016.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar