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.

Background The role of vitamin D status in COVID-19 patients is a matter of debate. Objectives To assess serum 25-hydroxyvitamin D (25OHD) levels in hospitalized patients with COVID-19 and to analyze the possible influence of vitamin D status on disease severity. Methods Retrospective case–control study of 216 COVID-19 patients and 197 population-based controls. Serum 25OHD levels were measured in both groups. The association of serum 25OHD levels with COVID-19 severity (admission to the intensive care unit, requirements for mechanical ventilation, or mortality) was also evaluated. Results Of the 216 patients, 19 were on vitamin D supplements and were analyzed separately. In COVID-19 patients, mean ± standard deviation 25OHD levels were 13.8 ± 7.2 ng/mL, compared with 20.9 ± 7.4 ng/mL in controls (P &lt; .0001). 25OHD values were lower in men than in women. Vitamin D deficiency was found in 82.2% of COVID-19 cases and 47.2% of population-based controls (P &lt; .0001). 25OHD inversely correlates with serum ferritin (P = .013) and D-dimer levels (P = .027). Vitamin D-deficient COVID-19 patients had a greater prevalence of hypertension and cardiovascular diseases, raised serum ferritin and troponin levels, as well as a longer length of hospital stay than those with serum 25OHD levels ≥20 ng/mL. No causal relationship was found between vitamin D deficiency and COVID-19 severity as a combined endpoint or as its separate components. Conclusions 25OHD levels are lower in hospitalized COVID-19 patients than in population-based controls and these patients had a higher prevalence of deficiency. We did not find any relationship between vitamin D concentrations or vitamin deficiency and the severity of the disease.