葛雅婷,蘇美玲,宋佳治，熊祖應(yīng)△,侯 霜▲

(1.安徽醫(yī)科大學(xué)北京大學(xué)深圳醫(yī)院臨床學(xué)院，廣東深圳 518000; 2.北京大學(xué)深圳醫(yī)院腎內(nèi)科，廣東深圳 518000)

Megsin基因3個多態(tài)性位點與亞洲人群IgA腎病易感性關(guān)聯(lián)的Meta分析*

葛雅婷1，2,蘇美玲2,宋佳治1，2，熊祖應(yīng)1，2△,侯 霜2▲

(1.安徽醫(yī)科大學(xué)北京大學(xué)深圳醫(yī)院臨床學(xué)院，廣東深圳 518000; 2.北京大學(xué)深圳醫(yī)院腎內(nèi)科，廣東深圳 518000)

目的 通過Meta分析評估亞洲人群Megsin基因3個多態(tài)性位點(rs1055901、rs1055902和rs2689399)與IgA腎病易感性的關(guān)聯(lián)。方法 通過中國知網(wǎng)、維普、萬方數(shù)據(jù)庫、中國生物醫(yī)學(xué)文獻(xiàn)數(shù)據(jù)庫、PubMed、Web of Science和Google學(xué)術(shù)數(shù)據(jù)庫數(shù)據(jù)庫，電子檢索關(guān)于亞洲人群Megsin基因rs1055901、rs1055902和rs2689399多態(tài)性位點與IgA腎病相關(guān)性的研究，檢索年限為1960年1月1日開始至2016年5月2日。采用Stata12.0軟件計算合并比值比(OR)和95%置信區(qū)間(CI)，并進(jìn)行敏感性及發(fā)表偏倚分析。結(jié)果 最終納入6篇文獻(xiàn)(含9項研究)，包含2 179例IgA腎病患者和1 769例健康對照。Meta分析結(jié)果顯示，Megsin基因rs1055901、rs1055902多態(tài)性位點與亞洲人群IgA腎病易感性無明顯關(guān)聯(lián)。然而，Megsin基因rs2689399位點與亞洲人群IgA腎病易感性呈顯著負(fù)相關(guān)(G和C：OR=0.754，95%CI0.592～0.961，P=0.022；GG和CC：OR=0.506，95%CI0.287～0.892，P=0.019；GG和GC+CC：OR=0.551，95%CI0.316～0.961，P=0.036)。結(jié)論 Megsin基因的rs2689399位點的G等位基因、GG基因型可能是亞洲人群IgA腎病的保護(hù)因子。

腎小球腎炎,IGA；疾病易感性；連鎖不平衡；Megsin；多態(tài)性；Meta分析

IgA腎病是全世界范圍內(nèi)常見的原發(fā)性腎小球疾病之一，也是導(dǎo)致終末期腎病的主要病因[1]。成年人中每年IgA發(fā)病率約2.5/10萬[2]。在初次腎活檢后20年內(nèi)，約30%～50%患者會進(jìn)展為終末期腎病[3]。腎活檢是診斷IgA腎病的金標(biāo)準(zhǔn)，其病理學(xué)特點主要是IgA為主沉積在系膜區(qū)伴系膜細(xì)胞和系膜基質(zhì)增生[4]。近年來，大量研究關(guān)注遺傳易感基因與IgA腎病的關(guān)聯(lián)，提示遺傳因素與環(huán)境的交互作用可能作用于該病的發(fā)生。

Megsin基因位于染色體18q21.3，主要是絲氨酸蛋白酶(Serpin)超家族成員，在諸多伴有系膜損傷的腎小球疾病的人類及動物模型中，證實Megsin在腎小球系膜細(xì)胞中表達(dá)上調(diào)，介導(dǎo)了系膜基質(zhì)增多和系膜細(xì)胞數(shù)量增加，這很可能是IgA腎病的發(fā)病機制之一[5-6]。一些研究已經(jīng)報道了Megsin基因3′段非編碼區(qū)rs1055901和rs1055902多態(tài)性位點與IgA腎病相關(guān)性[6-9]；也有研究報道C25663G多態(tài)性位點與漢族人群IgA腎病易感性的關(guān)系[10-11]，然而這些研究主要是小樣本資料的研究，且結(jié)果、結(jié)論尚未統(tǒng)一。

本研究利用系統(tǒng)性綜述的方法，通過增大樣本量，提高檢驗效能，進(jìn)一步揭示Megsin基因rs1055901、rs1055902及rs2689399 3個多態(tài)性位點與亞洲人群IgA腎病的關(guān)聯(lián)，為IgA腎病遺傳學(xué)研究提供依據(jù)。

1 資料與方法

1.1 檢索策略 用計算機檢索中國知網(wǎng)、維普、萬方數(shù)據(jù)庫、中國生物醫(yī)學(xué)文獻(xiàn)數(shù)據(jù)庫、PubMed、Web of Science和Google學(xué)術(shù)數(shù)據(jù)庫，中文以IgA腎病、Megsin、多態(tài)性為檢索詞，英文以“IgA nephropathy OR Immunoglobulin A nephropathy” AND “megsin” AND “variant OR mutation OR polymorphism OR SNP”為檢索詞。手工檢索國內(nèi)外有關(guān)Megsin基因多態(tài)性位點與亞洲人群IgA腎病發(fā)病風(fēng)險的文獻(xiàn)。所有文獻(xiàn)語種只限中、英文，檢索時間1960年1月1日至2016年5月2號。

1.2 納入和排除標(biāo)準(zhǔn) 納入標(biāo)準(zhǔn)：(1)納入研究均為病例-對照研究，研究人群為亞洲人；(2)符合IgA腎病診斷標(biāo)準(zhǔn)[12]；(3)關(guān)于Megsin基因rs1055901(核苷酸位點C2093T)、rs1055902(核苷酸位點C2180T)及rs2689399(核苷酸位點為C25663G)3個多態(tài)性位點與亞洲人群IgA腎病相關(guān)性研究；(4)排除基礎(chǔ)數(shù)據(jù)不完整及數(shù)據(jù)錯誤的研究，排除非病例對照研究、多態(tài)性位點非rs1055901、rs1055902或rs2689399的文獻(xiàn)；(5)對于重復(fù)發(fā)表或資料雷同的研究只保留其中質(zhì)量最好者。

1.3 文獻(xiàn)質(zhì)量評價及資料提取 對每個納入的研究參考紐卡斯?fàn)?渥太華量表(Newcastle-Ottawa scale,NOS)[13]標(biāo)準(zhǔn)，從病例-對照對象選擇、可比性、結(jié)局和暴露3個方面對文獻(xiàn)質(zhì)量評分(其中5～9分為相對高質(zhì)量的研究)。本分析提取的數(shù)據(jù)包括納入研究的作者、發(fā)表年限、種族、基因分型方法、對照組來源、Megsin基因病例和對照的例數(shù)及其3個位點基因頻率。由兩位評價者獨立提取數(shù)據(jù)后交叉核對及文獻(xiàn)質(zhì)量評價，不同意見與第三者討論解決。

1.4 統(tǒng)計學(xué)處理 采用STATA12.0軟件進(jìn)行統(tǒng)計學(xué)分析。利用5種遺傳學(xué)模型檢測Megsin基因多態(tài)性位點和IgA腎病的相關(guān)性，包括等位基因模型B和A，顯性模型BA+BB和AA，隱性模型BB和BA+A，純合子模型BB和AA，雜合子模型BA和AA(其中A為野生型位點，B為突變型位點)。使用Q或I2檢驗分析各研究間的異質(zhì)性，若P>0.1和I2<50%，則采用固定效應(yīng)模型；否則采用隨機效應(yīng)模型。采用比值比(OR)值及95%置信區(qū)間(CI)作為效應(yīng)量。各研究的對照組均檢測哈迪溫伯格平衡(Hardy Weinberg Equilibrium，HWE)。使用敏感性分析評估結(jié)果的穩(wěn)定性及判斷每個獨立研究對總體OR值有無影響，通過Begg′s及Egger′s檢驗評估發(fā)表偏倚。通過Hapmap數(shù)據(jù)庫，利用Haploview4.0軟件，分析Megsin基因rs1055901、rs1055902和rs2689399多態(tài)性位點在亞洲人群包括北京漢族人(CHB)及東京日本人(JPT)、非洲人群(伊巴丹約魯巴人，YRI)及歐洲人群(祖先來自西北歐的猶他州白人，CEU)中的連鎖不平衡關(guān)系。

2 結(jié) 果

2.1 文獻(xiàn)檢索結(jié)果及納入研究的基本特征 文獻(xiàn)篩選流程及結(jié)果見圖1。最終納入6篇文獻(xiàn)(含9個研究)[6-11]，包含2 179例IgA腎病患者(病例組)和1 769例健康對照(對照組)。研究對象均屬于亞洲人群，且對照組與病例組均具有可比性。對于rs1055901多態(tài)性位點，共納入4項研究，包含病例組856例，對照組786例。對于rs1055902多態(tài)性位點，共納入3項研究，包含病例組808例，對照組686例。對于rs2689399多態(tài)性位點，共納入2項研究，包含病例組515例，對照組297例。基因型檢測方法包括polymerase chain reaction-restriction fragment length polymorphism、sequencing、TaqMan、MassARRAY。納入研究的基本特征見表1。對納入研究對照組的基因型頻率進(jìn)行HWE檢測，發(fā)現(xiàn)有1項研究[11]不符合HWE檢驗(P<0.05)。

圖1 納入文獻(xiàn)流程及結(jié)果

SNP作者發(fā)表年限基因分型方法TT(CC)基因型例數(shù)(病例組/對照組)TC(CG)基因型例數(shù)(病例組/對照組)CC(GG)基因型例數(shù)(病例組/對照組)HWENOS評分(分)rs1055901Maixnerova等[6]2008sequencing56/1695/3546/100.2159Lim等[7]2008TaqMan33/35106/143121/1370.8008李露等[8]2013PCR-RFLP5/1017/5426/360.1138Wei等[9]2016MassARRAY27/19138/114186/1770.9107rs1055902Maixnerova等[6]2008sequencing53/1192/3352/170.4709Lim等[7]2008TaqMan123/135107/14730/330.4498Wei等[9]2016MassARRAY201/187131/11319/100.1507rs2689399夏運風(fēng)等[10]2006sequencing277/108169/7225/170.3217李露等[11]2011PCR-RFLP28/4612/364/180.0298

PCR-RFLP：聚合酶鏈反應(yīng)-限制性片段長度多態(tài)性；HWE：哈迪-溫伯格平衡；NOS：紐卡斯?fàn)?渥太華量表。

2.2 文章質(zhì)量評價 所納入的研究均為病例-對照研究，病例均來自醫(yī)院，納入方法合適，具有統(tǒng)一的診斷標(biāo)準(zhǔn)；對照組的選取標(biāo)準(zhǔn)合適，為體檢健康人群?；蛐偷臏y量方法準(zhǔn)確可靠，研究結(jié)果合理可信，且經(jīng)NOS評分后所納入的研究也波動在7～9分。因此，所有納入研究的方法學(xué)質(zhì)量均較高。

2.3 Meta分析結(jié)果

2.3.1 rs1055901多態(tài)性位點與IgA腎病易感性關(guān)系 對于rs1055901多態(tài)性位點，共納入4項研究。Meta分析結(jié)果顯示rs1055901多態(tài)性位點與IgA腎病易感性之間無明顯關(guān)聯(lián)(C和T：OR=1.004，95%CI：0.862～1.170，P=0.955；CC和TT：OR=0.953，95%CI：0.672～1.352，P=0.788；CT和TT：OR=0.788，95%CI0.562～1.105，P=0.168；CT+CC和TT：OR=0.853，95%CI0.620～1.175，P=0.331；CC和CT+TT：OR=1.187，95%CI0.836～1.686，P=0.338)，見表2。

2.3.2 rs1055902多態(tài)性位點與IgA腎病易感性關(guān)系 對于rs1055902多態(tài)性位點，共納入3項研究。結(jié)果顯示，rs1055902多態(tài)性位點與IgA腎病易感性無明顯關(guān)聯(lián)(C和T：OR=0.993，95%CI0.843～1.169，P=0.933；CC和TT：OR=1.049，95%CI0.708～1.555，P=0.810；CT和TT：OR=0.898，95%CI0.718～1.125，P=0.350；CT+CC和TT：OR=0.937，95%CI0.756～1.161，P=0.552；CC和CT+TT：OR=1.159，95%CI0.809～1.661，P=0.420)，見表2。

2.3.3 rs2689399多態(tài)性位點與IgA腎病易感性關(guān)系 對于rs2689399多態(tài)性位點，共納入2個研究。結(jié)果顯示，rs2689399等位基因模型、純合子模型及隱性模型與IgA腎病易感性呈顯著負(fù)相關(guān)(G和C：OR=0.754，95%CI0.592～0.961，P=0.022；GG和CC：OR=0.506，95%CI0.287～0.892，P=0.019；GG和GC+CC：OR=0.551，95%CI0.316～0.961，P=0.036)。而雜合子模型及顯性模型與IgA腎病易感性無明顯關(guān)聯(lián)(GC和CC：OR=0.839，95%CI0.608～1.159，P=0.287；GC+GG和CC：OR=0.769，95%CI0.568～1.042，P=0.090)。

表2 Megsin基因rs1055901、rs1055902和rs2689399多態(tài)性與IgA易感性Meta分析結(jié)果

OR：比值比；CI：置信區(qū)間；PH：異質(zhì)性P值；PZ：顯著性P值。

2.4 敏感性分析和發(fā)表偏倚 采用逐一排除的方法進(jìn)行敏感性分析，排除任一項研究后，Meta分析的總效應(yīng)量在剔除前后未發(fā)生明顯變化，提示結(jié)果較穩(wěn)定。此外，利用Begg′s和Egger′s檢驗評估發(fā)表偏倚，漏斗圖提示所有研究在圖中基本對稱，且Egger′s檢驗提示無明顯發(fā)表偏倚(rs1055902：P>|t|=0.424)。

圖2 通過Haploview軟件分析Megsin基因rs1055901、rs1055902和 rs2689399在CEU人群和YRI人群中分別構(gòu)成的單體型

2.5 Megsin rs1055901、rs1055902和rs2689399多態(tài)性位點的連鎖不平衡關(guān)系 通過Hapmap數(shù)據(jù)庫，利用Haploview4.0軟件分析，發(fā)現(xiàn)在CHB及JPT人群中未發(fā)現(xiàn)rs1055901、rs1055902和rs2689399存在連鎖關(guān)系；在CEU人群中發(fā)現(xiàn)Megsin rs1055901、rs1055902和rs2689399多態(tài)性位點之間標(biāo)準(zhǔn)化的連鎖不平衡系數(shù)D′=1，相關(guān)系數(shù)r2=1(圖2)，表明兩者呈完全連鎖；在YRI人群中，rs2689399和rs1055901兩位點之間呈完全連鎖(D′=1，r2=1)，而rs1055901和rs1055902兩個位點、rs2689399和rs1055901兩位點之間分別也存在一定程度連鎖，且連鎖程度一致(D′=1，r2=0.667，圖2)。

3 討 論

隨著分子生物學(xué)研究手段日益成熟與發(fā)展，國內(nèi)外學(xué)者發(fā)現(xiàn)IgA腎病的發(fā)生與發(fā)展與遺傳因素密切相關(guān)。近年諸多研究報道了一些候選基因如編碼人類HLA抗原、Megsin基因、RAS系統(tǒng)基因、MCP基因等的多態(tài)性位點與IgA腎病易感性相關(guān)[14]。在此僅探索Megsin基因與IgA腎病的相關(guān)性。

Megsin屬于Serpin超家族成員之一，有研究發(fā)現(xiàn)在IgA腎病患者中，其RNA及蛋白質(zhì)水平顯著上調(diào)，且發(fā)現(xiàn)Megsin表達(dá)與系膜細(xì)胞數(shù)量增多及系膜基質(zhì)增多有關(guān)。國內(nèi)外學(xué)者對Megsin rs1055901、rs1055902和rs2689399多態(tài)性位點與IgA腎病易感性進(jìn)行大量研究，部分學(xué)者對此進(jìn)行Meta分析，但是現(xiàn)有的關(guān)于Megsin rs1055901、rs1055902及rs2689399與IgA腎病易感性的Meta分析存在一定的缺陷。如Zhou等[15]在2015年完成的納入7篇文獻(xiàn)的Meta分析中，其中所納入的1篇與Megsin基因無關(guān)[16]，1篇數(shù)據(jù)不全[17]，還有1篇數(shù)據(jù)提取錯誤[6]，且上述Meta分析針對的是總體人群。

因此，本研究針對Megsin基因的3個位點的多態(tài)性位點與亞洲人群的IgA腎病易感性進(jìn)行更新Meta分析。結(jié)果發(fā)現(xiàn)，rs1055901和rs1055902多態(tài)性位點的5種模型(C和T，CC和TT，CT和TT，CT+CC和TT、CC和CT+TT)均與在亞洲人種IgA腎病易感性無關(guān)。但因研究及樣本量有限，未來仍需大量樣本去證實rs1055901和rs1055902多態(tài)性位點與IgA腎病易感性的關(guān)系。同時，本文也分析了rs2689399多態(tài)性位點與IgA腎病患病風(fēng)險的關(guān)系，發(fā)現(xiàn)在等位基因模型、純合子模型及隱性模型中，該基因多態(tài)性位點與IgA腎病易感性有統(tǒng)計學(xué)關(guān)聯(lián)，Megsin基因rs2689399的G等位基因和GG基因型可能是IgA腎病的一個保護(hù)因子。為了進(jìn)一步探討Megsin基因rs1055901、rs1055902和rs2689399之間的關(guān)系，通過Hapmap數(shù)據(jù)庫，利用Haploview4.0軟件，對Megsin基因連鎖不平衡分析，結(jié)果發(fā)現(xiàn)rs1055901、rs1055902和rs2689399在CEU人群中存在完全連鎖不平衡，在YRI人群中三者存在一定程度連鎖。但這個結(jié)果在亞洲人群中并未得到證實，這意味著更需要大量研究去證明這3個位點與亞洲人群IgA腎病易感性之間的關(guān)系。本篇Meta分析納入的數(shù)據(jù)嚴(yán)格遵循納入、剔除標(biāo)準(zhǔn)，病例和對照均較一致，且在年齡、性別、種族方面具有可比性。然而，本文也存在缺陷，如樣本量較小，不確定對照組有無篩查具有相關(guān)的亞臨床IgA腎病，且有一個研究的對照組不符合HWE檢驗，其多態(tài)性檢測方法也不一致。

綜上所述，本研究認(rèn)為Megsin基因rs2689399位點與IgA腎病易感性相關(guān)，rs2689399的G等位基因和GG基因型可能是IgA腎病的一個保護(hù)因子，但仍需要未來更多合理設(shè)計的大樣本研究進(jìn)一步證實本研究的結(jié)論。

[1]Donadio JV,Grande JP.IgA nephropathy[J].N Engl J Med,2002,347(10):738-748.

[2]Floege J,Amann K.Primary glomerulonephritides[J].Lancet,2016,387(10032):2036-2048.

[3]Schena FP.A retrospective analysis of the natural history of primary IgA nephropathy worldwide[J].Am J Med,1990,89(2):209-215.

[4]Wyatt RJ,Julian BA.IgA nephropathy[J].N Engl J Med,2013,368(25):2402-2414.

[5]陳楠,王偉銘,俞海瑾，等.腎小球腎炎腎組織megsin的表達(dá)及其意義[J].腎臟病與透析腎移植雜志,2003,12(5):412-414.

[6]Maixnerova D,Merta M,Reiterova J,et al.The influence of two megsin polymorphisms on the progression of IgA nephropathy[J].Folia Biol(Praha),2008,54(2):40-45.

[7]Lim CS,Kim SM,Oh YK,et al.Megsin 2093T-2180C haplotype at the 3′ untranslated region is associated with poor renal survival in Korean IgA nephropathy patients[J].Clin Nephrol,2008,70(2):101-109.

[8]李露,張?zhí)煸?崔麗瓊.漢族人群Megsin基因C-2093T多態(tài)性與IgA腎病的關(guān)系[J].醫(yī)學(xué)信息,2013(24):24-25.

[9]Wei LT,Fu RG,Gao J,et al.Association of Megsin gene variants with IgA nephropathy in northwest Chinese population:a STROBE-compliant observational study[J].Medicine(Baltimore),2016,95(6):2694-2699.

[10]夏運風(fēng),黃霜,李采霞，等.Megsin基因C25663G多態(tài)性與我國漢族人群IgA腎病的關(guān)系[J].中國中西醫(yī)結(jié)合腎病雜志,2006,7(2):91-93.

[11]李露,張?zhí)煸?崔麗瓊.漢族人群Megsin基因C-25663G多態(tài)性與IgA腎病的關(guān)系[J].中外醫(yī)療,2011,30(30):2-3.

[12]王德潤,解汝娟,禹程遠(yuǎn)，等.不同性別原發(fā)性IgA腎病患者臨床與病理分析[J].實用醫(yī)學(xué)雜志,2011,27(20):3744-3746.

[13]曾憲濤,劉慧,陳曦，等.Meta分析系列之四:觀察性研究的質(zhì)量評價工具[J].中國循證心血管醫(yī)學(xué)雜志,2012,4(4):297-299.

[14]賈妮亞,王彩麗.IgA腎病基因多態(tài)性研究進(jìn)展[J].中國中西醫(yī)結(jié)合腎病雜志,2009,10(9):833-835.

[15]Zhou TB,Mo XN,Kou QY,et al.Association of Megsin gene polymorphism with IgA nephropathy risk[J].J Recept Signal Transduct Res,2015,35(2):122-126.

[16]Szelestei T,Bahring S,Kovacs T,et al.Association of a uteroglobin polymorphism with rate of progression in patients with IgA nephropathy[J].Am J Kidney Dis,2000,36(3):468-473.

[17]Li YJ,Du Y,Li CX,et al.Family-based association study showing that immunoglobulin A nephropathy is associated with the polymorphisms 2093C and 2180T in the 3′ untranslated region of the Megsin gene[J].J Am Soc Nephrol,2004,15(7):1739-1743.

Meta-analysis of the relevance between Megsin rs1055901,rs1055902 and rs2689399 polymorphism and susceptibility of IgA nephrology in Asian population*

GeYating1，2,SuMeiling2,SongJiazhi1，2,XiongZuying1,2△,HouShuang2▲

(1.ClinicalCollege,ShenzhenHospital,AnhuiMedicalUniversity&PekingUniversity,Shenzhen,Guangdong518000,China; 2.DepartmentofNephrology,ShenzhenHospital,PekingUniversity,Shenzhen,Guangdong518000,China)

Objective To assess the association of three polymorphisms in Megsin (rs1055901,rs1055902 and rs2689399)and susceptibility of IgA nephropathy in Asian population.Methods We conducted a comprehensive search of electronic CNKI,VIP,WangFang Data,CBM,Pubmed,Web of Science and Google Scholar database on the association between Megsin rs1055901,rs1055902 and rs2689399 polymorphism and susceptibility of IgA nephrology in Asian population (last search update on 2 May 2016).Stata 12.0 software was used to calculate the odds ratio(OR)and 95%CI(confidence interval),as well as sensitivity and publication bias analyses.Results Six publications encompassing nine case-control studies were finally included,including 2 179 cases and 1 769 controls.Finally,no significant association between Megsin rs1055901 and rs1055902 polymorphism and IgA nephrology in Asian population was identified,while a significantly decreased risk of IgA nephrology for rs2689399 polymorphism,was identified in Asian population (G and C:OR=0.754,95%CI0.592-0.961,P=0.022;GG and CC:OR=0.506,95%CI0.287-0.892,P=0.019;GG and GC+CC:OR=0.551,95%CI0.316-0.961,P=0.036).Conclusion Rs2689399 G allele and GG genotype of Megsin may be the protective factors for IgA nephropathy in Asian population.

glomeruloneohritis,IgA;disease susceptibility;linkage disequilibrium;Megsin;polymorphism;Meta-analysis

2016年度重慶市出版專項資金資助項目

?證醫(yī)學(xué)·

10.3969/j.issn.1671-8348.2017.05.022

深圳市科技研發(fā)資金(JCYJ20140415162543010)。 作者簡介：葛雅婷(1989-)，在讀碩士，主要從事原發(fā)性腎小球疾病的發(fā)病機制的研究?！?/p>

，E-mail:xiongzy2005@163.com。▲共同通信作者，E-mail:hssz2007@126.com。

R692.3

A

1671-8348(2017)05-0648-03

2016-09-21

2016-11-22)