中國人群miR-146a 基因rs2910164多態(tài)性與缺血性腦卒中易感性的meta 分析

黃甫高，王春芳，韋 穎，雷 茗，譚 曇，羅宏成

右江民族醫(yī)學(xué)院附屬醫(yī)院1超聲科，2檢驗科，廣西 百色 533000

中國人群miR-146a 基因rs2910164多態(tài)性與缺血性腦卒中易感性的meta 分析

黃甫高1，王春芳2，韋 穎2，雷 茗2，譚 曇2，羅宏成2

右江民族醫(yī)學(xué)院附屬醫(yī)院1超聲科，2檢驗科，廣西 百色 533000

目的探討中國人群miR-146aC＞G（rs2910164）位點單核苷酸多態(tài)性與缺血性腦卒中遺傳易感性。方法檢索2015年12月前中國生物醫(yī)學(xué)文獻數(shù)據(jù)庫（CBM）、中國期刊全文數(shù)據(jù)庫（CNKI）、維普中文科技期刊數(shù)據(jù)庫（VIP）、萬方數(shù)據(jù)庫及PubMed、Ovid、Cochrane Library、EMBASE等數(shù)據(jù)庫，搜集有關(guān)中國人群miR-146aC＞G（rs2910164）位點多態(tài)性與缺血性腦卒中（IS）研究，采用NOS工具評價納入相關(guān)性研究的質(zhì)量，提取科學(xué)合理有效數(shù)據(jù)，采用Review Manager 5.0軟件進行Meta分析。結(jié)果共納入6篇文獻，共有IS患者1945例，健康對照者2456例，Meta分析尚未發(fā)現(xiàn)miR146a rs2910164 G/C基因多態(tài)性與缺血性腦卒中易感性具有相關(guān)性：基因型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。結(jié)論本研究尚未發(fā)現(xiàn)中國人群miR-146aC＞G（rs2910164）位點多態(tài)性與缺血性腦卒中易感性之間具有相關(guān)性。

缺血性腦卒中；miR-146a；基因多態(tài)性；Meta分析；中國人群

缺血性腦卒中（IS）是一種高發(fā)病率、高致殘率及高死亡率的腦血管疾病。在中國，每年每10萬人就有120～180人死于IS，在美國，IS已超越心臟病成為致死的第3大死因[1-2]。研究表明，遺傳因素在IS的發(fā)病機制中發(fā)揮重要作用。不同基因型的個體，患IS的風(fēng)險完全不同。微小RNA是一類長約22核苷酸非編碼小分子RNA，通過結(jié)合到目的基因信使RNA（mRNA）3’非編碼端，在轉(zhuǎn)錄后水平調(diào)節(jié)基因的表達，從而參與疾病的發(fā)生、發(fā)展和轉(zhuǎn)歸[3-4]。最近有一些研究探討了前體miR-146a的單核苷酸多態(tài)性（SNP）位點C＞G（rs2910164）與中國人群IS的關(guān)系，但不同地區(qū)、不同人群結(jié)果存在較大差異和不確定性[5-10]。因此，本文檢索并系統(tǒng)回顧了有關(guān)中國人群miR-146aC＞G（rs2910164）位點SNP與IS的國內(nèi)外文獻報道，采用Meta分析的方式探討二者關(guān)系。

1 資料與方法

1.1 一般資料

文獻入選標(biāo)準(zhǔn)：（1）以中國人群為研究對象，民族不限；（2）國內(nèi)外公開發(fā)表有關(guān)miR-146aC＞G（rs 2910164）位點多態(tài)性與IS的風(fēng)險的病例-對照研究，限于中英文；（3）研究方法、診斷標(biāo)準(zhǔn)及對照人群選擇標(biāo)準(zhǔn)相似；（4）對照組所研究對象標(biāo)本來源符合溫伯格（HWE）定律；（5）統(tǒng)計數(shù)據(jù)合理有效、統(tǒng)計學(xué)方法選擇恰當(dāng)、實驗原理和技術(shù)路線設(shè)計合理，統(tǒng)計結(jié)果表達準(zhǔn)確可靠，數(shù)據(jù)來源相同的同一地方人群，選擇樣本量最大的研究文獻。排除標(biāo)準(zhǔn)：數(shù)據(jù)來源重復(fù)不可靠或無法獲取數(shù)據(jù)的文獻、實驗設(shè)計缺乏科學(xué)性。

1.2 方法

1.2.1 文獻檢索 以miR-146a、基因多態(tài)性、腦卒中、IS、SNP等為關(guān)鍵檢索詞，從檢索中文數(shù)據(jù)庫有：萬方數(shù)據(jù)庫、維普中文科技期刊數(shù)據(jù)庫（VIP）、中國生物醫(yī)學(xué)文獻數(shù)據(jù)庫（CBM）、中國期刊全文數(shù)據(jù)庫（CNKI）；而檢索的英文數(shù)據(jù)庫包括：PubMed、Ovid、Cochrane Library、EMBASE等數(shù)據(jù)庫。檢索截止日期為2015年12月31日。說明是否是由2個系統(tǒng)評價員進行文獻檢索的。

1.2.2 入選文獻資料篩選、質(zhì)量評價標(biāo)準(zhǔn) 由2位獨立研究者根據(jù)以下納入與排除標(biāo)準(zhǔn)提取入選資料：（1）樣本量是否充足；（2）統(tǒng)計學(xué)設(shè)計及方法選擇是否合理及準(zhǔn)確；（3）疾病診斷標(biāo)準(zhǔn)是否明確；（4）等位基因及基因分型方法是否準(zhǔn)確可靠；（5）運用H-W定律檢測對照組基因分布是否符合要求（6）缺乏完整數(shù)據(jù)的資料由兩名不同研究者分別提取后統(tǒng)計分析比較，若存在數(shù)據(jù)納入標(biāo)準(zhǔn)異議，由第3位研究者參與共同討論決定。采用Newcastle-Ottawa Scale（NOS）系統(tǒng)對入選資料進行質(zhì)量評價。

1.3 統(tǒng)計學(xué)處理

選擇Revman 5.3進行統(tǒng)計處理。H-W遺傳平衡檢驗:H-W遺傳平衡檢測各入選研究資料中對照組基因型分布是否符合要求。miR146a rs2910164基因多態(tài)性與IS之間的相關(guān)性，以O(shè)R值及95%CI表示，各研究資料的異質(zhì)性以Q檢驗和I2檢驗表示，顯著性異質(zhì)性用隨機效應(yīng)模型進行數(shù)據(jù)合并，反之則采用固定效應(yīng)模型。采用逐一排除研究的方法進行敏感性分析，Egger’檢驗和Begg檢驗評估文獻是否存在發(fā)表偏移。

2 結(jié)果

2.1 納入文獻基本資料

中英文數(shù)據(jù)庫共獲取文獻36篇，排除研究基因非rs2910164（C＞G）位點及數(shù)據(jù)不全文獻，篩選到符合標(biāo)準(zhǔn)文獻6篇，最終入選文獻研究結(jié)果都符合HWE定律。納入研究結(jié)果NOS評分均＞6分，質(zhì)量評價結(jié)果表明，入選6組研究疾病診斷標(biāo)準(zhǔn)恰當(dāng)明確，對照組與病例組具有可比性，統(tǒng)計學(xué)設(shè)計及方法選擇合理和準(zhǔn)確，基因及基因型檢測方法正確；結(jié)果數(shù)據(jù)準(zhǔn)確可靠。入選6組中國人群miR-146aC＞G（rs2910164）位點多態(tài)性與IS易感性相關(guān)病例對照資料中，健康對照組2456例，患者1945例；文獻納入基本情況基因及基因型分布情況見表1，質(zhì)量評價結(jié)果見表2。

表1 入選文獻基本情況及miR-146aC＞G（rs2910164）基因型分布及等位基因頻率

2.2 Meta分析結(jié)果

2.2.1 G等位基因vsC等位基因與缺血性腦卒中的關(guān)系 G等位基因為暴露因素，C等位基因為非暴露因素，各研究異質(zhì)性檢驗結(jié)果有統(tǒng)計學(xué)意義（χ2=13.65，P=0.02），故采用隨機效應(yīng)模型，Meta分析結(jié)果顯示：中國人群中具有G等位基因與C等位基因人群缺血性腦卒中發(fā)病風(fēng)險的差異無統(tǒng)計學(xué)意義，OR=0.97，95%CI（0.89～1.06），P=0.056（圖1）。

表2 入選資料研究方法學(xué)質(zhì)量評價（分）

圖1 中國人群miR-146a（rs2910164）位點G等位基因 vs C等位基因與缺血性腦卒中相關(guān)性的Meta分析

2.2.2 基因型GGvs基因型CC與缺血性腦卒中的關(guān)系 GG基因型為暴露因素，CC基因型為非暴露因素，各研究異質(zhì)性檢驗結(jié)果有統(tǒng)計學(xué)意義（χ2=19.12，P=0.002），故采用隨機效應(yīng)模型，Meta分析結(jié)果顯示：中國人群具有GG基因型與CC基因型人群缺血性腦卒中發(fā)病風(fēng)險無統(tǒng)計學(xué)差異，OR=0.99，95%CI（0.82～1.18），P=0.88（圖2）。

2.2.3 基因型GCvs基因型CC與缺血性腦卒中的關(guān)系 GC基因型為暴露因素，CC基因型為非暴露因素，各研究異質(zhì)性檢驗結(jié)果有統(tǒng)計學(xué)意義（χ2=5.76，P=0.33），故采用固定效應(yīng)模型，Meta分析結(jié)果顯示：中國人群具有GC基因型與具有CC基因型人群缺血性腦卒中發(fā)病風(fēng)險差異無統(tǒng)計學(xué)意義，OR=1.02，95%CI（0.90～1.17），P=0.75（圖3）。

圖2 中國人群miR-146a（rs2910164）位點GG基因型 vs CC基因型與缺血性腦卒中相關(guān)性的Meta分析

圖3 中國人群miR-146a（rs2910164）位點GC基因型 vs CC基因型與缺血性腦卒中相關(guān)性的Meta分析

2.2.4 基因型GG+GC vs基因型CC與缺血性腦卒中的關(guān)系 GG+GC基因型為暴露因素，CC基因型為非暴露因素，各研究異質(zhì)性檢驗結(jié)果有統(tǒng)計學(xué)意義（χ2=10.98，P=0.05），故采用固定效應(yīng)模型，Meta分析結(jié)果顯示：中國人群具有GG和GC基因型與具有CC基因型人群缺血性腦卒中發(fā)病風(fēng)險無統(tǒng)計學(xué)差異，OR=1.01，95%CI（0.90～1.15），P=0.82（圖4）。

圖4 中國人群miR-146a（rs2910164）位點GC+GG基因型 vs CC基因型與缺血性腦卒中相關(guān)性的Meta分析

2.3 敏感性分析

采用逐一排除各研究的方法進行敏感性分析，重新估計合并效應(yīng)量，并與排除前的合并效應(yīng)量進行分析比較，發(fā)現(xiàn)效應(yīng)量無明顯改變，表明此次研究結(jié)果穩(wěn)定可靠。

2.4 發(fā)表偏倚分析

以O(shè)R值為橫坐標(biāo)，log（RR）的標(biāo)準(zhǔn)差為縱坐標(biāo)繪制漏斗圖，從漏斗圖可以看出，4種遺傳學(xué)模型研究基本對稱（圖5），采用Egger′s檢驗進一步分析，結(jié)果分別為GvsC：P＞|t|=0.167，GGvsCC：P＞|t|=0.576，GCvsCC：P＞|t|=0.312，GG+GCvsCC：P＞|t|=0.255，結(jié)果表明無明顯發(fā)表偏倚分析存在，但本研究符合納入標(biāo)準(zhǔn)文獻少，失安全系數(shù)小，會存在一定的發(fā)表偏倚。

3 討論

IS是當(dāng)今嚴(yán)重威脅人類生命的一類常見的腦血管疾病，由遺傳和環(huán)境等多種因素共同作用而導(dǎo)致的一類極其復(fù)雜的疾病[11-15]。傳統(tǒng)觀點認為動脈粥樣硬化、高血壓、血脂異常、高同型半胱氨酸血癥等是IS的危險因子。然而，這些危險因子僅能解釋部分IS的病因，表明除了環(huán)境因素外，遺傳因素在IS發(fā)病中發(fā)揮著不可忽視的作用。SNP主要由于單個堿基改變而引起分子水平上DNA表達序列多態(tài)性，是人類可遺傳變異中最常見的一種，人類個體對疾病的易感性及對藥物的反應(yīng)性差異與SNP密切相關(guān)。人體基因組大約30%轉(zhuǎn)錄體受miRNAs的調(diào)控[16]，miRNA參與細胞的分化、增殖、轉(zhuǎn)化等生理病理過程[17]。MiR-146a位于人類基因組的5q33，存在于前體miR-146a（pre-miRNA ）G＞C SNP位點通過調(diào)控成熟miR-146a的表達水平，進而參與IS的病理過程[5-10]。近年來有關(guān)miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性的單個研究結(jié)果不盡一致，且同一種族不人群存在不確定性。為此，本研究就中國人群miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性進行Meta分析，旨在揭示miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性之間的相關(guān)性，為遺傳學(xué)上對IS的分子診治具有重要意義。

圖5 4種遺傳學(xué)模型發(fā)表偏倚漏斗圖A: G vs C; B: GG vs CC; C: GC vs CC; D: GG+GC vs CC.

本文首次搜集中國人群miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性之間病例對照研究資料。通過研究尚未發(fā)現(xiàn)中國人群miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性之間存在相關(guān)性，國內(nèi)外學(xué)者有關(guān)miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性之間的研究結(jié)果不盡相同。國內(nèi)研究結(jié)果[5-9]認為致中國人群miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性之間存在相關(guān)性。而Liu[10]研究認為中國人群miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性沒有相關(guān)性。國內(nèi)外學(xué)者[18-21]研究認為miR-146aC＞G（rs2910164）可能增加缺血性腦卒中的發(fā)病風(fēng)險。研究結(jié)果的差異可能與人類個體種族地域差異、物理、化學(xué)環(huán)境、生物因素誘導(dǎo)作用差異相關(guān)，也可能由于納入研究對象標(biāo)準(zhǔn)及樣本量大小不同的原因。

本研究納入文獻僅限于中英文符合要求的文獻數(shù)量偏少；研究人群僅限于中國少數(shù)地區(qū)人群，沒有體現(xiàn)群體研究的廣泛性和代表性；病例對照研究組性別、年齡、民族等未能按要求完全匹配。中國人群miR-146aC＞G（rs2910164）位點基因多態(tài)性與IS易感性之間相關(guān)性更合理、可靠的結(jié)論有賴于開展大樣本、多中心、同質(zhì)性的病例對照研究。

[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基因多態(tài)性與長沙地區(qū)漢族人群缺血性腦卒中遺傳易感性的關(guān)聯(lián)研究[D]. 長沙: 中南大學(xué), 2011.

[6]李 友, 朱 菁, 馬國達, 等. miR-146a基因多態(tài)性與腔隙性腦梗死的關(guān)系[J]. 山東醫(yī)藥, 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基因多態(tài)性與缺血性腦卒中易感性的研究[J]. 現(xiàn)代生物醫(yī)學(xué)進展,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與缺血性腦卒中的關(guān)系[J]. 新醫(yī)學(xué), 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

黃甫高，主治醫(yī)師，E-mail: 513223727@qq.com

羅宏成，碩士，主管技師，E-mail: weizhuren999@163.com