(姜鈺超　匡樂樂　王紅艷　段文元　蔡春泉

(1復(fù)旦大學(xué)生命科學(xué)學(xué)院遺傳學(xué)系　上海　200433; 2 天津市兒童醫(yī)院外科　天津　300074;3 濟南軍區(qū)總院心血管病研究所 濟南　250031)

?

SHMT1基因啟動子的rs638416位點與先天性心臟病的相關(guān)性

(姜鈺超1匡樂樂1王紅艷1段文元3蔡春泉2△

(1復(fù)旦大學(xué)生命科學(xué)學(xué)院遺傳學(xué)系上海200433;2天津市兒童醫(yī)院外科天津300074;3濟南軍區(qū)總院心血管病研究所 濟南250031)

目的通過研究絲氨酸羥甲基轉(zhuǎn)移酶(serine hydroxy methyl transferase,SHMT1)基因啟動子區(qū)的單核苷酸多態(tài)性(single nucleotide polymorphism,SNP)來探討該基因與先天性心臟病 (congenital heart disease,CHD)的相關(guān)性。方法選取來自山東省的201例CHD患兒(病例組)和192例正常兒童(對照組),采用Sanger測序法對rs638416和rs117940726兩個位點的基因分型進行病例-對照研究,進一步用雙熒光素酶報告基因分析其是否為功能性SNP。結(jié)果兩個SNP位點的基因型分布在病例組和對照組中均符合Hardy-Weinberg平衡;單位點分析顯示rs638416位點G等位基因在病例組中分布頻率顯著高于對照組 (P=0.009)。經(jīng)非條件Logistic回歸分析確定,Log-additive模型最適用于分析rs638416,且該模型下統(tǒng)計學(xué)差異最顯著 (OR=1.49,95%CI:1.11～ 2.01,P=0.007 4);雙熒光素酶報告基因分析顯示帶有G等位基因的啟動子序列轉(zhuǎn)錄效率較帶有C等位基因的啟動子序列低36% (P=0.013)。結(jié)論rs638416位點能影響SHMT1轉(zhuǎn)錄效率,并與心臟發(fā)育異常相關(guān),提示rs638416位點G等位基因是山東人群CHD發(fā)生的遺傳風(fēng)險因子。

先天性心臟病;SHMT1;胸腺嘧啶合成;葉酸

先天性心臟病 (congenital heart disease,CHD)在所有出生缺陷中發(fā)生率最高,世界范圍內(nèi)新生兒的發(fā)生率接近1%[1]。目前在我國CHD總體發(fā)病率逐年升高,而且呈現(xiàn)南方高于北方、城市高于農(nóng)村的特點,為提高國民出生人口素質(zhì),探索CHD致病機制已刻不容緩[2]。研究表明,環(huán)境因素和遺傳因素與CHD的發(fā)生有關(guān),該疾病的遺傳度為55%～65%[3]。

葉酸是人體必需的營養(yǎng)元素,其作為一碳單位供體參與核苷酸合成以及甲基轉(zhuǎn)移反應(yīng)。女性孕期補充葉酸可以預(yù)防包括CHD在內(nèi)的多種出生缺陷[4-6]。葉酸代謝途徑中關(guān)鍵基因多態(tài)性與CHD的關(guān)聯(lián)研究有助于揭示葉酸預(yù)防CHD發(fā)生的機制。

在葉酸通路中涉及葉酸中間代謝途徑的多個基因已經(jīng)被證明與CHD的發(fā)生有密切關(guān)系,但究竟是產(chǎn)生甲基化供體受阻,還是核苷酸從頭合成受限的效應(yīng)影響了心臟的早期發(fā)育,原因尚不明確[7-8]。本研究關(guān)注絲氨酸羥甲基轉(zhuǎn)移酶基因(serine hydroxy methyl transferase,SHMT1)。研究證實,細胞質(zhì)SHMT在脫氧胸苷酸的從頭合成過程中起關(guān)鍵作用,該酶能夠?qū)⒔z氨酸上的一碳單位傳遞至亞甲基四氫葉酸上,后者作為一碳單位供體在胸苷酸合酶 (thymidylate synthase,TYMS)催化下將dUMP轉(zhuǎn)化為dTMP[9-11]。亞甲基四氫葉酸除了向dUMP提供甲基外,同時也能提供甲硫氨酸合成所需的甲基,同位素示蹤研究顯示N5,N10-亞甲基四氫葉酸上的甲基主要用于合成dTMP[12]。SHMT1在參與葉酸代謝途徑的基因中相對重要,許多涉及葉酸代謝的疾病關(guān)聯(lián)研究將該基因納入到候選基因中。以往針對SHMT1的研究主要集中在各類腫瘤和神經(jīng)管畸形等,即使有針對CHD的研究因涉及葉酸代謝途徑而將SHMT1基因納入研究范圍,也只是將常見的rs1979277 (C1420T)、rs12952556等幾個位點作為標(biāo)簽SNP[13-15]。本研究選擇了SHMT1啟動子區(qū)的rs638416和rs117940726位點作為標(biāo)簽SNP,其中rs638416位點基因型為G/C,而rs117940726位點基因型為G/A。近年已有研究發(fā)現(xiàn)rs638416位點與某些疾病相關(guān),而尚未見rs117940726位點與疾病相關(guān)的報道。

本研究試圖通過探索dTMP從頭合成途徑的關(guān)鍵基因SHMT1與CHD的關(guān)聯(lián)性來闡述CHD的部分致病機制,利用山東人群樣本,以SHMT1基因啟動子區(qū)rs638416和rs117940726為標(biāo)簽SNP位點進行病例-對照關(guān)聯(lián)研究檢測這些多態(tài)位點是否具有CHD易感性,同時分析這些SNP位點是否影響SHMT1基因表達。

資 料 和 方 法

病例和對照樣本來源濟南軍區(qū)心血管?？漆t(yī)院在2009年6月至2012年6月期間收集來自山東地區(qū)的CHD病例201例和對照192例,并進行病理診斷。對照樣本和病例樣本年齡與性別盡量匹配,對照樣本無CHD家族史。201例CHD病例中包含129例間隔缺損類CHD。本研究符合倫理學(xué)各項規(guī)定,經(jīng)復(fù)旦大學(xué)倫理委員會批準,流行病調(diào)查資料收集及血樣采集均經(jīng)患兒監(jiān)護人的知情同意。

SHMT1基因啟動子區(qū)2個SNP位點分型

DNA提取 取外周血3 mL置于EDTA抗凝管中,冰箱保存;采用德國Qiagen公司血液基因組DNA提取試劑盒進行DNA的提取;抽提的DNA采用美國Thermo公司的Nanodrop 2000微量紫外分光光度計測定DNA濃度。

標(biāo)簽SNP選擇參考千人基因組中國北京人群數(shù)據(jù),挑選SHMT1起始密碼子上游啟動子區(qū)2.5 kbp范圍內(nèi)MAF>10%的2個SNP位點 (rs638416、rs117940726)納入本文研究范圍。

Sanger測序法分型針對兩個SNP位點設(shè)計PCR擴增引物,引物由上海生工生物技術(shù)有限公司合成。使用日本Takara公司的Taq酶、dNTP及反應(yīng)緩沖液,利用說明書提供的標(biāo)準體系和擴增條件進行PCR。使用美國ABI公司的Bigdye試劑盒及3730測序儀針對PCR擴增產(chǎn)物進行標(biāo)準化測序反應(yīng)并檢測。測序結(jié)果用Seqman軟件讀取并進行基因型統(tǒng)計。

雙熒光素酶報告基因檢測

報告基因載體構(gòu)建使用美國Promega公司的PGL3 Basic質(zhì)粒構(gòu)建報告基因載體,利用標(biāo)準的克隆步驟從基因組DNA中擴增出包含rs638416位點的DNA序列 (1.5 kbp)并構(gòu)建到報告基因載體上。已構(gòu)建成功的報告基因rs638416位點為G基因型,通過G>C的定點突變構(gòu)建C基因型報告基因。

報告基因熒光信號檢測報告基因檢測基于293T細胞,按照標(biāo)準的細胞培養(yǎng)流程及美國Promega公司雙熒光素酶報告基因系統(tǒng)說明書進行實驗操作。以PGL3報告基因載體/pRL-TK內(nèi)參載體按10∶1的比例進行細胞共轉(zhuǎn)染。利用美國Promega公司GloMax 96 微孔板發(fā)光檢測儀檢測G/A基因型的PGL3載體報告基因表達。

統(tǒng)計學(xué)分析通過χ2檢驗計算正常對照的Hardy-Weinberg平衡 (HWE)的理論值。通過Haploview軟件對兩個SNP位點進行連鎖不平衡分析。相對OR和95%CI在校正年齡和性別之后由非條件Logistic回歸模型計算獲得。利用在線工具SNPstats的Codominant、Dominant、Recessive、Overdominant以及Log-additive 5種模型對SNP進行分析。采用SPSS 19.0軟件進行統(tǒng)計分析,P<0.05為差異有統(tǒng)計學(xué)意義。

結(jié)　　　果

SHMT1基因啟動子區(qū)SNP連鎖不平衡分析軟件計算結(jié)果顯示,rs638416和rs117940726之間SHMT1基因啟動子區(qū)的2個SNP位點不存在強連鎖現(xiàn)象(圖1)。

圖1SHMT1基因啟動子區(qū)SNP位點的連鎖不平衡分析

Fig 1Linkage disequilibrium analysis of SNPs in the promoter region ofSHMT1

SHMT1基因啟動子區(qū)單個位點SNP與CHD的關(guān)聯(lián)分析SHMT1基因啟動子區(qū)2個SNP位點rs638416和rs117940726的各基因型在對照組和病例組中的分布情況如表1所示。HWE檢驗結(jié)果顯示,在對照組中rs638416和rs117940726的χ2檢驗P值分別為0.46和0.16,基因型分布符合遺傳平衡。在經(jīng)過年齡和性別因素校正后,根據(jù)遺傳學(xué)模型進行Logistic回歸分析,AIC和BIC值越小說明模型越精確。結(jié)果顯示:除了Overdominant模型外,其余4種模型中rs638416位點都與CHD顯著關(guān)聯(lián),且認為Log-additive模型最適用于分析該SNP位點 (表2);而在5組模型中均未觀察到rs117940726位點與CHD的關(guān)聯(lián)性 (表3)。

雙熒光素酶報告基因檢測 如圖2所示,rs638416位點的G等位基因與C等位基因的報告基因相比熒光強度比值下降了36% (P=0.013)。

討　　　論

為探索孕前補充葉酸干預(yù)CHD發(fā)生的機制,本研究選取葉酸代謝網(wǎng)絡(luò)中與dTMP合成相關(guān)的關(guān)鍵酶SHMT1作為研究對象,在山東人群中通過選擇SHMT1基因上游啟動子區(qū)2個SNP進行病例-對照研究。結(jié)果顯示rs638416位點與CHD有關(guān),且通過雙熒光素酶報告基因分析確定rs638416位點的G等位基因能夠顯著降低啟動子的轉(zhuǎn)錄效率。本研究首次在人群中發(fā)現(xiàn)SHMT1與CHD的關(guān)聯(lián)性,研究結(jié)果為“dTMP合成受阻是CHD病因之一”的觀點提供了有力證據(jù)。

表1　在病例組和對照組中rs638416和rs117940726基因型的頻率分布

表2　rs638416位點基因型分布的統(tǒng)計分析

AIC:Akaike′s information criterion;BIC:Bayes information criterion.

表3　rs117940726位點基因型分布的統(tǒng)計分析

AIC:Akaike′s information criterion;BIC:Bayes information criterion.

SHMT1是一個dTMP合成通路上游的限速酶,一旦SHMT1被破壞,dTMP合成效率大為降低。在SHMT1純合敲除小鼠模型構(gòu)建的研究中發(fā)現(xiàn),雖有部分小鼠仍可無表型存活,但小鼠胚胎有較高的概率出現(xiàn)嚴重的神經(jīng)管畸形表型,同時觀察到由于dTMP的相對缺乏導(dǎo)致dUMP混入基因組序列中,此外SHMT1純合敲除小鼠胚胎的體質(zhì)量顯著小于雜合敲除及野生型的胚胎[16]。小鼠模型顯示了dTMP合成途徑對于細胞乃至生物個體生長具有非常重要的作用,尤其是在胚胎和器官發(fā)育的早期階段,然而在SHMT1敲除的小鼠中并未觀察到心臟發(fā)育異常明顯的表型,這也提示了人體心臟發(fā)育過程中可能對dTMP劑量需求更加敏感。

圖2rs638416位點C等位基因與G等位基因的熒光素酶報告基因檢測結(jié)果

Fig 2The luciferase reporter gene assays results of C allele and G allele of rs638416

本研究利用293T細胞驗證rs638416位點對啟動子效率的影響,而此前有研究在HepG2細胞中對rs638416位點開展雙熒光素酶報告基因分析,兩項研究得出的結(jié)論一致,G位點均可以降低啟動子的轉(zhuǎn)錄效率[17]。我們推測rs638416位點所在區(qū)域可能是一個轉(zhuǎn)錄因子的結(jié)合位點,rs638416位點的G等位基因可能會破壞這個轉(zhuǎn)錄因子結(jié)合位點與轉(zhuǎn)錄因子的結(jié)合能力,從而降低啟動子的轉(zhuǎn)錄效率。此前查證的文獻中rs638416位點只與顳下頜關(guān)節(jié)紊亂、卵巢癌及肺癌的發(fā)生明確相關(guān)[18-20]。本研究解釋了rs638416位點的部分機制,作為功能性SNP位點,其在今后有關(guān)葉酸代謝通路的研究中應(yīng)更受關(guān)注。

葉酸在人體內(nèi)的主要作用除了合成嘌呤及胸腺嘧啶外,同時也生成甲基化的S-腺苷甲硫氨酸。因此補充葉酸的作用主要是促進嘌呤、胸腺嘧啶及S-腺苷甲硫氨酸的合成。本課題組此前進行的一項針對合成胸腺嘧啶的關(guān)鍵基因TYMS多態(tài)性與CHD的關(guān)聯(lián)研究未發(fā)現(xiàn)TYMS與CHD有關(guān)聯(lián)性[10]。我們認為前期研究結(jié)果與本研究結(jié)論并不矛盾:TYMS是從dUMP合成dTMP的唯一酶,因而在進化角度上必須非常保守,影響酶活性的變異均可能被淘汰,該證據(jù)顯示了胸腺嘧啶合成通路的重要性。生物體演化出SHMT1及SHMT2這兩個同源酶來共同催化產(chǎn)生N5,N10-亞甲基四氫葉酸,此外四氫葉酸也可以被催化為N5,N10-亞甲基四氫葉酸,此種多重補償機制同樣體現(xiàn)了dTMP合成通路的重要性。我們推測在心臟發(fā)育早期,細胞增殖對N5,N10-亞甲基四氫葉酸劑量非常敏感,一旦N5,N10-亞甲基四氫葉酸劑量不足引發(fā)胸腺嘧啶合成受阻,很可能誘發(fā)CHD。

本研究結(jié)果在遺傳學(xué)的層面上揭示了葉酸對預(yù)防CHD的部分機制,而針對dTMP合成通路與CHD病的關(guān)系需要進一步在不同人群中進行驗證并開展相關(guān)功能研究,以完善對dTMP合成通路參與CHD發(fā)生機制的認識。

[1]VAN DER LINDE D,KONINGS EE,SLAGER MA,etal.Birth prevalence of congenital heart disease worldwide:a systematic review and meta-analysis[J].JAmCollCardiol,2011,58 (21):2241-2247.

[2]朱軍,周光萱,代禮,等.1996～2000年全國圍產(chǎn)期先天性心臟病發(fā)生率的分析[J].四川大學(xué)學(xué)報(醫(yī)學(xué)版),2005,35 (6):875-877.

[3]VAN DER BOM T,ZOMER AC,ZWINDERMAN AH,etal.The changing epidemiology of congenital heart disease[J].NatRevCardiol,2011,8 (1):50-60.

[4]BAILEY LB,BERRY RJ.Folic acid supplementation and the occurrence of congenital heart defects,orofacial clefts,multiple births,and miscarriage[J].AmJClinNutr,2005,81 (5):1213S-1217S.

[5]HERNANDEZ-DIAZ S,WERLER MM,WALKER AM,etal.Folic acid antagonists during pregnancy and the risk of birth defects[J].NEnglJMed,2000,343 (22):1608-1614.

[6]IOoNESCU-ITTU R,MARELLI AJ,MACKIE AS,etal.Prevalence of severe congenital heart disease after folic acid fortification of grain products:time trend analysis in Quebec,Canada[J].BMJ,2009,338:b1673.

[7]ZHAO J,YANG XY,GONG XH,etal.Functional variant in methionine synthase reductase intron-1 significantly increases the risk of congenital heart disease in the Han Chinese population[J].Circulation,2012,125 (3):482-490.

[8]ZHAO JY,QIAO B,DUAN WY,etal.Genetic variants reducing MTR gene expression increase the risk of congenital heart disease in Han Chinese populations[J].EurHeartJ,2014,35 (11):733-742.

[9]LEDUC D,GRAZIANI S,MESLET-CLADIERE L,etal.Two distinct pathways for thymidylate (dTMP) synthesis in (hyper) thermophilic Bacteria and Archaea[J].BiochemSocTrans,2004,32 (Pt 2):231-235.

[10]ZHAO JY,SUN JW,GU ZY,etal.Genetic polymorphisms of the TYMS gene are not associated with congenital cardiac septal defects in a Han Chinese population[J].PLoSOne,2012,7 (2):e31644.

[11]ERMENS AA,SCHOESTER M,LINDEMANS J,etal.Effect of nitrous oxide on folate coenzyme distribution and de novo synthesis of thymidylate in human bone marrow cells[J].ToxicolInVitro,1992,6 (2):133-137.

[12]HERBIG K,CHIANG EP,LEE LR,etal.Cytoplasmic serine hydroxy methyl transferase mediates competition between folate-dependent deoxyribonucleotide and S-adenosylmethionine biosyntheses[J].JBiolChem,2002,277 (41):38381-3839.

[13]PANDEY S,GARG P,LEE S,etal.Nucleotide biosynthesis arrest by silencing SHMT1 function via vitamin B6-coupled vector and effects on tumor growth inhibition[J].Biomaterials,2014,35 (34):9332-9342.

[14]LI L,JIN C,ZHENG Q,etal.Polymorphism of cytosolic serine hydroxy methyl transferase and breast cancer risk:evidence from a meta-analysis[J].TumourBiol,2014,35 (8):7361-7367.

[15]BEAUDIN AE,ABARINOV EV,MALYSHEVA O,etal.Dietary folate,but not choline,modifies neural tube defect risk in Shmt1 knockout mice[J].AmJClinNutr,2012,95 (1):109-114.

[16]BEAUDIN AE,ABARINOV EV,NODEN DM,etal.Shmt1 and de novo thymidylate biosynthesis underlie folate-responsive neural tube defects in mice[J].AmJClinNutrition,2011,93 (4):789-798.

[17]HEBBRING SJ,CHAI Y,JI Y,etal.Serine hydroxy methyl transferase 1 and 2:gene sequence variation and functional genomic characterization[J].JNeurochem,2012,120 (6):881-890.

[18]WANG L,LU J,AN J,etal.Polymorphisms of cytosolic serine hydroxy methyl transferase and risk of lung cancer:a case-control analysis[J].LungCancer,2007,57 (2):143-151.

[19]KELEMEN LE,SELLERS TA,SCHILDKRAUT JM,etal.Genetic variation in the one-carbon transfer pathway and ovarian cancer risk[J].CancerRes,2008,68 (7):2498-2506.

[20]ANEIROS-GUERRERO A,LENDINEZ AM,PALOMARES AR,etal.Genetic polymorphisms in folate pathway enzymes,DRD4 and GSTM1 are related to temporomandibular disorder[J].BMCMedGenet,2011,12:75.

E-mail:tjpns@126.com

Association of rs638416 in SHMT1 promoter and congenital heart diseases

JIANG Yu-chao1, KUANG Le-le1, WANG Hong-yan1, DUAN Wen-yuan3, CAI Chun-quan2△

(1DepartmentofGenetics,SchoolofLifeSciences,FudanUniversity,Shanghai200433,China;2DepartmentofSurgery,Children′sHospital,Tianjing300074,China;3InstituteofCardiovascularDissease,GeneralHospitalofJinanMilitaryRegion,Jinan250031,China)

ObjectiveTo explore the association between single nucleotide polymorphism (SNPs) in the promoter region of cytoplasmic serine hydroxy methyl transferase (SHMT1) and congenital heart disease (CHD).MethodsWe performed a case-control study including 201 CHD patients and 192 control participants from Shandong province,China through Sanger sequencing of rs638416 and rs117940726.Then we conducted luciferase reporter gene assays to confirm if they were functional SNP.ResultsThese two SNPs were genotyped and they were in Hardy-Weinberg equilibrium both in control and case groups.The single locus analysis revealed the genotype distribution of rs638416 was significant different between case patients and control subjects (P=0.009).In the unconditional Logistic regression analysis,Log-additive model was the most suitable model to analysis rs638416,and in this model association data showed the most significant statistical different (OR=1.49,95%CI:1.11-2.01,P=0.007 4).The result of luciferase reporter gene assays also showed 36% transcription activity reduction (P=0.0013) of the promoter region containing G allele when compared to the promoter region containing G allele.ConclusionsWe confirmed rs638416 reduced the transcription activity ofSHMT1 and may involve with abnormal heart development.The G allele of rs638416 is a genetic risk factor of CHD in Shandong population.

congenital heart diseases;SHMT1;purine de novo synthesis;folic acid

Q39, R725.4

Adoi: 10.3969/j.issn.1672-8467.2016.04.008

2015-07-28;編輯:段佳)

天津市衛(wèi)生行業(yè)重點攻關(guān)項目 (12KG116);天津市自然科學(xué)基金 (14JCYBJC25000)

*This work was supported by the Key Project of Tianjin Health Care Professionals (12KG116) and the Natural Science Foundation of Tianjin City,China (14JCYBJC25000).