IL-10基因多态性与复发性流产相关性的Meta分析*

袁恩武，闫广伟，代延朋，邢金芳

(郑州大学第三附属医院：1.检验科;2.病理科， 郑州 450052)

复发性流产(RSA)是指一对夫妇连续发生两次或两次以上的自然流产[1]。RSA是妊娠常见并发症且病因学复杂[2-3]。白细胞介素(IL)-10主要由外周血辅助T细胞(Th)2细胞产生，其可以下调Th1细胞产生IL-12和干扰素(INF)-α等细胞因子来调节Th2细胞的免疫反应[4]。研究表明， IL-10 基因多态性直接影响体内 IL-10的表达水平，从而和疾病的易感、严重性及疾病的发展进程密切相关[5]。IL-10基因位于染色体1q31-32上，由5个外显子和 4个内含子组成[6]，IL-10在启动子区域存在3个单核苷酸多态性(SNP)：1082G>A(rs：1800896)、819C>T(rs：1800871)和592C>A(rs：1800872)。 虽然有研究显示IL-10基因多态性和RSA的发生有关，但不同的研究得出的结果并不一致，现有相关Meta分析文章也存在分歧。因此本研究通过检索国内外已发表的IL-10基因多态性与RSA相关文章，采用Meta分析方法进行综合评价，探索两者的相关性，为RSA的临床诊断和相关机制的研究提供依据。

1 材料与方法

1.1文献检索 检索PubMed、Embase、Web of Science、中国知网、万方知识服务平台和维普数据库，文献检索时间为建库至2016年10月1日。英文检索词：RSA,recurrent spontaneous abortion,recurrent miscarriage,recurrent pregnancy loss,IL-10,interleukin-10,gene,polymorphism,SNP；中文检索词为复发性流产、基因多态性、IL-10、白细胞介素-10。采用主题词结合自由词检索的方式，并人工检索纳入研究的参考文献。

1.2纳入和排除标准

1.2.1纳入标准 (1)入选文献必须是IL-10基因多态性与RSA相关研究且RSA定义明确(即两次或两次以上的流产)；(2)研究方案为病例对照研究，并且文献中明确给出病例组和对照组不同基因型的具体人数；(3)病例组为患有RSA的妇女，对照组为至少1次成功妊娠且无自然流产史的健康妇女，研究中病例组和健康对照来源于年龄、民族等一般信息均衡可比的统一人群。

1.2.2排除标准 (1)未给出RSA的明确定义；(2)数据不完整或不能计算出不同基因型的基因频数；(3)重复或有重叠数据文章。

1.3文献质量评价和资料提取 由两名评价员独立按照纳入和排除标准进行文献筛选，如遇分歧进行商讨或由第3名评价员裁决。对筛选出的文献提取主要内容，包括作者、发表时间、国家地区、病例组和对照组3种基因型的例数等。同时由两位评价员根据纽尔卡斯-渥太华(NOS)评分量表对纳入文献的质量进行评价，满分为9分，若总分大于或等于6分为高质量研究。评价标准包括病例组与对照组的选择(4分)，病例组与对照组的可比性(2分)，病例组与对照组暴露因素测量(3分)。

1.4统计学处理 采用RevMan 5.3软件进行Meta分析，选择等位基因模型、显性基因模型、隐性基因模型、纯合子模型4种遗传模型进行Meta分析。通过Q检验对纳入研究进行异质性检验，若各研究结果I2≤50%表明无明显异质性，采用固定效应模型进行分析，否则采用随机效应模型进行分析，并根据国家和地区进行亚组分析。以OR和95%CI作为统计指标，以P<0.05为差异有统计学意义。用χ2检验来分析对照组的Hardy-Weinberg平衡，通过绘制漏斗图评估发表偏倚。

2 结 果

2.1文献检索结果 初步筛选得到551篇相关文献，按照纳入和排除标准进一步阅览标题和摘要剩余49篇，详细查看全篇，排除缺失数据和综述类等不满足纳入标准的文献，最后纳入文献18篇[7-24],NOS评分均大于或等于6分，累计病例2 688例，对照2 522例，12篇文献研究-592A/C位点，8篇文献研究-819C/T位点，14篇文献研究-1082G/C位点。文献筛选流程见图1，纳入各项研究的相关信息见表1～3。

图1 文献筛选流程及结果

纳入文献年份国家地区病例组基因(n)AAACCC对照组基因(n)AAACCCNOS评分(分)HWEP值杨秋红等[7]2010年中国亚洲53642534853770.260ZAMMITI等[8]2006年突尼斯非洲5193206254113480.000

续表1 关于IL-10-592A/C位点的纳入文献的基本特征

HWE：Hard-Weinberg平衡，P<0.05表示对照组人群的基因型频率分布不符合HWE平衡

表2 关于IL-10-819C/T位点的纳入文献的基本特征

表3 关于IL-10-1082G/A位点的纳入文献的基本特征

2.2Meta分析结果

2.2.1IL-10-592A/C基因多态性和RSA相关性的Meta分析 共有12篇文献研究了IL-10-592A/C基因多态性和RSA发生风险的关系。异质性检验显示均有明显性，采用随机效应模型进行分析。在总人群中，结果显示IL-10-592A/C多态性与RSA无明显相关性(Avs.C:OR=1.02，95%CI=0.81～1.28，P=0.87；CC+CAvs.AA:OR=0.94，95%CI=0.67～1.32，P=0.74)。对不同地区人群亚组分析未见其相关性，见表4。

2.2.2IL-10-819C/T基因多态性和RSA相关性的Meta分析 共有8篇文献研究了IL-10-819C/T基因多态性和RSA发生风险的关系。异质性检验显示等位基因(Cvs.T)和显性基因模型(TT+TCvs.CC)存在较大异质性，采用随机效应模型进行分析；隐性基因模型(TTvs.TC+CC)和等位基因(TTvs.CC)无明显异质性，采用固定效应模型进行分析。在总体人群中，各项研究间等位基因Cvs.T 和TT+TCvs.CC基因差异无统计学意义(Tvs.C:OR=1.30，95%CI=0.99～1.70，P=0.06；TT+TCvs.CC:OR=1.40，95%CI=0.93～2.11，P=0.11)，但在隐性基因模型和纯合子模型下发现IL-10-819C/T基因多态性与RSA具有明显相关性(TTvs.TC+CC:OR=1.40，95%CI=1.10～1.78，P<0.01；TTvs.CC:OR=1.57，95%CI=1.21～2.04，P<0.05)。对不同地区人群亚组分析显示在亚洲人群中，在纯合子模型(TTvs.CC)下发现IL-10-819C/T基因多态性与RSA发生具有关(OR=1.61，95%CI=1.01～2.57，P=0.05)，其他结果均无统计学意义，见表5。

2.2.3IL-10-1082G/A基因多态性和RSA相关性的Meta分析 共有14篇文献研究了IL-10-1082G/A基因多态性和RSA发生风险的关系。异质性检验显示除隐性基因模型(GGvs.GA+AA)和亚洲亚组具有较大异质性，其余研究均无明显异质性，分别采用随机效应模型和固定效应模型进行分析。在总体人群中，等位基因G与RSA发生具有相关性，携带等位基因G较携带等位基因A的发病风险可能升高(OR=1.13，95%CI=1.03～1.24，P=0.01)，同时在纯合子模型下也显示出和RSA的相关性(OR=1.39，95%CI=1.14～1.71，P<0.05)。各项研究间GG+GAvs.CC基因和GGvs.GA+AA基因差异无统计学意义(P>0.05)。对不同地区人群亚组分析发现在欧洲人群中IL-10-1082G/A基因多态性和RSA发生具有相关性(Gvs.A:OR=1.31，95%CI=1.02～1.68，P=0.03；GGvs.AA:OR=1.77，95%CI=1.08～2.88，P=0.02)，见表6、图2。

表4 IL-10-592A/C基因多态性与RSA相关性

表5 IL-10-819C/T基因多态性与RSA相关性

表6 IL-10-1082G/A基因多态性与RSA相关性

图2 IL-10-1082G/A等位基因模型亚组分析结果

2.3敏感性和发表偏倚分析 HWE是保证SNP位点有效性的基本遗传定律，去掉不符合HWE的文献后分析结果基本没有发生变化，表明此分析结果相对稳定可靠。由于纳入文献过少不适合绘制漏斗图，因此仅对总体人群不同基因模型进行漏斗图分析，IL-10-819C/T在隐性模型下漏斗图不对称，可能是由于纳入了低质量小样本研究，删除3篇样本量小于100的文献，漏斗图显示对称，其余漏斗图散点分布均匀，基本对称，发表偏倚较小。

3 讨 论

随着中国“二胎”政策的实施，许多大龄女性(≥35岁)选择再生育，高龄孕妇也随之增多，而随着年龄增加，女性妊娠不良结局和并发症也会增加。高龄产妇发生流产的危险性较高,这是因为染色体异常胎儿的发生率增加,从而使流产发生率增加[25-26]。RSA是自然流产中常见的，也是目前人们研究较多的，本文采用Meta分析方法，分别采用等位基因模型，显性基因模型，隐性基因模型和纯合子模型对IL-10基因多态性和RSA相关性进行评价。

本研究结果显示，IL-10-592A/C基因多态性和RSA无直接相关性，亚组分析也未发现其相关性，但国内相关研究均显示IL-10-592A/C基因多态性可能与RSA发病有关，詹福寿等[24]研究发现IL-10-592A/C多态性位点可能是宁夏地区汉族女性复发性自然流产的易感SNP，但这可能和不同民族遗传差异有关,并且纳入的相关文献之间均存在较大异质性。对IL-10-819C/T位点的研究发现，基因型为CC者女性RSA发病风险高于基因型TT者。由于纳入文献亚洲地区占多数，笔者按地区进行亚组分析，亚洲人群分析结果和总体人群结果一致，而欧洲人群未发现这种相关性。IL-10-1082G/A与RSA相关性分析表明，等位基因G和GG基因型均与RSA发生相关，DAHER等[18]研究发现GG基因型在RSA病例组分布频率明显高于健康对照组；PENG等[27]对IL-10-1082G/A多态性位点与RSA相关性进行Meta分析发现，IL-1082G等位基因和GG基因型会增加RSA的发生风险，与本次研究结果一致。IL-10-819C/T和IL-10-1082G/A两个多态性位点与RSA之间的关系目前国内未见相关报道，还需进一步设计严谨、大样本量的试验对IL-10基因多位点检测进行确认。

尽管漏斗图显示发表偏倚较小，但本次Meta分析也存在一些局限性。(1)IL-10-592A/C多态性分析中纳入研究之间存在较大异质性，经过亚组分析异质性仍然存在，异质性可能和地区、样本量、检测方法等因素有关，较大的异质性在一定程度上影响IL-10-592A/C位点结果可靠性；(2)在文献搜索过程中不免遗漏个别相关文献，而且有的文献无法获得各个基因型相关数据被排除[5,22,23];(3)本次纳入文献大多数来源于亚洲地区，具有一定地区局限性。

综上所述，本研究结果显示，IL-10-819C/T、IL-10-1082G/A基因多态性和RSA的发病相关，IL-10-592A/C基因多态性和RSA无明显相关性，对今后研究基因与RSA发生的易感性提供理论基础，但本研究存在一定局限性，还需更多大样本同质性的研究证明。