中國人群PARK16位點rs947211基因多態(tài)性與帕金森病易感性的Meta分析

(青島大學(xué)附屬醫(yī)院，山東 青島266003 1 神經(jīng)外科； 2 神經(jīng)外科監(jiān)護室)

帕金森病(PD)是一種常見神經(jīng)退行性疾病，其特點是分布在黑質(zhì)致密部和路易小體的多巴胺能神經(jīng)元的損失[1]。PD的確切病因和發(fā)病機制尚未完全明確,遺傳和環(huán)境因素可能在PD的發(fā)生和進展過程中起著至關(guān)重要的作用[2]。PARK16位點最先是由SATAKE等[3]通過全基因組關(guān)聯(lián)研究(GWAS)在日本人群中得到確定，證明了PARK16位點和PD易感性有很強的關(guān)聯(lián)性。隨后，在歐洲和北美PD人群中也分別進行了GWAS,結(jié)果顯示PARK16位點單核苷酸多態(tài)性(SNP)與PD的發(fā)病風(fēng)險無關(guān)[4-5]。目前已有相關(guān)的病例-對照研究文獻報道了PARK16位點rs947211基因的多態(tài)性與中國人群PD發(fā)病的關(guān)系[6-13]。本研究擬對中國人群中PARK16位點rs947211基因多態(tài)性和PD易感性的關(guān)系進行Meta分析，從而為臨床研究提供理論依據(jù)。

1 材料和方法

1.1 文獻檢索策略

由兩位研究人員系統(tǒng)地對PubMed、Google Scholar、EMBASE, Web of Science、中國知網(wǎng)、維普網(wǎng)以及萬方數(shù)據(jù)知識服務(wù)平臺上有關(guān)中國人群PARK16位點rs947211基因多態(tài)性和PD易感性關(guān)系研究的文獻進行檢索,檢索相關(guān)關(guān)鍵詞如下：“PARK16 or rs947211”、“Parkinson’s disease or PD or Parkinsonism”、“PARK16 或者rs947211”與“帕金森病”。同時兩位研究人員采用手工搜索及文獻追溯的方法以確定其他潛在可利用的研究文獻。

1.2 文獻納入及排除標(biāo)準(zhǔn)

納入標(biāo)準(zhǔn):(1)評價rs947211基因多態(tài)性與PD易感性關(guān)系的病例-對照或隊列研究；(2)文獻原始資料完整，均可獲得各基因型頻數(shù)；(3)研究樣本重疊時,優(yōu)先納入較大樣本或最近的研究；(4)對照組的基因型分布必須符合Hardy-Weinberg(H-W)遺傳平衡定律。排除標(biāo)準(zhǔn):動物實驗研究文獻、非病例對照或隊列研究文獻、基因型數(shù)據(jù)不完整的研究文獻，以及病例報告、綜述、摘要或者Meta分析研究文獻。

1.3 數(shù)據(jù)提取與質(zhì)量評價

兩位研究者分別閱讀納入的文獻，然后進行資料提取。提取內(nèi)容如下:第一作者、出版時間、國家以及種族、基因分型方法、PD病例及對照組的基因型分布。意見不一致處與第3位研究者進行討論解決。采用Newcastle-Ottawa Scale(NOS)量表對納入的所有文獻進行質(zhì)量評價?？偡?分，評分達6分以上為高質(zhì)量文獻。

1.4 統(tǒng)計學(xué)方法

根據(jù)WOOLF等[14]的方法分別計算5種模型(等位基因模型: A vs.G，顯性模型：AA+AG vs.GG，純合子模型：AA vs.GG，雜合子模型：AG vs.GG，隱性模型：AA vs.AG+GG)的ORs和相應(yīng)的95%置信區(qū)間(95%CI)評價rs947211基因多態(tài)性與PD易感性的關(guān)系。同樣，根據(jù)不同民族本研究進行了亞組分析。應(yīng)用I2對各研究間異質(zhì)性進行定量分析。如各研究間具有同質(zhì)性(I2<50%)，則采用固定效應(yīng)模型對統(tǒng)計量進行合并；如各研究間具有異質(zhì)性(I2≥50%)，則采用隨機效應(yīng)模型進行分析[15]。同時采用Begg’s 檢驗來評估發(fā)表偏倚，并繪制漏斗圖；P<0.05和(或)不對稱的漏斗圖表明可能存在潛在的發(fā)表偏倚[16]。為了評估Meta分析結(jié)果穩(wěn)定性，應(yīng)用依次從納入的所有文獻中去除1篇文獻然后合并剩余的文獻的方法對文獻進行影響性分析。所有的統(tǒng)計學(xué)檢驗均采用Stata12.0軟件(Stata Corp, College Station, Texas 77845, USA)進行分析。

2 結(jié) 果

2.1 納入文獻的基本特點及質(zhì)量評價

圖1為文獻篩選過程的流程圖。初檢文獻為302篇。按照上述文獻納入標(biāo)準(zhǔn)和排除標(biāo)準(zhǔn)對文獻進行篩選，最終納入本次Meta分析的文獻總共8篇(5篇英文和3篇中文文獻)。各對照組基因型均符合H-W遺傳平衡定律。納入文獻NOS量表評分6～9分，文獻質(zhì)量總體較高。文獻的基本特征見表1。

圖1 文獻篩選流程圖

2.2 Meta分析結(jié)果

通過5種模型對納入的文獻進行異質(zhì)性檢驗，結(jié)果顯示本次納入的文獻之間異質(zhì)性均小于50%(見表2)，采用固定模型進行分析。研究結(jié)果顯示，PARK16位點rs947211基因多態(tài)性與中國人群PD的發(fā)病風(fēng)險呈負(fù)相關(guān)的關(guān)系。亞組分析結(jié)果顯示，rs947211基因多態(tài)性不僅與降低漢族人群PD發(fā)病風(fēng)險有關(guān)，也與降低非漢族人群PD發(fā)病風(fēng)險有關(guān)(見表2、圖2)。

2.3 影響性分析和發(fā)表偏倚評價

影響性分析結(jié)果顯示，在每一次刪掉1篇然后合并剩余納入文獻的情況下，ORs以及95%CI并沒有明顯改變(圖3)，表明本研究結(jié)果可靠。發(fā)表偏倚通過漏斗圖的形狀以及P值來評價。漏斗圖見圖4，可見基本對稱；此外，各模型P值大于0.05(A vs. G:P=0.062；AA+AG vs.GG:P=0.213；AA vs.GG:P=0.119；AG vs.GG:P=0.213；AA vs.AG+GG:P=0.276)。以上表明本研究不存在明顯的發(fā)表偏倚。

表1 Meta分析納入研究文獻的基本信息

表2 PARK16位點rs947211基因多態(tài)性與PD易感性Meta分析結(jié)果

圖2 PARK16位點rs94721基因多態(tài)性與PD相關(guān)性森林圖(A vs.G)

圖3 PARK16位點rs947211基因多態(tài)性敏感性分析(A vs.G)

3 討 論

PD是一種呈進行性發(fā)展的神經(jīng)退行性疾病，全球85歲以上人群中約有2%以上的人群受到了該病的影響[17]。為了解其發(fā)病相關(guān)因素，本研究通過Meta分析評估了PARK16位點rs947211基因的多態(tài)性和PD發(fā)病風(fēng)險之間的關(guān)系，結(jié)果顯示PARK16位點rs947211基因的多態(tài)性與中國人群PD的發(fā)病風(fēng)險呈負(fù)相關(guān)關(guān)系。PARK16位點位于染色體1q32上，它包括5個PD易感候選基因：NUCKSl、SLC45A3、RAB7Ll、PM20Dl和 SLC4lAl。相對于家族性PD來說，PARK16對于散發(fā)性老年P(guān)D病人的保護性更為顯著[18]。最初在對日本人群GWAS中被發(fā)現(xiàn)其與PD相關(guān)聯(lián)，雖然隨后對亞洲其他地區(qū)以及高加索地區(qū)人群中進行了重復(fù)研究[19]，但其發(fā)病機制還不清楚。rs947211基因多態(tài)性位于SLC41A1和RAB7L1基因之間，而且其與NUCKS1基因的轉(zhuǎn)錄水平密切相關(guān)[20]。周期蛋白依賴性激酶由NUCKS1基因進行編碼，rs947211多態(tài)性在細(xì)胞周期過程中也扮演著重要的角色[21]?？紤]到細(xì)胞周期蛋白在多巴胺能神經(jīng)元凋亡過程中起到重要作用[22-24]，NUCKS1基因轉(zhuǎn)錄及翻譯可能與PD的發(fā)病機制有關(guān)?？傊?該研究結(jié)果顯示，PARK16位點rs947211基因的多態(tài)性與中國人群PD的發(fā)病風(fēng)險呈負(fù)相關(guān)關(guān)系。亞組分析結(jié)果顯示，rs947211基因多態(tài)性不僅與降低漢族人群PD發(fā)病風(fēng)險有關(guān)，而且也與降低非漢族人群PD發(fā)病風(fēng)險有關(guān)。該項研究的優(yōu)勢在于針對中國地區(qū)人群進行，且對漢族人群與非漢族人群進行了亞組分析，受地區(qū)及基因差異影響相對較?。淮送?，雖然在所有5個模型中沒有發(fā)現(xiàn)顯著的異質(zhì)性,但該項研究中存在的部分異質(zhì)性來源可能與所納入的樣本量仍不充足、環(huán)境因素和基因因素干擾、種族及地理位置差異等有關(guān)。

圖4PARK16位點rs947211基因多態(tài)性Begg’s漏斗圖發(fā)表偏倚評價(Avs.G)

本研究中存在如下一些局限性：首先,該Meta分析中僅包含了已經(jīng)發(fā)表的研究。因此,該Meta分析可能存在發(fā)表偏倚。其次，由于缺乏晚發(fā)型PD(LOPD)和早發(fā)型PD(EOPD)基因型分布的數(shù)據(jù),我們只對PARK16位點rs947211基因多態(tài)性與PD的發(fā)病風(fēng)險之間的關(guān)系進行了Meta分析，而沒有進一步對LOPD和EOPD進行Meta分析。第三,PD的發(fā)病原因比較復(fù)雜，除研究遺傳方面的因素外，還應(yīng)該進一步分析基因與環(huán)境及基因與基因的交互作用。

綜上所述，該Meta分析證明了PARK16位點rs947211基因多態(tài)性可以降低中國人群PD的發(fā)病風(fēng)險。此外,考慮到上述提到的局限性,仍需更大樣本量的病例-對照研究來驗證PARK16位點rs947211基因多態(tài)性和PD發(fā)病風(fēng)險之間的關(guān)系。