Genetic Variants in WNT16 and PKD2L1 Locus Affect Heel Ultrasound Bone Stiffness: Analyses from the General Population and Patients Evaluated for Osteoporosis

Osteoporosis, a complex chronic disease with increasing prevalence, is characterised by reduced bone mineral density (BMD) and increased fracture risk. The high heritability of BMD suggests substantial impact of the individual genetic disposition on bone phenotypes and the development of osteoporosis. In the past years, genome-wide association studies (GWAS) identified hundreds of genetic variants associated with BMD or osteoporosis. Here, we analysed 1103 single nucleotide polymorphisms (SNPs), previously identified as associated with estimated BMD (eBMD) in the UK Biobank. We assessed whether these SNPs are related to heel stiffness index obtained by quantitative ultrasound in 5665 adult participants of the Study of Health in Pomerania (SHIP). We confirmed 45 significant associations after correction for multiple testing. Next, we analysed six selected SNPs in 631 patients evaluated for osteoporosis [rs2707518 (CPED1/WNT16), rs3779381 (WNT16), rs115242848 (LOC101927709/EN1), rs10239787 (JAZF1), rs603424 (PKD2L1) and rs6968704 (JAZF1)]. Differences in minor allele frequencies (MAF) of rs2707518 and rs3779381 between SHIP participants (higher MAF) and patients evaluated for osteoporosis (lower MAF) indicated a protective effect of the minor allele on bone integrity. In contrast, differences in MAF of rs603424 indicated a harmful effect. Co-localisation analyses indicated that the rs603424 effect may be mediated via stearoyl-CoA desaturase (SCD) expression, an enzyme highly expressed in adipose tissue with a crucial role in lipogenesis. Taken together, our results support the role of the WNT16 pathway in the regulation of bone properties and indicate a novel causal role of SCD expression in adipose tissue on bone integrity.