骨质疏松症易感基因BDNF的遗传学关联分析及功能研究

doi:10.16288/j.yczz.17-145

遗传 ›› 2017, Vol. 39 ›› Issue (8): 726-736.doi: 10.16288/j.yczz.17-145

骨质疏松症易感基因BDNF的遗传学关联分析及功能研究

杨熳(),卢冰婕,段媛媛,陈晓峰,马建岗,郭燕()

西安交通大学生命科学与技术学院,生物医学信息工程教育部重点实验室,西安 710049

收稿日期:2017-04-20 修回日期:2017-06-08 出版日期:2017-08-20 发布日期:2017-12-25
作者简介:杨熳,硕士研究生,研究方向：遗传与分子生物学。E-mail: ymx6800@163.com|郭燕,博士,副教授,研究方向：骨质疏松症遗传与代谢。E-mail: guoyan253@mail.xjtu.edu.cn
基金资助:
国家自然科学基金项目(31471188, 81573241);中国博士后科学基金项目(2016M602797);陕西省自然科学基础研究项目(2016JQ3026)

Genetics association study and functional analysis on osteoporosis susceptibility gene BDNF

Yang Man(),Lu Bingjie,Duan Yuanyuan,Chen Xiaofeng,Ma Jiangang,Guo Yan()

Key Laboratory of Biomedical Information Engineering of Ministry of Education, and Institute of Molecular Genetics, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an 710049, China

Received:2017-04-20 Revised:2017-06-08 Online:2017-08-20 Published:2017-12-25
Supported by:
the National Natural Science Foundation of China(31471188, 81573241);China Postdoctoral Science Foundation(2016M602797);the Natural Science Basic Research Plan in Shaanxi Province of China(2016JQ3026)

摘要/Abstract

摘要：

为探索脑源性神经营养因子BDNF基因与中国汉族人群骨密度(bone mineral density, BMD)的关系,解析该基因调节骨质疏松症的功能机制,进而为中国汉族人群的骨质疏松症的防治提供依据。本研究从我国陕西地区征集了1300例汉族样本,并测量髋部/脊椎骨密度。选取BDNF基因上的14个标签SNPs进行基因分型,与1300例样本BMD进行关联分析,发现8个SNPs与髋部/脊椎BMD显著关联(P < 0.05)。其中,SNP rs16917237同时与髋部和脊椎骨密度关联,经Bonferroni校正后仍表现出显著性(0.05/14 = 0.0036)。整合连锁不平衡和单体型分析、表观功能注释、表达数量基因座分析、代谢通路分析进一步探索BDNF基因调节骨质疏松症的机制。构建小鼠前成骨细胞系(MC3T3-E1)人骨形态形成蛋白(rhBMP-2)诱导分化模型,利用小干扰RNA(siRNA)敲除BDNF。结果显示：14个SNPs位于同一单体型内;rs16917237在成骨细胞中表现出较强的激活型组蛋白H3K4me1、H3K4me3、H3K27ac修饰信号以及P300结合信号,表明其在成骨细胞中可能具有调控活性;rs16917237在11个组织中均能显著影响BDNF基因表达;BDNF基因位于与成骨细胞增殖分化相关的MAPK通路中;BDNF敲低能够显著降低MAPK通路中与成骨分化相关的CREB 基因mRNA和蛋白表达水平,提示其可能通过调控CREB表达进而影响成骨分化。生物信息学分析和功能实验结果一致,表明BDNF基因可能是影响骨质疏松症的重要功能基因。

关键词: 关联分析, 单核苷酸多态性, BDNF, 骨质疏松症, 成骨分化

Abstract:

To explore the relationship between brain-derived neurotrophic factor (BDNF) gene and bone mineral density (BMD) in Chinese Han population, we performed association analysis of 14 tag SNPs on BDNF gene with hip/spine BMD in 1300 Han Chinese samples from Shaanxi Province. We found that 8 of the 14 SNPs were significantly associated with hip or spine BMD (P < 0.05). Moreover, the SNP rs16917237 was significantly associated with both hip and spine BMD, with significant Bonferroni correlation (P value 0.05/14 = 0.0036) in hip BMD. To further explore the regulatory mechanism of BDNF gene in osteoporosis, we further performed a set of data analyses, including linkage disequilibrium and haplotype analysis, epigenetic annotation, expression quantitative trait locus (eQTL) analysis and metabolic pathway analysis. Further, we have established a mouse pre-osteoblasts differentiation cell model (MC3T3-E1) by recombination human bone morphogenetic protein (rh-BMP2) induction. siRNA- mediated knock down of BDNF in this cell model showed that all 14 SNPs are in the same haplotype block. Strong signals of active histone H3K4me1, H3K4me3, H3K27ac modifications and P300 binding were observed in osteoblasts, in the region surrounding the most significant SNP rs16917237, suggesting that this SNP might have a regulatory function in osteoblasts. Furthermore, analysis of genotype data of rs16917237 and BDNF expression in multiple tissues from GTEx showed that rs16917237 SNP could significantly affect the expression of BDNF in 11 tissues. Through analysis of the various BDNF pathways, we showed that BDNF participates in the MAPK pathway, which is a vital and well-established pathway affecting osteoblasts proliferation and differentiation. siRNA knock down of BDNF significantly decreased the mRNA and protein levels of CREB, which is important in the MAPK pathway in osteoblast differentiation. These findings suggest that BDNF might affect osteoblast differentiation via regulation of CREB expression. In conclusion, our results from combined genetic association and functional analyses show that BDNF is a vital osteoporosis susceptibility gene, which can affect BMD not only in Chinese Han but also likely in other populations.

Key words: candidate gene association analysis, SNP, BDNF, osteoporosis, osteoblast differentiation

杨熳,卢冰婕,段媛媛,陈晓峰,马建岗,郭燕. 骨质疏松症易感基因BDNF的遗传学关联分析及功能研究[J]. 遗传, 2017, 39(8): 726-736.

Yang Man,Lu Bingjie,Duan Yuanyuan,Chen Xiaofeng,Ma Jiangang,Guo Yan. Genetics association study and functional analysis on osteoporosis susceptibility gene BDNF[J]. Hereditas(Beijing), 2017, 39(8): 726-736.

参考文献

[1]	Harvey NC, Matthews P, Collins R, Cooper C.Osteoporosis epidemiology in UK Biobank: a unique opportunity for international researchers. Osteoporos Int, 2013, 24(12): 2903-2905.
[2]	Hadjidakis DJ, Androulakis II.Bone remodeling. Ann N Y Acad Sci, 2006, 1092: 385-396.
[3]	Peacock M, Turner CH, Econs MJ, Foroud T.Genetics of osteoporosis. Endocr Rev, 2002, 23(3): 303-326.
[4]	Ralston SH, Uitterlinden AG.Genetics of osteoporosis. Endocr Rev, 2010, 31(5): 629-662.
[5]	Livshits G, Deng HW, Nguyen TV, Yakovenko K, Recker RR, Eisman JA.Genetics of bone mineral density: evidence for a major pleiotropic effect from an intercontinental study. J Bone Miner Res, 2004, 19(6): 914-923.
[6]	Xu XH, Dong SS, Guo Y, Yang TL, Lei SF, Papasian CJ, Zhao M, Deng HW.Molecular genetic studies of gene identification for osteoporosis: the 2009 update. Endocr Rev, 2010, 31(4): 447-505.
[7]	Cardon LR, Bell JI.Association study designs for complex diseases. Nat Rev Genet, 2001, 2(2): 91-99.
[8]	Geng L, Yao ZW, Luo JY, Han LL, Lu Q.Association between Val80 polymorphism of the CYP19 gene, A163G polymorphism of the OPG gene and bone mineral density in post-menopausal Chinese women. Hereditas (Beijing), 2007, 29(11): 1345-1350.
[8]	耿力,姚珍薇, 骆建云, 韩力力, 卢起. 绝经后女性CYP19基因Val80及OPG基因A163G多态性与骨密度的相关性研究. 遗传, 2007, 29(11): 1345-1350.
[9]	Ghosh A, Carnahan J, Greenberg ME.Requirement for BDNF in activity-dependent survival of cortical neurons. Science, 1994, 263(5153): 1618-1623.
[10]	Acheson A, Conover JC, Fandl JP, DeChiara TM, Russell M, Thadani A, Squinto SP, Yancopoulos GD, Lindsay RM. A BDNF autocrine loop in adult sensory neurons prevents cell death. Nature, 1995, 374(6521): 450-453.
[11]	Yang YX, Chen YX, Wu JQ, Xing YM, Zeng FM, Huang YH, Cheng BW.Association study of CNR1, GAD1 and BDNF polymorphisms with male heroin dependence in the Dai population in Yunnan. Hereditas (Beijing), 2014, 36(9): 888-896.
[11]	杨亚鲜, 陈燕祥, 武健权, 邢豫明, 曾发明, 黄友浩, 程宝文. CNR1、GAD1和BDNF基因多态性与云南傣族男性海洛因依赖的相关性分析. 遗传, 2014, 36(9): 888-896.
[12]	Yang ZF, Ho DW, Lam CT, Luk JM, Lum CT, Yu WC, Poon RT, Fan ST.Identification of brain-derived neurotrophic factor as a novel functional protein in hepatocellular carcinoma. Cancer Res, 2005, 65(1): 219-225.
[13]	Hu Y, Wang YD, Guo T, Wei WN, Sun CY, Zhang L, Huang J.Identification of brain-derived neurotrophic factor as a novel angiogenic protein in multiple myeloma. Cancer Genet Cytogenet, 2007, 178(1): 1-10.
[14]	Sun CY, Chu ZB, She XM, Zhang L, Chen L, Ai LS, Hu Y.Brain-derived neurotrophic factor is a potential osteoclast stimulating factor in multiple myeloma. Int J Cancer, 2012, 130(4): 827-836.
[15]	Hutchison MR.BDNF alters ERK/p38 MAPK activity ratios to promote differentiation in growth plate chondrocytes. Mol Endocrinol, 2012, 26(8): 1406-1416.
[16]	Yamashiro T, Fukunaga T, Yamashita K, Kobashi N, Takano-Yamamoto T.Gene and protein expression of brain- derived neurotrophic factor and TrkB in bone and cartilage. Bone, 2001, 28(4): 404-409.
[17]	Deng FY, Tan LJ, Shen H, Liu YJ, Liu YZ, Li J, Zhu XZ, Chen XD, Tian Q, Zhao M, Deng HW.SNP rs6265 regulates protein phosphorylation and osteoblast differentiation and influences BMD in humans. J Bone Miner Res, 2013, 28(12): 2498-2507.
[18]	Deng HW, Shen H, Xu FH, Deng HY, Conway T, Zhang HT, Recker RR.Tests of linkage and/or association of genes for vitamin D receptor, osteocalcin, and parathyroid hormone with bone mineral density. J Bone Miner Res, 2002, 17(4): 678-686.
[19]	Guo Y, Tan LJ, Lei SF, Yang TL, Chen XD, Zhang F, Chen Y, Pan F, Yan H, Liu XG, Tian Q, Zhang ZX, Zhou Q, Qiu C, Dong SS, Xu XH, Guo YF, Zhu XZ, Liu SL, Wang XL, Li X, Luo Y, Zhang LS, Li M, Wang JT, Wen T, Drees B, Hamilton J, Papasian CJ, Recker RR, Song XP, Cheng J, Deng HW.Genome-wide association study identifies ALDH7A1 as a novel susceptibility gene for osteoporosis. PLoS Genet, 2010, 6(1): e1000806.
[20]	Zhu DL, Guo Y, Zhang Y, Dong SS, Xu W, Hao RH, Chen XF, Yan H, Yang SY, Yang TL.A functional SNP regulated by miR-196a-3p in the 3'UTR of FGF2 is associated with bone mineral density in the Chinese population. Hum Mutat, 2017, 38(6): 725-735.
[21]	Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, De Bakker PIW, Daly MJ, Sham PC.PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet, 2007, 81(3): 559-575.
[22]	Chen GQ, Deng CX, Li YP.TGF-β and BMP signaling in osteoblast differentiation and bone formation. Int J Biol Sci, 2012, 8(2): 272-288.
[23]	Zhang RR, Edwards JR, Ko SY, Dong SS, Liu HB, Oyajobi BO, Papasian C, Deng HW, Zhao M.Transcriptional regulation of BMP2 expression by the PTH-CREB signaling pathway in osteoblasts. PLoS One, 2011, 6(6): e20780.
[24]	Barroso I, Luan JA, Middelberg RPS, Harding AH, Franks PW, Jakes RW, Clayton D, Schafer AJ, O'Rahilly S, Wareham NJ. Candidate gene association study in type 2 diabetes indicates a role for genes involved in β-cell function as well as insulin action. PLoS Biol, 2003, 1(1): E20.
[25]	Tabor HK, Risch NJ, Myers RM.Candidate-gene approaches for studying complex genetic traits: practical considerations. Nat Rev Genet, 2002, 3(5): 391-397.
[26]	Estrada K, Styrkarsdottir U, Evangelou E, Hsu YH, Duncan EL, Ntzani EE, Oei L, Albagha OM, Amin N, Kemp JP, Koller DL, Li G, Liu CT, Minster RL, Moayyeri A, Vandenput L, Willner D, Xiao SM, Yerges-Armstrong LM, Zheng HF, Alonso N, Eriksson J, Kammerer CM, Kaptoge SK, Leo PJ, Thorleifsson G, Wilson SG, Wilson JF, Aalto V, Alen M, Aragaki AK, Aspelund T, Center JR, Dailiana Z, Duggan DJ, Garcia M, Garcia-Giralt N, Giroux S, Hallmans G, Hocking LJ, Husted LB, Jameson KA, Khusainova R, Kim GS, Kooperberg C, Koromila T, Kruk M, Laaksonen M, Lacroix AZ, Lee SH, Leung PC, Lewis JR, Masi L, Mencej-Bedrac S, Nguyen TV, Nogues X, Patel MS, Prezelj J, Rose LM, Scollen S, Siggeirsdottir K, Smith AV, Svensson O, Trompet S, Trummer O, Van Schoor NM, Woo J, Zhu K, Balcells S, Brandi ML, Buckley BM, Cheng S, Christiansen C, Cooper C, Dedoussis G, Ford I, Frost M, Goltzman D, González-Macías J, K?h?nen M, Karlsson M, Khusnutdinova E, Koh JM, Kollia P, Langdahl BL, Leslie WD, Lips P, Ljunggren ?, Lorenc RS, Marc J, Mellstr?m D, Obermayer-Pietsch B, Olmos JM, Pettersson-Kymmer U, Reid DM, Riancho JA, Ridker PM, Rousseau F, Slagboom PE, Tang NL, Urreizti R, Van Hul W, Viikari J, Zarrabeitia MT, Aulchenko YS, Castano-Betancourt M, Grundberg E, Herrera L, Ingvarsson T, Johannsdottir H, Kwan T, Li R, Luben R, Medina-Gómez C, Palsson ST, Reppe S, Rotter JI, Sigurdsson G, Van Meurs JB, Verlaan D, Williams FM, Wood AR, Zhou Y, Gautvik KM, Pastinen T, Raychaudhuri S, Cauley JA, Chasman DI, Clark GR, Cummings SR, Danoy P, Dennison EM, Eastell R, Eisman JA, Gudnason V, Hofman A, Jackson RD, Jones G, Jukema JW, Khaw KT, Lehtim?ki T, Liu Y, Lorentzon M, McCloskey E, Mitchell BD, Nandakumar K, Nicholson GC, Oostra BA, Peacock M, Pols HA, Prince RL, Raitakari O, Reid IR, Robbins J, Sambrook PN, Sham PC, Shuldiner AR, Tylavsky FA, Van Duijn CM, Wareham NJ, Cupples LA, Econs MJ, Evans DM, Harris TB, Kung AW, Psaty BM, Reeve J, Spector TD, Streeten EA, Zillikens MC, Thorsteinsdottir U, Ohlsson C, Karasik D, Richards JB, Brown MA, Stefansson K, Uitterlinden AG, Ralston SH, Ioannidis JP, Kiel DP, Rivadeneira F. Genome-wide meta-analysis identifyies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture. Nat Genet, 2012, 44(5): 491-501.
[27]	Guo Y, Dong SS, Chen XF, Jing YA, Yang M, Yan H, Shen H, Chen XD, Tan LJ, Tian Q, Deng HW, Yang TL.Integrating epigenomic elements and GWASs identifies BDNF gene affecting bone mineral density and osteoporotic fracture risk. Sci Rep, 2016, 6: 30558.
[28]	Asaumi K, Nakanishi T, Asahara H, Inoue H, Takigawa M.Expression of neurotrophins and their receptors (TRK) during fracture healing. Bone, 2000, 26(6): 625-633.
[29]	Ge J, Zhang Y, Zheng ZZ, Huang YW, Song HS.Research progress on CREB and the signal transduction pathways of its phosphorylation. J Anhui Agri Sci, 2010, 38(30): 16769-16771, 16774.
[29]	葛军, 张玉, 郑增长, 黄延旺, 宋红生. CREB转录因子及其磷酸化信号通路的研究进展. 安徽农业科学, 2010, 38(30): 16769-16771, 16774.
[30]	Sato K, Suematsu A, Nakashima T, Takemoto-Kimura S, Aoki K, Morishita Y, Asahara H, Ohya K, Yamaguchi A, Takai T, Kodama T, Chatila TA, Bito H, Takayanagi H.Regulation of osteoclast differentiation and function by the CaMK-CREB pathway. Nat Med, 2006, 12(12): 1410-1416.
[31]	Hutchison MR.Mice with a conditional deletion of the neurotrophin receptor TrkB are dwarfed, and are similar to mice with a MAPK14 deletion, PLoS One, 2013, 8(6): e66206.
[32]	Mikami Y, Suzuki S, Ishii Y, Watanabe N, Takahashi T, Isokawa K, Honda MJ.The p75 neurotrophin receptor regulates MC3T3-E1 osteoblastic differentiation. Differentiation, 2012, 84(5): 392-399.

编辑推荐

Metrics

www.chinagene.cn
备案号：京ICP备09063187号-4
总访问:,今日访问:,当前在线:

骨质疏松症易感基因BDNF的遗传学关联分析及功能研究

Genetics association study and functional analysis on osteoporosis susceptibility gene BDNF

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

[1]	梁文权,侯豫,赵存友. 精神分裂症相关单核苷酸多态性调控microRNA功能研究进展[J]. 遗传, 2019, 41(8): 677-685.
[2]	张晴,平杰,张昊翔,康波,李元丰,周钢桥. MKL1基因多态性与高原环境适应性的遗传关联研究[J]. 遗传, 2019, 41(7): 634-643.
[3]	文超良, 孙从佼, 杨宁. 从遗传力到肠菌力：概念及研究进展[J]. 遗传, 2019, 41(11): 1023-1040.
[4]	彭哲也,唐紫珺,谢民主. 机器学习方法在基因交互作用探测中的研究进展[J]. 遗传, 2018, 40(3): 218-226.
[5]	刘莉莉, 郭爱伟, 吴培福, 陈粉粉, 杨亚晋, 张勤. 敲降VPS28基因对中国荷斯坦奶牛乳脂合成的调控[J]. 遗传, 2018, 40(12): 1092-1100.
[6]	张统雨,朱才业,杜立新,赵福平. 羊重要性状全基因组关联分析研究进展[J]. 遗传, 2017, 39(6): 491-500.
[7]	胡鹏飞, 徐佳萍, 艾成, 邵秀娟, 王洪亮, 董依萌, 崔学哲, 杨福合, 邢秀梅. 利用全基因组重测序分析鹿茸重量相关基因[J]. 遗传, 2017, 39(11): 1090-1101.
[8]	邹文超, 沈林林, 沈建国, 蔡伟, 詹家绥, 高芳銮. 马铃薯Y病毒多基因系统发育分析及其在株系鉴定中的应用[J]. 遗传, 2017, 39(10): 918-929.
[9]	弓弦,张超,伊利亚斯·艾萨,时瑛,杨雪唯,努尔斯曼古丽奥斯曼,关亚群,徐书华. 2型糖尿病易感候选基因在世界不同人群中的多样性比较分析[J]. 遗传, 2016, 38(6): 543-559.
[10]	张涛, 王文浩, 张跟喜, 王金玉, 薛倩, 顾玉萍. 京海黄鸡体重性状全基因组关联分析[J]. 遗传, 2015, 37(8): 811-820.
[11]	陈开旭, 王为兰, 张富春, 郑秀芬. 人类身高的遗传学研究进展[J]. 遗传, 2015, 37(8): 741-755.
[12]	宋庆峰, 张红星, 马亦龙, 周钢桥. 复杂疾病的遗传易感基因区域的精细定位[J]. 遗传, 2014, 36(1): 2-10.
[13]	罗旭红刘志芳董长征. 基因水平的关联分析方法[J]. 遗传, 2013, 35(9): 1065-1071.
[14]	秦龙娟丁达霞崔璐璐黄青阳. SOST基因的表达调控[J]. 遗传, 2013, 35(8): 939-947.
[15]	李骞刘舒媛林克勤孙浩于亮黄小琴褚嘉祐杨昭庆. EGLN1基因6个单核苷酸多态性与高海拔低氧适应的相关性[J]. 遗传, 2013, 35(8): 992-998.