GWAS meta-analysis followed by Mendelian randomization revealed potential control mechanisms for circulating α-Klotho levels

Abstract The protein α-Klotho acts as transmembrane co-receptor for fibroblast growth factor 23 (FGF23) and is a key regulator of phosphate homeostasis. However, α-Klotho also exists in a circulating form, with pleiotropic, but incompletely understood functions and regulation. Therefore, we undertook a genome-wide association study (GWAS) meta-analysis followed by Mendelian randomization (MR) of circulating α-Klotho levels. Plasma α-Klotho levels were measured by enzyme-linked immunosorbent assay (ELISA) in the Ludwigshafen Risk and Cardiovascular Health and Avon Longitudinal Study of Parents and Children (mothers) cohorts, followed by a GWAS meta-analysis in 4376 individuals across the two cohorts. Six signals at five loci were associated with circulating α-Klotho levels at genome-wide significance (P < 5 × 10−8), namely ABO, KL, FGFR1, and two post-translational modification genes, B4GALNT3 and CHST9. Together, these loci explained >9% of the variation in circulating α-Klotho levels. MR analyses revealed no causal relationships between α-Klotho and renal function, FGF23-dependent factors such as vitamin D and phosphate levels, or bone mineral density. The screening for genetic correlations with other phenotypes followed by targeted MR suggested causal effects of liability of Crohn’s disease risk [Inverse variance weighted (IVW) beta = 0.059 (95% confidence interval 0.026, 0.093)] and low-density lipoprotein cholesterol levels [−0.198 (−0.332, −0.063)] on α-Klotho. Our GWAS findings suggest that two enzymes involved in post-translational modification, B4GALNT3 and CHST9, contribute to genetic influences on α-Klotho levels, presumably by affecting protein turnover and stability. Subsequent evidence from MR analyses on α-Klotho levels suggest regulation by mechanisms besides phosphate-homeostasis and raise the possibility of cross-talk with FGF19- and FGF21-dependent pathways, respectively. Significance statement: α-Klotho as a transmembrane protein is well investigated along the endocrine FGF23-α-Klotho pathway. However, the role of the circulating form of α-Klotho, which is generated by cleavage of transmembrane α-Klotho, remains incompletely understood. Genetic analyses might help to elucidate novel regulatory and functional mechanisms. The identification of genetic factors related to circulating α-Klotho further enables MR to examine causal relationships with other factors. The findings from the first GWAS meta-analysis of circulating α-Klotho levels identified six genome-wide significant signals across five genes. Given the function of two of the genes identified, B4GALNT3 and CHST9, it is tempting to speculate that post-translational modification significantly contributes to genetic influences on α-Klotho levels, presumably by affecting protein turnover and stability.

[1]  J. Erdmann,et al.  Genomewide Association Study of Severe Covid-19 with Respiratory Failure , 2020, The New England journal of medicine.

[2]  F. Martinez,et al.  Severe Covid-19. , 2020, The New England journal of medicine.

[3]  Tom M Palmer,et al.  Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: A multivariable Mendelian randomisation analysis , 2020, PLoS medicine.

[4]  James R. Apgar,et al.  Structure-function relationships of the soluble form of the antiaging protein Klotho have therapeutic implications for managing kidney disease , 2020, The Journal of Biological Chemistry.

[5]  J. Navarro-González,et al.  Association between serum levels of Klotho and inflammatory cytokines in cardiovascular disease: a case-control study , 2020, Aging.

[6]  Xiaokun Li The FGF metabolic axis , 2019, Frontiers of Medicine.

[7]  David A. Drew,et al.  Performance of soluble Klotho assays in clinical samples of kidney disease , 2019, Clinical kidney journal.

[8]  Mark R. Hurle,et al.  Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases , 2019, bioRxiv.

[9]  Weiping Xia,et al.  Klotho gene polymorphisms are associated with healthy aging and longevity: Evidence from a meta-analysis , 2019, Mechanisms of Ageing and Development.

[10]  E. Zeggini,et al.  Mendelian Randomization Analysis Reveals a Causal Influence of Circulating Sclerostin Levels on Bone Mineral Density and Fractures , 2019, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[11]  Di Zou,et al.  The role of klotho in chronic kidney disease , 2018, BMC Nephrology.

[12]  M. Kuro-o The Klotho proteins in health and disease , 2018, Nature Reviews Nephrology.

[13]  Sina A. Gharib,et al.  Unraveling the polygenic architecture of complex traits using blood eQTL metaanalysis , 2018, bioRxiv.

[14]  Rong Wang,et al.  Correlation between Soluble α-Klotho and Renal Function in Patients with Chronic Kidney Disease: A Review and Meta-Analysis , 2018, BioMed research international.

[15]  M. Gawaz,et al.  Therapeutic Interference With Vascular Calcification—Lessons From Klotho-Hypomorphic Mice and Beyond , 2018, Front. Endocrinol..

[16]  F. Windmeijer,et al.  An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings , 2018, bioRxiv.

[17]  J. Xie,et al.  New Insights into the Mechanism of Action of Soluble Klotho , 2017, Front. Endocrinol..

[18]  J. Michael Cherry,et al.  The Encyclopedia of DNA elements (ENCODE): data portal update , 2017, Nucleic Acids Res..

[19]  Michael H. Guo,et al.  Epigenetic profiling of growth plate chondrocytes sheds insight into regulatory genetic variation influencing height , 2017, eLife.

[20]  Nicola J. Rinaldi,et al.  Genetic effects on gene expression across human tissues , 2017, Nature.

[21]  Tom R. Gaunt,et al.  PhenoSpD: an integrated toolkit for phenotypic correlation estimation and multiple testing correction using GWAS summary statistics , 2017, bioRxiv.

[22]  Fernando Pires Hartwig,et al.  Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption , 2017, bioRxiv.

[23]  N. Sheehan,et al.  Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic , 2016, International journal of epidemiology.

[24]  Christian Gieger,et al.  Novel multiple sclerosis susceptibility loci implicated in epigenetic regulation , 2016, Science Advances.

[25]  Tom R. Gaunt,et al.  LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis , 2016, bioRxiv.

[26]  David M. Evans,et al.  Mendelian Randomization: New Applications in the Coming Age of Hypothesis-Free Causality. , 2015, Annual review of genomics and human genetics.

[27]  M. Daly,et al.  An Atlas of Genetic Correlations across Human Diseases and Traits , 2015, Nature Genetics.

[28]  G. Davey Smith,et al.  Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression , 2015, International journal of epidemiology.

[29]  M. Mohammadi,et al.  The demonstration of αKlotho deficiency in human chronic kidney disease with a novel synthetic antibody. , 2015, Nephrology, dialysis, transplantation : official publication of the European Dialysis and Transplant Association - European Renal Association.

[30]  G. Davey Smith,et al.  Mendelian randomization: genetic anchors for causal inference in epidemiological studies , 2014, Human molecular genetics.

[31]  M. Daly,et al.  LD Score regression distinguishes confounding from polygenicity in genome-wide association studies , 2014, Nature Genetics.

[32]  A. Ortiz,et al.  Hyperlipidemia-Associated Renal Damage Decreases Klotho Expression in Kidneys from ApoE Knockout Mice , 2013, PloS one.

[33]  C. Wallace,et al.  Bayesian Test for Colocalisation between Pairs of Genetic Association Studies Using Summary Statistics , 2013, PLoS genetics.

[34]  O. Delaneau,et al.  Supplementary Information for ‘ Improved whole chromosome phasing for disease and population genetic studies ’ , 2012 .

[35]  Eurie L. Hong,et al.  Annotation of functional variation in personal genomes using RegulomeDB , 2012, Genome research.

[36]  David Haussler,et al.  The UCSC genome browser and associated tools , 2012, Briefings Bioinform..

[37]  J. Marchini,et al.  Fast and accurate genotype imputation in genome-wide association studies through pre-phasing , 2012, Nature Genetics.

[38]  D. Fiete,et al.  Peptide-specific Transfer of N-Acetylgalactosamine to O-Linked Glycans by the Glycosyltransferases β1,4-N-Acetylgalactosaminyl Transferase 3 (β4GalNAc-T3) and β4GalNAc-T4* , 2012, The Journal of Biological Chemistry.

[39]  D. Lawlor,et al.  Cohort Profile: The Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort , 2012, International journal of epidemiology.

[40]  D. Lawlor,et al.  Cohort Profile: The ‘Children of the 90s’—the index offspring of the Avon Longitudinal Study of Parents and Children , 2012, International journal of epidemiology.

[41]  P. Visscher,et al.  Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits , 2012, Nature Genetics.

[42]  K. Rosenblatt,et al.  Klotho Depletion Contributes to Increased Inflammation in Kidney of the db/db Mouse Model of Diabetes via RelA (Serine)536 Phosphorylation , 2011, Diabetes.

[43]  M. Econs,et al.  Establishment of sandwich ELISA for soluble alpha-Klotho measurement: Age-dependent change of soluble alpha-Klotho levels in healthy subjects. , 2010, Biochemical and biophysical research communications.

[44]  Yun Li,et al.  METAL: fast and efficient meta-analysis of genomewide association scans , 2010, Bioinform..

[45]  D. Besselsen,et al.  Tumor necrosis factor and interferon-gamma down-regulate Klotho in mice with colitis. , 2010, Gastroenterology.

[46]  Aaron R. Quinlan,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2022 .

[47]  P. Donnelly,et al.  A Flexible and Accurate Genotype Imputation Method for the Next Generation of Genome-Wide Association Studies , 2009, PLoS genetics.

[48]  S. Leeman,et al.  Insulin stimulates the cleavage and release of the extracellular domain of Klotho by ADAM10 and ADAM17 , 2007, Proceedings of the National Academy of Sciences.

[49]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[50]  Daniel L. Koller,et al.  The role of mutant UDP-N-acetyl-alpha-D-galactosamine-polypeptide N-acetylgalactosaminyltransferase 3 in regulating serum intact fibroblast growth factor 23 and matrix extracellular phosphoglycoprotein in heritable tumoral calcinosis. , 2006, The Journal of clinical endocrinology and metabolism.

[51]  S. Ebrahim,et al.  'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? , 2003, International journal of epidemiology.

[52]  Monica M. Palcic,et al.  The structural basis for specificity in human ABO(H) blood group biosynthesis , 2002, Nature Structural Biology.

[53]  Daniel S. Roseman,et al.  The Mannose/N-Acetylgalactosamine-4-SO4Receptor Displays Greater Specificity for Multivalent than Monovalent Ligands* , 2001, The Journal of Biological Chemistry.

[54]  W. März,et al.  Rationale and design of the LURIC study--a resource for functional genomics, pharmacogenomics and long-term prognosis of cardiovascular disease. , 2001, Pharmacogenomics.

[55]  R. Nagai,et al.  Endothelial dysfunction in the klotho mouse and downregulation of klotho gene expression in various animal models of vascular and metabolic diseases , 2000, Cellular and Molecular Life Sciences CMLS.

[56]  V. Luria,et al.  Maternally expressed PGK-Cre transgene as a tool for early and uniform activation of the Cre site-specific recombinase , 1998, Transgenic Research.

[57]  D. Fiete,et al.  The macrophage/endothelial cell mannose receptor cDNA encodes a protein that binds oligosaccharides terminating with SO4-4-GalNAcbeta1,4GlcNAcbeta or Man at independent sites. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[58]  C. Eun A Genome-Wide Association Study Identifies IL23R as an Inflammatory Bowel Disease Gene. , 2007 .

[59]  S. Kumar,et al.  Circulatory half-life but not interaction with the lutropin/chorionic gonadotropin receptor is modulated by sulfation of bovine lutropin oligosaccharides. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[60]  O. Hindsgaul,et al.  A hepatic reticuloendothelial cell receptor specific for SO4-4GalNAc beta 1,4GlcNAc beta 1,2Man alpha that mediates rapid clearance of lutropin. , 1991, Cell.