Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2019 Feb;51(2):258-266.
doi: 10.1038/s41588-018-0302-x. Epub 2018 Dec 31.

An atlas of genetic influences on osteoporosis in humans and mice

John A Morris  1   2 John P Kemp  3   4 Scott E Youlten  5 Laetitia Laurent  2 John G Logan  6 Ryan C Chai  5 Nicholas A Vulpescu  7 Vincenzo Forgetta  2 Aaron Kleinman  8 Sindhu T Mohanty  5 C Marcelo Sergio  5 Julian Quinn  5 Loan Nguyen-Yamamoto  9 Aimee-Lee Luco  9 Jinchu Vijay  10 Marie-Michelle Simon  10 Albena Pramatarova  10 Carolina Medina-Gomez  11 Katerina Trajanoska  11 Elena J Ghirardello  6 Natalie C Butterfield  6 Katharine F Curry  6 Victoria D Leitch  6 Penny C Sparkes  6 Anne-Tounsia Adoum  6 Naila S Mannan  6 Davide S K Komla-Ebri  6 Andrea S Pollard  6 Hannah F Dewhurst  6 Thomas A D Hassall  3 Michael-John G Beltejar  12 23andMe Research TeamDouglas J Adams  13 Suzanne M Vaillancourt  14 Stephen Kaptoge  15 Paul Baldock  5 Cyrus Cooper  16   17   18 Jonathan Reeve  18 Evangelia E Ntzani  19   20 Evangelos Evangelou  19   21 Claes Ohlsson  22 David Karasik  23 Fernando Rivadeneira  11 Douglas P Kiel  23   24   25   26 Jonathan H Tobias  27 Celia L Gregson  27 Nicholas C Harvey  16   17 Elin Grundberg  10   28 David Goltzman  9 David J Adams  29 Christopher J Lelliott  29 David A Hinds  8 Cheryl L Ackert-Bicknell  30 Yi-Hsiang Hsu  23   24   25   26 Matthew T Maurano  7 Peter I Croucher  5 Graham R Williams  6 J H Duncan Bassett  6 David M Evans  31   32 J Brent Richards  33   34   35   36   37
Collaborators, Affiliations

An atlas of genetic influences on osteoporosis in humans and mice

John A Morris et al. Nat Genet. 2019 Feb.

Erratum in

  • Author Correction: An atlas of genetic influences on osteoporosis in humans and mice.
    Morris JA, Kemp JP, Youlten SE, Laurent L, Logan JG, Chai RC, Vulpescu NA, Forgetta V, Kleinman A, Mohanty ST, Sergio CM, Quinn J, Nguyen-Yamamoto L, Luco AL, Vijay J, Simon MM, Pramatarova A, Medina-Gomez C, Trajanoska K, Ghirardello EJ, Butterfield NC, Curry KF, Leitch VD, Sparkes PC, Adoum AT, Mannan NS, Komla-Ebri DSK, Pollard AS, Dewhurst HF, Hassall TAD, Beltejar MG; 23andMe Research Team; Adams DJ, Vaillancourt SM, Kaptoge S, Baldock P, Cooper C, Reeve J, Ntzani EE, Evangelou E, Ohlsson C, Karasik D, Rivadeneira F, Kiel DP, Tobias JH, Gregson CL, Harvey NC, Grundberg E, Goltzman D, Adams DJ, Lelliott CJ, Hinds DA, Ackert-Bicknell CL, Hsu YH, Maurano MT, Croucher PI, Williams GR, Bassett JHD, Evans DM, Richards JB. Morris JA, et al. Nat Genet. 2019 May;51(5):920. doi: 10.1038/s41588-019-0415-x. Nat Genet. 2019. PMID: 30988516

Abstract

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 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 estimated BMD (eBMD), in ~1.2 million individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds ratio (OR) = 58, P = 1 × 10-75) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice with disruptions in predicted 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 linking associated SNPs to causal genes, offers new insight into osteoporosis pathophysiology, and highlights opportunities for drug development.

PubMed Disclaimer

Conflict of interest statement

Competing Interests Statement

A.K. and D.A.H. are employees of 23andMe, Inc.

Figures

Figure 1.
Figure 1.. Manhattan plot of genome-wide association results for eBMD in the UK Biobank.
The dashed red line denotes the threshold for declaring genome-wide significance (6.6×10−9). 1,103 conditionally independent SNPs at 515 loci passed the criteria for genome-wide significance in n=426,824 UK Biobank participants. 301 novel loci (defined as > 1 Mbp from previously reported genome-wide significant BMD variants) reaching genome-wide significance are displayed in blue. Previously reported loci that reached genome-wide significance are displayed in red, and previously reported loci failing to reach genome-wide significance in our study are shown in black.
Figure 2.
Figure 2.. Fine-mapping SNPs and target gene selection diagram.
a) For each 500 Mbp region around a conditionally independent lead SNP (p<6.6×10−9 after conditional independence testing; n=426,824 UK Biobank participants) we applied statistical fine-mapping to calculate log10 Bayes factors for each SNP as a measure of their posterior probability for causality. Conditional independence testing was implemented using GCTA-COJO, and log10 Bayes factors were estimated using FINEMAP. SNPs that were conditionally independent lead SNPs or that had log10 Bayes factors > 3 were considered our fine-mapped SNPs that we then used for target gene identification. b) Target Genes were identified if: 1) It was the gene closest to a fine-mapped SNP. 2) A fine-mapped SNP was in its gene body. 3) A fine-mapped SNP was coding. 4) The gene mapped closest to a fine-mapped SNP which resided in an SaOS-2 ATAC-seq peak. 5) A fine-mapped SNP was present in a Hi-C osteoblast or osteocyte promoter interaction peak, therefore being closer to a target gene in three-dimensions than linearly on the genome.
Figure 3.
Figure 3.. SNPs at genome-wide significant loci are enriched for bone-relevant open chromatin sites.
Comparison of eBMD-associated SNPs in terms of enrichment for DHSs from primary osteoblasts, and ATAC-seq peaks from SaOS-2 osteosarcoma cells. Odds ratios were computed relative to all SNPs at genome-wide significant regions. Enrichments for missense protein coding SNPs are shown as baselines. a) Enrichments for conditionally independent (COJO) or log10 Bayes factor >3 (FINEMAP); note the latter set contains nearly twice the number of SNPs. b) Ranking SNPs by log10 Bayes factor (x-axis) showed increasing enrichment. 95% confidence interval (shaded region) was calculated by a two-sided Fisher’s Exact Test.
Figure 4.
Figure 4.
Target Gene Identification Workflow.
Figure 5.
Figure 5.. Reduction of DAAM2 protein resulted in reduced mineralization in SaOS-2 cells.
Mineralization quantification in control cells and DAAM2 exon 2 double-stranded break (DSB) induced cells in either the presence of osteogenic factors (treated) or absence (untreated). a) Dot plot of n=6 independent experiments ± standard error of the mean (SEM) from Alizarin red staining in (b) to quantify mineralization; Bar=5mm. ***p=1.3×10−15 compared to untreated control cells and &&&p=9.3×10−15 (left) and 8.2×10−13 (right) compared to treated control cells determined by one-way ANOVA (F=49.7, df=5) and Bonferroni post-hoc tests.
Figure 6.
Figure 6.. Biomechanical analyses of mice with Daam2 knockdown.
a) Femur biomechanical analysis. Destructive 3-point bend testing (Instron 5543 load frame) of femurs from wild-type (WT, nfemale=3, nmale=4), Daam2+/tm1a (nfemale=6, nmale=4) and Daam2tm1a/tm1a (nfemale=5, nmale=9) mice. Graphs show yield load, maximum load, fracture load, stiffness (gradient of the linear elastic phase) and toughness (energy dissipated prior to fracture). Female data are shown on the left and male data on the right. Data are shown as mean ± standard error of the mean (SEM). Female maximum load analyses for WT versus Daam2tm1a/tm1a (**) and Daam2+/tm1a versus Daam2tm1a/tm1a (#) had statistically significant differences (one-way ANOVA p=3.0×10−3, F=10.29, df=13, Tukey’s post-hoc test **p<0.01 and #p<0.05). Male maximum load analyses for WT versus Daam2tm1a/tm1a (***) and Daam2+/tm1a versus Daam2tm1a/tm1a had statistically significant differences [one-way ANOVA p<1.0×10−4 (GraphPad Prism does not report smaller p-values), F=50.11, df=16, Tukey’s post-hoc test ***p<1.0×10−3 and ###p<1.0×10−3]. Male fracture load analyses for WT vs Daam2tm1a/tm1a (***) and Daam2+/tm1a vs Daam2tm1a/tm1 (##) had statistically significant differences (one-way ANOVA p=3.0×10−4, F=15.49, df=16, Tukey’s post-hoc test ***p<1.0×10−3 and ##p<0.01). b) Vertebra biomechanical analyses. Destructive compression testing (Instron 5543 load frame) of caudal vertebrae from WT (nfemale=3, nmale=4), Daam2+/tm1a (nfemale=6, nmale=4) and Daam2tm1a/tm1a (nfemale=5, nmale=9) mice. Graphs show yield load, maximum load and stiffness. Data are shown as mean ± SEM. Female yield load analysis for WT versus Daam2tm1a/tm1a (**) had a statistically significant difference (one-way ANOVA p=6.5×10−3, F=8.26, df=13, Tukey’s post-hoc test **p<0.01). Female maximum load analyses for WT versus Daam2tm1a/tm1a (**) and WT versus Daam2+/tm1a (*) had statistically significant differences (one-way ANOVA p=2.9×10−3, F=10.45, df=13, Tukey’s post-hoc test **p<0.01 and *p<0.05). Male maximum load analysis for WT vs Daam2tm1a/tm1a (*) had a statistically significant difference (one-way ANOVA p=0.04, F=4.10, df=16, Tukey’s post-hoc test *p<0.05). c) Bone quality analysis from rapid throughput screening mouse knockouts. The graph demonstrates the physiological relationship between bone mineral content and stiffness in caudal vertebrae from P112 female WT mice (n=320). The blue line shows the linear regression (Pearson’s r=0.21, p=1.2×10−4) and the grey box indicates ± 2 standard deviations (SD). The mean value for female Daam2tm1a/tm1a [n=2 from initial OBCD screen (Supplementary Note)] mice is shown in orange (−2.14 SD).

Comment in

  • Unravelling the genetics of osteoporosis.
    Greenhill C. Greenhill C. Nat Rev Endocrinol. 2019 Mar;15(3):129. doi: 10.1038/s41574-019-0158-x. Nat Rev Endocrinol. 2019. PMID: 30647468 No abstract available.
  • GWAS cracks fracture risk.
    Trenkmann M. Trenkmann M. Nat Rev Rheumatol. 2019 Mar;15(3):126. doi: 10.1038/s41584-019-0173-2. Nat Rev Rheumatol. 2019. PMID: 30697000 No abstract available.

Similar articles

  • Twelve New Genomic Loci Associated With Bone Mineral Density.
    Liu L, Zhao M, Xie ZG, Liu J, Peng HP, Pei YF, Sun HP, Zhang L. Liu L, et al. Front Endocrinol (Lausanne). 2020 Apr 22;11:243. doi: 10.3389/fendo.2020.00243. eCollection 2020. Front Endocrinol (Lausanne). 2020. PMID: 32390946 Free PMC article.
  • Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium.
    Moayyeri A, Hsu YH, Karasik D, Estrada K, Xiao SM, Nielson C, Srikanth P, Giroux S, Wilson SG, Zheng HF, Smith AV, Pye SR, Leo PJ, Teumer A, Hwang JY, Ohlsson C, McGuigan F, Minster RL, Hayward C, Olmos JM, Lyytikäinen LP, Lewis JR, Swart KM, Masi L, Oldmeadow C, Holliday EG, Cheng S, van Schoor NM, Harvey NC, Kruk M, del Greco M F, Igl W, Trummer O, Grigoriou E, Luben R, Liu CT, Zhou Y, Oei L, Medina-Gomez C, Zmuda J, Tranah G, Brown SJ, Williams FM, Soranzo N, Jakobsdottir J, Siggeirsdottir K, Holliday KL, Hannemann A, Go MJ, Garcia M, Polasek O, Laaksonen M, Zhu K, Enneman AW, McEvoy M, Peel R, Sham PC, Jaworski M, Johansson Å, Hicks AA, Pludowski P, Scott R, Dhonukshe-Rutten RA, van der Velde N, Kähönen M, Viikari JS, Sievänen H, Raitakari OT, González-Macías J, Hernández JL, Mellström D, Ljunggren O, Cho YS, Völker U, Nauck M, Homuth G, Völzke H, Haring R, Brown MA, McCloskey E, Nicholson GC, Eastell R, Eisman JA, Jones G, Reid IR, Dennison EM, Wark J, Boonen S, Vanderschueren D, Wu FC, Aspelund T, Richards JB, Bauer D, Hofman A, Khaw KT, Dedoussis G, Obermayer-Pietsch B, Gyllensten U, Pramstaller PP, Lorenc RS, Cooper C, Kung AW, Lips P, Alen M, Attia J, Brandi ML, de Groot … See abstract for full author list ➔ Moayyeri A, et al. Hum Mol Genet. 2014 Jun 1;23(11):3054-68. doi: 10.1093/hmg/ddt675. Epub 2014 Jan 14. Hum Mol Genet. 2014. PMID: 24430505 Free PMC article.
  • Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis.
    Kemp JP, Morris JA, Medina-Gomez C, Forgetta V, Warrington NM, Youlten SE, Zheng J, Gregson CL, Grundberg E, Trajanoska K, Logan JG, Pollard AS, Sparkes PC, Ghirardello EJ, Allen R, Leitch VD, Butterfield NC, Komla-Ebri D, Adoum AT, Curry KF, White JK, Kussy F, Greenlaw KM, Xu C, Harvey NC, Cooper C, Adams DJ, Greenwood CMT, Maurano MT, Kaptoge S, Rivadeneira F, Tobias JH, Croucher PI, Ackert-Bicknell CL, Bassett JHD, Williams GR, Richards JB, Evans DM. Kemp JP, et al. Nat Genet. 2017 Oct;49(10):1468-1475. doi: 10.1038/ng.3949. Epub 2017 Sep 4. Nat Genet. 2017. PMID: 28869591 Free PMC article.
  • Dissecting the Genetics of Osteoporosis using Systems Approaches.
    Al-Barghouthi BM, Farber CR. Al-Barghouthi BM, et al. Trends Genet. 2019 Jan;35(1):55-67. doi: 10.1016/j.tig.2018.10.004. Epub 2018 Nov 20. Trends Genet. 2019. PMID: 30470485 Free PMC article. Review.
  • Associations of Smoking and Alcohol and Coffee Intake with Fracture and Bone Mineral Density: A Mendelian Randomization Study.
    Yuan S, Michaëlsson K, Wan Z, Larsson SC. Yuan S, et al. Calcif Tissue Int. 2019 Dec;105(6):582-588. doi: 10.1007/s00223-019-00606-0. Epub 2019 Sep 4. Calcif Tissue Int. 2019. PMID: 31482193 Review.

Cited by

References

    1. World Health Organization. Consensus development conference: Prophylaxis and treatment of osteoporosis. Osteoporos. Int. 1, 114–117 (1991). - PubMed
    1. Richards JB, Zheng H-F & Spector TD Genetics of osteoporosis from genome-wide association studies: advances and challenges. Nat. Rev. Genet. 13, 576–588 (2012). - PubMed
    1. Johnell O et al. Predictive value of BMD for hip and other fractures. J. Bone Miner. Res. 20, 1185–1194 (2005). - PubMed
    1. Kemp JP et al. Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis. Nat. Genet. 49, 1468–1475 (2017). - PMC - PubMed
    1. Arden NK, Baker J, Hogg C, Baan K & Spector TD The heritability of bone mineral density, ultrasound of the calcaneus and hip axis length: a study of postmenopausal twins. J. Bone Miner. Res. 11, 530–534 (1996). - PubMed

Publication types

-