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

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

1 材料和方法

1.1 文獻(xiàn)檢索策略

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

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

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

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

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

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

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

2 結(jié) 果

2.1 納入文獻(xiàn)的基本特點(diǎn)及質(zhì)量評價(jià)

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

圖1 文獻(xiàn)篩選流程圖

2.2 Meta分析結(jié)果

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

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

影響性分析結(jié)果顯示，在每一次刪掉1篇然后合并剩余納入文獻(xiàn)的情況下，ORs以及95%CI并沒有明顯改變(圖3)，表明本研究結(jié)果可靠。發(fā)表偏倚通過漏斗圖的形狀以及P值來評價(jià)。漏斗圖見圖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分析納入研究文獻(xiàn)的基本信息

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

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

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

3 討 論

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

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

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

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