黃晶，呂學鋒，張晉欣，張春艷

卒中是導致成年人殘疾的最常見疾病，也是導致成年人死亡的第二大疾病[1]。其中，缺血性卒中（ischemic stroke，IS）患者占所有卒中患者的70%～80%[2]。傳統(tǒng)的IS高危因素包括高血壓、高脂血癥、糖尿病、吸煙以及肥胖，近些年，有研究發(fā)現(xiàn)，IL-10水平較低的人群IS發(fā)病風險明顯更高[3]。IL-10基因位于1號染色體1q31-q32的位置，由5個外顯子和4個內(nèi)含子組成[4]，其啟動子區(qū)的2種單核苷酸多態(tài)性：IL-10-1082A/G和IL-10-819T/C，對IL-10的mRNA表達有影響，并與IS的發(fā)病有關[5]，但上述結論尚存爭議。本研究采用系統(tǒng)評價的方法，明確IL-10啟動子區(qū)基因（IL-10-1082A/G和IL-10-819T/C）多態(tài)性與IS易感性的關系，為IS的預防及治療提供依據(jù)。

1 研究對象與方法

1.1 文獻檢索策略 計算機檢索Pubmed、Embase、Web of Science、萬方及中國知網(wǎng)數(shù)據(jù)庫，檢索時限為建庫至2019年2月，英文檢索詞：“interleukin-10”or“IL-10”or“IL 10”AND“stroke”or“cerebral infarction”or“cerebral thrombosis”。中文檢索詞：“白細胞介素10”“白介素10”“IL 10”“腦梗死”“卒中”“腦血栓形成”。

1.2 納入標準和排除標準 納入標準：①研究內(nèi)容為IL-10基因多態(tài)性與IS發(fā)病風險相關性的研究；②研究類型為病例-對照研究；③病例組研究對象為確診的IS患者，對照組為健康人群；④病例組和對照組均有詳細或可以計算出的等位基因和基因型分布數(shù)據(jù)；⑤若相同的數(shù)據(jù)發(fā)表多篇文章，納入最近發(fā)表的研究。

排除標準：①重復發(fā)表的研究；②綜述或評論性文章；③動物實驗；④對照組基因型分布不符合哈代-溫伯格平衡定律（Hardy-Weinberg equilibrium，HWE）。

1.3 資料提取 由2名評價者根據(jù)納入標準和排除標準獨立篩選文獻、提取資料，并交叉核對，若意見不同，則重新審閱全文并相互討論，難以決定時，交第三方裁決。需提取的數(shù)據(jù)如下：第一作者姓名；研究對象種族和國籍；文獻發(fā)表年份；對照組來源；基因分型方法；兩組間年齡和性別是否匹配；病例組和對照組的總例數(shù)；病例組和對照組的等位基因及基因型分布的數(shù)據(jù)。

1.4 質量評價 由2名評價者參照Cochrane手冊的質量評價標準對納入研究進行質量評價。文獻質量評分使用紐卡斯爾-渥太華量表（Newcastle-Ottawa Scale，NOS），其中包括三個方面：選擇、可比性和暴露，總評分為9分。

1.5 統(tǒng)計學方法 使用RevMan 5.3軟件進行Meta分析。采用以下5組基因模型分析IL-10-1082A/G和IL-10-819T/C基因多態(tài)性與IS患病風險間的相關性：等位基因模型（GvsA），（CvsT）；顯性模型（GG+AGvsAA），（CC+CTvsTT）；隱性模型（GGvsAG+AA），（CCvsCT+TT）；純合子模型（GGvsAA），（CCvsTT）和雜合子模型（AGvsAA），（CTvsTT）。結果采用OR及95%CI作為效應量表示，根據(jù)種族劃分，對中國人群進行亞組分析，P＜0.05為差異有統(tǒng)計學意義。采用χ2檢驗對研究結果進行異質性檢驗，P＜0.10或I2＞50%時，認為各研究間有異質性，選用隨機效應模型進行分析；P＞0.10且I2＜50%，則選用固定效應模型。根據(jù)對照組的基因型分布進行HWE檢驗，定義P＞0.05為對照組符合HWE。最后，繪制漏斗圖，并使用State 14.0軟件進行Egger’s回歸分析，評估發(fā)表偏倚。

2 結果

2.1 文獻納入情況 本文最終納入13篇文獻[6-18]，均為病例對照研究，12篇文獻[6-13，15-18]研究了IL-10-1082A/G基因多態(tài)性與IS發(fā)病風險的相關性，6篇文獻[7，11-12，14，17-18]研究了IL-10-819T/C基因多態(tài)性與IS發(fā)病風險的相關性。共5328例受試者，其中病例組2822例，對照組2506例。文獻篩選流程見圖1，納入文獻基線資料及文獻質量評分見表1，基因型的分布情況以及HWE檢驗結果見表2。

2.2 薈萃分析結果

2.2.1IL-10-1082A/G基因多態(tài)性與缺血性卒中發(fā)病風險的相關性 5組基因模型Meta分析匯總見表3。在總體人群中，GGvsAA基因模型中各研究間無明顯異質性，采用固定效應模型；其余4組基因模型中各研究間均存在異質性，采用隨機效應模型。結果顯示，在各個遺傳模型中，IL-10-1082A/G基因多態(tài)性與IS發(fā)病風險明顯相關（GvsA：OR0.71，95%CI0.59～0.86，P＜0.001；GGvsAA：OR0.61，95%CI0.49～0.76，P＜0.001；AGvsAA：OR0.72，95%CI0.55～0.94，P=0.020；GG+AGvsAA：OR0.68，95%CI0.53～0.87，P=0.002；GGvsAG+AA：OR0.68，95%CI0.52～0.89，P=0.005）（圖2）。按種族劃分對中國人群進行亞組分析同樣發(fā)現(xiàn)，IL-10-1082A/G基因多態(tài)性與IS易感性明顯相關。

表1 納入文獻基線資料及文獻評分

表2 基因型分布及HWE檢驗

表3 IL-10啟動子基因（IL-10-1082A/G）多態(tài)性與缺血性卒中發(fā)病風險相關性的Meta分析

2.2.2IL-10-819T/C基因多態(tài)性與缺血性卒中發(fā)病風險的相關性 5組基因模型Meta分析匯總見表4。所有遺傳模型中均未發(fā)現(xiàn)異質性，采用固定效應模型。結果顯示：在各個遺傳模型中，IL-10-819T/C基因多態(tài)性與IS發(fā)病風險無明顯關聯(lián)（CvsT：OR1.11，95%CI0.97～1.29，P=0.140；CCvsTT：OR1.17，95%CI0.87～1.57，P=0.300；CTvsTT：OR1.09，95%CI0.88～1.36，P=0.420；GG+AGvsAA：OR1.13，95%CI0.93～1.39，P=0.230；GGvsAG+AA：OR1.16，95%CI0.89～1.50，P=0.270）（圖3）。按種族劃分對中國人群進行亞組分析同樣發(fā)現(xiàn)，IL-10-819T/C基因多態(tài)性與IS易感性無明顯相關性。

2.3 敏感性分析 敏感性分析采用逐一剔除各項研究后，重新進行Meta分析的方法。結果發(fā)現(xiàn)在雜合子模型（AGvsAA）中，剔除Tuttolomondo等[7]的研究后，與原結果（OR0.72，95%CI0.55～0.94，P=0.02）相比，發(fā)生了改變（OR0.76，95%CI0.59～0.99，P=0.05），尚不能下肯定結論。在其他模型中，剔除任一研究，未發(fā)現(xiàn)結果產(chǎn)生明顯改變。

2.4 發(fā)表偏倚 采用漏斗圖及Egger’s回歸分析評估研究的發(fā)表偏倚。漏斗圖結果顯示，納入研究基本對稱分布于合并OR值的兩側（圖4）。Egger’s檢驗結果顯示，所有基因模型中，僅有純合子模型（CCvsTT）存在明顯發(fā)表偏倚（P＜0.001）（表3、表4）。

圖2 總體人群IL-10-1082A/G基因多態(tài)性與缺血性卒中發(fā)病風險相關性的森林圖

表4 IL-10啟動子基因（IL-10-819T/C）多態(tài)性與缺血性卒中發(fā)病風險相關性的Meta分析

圖3 總體人群IL-10-819T/C基因多態(tài)性與缺血性卒中發(fā)病風險相關性的森林圖

圖4 IL-10基因多態(tài)性與缺血性卒中發(fā)病風險相關性的漏斗圖

3 討論

動脈粥樣硬化是動脈粥樣硬化性IS病理生理過程中必不可少的一環(huán)，而慢性炎癥貫穿動脈粥樣硬化發(fā)生發(fā)展過程的始終[19]。IL-10是一種具有免疫調節(jié)作用的細胞因子，主要由活化的單核細胞、淋巴細胞和巨噬細胞分泌[20]，可延緩血管內(nèi)動脈粥樣硬化進展[22]，從而具有腦缺血保護作用，其潛在機制可能是：①抑制一些炎癥因子的生成，如抑制巨噬細胞產(chǎn)生的IL-6、IL-1β、IL-1α、TNF-α，以及T細胞產(chǎn)生的γ干擾素（interferon-γ，IFN-γ）[21]；②促進細胞因子拮抗劑[如可溶性p55、p75腫瘤壞死因子受體（tumor necrosis factor receptor，TNFR）]的表達；③阻斷細胞因子的信號轉導通路。此外，在IS患者的血清、腦脊液中檢測到IL-10濃度的升高[23]，提示IL-10在IS發(fā)病中起著重要作用。

IL-10啟動子區(qū)基因（-1082、-819、-652、-592、-127、-41）存在6種不同的多態(tài)性[24]，多項研究表明，IL-10-1082A/G、IL-10-819T/C多態(tài)性可通過影響轉錄、翻譯從而影響IL-10的合成。目前國內(nèi)外有關IL-10基因多態(tài)性與IS發(fā)病風險的相關研究，結論不盡相同，進一步評估IL-10基因多態(tài)性與IS的關系，可為IS的防治及高危人群的篩查提供新的思路。本文收集現(xiàn)有研究進行薈萃分析，探討IL-10基因多態(tài)性與IS發(fā)病風險間的關系。

對IL-10-1082A/G基因多態(tài)性與IS易感性的Meta分析發(fā)現(xiàn)，以等位基因G為突變型，等位基因A為野生型，無論在哪個基因模型中，對照組暴露于G基因的人數(shù)均多于病例組，提示G基因可明顯減少人群中IS的發(fā)病風險。但除純合子模型（GGvsAA）外，其余基因模型各研究間均存在一定異質性。對中國人群進行亞組分析后，各組異質性均有不同程度的下降，但結果與總體人群一致。而所有基因模型中，IL-10-819T/C基因多態(tài)性在病例組和對照組中的差異無統(tǒng)計學意義，說明IL-10-819T/C基因多態(tài)性可能與IS的發(fā)病風險無明顯關聯(lián)。

此前已有Meta分析探討過IL-10啟動子區(qū)基因多態(tài)性與IS易感性的關系，但結論各不相同。Jin等[25]納入了7項病例對照研究進行薈萃分析，結果發(fā)現(xiàn)，在所有遺傳模型中，IL-10-1082A/G基因多態(tài)性與IS易感性均無明顯關聯(lián)，按照種族進行亞組分析后仍未發(fā)現(xiàn)二者的相關性。Yin等[26]的薈萃分析表明，在顯性基因模型（GG+GAvsAA）和等位基因模型（GvsA）中，IL-10-1082A/G多態(tài)性與IS發(fā)病風險顯著相關，但在隱性基因模型中，二者無明顯相關性；在所有基因模型中，IL-10-819T/C基因的多態(tài)性與IS發(fā)病風險無明顯關聯(lián)。導致本篇薈萃分析與上述研究結果不一致的原因可能包括種族差異、樣本量的大小、試驗設計不同。本文綜合現(xiàn)有研究進行薈萃分析，與上述研究相比，納入研究最多，樣本量最大，且針對中國人群進行了亞組分析，以最大程度減小異質性來源。同時，本文排除了所有不符合HWE檢驗的研究，增加了結果的可信度。

本文仍存在一些局限性：①在IL-10-1082A/G基因多態(tài)性與IS易感性的對比中，除純合子模型（GGvsAA）外，其余各基因模型中研究間均存在異質性，按照不同種族進行亞組分析后也未完全消除異質性，說明除外種族原因，還存在研究設計不同、基因分型方式不同等原因，導致了研究間的異質性；②IS的發(fā)病包含多種因素、多種基因的作用，本文只分析了一種因素，未考慮不同生存環(huán)境等混雜因素的影響；③本文納入的國內(nèi)研究較多，還需要其他國家的大樣本、高質量的研究來證實。

綜上所述，IL-10-1082A/G基因多態(tài)性與IS易感性顯著相關，而IL-10-819T/C基因多態(tài)性與IS發(fā)病風險無明顯關聯(lián)?？紤]到本研究尚存不足，未來還需大樣本、多中心、高質量的研究，進一步去除諸如肥胖、吸煙等混雜因素的影響，來進一步明確IL-10基因多態(tài)性與IS易感性的關聯(lián)。

【點睛】本文明確了IL-10啟動子基因多態(tài)性與卒中發(fā)病的相關性，對卒中高危人群篩查提供了依據(jù)。