HLA-DRA rs3129882多態性與帕金森病關系的meta分析*

曹美群，馮毅菲，蔣 屏，范大華

(1.廣東省深圳市第二人民醫院老年醫學研究所 518035；2.廣西中醫藥大學研究生學院，南寧 530200；3.湖南中醫藥大學第一附屬醫院兒科，長沙 410007；4.深圳金普邁生物科技有限公司，廣東 518101)

帕金森病(PD)是僅次于阿爾茨海默病的第二大神經退行性疾病，全世界有600多萬人受到影響[1]。目前，其病因仍不十分清楚，但越來越多的證據表明，多個基因與PD發病有關[1-2]。最近一項研究表明，高加索人PD與Ⅱ類主要組織相容性復合體DRα(HLA-DRA)內含子1上的一個非編碼區單核苷酸多態性位點 rs3129882關聯性很高[3]。然而，后續系列研究試圖驗證PD與HLA-DRA rs3129882的遺傳關聯卻得到了不一致的結果[4-6]。為探究以上研究結論矛盾的原因，本研究通過對15篇文獻包括28項獨立研究進行meta分析，評估了HLA-DRA rs3129882與PD的相關性。

1 資料與方法

1.1 研究對象

在Medline數據庫和Google Scholar檢索，將研究HLA-DRA rs3129882多態性與PD關系的英文文獻作為研究對象。納入標準：(1)評估rs3129882多態性與PD的關聯性；(2)代表普通人群的隊列研究或診斷為PD的病例對照研究；(3)提供優勢比(OR)及其95%可信區間(95%CI)或提供足夠的統計數據。排除標準：(1)非英文形式發表的文獻；(2)綜述、評論、社論、通訊、meta分析和摘要；(3)重復發表的文獻(當兩個或更多的研究使用了相同的樣本或大部分重疊的樣本以報告關聯結果時只選擇樣本量最大的研究)。

1.2 方法

1.2.1文獻檢索

在Medline數據庫和Google Scholar檢索HLA-DRA rs3129882多態性與PD關系的相關研究文獻。關鍵詞為帕金森病、人類白細胞抗原或HLA。檢索時間為建庫至2020年12月，同時追溯檢索原始文章和綜述的參考文獻使文獻檢索盡可能詳盡。

1.2.2資料提取

從符合條件的研究中提取以下信息：(1)第一作者姓名；(2)發表年份；(3)人群和種族；(4)PD患者組和對照組rs3129882基因多態性的次要/主要等位基因；(5)PD病例和對照者數量；(6)基因分型平臺；(7)OR和95%CI，如報道了潛在混雜因素的調整，提取調整后的OR和95%CI，在未報道調整OR或95%CI的情況下用報道的基因型/等位基因數據計算原始OR和95%CI。

1.2.3文獻質量評價標準

采用畢卡斯爾-渥太華量表(NOS量表)評價文獻質量,病例及對照的選擇計4分(包括病例的定義是否充分、病例的代表性、對照的選擇和對照的定義)、基于設計或分析所得的病例與對照的可比性計2分、暴露計3分(包括暴露的確定、病例及對照的暴露是否采用了相同的確定方法和無應答率)共3個大項9個小項，滿分為9分，≥5分為高質量文獻[7]。

1.3 統計學處理

采用Review Manager(Revman)5.3、Stata11.0軟件進行數據分析，以P<0.05為差異有統計學意義。各納入研究異質性檢驗采用CochraneQ檢驗，計算I2值，I2為0～<25%、25%～<50%、50%～<75%、75%～100%分別表示低、中、高、極不均質性。若研究間無統計學異質性(I2≤50%，P≥0.10)，選擇固定效應模型進行meta分析[8-9]；若存在統計學異質性(I2>50%，P<0.10)則采用隨機效應模型進行meta分析。使用逆方差法計算合并OR。敏感性分析通過一一剔除文獻重新進行meta 分析，用漏斗圖、Begg’s和Egger’s檢驗評估潛在的發表偏倚[10]。

2 結 果

2.1 文獻篩選流程

在Medline上搜索得到171篇文獻，其中12篇符合納入標準；手動檢索已納入文獻所引用的參考文獻，補充納入2篇文獻[6，11]；檢索Google Scholar，補充納入1篇文章[12]，共納入文獻15篇[5,11-24]，包括28項獨立研究。文獻篩選流程見圖1。納入文獻基本特征見表1。

圖1 文獻篩選流程

表1 納入文獻基本特征

2.2 meta分析結果

15篇文獻包括28項獨立研究(21項高加索人群研究[5,11-17,23]和7項亞洲人群研究[18-24])納入meta分析，共74 715名受試者(24 269名PD患者和50 446名健康對照者)，其中63 778名受試者(19 509名PD患者和44 269名健康對照者)為高加索人，10 937名受試者(4 760名PD患者和6 177名健康對照者)為亞洲人。15篇文獻[5,11-24]為混合人群的研究，各研究間存在高度統計學異質性(I2=56%，P=0.000 2),采用隨機效應模型進行meta分析。混合人群rs3129882變異與PD無關(OR=1.04,95%CI：0.99～1.09，P=0.120 0)。7篇文獻[18-24]為亞洲人群的研究，各研究間存在高度統計學異質性(I2=68%，P=0.004 0)，采用隨機效應模型進行meta分析。亞洲人群rs3129882變異與PD無關(OR=0.99，95%CI：0.87～1.12，P=0.820 0)。9篇文獻[5,11-17,23]為高加索人群的研究，各研究間無統計學異質性(I2=46%，P=0.010 0)，采用固定效應模型進行meta分析。高加索人群rs3129882變異與PD明顯相關(OR=1.04,95%CI：1.01～1.08，P=0.005 0)。見表2。

表2 rs3129882多態性與PD關系的亞組分析

2.3 敏感度分析與發表偏倚

通過依次剔除高加索人群或亞洲人群中的研究進行敏感性分析發現結果和結論一致。使用Revman創建漏斗圖評估兩個亞組研究的發表偏倚風險見圖2。漏斗圖沒有明顯的不對稱性，高加索人群的Begg’s檢驗P=0.291，Egger’s檢驗P=0.186，亞洲人群的Begg’s檢驗P=0.230，Egger’s檢驗P=0.065。

圖2 評價納入研究發表偏倚的漏斗圖

3 討 論

研究表明，很多基因與PD的發病有關，目前，已經在家族性PD中發現了18個相關位點，在散發性PD中也發現了多個相關位點[1-2]。最近一項對美國歐洲高加索人的研究發現，HLA 6號染色體短臂2區1帶3亞帶1次亞帶(6p21.3)區域與PD相關，且與HLA-DRA rs3129882關聯性最高[3]。然而，后續系列PD與HLA-DRA rs3129882的遺傳關聯研究卻得到了不一致的結果[4-6]。DO等[18]在高加索人群中進行的一項大規模全基因組關聯，研究結果顯示，rs3129882與PD存在關聯，而其他樣本量較小的高加索人研究卻沒有發現這種關聯[19-21]。PUSCHMANN等[22]在愛爾蘭和波蘭混合人群中研究，結果甚至顯示rs3129882的“G”等位基因與PD呈負相關。同時亞洲人群的幾項rs3129882與PD關聯的病例對照研究結果也不一致。中國大陸人群的兩項研究顯示，rs3129882是PD的易感基因[23-24]，而中國臺灣人群中的兩項研究卻未顯示出這樣的關聯[4，8]，甚至還有研究發現，rs3129882多態性與新加坡和馬來西亞華裔人群患PD的風險降低有關[5]。最近一項包括中國、新加坡、韓國和日本等地區的亞洲混合人群的大規模病例對照研究也未能顯示出這種關聯[6]。

本研究對15篇文獻包括28項獨立研究進行了meta分析，以評估HLA-DRA rs3129882在PD發病中的作用，結果顯示，在混合人群或亞洲人群樣本中未發現這種關聯，但在高加索人中rs3129882與PD明顯相關。這種在高加索人群和亞洲人群中rs3129882與PD的關聯不一致很可能與不同種族和地理群體之間基因多態性的遺傳變異(也稱為群體亞結構)有關，在同一種族但不同地域人群中甚至存在重要的亞結構[3]。HAMZA等[3]在高加索人群中觀察到，rs3129882風險等位基因的頻率梯度范圍在北歐受試者的0.35至南歐受試者的0.46之間，同樣，在東亞華人人群中也有顯著變動，其水平從中國大陸的0.28升高至中國臺灣的0.37[4-5,8,23-24]。此外，值得注意的是，rs3129882的次要等位基因從高加索人的鳥嘌呤“G”轉換為亞洲人的腺嘌呤“A”[3，18，23-24]。因此，如不考慮混合地域來源樣本的種群亞結構，rs3129882與PD無關，甚至是反向關聯[5]。未來對HLA區域的遺傳變異和PD易感性研究應該按種族和地理來源進行仔細分層。

本研究雖已證實rs3129882與高加索人PD的風險有關聯，但其潛在的分子機制仍有待于闡明。有研究表明，rs3129882位點突變會影響HLA-DR和HLA-DQ的基因轉錄和選擇性剪接，最終影響HLA-DRA和HLA-DRB抗原表達水平[10]。另有研究表明，自發性PD或神經毒素1-甲基-4-苯基-1,2,3,6-四氫吡啶(MPTP)誘導的PD大腦HLA-DR抗原表達均顯著上調[25]，長期使用非甾體抗炎藥可顯著降低PD的風險[26-27]。表明神經炎癥在PD發病機制中具有重要作用[28]。因此，本研究發現并確認HLA-DRA rs3129882與高加索人PD發病相關，進一步證實了免疫系統在PD發病中的作用，為高加索人PD的治療提供了新靶點。

但本研究仍存在局限性：(1)僅使用全基因組關聯研究和meta分析所采用的可加性遺傳模型[29]，沒有使用顯性模型或隱性模型進行分析，無法獲得所有研究中3種基因型的確切數量。今后應更多地關注顯性和隱性模型下不同種族和地理來源人群中HLA-DRA rs3129882與PD的相關性。(2)僅關注了HLA-DRA rs3129882一個基因多態性，越來越多的證據表明，在HLA區域內存在多個與PD相關的變異[30]。今后應對HLA區域內受人群亞結構調整的PD相關變異進行全面系統分析,將有助于進一步解釋關聯研究結果的差異。