中國(guó)人群 ALOX5AP基因 SG13S114A/T多態(tài)性與缺血性卒中相關(guān)性的 Meta分析

徐佳亮, 藺慕會(huì), 陳曉虹, 何志義, 白 璇, 呂 姣

缺血性卒中(ischemic stroke,IS)是一個(gè)復(fù)雜性狀的多基因遺傳病,遺傳變異決定了其遺傳易感性,但其確切的遺傳機(jī)制尚不清楚[1]。近年來(lái),IS在遺傳學(xué)方面新的進(jìn)展主要集中在一個(gè)新的炎癥相關(guān)基因,即 ALOX5AP基因的研究上。ALOX5AP基因作為 IS的易感基因,是 2004年 DECODE[2]研究組首先鑒定出來(lái)的,此后 ALOX5AP基因遺傳變異與 IS的相關(guān)性在其他人群中也得到了驗(yàn)證。國(guó)內(nèi)的研究顯示 ALOX5AP基因第 1內(nèi)含子 SG13S114A/T多態(tài)性位點(diǎn),A等位基因與 IS的發(fā)病密切相關(guān),但另一部分研究結(jié)果則不一致。推測(cè)這種相互矛盾結(jié)果的得出與樣本含量不足有關(guān)。為此,本研究采用 Meta分析的方法,對(duì)目前國(guó)內(nèi)進(jìn)行的相關(guān)研究進(jìn)匯總分析,以期更有效地評(píng)價(jià)兩者之間的關(guān)系,提供循證醫(yī)學(xué)資料。

1 資料與方法

1.1 研究對(duì)象 通過(guò)收集公開(kāi)發(fā)表的有關(guān)中國(guó)人群 ALOX5AP基因 SG13S114A/T多態(tài)性與 IS相關(guān)性的病例-對(duì)照研究資料,以其中的 IS患者為研究病例,以健康人群為對(duì)照組,觀察指標(biāo)為兩組基因型頻數(shù)為優(yōu)勢(shì)比值(OR),及等位基因頻數(shù)的 OR值。

1.2 資料來(lái)源 以缺血性卒中或腦梗死、5脂氧合酶激活蛋白基因或 ALOX5AP基因?yàn)橹黝}詞,聯(lián)合檢索 2004年至今的中國(guó)生物醫(yī)學(xué)文獻(xiàn)數(shù)據(jù)庫(kù),中國(guó)學(xué)術(shù)期刊全文數(shù)據(jù)庫(kù),萬(wàn)方數(shù)據(jù)庫(kù),及谷歌學(xué)術(shù)搜索引擎獲得全部中文文獻(xiàn)。并檢索同期 Medline在外文期刊發(fā)表的中國(guó)人資料,同時(shí)對(duì)萬(wàn)方資料庫(kù)博士和碩士學(xué)位論文進(jìn)行檢索,將未公開(kāi)發(fā)表的學(xué)位論資料納入本研究范圍。為避免上述電子數(shù)據(jù)庫(kù)未收錄和未發(fā)表的文獻(xiàn)資料,還對(duì)國(guó)內(nèi)一些大型學(xué)術(shù)會(huì)議的論文匯編進(jìn)行檢索。通過(guò)以上方法獲得有關(guān)中國(guó)人群 ALOX5AP基因 SG13S114A/T多態(tài)性與IS相關(guān)性的病例-對(duì)照研究文獻(xiàn)。

1.3 納入標(biāo)準(zhǔn)及研究質(zhì)量控制 由神經(jīng)內(nèi)科醫(yī)師和流行病學(xué)工作者采用盲法獨(dú)立進(jìn)行文獻(xiàn)提取工作,并共同對(duì)文獻(xiàn)進(jìn)行質(zhì)量評(píng)估取得一致性意見(jiàn):(1)研究對(duì)象均為漢族人或以漢族人為主體,不包括少數(shù)民族為主體的研究;(2)均為獨(dú)立的研究,同一組資料重復(fù)報(bào)道者,取其資料完整的文獻(xiàn);(3)各文獻(xiàn)研究方法、診斷標(biāo)準(zhǔn)及對(duì)照組人群的納入標(biāo)準(zhǔn)相同;(4)符合 Hardy-Weinberg(H-W)遺傳平衡定律。

1.4 數(shù)據(jù)處理與統(tǒng)計(jì) (1)對(duì)納入研究的基因型分布進(jìn)行 H-W平衡檢驗(yàn);(2)應(yīng)用漏斗圖初步分析納入的文獻(xiàn)是否存在明顯的發(fā)表偏倚,然后根據(jù)Egger法進(jìn)行偏倚檢驗(yàn)[3];(3)應(yīng)用 RVEMAN4.2(國(guó)際循證醫(yī)學(xué)協(xié)會(huì)組提供)統(tǒng)計(jì)軟件計(jì)算各個(gè)研究的SG13S114位點(diǎn)(AA+AT)/TT的 OR值,等位基因頻數(shù) A/T的 OR值,及其 95%CI;(4)以 Q檢驗(yàn)對(duì)各研究結(jié)果進(jìn)行同質(zhì)性檢驗(yàn),若無(wú)異質(zhì)性,則采用 Peto固定效應(yīng)模型進(jìn)行數(shù)據(jù)合并,計(jì)算總 OR值,若存在顯著異質(zhì)性,則采用 D-L隨機(jī)效應(yīng)模型法進(jìn)行數(shù)據(jù)合并,檢驗(yàn)水準(zhǔn)取 0.05。

2 結(jié) 果

2.1 文獻(xiàn)檢索結(jié)果 共檢出 10篇文獻(xiàn),有 4篇文獻(xiàn)是同組資料的重復(fù)性研究,取其中一篇資料完整的文獻(xiàn),其他 3篇被剔除;另有 2篇文獻(xiàn)報(bào)道的多態(tài)性位點(diǎn)為 SG13S100,也被剔除。最后納入本研究的為 5篇文獻(xiàn),包括 2篇學(xué)位論文,其余 3篇發(fā)表于中國(guó)科技核心期刊。5項(xiàng)研究中,IS組為 2620例,對(duì)照組為 2411例;所有入選病例均符合國(guó)際疾病分類(lèi)第 9版(ICD9)的診斷標(biāo)準(zhǔn),且都符合 H-W遺傳平衡定律;有 4個(gè)研究結(jié)果顯示 SG13S114A/T多態(tài)性與 IS具有相關(guān)性,1個(gè)研究結(jié)果顯示沒(méi)有相關(guān)性。入選文獻(xiàn)的資料及基因型分布情況見(jiàn)表1,等位基因頻數(shù)分布見(jiàn)表2。

表1 入選文獻(xiàn)一般情況及 SG13S114A/T基因型分布

2.2 入選文獻(xiàn)的發(fā)表偏倚 對(duì)漏斗圖目測(cè)結(jié)果顯示,未見(jiàn)明顯偏倚(見(jiàn)圖1,圖2)。按 Egger法,(AA+AT)/TT的 OR值除以其標(biāo)準(zhǔn)誤與標(biāo)準(zhǔn)誤倒數(shù)所建立的直線(xiàn)回歸方程的 Y軸截距別為 0.042,95%CI為-2.248～2.333,包含 0,P=0.957;A/T的OR值除以其標(biāo)準(zhǔn)誤與標(biāo)準(zhǔn)誤倒數(shù)所建立的直線(xiàn)回歸方程的 Y軸截距別為-0.273,95%CI為-3.225～2.697,包含 0,P=0.788,因此可以認(rèn)為沒(méi)有顯著性偏倚存在。

表2 入選文獻(xiàn)SG13S114A/T等位基因頻數(shù)分布

圖1 SG13S114A/T基因型(AA+AT)/TT漏斗圖,橫軸為 OR值,縱軸為 OR值自然對(duì)數(shù)的標(biāo)準(zhǔn)誤

圖2 SG13S114A/T等位基因 A/T漏斗圖,橫軸為 OR值,縱軸為 OR值自然對(duì)數(shù)的標(biāo)準(zhǔn)誤

2.3 異質(zhì)性評(píng)估 同質(zhì)性檢驗(yàn)顯示,SG13S114A/T位點(diǎn)基因型(AA+AT)/TT的 OR值統(tǒng)計(jì)量 χ2=1.00,P=0.91;等位基因頻數(shù) A/T的OR值統(tǒng)計(jì)量 χ2=8.18,P=0.09,按照 α=0.05水準(zhǔn)認(rèn)為各研究的基因型結(jié)果和等位基因頻數(shù)結(jié)果之間的一致性均好。

2.4 研究對(duì)象基因型分析 根據(jù)同質(zhì)性檢驗(yàn)結(jié)果,應(yīng)用 Peto固定模型模型法合并數(shù)據(jù),計(jì)算 OR值,并繪制森林圖。5個(gè)研究合并后,SG13S114A/T位點(diǎn)(AA+AT)/TT合并的 OR值為 1.22,95%CI為1.09～1.36,OR值顯著性檢驗(yàn) Z值為 3.5,P＜0.05(見(jiàn)圖3);等位基因頻數(shù) A/T的 OR值為 1.05,95%CI為 0.96～1.14,OR值顯著性檢驗(yàn) Z值為 1.09,P=0.28(見(jiàn)圖4)。

3 討 論

目前針對(duì) IS遺傳易感性獨(dú)立的關(guān)聯(lián)性研究均存在樣本量不足,統(tǒng)計(jì)校能低的缺點(diǎn)。通過(guò) Meta分析,可以利用現(xiàn)有研究的數(shù)據(jù)彌補(bǔ)樣本量的不足提高統(tǒng)計(jì)效力,鑒定 IS有效的易感基因,同時(shí)也可以對(duì) IS的遺傳危險(xiǎn)因素行定量地研究。Meta分析中存在著各種偏倚,但以發(fā)表偏倚最為常見(jiàn)。漏斗形圖是初步判斷發(fā)表偏倚簡(jiǎn)便易行的方法,但其缺點(diǎn)是只能對(duì)結(jié)果作定性的判斷。而 Egger法[3]可對(duì)發(fā)表偏倚作定量的評(píng)估,是漏斗形圖有效的補(bǔ)充。目前多將兩者聯(lián)合應(yīng)用檢驗(yàn)發(fā)表偏倚。本研究結(jié)果所示,根據(jù) Egger的線(xiàn)性回歸法,所建立的基因型和等位基因頻數(shù)統(tǒng)計(jì)效應(yīng)量的線(xiàn)性回歸方程 Y軸截距的95%CI均包含 0,且 P＞0.1,結(jié)合漏斗圖的目測(cè)結(jié)果,可以認(rèn)為沒(méi)有明顯的發(fā)表偏倚。

ALOX5AP基因編碼的產(chǎn)物是 5-脂氧合激活蛋白(5-Lipoxygenase activiting protein,FLAP),5-脂氧合酶(5-Lipoxygenase,5-LO)被 FLAP激活后可以將花生四烯酸轉(zhuǎn)化成白三烯(leukotriene,LT),LT被認(rèn)為是典型的炎癥細(xì)胞因子,因此 ALOX5AP基因也可以被為是炎癥因子的關(guān)鍵調(diào)控基因。2004年 DECODE研究組通過(guò)應(yīng)用全基因組掃描的方法發(fā)現(xiàn)ALOX5AP基因 HapA單倍型可以明顯增加北歐人群IS的遺傳易感性。這一結(jié)果引起了各國(guó)學(xué)者的極大關(guān)注,我國(guó)學(xué)者也相繼開(kāi)展了 ALOX5AP基因多態(tài)性與 IS相關(guān)性的研究,但研究結(jié)果并不一致。

本研究在入選的有關(guān)第 1內(nèi)含子 SG13S114A/T多態(tài)性與中國(guó)人 IS相關(guān)性的 5個(gè)病例對(duì)照組比較中,4項(xiàng)研究報(bào)告該位點(diǎn)的遺傳變異與 IS密切相關(guān),Meta分析數(shù)據(jù)合并后結(jié)果顯示 SG13S114A/T位點(diǎn)AA基因型可以顯著增加我國(guó)人群 IS的遺傳易感性。該位點(diǎn)多態(tài)性增加 IS易感性的機(jī)制目前尚未完全明確,但多數(shù)研究提示,ALOX5AP基因在轉(zhuǎn)錄過(guò)程中 SG13S114位點(diǎn)的變異可能具有某種基因調(diào)控作用,導(dǎo)致 FLAP產(chǎn)生增多從而上調(diào) LT活性,從而加重動(dòng)脈粥樣硬化的炎癥反應(yīng)導(dǎo)致 IS的形成[9]。針對(duì)A等位基因研究顯示,4個(gè)獨(dú)立的研究均發(fā)現(xiàn)A等位基因與 IS存在相關(guān)性,但數(shù)據(jù)合并后,A/T的OR值為 1.05,P=0.28,沒(méi)有發(fā)現(xiàn)相關(guān)性。究其原因可能是:(1)SG13S114A等位基因在人群中出現(xiàn)的頻率較低,而且納入本研究的文獻(xiàn)較少;(2)納入文獻(xiàn)的樣本量較小,導(dǎo)致 Meta分析統(tǒng)計(jì)校能較低;(3)在 ALOX5AP基因中存在其他的導(dǎo)致 IS易感性增加的遺傳變異。一項(xiàng)由 Zintzaras等[10]進(jìn)行的 Meta分析結(jié)果顯示,SG13S114A/T位點(diǎn)合并的 OR值為1.03,P=0.39與 IS無(wú)相關(guān)性,而 SG13S89G/A位點(diǎn)合并的 OR值為 1.25,P＜0.01與 IS相關(guān)。導(dǎo)致這種差異的形成,可能與人群遺傳背景的不同,研究方法,入選對(duì)象入選標(biāo)準(zhǔn),及性別構(gòu)成等有關(guān)。

圖3 SG13S114A/T基因型(AA+AT)/TT森林圖

圖4 SG13S114A/T等位基因頻數(shù)A/T森林圖

Meta分析作為一種現(xiàn)代科學(xué)研究中定量的文獻(xiàn)分析方法,在 IS遺傳學(xué)研究中的應(yīng)用越來(lái)越廣泛。相信隨著 ALOX5AP基因多態(tài)性和 IS相關(guān)性研究的深入和研究資料的積累,我們可以通過(guò) Meta分析進(jìn)一步明確 ALOX5AP基因的遺傳變異對(duì) IS發(fā)病機(jī)制的影響,及為 IS的早期診斷和綜合預(yù)防提供更多的依據(jù)。

[1] Dichgans M.Genetics of ischaemic stroke[J].Lancet Neurol,2007,6:149-161.

[2] Helgadottir A,Manolescu A,Thorleifsson G,et al.The gene encoding 5-lipoxygenase activating protein confers risk of myocardial infarction and stroke[J].Nat Genet,2004,36:233-238.

[3] Egger M,Davey SG,Schneider M,et al.Bias in Meta-analysis detected by a simple,graphical test[J].BMJ,1997,315:629-634.

[4] 張雙彥,鄔英全.腦梗死與 ALOX5AP的 SNPs及 HLA-DRB基因多態(tài)性的關(guān)聯(lián)研究[J].吉林大學(xué)學(xué)報(bào),2006,38:43.

[5] Zhang WL,Yang XM,Shi J,et al.Polymorphism of SG13S114T/A in the ALOX5APGeneand the Risk for Stroke in a Large Chinese Cohort[J].Acta Genetica Sinica,2006,33:678-684.

[6] 賀 穎,董子明.五個(gè)缺血性腦血管病候選基因的 SNP分析及遺傳易感性研究[J].鄭州大學(xué)學(xué)報(bào),2007,32:54.

[7] 高志嵩,杜秦川.缺血性腦卒中 ALOX5AP基因 SG13S114T/A和SG13S89G/A多態(tài)性研究[J].寧夏醫(yī)學(xué)雜志,2008,30:388-390.

[8] 徐佳亮,何志義,李 楠,等.5-脂氧合酶激活蛋白基因和白介素-1A基因多態(tài)與缺血性腦卒中的相關(guān)性[J].中華神經(jīng)科雜志,2009,42:803-807.

[9] Kaushal R,Pal P,Alwell K,et al.Association of ALOX5APwith ischemic stroke:a population-based case-control study[J].Hum Genet,2007,121:601-608.

[10] Zintzaras E,Rodopoulou P,Sakellaridis N.Variantsof the Arachidonate 5-Lipoxygenase-Activating Protein(ALOX5AP)Gene and Risk of Stroke:A HuGEGene-Disease Association Review and Meta-Analysis[J].Am J Epidemiol,2009,169:523-532.