Pathway analysis of genome-wide association study for bone mineral density
Authors:
- 25 April 2012
- journal article
- research article
- Published by Springer Nature in Molecular Biology Reports
- Vol. 39 (8) , 8099-8106
Abstract
The aim of this study was to identify the candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms that contribute to bone mineral density (BMD) and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach. We used hip BMD GWAS data of the genotypes of 301,019 SNPs in 5,715 Europeans. ICSNPathway (identify candidate causal SNPs and pathways) analysis was applied to the BMD GWAS dataset. The first stage involved the pre-selection of candidate causal SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate causal SNPs using improved-gene set enrichment analysis. ICSNPathway analysis identified seven candidate SNPs, eight candidate pathways, and seven hypothetical biological mechanisms. Eight pathways are as follows; gamma-hexachlorocyclohexane degradation (nominal p-value < 0.001, false discovery rate (FDR) p-value < 0.001, FDR = 0.016), TACI and BCMA stimulation of B cell immune response (nominal p-value < 0.001, FDR = 0.021), endonuclease activity (nominal p-value = 0.001, FDR = 0,026), regulation of defense response to virus (nominal p-value = 0.001, FDR = 0.028), serine_type_endopeptidase_inhibitor_activity (nominal p-value = 0.001, FDR = 0.044), endoribonuclease activity (nominal p-value = 0.002, FDR = 0.045), and myeloid leukocyte differentiation (nominal p-value = 0.001, FDR = 0.050). The most significant causal pathway was gamma-hexachlorocyclohexane degradation. CYP3A5, PON2, PON3, CMBL, PON1, ALPL, CYP3A43, CYP3A7, ACP6, ACPP, and ALPI (p < 0.05) are involved in the pathway of gamma-hexachlorocyclohexane degradation. Further examination of the gene contents revealed that DBR1, DICER1, EXO1, FEN1, POP1, POP4, RPP30, and RPP38 were involved in 2 of the 8 pathways (p < 0.05). By applying ICSNPathway analysis to BMD GWAS data, we identified seven candidate SNPs and eight pathways involving gamma-hexachlorocyclohexane degradation, which may contribute to low BMD.Keywords
This publication has 23 references indexed in Scilit:
- Pathway-based analysis of GWAS datasets: effective but caution requiredInternational Journal of Neuropsychopharmacology, 2011
- Analysing biological pathways in genome-wide association studiesNature Reviews Genetics, 2010
- Molecular networks as sensors and drivers of common human diseasesNature, 2009
- Associations between osteoprotegerin polymorphisms and bone mineral density: a meta-analysisMolecular Biology Reports, 2009
- Strategies and issues in the detection of pathway enrichment in genome-wide association studiesHuman Genetics, 2009
- SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMapBioinformatics, 2008
- A Gene Wiki for Community Annotation of Gene FunctionPLOS Biology, 2008
- Association of Systemic Lupus Erythematosus withC8orf13–BLKandITGAM–ITGAXNew England Journal of Medicine, 2008
- Searching for genes underlying susceptibility to osteoporotic fracture: current progress and future prospectOsteoporosis International, 2007
- Genetic Epidemiological Approaches to the Search for Osteoporosis GenesJournal of Bone and Mineral Research, 2000