Targeted sequencing of genome wide significant loci associated with bone mineral density (BMD) reveals significant novel and rare variants: the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) targeted sequencing study.

Background Bone mineral density (BMD) is a heritable phenotype that predicts fracture risk. We performed fine-mapping by targeted sequencing at WLS, MEF2C, ARHGAP1/F2 and JAG1 loci prioritized by eQTL and bioinformatic approaches among 56 BMD loci from our previous GWAS meta-analysis. Targeted sequencing was conducted in 1,291 Caucasians from the Framingham Heart Study ( n =  925) and Cardiovascular Health Study ( n =  366), including 206 women and men with extreme low femoral neck (FN) BMD. A total of 4,964 sequence variants (SNVs) were observed and 80% were rare with MAF <1%. The associations between previously identified SNPs in these loci and BMD, while nominally significant in sequenced participants, were no longer significant after multiple testing corrections. Conditional analyses did not find protein-coding variants that may be responsible for GWAS signals. On the other hand, in the sequenced subjects, we identified novel associations in WLS , ARHGAP1 , and 5' of MEF2C ( P- values < 8x10  -   5 ; false discovery rate (FDR) q-values < 0.01) that were much more strongly associated with BMD compared to the GWAS SNPs. These associated SNVs are less-common; independent from previous GWAS signals in the same loci; and located in gene regulatory elements. Our findings suggest that protein-coding variants in selected GWAS loci did not contribute to GWAS signals. By performing targeted sequencing in GWAS loci, we identified less-common and rare non-coding SNVs associated with BMD independently from GWAS common SNPs, suggesting both common and less-common variants may associate with disease risks and phenotypes in the same loci.

[1]  Jason H. Moore,et al.  Chapter 11: Genome-Wide Association Studies , 2012, PLoS Comput. Biol..

[2]  Issei Komuro,et al.  Wnt signaling and aging-related heart disorders. , 2010, Circulation research.

[3]  S. Nelson,et al.  BFAST: An Alignment Tool for Large Scale Genome Resequencing , 2009, PloS one.

[4]  K. Lunetta,et al.  Methods in Genetics and Clinical Interpretation Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Design of Prospective Meta-Analyses of Genome-Wide Association Studies From 5 Cohorts , 2010 .

[5]  Claus Christiansen,et al.  Diagnosis of Osteoporosis , 1992, Southern medical journal.

[6]  O Johnell,et al.  Meta-analysis of how well measures of bone mineral density predict occurrence of osteoporotic fractures. , 1996, BMJ.

[7]  Jacob A. Tennessen,et al.  Evolution and Functional Impact of Rare Coding Variation from Deep Sequencing of Human Exomes , 2012, Science.

[8]  C. Antoniades,et al.  Wnt signaling in cardiovascular physiology , 2012, Trends in Endocrinology & Metabolism.

[9]  Steven R Goldring,et al.  Eating bone or adding it: the Wnt pathway decides , 2007, Nature Medicine.

[10]  David Z. Chen,et al.  Architecture of the human regulatory network derived from ENCODE data , 2012, Nature.

[11]  Christian Gieger,et al.  New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk , 2010, Nature Genetics.

[12]  W. Hsu,et al.  Gpr177, a novel locus for bone mineral density and osteoporosis, regulates osteogenesis and chondrogenesis in skeletal development , 2013, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[13]  Daniel L. Koller,et al.  Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture , 2012, Nature Genetics.

[14]  Tom R. Gaunt,et al.  Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci. , 2012, American journal of human genetics.

[15]  D. Kiel,et al.  Clinical review: Genome-wide association studies of skeletal phenotypes: what we have learned and where we are headed. , 2012, The Journal of clinical endocrinology and metabolism.

[16]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[17]  R. Nusse,et al.  The Wnt signaling pathway in development and disease. , 2004, Annual review of cell and developmental biology.

[18]  J. Pritchard,et al.  Genomics: ENCODE explained , 2012, Nature.

[19]  Data production leads,et al.  An integrated encyclopedia of DNA elements in the human genome , 2012 .

[20]  Tanya M. Teslovich,et al.  Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index , 2010 .

[21]  Shane J. Neph,et al.  Systematic Localization of Common Disease-Associated Variation in Regulatory DNA , 2012, Science.

[22]  Kenny Q. Ye,et al.  An integrated map of genetic variation from 1,092 human genomes , 2012, Nature.

[23]  ENCODEConsortium,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[24]  M. Marazita,et al.  Genome-wide Association Studies , 2012, Journal of dental research.

[25]  Xihong Lin,et al.  Rare-variant association testing for sequencing data with the sequence kernel association test. , 2011, American journal of human genetics.

[26]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[27]  Eric Boerwinkle,et al.  Strategies to Design and Analyze Targeted Sequencing Data: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study , 2014, Circulation. Cardiovascular genetics.

[28]  Bjarni V. Halldórsson,et al.  Association of JAG1 with bone mineral density and osteoporotic fractures: a genome-wide association study and follow-up replication studies. , 2010, American journal of human genetics.

[29]  F. Natt,et al.  Control of the SOST Bone Enhancer by PTH Using MEF2 Transcription Factors , 2007, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[30]  E. Creemers,et al.  Myocardin is a direct transcriptional target of Mef2, Tead and Foxo proteins during cardiovascular development , 2006, Development.

[31]  Nicole Soranzo,et al.  An Integration of Genome-Wide Association Study and Gene Expression Profiling to Prioritize the Discovery of Novel Susceptibility Loci for Osteoporosis-Related Traits , 2010, PLoS genetics.

[32]  C. Bucana,et al.  Control of mouse cardiac morphogenesis and myogenesis by transcription factor MEF2C. , 1997, Science.

[33]  John J. Andreucci,et al.  Protein Kinase A-regulated Assembly of a MEF2·HDAC4 Repressor Complex Controls c-Jun Expression in Vascular Smooth Muscle Cells* , 2009, The Journal of Biological Chemistry.

[34]  Eurie L. Hong,et al.  Annotation of functional variation in personal genomes using RegulomeDB , 2012, Genome research.

[35]  J. Stamatoyannopoulos,et al.  Power of deep, all-exon resequencing for discovery of human trait genes , 2009, Proceedings of the National Academy of Sciences.

[36]  Christian Gieger,et al.  Genetic Predictors of Fibrin D-Dimer Levels in Healthy Adults , 2011, Circulation.

[37]  S. Lipton,et al.  Assignment of human myocyte-specific enhancer binding factor 2C (hMEF2C) to human chromosome 5q14 and evidence that MEF2C is evolutionarily conserved. , 1995, Genomics.

[38]  Manolis Kellis,et al.  HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants , 2011, Nucleic Acids Res..

[39]  E. Oetjen,et al.  Genome-Wide Association Identifies Nine Common Variants Associated With Fasting Proinsulin Levels and Provides New Insights Into the Pathophysiology of Type 2 Diabetes , 2011, Diabetes.

[40]  S. Hannenhalli,et al.  Evidence for Coregulation of Myocardial Gene Expression by MEF2 and NFAT in Human Heart Failure , 2009, Circulation. Cardiovascular genetics.