李 娜，安 宇，張永錄，丁志杰，吳柏林，蘭小平

1. 隴東學院生命科學與技術(shù)學院，慶陽 745000；2. 復(fù)旦大學人類表型組研究院，上海 201203；3. 甘肅省天水市第三人民醫(yī)院精神科，天水741000；4. 復(fù)旦大學生物醫(yī)學研究院，上海200032；5. 天水師范學院生物工程與技術(shù)學院，天水741000；6. 上海市兒童醫(yī)院，上海交通大學附屬兒童醫(yī)院檢驗科，上海200040

精神分裂癥（schizophrenia，SZ）是一類病因未明的精神類疾病，以妄想、幻想及認知障礙為典型特征，終身患病率為1%[1]，遺傳度高達79%[2]。目前認為SZ具有高度遺傳異質(zhì)性，不同性別SZ患者的發(fā)病年齡和SZ亞型不同[3-4]。

隨著分子遺傳學技術(shù)的發(fā)展，已有大量研究發(fā)現(xiàn)γ-氨基丁酸（γ-aminobutyric acid，GABA）與 SZ的發(fā)病具有相關(guān)性。GABA是中樞神經(jīng)系統(tǒng)中主要的抑制性神經(jīng)遞質(zhì)，由腦內(nèi)興奮性神經(jīng)遞質(zhì)谷氨酸在谷氨酸脫羧酶（glutamicacid decarboxylase，GAD）的催化下生成。研究[5-7]發(fā)現(xiàn)SZ患者死亡后腦組織中GABA活性水平及谷氨酸脫羧酶 67（glutamate decarboxylase 67，GAD67）亞型的mRNA的表達量下降。GAD67由位于2q31.1區(qū)的谷氨酸脫羧酶1（glutamate decarboxylase 1，GAD1）基因編碼，該基因包含17個外顯子。除了合成GABA的限速酶GAD1基因表達量改變引起GABA活性水平下降外，GABA受體也參與SZ的發(fā)病。γ-氨基丁酸A受體β3亞基 [γ-aminobutyric acid （GABA）Areceptor β-3，GABRB3]基因、γ-氨基丁酸 A 受體 α5亞基 [γ-aminobutyric acid（GABA）Areceptor α-5，GABRA5] 基因和 γ-氨基丁酸 A受體γ3亞基[γ-aminobutyric acid （GABA）Areceptor γ-3，GABRG3] 基因構(gòu)成 γ-氨基丁酸 A 受體 [γ-aminobutyric acid （GABA）Areceptor，GABAA]的基因簇，分別編碼 GABAA受體的 β3、α5和 γ3亞基。已有研究[8-10]表明GABRB3基因是SZ的易感基因。然而，以往的研究所選的相關(guān)基因的單核苷酸多態(tài)性（single nucleotide polymorphisms，SNPs）位點大多數(shù)來自于基因編碼區(qū)，主要關(guān)注因氨基酸序列改變導(dǎo)致蛋白質(zhì)功能異常對SZ的影響，很少對基因非編碼區(qū)SNPs位點進行研究。本文選取GAD1和GABRB3基因啟動子區(qū)SNPs位點為研究對象，對中國漢族人群中GAD1基因rs3791878、rs3749034位點和GABRB3基因rs4906902位點的多態(tài)性與SZ的相關(guān)性進行分析。統(tǒng)計手冊》（第 4 版）（Diagnostic and Statistical Manual of Mental Disorders，fourth edition，DSM-Ⅳ）診斷標準明確診斷，且樣本間不存在親緣關(guān)系。624例正常個體樣本（對照組）來自天水市健康志愿獻血者，其中男性350名，女性274名，平均年齡（52.22±19.24）歲，所有正常個體均由2名以上精神科副主任醫(yī)師確認無精神類疾病和陽性家族史，樣本間不存在親緣關(guān)系。

1 對象與方法

1.1 研究對象

本研究所有對象均來源于中國漢族人群，545例SZ患者樣本（病例組）來自2012年10月至2013年10月在天水市第三人民醫(yī)院就診的患者，其中男性310例，女性235例，平均年齡（37.35±12.82）歲。所有SZ患者均由2名以上精神科副主任醫(yī)師參照《美國精神障礙診斷與

1.2 研究方法

1.2.1 基因選擇、啟動子區(qū)預(yù)測及Tag-SNPs的選擇 根據(jù)已有文獻報道[11-12]和蛋白質(zhì)互作（protein-protein interaction，PPI）在線數(shù)據(jù)庫STRING （http://string-db.org/）的分析結(jié)果，選擇了在GABA通路中存在互作關(guān)系的GAD1和GABRB3基因作為目標基因。通過序列保守性分析預(yù)測了GAD1和GABRB3基因的啟動子區(qū)域，結(jié)合1000 GENOMES數(shù)據(jù)庫中的中國北京漢族人群（Han Chinese in Beijing，CHB）和中國南方漢族人群（Southern Han Chinese，CHS）（CHB+CHS=197）的 SNPs數(shù)據(jù)（dbSNP137），利用Haploview 4.2軟件，以條件決定系數(shù)r2＞0.8、最小等位基因頻率（minor allele frequency，MAF） ＞0.05為標準進行選擇，共篩選出 3個標簽SNPs（Tag-SNPs），分別是位于GAD1基因啟動子區(qū)的rs3791878和rs3749034位點及位于GABRB3基因啟動子區(qū)的rs4906902位點。

1.2.2 樣本采集及基因組DNA提取 本研究經(jīng)天水市第三人民醫(yī)院倫理委員會批準。受試者均由本人或監(jiān)護人簽署知情同意書后，采集EDTA抗凝外周靜脈血3～5 mL，采用QIAamp DNA Blood Mini Kit（QIAGEN，德國）試劑盒提取基因組DNA。

1.2.3 SNaPshot基因分型 采用SNaPshot分型技術(shù)檢測Tag-SNPs的多態(tài)性。在T100TMPCR 儀（Bio-Rad，美國）上 利 用 Qiagen multiplex PCR kit（100）（QIAGEN， 德國）試劑構(gòu)建多重PCR反應(yīng)體系，50 μL PCR反應(yīng)物 包含 0.2 μmol/L PCR 引物（表 1）和 1 μL 模板 DNA，PCR反應(yīng)條件為：97 ℃ 5 min，35個循環(huán)，95 ℃ 30 s，58 ℃90 s，72 ℃ 60 s，72℃ 10 min。經(jīng)純化和 SNaPshot單堿基延伸反應(yīng)后，利用ABI 3130 DNA 測序儀（Applied Biosystems，美國）檢測等位基因分布。SNaPshot試劑盒為 SNaPshot master mix（Applied Biosystems，美國）。應(yīng)用Soft Genetics Gene Marker v 2.2.0 軟件進行基因型分析。隨機抽取 10%的樣本使用 ABI 3130 DNA 測序儀進行雙向Sanger測序驗證，基因分型一致率為100%。

表1 SNaPshot分型引物Tab 1 Primers for SNaPshot genotyping

1.3 統(tǒng)計學分析

使用 R64 3.5.3 軟件分析病例組與對照組間年齡和性別的差異，年齡分析采用兩獨立樣本t檢驗、性別分析采用χ2檢驗。應(yīng)用SNPstats（https://www.snpstats.net/snpstats/start.htm）在線軟件進行哈迪 - 溫伯格平衡 （Hardy-Weinberg Equilibrium，HWE） 檢驗、等位基因和不同遺傳模式下基因型分布分析。采用Kaplan-Meier統(tǒng)計分析檢驗rs3791878位點風險基因型G/T與男性SZ患者發(fā)病年齡的相關(guān)性。樣本量統(tǒng)計效能檢驗使用Quanto 1.2.4軟件計算。所有檢驗均為雙側(cè)檢驗，P＜0.05 為差異有統(tǒng)計學意義。

2 結(jié)果

2.1 一般人口學資料

病例組共545例，男性310例（56.9%），女性235例（43.1%），年齡14～76歲，平均年齡（37.35±12.82）歲。對照組共624例，男性350例（56.1%），女性274例（43.9%），年齡19～88歲，平均年齡（52.22±19.24）歲。病例組與對照組相比，性別構(gòu)成差異無統(tǒng)計學意義（P=0.831），年齡差異有統(tǒng)計學意義（P=0.000）。

2.2 統(tǒng)計效能檢驗

Tag-SNPs中rs3749034 的MAF最?。∕AF=0.233）。假設(shè)OR=1.5，加性模式，疾病一般人群患病率為1%[1]，根據(jù)Quanto1.2.4軟件計算，結(jié)果顯示統(tǒng)計效能為0.99。

2.3 Tag-SNPs等位基因分布比較

對照組中3個Tag-SNPs位點基因型分布符合HWE檢驗。由于樣本損失，有18例樣本的rs3791878位點數(shù)據(jù)缺失，12例樣本的rs3749034位點數(shù)據(jù)缺失，3例樣本的rs4906902位點數(shù)據(jù)缺失。分別在總體樣本和性別分層后樣本的病例組與對照組之間對3個Tag-SNPs位點的等位基因分布進行分析。結(jié)果顯示：總體樣本病例組與對照組相比，3個Tag-SNPs位點等位基因分布差異均無統(tǒng)計學意義（P＞0.05）（表2）；男性病例組與對照組相比，rs3791878位點等位基因分布差異有統(tǒng)計學意義（P=0.039），P值經(jīng) Bonferroni校正后，差異無統(tǒng)計學意義（P=0.117）；女性病例組與對照組相比，3個Tag-SNPs位點等位基因分布差異均無統(tǒng)計學意義（P＞0.05）（表 3）。

表2 病例組與對照組Tag-SNPs 等位基因分布的比較Tab 2 Comparison of allele distribution of Tag-SNPs between case group and control group

表3 病例組與對照組Tag-SNPs等位基因分布的比較（按性別分層）Tab 3 Comparison of allele distribution of Tag-SNPs between case group and control group （stratified by sex）

2.4 Tag-SNPs基因型分布比較

分別在共顯性（co-dominant）、顯性（dominant）、隱性（recessive）、超顯性（over-dominant）和加性化（logadditive）5種不同的遺傳模式下對3個Tag-SNPs位點基因型分布進行分析。結(jié)果顯示，總體樣本的病例組與對照組相比，rs3791878位點在共顯性、顯性和超顯性遺傳模式下基因型分布差異均有統(tǒng)計學意義（P=0.022，P=0.011，P=0.006）。挑選赤池信息準則（Akaike Information Criterion，AIC）和貝葉斯信息準則（Bayesian Information Criterion，BIC）數(shù)值最小的超顯性遺傳模式為rs3791878位點可能的遺傳模型，該模型下的P值經(jīng)Bonferroni 校正后，差異仍有統(tǒng)計學意義（P=0.018）。而rs3749034位點和rs4906902位點在5種不同遺傳模式下基因型分布差異均無統(tǒng)計學意義（P＞0.05）（表4）。在性別分層后樣本的病例組與對照組之間對rs3791878位點的基因型分布進行分析，結(jié)果顯示：男性病例組與對照組相比，rs3791878位點在超顯性遺傳模式下基因型分布差異有統(tǒng)計學意義（P=0.000），P值經(jīng)Bonferroni 校正后，差異仍有統(tǒng)計學意義（P=0.000）；女性病例組與對照組相比，rs3791878位點在5種不同遺傳模式下基因型分布差異均無統(tǒng)計學意義（P＞0.05）（表 5）。

表4 病例組與對照組 Tag-SNPs 不同遺傳模式下基因型分布的比較Tab 4 Comparison of genotype distribution of Tag-SNPs under different genetic models between case group and control group

Continued Tab

表5 病例組與對照組rs3791878位點不同遺傳模式下基因型分布的比較（按性別分層）Tab 5 Comparison of rs3791878 genotype distribution under different genetic models between case group and control group (stratified by sex)

Continued Tab

2.5 rs 3791878 位點風險基因型G/T與男性SZ發(fā)病年齡的相關(guān)性

310例男性SZ患者共有195例記錄了首次發(fā)病年齡，應(yīng)用Kaplan-Meier分析rs3791878位點風險基因型G/T與195例男性患者發(fā)病年齡的相關(guān)性，結(jié)果表明G/T基因型與男性SZ的發(fā)病年齡無相關(guān)性（χ2=0.271，P=0.603）（圖 1）。

圖1 rs3791878位點風險基因型G/T與男性SZ發(fā)病年齡的相關(guān)性Fig 1 Correlation between the risk genotype G/T of rs3791878 and the age of onset in male SZ

3 討論

本研究選擇GABA通路相關(guān)基因啟動子區(qū)SNPs為研究對象，探索GAD1基因rs3791878和rs3749034位點、GABRB3基因 rs4906902位點多態(tài)性與SZ的相關(guān)性。PPI分析顯示GAD1和GABRB3存在相互作用，且均與溶質(zhì)載體家族蛋白32成員A1（solute carrier family 32 member 1，SLC32A1）和突觸關(guān)聯(lián)蛋白25（synaptosomalassociated protein 25，SNAP25）有直接聯(lián)系；SLC32A1和SNAP25為GABA通路相關(guān)蛋白且與SZ發(fā)病相關(guān)[13-14]。關(guān)聯(lián)分析結(jié)果表明，rs3791878位點在超顯性遺傳模式下，G/T基因型會增加男性SZ的發(fā)病風險，且風險達1.91，但沒有得到G/T基因型提前男性SZ患者發(fā)病年齡的結(jié)果。由于本研究樣本量所限，這一結(jié)果還需要進一步在更大樣本量群體中進行驗證。本研究沒有發(fā)現(xiàn)rs3749034、rs4906902位點與SZ存在相關(guān)性。

Du等[15]在2008年對235個中國漢族SZ家系的研究結(jié)果也表明GAD1基因rs3791878位點與男性SZ的發(fā)生相關(guān)。結(jié)合本研究，我們認為rs3791878位點存在偏從性顯性遺傳的可能性，其雜合子的表達受到性別影響，從而增加SZ的發(fā)病風險。rs3791878位點在超顯性遺傳模式下與SZ顯著相關(guān)，其原因一方面可能與G/T基因型對轉(zhuǎn)錄起始因子與啟動子序列結(jié)合能力的改變有關(guān)，最終導(dǎo)致該基因表達量的改變；另一方面可能是G/T會啟動不同的蛋白質(zhì)或不同亞基的表達，導(dǎo)致谷氨酸脫羧酶GAD67亞型活性發(fā)生改變，促進SZ的發(fā)生。

Straub等[16]認為rs3749034位點突變會導(dǎo)致GAD1基因啟動子區(qū)產(chǎn)生鋅指結(jié)構(gòu)轉(zhuǎn)錄因子[ATP1A1 （Na+, K+ATPase α-1 subunit） regulatory element binding factor 6，AREB6］ 和成肌分化抗原（myogenic differentiation antigen，MYOD）的結(jié)合位點，這2種轉(zhuǎn)錄因子能夠影響GAD67mRNA的表達。但本次研究沒有發(fā)現(xiàn)該位點與SZ的相關(guān)性。這一結(jié)果與Addington等[17]發(fā)現(xiàn)rs3749034位點參與SZ早期發(fā)病的結(jié)果不一致。也有研究[18]發(fā)現(xiàn)兒茶酚-O-甲基轉(zhuǎn)移酶（catechol-O-methyltransferase，COMT） 基 因 Val158Met（rs4680）位點與rs3749034位點的相互作用會導(dǎo)致左側(cè)海馬旁回皮質(zhì)的厚度減少，且只有COMTVal158Met突變位點攜帶者才表現(xiàn)出rs3749034位點變異的陽性結(jié)果。這說明該位點可能參與SZ的發(fā)生，但并非主導(dǎo)因素。

Urak等[19]認為rs4906902位點的等位基因A突變?yōu)镚會改變神經(jīng)元特異性轉(zhuǎn)錄激活劑N-Oct-3的結(jié)合活性。rs4906902位點與rs8179184位點緊密連鎖，單倍型G-T頻率在SZ患者組的分布顯著高于對照組[20]，但本研究沒有發(fā)現(xiàn)rs4906902位點與SZ的相關(guān)性，這可能是因為SZ的異質(zhì)性所導(dǎo)致的。

本研究結(jié)果雖然表明rs3791878位點G/T基因型能夠顯著增加中國漢族男性人群SZ的發(fā)病風險，同時提示GAD1基因是中國漢族男性人群SZ的易感基因，但未涉及機制方面的研究工作。今后仍需大樣本研究進一步明確GAD1基因與SZ發(fā)生的相關(guān)性并開展相關(guān)的功能研究，以揭示其致病機制。