Bone mineral density (BMD) is the most important predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and East Asian ancestry. We tested the top-associated BMD markers for replication in 50,933 independent subjects and for risk of low-trauma fracture in 31,016 cases and 102,444 controls. We identified 56 loci (32 novel)associated with BMD atgenome-wide significant level (P<5×10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal-stem-cell differentiation, endochondral ossification and the Wnt signalling pathways. However, we also discovered loci containing genes not known to play a role in bone biology. Fourteen BMD loci were also associated with fracture risk (P<5×10−4, Bonferroni corrected), of which six reached P<5×10−8 including: 18p11.21 (C18orf19), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

SUMMARY The extent to which low-frequency (minor allele frequency [MAF] between 1–5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is largely unknown. Bone mineral density (BMD) is highly heritable, is a major predictor of osteoporotic fractures and has been previously associated with common genetic variants1–8, and rare, population-specific, coding variants9. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n=2,882 from UK10K), whole-exome sequencing (n= 3,549), deep imputation of genotyped samples using a combined UK10K/1000Genomes reference panel (n=26,534), and de-novo replication genotyping (n= 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size 4-fold larger than the mean of previously reported common variants for lumbar spine BMD8 (rs11692564[T], MAF = 1.7%, replication effect size = +0.20 standard deviations [SD], Pmeta = 2×10−14), which was also associated with a decreased risk of fracture (OR = 0.85; P = 2×10−11; ncases = 98,742 and ncontrols = 409,511). Using an En1Cre/flox mouse model, we observed that conditional loss of En1 results in low bone mass, likely as a consequence of high bone turn-over. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817[T], MAF = 1.1%, replication effect size = +0.39 SD, Pmeta = 1×10−11). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Paget’s disease of bone (PDB) is a common disorder with a strong genetic component characterised by focal increases in bone turnover which in some cases is caused by SQSTM1 mutations. To identify additional susceptibility genes we performed a genome wide association study in 750 PDB cases without SQSTM1 mutations and 1002 controls and identified three candidate loci for the disease which were replicated in an independent set of 500 cases and 535 controls. The strongest signal was with rs484959 on 1p13 close to the CSF1 gene (P = 5.38 × 10−24) and significant associations were also observed with rs1561570 on 10p13 within the OPTN gene (P = 6.09 × 10−13) and with rs3018362 on 18q21 close to the TNFRSF11A gene (P = 5.27 × 10−13). These studies provide new insights into the pathogenesis of PDB and identify OPTN, CSF1 and TNFRSF11A as novel candidate genes for disease susceptibility.

Paget's disease of bone (PDB) is a common disorder characterized by focal abnormalities of bone remodeling. We previously identified variants at the CSF1, OPTN and TNFRSF11A loci as risk factors for PDB by genome-wide association study. Here we extended this study, identified three new loci and confirmed their association with PDB in 2,215 affected individuals (cases) and 4,370 controls from seven independent populations. The new associations were with rs5742915 within PML on 15q24 (odds ratio (OR) = 1.34, P = 1.6 × 10(-14)), rs10498635 within RIN3 on 14q32 (OR = 1.44, P = 2.55 × 10(-11)) and rs4294134 within NUP205 on 7q33 (OR = 1.45, P = 8.45 × 10(-10)). Our data also confirmed the association of TM7SF4 (rs2458413, OR = 1.40, P = 7.38 × 10(-17)) with PDB. These seven loci explained ∼13% of the familial risk of PDB. These studies provide new insights into the genetic architecture and pathophysiology of PDB.