Integration of summary data from GWAS and eQTL studies identified novel causal BMD genes with functional predictions.
- Authors
- Meng, Xiang-He; Chen, Xiang-Ding; Greenbaum, Jonathan; Zeng, Qin; You, Sheng-Lan; Xiao, Hong-Mei; Tan, Li-Jun; Deng, Hong-Wen
- Year
- 2018
- Journal
- Bone
- PMID
- 29763751
- DOI
- 10.1016/j.bone.2018.05.012
- PMCID
- PMC6346739
PURPOSE: Osteoporosis is a common global health problem characterized by low bone mineral density (BMD) and increased risk of fracture. Genome-wide association studies (GWAS) have identified >100 genetic loci associated with BMD. However, the functional genes responsible for most associations remain largely unknown. We conducted an innovative summary statistic data-based Mendelian randomization (SMR) analysis to identify novel causal genes associated with BMD and explored their potential functional significance. METHODS: After quality control of the largest GWAS meta-analysis data of BMD and the largest expression quantitative trait loci (eQTL) meta-analysis data from peripheral blood samples, 5967 genes were tested using the SMR method. Another eQTL data was used to verify the results. Next we performed a fine-mapping association analysis to investigate the functional SNP in the identified loci. Weighted gene co-expression network analysis (WGCNA) was used to explore functional relationships for the identified novel genes with known putative osteoporosis genes. Further, we assessed functions of the identified genes through in vitro cellular study or previous functional studies. RESULTS: We identified two potentially causal genes (ASB16-AS1 and SYN2) associated with BMD. SYN2 was a novel osteoporosis candidate gene and ASB16-AS1 locus was known to be associated with BMD but was not the nearest gene to the top GWAS SNP. Fine-mapping association analysis showed that rs184478 and rs795000 was predicted to be possible causal SNPs in ASB16-AS1 and SYN2, respectively. ASB16-AS1 co-expressed with several known putative osteoporosis risk genes. In vitro cellular study showed that over-expressed ASB16-AS1 increased the expression of osteoblastogenesis related genes (BMP2 and ALPL), indicating its functional significance. CONCLUSION: Our findings support that ASB16-AS1 and SYN2 may represent two novel functional genes underlying BMD variation. The findings provide a basis for further functional mechanistic studies.
The SMR results at gene loci for BMD. (A) The SMR result at ASB16-AS1 locus for FN-BMD. (B) The SMR result at SYN2 locus for LS-BMD. In the top plot, black dots represent the p values for the SNPs from the latest GWAS meta-analysis for BMD (Y-axis), diamonds represent the p values for probes from the SMR test. In the bottom plot, the eQTL p values of the SNPs were from the eQTL study (Y-axis) for the ILMN_1676731 probe (or ILMN_1781060 probe) tagging ASB16-AS1 (or SYN2).
Estimated genetic associations and 95% confidence intervals with effect sizes in eQTL and GWAS studies for 16 genetic variants in the ASB16-AS1 gene region (A) and 15 genetic variants in the SYN2 gene region (B).
Regional association plot for ASB16-AS1 on chromosome 17. SNPs which were in this region were selected with their p values from the GWAS data of FN-BMD. r2 of pairwise LD is calculated between rs227580 and other SNPs. C17ORF65 is also known as ASB16-AS1.
The gene co-expression relationships for ASB16-AS1. ASB16-AS1 centered network provides a view of all edges and their corresponding nodes connected to ASB16-AS1 with a TOM > 0.15. We only selected those nodes proved to be associated with BMD before. Genes are color coded based on their correlation with ASB16-AS1, white (cor > 0) and grey (cor < 0).
ASB16-AS1 can promote the expression of osteoblastic genes. (A) The relative expressions of two variants of ASB16-AS1 were detected in hFOB1.19 cells. The relative expressions of ASB16-AS1 (B), ALPL (C), BMP2 (D) and RUNX2 (E) were detected after transfection with the pCEP4-ASB16-AS1 using the pCEP4 as control. Bars represented S.D. *p < 0.05, **p < 0.01.
No entities extracted from this document yet.
No uploaded files.
| Citation | PMID | DOI | Status |
|---|---|---|---|
| AmmannP, RizzoliR, Bone strength and its determinants, Osteoporos. Int 14 (Suppl. 3) (2003) S13β8.1273080010.1007/s00198-002-1345-4 | β | β | β |
| BaronR, KneisselM, WNT signaling in bone homeostasis and disease: from human mutations to treatments, Nat. Med 19 (2013) 179β192.2338961810.1038/nm.3074 | β | β | β |
| BoefAG, DekkersOM, le CessieS, Mendelian randomization studies: a review of the approaches used and the quality of reporting, Int. J. Epidemiol 44 (2015) 496β511.2595378410.1093/ije/dyv071 | β | β | β |
| BoyleAP, HongEL, HariharanM, , Annotation of functional variation in personal genomes using RegulomeDB, Genome Res. 22 (2012) 1790β1797.2295598910.1101/gr.137323.112PMC3431494 | β | β | β |
| BurgessS, DudbridgeF, ThompsonSG, Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods, Stat. Med 35 (2016) 1880β1906.2666190410.1002/sim.6835PMC4832315 | β | β | β |
| ChalmersTC, SmithHJr., BlackburnB, , A method for assessing the quality of a randomized control trial, Control. Clin. Trials 2 (1981) 31β49.726163810.1016/0197-2456(81)90056-8 | β | β | β |
| ChenG, WangC, WangJ, , Antiosteoporotic effect of icariin in ovariectomized rats is mediated via the Wnt/beta-catenin pathway, Exp. Ther. Med 12 (2016) 279β287.2734705010.3892/etm.2016.3333PMC4906828 | β | β | β |
| Davey SmithG, HemaniG, Mendelian randomization: genetic anchors for causal inference in epidemiological studies, Hum. Mol. Genet 23 (2014) R89β98.2506437310.1093/hmg/ddu328PMC4170722 | β | β | β |
| DengFY, LeiSF, ZhangY, , Peripheral blood monocyte-expressed ANXA2 gene is involved in pathogenesis of osteoporosis in humans, Mol. Cell. Proteomics 10 (2011) M111.011700.10.1074/mcp.M111.011700PMC322641121817168 | β | β | β |
| DengFY, TanLJ, ShenH, , SNP rs6265 regulates protein phosphorylation and osteoblast differentiation and influences BMD in humans, J. Bone Miner. Res 28 (2013) 2498β2507.2371240010.1002/jbmr.1997PMC4127979 | β | β | β |
| DoerksT, CopleyRR, SchultzJ, , Systematic identification of novel protein domain families associated with nuclear functions, Genome Res. 12 (2002) 47β56.1177983010.1101/gr.203201PMC155265 | β | β | β |
| EstradaK, StyrkarsdottirU, EvangelouE, , Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture, Nat. Genet 44 (2012) 491β501.2250442010.1038/ng.2249PMC3338864 | β | β | β |
| FarberCR, LusisAJ, Integrating Global Gene Expression Analysis and Genetics, 60 (2008), pp. 571β601.10.1016/S0065-2660(07)00420-8PMC310965718358333 | β | β | β |
| FarberCR, Systems genetics: a novel approach to dissect the genetic basis of osteoporosis, Curr. Osteoporos. Rep 10 (2012) 228β235.2280214610.1007/s11914-012-0112-5PMC3422742 | β | β | β |
| FarhKK-H, MarsonA, ZhuJ, , Genetic and epigenetic fine mapping of causal autoimmune disease variants, Nature 518 (2014) 337β343.2536377910.1038/nature13835PMC4336207 | β | β | β |
| GamazonER, WheelerHE, ShahKP, , A gene-based association method for mapping traits using reference transcriptome data, Nat. Genet 47 (2015) 1091β1098.2625884810.1038/ng.3367PMC4552594 | β | β | β |
| GeissmannF, ManzMG, JungS, , Development of monocytes, macrophages, and dendritic cells, Science 327 (2010) 656β661.2013356410.1126/science.1178331PMC2887389 | β | β | β |
| GhazalpourA, DossS, ZhangB, , Integrating genetic and network analysis to characterize genes related to mouse weight, PLoS Genet. preprint (2005) e130.10.1371/journal.pgen.0020130PMC155028316934000 | β | β | β |
| GiambartolomeiC, VukcevicD, SchadtEE, , Bayesian test for colocalisation between pairs of genetic association studies using summary statistics, PLoS Genet. 10 (2014) e1004383.2483039410.1371/journal.pgen.1004383PMC4022491 | β | β | β |
| GrundbergE, KwanT, GeB, , Population genomics in a disease targeted primary cell model, Genome Res. 19 (2009) 1942β1952.1965437010.1101/gr.095224.109PMC2775606 | β | β | β |
| GusevA, KoA, ShiH, , Integrative approaches for large-scale transcriptome-wide association studies, Nat. Genet 48 (2016) 245β252.2685491710.1038/ng.3506PMC4767558 | β | β | β |
| HindorffLA, SethupathyP, JunkinsHA, , Potential etiologic and functional implications of genome-wide association loci for human diseases and traits, Proc. Natl. Acad. Sci. U. S. A 106 (2009) 9362β9367.1947429410.1073/pnas.0903103106PMC2687147 | β | β | β |
| HsuYH, ZillikensMC, WilsonSG, , An integration of genome-wide association study and gene expression profiling to prioritize the discovery of novel susceptibility loci for osteoporosis-related traits, PLoS Genet. 6 (2010) e1000977.2054894410.1371/journal.pgen.1000977PMC2883588 | β | β | β |
| JiaP, ZhengS, LongJ, , dmGWAS: dense module searching for genome-wide association studies in protein-protein interaction networks, Bioinformatics 27 (2011) 95β102.2104507310.1093/bioinformatics/btq615PMC3008643 | β | β | β |
| JohnsonAD, HandsakerRE, PulitSL, , SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap, Bioinformatics 24 (2008) 2938β2939.1897417110.1093/bioinformatics/btn564PMC2720775 | β | β | β |
| KrallEA, Dawson-HughesB, Heritable and life-style determinants of bone mineral density, J. Bone Miner. Res 8 (1) (1993) 1β9.842704210.1002/jbmr.5650080102 | β | β | β |
| KrishnanV, BryantHU, MacdougaldOA, Regulation of bone mass by Wnt signaling, J. Clin. Invest 116 (2006) 1202β1209.1667076110.1172/JCI28551PMC1451219 | β | β | β |
| KwanT, GrundbergE, KokaV, , Tissue effect on genetic control of transcript isoform variation, PLoS Genet. 5 (2009) e1000608.1968054210.1371/journal.pgen.1000608PMC2719916 | β | β | β |
| LangfelderP, HorvathS, WGCNA: an R package for weighted correlation network analysis, BMC Bioinforma. 9 (2008) 559.10.1186/1471-2105-9-559PMC263148819114008 | β | β | β |
| LinX, XiongD, PengYQ, , Epidemiology and management of osteoporosis in the Peopleβs Republic of China: current perspectives, Clin. Interv. Aging 10 (2015) 1017β1033.2615070610.2147/CIA.S54613PMC4485798 | β | β | β |
| LiuYZ, DvornykV, LuY, , A novel pathophysiological mechanism for osteoporosis suggested by an in vivo gene expression study of circulating monocytes, J. Biol. Chem 280 (2005) 29011β29016.1596523510.1074/jbc.M501164200 | β | β | β |
| LiuYZ, ZhouY, ZhangL, , Attenuated monocyte apoptosis, a new mechanism for osteoporosis suggested by a transcriptome-wide expression study of monocytes, PLoS One 10 (2015) e0116792.2565907310.1371/journal.pone.0116792PMC4319757 | β | β | β |
| Lloyd-JonesLR, HollowayA, McRaeA, , The genetic architecture of gene expression in peripheral blood, Am. J. Hum. Genet 100 (2017) 228β237.2806546810.1016/j.ajhg.2016.12.008PMC5294670 | β | β | β |
| LynchMWB, Genetics and Analysis of Quantitative Traits, SinauerAssociates, 1998. | β | β | β |
| ManabeN, KawaguchiH, ChikudaH, , Connection between B lymphocyte and osteoclast differentiation pathways, J. Immunol 167 (2001) 2625β2631.1150960410.4049/jimmunol.167.5.2625 | β | β | β |
| MariePJ, DebiaisF, HayE, Regulation of human cranial osteoblast phenotype by FGF-2, FGFR-2 and BMP-2 signaling, Histol. Histopathol 17 (2002) 877β885.1216879910.14670/HH-17.877 | β | β | β |
| MeleM, FerreiraPG, ReverterF, , Human genomics. The human transcriptome across tissues and individuals, Science 348 (2015) 660β665.2595400210.1126/science.aaa0355PMC4547472 | β | β | β |
| NelsonMR, TipneyH, PainterJL, , The support of human genetic evidence for approved drug indications, Nat. Genet 47 (2015) 856β860.2612108810.1038/ng.3314 | β | β | β |
| PasaniucB, PriceAL, Dissecting the genetics of complex traits using summary association statistics, Nat. Rev. Genet 18 (2016) 117β127.2784042810.1038/nrg.2016.142PMC5449190 | β | β | β |
| PavlidesJMW, ZhuZ, GrattenJ, , Predicting gene targets from integrative analyses of summary data from GWAS and eQTL studies for 28 human complex traits, Genome Med. 8 (2016).10.1186/s13073-016-0338-4PMC497918527506385 | β | β | β |
| PocockNAEJ, HopperJL, YeatesMG, SambrookPN, EberlS, Genetic determinants of bone mass in adults.pdf, J. Clin. Invest 80 (3) (1987) 706β710.362448510.1172/JCI113125PMC442294 | β | β | β |
| RivadeneiraF, StyrkarsdottirU, EstradaK, , Twenty bone-mineral-density loci identified by large-scale meta-analysis of genome-wide association studies, Nat. Genet 41 (2009) 1199β1206.1980198210.1038/ng.446PMC2783489 | β | β | β |
| RivadeneiraF, StyrkarsdottirU, EstradaK, HalldorssonBV, , Twenty bone-mineral-density loci identified by large-scale meta-analysis of genome-wide association studies, Nat. Genet 41 (11) (2009) 1199β1206.1980198210.1038/ng.446PMC2783489 | β | β | β |
| SabikOL, FarberCR, Using GWAS to identify novel therapeutic targets for osteoporosis, Transl. Res 181 (2017) 15β26.2783764910.1016/j.trsl.2016.10.009PMC5357198 | β | β | β |
| SchadtEE, Molecular networks as sensors and drivers of common human diseases, Nature 461 (2009) 218β223.1974170310.1038/nature08454 | β | β | β |
| SchaubMA, BoyleAP, KundajeA, , Linking disease associations with regulatory information in the human genome, Genome Res. 22 (2012) 1748β1759.2295598610.1101/gr.136127.111PMC3431491 | β | β | β |
| SmithGD, EbrahimS, Mendelian randomizationβ: can genetic epidemiology contribute to understanding environmental determinants of disease? Int. J. Epidemiol 32 (2003) 1β22.1268999810.1093/ije/dyg070 | β | β | β |
| StyrkarsdottirU, HalldorssonBV, GretarsdottirS, , Multiple genetic loci for bone mineral density and fractures, N. Engl. J. Med 358 (2008) 2355β2365.1844577710.1056/NEJMoa0801197 | β | β | β |
| StyrkarsdottirU, HalldorssonBV, GretarsdottirS, , New sequence variants associated with bone mineral density, Nat. Genet 41 (2009) 15β17.1907926210.1038/ng.284 | β | β | β |
| StyrkarsdottirU, HalldorssonBV, GretarsdottirS, GudbjartssonDF, , Multiple genetic loci for bone mineral density and fractures, N. Engl. J. Med 358 (22) (2008) 2355β2365.1844577710.1056/NEJMoa0801197 | β | β | β |
| SwallowDM, PoveyS, ParkarM, , Mapping of the gene coding for the human liver/bone/kidney isozyme of alkaline phosphatase to chromosome 1, Ann. Hum. Genet 50 (1986) 229β235.344601110.1111/j.1469-1809.1986.tb01043.x | β | β | β |
| VanderWeeleTJ, Tchetgen TchetgenEJ, CornelisM, , Methodological challenges in Mendelian randomization, Epidemiology 25 (2014) 427β435.2468157610.1097/EDE.0000000000000081PMC3981897 | β | β | β |
| WeinMN, SpatzJ, NishimoriS, , HDAC5 controls MEF2C-driven sclerostin expression in osteocytes, J. Bone Miner. Res 30 (2015) 400β411.2527105510.1002/jbmr.2381PMC4342334 | β | β | β |
| WelterD, MacArthurJ, MoralesJ, , The NHGRI GWAS catalog, a curated resource of SNP-trait associations, Nucleic Acids Res. 42 (2014) D1001β6.2431657710.1093/nar/gkt1229PMC3965119 | β | β | β |
| WestraHJ, PetersMJ, EskoT, , Systematic identification of trans eQTLs as putative drivers of known disease associations, Nat. Genet 45 (2013) 1238β1243.2401363910.1038/ng.2756PMC3991562 | β | β | β |
| WittrantY, TheoleyreS, ChipoyC, , RANKL/RANK/OPG: new therapeutic targets in bone tumours and associated osteolysis, Biochim. Biophys. Acta 1704 (2004) 49β57.1536386010.1016/j.bbcan.2004.05.002 | β | β | β |
| WrightNC, LookerAC, SaagKG, , The recent prevalence of osteoporosis and low bone mass in the United States based on bone mineral density at the femoral neck or lumbar spine, J. Bone Miner. Res 29 (2014) 2520β2526.2477149210.1002/jbmr.2269PMC4757905 | β | β | β |
| YangJE, SongMS, ShenY, , The role of KV7.3 in regulating osteoblast maturation and mineralization, Int. J. Mol. Sci 17 (2016) 407.2699912810.3390/ijms17030407PMC4813262 | β | β | β |
| ZhangL, ChoiHJ, EstradaK, , Multistage genome-wide association meta-analyses identified two new loci for bone mineral density, Hum. Mol. Genet 23 (2014) 1923β1933.2424974010.1093/hmg/ddt575PMC3943521 | β | β | β |
| ZhengHF, ForgettaV, HsuYH, , Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture, Nature 526 (2015) 112β117.2636779410.1038/nature14878PMC4755714 | β | β | β |
| ZhuZ, ZhangF, HuH, , Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets, Nat. Genet 48 (2016) 481β487.2701911010.1038/ng.3538 | β | β | β |
In this knowledge base
| Title | Year | PMID |
|---|---|---|
| Multi-omics integration analysis identifies novel genes for alcoholism with potential overlap with neurodegenerative diseases. | 2021 | 34417470 |
External
| Title | Authors | Journal | Year | Link |
|---|---|---|---|---|
| Identification of novel causally related genes in adenomyosis: An integrated summary data-based Mendelian randomization study and bioinformatics analysis. | Yang Q et al. | β | 2025 | β |
| Integrative analysis of genome and transcriptome reveals a novel regulator for pork intramuscular fat content. | Zhao X et al. | β | 2025 | β |
| Causal relationship between key genes and metabolic dysfunction-associated fatty liver disease risk mediated by immune cells: A Mendelian randomization and mediation analysis. | Feng G et al. | β | 2024 | β |
| The Parkinson's disease risk gene cathepsin B promotes fibrillar alpha-synuclein clearance, lysosomal function and glucocerebrosidase activity in dopaminergic neurons. | Jones-Tabah J et al. | β | 2024 | β |
| The molecular pathology of schizophrenia: an overview of existing knowledge and new directions for future research. | Nakamura T et al. | β | 2023 | β |
| The Parkinsonβs disease risk gene cathepsin B promotes fibrillar alpha-synuclein clearance, lysosomal function and glucocerebrosidase activity in dopaminergic neurons | Jones-Tabah J et al. | β | 2023 | β |
| Ten challenges for clinical translation in psychiatric genetics. | Derks EM et al. | β | 2022 | β |
| Convergent lines of evidence support BIN1 as a risk gene of Alzheimer's disease. | Zhu J et al. | β | 2021 | β |
| Expression, prognosis value, and immune infiltration of lncRNA ASB16-AS1 identified by pan-cancer analysis. | Wu L et al. | β | 2021 | β |
| Identifying Thyroid Carcinoma-Related Genes by Integrating GWAS and eQTL Data. | Shen F et al. | β | 2021 | β |
| Multi-omics integration analysis identifies novel genes for alcoholism with potential overlap with neurodegenerative diseases. | Kapoor M et al. | β | 2021 | β |
| Translating genetic association of lipid levels for biological and clinical application. | Crone B et al. | β | 2021 | β |
| A novel computational strategy for DNA methylation imputation using mixture regression model (MRM). | Yu F et al. | β | 2020 | β |
| Approximation of bias and mean-squared error in two-sample Mendelian randomization analyses. | Deng L et al. | β | 2020 | β |
| A road map for understanding molecular and genetic determinants of osteoporosis. | Yang TL et al. | β | 2020 | β |
| Comprehensive Analysis of the Genetic and Epigenetic Mechanisms of Osteoporosis and Bone Mineral Density. | Dong H et al. | β | 2020 | β |
| Identification of Gingivitis-Related Genes Across Human Tissues Based on the Summary Mendelian Randomization. | Zhang J et al. | β | 2020 | β |
| Dissecting the Genetics of Osteoporosis using Systems Approaches. | Al-Barghouthi BM et al. | β | 2019 | β |
| Identification of crucial miRNAs and lncRNAs for ossification of ligamentum flavum. | Kong D et al. | β | 2019 | β |
| Integrate GWAS, eQTL, and mQTL Data to Identify Alzheimer's Disease-Related Genes. | Zhao T et al. | β | 2019 | β |
| Integrating Multi-Omics Data to Identify Novel Disease Genes and Single-Neucleotide Polymorphisms. | Zhao S et al. | β | 2019 | β |
| SNP eQTL status and eQTL density in the adjacent region of the SNP are associated with its statistical significance in GWA studies. | Gorlov I et al. | β | 2019 | β |
| Identification of Alzheimer's Disease-Related Genes Based on Data Integration Method. | Hu Y et al. | β | 2018 | β |