孫玲玲，董博翰，凌烈鋒

(皖南醫(yī)學院 生物化學教研室，安徽 蕪湖　241002)

?

臨床經(jīng)驗交流

IL 12B基因多態(tài)性與Graves病相關(guān)性的Meta分析

孫玲玲，董博翰，凌烈鋒

(皖南醫(yī)學院 生物化學教研室，安徽 蕪湖241002)

[摘要]目的 Meta分析評價IL 12B基因3′UTR 1188A/C多態(tài)性與Graves病的相關(guān)性。方法 設(shè)定入選條件，檢索PubMed、EMbase、中國知網(wǎng)、重慶維普、萬方數(shù)據(jù)庫，收集關(guān)于IL12B基因3′UTR 1188A/C多態(tài)性與Graves病相關(guān)性的文獻，對符合入選條件的研究進行異質(zhì)性檢驗Meta分析，計算不同遺傳模型下1188A/C多態(tài)性與Graves病發(fā)病的相關(guān)性。結(jié)果 5篇文獻納入研究，包括747名病例和698名對照。合并分析結(jié)果顯示，IL12B基因3′UTR 1188A/C多態(tài)性與Graves病發(fā)病風險性無統(tǒng)計學意義(顯性遺傳模型：OR=1.18，95% CI=0.24～2.17，P=0.59；隱性遺傳模型：OR=1.61，95% CI=0.97～2.67，P=0.06；等位基因模型：OR=1.32，95% CI=0.93～1.87，P=0.12)。對3個日本人群的研究進行亞組分析，結(jié)果亦未發(fā)現(xiàn)1188A/C多態(tài)位點的變異與GD發(fā)病有統(tǒng)計學相關(guān)性(顯性遺傳模型：OR=0.78，95% CI=0.45～1.37，P=0.39；隱性遺傳模型：OR=1.15，95% CI=0.87～1.52，P=0.31；等位基因模型：OR=0.99，95% CI=0.83～1.18，P=0.92)。結(jié)論 IL12B基因1188A/C多態(tài)性可能與GD發(fā)病相關(guān)性無關(guān)，但由于研究樣本數(shù)較少，尚需大樣本研究加以證實。

[關(guān)鍵詞]Graves?。籌L12B；基因多態(tài)性；Meta分析

Graves病(Graves’ disease，GD)又稱毒性彌漫性甲狀腺腫，是一種器官特異性自身免疫性疾病，其發(fā)病機制依賴于遺傳因素和環(huán)境因素的共同作用，屬于多基因遺傳病。連鎖分析研究發(fā)現(xiàn)5q31-33區(qū)域存在著東亞人群Graves病的易感位點[1-3]，白介素(interleukin，IL)12B位于此區(qū)域。關(guān)于IL12B基因多態(tài)性與GD的研究目前已有多篇報道，但均為小樣本研究，統(tǒng)計效力較低。本研究收集國內(nèi)外關(guān)于IL12B基因3′非翻譯區(qū)(3′UTR)1188A/C基因多態(tài)性與GD相關(guān)性的研究，并進行Meta分析，客觀評價IL12B基因1188A/C多態(tài)性在GD發(fā)病中的作用。

1資料與方法

1.1文獻收集及納入與排除標準以“Graves’ disease”、“GD”、“IL 12B”、“IL12 p40”、“polymorphism”、“Graves病”、“白介素12B”、“多態(tài)”作為關(guān)鍵詞聯(lián)合檢索PubMed、Embase、Cochrane英文數(shù)據(jù)庫和中國知網(wǎng)、重慶維普、萬方等中文數(shù)據(jù)庫(含碩博士學位論文)，收集檢索入庫—2015年3月30日以前有關(guān)IL12B基因多態(tài)性與GD相關(guān)性的文獻。文獻納入標準：①原始文獻為公開發(fā)表關(guān)于IL12B 基因多態(tài)性與GD關(guān)聯(lián)的病例對照研究；②病例組均為經(jīng)醫(yī)院確診的GD患者，對照組為健康人群，排除甲狀腺疾病史及甲狀腺疾病家族史；③排除有類風濕性關(guān)節(jié)炎、哮喘、系統(tǒng)性硬化等自身免疫性疾?。虎芪墨I數(shù)據(jù)完整或根據(jù)結(jié)果能推算出分析所需數(shù)據(jù)；⑤納入文獻各組對照基因型均需符合Hardy-Weinberg遺傳平衡定律。去除重復(fù)報道、無對照文獻和數(shù)據(jù)描述不清無法利用的文獻。

1.2資料提取與分析由2名研究者按照納入和排除標準獨立進行文獻篩選、資料提取、文獻質(zhì)量評價并交叉核對，如遇有爭議的文獻則由本小組成員共同討論決定。分析各個研究發(fā)表信息，摘錄研究的發(fā)表年代、第一作者姓名、研究人群、多態(tài)性檢測方法、基因型和等位基因頻率等文獻資料。

1.3統(tǒng)計學分析采用Stata 12.0和Revman 5.2統(tǒng)計軟件包進行數(shù)據(jù)分析。在顯性模型、隱性模型和等位基因模型下合并計算比值比(odds ratio，OR)和95%可信區(qū)間(confidence interval，CI)，評價IL12B 1188A/C多態(tài)性與GD的關(guān)聯(lián)性。合并計算前先進行異質(zhì)性檢驗， 若P>0.1，I2<50%，表示各研究間沒有顯著異質(zhì)性，采用固定效應(yīng)模型(Mantel-Haenszel法)合并計算效應(yīng)量；若，P≤0.10，I2≥50%，表示存在異質(zhì)性，則采用隨機效應(yīng)模型(Dersimonia-Laird法)進行合并分析[4]。采用Egger’s及Begg’s檢驗觀察各研究間是否存在發(fā)表偏倚，合并后的檢驗水準α=0.10。

2結(jié)果

2.1文獻檢索結(jié)果初次檢索獲得中文及英文文獻共9篇[5-13]，排除2篇未做1188A/C多態(tài)位點的文獻[5-6]，1篇無正常對照的文獻[7]，2篇數(shù)據(jù)重復(fù)的文獻選擇1篇[8-9]，篩選后最終納入5篇文獻，包括747例病例和698例對照。納入文獻的基本情況及多態(tài)位點基因型分布情況(見表1)。

表1納入文獻的基本情況

第一作者年份國別病例AAACCC對照AAACCC哈代平衡PIkeda等[10]Hiromatsu等[11]Yang等[12]孫玲玲等[13]Anvari等[8-9]20042006200620072010日本日本日本中國伊朗223929771629026604254345345320296430391206726612869232726080.2040.9031.3920.5000.4310.8060.3900.8230.9780.613

2.2納入研究的質(zhì)量評價采用紐爾斯卡-渥太華量表(the Newcastle-Ottawa Scale, NOS)對病例-對照研究的評價標準對納入的文獻進行質(zhì)量評價[14]。評價內(nèi)容包括3方面：研究人群選擇4分，組間可比性2分，暴露因素的測量3分，共9分。質(zhì)量評分結(jié)果顯示文獻中的Graves病臨床診斷明確，對照排除了自身免疫性疾病史。3個研究未提及病例-對照匹配，2個給出了病例組合對照組性別頻數(shù)分布，經(jīng)χ2檢驗顯示男女分布差異無統(tǒng)計學意義(P>0.05)。納入的5篇文獻4篇得分7分，1篇得分為6分。

2.3Meta分析結(jié)果5個研究的異質(zhì)性檢驗顯示各研究之間存在顯著的異質(zhì)性(P<0.10,I2>50%)，因此采用隨機效應(yīng)模型進行Meta分析。結(jié)果顯示IL12B 1188A/C多態(tài)位點在3種遺傳模式下，合并OR值經(jīng)檢驗均無統(tǒng)計學意義(P>0.05，見圖1)。

對納入文獻以研究人群的國別為分層因素進行了亞層分析。其中，對3個日本人群的研究進行亞組分析結(jié)果未發(fā)現(xiàn)1188A/C多態(tài)位點的變異與GD發(fā)病有統(tǒng)計學相關(guān)性(顯性遺傳模型：OR=0.78，95% CI=0.45～1.37，P=0.39；隱性遺傳模型：OR=1.15，95% CI=0.87～1.52，P=0.31；等位基因模型：OR=0.99，95% CI=0.83～1.18，P=0.92)。

A：顯性模型 ；B：隱性模型；C：等位基因模型。圖1　IL12B 1188A/C多態(tài)性與GD易感森林圖

2.4發(fā)表偏倚評估發(fā)表偏倚產(chǎn)生的主要原因是陰性結(jié)果的研究不易發(fā)表從而使文獻檢索向陽性結(jié)果的研究偏斜[15]，最終影響Meta分析的結(jié)果。Egger’s及Begg’s檢驗提示各研究無統(tǒng)計學意義的發(fā)表偏倚(P>0.10)。本研究只有5個研究，故未做漏斗圖[16]。

3討論

研究表明，Th1/Th2細胞失衡在GD發(fā)病中具有重要作用，IL12是決定Th細胞分化的重要細胞因子[17]。IL12又稱自然殺傷細胞刺激因子(NKSF)，是一種由2個亞基p35和p40通過二硫鍵連接的異源二聚體型細胞因子。IL12可刺激T細胞增殖，促進Th1細胞的生長、分化 。IL12的2個亞基基因定位于不同的染色體上，編碼p40的基因IL12B定位于5q31-33。據(jù)報道，IL12B基因多態(tài)性與糖尿病、哮喘、精神分裂癥等[17-19]自身免疫性疾病相關(guān)聯(lián)。動物模型研究發(fā)現(xiàn)IL12參與GD和HT發(fā)病，GD患者IL12水平顯著升高[17]。目前關(guān)于IL12B基因多態(tài)性與GD易感性的研究很多，絕大多數(shù)集中在3′UTR區(qū)能夠影響IL12B轉(zhuǎn)錄的功能SNP 1188A/C[19-21]，但該多態(tài)性與GD發(fā)病是否相關(guān)，研究結(jié)果尚不一致。Ikeda等[10]和Hiromatsu等[11]在日本人群的2個研究提示IL12B基因3′UTR 1188A/C多態(tài)與Graves病的發(fā)病沒有關(guān)聯(lián)，這與Yang等[12]的研究結(jié)果一致，但Anvari等[8-9]的研究卻發(fā)現(xiàn)1188A/C多態(tài)與GD、GO發(fā)病相關(guān)，1188C等位基因和CC基因型發(fā)病率顯著增高。

目前關(guān)于IL12B 基因多態(tài)性與Graves病的相關(guān)研究均為小樣本研究，統(tǒng)計效力較低，Meta分析通過對多個目的相同相互獨立的研究結(jié)果進行綜合評價分析，以增加樣本量提高檢驗效能，已成為探索疾病與基因關(guān)聯(lián)常用的研究方法。本研究采用Meta分析方法對IL12B基因1188A/C多態(tài)性與GD的易感性進行定量合并分析，納入符合標準的文獻5篇，所有研究的對照組基因型分布均符合Hardy-Weinberg平衡。因納入文獻只有5篇，故未做漏斗圖分析發(fā)表偏倚，Egger’s及Begg’s檢驗提示尚不能認為存在發(fā)表偏倚。

本研究結(jié)果顯示，IL12B基因1188A/C多態(tài)性在3種遺傳模型下與GD均無統(tǒng)計學意義相關(guān)，與GD總體發(fā)病風險關(guān)系不大。但對納入研究以國別為分層因素進行分析，發(fā)現(xiàn)在3個日本人群的研究中，3種遺傳模型下的差異均無統(tǒng)計學意義。但由于關(guān)于IL12B基因多態(tài)性與GD的關(guān)聯(lián)研究較少，且均為小樣本研究，5個研究總樣本量只有不足1 500人；在文獻數(shù)據(jù)收集中，部分文獻由于無法提供完整數(shù)據(jù)而損失；因此IL12B基因3′UTR 1188A/C多態(tài)性與GD易感性的關(guān)系仍需擴大樣本進行進一步研究。

[參考文獻]

[1] Jin Y, Teng W, Ben S, et al. Genome-wide scan of Graves' disease: evidence for linkage on chromosome 5q31 in Chinese Han pedigrees[J]. J Clin Endocrinol Metab, 2003, 88(4): 1798-1803.

[2] Akamizu T, Hiratani H, Ikegami S, et al. Association study of autoimmune thyroid disease at 5q23-q33 in Japanese patients[J]. J Hum Genet, 2003, 48(5): 236-242.

[3] Sakai K, Shirasawa S, Ishikawa N, et al. Identification of susceptibility loci for autoimmune thyroid disease to 5q31-q33 and Hashimoto's thyroiditis to 8q23-q24 by multipoint affected sib-pair linkage analysis in Japanese[J]. Hum Mol Genet, 2001, 10(13): 1379-1386.

[4] 胡世蕓, 陳妍雯, 潘宗琴, 等. 我國肝衰竭醫(yī)院感染危險因素的Meta分析[J]. 遵義醫(yī)學院學報, 2014, 37(5): 520-528.

[5] Liu S, Hou P, Morahan G, et al. No association between IL12B gene polymorphisms and Graves' disease in the Chinese population[J]. Endocr Res, 2012, 37(4): 182-187.

[6] Liu Y H, Chen C C, Liao L L, et al. Association of IL12B polymorphisms with susceptibility to Graves ophthalmopathy in a Taiwan Chinese population[J]. J Biomed Sci, 2012,19(1)：97-104.

[7] Walsh J P, Berry J, Liu S, et al. The clinical presentation of autoimmune thyroid disease in men is associated with IL12B genotype[J]. Clin Endocrinol (Oxf), 2011, 74(4): 508-512.

[8] Anvari M, Khalilzadeh O, Esteghamati A, et al. Graves’ disease and gene polymorphism of TNF-α, IL-2, IL-6, IL-12, and IFN-γ[J]. Endocrine, 2010, 37(2): 344-348.

[9] Anvari M, Khalilzadeh O, Esteghamati A, et al. Genetic susceptibility to Graves’ ophthalmopathy: the role of polymorphisms in proinflammatory cytokine genes[J]. Eye (Lond), 2010, 24(6): 1058-1063.

[10] Ikeda Y, Yoshida W, Noguchi T, et al. Lack of association between IL-12B gene polymorphism and autoimmune thyroid disease in Japanese patients[J]. Endocr J, 2004, 51(6): 609-613.

[11] Hiromatsu Y, Fukutani T, Ichimura M, et al. Interleukin-12B gene polymorphism does not confer susceptibility to Graves’ ophthalmopathy in Japanese population[J]. Endocr J, 2006, 53(6): 753-759.

[12] Yang J M, Nagasaka S, Yatagai T, et al. Interleukin-12p40 gene (IL-12B) polymorphism and Type 1 diabetes mellitus in Japanese: possible role in subjects without having high-risk HLA haplotypes[J]. Diabetes Res Clin Pract, 2006, 71(2): 164-169.

[13] 孫玲玲，褚迅，黃薇. 華東地區(qū)漢族人群IL 12B基因多態(tài)性與Graves病相關(guān)性研究[J]. 國際遺傳學雜志, 2007, 30(4): 258-260.

[14] Wells G, Shea B, O’Connell D, et al. NewCastle-Ottawa Quality Assessment Scale—Case Control Studies[EB/OL].[2012-06-15]. http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp.

[15] 楊書, 楊曉虹, 劉新. 發(fā)表性偏倚產(chǎn)生的原因及可行性論證[J]. 成都醫(yī)學院學報, 2008, 3(2): 132-135.

[16] Lau J, Ioannidis J P, Terrin N, et al. The case of the misleading funnel plot[J]. BMJ, 2006, 333(7568): 597-600.

[17] Fountoulakis S, Tsatsoulis A. On the pathogenesis of autoimmune thyroid disease: a unifying hypothesis[J]. Clin Endocrinol (Oxf), 2004, 60(4): 397-409.

[18] Morahan G, Huang D, Ymer S I, et al. Linkage disequilibrium of a type 1 diabetes susceptibility locus with a regulatory IL12B allele[J]. Nat Genet, 2001, 27(3): 218-221.

[19] Morahan G, Huang D, Wu M, et al. Association of IL12B promoter polymorphism with severity of atopic and non-atopic asthma in children[J]. Lancet, 2002, 360(9331): 455-459.

[20] Ozbey U, Tug E, Kara M, et al. The value of interleukin-12B (p40) gene promoter polymorphism in patients with schizophrenia in a region of East Turkey[J]. Psychiatry Clin Neurosci, 2008, 62(3): 307-312.

[21] Seegers D, Zwiers A, Strober W, et al. A TaqI polymorphism in the 3'UTR of the IL-12 p40 gene correlates with increased IL-12 secretion[J]. Genes Immun, 2002, 3(7): 419-423.

[收稿2015-12-29;修回2016-03-08]

(編輯：陳婭)

[基金項目]國家自然科學基金青年項目(NO：81402351)；安徽省重點實驗室項目(NO：1306C083008)。

[中圖法分類號]R581

[文獻標志碼]B

[文章編號]1000-2715(2016)03-0293-04

Meta-analysis of the association between IL 12B gene polymorphism and Graves’ disease

SunLingling,DongBohan,LingLiefeng

(Department of Biochemistry, Wannan Medical College, Wuhu Anhui 241002, China)

[Abstract]Objective To assess the association between 1188 A/C polymorphism of IL12B 3′UTR and Graves’ disease (GD) by meta-analysis.Methods The relevant studies about the correlation between IL12B 3′UTR 1188A/C polymorphism and GD were collected by searching PubMed, EMbase, CNKI, VIP and CBM according to the inclusion criteria. Heterogeneity tests and meta-analysis were conducted.Results Five studies consisting of 747 cases and 698 controls were included. The results of meta-analysis showed that there was no significant difference between the 1188A/C polrmorphism of IL12B 3′UTR and GD (for dominant genetic model: OR=1.18, 95% CI=0.24～2.17, P=0.59; for recessive genetic mode: OR=1.61, 95% CI=0.97～2.67, P=0.06; for allelic genetic model: OR=1.32, 95% CI=0.93～1.87, P=0.12). In the 3 Japanese studies of the correlation between 1188 A/C polymorphism and the susceptibility of GD, no statistically signicant correlation was found (for dominant genetic model: OR=0.78, 95% CI=0.45～1.37, P=0.39; for recessive genetic mode: OR=1.15, 95% CI=0.87～1.52, P=0.31; for allelic genetic model: OR=0.99, 95% CI=0.83～1.18, P=0.92).Conclusion There is no significant correlation between IL 12B 3′UTR 1188A/C polymorphism and GD.

[Key words]Graves disease; IL12B; gene polymorphism; Meta-analysis