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.

Mek is a dual-specificity kinase that activates the extracellular-signal-regulated (Erk) mitogen-activated protein (MAP) kinases upon agonist binding to receptors. The Erk MAP kinase cascade is involved in cell-fate determination in many organisms. In mammals, this pathway is proposed to regulate cell growth and differentiation. Genetic studies have shown that although a single mek gene is present in Caenorhabditis elegans, Drosophila and Xenopus, two mek homologs, Mek1 and Mek2, are present in the mammalian cascade. In the present study, we describe a mutant mouse line in which the mek1 gene has been disrupted by insertional mutagenesis. The null mutation was recessive lethal, as the homozygous mutant embryos died at 10.5 days of gestation. Histopathological analyses revealed a reduction in vascularization of the placenta that was due to a marked decrease of vascular endothelial cells in the labyrinthine region. The failure to establish a functional placenta probably explains the death of the mek1 -/embryos. Cell-migration assays indicated that mek1 -/fibroblasts could not be induced to migrate by fibronectin, although the levels of Mek2 protein and Erk activation were normal. Re-expression of Mek1 in the mutant mouse embryonic fibroblasts (MEFs) restored their ability to migrate. Our findings provide genetic evidence that establishes the unique role played by Mek1 in signal transduction. They also suggest that mek1 function is required for normal response to angiogenic signals that might promote vascularization of the labyrinthine region of the placenta.