Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We... Show moreSkull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases. Show less
A large-scale GWAS provides insight on diabetes-dependent genetic effects on the glomerular filtration rate, a common metric to monitor kidney health in disease.Reduced glomerular filtration rate ... Show moreA large-scale GWAS provides insight on diabetes-dependent genetic effects on the glomerular filtration rate, a common metric to monitor kidney health in disease.Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (n(DM) = 178,691, n(noDM) = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM. Show less
Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions... Show moreCommon single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries. Show less
Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here... Show moreGlycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 x 10(-8)), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.A trans-ancestry meta-analysis of GWAS of glycemic traits in up to 281,416 individuals identifies 99 novel loci, of which one quarter was found due to the multi-ancestry approach, which also improves fine-mapping of credible variant sets. Show less
Ahluwalia, T.S.; Prins, B.P.; Abdollahi, M.; Armstrong, N.J.; Aslibekyan, S.; Bain, L.; ... ; CHARGE Inflammation Working Grp 2021
Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been... Show moreInterleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67428 (n(discovery)=52654 and n(replication)=14774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (P-combined=1.8x10(-11)), HLA-DRB1/DRB5 rs660895 on Chr6p21 (P-combined=1.5x10(-10)) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (P-combined=1.2x10(-122)). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology. Show less
Context: There is a need for novel biomarkers and better understanding of the pathophysiology of diabetic kidney disease.Objective: To investigate associations between plasma metabolites and kidney... Show moreContext: There is a need for novel biomarkers and better understanding of the pathophysiology of diabetic kidney disease.Objective: To investigate associations between plasma metabolites and kidney function in people with type 2 diabetes (T2D).Design: 3089 samples from individuals with T2D, collected between 1999 and 2015, from 5 independent Dutch cohort studies were included. Up to 7 years follow-up was available in 1100 individuals from 2 of the cohorts.Main outcome measures: Plasma metabolites (n = 149) were measured by nuclear magnetic resonance spectroscopy. Associations between metabolites and estimated glomerular filtration rate (eGFR), urinary albumin-to-creatinine ratio (UACR), and eGFR slopes were investigated in each study followed by random effect meta-analysis. Adjustments included traditional cardiovascular risk factors and correction for multiple testing.Results: In total, 125 metabolites were significantly associated (P-FDR = 1.5x10(-32) - 0.046; beta = -11.98-2.17) with eGFR. Inverse associations with eGFR were demonstrated for branched-chain and aromatic amino acids (AAAs), glycoprotein acetyls, triglycerides (TGs), lipids in very low-density lipoproteins (VLDL) subclasses, and fatty acids (P-FDR < 0.03). We observed positive associations with cholesterol and phospholipids in high-density lipoproteins (HDL) and apolipoprotein A1 (P-FDR < 0.05). Albeit some metabolites were associated with UACR levels (P < 0.05), significance was lost after correction for multiple testing. Tyrosine and HDL-related metabolites were positively associated with eGFR slopes before adjustment for multiple testing (P-Tyr = 0.003; P-HDLrelated < 0.05), but not after.Conclusions: This study identified metabolites associated with impaired kidney function in T2D, implying involvement of lipid and amino acid metabolism in the pathogenesis. Whether these processes precede or are consequences of renal impairment needs further investigation. Show less
Epidemiology studies suggested that low birthweight was associated with a higher risk of hypertension in later life. However, little is known about the causality of such associations. In our study,... Show moreEpidemiology studies suggested that low birthweight was associated with a higher risk of hypertension in later life. However, little is known about the causality of such associations. In our study, we evaluated the causal association of low birthweight with adulthood hypertension following a standard analytic protocol using the study-level data of 183,433 participants from 60 studies (CHARGE-BIG consortium), as well as that with blood pressure using publicly available summary-level genome-wide association data from EGG consortium of 153,781 participants, ICBP consortium and UK Biobank cohort together of 757,601 participants. We used seven SNPs as the instrumental variable in the study-level analysis and 47 SNPs in the summary-level analysis. In the study-level analyses, decreased birthweight was associated with a higher risk of hypertension in adults (the odds ratio per 1 standard deviation (SD) lower birthweight, 1.22; 95% CI 1.16 to 1.28), while no association was found between genetically instrumented birthweight and hypertension risk (instrumental odds ratio for causal effect per 1 SD lower birthweight, 0.97; 95% CI 0.68 to 1.41). Such results were consistent with that from the summary-level analyses, where the genetically determined low birthweight was not associated with blood pressure measurements either. One SD lower genetically determined birthweight was not associated with systolic blood pressure (beta = - 0.76, 95% CI - 2.45 to 1.08 mmHg), 0.06 mmHg lower diastolic blood pressure (beta = - 0.06, 95% CI - 0.93 to 0.87 mmHg), or pulse pressure (beta = - 0.65, 95% CI - 1.38 to 0.69 mmHg, all p > 0.05). Our findings suggest that the inverse association of birthweight with hypertension risk from observational studies was not supported by large Mendelian randomization analyses. Show less
Merino, J.; Dashti, H.S.; Li, S.X.; Sarnowski, C.; Justice, A.E.; Graff, M.; ... ; Tanaka, T. 2019
BACKGROUNDOnly a handful of genetic discovery efforts in apparent treatment-resistant hypertension (aTRH) have been described.METHODSWe conducted a case-control genome-wide association study of... Show moreBACKGROUNDOnly a handful of genetic discovery efforts in apparent treatment-resistant hypertension (aTRH) have been described.METHODSWe conducted a case-control genome-wide association study of aTRH among persons treated for hypertension, using data from 10 cohorts of European ancestry (EA) and 5 cohorts of African ancestry (AA). Cases were treated with 3 different antihypertensive medication classes and had blood pressure (BP) above goal (systolic BP >= 140 mm Hg and/or diastolic BP >= 90 mm Hg) or 4 or more medication classes regardless of BP control (n(EA) = 931, n(AA) = 228). Both a normotensive control group and a treatment-responsive control group were considered in separate analyses. Normotensive controls were untreated (n(EA) = 14,210, n(AA) = 2,480) and had systolic BP/diastolic BP < 140/90 mm Hg. Treatment-responsive controls (n(EA) = 5,266, n(AA) = 1,817) had BP at goal (<140/90 mm Hg), while treated with one antihypertensive medication class. Individual cohorts used logistic regression with adjustment for age, sex, study site, and principal components for ancestry to examine the association of single-nucleotide polymorphisms with case-control status. Inverse variance-weighted fixed-effects meta-analyses were carried out using METAL.RESULTSThe known hypertension locus, CASZ1, was a top finding among EAs (P = 1.1 x 10(-8)) and in the race-combined analysis (P = 1.5 x 10(-9)) using the normotensive control group (rs12046278, odds ratio = 0.71 (95% confidence interval: 0.6-0.8)). Single-nucleotide polymorphisms in this locus were robustly replicated in the Million Veterans Program (MVP) study in consideration of a treatment-responsive control group. There were no statistically significant findings for the discovery analyses including treatment-responsive controls.CONCLUSIONThis genomic discovery effort for aTRH identified CASZ1 as an aTRH risk locus. Show less
Wuttke, M.; Li, Y.; Li, M.; Sieber, K.B.; Feitosa, M.F.; Gorski, M.; ... ; Waterwort 2019