Genetic Variation in and Its Association with Bone Mineral Density, Fractures and Osteoporosis in Children with Bone Fragility

Several genome-wide association studies (GWAS), GWAS meta-analyses, and mouse studies have demonstrated that wingless-related integration site 16 (WNT16) gene is associated with bone mineral density (BMD), cortical bone thickness, bone strength and fracture risk. Practically no data exist regarding the significance of WNT16 in childhood-onset osteoporosis and related fractures. We hypothesized that pathogenic variants and genetic variations in could explain skeletal fragility in affected children. We screened the gene by Sanger sequencing in three pediatric cohorts: 35 with primary osteoporosis, 59 with multiple fractures, and in 95 healthy controls. Altogether, we identified 12 variants in . Of them one was a rare 5'UTR variant rs1386898215 in genome aggregate and medical trans-omic databases (GnomAD, TOPMED; minor allele frequency (MAF) 0.00 and 0.000008, respectively). One variant rs1554366753, overrepresented in children with osteoporosis (MAF = 0.06 healthy controls MAF = 0.01), was significantly associated with lower BMD. This variant was found associated with increased gene expression at mRNA level in fibroblast cultures. None of the other identified variants were rare (MAF < 0.001) or deemed pathogenic by predictor programs. WNT16 may play a role in childhood osteoporosis but genetic variation is not a common cause of skeletal fragility in childhood.