白細胞介素2基因-330T/G和白細胞介素6基因-174G/C多態性與結核病易感性的Meta分析

李艷林 黎亮 唐秀芳 吳秀娟

結核病發病率僅以每年1.5%的速度低速下降，仍是全球死亡率最高的十大疾病之一，嚴重威脅著人類的健康[1]。據世界衛生組織的報道，2015年新發結核病患者為140萬例，結核病死亡患者數為140萬例[2]。流行病學研究顯示，全球約有1/3的人群感染過結核分枝桿菌(Mycobacteriumtuberculosis，MTB)，但僅有10%的感染者發展為活動性結核病。結核病的發病是MTB、宿主遺傳因素和環境因素共同影響的結果，宿主易感性在結核病的發生發展中起重要作用[3]。

目前，已有學者進行了白細胞介素2(IL-2)和IL-6基因啟動子區域中的多態性位點與結核病易感性的關聯研究，但結論不一。為了客觀評價IL-2基因-330T/G和IL-6基因-174G/C與結核病易感性的關系，筆者采用Meta分析的方法，通過增大樣本量，對既往研究結果進行系統、定量、綜合性分析，來增加結論的可信度，解決研究結果的不一致性。

材料和方法

1.文獻檢索策略：分別以“IL-2”或“Cytokine”或“Interleukine-2”和“TB”或“Tuberculosis”或“TB infection”，以及“Polymorphism”或“genotype”為關鍵詞，檢索數據庫PubMed、Medline、EMBase、Web of Science和Elsevier Science Direct中公開發表的IL-2與結核病易感性相關的病例-對照研究文獻。分別以“IL-6”或“Cytokine”或“Interleukine-6”和“TB”或“Tuberculosis”或“TB infection”，以及“Polymorphism”或“genotype”為關鍵詞，檢索公開發表的IL-6與結核病易感性相關的病例-對照研究文獻。檢索語種為“英語”。以“白細胞介素”“細胞因子”“結核”“多態性”和“易感性”為檢索詞，檢索1990年1月至2017年6月在萬方數據庫、中國知網(CNKI)和維普數據庫公開發表的IL-2、IL-6與結核病易感性相關的病例-對照研究文獻。檢索語種為“漢語”。初步檢索共獲得相關文獻573篇，根據納入排除標準剔除文獻561篇，最終納入文獻12篇[4-15]。

2.納入排除標準：(1)納入標準：①研究內容為結核病患者與健康對照者的IL-2基因-330T/G和(或)IL-6基因-174G/C位點多態性；②病例-對照研究，報告數據完整，文獻包含各個基因型的分布頻數；③有各基因型與結核病患病風險性的OR(95%CI)值。(2)排除標準：①非隨機對照研究；②系統評價、綜述、病例報告等無法利用的文獻；③具體數據描述不清，僅有摘要或會議匯編，無法得到原始數據的資料；④重復報道的文獻取其一；⑤僅有基因型頻率，未提及基因型頻數的研究。

3.資料提取：由2名研究者分別獨立進行文獻檢索，并評價所獲文獻，不同的檢索結果經討論解決，最終達成一致。提取的信息包括：作者、發表年份、樣本量、研究對象所屬地區、納入研究的多態性位點的基因型頻率及等位基因頻率等。

4.文獻質量評價：采用Begg檢驗與Egger檢驗評估納入研究文獻的發表偏倚，P<0.05提示可能存在發表偏倚。納入的所有研究根據Newcastle-Ottawa Scale(NOS)原則進行評估，評分≥7分提示文獻質量較好。

5.統計學分析：運用Stata 12.0軟件進行數據處理，采用綜合定量分析并繪制Meta分析的森林圖。對各研究結果的異質性進行檢驗，若各研究間異質性差異無統計學意義(P>0.05)，則采用固定效應模型進行數據合并；若各研究間異質性差異有統計學意義(P<0.05)，則采用隨機效應模型進行數據合并。計算OR(95%CI)值評價IL-2基因-330T/G 和IL-6基因-174G/C多態性與結核病易感性的關系，以P<0.05為差異有統計學意義。

結 果

1.基本情況：納入的12篇文獻中，累計病例組2416例，健康對照組3474名。納入文獻中研究對象為亞洲人群的有6篇，歐洲人群的有4篇，美洲人群的有2篇。納入文獻基本信息特點見表1。本研究納入的12篇文獻采用NOS系統進行評分，均≥7分。經Begg檢驗與Egger檢驗，未發現各項研究間存在明顯發表偏倚，因此認為所納入文獻質量較高。

2.異質性分析：將收集的文獻資料根據不同地區進行異質性分析。IL-2基因-330T/G位點各研究間均存在異質性，且差異有統計學意義(圖1顯示為-330T/G等位基因模型異質性分析：I2=70.3%，P<0.05)，采用隨機效應模型進行Meta分析。IL-6 基因-174G/C位點各研究間存在一定的異質性，但差異無統計學意義(圖2為-174G/C等位基因模型異質性分析：I2=18.2%，P>0.05)，故采用固定效應模型。

3. Meta分析結果：IL-2基因-330T/G位點納入研究共8篇，病例組971 例，對照組1519名。Meta分析結果顯示，IL-2基因-330T/G多態性可能與結核病患病風險降低有關(T對G：OR=0.75，95%CI=0.59～0.95，P=0.019)。亞組分析顯示，-330T等位基因與歐洲人群結核病患病風險降低有關(T對G：OR=0.55，95%CI=0.31～0.98，P=0.042)，見表2、圖1。

注NOS：Newcastle-Ottawa Scale量表；HWE：Hardy-Weinberg平衡

I2為異質性部分在總變異中所占的比重；分組分析時美洲組參考文獻僅1篇，故無I2和P值數據圖1 白細胞介素2基因-330T/G等位基因模型森林圖

地區等位基因模型(T對G)純合子模型(TT對GG)雜合子模型(TG對GG)顯性模型(TT+TG對GG)隱性模型(TT對TG+GG)OR(95%CI)值P值OR(95%CI)值P值OR(95%CI)值P值OR(95%CI)值P值OR(95%CI)值P值總體075(059～095)0019057(032～102)0056097(056～167)0902077(047～125)0290062(042～093)0020亞洲081(060～110)0174066(028～154)0331083(032～216)0702078(032～191)0581074(050～110)0134歐洲055(031～098)0042036(013～100)0049108(060～198)0791066(038～113)0127036(012～105)0061美洲102(074～140)0927113(056～227)0728123(061～250)0566118(060～229)0634096(063～147)0857

IL-6基因-174G/C位點納入研究共9篇，病例組1445例，對照組1955名。Meta分析結果顯示，IL-6基因-174G/C多態性可能與結核病患病風險增加相關(G對C：OR=1.40，95%CI=1.21～1.62，P<0.01；GG對GC：OR=1.98，95%CI=1.36～2.87，P<0.01；GG+GC對CC：OR=1.75，95%CI=1.22～2.52，P=0.002；GG對GC+CC：OR=1.44，95%CI=1.21～1.72，P<0.01)。亞組分析顯示，-174G等位基因與美國人群(G對C：OR=1.65，95%CI=1.28～2.13，P<0.01；GG對GC：OR=2.64，95%CI=1.39～5.01，P=0.003；GG+GC對CC：OR=2.20，95%CI=1.18～4.12，P=0.014；GG對GC+CC：OR=1.76，95%CI=1.28～2.41，P<0.01)和亞洲人群(G對C：OR=1.40，95%CI=1.13～1.72，P=0.002；GG對GC：OR=2.15，95%CI=1.17～3.92，P=0.013；GG+GC對CC：OR=1.85，95%CI=1.03～3.32，P=0.040；GG對GC+CC：OR=1.43，95%CI=1.12～1.83，P=0.005)結核病患病風險增加有關，見表3、圖2。

4.發表偏倚分析：各研究間均無統計學意義上的發表偏倚，結果具有一定的可信度(圖3，4)。

I2為異質性部分在總變異中所占的比重圖2 白細胞介素6基因-174G/C等位基因模型森林圖

地區等位基因模型(G對C)純合子模型(GG對CC)雜合子模型(GC對GG)顯性模型(GG+GC對CC)隱性模型(GG對GC+CC)OR(95%CI)值P值OR(95%CI)值P值OR(95%CI)值P值OR(95%CI)值P值OR(95%CI)值P值總體140(121～162)<001198(136～287)<001146(100～214)0049175(122～252)0002144(121～172)<001亞洲140(113～172)0002215(117～392)0013151(083～275)0179185(103～332)0040143(112～183)0005歐洲109(079～150)0594125(062～252)0528123(060～253)0573125(064～246)0518106(070～161)0773美洲165(128～213)<001264(139～501)0003163(084～315)0149220(118～412)0014176(128～241)<001

縱坐標為OR的對數值，橫坐標為log[OR]的標準誤圖3 白細胞介素2等位基因模型納入分析文獻Begg檢驗發表偏倚圖

縱坐標為OR的對數值，橫坐標為log[OR]的標準誤值圖4 白細胞介素6等位基因模型納入分析文獻Begg檢驗發表偏倚圖

討 論

大量研究證實，細胞因子在結核病的發展中發揮著重要的作用。肺結核的易感性及其疾病的進展可能是不同染色體的多個等位基因共同作用的結果[16]。細胞因子基因中的多態性影響基因的轉錄表達，從而導致了細胞因子之間的個體差異[17]。

IL-2和IL-6作為細胞因子共同作用以促進機體對MTB的免疫應答[18-19]。IL-2在感染初期通過效應細胞擴增產生的γ干擾素在抗結核免疫中起重要作用。此外，在感染后期IL-2通過控制炎癥反應促進調節性T細胞的生長，并延長其存活時間[20]。IL-6作為關鍵的促炎細胞因子，在活動性肺結核的發展中發揮作用，促進結核肉芽腫的形成，并誘導感染MTB的巨噬細胞凋亡[21]。結核病患者IL-2和IL-6的血清水平明顯升高，揭示IL-2和IL-6可能作為結核病的潛在生物標志物[22]。

從本研究結果可看出，歐洲人群的IL-2基因-330T/G等位基因模型中攜帶等位基因T者結核病患病率低于攜帶等位基因G者，提示IL-2基因-330T/G的T等位基因可能降低歐洲人群結核病患病風險。IL-6基因-174G/C與亞洲和美洲人群結核病患病風險增加有明顯相關性，提示IL-6基因-174G/C的G等位基因可能增加亞洲和美洲人群結核病患病風險。

Mao等[23]通過Meta分析發現，IL-6基因-174G/C多態性與亞洲人結核病風險的增加明顯相關。筆者經過Meta分析不僅發現-174G/C基因型的存在與亞洲人群結核病患病風險增加有關，同時還和美洲人群結核病患病風險增加有關。究其原因，可能為樣本量的大小及納入文獻量的不同導致的Meta分析結果的差異。

文獻質量評價結果顯示，本次Meta分析納入文獻的質量較高，其結果也具有一定的可信度，但本研究也存在一定的局限性。首先，本次分析中，4項研究[4-5,10-11]中對照的基因型分布不符合HWE，表明這些研究也可能導致異質性。其次，所收集的研究數量較少。第三，亞組分析和Meta回歸分析受到限制，因為一些研究沒有提供有關年齡、性別和其他可能與結核病患病風險相關的全面信息。因此，上述結論仍需更多高質量、大規模的研究進一步證實。

[1] Cazabon D, Alsdurf H, Satyanarayana S, et al. Quality of tuberculosis care in high burden countries: the urgent need to address gaps in the care cascade. Int J Infect Dis, 2017, 56: 111-116.

[2] World Health Organization. Global tuberculosis report 2016. Geneva: World Health Organization, 2016.

[3] 江道斌, 齊曼古力·吾守爾. 巨噬細胞功能相關的易感基因多態性與結核病相關性的研究進展. 中國防癆雜志, 2014, 36(5): 385-387.

[4] Oral HB, Budak F, Uzaslan EK, et al. Interleukin-10 (IL-10) gene polymorphism as a potential host susceptibility factor in tuberculosis. Cytokine, 2006, 35(3/4): 143-147.

[5] Amirzargar AA, Rezaei N, Jabbari H, et al. Cytokine single nucleotide polymorphisms in Iranian patients with pulmonary tuberculosis. Eur Cytokine Netw, 2006, 17(2): 84-89.

[6] Selvaraj P, Alagarasu K, Harishankar M, et al. Cytokine gene polymorphisms and cytokine levels in pulmonary tuberculosis. Cytokine, 2008, 43(1): 26-33.

[7] Naslednikova IO, Urazova OI, Voronkova OV, et al. Allelic polymorphism of cytokine genes during pulmonary tuberculosis. Bull Exp Biol Med, 2009, 148(2): 175-180.

[8] Trajkov D, Trajchevska M, Arsov T, et al. Association of 22 cytokine gene polymorphisms with tuberculosis in Macedo-nians. Indian J Tuberc, 2009, 56(3): 117-131.

[9] Sivangala R, Ponnana M, Thada S, et al. Association of cytokine gene polymorphisms in patients with tuberculosis and their household contacts. Scand J Immunol, 2014, 79(3): 197-205.

[10] Hu Y, Wu L, Li D, et al. Association between cytokine gene polymorphisms and tuberculosis in a Chinese population in Shanghai: a case-control study. BMC Immunol, 2015, 16: 8.

[11] Butov DO, Kuzhko MM, Makeeva NI, et al. Association of interleukins genes polymorphisms with multi-drug resistant tuberculosis in Ukrainian population. Pneumonol Alergol Pol, 2016, 84(3): 168-173.

[12] Milano M, Moraes MO, Rodenbusch R, et al. Single nucleotide polymorphisms in IL17A and IL6 are associated with decreased risk for pulmonary tuberculosis in southern Brazilian population. PLoS One, 2016, 11(2): e0147814.

[13] Henao MI, Montes C, París SC, et al. Cytokine gene polymorphisms in Colombian patients with different clinical presentations of tuberculosis. Tuberculosis (Edinb), 2006, 86(1): 11-19.

[14] Ansari A, Hasan Z, Dawood G, et al. Differential combination of cytokine and interferon-γ +874 T/A polymorphisms determines disease severity in pulmonary tuberculosis. PLoS One, 2011, 6(11): e27848.

[15] Zhang G, Zhou B, Wang W, et al. A functional single-nucleotide polymorphism in the promoter of the gene encoding interleukin 6 is associated with susceptibility to tuberculosis. J Infect Dis, 2012, 205(11): 1697-1704.

[16] Hill AV. The immunogenetics of human infectious diseases. Annu Rev Immunol, 1998, 16: 593-617.

[17] Pan F, Tian J, Pan YY, et al. Association of IL-10-1082 promoter polymorphism with susceptibility to gastric cancer: evidence from 22 case-control studies. Mol Biol Rep, 2012, 39(6): 7143-7154.

[18] Flynn JL, Chan J. Immunology of tuberculosis. Annu Rev Immunol, 2001, 19: 93-129.

[19] Khanna M, Srivastava LM, Kumar P. Defective interleukin-2 production and interleukin-2 receptor expression in pulmonary tuberculosis. J Commun Dis, 2003, 35(2): 65-70.

[20] Millington KA, Innes JA, Hackforth S, et al. Dynamic relationship between IFN-gamma and IL-2 profile ofMycobacteriumtuberculosis-specific T cells and antigen load. J Immunol, 2007, 178(8): 5217-5226.

[21] 陳興年, 劉奇棟, 王泓, 等. 活動性肺結核患者血清中白細胞介素6、8和腫瘤壞死因子α檢測的臨床意義. 中國防癆雜志, 2015, 37(4): 344-347.

[22] Chowdhury IH, Choudhuri S, Sen A, et al. Serum interleukin 6 (IL-6) as a potential biomarker of disease progression in active pulmonary tuberculosis following anti-tuberculosis drug therapy. Mol Immunol, 2015, 63(2): 601-602.

[23] Mao X, Ke Z, Liu S, et al. IL-1β +3953C/T, -511T/C and IL-6 -174C/G polymorphisms in association with tuberculosis susceptibility: a meta-analysis. Gene, 2015, 573(1): 75-83.