羅旭東，林舉達(dá)，鄒曉波，殷靜雯

(廣東醫(yī)學(xué)院附屬醫(yī)院精神心理科，廣東 湛江 524001)

·薈萃分析·

miR-137基因rs1625579多態(tài)性與精神分裂癥發(fā)病風(fēng)險(xiǎn)的Meta分析

羅旭東，林舉達(dá)，鄒曉波，殷靜雯

(廣東醫(yī)學(xué)院附屬醫(yī)院精神心理科，廣東 湛江 524001)

目的 系統(tǒng)評(píng)價(jià)miR-137 rs1625579位點(diǎn)基因多態(tài)性與精神分裂癥發(fā)病風(fēng)險(xiǎn)的關(guān)系。方法計(jì)算機(jī)檢索PubMed、Embase、Springer Link、Wiley Online Library、HighWire Press、Google Scholar和CNKI數(shù)據(jù)庫(kù)以獲取2015年8月之前發(fā)表的關(guān)于miR-137單核苷酸多態(tài)位點(diǎn)rs1625579基因多態(tài)性與精神分裂癥發(fā)病的病例對(duì)照研究文獻(xiàn)，采用Review Manager5.0和Stata12.0統(tǒng)計(jì)軟件進(jìn)行統(tǒng)計(jì)分析。結(jié)果最終納入8篇，miR-137基因rs1625579位點(diǎn)多態(tài)性與精神分裂癥在顯性模型(TT+GT vs GG)、隱性模型(TT vs GT+GG)、加性模型(TT vs GG)和等位基因模型(T vs G)分析中差異有統(tǒng)計(jì)學(xué)意義(P＜0.05)。合并OR值分別為1.47、1.21、1.53、1.17；P值分別為0.03、0.001、0.02、0.003。結(jié)論miR-137單核苷酸多態(tài)性位點(diǎn)rs1625579與精神分裂癥的易感性均具顯著相關(guān)，為miR-137多態(tài)性與精神分裂癥關(guān)聯(lián)性的研究提供了理論參考。

精神分裂癥；miR-137；rs1625579；多態(tài)性；Meta分析

精神分裂癥是神經(jīng)發(fā)育異常引起的復(fù)雜遺傳疾病[1]。微小RNA(microRNA，miRNA)通過(guò)多條信號(hào)通路調(diào)控多個(gè)基因表達(dá)對(duì)神經(jīng)發(fā)育起關(guān)鍵作用[2]。微小RNA-137(microRNA-137，miR-137)是在大腦中高豐度表達(dá)的miRNA，它對(duì)胚胎神經(jīng)干細(xì)胞決定、神經(jīng)細(xì)胞增殖與分化、突觸成熟等方面都具有重要作用，它的失調(diào)會(huì)引起神經(jīng)系統(tǒng)基因表達(dá)調(diào)控網(wǎng)絡(luò)的改變而誘發(fā)精神相關(guān)疾病，它已被確立為精神分裂癥易感基因[3]。全基因組關(guān)聯(lián)性掃描研究發(fā)現(xiàn)位于miR-137宿主基因上的單核苷酸多態(tài)性位點(diǎn)rs1625579，與精神分裂癥具有高相關(guān)性[4]。但是在不同人群中進(jìn)行的病例對(duì)照研究卻得出不一致的結(jié)果[5-9]，因此，本研究采用Meta分析的方法系統(tǒng)評(píng)價(jià)miR-137 rs1625579位點(diǎn)基因多態(tài)性與精神分裂癥發(fā)病風(fēng)險(xiǎn)的關(guān)系。

1 資料與方法

1.1 檢索策略 從PubMed(http://www.ncbi.nlm. nih.gov/pubmed)，Embase(http://www.elsevier.com/solutions/embase-biomedical-research)，Springer Link(http: //link.springer.com/)，Wiley Online Library(http://onlinelibrary.wiley.com/)，HighWire Press(http://highwire. stanford.edu)，Google Scholar(http://scholar.google. com)和 China NationalKnowledge Infrastructure (CNKI:http://www.cnki.net)數(shù)據(jù)庫(kù)以獲取2015年8月之前發(fā)表的關(guān)于miR-137多態(tài)性與精神分裂癥易感性關(guān)系的所有病例對(duì)照研究，使用檢索關(guān)鍵詞及主題詞包括：schizophrenia；miR-137；microRNA-137；rs1625579。

1.2 納入標(biāo)準(zhǔn) (1)語(yǔ)種為英語(yǔ)或中文；(2)病例組符合診斷標(biāo)準(zhǔn)(DSM-Ⅳ或者ICD-10)的精神分裂癥患者，對(duì)照組為健康人群；(3)病例組和對(duì)照組有明確的rs1625579基因型頻率分布數(shù)據(jù)，可以計(jì)算比值比(Odds ratio，OR)和95%可信區(qū)間(Confidence inter-val,CI)；(4)篩除數(shù)據(jù)不完整的文獻(xiàn)。

1.3 資料提取 對(duì)納入研究的文獻(xiàn)資料提取第一作者姓名、發(fā)表日期、研究對(duì)象種族來(lái)源、病例組和對(duì)照組的總?cè)藬?shù)、哈迪-溫伯格平衡(Hardy-Weinberg equilibrium,HWE)檢驗(yàn)結(jié)果、基因分型方法、基因型及等位基因頻率分布數(shù)據(jù)。文獻(xiàn)檢索、資料提取時(shí)由兩位研究者分別進(jìn)行，共同商議、匯總。無(wú)法達(dá)成一致，交由第3位研究者裁定。納入文獻(xiàn)采用NOS量表(Newcastle-Ottawa Scale)進(jìn)行質(zhì)量評(píng)價(jià)。

1.4 資料合并及統(tǒng)計(jì)學(xué)分析 采用OR作為檢驗(yàn)統(tǒng)計(jì)量和95%CI作為衡量miR-137多態(tài)性rs1625579與精神分裂癥發(fā)病風(fēng)險(xiǎn)之間關(guān)聯(lián)程度的效應(yīng)指標(biāo)。合并OR值時(shí)采用下述4種模型進(jìn)行計(jì)算：顯性模型(TT+GT vs GG)、隱性模型(TT vs GT+GG)、加性模型(TT vs GG)、等位基因模型(T vs G)。采用基于χ2檢驗(yàn)的校正Q檢驗(yàn)來(lái)衡量各研究間效應(yīng)量分布有無(wú)顯著異質(zhì)性，檢驗(yàn)水準(zhǔn)設(shè)為0.05，采用I2值評(píng)價(jià)各研究間變異占總變異的百分比，I2值在0～1之間，I2=0時(shí)各研究間無(wú)異質(zhì)性，I2值趨近于0時(shí)各研究間異質(zhì)性可歸因于偶然因素，低于25%認(rèn)為是輕度異質(zhì)性，25%～50%為中度異質(zhì)性，大于50%則認(rèn)為各研究間有高度異質(zhì)性。若P＞0.05或I2＜50%，各研究間效應(yīng)量分布異質(zhì)性顯著，采用固定效應(yīng)模型(Mantel-Haenszel)進(jìn)行合并，反之，則采用隨機(jī)效應(yīng)模型(DerSimonian-Laird)進(jìn)行合并運(yùn)算。采用Z檢驗(yàn)評(píng)價(jià)合并OR值有無(wú)統(tǒng)計(jì)學(xué)意義，檢驗(yàn)水準(zhǔn)設(shè)為0.05。根據(jù)Egger′s回歸檢驗(yàn)、Begg′s級(jí)聯(lián)檢驗(yàn)判斷納入文獻(xiàn)是否存在顯著發(fā)表偏倚。偏倚檢驗(yàn)P＜0.01則認(rèn)為存在顯著發(fā)表偏倚。所有P值均為雙側(cè)檢驗(yàn)，除異質(zhì)性評(píng)估及發(fā)表偏倚外，檢驗(yàn)水準(zhǔn)均設(shè)為P＜0.05時(shí)差異具有統(tǒng)計(jì)學(xué)意義。統(tǒng)計(jì)分析采用STATA 10.0(Stata Corporation，College Station，TX，USA)和Review Manager 5.0(Coehrane Collaboration,Software Update，Oxford，United Kingdom)軟件完成。

2 結(jié) 果

2.1 文獻(xiàn)檢索結(jié)果 共檢出相關(guān)文獻(xiàn)28篇，閱讀題目、摘要和全文，根據(jù)納入和排除標(biāo)準(zhǔn)以及結(jié)合文獻(xiàn)質(zhì)量評(píng)價(jià)結(jié)果最終對(duì)文獻(xiàn)做如下處理：剔除綜述1篇，非基因多態(tài)研究10篇，非精神分裂癥相關(guān)研究3篇，未報(bào)告rs1625579多態(tài)位點(diǎn)文獻(xiàn)1篇，未進(jìn)行基因型及等位基因分析3篇，非病例對(duì)照研究1篇，無(wú)完整基因型數(shù)據(jù)1篇，最終納入8篇文獻(xiàn)進(jìn)行Meta分析(圖1)。

2.2 納入文獻(xiàn)概要 依據(jù)納入、排除標(biāo)準(zhǔn)和NOS量表(Newcastle-Ottawa Scale)評(píng)價(jià)結(jié)果，最終納入8篇病例對(duì)照研究文獻(xiàn)，其中病例組包括3 802例精神分裂癥患者，對(duì)照組包括3 461名健康對(duì)照，納入文獻(xiàn)中的樣本均符合H-W遺傳平衡定律，見(jiàn)表1。

圖1 文獻(xiàn)納入和排除標(biāo)準(zhǔn)流程

表1 納入miR-137基因多態(tài)性與精神分裂癥病例對(duì)照研究基本信息

2.3 Meta分析結(jié)果 miR-137基因rs1625579位點(diǎn)多態(tài)性與精神分裂癥顯著相關(guān)性，見(jiàn)圖2。顯性模型(TT+GT vs GG)分析中，突變率異質(zhì)性檢驗(yàn)P=0.83 (P＞0.05)，I2=0%，各個(gè)文獻(xiàn)間不存在異質(zhì)性，故選用固定效應(yīng)模型進(jìn)行統(tǒng)計(jì)合并，OR合并值為1.47，95%CI為(1.03 ，2.11)，P=0.03。結(jié)果表明TT+GT基因型與GG基因型比較，rs1625579基因型多態(tài)性與精神分裂癥顯著相關(guān)(P＜0.05)。隱性模型(TT vs GT+GG)分析中，突變率異質(zhì)性檢驗(yàn)P=0.84(P＞0.05)，I2=0%，各個(gè)文獻(xiàn)間不存在異質(zhì)性，故選用固定效應(yīng)模型進(jìn)行統(tǒng)計(jì)合并，OR合并值為1.21，95%CI為(1.08，1.35)，P= 0.001。結(jié)果表明TT基因型與GT+GG基因型比較，rs1625579基因型多態(tài)性與精神分裂癥顯著相關(guān)(P＜0.05)。加性模型(TT vs GG)分析中，突變率異質(zhì)性檢驗(yàn)P=0.82(P＞0.05)，I2=0%，各個(gè)文獻(xiàn)間不存在異質(zhì)性，故選用固定效應(yīng)模型進(jìn)行統(tǒng)計(jì)合并，OR合并值為1.53，95%CI為(1.06，2.19），P=0.02。結(jié)果表明TT基因型與GG基因型比較rs1625579基因型多態(tài)性與精神分裂癥顯著相關(guān)(P＜0.05)。等位基因模型(T vs G)分析中，突變率異質(zhì)性檢驗(yàn)P=0.37(P＞0.05)，I2=7%，各個(gè)文獻(xiàn)間異質(zhì)性可歸因于偶然因素，故選用固定效應(yīng)模型進(jìn)行統(tǒng)計(jì)合并，OR合并值為1.17，95%CI為(1.06，1.30)，P=0.003。結(jié)果表明TT+GT基因型與GG基因型比較，rs1625579基因型多態(tài)性與精神分裂癥顯著相關(guān)(P＜0.05)。

2.4 發(fā)表偏倚分析 顯性模型(TT+GT vs GG)中，Begg′s檢測(cè)Z=0.00，P=1.00；Egger′s檢測(cè)t=0.17，P= 0.869。隱性模型(TT vs GT+GG)中，Begg’s檢測(cè)Z= 0.37，P=0.711；Egger′s檢測(cè)t=0.66，P=0.531。加性模型(TT vs GG)中，Begg′s檢測(cè)Z=0.12，P=1.00；Egger′s檢測(cè)t=0.25，P=0.811。等位基因模型(T vs G)分析中，Begg′s檢測(cè)Z=0.62，P=0.536；Egger′s檢測(cè)t=0.50，P=0.636。研究結(jié)果表明：本研究納入的文獻(xiàn)不存在發(fā)表偏倚。

圖2 miR-137基因多態(tài)性與精神分裂癥發(fā)病風(fēng)險(xiǎn)的Meta分析森林圖

3 討 論

近年來(lái)，miR-137對(duì)神經(jīng)細(xì)胞增殖和分化及對(duì)精神分裂癥的相關(guān)基因的調(diào)控可能使其成為精神分裂癥發(fā)病的易感基因[3]。在尸檢中發(fā)現(xiàn)miR-137在前額葉皮質(zhì)的表達(dá)在患者和對(duì)照組之間沒(méi)有顯著差異，但是在40名正常人群中，TT風(fēng)險(xiǎn)基因型的被試前額葉皮質(zhì)miR-137表達(dá)量顯著高于G等位基因攜帶者，精神分裂癥患者中則差異無(wú)統(tǒng)計(jì)學(xué)意義[13]。該研究首次提出了miR-137的表達(dá)量在正常人和精神分裂癥患者中都存在差異性的變化，特別是單核苷酸多態(tài)位點(diǎn)rs1625579的不同基因型miR-137表達(dá)量的差異，因此，提示rs1625579可能參與調(diào)控miR-137的表達(dá)，而與精神分裂癥的易感性相關(guān)。

rs1625579位于miR-137基因下游的8 kb區(qū)域，是迄今發(fā)現(xiàn)與精神分裂癥關(guān)聯(lián)性最強(qiáng)位點(diǎn)(P=1.6×10-11)[4]。但是，隨后在英國(guó)[8，14]、澳大利亞[5]、美國(guó)、加拿大[10]等高加索人群的相關(guān)研究中，miR-137與精神分裂癥病例對(duì)照研究得出了陰性的結(jié)果，在中國(guó)人群中rs1625579與精神分裂癥相關(guān)性的研究結(jié)果不一致[3，6，11]。差異的結(jié)果使miR-137在精神分裂癥中易感性的研究遇到了一定的困難。造成差異性結(jié)果的原因可能為：首先，由于精神分裂癥是一個(gè)高度異質(zhì)性的疾病，目前尚無(wú)客觀的實(shí)驗(yàn)室檢查作為診斷依據(jù)，造成不同研究的入組病例異質(zhì)性增高；其次，受到樣本量的影響，即使是對(duì)于同一多態(tài)位點(diǎn)在相同人群中的研究也得出不同甚至相反的結(jié)果；再次，由于在不同種族人群中等位基因頻率、單體型頻率不同。造成病例對(duì)照研究常出現(xiàn)矛盾性的結(jié)果。以rs1625579多態(tài)位點(diǎn)為例，筆者的檢索結(jié)果顯示rs1625579三種基因型頻率分布在亞洲人群、歐洲人群和非洲人群之間存在顯著的差異；同屬亞洲人群的日本和中國(guó)人群也具有差異性：中國(guó)人群的TT基因型的頻率最低，而非洲人群的TT基因型的頻率最高。我們查閱了miR-137 rs1625579多態(tài)性與精神分裂癥發(fā)病風(fēng)險(xiǎn)的病例對(duì)照研究進(jìn)行Meta分析，最終納入符合要求的8篇研究，合計(jì)患者3 802例，健康對(duì)照者3 461例。Meta分析的結(jié)果提示在顯性模型、隱性模型、加性模型和等位基因模型中miR-137單核苷酸多態(tài)性位點(diǎn)rs1625579與精神分裂癥的易感性均具顯著相關(guān)(分別為P=0.03，P=0.001，P=0.02，P=0.003)。進(jìn)一步證實(shí)了miR-137基因多態(tài)與精神分裂癥的相關(guān)性。

雖然Meta定量合成多個(gè)研究的結(jié)果，能夠提高檢驗(yàn)效能。但是仍存在一些可能影響結(jié)果的局限性：在研究中存在語(yǔ)言偏倚，檢索語(yǔ)種限定為中文和英文，排除了其他的語(yǔ)種，可能造成結(jié)果的偏倚。

綜上所述，miR-137單核苷酸多態(tài)性位點(diǎn)rs1625579與精神分裂癥的易感性均具顯著相關(guān)。我們的研究結(jié)果為miR-137多態(tài)性與精神分裂癥關(guān)聯(lián)性的研究提供理論參考。但是rs1625579基因多態(tài)如何影響miR-137的表達(dá)進(jìn)而調(diào)控其下游通路影響精神分裂癥的機(jī)制目前尚不清楚。rs1625579多態(tài)到底是發(fā)揮功能性作用，或者僅僅是與功能性多態(tài)位點(diǎn)相連鎖的遺傳標(biāo)記，這都需要進(jìn)一步研究證實(shí)。

[1]Cannon TD.Abnormalities of brain structure and function in schizophrenia:implications for aetiology and pathophysiology[J].Ann Med,1996,28(6):533-539.

[2]Miller BH,Wahlestedt C.MicroRNA dysregulation in psychiatric disease[J].Brain Research,2010,1338:89-99.

[3]Yin J,Lin J,Luo X,et al.miR-137:a new player in schizophrenia[J]. Int J Mol Sci,2014,15(2):3262-3271.

[4]Ripke S,Sanders A,Kendler K,et al.Genome-wide association study identifies five new schizophrenia loci[J].Nat Genet,2011,43(10): 969-976.

[5]Green MJ,Cairns MJ,Wu J,et al.Genome-wide supported variant MIR137 and severe negative symptoms predict membership of an impaired cognitive subtype of schizophrenia[J].Molecular Psychiatry, 2013,18(7):774-780.

[6]Guan F,Zhang B,Yan T,et al.MIR137 gene and target gene CACNA1C of miR-137 contribute to schizophrenia susceptibility in Han Chinese[J].Schizophr Res,2014,152(1):97-104.

[7]van Erp TG,Guella I,Vawter MP,et al.Schizophrenia miR-137 Locus Risk Genotype is Associated with Dorsolateral Prefrontal Cortex Hyperactivation[J].Biol Psychiatry,2014,75(5):398-405.

[8]Whalley HC,Papmeyer M,Romaniuk L,et al.Impact of a microRNA MIR137 Susceptibility Variant on Brain Function in People at High Genetic Risk of Schizophrenia or Bipolar Disorder[J].Neuropsychopharmacology,2012,37(12):2720-2729.

[9]Ma G,Yin J,Fu J,et al.Association of a miRNA-137 polymorphism with schizophrenia in a Southern Chinese Han population[J]. BioMed Research International,2014,2014:751267.

[10]Lett TA,Chakravarty MM,Felsky D,et al.The genome-wide supported microRNA-137 variant predicts phenotypic heterogeneity within schizophrenia[J].Mol Psychiatry,2013,18(4):443-450.

[11]Rose EJ,Morris DW,Fahey C,et al.The miR-137 schizophrenia susceptibility variant rs1625579 does not predict variability in brain volume in a sample of schizophrenic patients and healthy individuals[J]. Am J Med Genet B Neuropsychiatr Genet,2014,165B(6):467-471.

[12]Yuan J,Cheng Z,Zhang F,et al.Lack of association between microRNA-137 SNP rs1625579 and schizophrenia in a replication study of Han Chinese[J].Mol Genet Genomics,2015,290(1):297-301.

[13]Guella I,Sequeira A,Rollins B,et al.Analysis of miR-137 expression and rs1625579 in dorsolateral prefrontal cortex[J].J Psychiatr Res,2013,47(9):1215-1221.

[14]Cummings E,Donohoe G,Hargreaves A,et al.Mood congruent psychotic symptoms and specific cognitive deficits in carriers of the novel schizophrenia risk variant at MIR-137[J].Neurosci Lett,2013, 532:33-38.

Correlation between single nucleotide polymorphism rs1625579 in miR-137 gene and schizophrenia:a meta-analysis.

LUO Xu-dong,LIN Ju-da,ZOU Xiao-bo,YIN Jing-wen.Department of Psychiatry,the Affiliated Hospital of Guangdong Medical College,Zhanjiang 524001,Guangdong,CHINA

ObjectiveTo assay the association between the single nucleotide polymorphism(rs1625579)in miR-137 gene and schizophrenia.MethodsStudies about the effects of rs1625579 in miR-137 gene on the sensitivity in patients with schizophrenia until Aug.2015 were searched from PubMed,Embase,Springer Link,Wiley Online Library,HighWire Press,Google Scholar,and China National Knowledge Infrastructure(CNKI).Meta-analysis was performed by using Review Manager5.0 and Stata12.0.ResultsA total of 8 studies were included.A significant association between rs1625579 and schizophrenia was found using dominant[TT+GT vs GG,P=0.03,odds ratio(OR)=1.47],recessive(TT vs GT+GG,P=0.001,OR=1.21),additive(TT vs GG,P=0.02,OR=1.53),and allele(T vs G,P=0.003,OR= 1.17)models.ConclusionThe results of our study provide convincing evidence that rs1625579 is significantly associated with schizophrenia.

Schizophrenia;miR-137;rs1625579;Polymorphism;Meta-analysis

R749.3

A

1003—6350（2016）04—0656—04

10.3969/j.issn.1003-6350.2016.04.051

2015-10-22)

廣東省湛江市科技攻關(guān)計(jì)劃項(xiàng)目(編號(hào)：2013B01258)

通迅作者：殷靜雯。E-mail：amysays@live.com