IL-10基因-1082A/G多態性與糖尿病腎病易感性的Meta分析

胡 玲 王思思 周樂汀 王 涼 孫鑄興

約30%～40%的DM最終發展成為DN，DN是DM患者發生終末期腎病(end stage renal disease, ESRD)和死亡的主要原因[1]。研究表明，炎癥因素在DN發生及發展中起著重要作用[2，3]。Fathy等[4]研究發現，DN患者外周血IL-10水平顯著高于糖尿病非腎病患者。此外，Myliwska等[5]研究發現，DN患者IL-10水平的升高與蛋白尿水平呈正比，其水平與DN嚴重程度相關，這提示IL-10在DN病程中起著一定作用，但其機制仍不明確。

IL-10主要由Th2細胞、調節性T細胞產生，具有抗炎和免疫抑制等作用。IL-10基因位于染色體1q31-32上，-1082A/G多態性位于該基因的啟動子區域，其可改變IL-10基因的轉錄、調控mRNA的表達，影響IL-10的表達水平。目前有部分關于IL-10 -1082A/G多態性與DN發病風險的研究，但各研究結果并不一致[4,6～14]。本研究通過Meta分析合并多項研究結果，進一步探討IL-10基因-1082A/G多態位點與DN發生率之間的關系。

材料與方法

1.文獻檢索策略：通過對PubMed、Embase、Medline、中國知網、維普數據庫、萬方數據庫等全文數據庫檢索，檢索時間為建庫以來至2019年6月。英文關鍵詞主要包括：IL-10、IL-10 polymorphism、IL-10 gene、-1082A/G、rs1800896、SNP以及Diabetic nephropathy、diabetic kidney disease、DN、DKD。中文關鍵詞主要包括：IL-10、IL-10基因多態性、-1082A/G、rs1800896以及糖尿病腎病、糖尿病相關腎臟疾病、DN、DKD，將關鍵詞兩兩組合進行檢索。

2.納入與排除標準：納入標準：①均為基于人類的研究；②實驗組納入對象為DN患者，對照組為單純DM患者或健康人群；③根據文獻內容可獲得等位基因和基因型分布數據。排除標準：①不是基于人類的研究；②不是基于IL-10基因-1082A/G位點的研究；③重復報告或發表的研究；④等位基因及基因型分布數據無法獲得的研究。

3.數據提取與質量評價：文獻篩選和資料提取由兩名研究者獨立完成，意見不一致時通過討論或由第3位研究者協助解決。資料提取內容包括：第一作者、發表時間、地區、樣本量大小、等位基因分布情況、基因型分布情況等。通過NOS標準[15]對納入文獻進行質量評價，NOS標準包含3個部分：納入人群的選擇、群體可比性和暴露，總得分共9分。得分>7分被認為是高質量文獻，4～7分為中等質量文獻。

4.統計學方法：采用Stata12.0統計學軟件對數據進行統計分析，通過計算發病風險比值比(odds ratio, OR)及95%置信區間(confidence interval, CI)指標來評估IL-10基因-1082A/G多態位點與DN發生率的關系。使用的遺傳模型如下：等位基因模型(G vs A)、純合子模型(GG vs AA)、顯性模型(GG+GA vs AA)、隱性模型(GG vs GA+AA)、雜合子模型(GA vs AA)。本研究采用Q檢驗進行異質性分析，若P>0.1，提示研究之間異質性較小，采用固定效應模型進行分析，反之選用隨機效應模型。通過Begger檢驗進行發表偏倚評估，若P>0.05提示不存在發表偏倚。

結 果

1.文獻檢索結果：根據檢索策略初步檢出文獻91篇，經過嚴格篩選最終共納入10篇文獻，納入文獻均為病例對照研究，共包括DN患者1759例，單純DM患者1525例，健康對照者2209例。納入文獻中有6項研究來自亞洲人群，其余4項研究分別來自北美洲、歐洲、南美洲及非洲。其中有6篇文獻選擇DN患者為實驗組、單純DM患者為對照組；9篇文獻將DN患者與健康人群進行比較。納入文獻使用NOS標準進行質量評分，評分均在6～7分，納入文獻質量中等(表1、表2)。

表1 IL-10基因-1082A/G多態性在實驗組及對照組(DM)分布情況

表2 IL-10基因-1082A/G多態性在實驗組及對照組(健康人群)分布情況

2.Meta分析結果：與DM組比較，DN組IL-10基因-1082A/G多態性在5種遺傳模型中比較，差異無統計學意義：等位基因模型(G vs A)：OR=0.907，95% CI：0.673～1.224，P=0.524；純合子模型(GG vs AA)：OR=0.806，95% CI:0.379～1.714，P=0.576； 顯性遺傳模型(GG+GA vs AA)：OR=0.917，95% CI:0.588～1.430，P=0.702；隱性遺傳模型(GG vs GA+AA)：OR=0.888，95% CI:0.529～1.488，P=0.651； 雜合子模型(GA vs AA)：OR=1.110，95% CI:0.877～1.407，P=0.385(表3)。

與健康人群比較，DN患者在純合子模型下發現IL-10基因-1082A/G多態性與DN呈明顯相關性，純合子模型(GG vs AA)：OR=0.577, 95% CI:0.351～0.948，P<0.05，其他遺傳模型下差異無統計學意義。等位基因模型(G vs A)：OR=0.735，95% CI：0.487～1.110，P=0.143；顯性遺傳模型(GG+GA vs AA)：OR=0.761，95% CI：0.527～1.099，P=0.146；隱性遺傳模型(GG vs GA+AA)：OR=0.763，95% CI：0.566～1.030，P=0.077； 雜合子模型(GA vs AA)：OR=0.811，95% CI：0.563～1.169，P=0.262。以純合子模型GG vs AA為例，森林圖詳見圖1。納入文獻有6項研究來自亞洲，另外3項研究分別來自歐洲、北美洲、非洲。按地區進行亞組分析，結果顯示除雜合子模型外，IL-10基因-1082A/G多態性在其他遺傳模型上比較差異均有統計學意義。等位基因模型(G vs A)：OR=0.589, 95% CI：0.363～0.953，P=0.031；純合子模型(GG vs AA)：OR=0.349，95% CI：0.284～0.547，P=0.000；顯性模型(GG vs GA+AA)：OR=0.588，95% CI：0.382～0.904，P=0.015；隱性模型(GG+GA vs AA)：OR=0.629, 95% CI：0.513～0.771，P=0.000(表3)。

表3 IL-10基因-1082A/G多態性與DN易感性的Meta分析結果

3.敏感度分析：剔除不符合哈迪溫伯格定律的研究，再次計算OR及95% CI，得出結果較前相似，提示該研究結果穩定性較強(表4)。本研究應用Begger檢驗進行發表偏倚分析，結果顯示在5種遺傳模型中P均>0.05，提示該研究不存在發表偏倚。

表4 敏感度分析[OR(95% CI)]

討 論

近年來研究顯示，IL-10基因-1082A/G多態性與DN發生率密切相關。Fathy等[4]研究發現，DN組中-1082A/G多態位點A等位基因頻率明顯高于糖尿病非腎病組。Polina等[10]在瑞士人群中研究發現，G等位基因及GG、GA、AA基因型分布在DN、DM兩組間比較，差異無統計學意義。馬東紅等[14]研究發現，GG基因型與G等位基因在DN組分布顯著低于中國健康人群。Yin等[13]研究發現將DN與健康人群比較，攜帶AA基因型與GG基因型比較DN發生率增高；攜帶AA+GA基因型罹患DN的風險是攜帶GG基因型的2.38倍；攜帶AA基因型發生DN的風險是攜帶GG+GA基因型的2.08倍。Ma等[12]研究發現，攜帶A等位基因可顯著增加中國人群罹患DN的風險。Mtiraoui等[9]研究發現攜帶-1082G/-819T/-592A(GTA)或-1082G/-819T/-592C (GTC) 等位基因的人群更不易罹患DN。但Erdogan等[7]在基于土耳其地區的研究中發現，無論在健康人群還是DM群體中，-1082A/G多態性與DN發病均無明顯關聯。

本研究通過Meta分析發現等位基因分布及基因型分布，DN組與DM組比較，差異無統計學意義。在健康人群中，攜帶GG基因型的人群發生DN的風險是AA基因型人群的0.58倍。對地區進行亞組分析發現，亞洲人群中攜帶G等位基因、GG基因型罹患DN的風險明顯降低。該結果提示IL-10基因-1082A/G多態性與DN發生率存在相關性，亞洲人群中G等位基因、GG基因型是DN發病的保護因素。在合并效應量時，筆者研究發現5種遺傳模型上均存在一定的異質性，故嘗試將地區進一步亞組分析，結果發現亞組分析后各遺傳模型上異質性明顯降低，提示本研究異質性主要來源于不同地區。

2018年Naing等[16]對IL-10-1082A/G多態性與DN發生率進行Meta分析，共納入病例對照研究6項，對照組為DM患者或健康人群，該研究未將對照組DM或健康人群進一步分組，且僅在顯性及隱性2種遺傳模型上進行統計分析，結果發現在顯性及隱性模型上DN與對照組比較，差異無統計學意義；將地區進行亞組分析，發現亞洲人群僅在隱性模型上差異有統計學意義。而本研究納入文獻更多更全(n=10)，將對照組中DN或健康對照人群進一步細分，最終在5種遺傳模型上進行統計分析。

IL-10在DN發生、發展中起著一定作用，但其機制仍不明確。有研究表明，在腎臟中IL-10可增加多種生長因子、細胞因子和趨化因子的合成和分泌，誘導系膜細胞過度增殖。IL-10還可誘導系膜細胞外基質成分分泌增加，導致其積聚。與此同時，IL-10可誘導胱抑素C和轉化生長因子-β(transforming growth factor-β, TGF-β)的合成，胱抑素C可調節組織炎癥和誘導系膜細胞增殖；TGF-β與IL-10可同時促進腎小球硬化和間質纖維化的發生。此外，TGF-β誘導足細胞上皮向間充質轉化，導致蛋白尿的出現，蛋白尿的發生進一步加重了腎小球硬化和間質纖維化，最終導致ESRD。

本研究存在一定的局限性：(1)目前IL-10基因-1082A/G多態性及DN發病之間的文獻國內外仍較少，本研究共納入10篇病例對照研究，文章數目偏少，可能會降低研究的準確性。(2)對照組人群均來源于醫院，存在一定的選擇偏倚。(3)本研究按照地區進行亞組分析，亞洲以外的地區研究相對較少，且納入研究均為2型DM患者，缺乏其他類型DM相關文獻。(4)本研究探討-1082A/G多態性與DN發生率之間的關聯，未考慮基因與環境或基因與基因之間的相互作用。因此仍需開展大樣本、基于不同地區、不同種族、不同DM類型的相關研究，進一步明確IL-10基因-1082A/G多態位點與DN發生率之間的關聯，這對于DN的早期篩查、診斷、治療及預后均具有重要的指導意義。