Key Points

Multitrait analysis of genome-wide association study boosts the statistical power to identify novel genetic traits for fibroblast growth factor 23.A functional genomics approach aided network discovery to identify histidine-rich glycoprotein (HRG) and high-mobility group protein box 1 (HMGB1) as key regulators of cardiac disease in CKD.Integration of clinical and genetic data enhances the discovery power and is crucial for understanding the genetic underpinnings of mineral bone disorder related to CKD.

Background

Genome-wide association studies (GWAS) have identified numerous genetic loci associated with mineral metabolism markers but have exclusively focused on single-trait analysis. In this study, we performed a multitrait analysis of GWAS (MTAG) of mineral metabolism, exploring overlapping genetic architecture between traits to identify novel genetic associations for fibroblast growth factor 23 (FGF23).

Methods

We applied MTAG to variants common to GWAS of five genetically correlated mineral metabolism markers in participants of European ancestry. We integrated UK Biobank GWAS for blood levels for phosphate, 25-hydroxyvitamin D, and calcium (n=366,484) and Cohorts for Heart and Aging Research in Genetic Epidemiology GWAS for parathyroid hormone (n=29,155) and FGF23 (n=13,716). We then used supervised and unsupervised deep machine learning to identify novel associations between genetic traits and FGF23.

Results

MTAG increased the effective sample size for mineral metabolism markers to n=50,325 for FGF23. After clumping, MTAG identified independent genome-wide significant single-nucleotide polymorphisms for all traits, including 62 loci for FGF23. Many of these loci have not been previously reported in single-trait analyses. Through a functional genomics approach, we identified histidine-rich glycoprotein (HRG) and high-mobility group box 1 (HMGB1) as master regulators of downstream canonical pathways associated with circulating FGF23, and both genes were highly enriched in hypertrophied cardiac tissue of deceased hemodialysis patients. In addition, we found that DNMT3A was associated with uremic toxin, 8-hydroxy-2-deoxyguanosine, a biomarker of DNA damage. In silico gene perturbation analysis revealed that DNMT3A is protective in patients with heart failure caused by hypertrophied or dilated cardiomyopathy.

Conclusions

Our findings highlight the importance of MTAG analysis of mineral metabolism markers to boost the number of genome-wide significant loci for FGF23 to identify novel genetic traits. Functional genomics revealed novel networks that inform unique cellular functions and identified HRG and HMGB1 as key master regulators of FGF23 and cardiovascular disease in CKD.