ADAMTS14基因单核苷酸多态性与汉族女性颞下颌关节骨关节炎的相关性研究
王丹丹1, 甘业华2, 马绪臣2, 孟娟红1,
1. 北京大学口腔医学院·口腔医院,口腔颌面外科,北京 100081
2. 北京大学口腔医学院·口腔医院,颞下颌关节病及口颌面痛诊治中心 口腔数字化医疗技术和材料国家工程实验室 口腔数字医学北京市重点实验室, 北京 100081
摘要
目的:探讨 ADAMTS14基因rs4747096位点单核苷酸多态性与中国汉族女性颞下颌关节骨关节炎之间的相关性。方法:采用病例-对照研究,纳入213例受试者,均为中国汉族女性,其中病例组103例,健康对照组110例。从采集的外周静脉血中提取基因组DNA,通过聚合酶链式反应扩增目的片段,对反应产物测序确定基因分型并进行统计分析。病例组及实验组的年龄分布采用 t检验,基因型及等位基因频率分布采用卡方检验。结果:病例组103例中基因型频率为38.8%(AA)、55.4%(AG)、5.8%(GG),对照组110例中基因型频率为40.9%(AA)、43.6%(AG)、15.5%(GG)。基因型频率在两组中的差异有统计学意义( P=0.047),等位基因频率在两组中差异无统计学意义( P>0.05)。显性遗传模式中,病例组基因型频率为94.2%(AA或AG)、5.8%(GG),对照组基因型频率为84.5%(AA或AG)、15.5%(GG),病例组中AA、AG型明显高于对照组,是发病的风险因素( OR=1.114,95% CI:1.015~1.223, P=0.028)。结论: ADAMTS14基因rs4747096位点的单核苷酸多态性与中国汉族女性人群颞下颌关节骨关节炎之间存在明显相关性,rs4747096位点的单核苷酸多态性在颞下颌关节骨关节炎患者及健康人群中的分布存在差异。
关键词: 颞下颌关节; 骨关节炎; ADAMTS14; 单核苷酸多态性; 女性
中图分类号:R782.6 文献标志码:A 文章编号:1671-167X(2018)02-0279-05
Association between ADAMTS14 gene polymorphism and the temporomandibular joint osteoarthritis in Chinese Han females
WANG Dan-dan1, GAN Ye-hua2, MA Xu-chen2, MENG Juan-hong1,
1. Department of Oral and Maxillofacial Surgery, Beijing 100081, China
2. Center for Temporomandibular Joint Disorder and Orofacial Pain, Peking University School and Hospital of Stomatology & National Engineering Laboratory for Digital and Material Technology of Stomatology & Beijing Key Laboratory of Digital Stomatology, Beijing 100081, China
△ Corresponding author’s e-mail, jhmeng@263.net
Abstract
Objective:To investigate the association between single nucleotide polymorphisms (SNP) of ADAMTS14 gene rs4747096 and osteoarthritis of the temporomandibular joint in Chinese Han females.Methods:As a case-control study, a total of 213 Chinese Han females were involved in the present stu-dy, which contained 103 temporomandibular joint osteoarthritis patients and 110 healthy people who had no symptoms or signs of temporomandibular joint osteoarthritis as control. Peripheral blood samples were collected from each participant. Genomic DNAs of temporomandibular joint osteoarthritis patients and healthy control were extracted from peripheral venous blood, which were stored in -80 ℃ refrigerator by using DNA extraction kits. The designed primers were used for polymerase chain reaction (PCR) amplification of specific DNA fragments. Genotype was determined by sequencing the PCR products. The software Chromas 2.22 was used to analyze the genotype. The genotype distributions, allele frequencies and genetic models between the patients and controls were compared. The age distribution was checked by t-test. Genotype and allele frequency were detected by Chi-square test.Results:In the present study, there were no significant differences between the osteoarthritis patients and healthy controls in terms of age. The genotype distribution was in accordance with Hardy-Weinberg equilibrium in the two groups. The genotype frequency of the ADAMTS14 (rs4747096) in the experimental group was 38.8% (AA), 55.4% (AG), and 5.8% (GG), respectively. The genotype frequency in the control group was 40.9% (AA), 43.6% (AG), and 15.5% (GG), respectively. The difference of genotype frequency of the ADAMTS14 (rs4747096) was significant between the experimental group and the control group ( P=0.047). There was no significant difference in allele frequency between the two groups ( P=0.415). AA and AG genotypes significantly increased the risk of the disease compared with GG in dominant model ( OR=1.114, 95% CI: 1.015-1.223, P=0.028).Conclusion:A significant correlationship was found between the ADAMTS14 (rs4747096) SNP and the temporomandibular joint osteo-arthritis in Chinese Han females. The distribution of rs4747096 may be different between temporomandi-bular joint osteoarthritis and healthy population.
Key words: Temporomandibular joint; Osteoarthritis; ADAMTS14; Single nucleotide polymorphism; Female

颞下颌关节骨关节炎(temporomandibular joint osteoarthritis, TMJOA)是以关节软骨的进行性破坏为特征性表现的关节退行性疾病[1]。确切病因及发病机制尚不明确, 目前普遍认为环境因素与多基因协同作用使TMJOA易感性增加[2]。单核苷酸多态性(single nucleotide polymorphism, SNP)约占人类基因多态性的90%, 是人类DNA变异中最常见、最稳定的变异, 目前广泛用于基因多态性与疾病的关联性研究[3]。位于编码区内的SNP(coding SNP, cSNP)包括同义cSNP(synonymous cSNP)和非同义cSNP(non-synonymous cSNP)。

负荷过重、蛋白酶及细胞因子表达异常等因素可使软骨的细胞外基质(extracellular matrix, ECM)分解和合成代谢失衡, 最终导致软骨破坏[4]。带有血小板凝血酶敏感蛋白结构域的解聚素与金属蛋白酶(a disintegrin and metalloproteinase with thrombospondin motifs, ADAMTS)在软骨ECM降解中发挥至关重要的作用[5, 6]。Colige等[7]的研究发现, ADAMTS14参与Ⅱ 型胶原纤维的形成。Ⅱ 型胶原是构成关节软骨ECM的主要成分之一, 其代谢异常可导致关节软骨机械强度减低, 是骨关节炎(osteoarthritis, OA)发病过程中的重要因素[8]

ADAMTS14基因rs4747096是位于编码区的非同义SNP, 可使翻译的氨基酸序列发生改变, 最终可能影响生物学性状。Rodriguez-Lopez等[9]的研究表明, 需行膝关节置换的女性OA患者中, ADAMTS14基因rs4747096非同义cSNP等位基因G频率明显高于对照组, ADAMTS14基因多态性与膝关节OA的发病有相关性。然而, ADAMTS14基因多态性与TMJOA的相关性目前尚未见研究报道。本研究旨在通过对汉族女性TMJOA病例组及健康对照组ADAMTS14基因rs4747096非同义cSNP进行检测, 探讨ADAMTS14基因多态性与汉族女性TMJOA之间的相关性。

1 资料与方法
1.1 研究对象

本研究的纳入对象均为中国汉族女性, 采用病例-对照研究。

病例组为在北京大学口腔医院颞下颌关节病及口颌面痛诊治中心就诊, 根据RDC/TMD (Research Diagnostic Criteria for Temporomandibular Disorders)诊断为TMJOA, 且影像学颞下颌关节有髁突骨质破坏、增生硬化、骨赘形成或囊样变的患者, 排除颞下颌关节感染性、创伤性、类风湿性等特殊类型的关节炎, 以及伴有髁突畸形、肿瘤的患者, 最终共纳入103例。病例组患者平均年龄为(33.74± 13.41)岁, 单侧发病44例(42.7%), 双侧发病59例(57.3%)。

对照组为无TMJOA及其他部位关节炎表现者, 无颞下颌关节紊乱病病史, 且影像学表现双侧髁突骨皮质光滑、连续, 共纳入110例。对照组平均年龄为(30.88± 9.57)岁。

本研究经北京大学口腔医院生物医学伦理委员会批准, 受试者均充分告知并签署知情同意书。

1.2 基因组DNA提取

病例组与对照组均抽取2 mL外周静脉血, 置于EDTA抗凝管中, -80 ℃保存。使用血液基因组DNA提取试剂盒(QIAamp DNA Blood Mini Kits, QIAGEN, 德国)提取基因组DNA, 具体方法根据产品操作手册进行。采用分光光度计(721分光光度计, 上海第三分析仪器厂)测光密度值(D)和DNA浓度, -20 ℃保存。

1.3 PCR

采用Generunner 6.2.07 Beta软件进行引物设计, 由北京六合华大基因公司进行引物合成。正向引物序列:5'-TGTGCAGGACCAACGCCAACAG-3', 反向引物序列:5'-GGAATTGCAGGTAACGGCTCATG-3'。PCR反应体系内含2× Es Taq MasterMix(康为世纪)12 μ L, 10 μ mol/L正、反向引物各2 μ L, DNA模版2 μ L, 灭菌双蒸水9 μ L。采用PCR仪(Biometra, 华粤)进行扩增, 反应条件:94 ℃预变性2 min后进入30个循环, 每个循环包括98 ℃变性10 s, 61 ℃退火30 s, 72 ℃延伸30 s, 最终72 ℃延伸2 min, 4 ℃保存。

1.4 测序分析

将PCR产物送往北京六合华大基因公司进行DNA测序, 采用Chromas 2.22软件分析测序峰图, 确定基因分型。

1.5 统计学分析

采用SPSS 21.0软件进行统计学分析。年龄分布采用t检验, 基因型及等位基因频率分布采用卡方检验, 计算比值比(odds ratio, OR) 及95%可信区间(confidence interval, CI), P< 0.05为差异有统计学意义。

2 结果

病例组平均年龄(33.74± 13.41)岁, 对照组平均年龄(30.88± 9.57)岁, 病例组与对照组平均年龄差异无统计学意义(t=1.798, P=0.074)。

rs4747096位点AA、AG、GG基因型测序峰图见图1。基因型和等位基因频率满足Hardy-Weinberg平衡定律。ADAMTS14基因rs4747096位点SNP基因型频率在病例组和对照组中差异有统计学意义(χ 2=6.103, P=0.047, 表1), 等位基因型频率在两组中差异无统计学意义(χ 2=0.663, P=0.415, 表2)。显性模型中AA、AG为患病风险因素, 相对于GG明显增加了患TMJOA的风险(OR=1.114, 95% CI: 1.015~1.223, P=0.028, 表2)。

图1 rs4747096位点DNA测序峰图(AA、AG、GG)Figure 1 DNA sequencing chromatograms of rs474709 (AA, AG, GG)

表1 ADAMTS14基因rs4747096非同义编码区SNP基因型频率分布 Table 1 Distribution of genotype frequency of non-synonymous coding SNP for rs4747096 of ADAMTS14 gene
表2 rs4747096等位基因频率, 显性模型、隐形模型基因型频率分布 Table 2 Genotype frequency distribution of alleles, dominant model and recessive model for rs4747096
3 讨论

TMJOA可引起颞下颌关节疼痛、下颌运动受限以及咀嚼、言语等功能障碍, 长期慢性疼痛的症状甚至可引发心理问题。颞下颌关节髁突软骨破坏、骨质吸收严重者可导致错牙合、面部畸形, 影响患者的生活质量。不同个体对TMJOA的易感性存在差

[10], 目前普遍认为是多基因协同作用的结果[11], 但确切的分子机制仍不明确。

关节软骨的进行性破坏是TMJOA的特征性表现之一。颞下颌关节的关节软骨为纤维软骨, 软骨基质中的胶原蛋白占软骨干重的50%~60%, 以Ⅱ 型胶原(90%)为主[4]。Ⅱ 型胶原蛋白可抵抗软骨承受的张力, 其正常结构对维护软骨的机械强度至关重要。关节软骨浅表区的胶原纤维破坏是早期OA的组织病理学标志[12]。动物实验表明, 早期OA关节软骨中胶原纤维直径增大, 此类软骨微结构的改变比软骨下骨破坏出现得更早[13, 14]

ADAMTSs是一类具有金属蛋白酶结构域、解聚素结构域及血小板反应蛋白结构域的分泌型蛋白酶[15], 参与软骨基质中的胶原、蛋白多糖降解。针对ADAMTSs家族成员基因多态性、蛋白表达水平等的相关研究, 已成为OA发病机制研究的焦点。Li等[16]发现早期膝关节OA患者血清中ADAMTS4水平要明显高于健康人群, 而中、晚期患者血清ADAMTS5水平明显高于健康人群, 并推测ADAMTS4有可能成为早期OA的血清学标记物。Majumdar等[17]的研究发现敲除ADAMTS4/5基因的小鼠OA的严重程度较野生型明显减轻。Miller等[18]研究发现应用ADAMTS5抗体可减缓小鼠膝关节OA病程中的软骨破坏。Gok等[19]发现ADAMTS9基因CA重复序列长度≥ 20者膝关节OA更严重。ADAMTS1、4、5、8、9、15可促进关节软骨中的蛋白聚糖降解, ADAMTS7、12参与软骨基质中寡聚蛋白的降解[20]。国内外的研究成果显示ADAMTSs与OA密切相关, 而进一步明确这种相关性则有助于从根本上预防TMJOA的发生, 并有助于TMJOA的早期诊断和治疗。

ADAMTS14基因位于染色体10q22.1, 目前, ADAMTS14的功能尚不完全明确。由于ADAMTS14基因和ADAMTS3基因(63%)具有同源性, Fernandes等[21]推测ADAMTS14与ADAMTS3功能相似, 可裂解胶原单体N-端前肽, 促进Ⅱ 型胶原单体聚合成Ⅱ 型胶原纤维, 参与软骨损伤后的修复过程。

ADAMTS14基因编码的蛋白质包含1 223个氨基酸, 由信号肽、前结构域、金属蛋白酶结构域、解聚素结构域、4个凝血酶敏感蛋白结构域和羧基末端结构域组成, 4个凝血酶敏感蛋白结构域被富含半胱氨酸的氨基酸片段间隔开。rs4747096位点A> G, 等位基因G是祖先型等位基因, 但在欧、亚、非洲人中等位基因A所占比例更大, 因此等位基因A可能是自然选择过程中的优势基因。碱基序列的改变导致编码的蛋白质分子羧基端结构域中的氨基酸由谷氨酸(GAA)转变为甘氨酸(GGA), 但目前由此引发的功能改变尚不明确。我们推测ADAMTS14基因rs4747096位点的多态性, 可能会引起Ⅱ 型胶原蛋白成熟过程异常, 从而导致软骨破坏, 但这一过程仍需进一步的研究验证。

目前国际上已有多项关于ADAMTS14基因rs4747096非同义cSNP与OA相关性的研究。有一项涉及3 000例英国女性的全基因组关联分析(genome-wide association study, GWAS)研究表明, ADAMTS14基因rs4747096非同义cSNP与膝关节OA有明显相关性。Poonpet等[22]的研究发现, 泰国女性膝关节OA患者的基因型AA及等位基因A明显高于对照组, AA和AG增加了患病风险, 这一结果与我们的研究结果相符。除与OA易感性有关外, El Khoury等[23]的研究发现在跟腱病患者中, ADAMTS14基因rs4747096位点GG基因型可能是一个保护性因素, 可延缓病程的进展。

本研究通过对中国汉族女性病例组与对照组的ADAMTS14基因rs4747096非同义cSNP的基因型检测分析发现, 两组中ADAMTS14基因rs4747096非同义cSNP基因型频率的差异有统计学意义。显性模型中, 病例组的AA和AG明显高于对照组, 是TMJOA的风险因素(OR=1.114, 95% CI: 1.015~1.223, P=0.028), GG则是TMJOA的保护性因素(OR=0.377, 95% CI: 0.155~0.919, P=0.028)。

通过分析本研究结果以及回顾以往文献, 我们发现ADAMTS14基因rs4747096 非同义cSNP与高加索人种、蒙古人种的女性OA都有不同程度的相关性, 但在男性OA患者中未发现明显相关性[22], 我们推测这可能与女性特殊的雌激素水平有关, 今后可进一步对ADAMTS14、雌激素及其受体的基因多态性在OA发病中是否具有协同作用进行分析。

现有的TMJOA治疗方法的主要作用为缓解症状, 恢复功能, 防止软骨及软骨下骨进一步破坏, 较难从根本上阻止或逆转已发生的软骨病损[24]。然而, 进一步明确ADAMTS14基因rs4747096非同义cSNP与TMJOA的相关性, 明确其在发病过程中的作用, 有助于了解OA发病的分子机制, 探索新的药物靶点[2]。本研究样本量较小且仅限于女性, 今后还需要进一步扩大样本量进行验证, 并对更多的ADAMTS14基因SNP位点进行检测, 同时也应增进ADAMTS14基因的功能研究。

The authors have declared that no competing interests exist.

参考文献
[1] 孟娟红, 甘业华, 马绪臣. 颞下颌关节骨关节炎发病的分子机制及相关治疗的实验研究[J]. 北京大学学报(医学版), 2013, 45(1): 5-8. [本文引用:1]
[2] Meng J, Ma X, Ma D, et al. Microarray analysis of differential gene expression in temporomand ibular joint condylar cartilage after experimentally induced osteoarthritis[J]. Osteoarthritis Cartilage, 2005, 13(12): 1115-1125. [本文引用:2]
[3] Collins FS, Brooks LD, Chakravarti A. A DNA polymorphism discovery resource for research on human genetic variation[J]. Genome Res, 1998, 8(12): 1229-1231. [本文引用:1]
[4] Martel-Pelletier J, Boileau C, Pelletier J P, et al. Cartilage in normal and osteoarthritis conditions[J]. Best Pract Res Clin Rheumatol, 2008, 22(2): 351-384. [本文引用:2]
[5] Verma P, Dalal K. ADAMTS-4 and ADAMTS-5: key enzymes in osteoarthritis[J]. J Cell Biochem, 2011, 112(12): 3507-3514. [本文引用:1]
[6] Apte SS. A disintegrin-like and metalloprotease (reprolysin-type) with thrombospondin type 1 motif (ADAMTS) superfamily: functions and mechanisms[J]. J Biol Chem, 2009, 284(46): 31493-31497. [本文引用:1]
[7] Colige A, Vand enberghe I, Thiry M, et al. Cloning and characte-rization of ADAMTS-14, a novel ADAMTS displaying high homology with ADAMTS-2 and ADAMTS-3[J]. J Biol Chem, 2002, 277(8): 5756-5766. [本文引用:1]
[8] Nagase H, Kashiwagi M. Aggrecanases and cartilage matrix degradation[J]. Arthritis Res Ther, 2003, 5(2): 94-103. [本文引用:1]
[9] Rodriguez-Lopez J, Pombo-Suarez M, Loughlin J, et al. Association of a nsSNP in ADAMTS14 to some osteoarthritis phenotypes[J]. Osteoarthritis Cartilage, 2009, 17(3): 321-327. [本文引用:1]
[10] Tanaka E, Detamore MS, Mercuri LG. Degenerative disorders of the temporomand ibular joint: etiology, diagnosis, and treatment[J]. J Dent Res, 2008, 87(4): 296-307. [本文引用:1]
[11] Xiao JL, Meng JH, Gan YH, et al. Association of GDF5, SMAD3 and RUNX2 polymorphisms with temporomand ibular joint osteoarthritis in female Han Chinese[J]. J Oral Rehabil, 2015, 42(7): 529-536. [本文引用:1]
[12] Poole AR, Kobayashi M, Yasuda T, et al. Type II collagen degradation and its regulation in articular cartilage in osteoarthritis [J]. Ann Rheum Dis, 2002, 61(Suppl 2): ii78-81. [本文引用:1]
[13] Silbermann M, Livne E. Age-related degenerative changes in the mouse mand ibular joint[J]. J Anat, 1979, 129(Pt 3): 507-520. [本文引用:1]
[14] Chen J, Gupta T, Barasz JA, et al. Analysis of microarchitectural changes in a mouse temporomand ibular joint osteoarthritis model[J]. Arch Oral Biol, 2009, 54(12): 1091-1098. [本文引用:1]
[15] Bolz H, Ramirez A, von Brederlow B, et al. Characterization of ADAMTS14, a novel member of the ADAMTS metalloproteinase family[J]. Biochim Biophys Acta, 2001, 1522(3): 221-225. [本文引用:1]
[16] Li W, Du C, Wang H, et al. Increased serum ADAMTS-4 in knee osteoarthritis: a potential indicator for the diagnosis of osteoarthritis in early stages[J]. Genet Mol Res, 2014, 13(4): 9642-9649. [本文引用:1]
[17] Majumdar MK, Askew R, Schelling S, et al. Double-knockout of ADAMTS-4 and ADAMTS-5 in mice results in physiologically normal animals and prevents the progression of osteoarthritis[J]. Arthritis Rheum, 2007, 56(11): 3670-3674. [本文引用:1]
[18] Miller RE, Tran PB, Ishihara S, et al. Therapeutic effects of an anti-ADAMTS-5 antibody on joint damage and mechanical allody-nia in a murine model of osteoarthritis[J]. Osteoarthritis Cartilage, 2016, 24(2): 299-306. [本文引用:1]
[19] Gok K, Cemeroglu O, Cakirbay H, et al. Relationship between cytosine-adenine repeat polymorphism of ADAMTS9 gene and cli-nical and radiologic severity of knee osteoarthritis[J]. Int J Rheum Dis, 2016, 19(5): 498-505. [本文引用:1]
[20] Yang CY, Chanalaris A, Troeberg L. ADAMTS and ADAM metalloproteinases in osteoarthritis—looking beyond the “usual suspects”[J]. Osteoarthritis Cartilage, 2017, 25(7): 1000-1009. [本文引用:1]
[21] Fernand es RJ, Hirohata S, Engle JM, et al. Procollagen II amino propeptide processing by ADAMTS-3. Insights on dermatosparaxis[J]. J Biol Chem, 2001, 276(34): 31502-31509. [本文引用:1]
[22] Poonpet T, Honsawek S, Tammachote N, et al. ADAMTS14 gene polymorphism associated with knee osteoarthritis in Thai women[J]. Genet Mol Res, 2013, 12(4): 5301-5309. [本文引用:2]
[23] El Khoury L, Posthumus M, Collins M, et al. Polymorphic variation within the ADAMTS2, ADAMTS14, ADAMTS5, ADAM12 and TIMP2 genes and the risk of Achilles tendon pathology: a genetic association study[J]. J Sci Med Sport, 2013, 16(6): 493-498. [本文引用:1]
[24] Wang XD, Zhang JN, Gan YH, et al. Current understand ing of pathogenesis and treatment of TMJ osteoarthritis[J]. J Dent Res, 2015, 94(5): 666-673. [本文引用:1]