中国人群miR-146a 基因rs2910164多态性与缺血性脑卒中易感性的meta 分析

黄甫高，王春芳，韦 颖，雷 茗，谭 昙，罗宏成

右江民族医学院附属医院1超声科，2检验科，广西 百色 533000

中国人群miR-146a 基因rs2910164多态性与缺血性脑卒中易感性的meta 分析

黄甫高1，王春芳2，韦 颖2，雷 茗2，谭 昙2，罗宏成2

右江民族医学院附属医院1超声科，2检验科，广西 百色 533000

目的探讨中国人群miR-146aC＞G（rs2910164）位点单核苷酸多态性与缺血性脑卒中遗传易感性。方法检索2015年12月前中国生物医学文献数据库（CBM）、中国期刊全文数据库（CNKI）、维普中文科技期刊数据库（VIP）、万方数据库及PubMed、Ovid、Cochrane Library、EMBASE等数据库，搜集有关中国人群miR-146aC＞G（rs2910164）位点多态性与缺血性脑卒中（IS）研究，采用NOS工具评价纳入相关性研究的质量，提取科学合理有效数据，采用Review Manager 5.0软件进行Meta分析。结果共纳入6篇文献，共有IS患者1945例，健康对照者2456例，Meta分析尚未发现miR146a rs2910164 G/C基因多态性与缺血性脑卒中易感性具有相关性：基因型GvsC：OR=0.97，95%CI（0.89～1.06），P=0.06；基因型GGvsCC：OR=0.99，95%CI（0.82～1.18），P=0.88；基因型GCvsCC：OR=1.02，95%CI（0.90～1.17），P=0.75；基因型GG+GCvsCC：OR=1.01，95%CI（0.90～1.15），P=0.82。结论本研究尚未发现中国人群miR-146aC＞G（rs2910164）位点多态性与缺血性脑卒中易感性之间具有相关性。

缺血性脑卒中；miR-146a；基因多态性；Meta分析；中国人群

缺血性脑卒中（IS）是一种高发病率、高致残率及高死亡率的脑血管疾病。在中国，每年每10万人就有120～180人死于IS，在美国，IS已超越心脏病成为致死的第3大死因[1-2]。研究表明，遗传因素在IS的发病机制中发挥重要作用。不同基因型的个体，患IS的风险完全不同。微小RNA是一类长约22核苷酸非编码小分子RNA，通过结合到目的基因信使RNA（mRNA）3’非编码端，在转录后水平调节基因的表达，从而参与疾病的发生、发展和转归[3-4]。最近有一些研究探讨了前体miR-146a的单核苷酸多态性（SNP）位点C＞G（rs2910164）与中国人群IS的关系，但不同地区、不同人群结果存在较大差异和不确定性[5-10]。因此，本文检索并系统回顾了有关中国人群miR-146aC＞G（rs2910164）位点SNP与IS的国内外文献报道，采用Meta分析的方式探讨二者关系。

1 资料与方法

1.1 一般资料

文献入选标准：（1）以中国人群为研究对象，民族不限；（2）国内外公开发表有关miR-146aC＞G（rs 2910164）位点多态性与IS的风险的病例-对照研究，限于中英文；（3）研究方法、诊断标准及对照人群选择标准相似；（4）对照组所研究对象标本来源符合温伯格（HWE）定律；（5）统计数据合理有效、统计学方法选择恰当、实验原理和技术路线设计合理，统计结果表达准确可靠，数据来源相同的同一地方人群，选择样本量最大的研究文献。排除标准：数据来源重复不可靠或无法获取数据的文献、实验设计缺乏科学性。

1.2 方法

1.2.1 文献检索 以miR-146a、基因多态性、脑卒中、IS、SNP等为关键检索词，从检索中文数据库有：万方数据库、维普中文科技期刊数据库（VIP）、中国生物医学文献数据库（CBM）、中国期刊全文数据库（CNKI）；而检索的英文数据库包括：PubMed、Ovid、Cochrane Library、EMBASE等数据库。检索截止日期为2015年12月31日。说明是否是由2个系统评价员进行文献检索的。

1.2.2 入选文献资料筛选、质量评价标准 由2位独立研究者根据以下纳入与排除标准提取入选资料：（1）样本量是否充足；（2）统计学设计及方法选择是否合理及准确；（3）疾病诊断标准是否明确；（4）等位基因及基因分型方法是否准确可靠；（5）运用H-W定律检测对照组基因分布是否符合要求（6）缺乏完整数据的资料由两名不同研究者分别提取后统计分析比较，若存在数据纳入标准异议，由第3位研究者参与共同讨论决定。采用Newcastle-Ottawa Scale（NOS）系统对入选资料进行质量评价。

1.3 统计学处理

选择Revman 5.3进行统计处理。H-W遗传平衡检验:H-W遗传平衡检测各入选研究资料中对照组基因型分布是否符合要求。miR146a rs2910164基因多态性与IS之间的相关性，以OR值及95%CI表示，各研究资料的异质性以Q检验和I2检验表示，显著性异质性用随机效应模型进行数据合并，反之则采用固定效应模型。采用逐一排除研究的方法进行敏感性分析，Egger’检验和Begg检验评估文献是否存在发表偏移。

2 结果

2.1 纳入文献基本资料

中英文数据库共获取文献36篇，排除研究基因非rs2910164（C＞G）位点及数据不全文献，筛选到符合标准文献6篇，最终入选文献研究结果都符合HWE定律。纳入研究结果NOS评分均＞6分，质量评价结果表明，入选6组研究疾病诊断标准恰当明确，对照组与病例组具有可比性，统计学设计及方法选择合理和准确，基因及基因型检测方法正确；结果数据准确可靠。入选6组中国人群miR-146aC＞G（rs2910164）位点多态性与IS易感性相关病例对照资料中，健康对照组2456例，患者1945例；文献纳入基本情况基因及基因型分布情况见表1，质量评价结果见表2。

表1 入选文献基本情况及miR-146aC＞G（rs2910164）基因型分布及等位基因频率

2.2 Meta分析结果

2.2.1 G等位基因vsC等位基因与缺血性脑卒中的关系 G等位基因为暴露因素，C等位基因为非暴露因素，各研究异质性检验结果有统计学意义（χ2=13.65，P=0.02），故采用随机效应模型，Meta分析结果显示：中国人群中具有G等位基因与C等位基因人群缺血性脑卒中发病风险的差异无统计学意义，OR=0.97，95%CI（0.89～1.06），P=0.056（图1）。

表2 入选资料研究方法学质量评价（分）

图1 中国人群miR-146a（rs2910164）位点G等位基因 vs C等位基因与缺血性脑卒中相关性的Meta分析

2.2.2 基因型GGvs基因型CC与缺血性脑卒中的关系 GG基因型为暴露因素，CC基因型为非暴露因素，各研究异质性检验结果有统计学意义（χ2=19.12，P=0.002），故采用随机效应模型，Meta分析结果显示：中国人群具有GG基因型与CC基因型人群缺血性脑卒中发病风险无统计学差异，OR=0.99，95%CI（0.82～1.18），P=0.88（图2）。

2.2.3 基因型GCvs基因型CC与缺血性脑卒中的关系 GC基因型为暴露因素，CC基因型为非暴露因素，各研究异质性检验结果有统计学意义（χ2=5.76，P=0.33），故采用固定效应模型，Meta分析结果显示：中国人群具有GC基因型与具有CC基因型人群缺血性脑卒中发病风险差异无统计学意义，OR=1.02，95%CI（0.90～1.17），P=0.75（图3）。

图2 中国人群miR-146a（rs2910164）位点GG基因型 vs CC基因型与缺血性脑卒中相关性的Meta分析

图3 中国人群miR-146a（rs2910164）位点GC基因型 vs CC基因型与缺血性脑卒中相关性的Meta分析

2.2.4 基因型GG+GC vs基因型CC与缺血性脑卒中的关系 GG+GC基因型为暴露因素，CC基因型为非暴露因素，各研究异质性检验结果有统计学意义（χ2=10.98，P=0.05），故采用固定效应模型，Meta分析结果显示：中国人群具有GG和GC基因型与具有CC基因型人群缺血性脑卒中发病风险无统计学差异，OR=1.01，95%CI（0.90～1.15），P=0.82（图4）。

图4 中国人群miR-146a（rs2910164）位点GC+GG基因型 vs CC基因型与缺血性脑卒中相关性的Meta分析

2.3 敏感性分析

采用逐一排除各研究的方法进行敏感性分析，重新估计合并效应量，并与排除前的合并效应量进行分析比较，发现效应量无明显改变，表明此次研究结果稳定可靠。

2.4 发表偏倚分析

以OR值为横坐标，log（RR）的标准差为纵坐标绘制漏斗图，从漏斗图可以看出，4种遗传学模型研究基本对称（图5），采用Egger′s检验进一步分析，结果分别为GvsC：P＞|t|=0.167，GGvsCC：P＞|t|=0.576，GCvsCC：P＞|t|=0.312，GG+GCvsCC：P＞|t|=0.255，结果表明无明显发表偏倚分析存在，但本研究符合纳入标准文献少，失安全系数小，会存在一定的发表偏倚。

3 讨论

IS是当今严重威胁人类生命的一类常见的脑血管疾病，由遗传和环境等多种因素共同作用而导致的一类极其复杂的疾病[11-15]。传统观点认为动脉粥样硬化、高血压、血脂异常、高同型半胱氨酸血症等是IS的危险因子。然而，这些危险因子仅能解释部分IS的病因，表明除了环境因素外，遗传因素在IS发病中发挥着不可忽视的作用。SNP主要由于单个碱基改变而引起分子水平上DNA表达序列多态性，是人类可遗传变异中最常见的一种，人类个体对疾病的易感性及对药物的反应性差异与SNP密切相关。人体基因组大约30%转录体受miRNAs的调控[16]，miRNA参与细胞的分化、增殖、转化等生理病理过程[17]。MiR-146a位于人类基因组的5q33，存在于前体miR-146a（pre-miRNA ）G＞C SNP位点通过调控成熟miR-146a的表达水平，进而参与IS的病理过程[5-10]。近年来有关miR-146aC＞G（rs2910164）位点基因多态性与IS易感性的单个研究结果不尽一致，且同一种族不人群存在不确定性。为此，本研究就中国人群miR-146aC＞G（rs2910164）位点基因多态性与IS易感性进行Meta分析，旨在揭示miR-146aC＞G（rs2910164）位点基因多态性与IS易感性之间的相关性，为遗传学上对IS的分子诊治具有重要意义。

图5 4种遗传学模型发表偏倚漏斗图A: G vs C; B: GG vs CC; C: GC vs CC; D: GG+GC vs CC.

本文首次搜集中国人群miR-146aC＞G（rs2910164）位点基因多态性与IS易感性之间病例对照研究资料。通过研究尚未发现中国人群miR-146aC＞G（rs2910164）位点基因多态性与IS易感性之间存在相关性，国内外学者有关miR-146aC＞G（rs2910164）位点基因多态性与IS易感性之间的研究结果不尽相同。国内研究结果[5-9]认为致中国人群miR-146aC＞G（rs2910164）位点基因多态性与IS易感性之间存在相关性。而Liu[10]研究认为中国人群miR-146aC＞G（rs2910164）位点基因多态性与IS易感性没有相关性。国内外学者[18-21]研究认为miR-146aC＞G（rs2910164）可能增加缺血性脑卒中的发病风险。研究结果的差异可能与人类个体种族地域差异、物理、化学环境、生物因素诱导作用差异相关，也可能由于纳入研究对象标准及样本量大小不同的原因。

本研究纳入文献仅限于中英文符合要求的文献数量偏少；研究人群仅限于中国少数地区人群，没有体现群体研究的广泛性和代表性；病例对照研究组性别、年龄、民族等未能按要求完全匹配。中国人群miR-146aC＞G（rs2910164）位点基因多态性与IS易感性之间相关性更合理、可靠的结论有赖于开展大样本、多中心、同质性的病例对照研究。

[1]Liu L, Wang D, Wong K. Stroke and stroke care in China: Huge burden, significant workload and a national priority[J]. Stroke,2011, 42(5): 3651-4.

[2]Goldstein LB, Adams R, Becker K, et al. Primary prevention of ischemic stroke: A statement for healthcare professionals from the Stroke Council of the American Heart Association[J]. Circulation,2001, 103(1): 163-82.

[3]Bartel DP. MicroRNAs: genomics, biogenesis, mechanism and function[J]. Cell, 2004, 116(2): 281-97.

[4]Pushparaj PN, Aarthi JJ, Kumar SD, et al. RNAi and RNAa--the yin and yang of RNAome[J]. Bioinformatics, 2008, 2(6): 235-7.

[5]孙 吉. miRNA-146a和EPHX2基因多态性与长沙地区汉族人群缺血性脑卒中遗传易感性的关联研究[D]. 长沙: 中南大学, 2011.

[6]李 友, 朱 菁, 马国达, 等. miR-146a基因多态性与腔隙性脑梗死的关系[J]. 山东医药, 2014, 54(16): 1-3.

[7]Zhu R, Liu X, He Z, et al. miR-146a and miR-196a2 polymorphisms in patients with ischemic stroke in the northern Chinese Han population[J]. Neurochem Res, 2014, 39(9): 1709-16.

[8]胡亚梅, 李书剑, 姜晓峰, 等. miR-146aC＞G, miR-149T＞C基因多态性与缺血性脑卒中易感性的研究[J]. 现代生物医学进展,2014, 14(29): 5643-8.

[9]Huang S, Zhou S, Zhang Y, et al. Association of the genetic polymorphisms in pre-microRNAs with risk of ischemic stroke in a Chinese population[J]. PLoS One, 2015, 10(2): e0117007-12.

[10]Liu Y, Ma Y, Zhang B, et al. Genetic polymorphisms in premicroRNAs and risk of ischemic stroke in a Chinese population[J].J Mol Neurosci, 2014, 52(4): 473-80.

[11]Falcone GJ, Malik R, Dichgans M, et al. Current concepts and clinical applications of stroke genetics[J]. Lancet Neurology, 2014,13(4): 405-8.

[12]Meschia JF, Arnett DK, Ay H, et al. Stroke genetics network(SiGN) study: design and rationale for a genome-wide association study of ischemic stroke subtypes[J]. Stroke, 2013, 44(10): 2694-702.

[13]Sharma P, Yadav S, Meschia JF. Genetics of ischaemic stroke[J]. J Neurol Neurosurg Psychiatry, 2013, 84(12): 1302-8.

[14]Anderson V, Gomes A, Greenham M, et al. Social competence following pediatric stroke: contributions of brain insult and family environment[J]. Soc Neurosci, 2014, 9(5): 471-83.

[15]Hamano T, Kawakami N, Li X, et al. Neighbourhood environment and stroke: a follow-up study in Sweden[J]. PLoS One, 2013, 8(2):e56680-4.

[16]Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets[J]. Cell, 2005, 120(1): 15-20.

[17]Ambros V. The functions of animal microRNAs[J]. Nature, 2004,431(76): 350-5.

[18]吕国天, 王 爽, 王庆坤. 微小RNA-146a rs2910164及微小RNA-499 rs3746444与缺血性脑卒中的关系[J]. 新医学, 2016, 47(04):257-60.

[19]Zhu H, Zhang H, Bao L, et al. Analysis of association of genetic polymorphisms of microRNAs with ischemic stroke[J]. Chin J Med Genet, 2017, 34(2): 261-5.

[20]Du J, Cui C, Zhang S, et al. Association of MicroRNA-146a and MicroRNA-149 polymorphisms with strokes in Asian populations:an updated Meta-analysis[J]. Angiology, 2017, 68(10): 863-70.

[21]Su ZF, Sun ZW, Zhang Y, et al. Regulatory effects of miR-146a/b on the function of endothelial progenitor cells in acute ischemic stroke in mice[J]. J Med Sci, 2017, 33(8): 369-78.

Association between miR-146a (rs2910164) gene polymorphism and susceptibility to ischemic stroke in Chinese population: a meta-analysis

HUANG Fugao1, Wang Chunfang2, Wei Ying2, Lei Ming2, Tan Tan2, LUO Hongcheng2

1Department of Ultrasonography,2Department of Clinical Laboratory, The Affiliated Hospital of Youjiang Medical College for Nationalities,Baise 533000, China

ObjectiveTo investigate the association between miR-146a gene polymorphism and Ischemic stroke (IS) in Chinese population.MethodsPublications in PubMed, Ovid, Cochrane Library EMBASE, CNKI, VIP and Wangfang Datebase were retrieved up to December 2015. Studies of investigating association between miR-146aC＞G (rs2910164) and IS in Chinese population were included. Quality of inclubed studies were evaluated by a Newcastle-Ottawa Scale (NOS) tool.The valied data were extracted and analyzed by Review Manager5.0 software.ResultsSix studies involving 1945 patients of IS and 2456 controls were included. No association between miR-146a rs2910164 polymorphism and the risk of IS was observed [GvsC:OR=0.97, 95%CI: (0.89-1.06),P=0.06]; [GGvsCC: OR=0.99, 95%CI: (0.82-1.18),P=0.88]; [ GCvsCC: OR=1.02, 95%CI (0.90-1.17),P=0.75]; [GG+GCvsCC: OR=1.01, 95%CI (0.90-1.15),P=0.82].ConclusionThe Meta-analysis suggested that the rs2910164 polymorphism in miR-146a is not associated with IS in Chinese population.

ischemic stroke; miR-146a; gene polymorphism; Meta-analysis; Chinese

2017-07-06

黄甫高，主治医师，E-mail: 513223727@qq.com

罗宏成，硕士，主管技师，E-mail: weizhuren999@163.com