miR-137基因rs1625579多態性與精神分裂癥發病風險的Meta分析

羅旭東，林舉達，鄒曉波，殷靜雯

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

·薈萃分析·

miR-137基因rs1625579多態性與精神分裂癥發病風險的Meta分析

羅旭東，林舉達，鄒曉波，殷靜雯

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

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

精神分裂癥；miR-137；rs1625579；多態性；Meta分析

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

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)數據庫以獲取2015年8月之前發表的關于miR-137多態性與精神分裂癥易感性關系的所有病例對照研究，使用檢索關鍵詞及主題詞包括：schizophrenia；miR-137；microRNA-137；rs1625579。

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

1.3 資料提取 對納入研究的文獻資料提取第一作者姓名、發表日期、研究對象種族來源、病例組和對照組的總人數、哈迪-溫伯格平衡(Hardy-Weinberg equilibrium,HWE)檢驗結果、基因分型方法、基因型及等位基因頻率分布數據。文獻檢索、資料提取時由兩位研究者分別進行，共同商議、匯總。無法達成一致，交由第3位研究者裁定。納入文獻采用NOS量表(Newcastle-Ottawa Scale)進行質量評價。

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

2 結 果

2.1 文獻檢索結果 共檢出相關文獻28篇，閱讀題目、摘要和全文，根據納入和排除標準以及結合文獻質量評價結果最終對文獻做如下處理：剔除綜述1篇，非基因多態研究10篇，非精神分裂癥相關研究3篇，未報告rs1625579多態位點文獻1篇，未進行基因型及等位基因分析3篇，非病例對照研究1篇，無完整基因型數據1篇，最終納入8篇文獻進行Meta分析(圖1)。

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

圖1 文獻納入和排除標準流程

表1 納入miR-137基因多態性與精神分裂癥病例對照研究基本信息

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

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

圖2 miR-137基因多態性與精神分裂癥發病風險的Meta分析森林圖

3 討 論

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

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

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

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

[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)

廣東省湛江市科技攻關計劃項目(編號：2013B01258)

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