Osteoporosis is a common aging-related disease diagnosed primarily using bone mineral density (BMD). We assessed genetic determinants of BMD as estimated by heel quantitative ultrasound (eBMD) in 426,824 individuals, identifying 518 genome-wide significant loci (301 novel), explaining 20% of its variance. We identified 13 bone fracture loci, all associated with eBMD, in ~1.2M individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds-ratio=58, p=10 −75 ) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice lacking target genes and found an increased abnormal skeletal phenotype frequency compared to 526 unselected lines (p<0.0001). In-depth analysis of one gene, DAAM2 , showed a disproportionate decrease in bone strength relative to mineralization. This genetic atlas provides evidence testing how to link associated-SNPs to causal genes, offers new insights into osteoporosis pathophysiology and highlights opportunities for drug development.

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study in 142,487 individuals from the UK Biobank to identify loci associated with BMD estimated by quantitative ultrasound of the heel (“eBMD”). We identified 307 conditionally independent SNPs attaining genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 novel loci, and several rare variants with large effect sizes. To investigate underlying mechanisms we undertook: 1) bioinformatic, functional genomic annotation and human osteoblast expression studies; 2) gene function prediction; 3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and 4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. These studies strongly implicate GPC6 as a novel determinant of BMD and also identify abnormal skeletal phenotypes in knockout mice for a further 100 prioritized genes.

Summary Osteoclasts are large multinucleated bone-resorbing cells formed by the fusion of monocyte/macrophage-derived precursors that are thought to undergo apoptosis once resorption is complete. Here, by intravital imaging, we reveal that RANKL-stimulated osteoclasts have an alternative cell fate in which they fission into daughter cells called osteomorphs. Inhibiting RANKL blocked this cellular recycling and resulted in osteomorph accumulation. Single-cell RNA sequencing showed that osteomorphs are transcriptionally distinct from osteoclasts and macrophages and express a number of non-canonical osteoclast genes that are associated with structural and functional bone phenotypes when deleted in mice. Furthermore, genetic variation in human orthologs of osteomorph genes causes monogenic skeletal disorders and associates with bone mineral density, a polygenetic skeletal trait. Thus, osteoclasts recycle via osteomorphs, a cell type involved in the regulation of bone resorption that may be targeted for the treatment of skeletal diseases.

Achondroplasia (ACH) is the most frequent form of dwarfism and is caused by gain-of-function mutations in the fibroblast growth factor receptor 3-encoding (FGFR3-encoding) gene. Although potential therapeutic strategies for ACH, which aim to reduce excessive FGFR3 activation, have emerged over many years, the use of tyrosine kinase inhibitor (TKI) to counteract FGFR3 hyperactivity has yet to be evaluated. Here, we have reported that the pan-FGFR TKI, NVP-BGJ398, reduces FGFR3 phosphorylation and corrects the abnormal femoral growth plate and calvaria in organ cultures from embryos of the Fgfr3Y367C/+ mouse model of ACH. Moreover, we demonstrated that a low dose of NVP-BGJ398, injected subcutaneously, was able to penetrate into the growth plate of Fgfr3Y367C/+ mice and modify its organization. Improvements to the axial and appendicular skeletons were noticeable after 10 days of treatment and were more extensive after 15 days of treatment that started from postnatal day 1. Low-dose NVP-BGJ398 treatment reduced intervertebral disc defects of lumbar vertebrae, loss of synchondroses, and foramen-magnum shape anomalies. NVP-BGJ398 inhibited FGFR3 downstream signaling pathways, including MAPK, SOX9, STAT1, and PLCγ, in the growth plates of Fgfr3Y367C/+ mice and in cultured chondrocyte models of ACH. Together, our data demonstrate that NVP-BGJ398 corrects pathological hallmarks of ACH and support TKIs as a potential therapeutic approach for ACH.

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10−22) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10−13) and osteoarthritis (P = 1.6 × 10−7). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease.