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.

Background: Some years ago we established an N-ethyl-N-nitrosourea screen for modifiers of transgene variegation in the mouse and a preliminary description of the first six mutant lines, named MommeD1-D6, has been published. We have reported the underlying genes in three cases: MommeD1 is a mutation in SMC hinge domain containing 1 (Smchd1), a novel modifier of epigenetic gene silencing; MommeD2 is a mutation in DNA methyltransferase 1 (Dnmt1); and MommeD4 is a mutation in Smarca 5 (Snf2h), a known chromatin remodeler. The identification of Dnmt1 and Smarca5 attest to the effectiveness of the screen design.

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.