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

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×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×10−13; SOX6, p = 6.4×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.

[1]  F. Dudbridge,et al.  Estimation of significance thresholds for genomewide association scans , 2008, Genetic epidemiology.

[2]  P. Allhoff,et al.  The Framingham Offspring Study , 1991 .

[3]  M. Olivier.  A haplotype map of the human genome , 2003, Nature.

[4]  Bjarni V. Halldórsson,et al.  New sequence variants associated with bone mineral density , 2009, Nature Genetics.

[5]  E. Schadt,et al.  Identification of Abcc6 as the major causal gene for dystrophic cardiac calcification in mice through integrative genomics , 2007, Proceedings of the National Academy of Sciences.

[6]  Joshua T. Burdick,et al.  In silico method for inferring genotypes in pedigrees , 2006, Nature Genetics.

[7]  Yurii S. Aulchenko,et al.  Twenty bone mineral density loci identified by large-scale meta-analysis of genome-wide association studies , 2009, Nature Genetics.

[8]  H. Barden,et al.  Femur strength index predicts hip fracture independent of bone density and hip axis length , 2005, Osteoporosis International.

[9]  S. Horvath,et al.  Variations in DNA elucidate molecular networks that cause disease , 2008, Nature.

[10]  K. Roeder,et al.  Genomic Control for Association Studies , 1999, Biometrics.

[11]  J. Seidman,et al.  A mouse model of human familial hypocalciuric hypercalcemia and neonatal severe hyperparathyroidism , 1995, Nature Genetics.

[12]  D. Kiel,et al.  Collaborative Meta-analysis: Associations of 150 Candidate Genes With Osteoporosis and Osteoporotic Fracture , 2009, Annals of Internal Medicine.

[13]  W. Kannel,et al.  An investigation of coronary heart disease in families. The Framingham offspring study. , 1979, American journal of epidemiology.

[14]  Robert W. Williams,et al.  Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function , 2005, Nature Genetics.

[15]  M. Almeida,et al.  Kinase-mediated regulation of common transcription factors accounts for the bone-protective effects of sex steroids. , 2003, The Journal of clinical investigation.

[16]  H. Stefánsson,et al.  Genetics of gene expression and its effect on disease , 2008, Nature.

[17]  J. Castle,et al.  An integrative genomics approach to infer causal associations between gene expression and disease , 2005, Nature Genetics.

[18]  M. Zaidi,et al.  Impaired osteoblastic differentiation, reduced bone formation, and severe osteoporosis in noggin-overexpressing mice. , 2003, The Journal of clinical investigation.

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

[20]  A. Uitterlinden,et al.  Differential genetic effects of ESR1 gene polymorphisms on osteoporosis outcomes. , 2004, JAMA.

[21]  A Hofman,et al.  Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study , 2008, The Lancet.

[22]  Christopher B. Ruff,et al.  Sex differences in geometry of the femoral neck with aging: A structural analysis of bone mineral data , 2004, Calcified Tissue International.

[23]  S. Ralston,et al.  Genetic regulation of bone mass and susceptibility to osteoporosis. , 2006, Genes & development.

[24]  Eric E Schadt,et al.  An Integrative Genetics Approach to Identify Candidate Genes Regulating BMD: Combining Linkage, Gene Expression, and Association , 2009, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[25]  T. Niu,et al.  Large‐Scale Genome‐Wide Linkage Analysis for Loci Linked to BMD at Different Skeletal Sites in Extreme Selected Sibships , 2006, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[26]  David M. Evans,et al.  Common variants in the region around Osterix are associated with bone mineral density and growth in childhood , 2009, Human molecular genetics.

[27]  John D. Storey.  The positive false discovery rate: a Bayesian interpretation and the q-value , 2003 .

[28]  John D. Storey,et al.  Mapping the Genetic Architecture of Gene Expression in Human Liver , 2008, PLoS biology.

[29]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[30]  K. Nakashima,et al.  [The Rotterdam study]. , 2011, Nihon rinsho. Japanese journal of clinical medicine.

[31]  C. Hoggart,et al.  Genome-wide association analysis of metabolic traits in a birth cohort from a founder population , 2008, Nature Genetics.

[32]  Tim D Spector,et al.  The UK Adult Twin Registry (TwinsUK) , 2006, Twin Research and Human Genetics.

[33]  D. Karasik,et al.  Contribution of Gender‐Specific Genetic Factors to Osteoporosis Risk , 2008, Annals of human genetics.

[34]  Gregor Eichele,et al.  GenePaint.org: an atlas of gene expression patterns in the mouse embryo , 2004, Nucleic Acids Res..

[35]  A. Hinnebusch,et al.  IMPACT, a Protein Preferentially Expressed in the Mouse Brain, Binds GCN1 and Inhibits GCN2 Activation* , 2005, Journal of Biological Chemistry.

[36]  G. Abecasis,et al.  Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies , 2006, Nature Genetics.

[37]  S. Ferrari,et al.  Beta-arrestin2 regulates parathyroid hormone effects on a p38 MAPK and NFkappaB gene expression network in osteoblasts. , 2009, Bone.

[38]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[39]  Benjamin M. Bolstad,et al.  affy - analysis of Affymetrix GeneChip data at the probe level , 2004, Bioinform..

[40]  S. K. Lee,et al.  Printed in U.S.A. Copyright © 1999 by The Endocrine Society Parathyroid Hormone Stimulates TRANCE and Inhibits Osteoprotegerin Messenger Ribonucleic Acid Expression in Murine Bone Marrow Cultures: Correlation with , 2022 .

[41]  A. Masuda,et al.  JNK Activity Is Essential for Atf4 Expression and Late‐Stage Osteoblast Differentiation , 2009, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[42]  R. Collins,et al.  Common variants at 30 loci contribute to polygenic dyslipidemia , 2009, Nature Genetics.

[43]  T. Manolio,et al.  How to Interpret a Genome-wide Association Study Topic Collections , 2022 .

[44]  K. Mossman.  The Wellcome Trust Case Control Consortium, U.K. , 2008 .

[45]  T. Zhu,et al.  Src-CrkII-C3G-dependent Activation of Rap1 Switches Growth Hormone-stimulated p44/42 MAP Kinase and JNK/SAPK Activities* , 2003, Journal of Biological Chemistry.

[46]  Richard Baldock,et al.  Computational genetics: Bioinformatics beyond sequence: mapping gene function in the embryo , 2001, Nature Reviews Genetics.

[47]  Scott A. Rifkin,et al.  Revealing the architecture of gene regulation: the promise of eQTL studies. , 2008, Trends in genetics : TIG.

[48]  A. Uitterlinden,et al.  The Association between Common Vitamin D Receptor Gene Variations and Osteoporosis: A Participant-Level Meta-Analysis , 2006, Annals of Internal Medicine.

[49]  Michael Boehnke,et al.  Efficient study designs for test of genetic association using sibship data and unrelated cases and controls. , 2006, American journal of human genetics.

[50]  P. Deloukas,et al.  Common Regulatory Variation Impacts Gene Expression in a Cell Type–Dependent Manner , 2009, Science.

[51]  Simon C. Potter,et al.  Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls , 2007, Nature.

[52]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[53]  T. Dawber,et al.  Epidemiological approaches to heart disease: the Framingham Study. , 1951, American journal of public health and the nation's health.

[54]  Monique M. B. Breteler,et al.  The Rotterdam Study: 2016 objectives and design update , 2015, European Journal of Epidemiology.

[55]  M. Daly,et al.  Identifying Relationships among Genomic Disease Regions: Predicting Genes at Pathogenic SNP Associations and Rare Deletions , 2009, PLoS genetics.

[56]  Kari Stefansson,et al.  Multiple genetic loci for bone mineral density and fractures. , 2008, The New England journal of medicine.

[57]  L. Liang,et al.  A genome-wide association study of global gene expression , 2007, Nature Genetics.

[58]  L. Liang,et al.  Mapping complex disease traits with global gene expression , 2009, Nature Reviews Genetics.

[59]  K. Sirotkin,et al.  The NCBI dbGaP database of genotypes and phenotypes , 2007, Nature Genetics.

[60]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[61]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[62]  David Botstein,et al.  GO: : TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes , 2004, Bioinform..

[63]  T. Jämsä,et al.  Experimental hip fracture load can be predicted from plain radiography by combined analysis of trabecular bone structure and bone geometry , 2008, Osteoporosis International.

[64]  Judy H. Cho,et al.  Finding the missing heritability of complex diseases , 2009, Nature.

[65]  Josée Dupuis,et al.  Genome-wide association with bone mass and geometry in the Framingham Heart Study , 2007, BMC Medical Genetics.

[66]  G. Abecasis,et al.  Merlin—rapid analysis of dense genetic maps using sparse gene flow trees , 2002, Nature Genetics.

[67]  D. Kiel,et al.  Risk Factors for Longitudinal Bone Loss in Elderly Men and Women: The Framingham Osteoporosis Study , 2000, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[68]  D. Kiel,et al.  Large-scale analysis of association between LRP5 and LRP6 variants and osteoporosis. , 2008, JAMA.

[69]  A. Uitterlinden,et al.  Tissue Effect on Genetic Control of Transcript Isoform Variation , 2009, PLoS genetics.

[70]  Daniel L. Koller,et al.  Meta‐Analysis of Genome‐Wide Scans Provides Evidence for Sex‐ and Site‐Specific Regulation of Bone Mass , 2006, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[71]  M. Olivier.  A haplotype map of the human genome. , 2003, Nature.

[72]  T. Niu,et al.  Variation in genes involved in the RANKL/RANK/OPG bone remodeling pathway are associated with bone mineral density at different skeletal sites in men , 2005, Human Genetics.

[73]  K Yano,et al.  Osteoclast differentiation factor is a ligand for osteoprotegerin/osteoclastogenesis-inhibitory factor and is identical to TRANCE/RANKL. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[74]  G. Abecasis,et al.  Association Analysis in a Variance Components Framework , 2001, Genetic epidemiology.

[75]  Jehyuk Lee,et al.  A Robust Approach to Identifying Tissue-Specific Gene Expression Regulatory Variants Using Personalized Human Induced Pluripotent Stem Cells , 2009, PLoS genetics.

[76]  S. Cummings,et al.  Bone mass measurements and risk of fracture in Caucasian women: a review of findings from prospective studies. , 1995, The American journal of medicine.

[77]  Yan Guo,et al.  Genome-wide association and follow-up replication studies identified ADAMTS18 and TGFBR3 as bone mass candidate genes in different ethnic groups. , 2009, American journal of human genetics.

[78]  Richard L Prince,et al.  In vivo short-term precision of hip structure analysis variables in comparison with bone mineral density using paired dual-energy X-ray absorptiometry scans from multi-center clinical trials. , 2005, Bone.

[79]  Fernando Rivadeneira,et al.  Femoral Neck BMD Is a Strong Predictor of Hip Fracture Susceptibility in Elderly Men and Women Because It Detects Cortical Bone Instability: The Rotterdam Study , 2007, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[80]  E. Seeman,et al.  Clinical review 137: Sexual dimorphism in skeletal size, density, and strength. , 2001, The Journal of clinical endocrinology and metabolism.

[81]  Jehyuk Lee,et al.  Digital RNA Allelotyping Reveals Tissue-specific and Allele-specific Gene Expression in Human , 2009, Nature Methods.

[82]  R. Pereira,et al.  CCAAT/enhancer binding protein homologous protein (DDIT3) induces osteoblastic cell differentiation. , 2004, Endocrinology.

[83]  P. Donnelly,et al.  A new multipoint method for genome-wide association studies by imputation of genotypes , 2007, Nature Genetics.

[84]  M. Daly,et al.  Evaluating and improving power in whole-genome association studies using fixed marker sets , 2006, Nature Genetics.

[85]  Enrico Petretto,et al.  Heritability and Tissue Specificity of Expression Quantitative Trait Loci , 2006, PLoS genetics.

[86]  J. Stockman.  Genetic Variants Associated With Cardiac Structure and Function: A Meta-analysis and Replication of Genome-wide Association Data , 2011 .

[87]  G. Stein,et al.  Relationship of cell growth to the regulation of tissue‐specific gene expression during osteoblast differentiation , 1990, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[88]  John P A Ioannidis,et al.  Required sample size and nonreplicability thresholds for heterogeneous genetic associations , 2008, Proceedings of the National Academy of Sciences.

[89]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[90]  J. Malvehy,et al.  Genome-wide association study identifies three loci associated with melanoma risk , 2009, Nature Genetics.