An Atlas of Human and Murine Genetic Influences on Osteoporosis

Osteoporosis is a common debilitating chronic disease diagnosed primarily using bone mineral density (BMD). We undertook a comprehensive assessment of human genetic determinants of bone density in 426,824 individuals, identifying a total of 518 genome-wide significant loci, (301 novel), explaining 20% of the total variance in BMD—as estimated by heel quantitative ultrasound (eBMD). Next, meta-analysis identified 13 bone fracture loci in ~1.2M individuals, which were also associated with BMD. We then identified target genes from cell-specific genomic landscape features, including chromatin conformation and accessible chromatin sites, that were strongly enriched for genes known to influence bone density and strength (maximum odds ratio = 58, P = 10−75). We next performed rapid throughput skeletal phenotyping of 126 knockout mice lacking eBMD Target Genes and showed that these mice had an increased frequency of abnormal skeletal phenotypes compared to 526 unselected lines (P < 0.0001). In-depth analysis of one such Target Gene, DAAM2, showed a disproportionate decrease in bone strength relative to mineralization. This comprehensive human and murine genetic atlas provides empirical evidence testing how to link associated SNPs to causal genes, offers new insights into osteoporosis pathophysiology and highlights opportunities for drug development.

[1]  C. Greenwood,et al.  Genetic architecture: the shape of the genetic contribution to human traits and disease , 2017, Nature Reviews Genetics.

[2]  D. Kiel,et al.  Genetic Determinants of Circulating Estrogen Levels and Evidence of a Causal Effect of Estradiol on Bone Density in Men , 2018, The Journal of clinical endocrinology and metabolism.

[3]  David M. Evans,et al.  Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects. , 2018, American journal of human genetics.

[4]  Erdogan Taskesen,et al.  Functional mapping and annotation of genetic associations with FUMA , 2017, Nature Communications.

[5]  M. Pirinen,et al.  Prospects of Fine-Mapping Trait-Associated Genomic Regions by Using Summary Statistics from Genome-wide Association Studies. , 2017, American journal of human genetics.

[6]  Jacqueline K. White,et al.  Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis , 2017, Nature Genetics.

[7]  P. Donnelly,et al.  Genome-wide genetic data on ~500,000 UK Biobank participants , 2017, bioRxiv.

[8]  P. Visscher,et al.  10 Years of GWAS Discovery: Biology, Function, and Translation. , 2017, American journal of human genetics.

[9]  Job Dekker,et al.  Hi-C 2.0: An optimized Hi-C procedure for high-resolution genome-wide mapping of chromosome conformation. , 2017, Methods.

[10]  Jacques P. Brown,et al.  10 years of denosumab treatment in postmenopausal women with osteoporosis: results from the phase 3 randomised FREEDOM trial and open-label extension. , 2017, The lancet. Diabetes & endocrinology.

[11]  Borbala Mifsud,et al.  GOTHiC, a probabilistic model to resolve complex biases and to identify real interactions in Hi-C data , 2017, PloS one.

[12]  Debbie A Lawlor,et al.  Triangulation in aetiological epidemiology , 2016, International journal of epidemiology.

[13]  Houda Belaghzal,et al.  Hi-C 2.0: An Optimized Hi-C Procedure for High-Resolution Genome-Wide Mapping of Chromosome Conformation , 2016, bioRxiv.

[14]  Anthony D. Schmitt,et al.  A Compendium of Chromatin Contact Maps Reveals Spatially Active Regions in the Human Genome. , 2016, Cell reports.

[15]  J. Kent,et al.  Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR , 2016, Genome Biology.

[16]  Neva C. Durand,et al.  Juicer Provides a One-Click System for Analyzing Loop-Resolution Hi-C Experiments. , 2016, Cell systems.

[17]  F. Rivadeneira,et al.  Osteoporosis and Bone Mass Disorders: From Gene Pathways to Treatments , 2016, Trends in Endocrinology & Metabolism.

[18]  F. Cunningham,et al.  The Ensembl Variant Effect Predictor , 2016, Genome Biology.

[19]  Sayan Mukherjee,et al.  Fast Principal-Component Analysis Reveals Convergent Evolution of ADH1B in Europe and East Asia. , 2016, American journal of human genetics.

[20]  Norma E. Conner,et al.  Advances and Challenges , 2016, The American journal of hospice & palliative care.

[21]  Meagan E. Sullender,et al.  Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9 , 2015, Nature Biotechnology.

[22]  Ryan Miller,et al.  WikiPathways: capturing the full diversity of pathway knowledge , 2015, Nucleic Acids Res..

[23]  Jean-Philippe Vert,et al.  HiC-Pro: an optimized and flexible pipeline for Hi-C data processing , 2015, Genome Biology.

[24]  Wei-Chen Chen,et al.  EM Algorithm for Model-Based Clustering of Finite MixtureGaussian Distribution , 2015 .

[25]  Gabor T. Marth,et al.  A global reference for human genetic variation , 2015, Nature.

[26]  Matti Pirinen,et al.  FINEMAP: efficient variable selection using summary data from genome-wide association studies , 2015, bioRxiv.

[27]  Yakir A Reshef,et al.  Partitioning heritability by functional annotation using genome-wide association summary statistics , 2015, Nature Genetics.

[28]  F. Berenbaum,et al.  The pro-inflammatory cytokine 14-3-3ε is a ligand of CD13 in cartilage , 2015, Journal of Cell Science.

[29]  Z. Niu,et al.  WAC loss-of-function mutations cause a recognisable syndrome characterised by dysmorphic features, developmental delay and hypotonia and recapitulate 10p11.23 microdeletion syndrome , 2015, Journal of Medical Genetics.

[30]  Mulin Jun Li,et al.  Nature Genetics Advance Online Publication a N a Ly S I S the Support of Human Genetic Evidence for Approved Drug Indications , 2022 .

[31]  A. Boyde,et al.  Thyrostimulin Regulates Osteoblastic Bone Formation During Early Skeletal Development , 2015, Endocrinology.

[32]  M. Daly,et al.  An Atlas of Genetic Correlations across Human Diseases and Traits , 2015, Nature Genetics.

[33]  Lesley S. Chaboub,et al.  Daam2-PIP5K Is a Regulatory Pathway for Wnt Signaling and Therapeutic Target for Remyelination in the CNS , 2015, Neuron.

[34]  P. Elliott,et al.  UK Biobank: An Open Access Resource for Identifying the Causes of a Wide Range of Complex Diseases of Middle and Old Age , 2015, PLoS medicine.

[35]  A. Moayyeri,et al.  Predictive ability of heel quantitative ultrasound for incident fractures: an individual-level meta-analysis , 2015, Osteoporosis International.

[36]  Florian Kronenberg,et al.  EasyStrata: evaluation and visualization of stratified genome-wide association meta-analysis data , 2015, Bioinform..

[37]  Carson C Chow,et al.  Second-generation PLINK: rising to the challenge of larger and richer datasets , 2014, GigaScience.

[38]  Paul Theodor Pyl,et al.  HTSeq—a Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[39]  B. Berger,et al.  Efficient Bayesian mixed model analysis increases association power in large cohorts , 2014, Nature Genetics.

[40]  M. Daly,et al.  LD Score regression distinguishes confounding from polygenicity in genome-wide association studies , 2014, Nature Genetics.

[41]  Neva C. Durand,et al.  A 3D Map of the Human Genome at Kilobase Resolution Reveals Principles of Chromatin Looping , 2014, Cell.

[42]  J. Witte,et al.  The contribution of genetic variants to disease depends on the ruler , 2014, Nature Reviews Genetics.

[43]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[44]  Howard Y. Chang,et al.  Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position , 2013, Nature Methods.

[45]  Traver Hart,et al.  Finding the active genes in deep RNA-seq gene expression studies , 2013, BMC Genomics.

[46]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[47]  C. Esapa,et al.  Bone Mineral Content and Density. , 2012, Current protocols in mouse biology.

[48]  Nathan C. Sheffield,et al.  The accessible chromatin landscape of the human genome , 2012, Nature.

[49]  Tim D. Spector,et al.  Genetics of osteoporosis from genome-wide association studies: advances and challenges , 2012, Nature Reviews Genetics.

[50]  Alan Boyde,et al.  Rapid-Throughput Skeletal Phenotyping of 100 Knockout Mice Identifies 9 New Genes That Determine Bone Strength , 2012, PLoS genetics.

[51]  P. Visscher,et al.  Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits , 2012, Nature Genetics.

[52]  B. Deneen,et al.  Daam2 is required for dorsal patterning via modulation of canonical Wnt signaling in the developing spinal cord. , 2012, Developmental cell.

[53]  J. Bassett,et al.  Quantitative X-ray imaging of rodent bone by Faxitron. , 2012, Methods in molecular biology.

[54]  D. Kiel,et al.  Genetic Determinants of Serum Testosterone Concentrations in Men , 2011, PLoS genetics.

[55]  P. Visscher,et al.  GCTA: a tool for genome-wide complex trait analysis. , 2011, American journal of human genetics.

[56]  A. Moayyeri,et al.  Quantitative ultrasound of the heel and fracture risk assessment: an updated meta-analysis , 2011, Osteoporosis International.

[57]  Sharon R Grossman,et al.  Integrating common and rare genetic variation in diverse human populations , 2010, Nature.

[58]  A. Boskey,et al.  Contribution of Mineral to Bone Structural Behavior and Tissue Mechanical Properties , 2010, Calcified Tissue International.

[59]  P. Baldock,et al.  Critical role for Y1 receptors in mesenchymal progenitor cell differentiation and osteoblast activity , 2010, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[60]  Yun Li,et al.  METAL: fast and efficient meta-analysis of genomewide association scans , 2010, Bioinform..

[61]  C. Glass,et al.  Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. , 2010, Molecular cell.

[62]  R. H. Bassett,et al.  Optimal bone strength and mineralization requires the type 2 iodothyronine deiodinase in osteoblasts , 2010, Proceedings of the National Academy of Sciences.

[63]  L. Bonewald,et al.  Dynamics of the transition from osteoblast to osteocyte , 2010, Annals of the New York Academy of Sciences.

[64]  Aaron R. Quinlan,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2022 .

[65]  Daniel Sinnett,et al.  Population genomics in a disease targeted primary cell model. , 2009, Genome research.

[66]  Z. Ou,et al.  Microdeletions including YWHAE in the Miller–Dieker syndrome region on chromosome 17p13.3 result in facial dysmorphisms, growth restriction, and cognitive impairment , 2009, Journal of Medical Genetics.

[67]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[68]  P. Ridker,et al.  Genome-wide association studies identify novel loci associated with age at menarche and age at natural menopause , 2009, Nature Genetics.

[69]  P. Donnelly,et al.  Designing Genome-Wide Association Studies: Sample Size, Power, Imputation, and the Choice of Genotyping Chip , 2009, PLoS genetics.

[70]  A. Boskey,et al.  Fourier transform-infrared microspectroscopy and microscopic imaging. , 2008, Methods in molecular biology.

[71]  C. Tanikawa,et al.  RGC32, a novel p53-inducible gene, is located on centrosomes during mitosis and results in G2/M arrest , 2007, Oncogene.

[72]  M. Cohn,et al.  Essential Role for Co-chaperone Fkbp52 but Not Fkbp51 in Androgen Receptor-mediated Signaling and Physiology* , 2007, Journal of Biological Chemistry.

[73]  D. C. Bauer,et al.  Quantitative ultrasound predicts hip and non-spine fracture in men: the MrOS study , 2007, Osteoporosis International.

[74]  A. Silman,et al.  Predictive Value of BMD for Hip and Other Fractures , 2005, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[75]  Luke E. K. Achenie,et al.  Expression Profile of Osteoblast Lineage at Defined Stages of Differentiation* , 2005, Journal of Biological Chemistry.

[76]  Claus Christiansen,et al.  Consensus development conference: Prophylaxis and treatment of osteoporosis , 2005, Osteoporosis International.

[77]  Juliet M Chapman,et al.  Detecting Disease Associations due to Linkage Disequilibrium Using Haplotype Tags: A Class of Tests and the Determinants of Statistical Power , 2003, Human Heredity.

[78]  Kyoko Yokomori,et al.  The alternative Ctf18-Dcc1-Ctf8-replication factor C complex required for sister chromatid cohesion loads proliferating cell nuclear antigen onto DNA , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[79]  D. Kiel,et al.  Mapping of Quantitative Ultrasound of the Calcaneus Bone to Chromosome 1 by Genome-Wide Linkage Analysis , 2002, Osteoporosis International.

[80]  Andrew J. Bannister,et al.  Selective recognition of methylated lysine 9 on histone H3 by the HP1 chromo domain , 2001, Nature.

[81]  T. Spector,et al.  Genetic Variation in Bone Mineral Density and Calcaneal Ultrasound: A Study of the Influence of Menopause Using Female Twins , 2001, Osteoporosis International.

[82]  P. Sexton,et al.  Calcitonin receptor antibodies in the identification of osteoclasts. , 1999, Bone.

[83]  Harry K. Genant,et al.  Broadband ultrasound attenuation predicts fractures strongly and independently of densitometry in older women. A prospective study. Study of Osteoporotic Fractures Research Group. , 1997, Archives of internal medicine.

[84]  T. Spector,et al.  The heritability of bone mineral density, ultrasound of the calcaneus and hip axis length: A study of postmenopausal twins , 1996, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[85]  W. G. Cochran.  The combination of estimates from different experiments. , 1954 .