An Integration of Genome-Wide Association Study and Gene Expression Profiling to Prioritize the Discovery of Novel Susceptibility Loci for Osteoporosis-Related Traits

YH Hsu, M.C. Zillikens, SG Wilson, CR Farber, S Demissie, N Soranzo, EN Bianchi, E Grundberg, LM Liang, JB Richards, Karol Estrada Gil, Y (Yanhua) Zhou, A van Nas, MF Moffatt, G Zhai, Bert Hofman, Joyce van Meurs, Huib Pols, RI Price, O NilssonT Pastinen, LA Cupples, AJ Lusis, EE Schadt, S Ferrari, André Uitterlinden, Fernando Rivadeneira, TD Spector, D Karasik, DP Kiel

Research output: Contribution to journalArticleAcademicpeer-review

164 Citations (Scopus)
1 Downloads (Pure)

Abstract

Osteoporosis is a complex disorder and commonly leads to fractures in elderly persons. Genome-wide association studies (GWAS) have become an unbiased approach to identify variations in the genome that potentially affect health. However, the genetic variants identified so far only explain a small proportion of the heritability for complex traits. Due to the modest genetic effect size and inadequate power, true association signals may not be revealed based on a stringent genome-wide significance threshold. Here, we take advantage of SNP and transcript arrays and integrate GWAS and expression signature profiling relevant to the skeletal system in cellular and animal models to prioritize the discovery of novel candidate genes for osteoporosis-related traits, including bone mineral density (BMD) at the lumbar spine (LS) and femoral neck (FN), as well as geometric indices of the hip (femoral neck-shaft angle, NSA; femoral neck length, NL; and narrow-neck width, NW). A two-stage meta-analysis of GWAS from 7,633 Caucasian women and 3,657 men, revealed three novel loci associated with osteoporosis-related traits, including chromosome 1p13.2 (RAP1A, p = 3.6 x 10(-8)), 2q11.2 (TBC1D8), and 18q11.2 (OSBPL1A), and confirmed a previously reported region near TNFRSF11B/OPG gene. We also prioritized 16 suggestive genome-wide significant candidate genes based on their potential involvement in skeletal metabolism. Among them, 3 candidate genes were associated with BMD in women. Notably, 2 out of these 3 genes (GPR177, p = 2.6 x 10(-13); SOX6, p = 6.4 x 10(-10)) associated with BMD in women have been successfully replicated in a large-scale meta-analysis of BMD, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of our prioritization strategy. In the absence of direct biological support for identified genes, we highlighted the efficiency of subsequent functional characterization using publicly available expression profiling relevant to the skeletal system in cellular or whole animal models to prioritize candidate genes for further functional validation.
Original languageUndefined/Unknown
Article numbere1000977
JournalPLoS Genetics (print)
Volume6
Issue number6
DOIs
Publication statusPublished - 2010

Cite this