Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Ling Oei - Oei, Karol Estrada Gil, EL Duncan, C Christiansen, CT Liu, BL Langdahl, B Obermayer-Pietsch, JA Riancho, RL Prince, NM Schoor, E McCloskey, YH Hsu, E Evangelou, E Ntzani, DM Evans, N Alonso, LB Husted, C Valero, JL Hernandez, JR LewisSK Kaptoge, K Zhu, LA Cupples, C Medina-Gomez, L Vandenput, GS Kim, SH Lee, MC Castano-Betancourt, Edwin Oei, J Martinez, A Daroszewska, M Klift, D Mellstrom, L Herrera, MK Karlsson, Bert Hofman, O Ljunggren, Huib Pols, Lisette Stolk, Joyce van Meurs, JPA Ioannidis, M.C. Zillikens, P Lips, D Karasik, André Uitterlinden, U Styrkarsdottir, MA Brown, JM Koh, JB Richards, J Reeve, C Ohlsson, SH Ralston, DP Kiel, Fernando Rivadeneira

Research output: Contribution to journalArticleAcademicpeer-review

27 Citations (Scopus)


Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fracture applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55 years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5 x 10(-8). In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p = 4.6 x 10(-8). However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% Cl: 0.98-1.14; p = 0.17), displaying high degree of heterogeneity (I-2= 57%; Q(het)p = 0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p = 0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures. (C) 2013 Elsevier Inc. All rights reserved.
Original languageUndefined/Unknown
Pages (from-to)20-27
Number of pages8
Publication statusPublished - 2014

Cite this