A longitudinal genome-wide association study of bone mineral density mean and variability in the UK Biobank

Summary Bone mineral density (BMD) is an essential predictor of osteoporosis and fracture. We conducted a genome-wide trajectory analysis of BMD and analyzed the BMD change. Purpose This study aimed to identify the genetic architecture and potential biomarkers of BMD. Methods Our analysis included 141,261 white participants from the UK Biobank with heel BMD phenotype data. We used a genome-wide trajectory analysis tool, TrajGWAS, to conduct a genome-wide association study (GWAS) of BMD. Then, we validated our findings in previously reported BMD genetic associations and performed replication analysis in the Asian participants. Finally, gene-set enrichment analysis (GSEA) of the identified candidate genes was conducted using the FUMA platform. Results A total of 52 genes associated with BMD trajectory mean were identified, of which the top three significant genes were WNT16 ( P = 1.31 × 10 −126 ), FAM3C ( P = 4.18 × 10 −108 ), and CPED1 ( P = 8.48 × 10 −106 ). In addition, 114 genes associated with BMD within-subject variability were also identified, such as AC092079.1 ( P = 2.72 × 10 −13 ) and RGS7 ( P = 4.72 × 10 −10 ). The associations for these candidate genes were confirmed in the previous GWASs and replicated successfully in the Asian participants. GSEA results of BMD change identified multiple GO terms related to skeletal development, such as SKELETAL SYSTEM DEVELOPMENT ( P adjusted = 2.45 × 10 −3 ) and REGULATION OF OSSIFICATION ( P adjusted = 2.45 × 10 −3 ). KEGG enrichment analysis showed that these genes were mainly enriched in WNT SIGNALING PATHWAY. Conclusions Our findings indicated that the CPED1 - WNT16 - FAM3C locus plays a significant role in BMD mean trajectories and identified several novel candidate genes contributing to BMD within-subject variability, facilitating the understanding of the genetic architecture of BMD.