山羊STAT3基因克隆、生物信息學(xué)分析及甲基化修飾研究

賈文超，劉亮亮，吳賢鋒，趙釗艷，王　珂，王　毛，高　靜，

王立強(qiáng)1，陳　宏1，潘傳英1，藍(lán)賢勇1*

(1.西北農(nóng)林科技大學(xué) 動(dòng)物科技學(xué)院/生命科學(xué)學(xué)院，陜西省農(nóng)業(yè)分子生物學(xué)重點(diǎn)實(shí)驗(yàn)室，楊凌 712100；2.江蘇科技大學(xué) 計(jì)算機(jī)科學(xué)與工程學(xué)院，鎮(zhèn)江 212003)

?

山羊STAT3基因克隆、生物信息學(xué)分析及甲基化修飾研究

賈文超1#，劉亮亮2#，吳賢鋒1，趙釗艷1，王珂1，王毛1，高靜1，

王立強(qiáng)1，陳宏1，潘傳英1，藍(lán)賢勇1*

(1.西北農(nóng)林科技大學(xué) 動(dòng)物科技學(xué)院/生命科學(xué)學(xué)院，陜西省農(nóng)業(yè)分子生物學(xué)重點(diǎn)實(shí)驗(yàn)室，楊凌 712100；2.江蘇科技大學(xué) 計(jì)算機(jī)科學(xué)與工程學(xué)院，鎮(zhèn)江 212003)

摘要：本研究利用RT-PCR、TA克隆和重亞硫酸鹽測(cè)序(BSP)等方法，旨在克隆、分析山羊STAT3基因完整編碼區(qū)，并解析該基因甲基化修飾水平及其與體重的關(guān)系。結(jié)果表明，山羊STAT3基因編碼區(qū)全長(zhǎng)2 313 bp，編碼770個(gè)氨基酸；與綿羊、牛、野豬、小鼠、大鼠和人的氨基酸序列一致性分別為99.7%、99.6%、99.4%、99.2%、99.4%和99.5%。山羊高體重組與低體重組之間甲基化差異研究發(fā)現(xiàn)，高體重組STAT3基因啟動(dòng)子區(qū)甲基化水平顯著高于低體重組(P=0.020)，提示該基因甲基化修飾對(duì)體重具有顯著影響，可作為標(biāo)記輔助選擇(MAS)的有效表觀遺傳標(biāo)記。這些結(jié)果為研究山羊肉用性能的遺傳與表觀遺傳機(jī)制提供了科學(xué)資料。

關(guān)鍵詞：山羊；STAT3基因；分子克隆；序列分析；DNA甲基化

羊肉較豬肉更細(xì)嫩，而且脂肪與膽固醇含量更少，為此，近年來(lái)，羊肉越來(lái)越受到消費(fèi)者的青睞。當(dāng)前，進(jìn)一步提升山羊產(chǎn)肉性能已成為現(xiàn)代養(yǎng)羊業(yè)發(fā)展的關(guān)鍵問題，故對(duì)山羊產(chǎn)肉調(diào)控機(jī)制進(jìn)行研究顯得極其重要。目前，關(guān)于山羊產(chǎn)肉調(diào)控的表觀遺傳機(jī)制的研究很少，其中山羊重要基因表觀遺傳修飾的研究則更是鮮見報(bào)道[1]。

信號(hào)轉(zhuǎn)導(dǎo)與轉(zhuǎn)錄激活因子(STATs)是一類活化后能轉(zhuǎn)入細(xì)胞核、并與相應(yīng)DNA結(jié)合的蛋白家族，具有信號(hào)轉(zhuǎn)導(dǎo)和轉(zhuǎn)錄調(diào)控的雙重功能[2]。STATs家族包括7個(gè)成員，即STAT1-STAT4、STAT5a、STAT5b和STAT6，均含有SH2同源結(jié)構(gòu)域和C端酪氨酸磷酸化位點(diǎn)[2]。作為STATs家族成員之一，STAT3 于1994年作為白細(xì)胞介素-6(IL-6)信號(hào)傳遞中的急性期反應(yīng)因子(APRF) 被純化，它在調(diào)節(jié)肝急性期蛋白和抑制單核細(xì)胞生長(zhǎng)等方面具有重要作用[3-5]，不僅能調(diào)節(jié)細(xì)胞自噬[6]，還能影響胚胎早期發(fā)育[7]。動(dòng)物STAT3通過調(diào)控催乳素(PRL)和生長(zhǎng)激素(GH)的表達(dá)來(lái)影響體重等生產(chǎn)性能[8-9]。據(jù)此推測(cè)，STAT3基因?qū)w重等生產(chǎn)性狀有重要遺傳效應(yīng)，故該基因是研究現(xiàn)代山羊產(chǎn)業(yè)中影響產(chǎn)肉性能的重要候選基因。

所謂表觀遺傳就是指在DNA 序列不改變的情況下基因表達(dá)與調(diào)控的可遺傳變化，主要包括DNA甲基化、組蛋白修飾、染色質(zhì)重塑和非編碼RNA等。作為最主要的表觀遺傳修飾方式之一，DNA甲基化具有時(shí)空特異性的特點(diǎn)，常見于基因CpG島的CG二核苷酸位點(diǎn)[1，10]。盡管DNA甲基化不涉及DNA序列的改變，但它卻能影響基因的表達(dá)調(diào)控，它可以通過有絲分裂和減數(shù)分裂穩(wěn)定地將DNA甲基化信息傳遞給后代，是一種新的表觀遺傳標(biāo)記[10]。目前，關(guān)于山羊STAT3基因甲基化修飾及其與產(chǎn)肉性能關(guān)系的研究尚未見報(bào)道。

因山羊STAT3基因編碼區(qū)未被克隆，其氨基酸信息未知，故不利于其功能的深入解析；同時(shí)，因山羊STAT3基因甲基化模式未知，故其甲基化功能研究也無(wú)法深入。為此，本研究旨在克隆STAT3基因編碼區(qū)(CDS)，并分析其序列特征，同時(shí)還解析該基因啟動(dòng)子區(qū)甲基化修飾對(duì)體重的影響，期望為山羊產(chǎn)肉調(diào)控的遺傳與表觀遺傳機(jī)制研究提供科學(xué)資料。

1材料與方法

1.1試驗(yàn)材料

山羊樣品包括兩類：一類是用于基因克隆的組織樣品，另一類是用于DNA甲基化分析的耳組織樣品。其中，組織樣品來(lái)源于6只健康的2～3歲左右布爾山羊和海南黑山羊肝組織，用于山羊STAT3基因的克隆研究；耳組織樣品來(lái)源于284只成年海南黑山羊，用于DNA甲基化修飾研究。在DNA甲基化研究中，基于數(shù)理統(tǒng)計(jì)的概率知識(shí)，參考C.Y.Pan等[11]和趙海諭[12]的方法，對(duì)收集的海南黑山羊群體的體重進(jìn)行分類，將高于群體均值1.96倍標(biāo)準(zhǔn)差(Mean+1.96SD)的個(gè)體設(shè)為高體重組(H)，而低于群體均值1.96倍標(biāo)準(zhǔn)差(Mean-1.96SD)的個(gè)體設(shè)為低體重組(L)。

1.2主要試劑

Trizol 購(gòu)自 Invitrogen 公司(USA)；RevertAid First Strand cDNA Synthesis Kit、DNA重亞硫酸氫鹽處理試劑盒EZ?DNA Methylation-GoldTM Kit、ZYMO RESEARCH試劑盒和DNA柱式膠回收試劑盒購(gòu)于康寧生命科學(xué)有限公司；T-easy載體購(gòu)自Promega 公司(USA)；iQTMSYBR?Green Supermix 購(gòu)自 Bio-Rad公司(USA)；2×TaqPCR Master Mix 購(gòu)自天根生化科技有限公司(中國(guó)北京)；大腸桿菌(Escherichiacoli)DH5α感受態(tài)細(xì)胞(Trans-5α)為本實(shí)驗(yàn)室保存。

1.3方法

1.3.1山羊STAT3基因的分子克隆

1.3.1.1總RNA的提取和cDNA第一鏈的合成：采用Trizol 法提取肝組織的總RNA[12-13]；利用紫外分光光度計(jì)和瓊脂糖凝膠電泳法檢測(cè)總RNA質(zhì)量；按照反轉(zhuǎn)錄試劑盒說(shuō)明書獲得cDNA：5×PrimeScript Buffer 2(for Real Time)2.0 μL，PrimeScript RT Enzyme Mix I 0.5 μL，RT Primer Mix 0.5 μL，除基因組DNA的反應(yīng)液 1.0 μL，RNase Free dH2O 10.0 μL。反應(yīng)條件：37 ℃反轉(zhuǎn)錄15 min→85 ℃反轉(zhuǎn)錄酶失活5 s →反轉(zhuǎn)錄獲得cDNA，-20 ℃?zhèn)溆?；根?jù)反轉(zhuǎn)錄cDNA的質(zhì)量和濃度構(gòu)建cDNA池，并以此cDNA池為模板克隆STAT3基因。1.3.1.2克隆引物的設(shè)計(jì)與PCR擴(kuò)增：參考NCBI的基因序列(NM_001012671)，利用Primer Premier軟件(5.0)設(shè)計(jì)克隆山羊STAT3基因CDS區(qū)的PCR引物1對(duì)(上游：5′-AGCAGCAGTTGTGACCCCT-GATTC-3′；下游:5′-TGTAAGGGCTGGTGGGACATGGGT-3′)，其預(yù)期擴(kuò)增片段長(zhǎng)度為2 410 bp。在 PCR過程中，總反應(yīng)體系20 μL：上、下游引物各0.4 μL，RNase-Free ddH2O 7.4 μL，2×PowerTaqPCR MasterMix 11 μL，cDNA模板0.8 μL；PCR反應(yīng)條件：94 ℃預(yù)變性5 min；94 ℃變性30 s，60 ℃退火30 s，72 ℃延伸150 s，36個(gè)循環(huán)；72 ℃終延伸10 min，10 ℃保存。

1.3.1.3PCR產(chǎn)物的回收與單克隆的測(cè)序：利用0.8%瓊脂糖凝膠電泳檢測(cè)PCR產(chǎn)物，再根據(jù)凝膠回收試劑盒回收目的片段；質(zhì)檢合格純化后PCR產(chǎn)物與 T-easy載體4 ℃過夜連接；連接產(chǎn)物轉(zhuǎn)化至25 μL Trans-5α感受態(tài)細(xì)胞，經(jīng)復(fù)蘇、平板培養(yǎng)后，通過藍(lán)、白斑篩選挑取白色單克隆菌落，用LB(Amp+)培養(yǎng)液37 ℃擴(kuò)大培養(yǎng)，之后利用菌液PCR 鑒定陽(yáng)性克隆，并送南京金斯瑞公司測(cè)序。

1.3.2生物信息學(xué)分析利用DNAMAN 軟件(6.0)對(duì)測(cè)序所獲得的山羊STAT3序列進(jìn)行開放閱讀框(ORF)和氨基酸(aa)序列預(yù)測(cè)；采用 MEGA軟件(5.0)構(gòu)建系統(tǒng)進(jìn)化樹；利用 EXPASY(http://www.expasy.org)對(duì)山羊STAT3蛋白進(jìn)行理化性質(zhì)研究；跨膜結(jié)構(gòu)預(yù)測(cè)采用 TMHMM(http://www.cbs.dtu.dk/services/TMHMM)進(jìn)行；利用NetPhos 2.0(http://cbs.dtu.dk/services/NetPhos/)對(duì)STAT3 進(jìn)行磷酸化位點(diǎn)預(yù)測(cè)，利用NetNGlyc(http://www.cbs.dtu.dk/services/NetNGlyc/)分析蛋白質(zhì)糖基化位點(diǎn)；利用在 線 程序Hopfield(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopm.html) 對(duì)蛋白質(zhì)二級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)；亞細(xì)胞定位采用 PsortⅡ(http://psort.hgc.jp/)進(jìn)行[14-21]。

1.3.3山羊STAT3基因甲基化修飾研究

1.3.3.1基因組DNA的提取與重亞硫酸鹽處理：采用標(biāo)準(zhǔn)的苯酚-氯仿方法提取來(lái)源于高體重組和低體重組個(gè)體的耳組織基因組DNA[11-12，18-19]；嚴(yán)格按照EZ DNA Methylation-GoldTM Kit試劑盒說(shuō)明書對(duì)基因組DNA進(jìn)行重亞硫酸鹽處理。

1.3.3.2山羊STAT3基因甲基化引物的設(shè)計(jì)與擴(kuò)增：利用在線軟件Methprimer(http://www.urogene.org/methprimer/)篩選山羊STAT3基因CpG島，并針對(duì)CpG島設(shè)計(jì)了1對(duì)位于啟動(dòng)子區(qū)的甲基化引物(上游Methprimer-STAT3-M1-F：5′-GTATTGTTTTTTTATTTTTTTGAAT-3′；下游STAT3-M1-R:5′-CCTACTTTAAACTTCAATT-TCTAC-3′)，擴(kuò)增片段為191 bp(-807～-617 bp)，包括10個(gè)CpG位點(diǎn)。25 μL PCR反應(yīng)體系：50 ng·μL-1的重亞硫酸鹽處理后的基因組DNA 1.0 μL，10 μmol·L-1的上、下游引物各0.5 μL，2×Reaction Mix(含Buffer、TaqDNA聚合酶、Mg2+、dNTPs等) 12.5 μL，滅菌超純水10.5 μL。Touch-Down PCR反應(yīng)程序：94 ℃預(yù)變性5 min；94 ℃變性30 s，68～54 ℃(每?jī)裳h(huán)降2 ℃)30 s，72 ℃延伸30 s，14個(gè)循環(huán)；94 ℃變性30 s，51 ℃復(fù)性30 s 72 ℃延伸30 s，40個(gè)循環(huán)；72 ℃延伸10 min；4 ℃保存。

1.3.3.3重亞硫酸鹽處理后測(cè)序(BSP)和DNA甲基化模式分析：參考1.3.1.3中描述的方法進(jìn)行TA克隆；隨機(jī)挑選10～15個(gè)陽(yáng)性克隆送南京金斯瑞公司測(cè)序，并利用T7正向或M13反向測(cè)序；根據(jù)DNA甲基化測(cè)序結(jié)果，分析陽(yáng)性克隆中甲基化的CpG二核苷酸位點(diǎn)的百分?jǐn)?shù)[11-13，15，17-19]；分析高體重組與低體重組STAT3基因啟動(dòng)子區(qū)CpG島的整體甲基化率[17-19，22-24]；分析高體重組與低體重組STAT3基因啟動(dòng)子區(qū)甲基化差異與體重的關(guān)系：根據(jù)甲基化位點(diǎn)的甲基化差異數(shù)據(jù)特征，利用SPSS(18.0)軟件中Fisher精確概率檢驗(yàn)法分析高體重組與低體重組的DNA甲基化差異[11-12，16]。

2結(jié)果

2.1山羊STAT3基因CDS的分子克隆及氨基酸序列分析

利用RT-PCR方法獲得一條長(zhǎng)度為2 410 bp左右的條帶(圖1)，與預(yù)期結(jié)果一致。經(jīng)TA克隆、測(cè)序，成功地獲得了山羊STAT3基因CDS序列；ORF finder分析表明，山羊STAT3基因CDS全長(zhǎng)2 313 bp，編碼770個(gè)氨基酸(圖2)。

圖1　山羊STAT3基因的RT-PCR電泳結(jié)果Fig.1　The RT-PCR electrophoresis results of the goat STAT3 gene

2.2山羊STAT3氨基酸序列的系統(tǒng)進(jìn)化樹構(gòu)建

將山羊STAT3氨基酸序列與牛(NP_001012689)、人(AAS66986)、小鼠(AAA19452)、大鼠(AAH87025)、野豬(ABF21150)、綿羊(XP_004012974)氨基酸序列進(jìn)行同源性比對(duì)，發(fā)現(xiàn)它們高度同源；在氨基酸同源性分析中(圖3)，山羊STAT3基因氨基酸序列與牛、綿羊、小鼠、大鼠、野豬和人氨基酸序列的相似性分別為99.6%、99.7%、99.2%、99.4%、99.4%和99.5%。系統(tǒng)發(fā)育樹顯示(圖4)，山羊與牛的親緣關(guān)系最為接近，與人類、大鼠、小鼠和綿羊親緣關(guān)系較遠(yuǎn)。

圖2　山羊STAT3基因CDS與氨基酸序列Fig.2　The CDS and amino acid sequence of the goat STAT3 gene

2.3山羊STAT3蛋白結(jié)構(gòu)與功能預(yù)測(cè)

2.3.1山羊STAT3蛋白理化性質(zhì)分析山羊STAT3蛋白由770個(gè)氨基酸殘基組成，其分子式為：C3909H6161N1063O1171S41，其分子量為88 099.8D，其理論等電點(diǎn)(pI)為5.94；該蛋白的氨基酸組成及其比例分別為：Ala(5.8%)、Arg(4.7%)、Asn(4.9%)、Asp(4.2%)、Cys(1.8%)、Gln(7.8%)、Glu(7.7%)、Gly(4.9%)、His(1.7%)、Ile(5.3%)、Leu(10.6%)、Lys(6.0%)、Met(3.5%)、Phe(3.9%)、Pro(4.3%)、Ser(6.8%)、Thr(6.1%)、Trp(1.9%)、Tyr(2.7%) 和Val(5.3%)。其不穩(wěn)定系數(shù)為49.58，且預(yù)計(jì)在哺乳動(dòng)物網(wǎng)織紅細(xì)胞內(nèi)半衰期為30h；脂溶指數(shù)為83.58，總平均疏水指數(shù)-0.400，屬親水蛋白。

圖4　不同物種STAT3基因的系統(tǒng)發(fā)育樹Fig.4　Phylogenetic tree of STAT3 gene among various species

2.3.2山羊STAT3蛋白結(jié)構(gòu)預(yù)測(cè)山羊STAT3蛋白主要由α-螺旋(46.62%)、β折疊(18.18%)、β轉(zhuǎn)角(7.40%)及無(wú)規(guī)則卷曲(27.79%)構(gòu)成。采用SWISS-MODEL和PyMol軟件對(duì)山羊STAT3建模獲得二級(jí)結(jié)構(gòu)模型及結(jié)構(gòu)域的預(yù)測(cè)(圖5，圖6)。

圖5　山羊STAT3蛋白質(zhì)二級(jí)結(jié)構(gòu)分析Fig.5　The secondary structure analysis of goat STAT3 protein

圖6　山羊STAT3蛋白質(zhì)結(jié)構(gòu)域Fig.6　Protein domains of goat STAT3

2.3.3山羊STAT3亞細(xì)胞定位、蛋白磷酸化和糖基化位點(diǎn)預(yù)測(cè)山羊STAT3蛋白在細(xì)胞質(zhì)、細(xì)胞核、線粒體基質(zhì)、內(nèi)質(zhì)網(wǎng)、膜泡分泌系統(tǒng)及質(zhì)膜中的分布比例分別為34.8%、30.4%、17.4%、8.7%、4.3%和4.3%，在細(xì)胞質(zhì)中分布最多。磷酸化分析顯示(圖7a)，山羊STAT3蛋白存在較多的絲氨酸(Ser)磷酸化位點(diǎn)，而酪氨酸(Tyr)和蘇氨酸(Thr)磷酸化位點(diǎn)較少。糖基化位點(diǎn)分析表明，山羊STAT3蛋白存在3個(gè)糖基化修飾位點(diǎn)，分別位于氨基酸序列的第192、401和538位(圖7b)。

2.4山羊STAT3基因甲基化修飾分析

山羊STAT3基因啟動(dòng)子區(qū)CpG島的PCR片段長(zhǎng)度為191bp(圖8)，共包括了10個(gè)CpG二核苷酸位點(diǎn)(圖9，圖10)。經(jīng)研究發(fā)現(xiàn)，高體重組(H)和低體重組(L)的整體甲基化率分別為42.86%和0。利用獨(dú)立卡方χ2進(jìn)行檢驗(yàn)，發(fā)現(xiàn)H組和L組之間DNA甲基化修飾程度存在顯著差異(P=0.020)(圖11)。

3討論

圖7　山羊STAT3蛋白的磷酸化位點(diǎn)(a)和糖基化位點(diǎn)(b)Fig.7　Phosphorylation site(a) and glycosylation site(b) analyses of goat STAT3 protein

圖8　山羊STAT3-M1 PCR電泳結(jié)果Fig.8　PCR electrophoresis results of the goat STAT3-M1

STAT3作為一種重要轉(zhuǎn)錄因子，是胞漿中信息的傳遞者，是胞核中轉(zhuǎn)錄因子的激活者，還是乳腺癌細(xì)胞中多種細(xì)胞因子信號(hào)傳導(dǎo)途徑的匯聚點(diǎn)，參與腫瘤發(fā)生發(fā)展的多個(gè)環(huán)節(jié)[20]。研究表明，STATs家族基因?qū)w重具有顯著遺傳效應(yīng)，其中STAT5a、STAT3和STAT1等蛋白參與GH與PRL的信號(hào)轉(zhuǎn)導(dǎo)，從而對(duì)山羊體重均有重要影響；H.W.Davey等[21]發(fā)現(xiàn)，STAT5a與脂肪細(xì)胞增生密切相關(guān)；過表達(dá)動(dòng)物IL-6將導(dǎo)致脂肪分解基因的mRNA水平升高，從而觸發(fā)作為IL-6傳遞中的STAT3磷酸化水平，從而導(dǎo)致體重下降。本課題組最新研究表明，山羊STAT5a和STAT3基因遺傳變異對(duì)生長(zhǎng)性狀有顯著影響[22-23]；STAT5a和STAT3為同一基因家族成員，且同源性較高、功能相似。因此，研究STAT3基因的功能對(duì)理解山羊體重和產(chǎn)肉性能具有重要意義。

本研究成功克隆了山羊STAT3基因CDS全長(zhǎng)(2 313bp，編碼770個(gè)氨基酸)，為該基因的功能研究奠定了基礎(chǔ)。不同各物種STAT3氨基酸序列比對(duì)結(jié)果提示該基因較為保守。生物信息學(xué)預(yù)測(cè)表明，第600～700位氨基酸具有SH2結(jié)構(gòu)域，而且該區(qū)域中還存在較多的絲氨酸、蘇氨酸(Thr)和酪氨酸(Tyr)的磷酸化位點(diǎn)，這些為該基因深入的功能驗(yàn)證提供了科學(xué)參考。

研究表明，基因啟動(dòng)子區(qū)DNA甲基化修飾會(huì)抑制基因的表達(dá)[24]，當(dāng)基因啟動(dòng)子區(qū)被甲基化后，甲基化結(jié)合蛋白會(huì)抑制SP1轉(zhuǎn)錄因子介導(dǎo)的轉(zhuǎn)錄激活[25]。STAT3基因在JAK-STAT信號(hào)途徑中發(fā)揮著重要功能[26]。E.Kim等[27]指出，組蛋白甲基化轉(zhuǎn)移酶(EZH2)與STAT3的結(jié)合會(huì)使STAT3甲基化，進(jìn)而增強(qiáng)STAT3的活性。為此，本研究還分析了STAT3基因啟動(dòng)子區(qū)甲基化狀態(tài)及其與山羊體重的關(guān)系。結(jié)果發(fā)現(xiàn)，低體重組STAT3基因啟動(dòng)子區(qū)是非甲基化的，而高體重組STAT3基因啟動(dòng)子區(qū)為中度甲基化水平；統(tǒng)計(jì)分析揭示，低體重組與高體重組的STAT3基因甲基化修飾程度之間存在顯著差異，提示該基因甲基化修飾對(duì)體重有顯著影響，同時(shí)，還提示該基因的DNA甲基化修飾可以作為體重性狀篩選的有效表觀遺傳標(biāo)記。為此，建議將STAT3基因作為影響山羊產(chǎn)肉性能的重要候選基因，并有必要進(jìn)行深入研究。

上方數(shù)字代表源序列，下方數(shù)字代表重亞硫酸鹽處理后序列比對(duì)所用的模式序列。CG或YG代表各CpG二核苷酸位點(diǎn)，橫線指示引物位置The upward side sequence indicated origin sequence，while the downward side sequence indicated the bisulfite treatment sequence；CG or YG represented the location of CpG dinucleotide；the line indicated the location of primers圖9　STAT3-M1的 CpG位點(diǎn)的甲基化模式圖Fig.9　Methylated CpG sites mode of STAT3-M1 region

測(cè)序圖中帶框的字母CG或TG表示各CpG二核苷酸位點(diǎn)的位置；橫線指示引物位置CG or TG represented the location of CpG dinucleotide in sequence map；the line indicated the location of primers圖10　山羊STAT3基因的重亞硫酸鹽處理后測(cè)序圖Fig.10　The bisulfite sequencing results of goat STAT3 gene

圖11　不同體重組的山羊STAT3基因啟動(dòng)子區(qū)CpG島整體甲基化率比較Fig.11　Overall methylation percentage in the CpG island of the STAT3 gene promoter region between the high body weight group and the low body weight group

總之，本研究成功克隆了山羊STAT3基因CDS序列，并發(fā)現(xiàn)該基因啟動(dòng)子區(qū)甲基化修飾顯著影響山羊體重性狀，可作為體重性狀選擇的有效表觀遺傳標(biāo)記，這將有助于進(jìn)一步解析山羊產(chǎn)肉調(diào)控的遺傳與表觀遺傳機(jī)制，為未來(lái)山羊產(chǎn)業(yè)的發(fā)展提供科學(xué)資料。

參考文獻(xiàn)(References)：

[1]趙海諭，藍(lán)賢勇，雷初朝，等.表觀遺傳學(xué)：家畜遺傳育種的新挑戰(zhàn)[J].家畜生態(tài)學(xué)報(bào)，2014，35(8)：1-5.

ZHAOHY，LANXY，LEICZ，etal.Epigenetics：Newchallengesofbreedingandgeneticsinlivestock[J].Acta Ecologiae Animalis Domastici，2014，35(8)：1-5.(inChinese)

[2]QINH，YEHWI，DESARNOP，etal.Signaltransducerandactivatoroftranscription-3/suppressorofcytokinesignaling-3(STAT3/SOCS3)axisinmyeloidcellsregulatesneuroinflammation[J].Proc Natl Acad Sci U S A，2012，109(13)：5004-5009.

[3]ZHUANGS.RegulationofSTATsignalingbyacetylation[J].Cell Signal，2013，25(9)：1924-1931.

[4]DEMARIAM，CAMPOREALEA，POLIV.STAT3andmetabolism：Howmanywaystouseasinglemolecule?[J].Int J Cancer，2014，135(9)：1997-2003.

[5]SUBRAMANIAMA，SHANMUGAMMK，PERUMALE，etal.Potentialroleofsignaltransducerandactivatoroftranscription(STAT) 3signalingpathwayininflammation，survival，proliferationandinvasionofhepatocellularcarcinoma[J].Biochim Biophys Acta，2013，1835(1)：46-60.[6]KAMRANMZ，PATILP，GUDERP.RoleofSTAT3incancermetastasisandtranslationaladvances[J].Biomed Res Int，2013，2013：421821.

[7]OHKUBON，SUZUKIY，AOTOM，etal.AccelerateddestructionoferythrocytesinTie2promoter-drivenSTAT3conditionalknockoutmice[J]. Life Sci，2013，93(9-11)：380-387.

[8]RESEMANNHK，WATSONCJ，LLOYD-LEWISB.TheStat3paradox：Akillerandanoncogene[J].Mol Cell Endocrinol，2014，382(1)：603-611.

[9]RICHARDAJ，STEPHENSJM.TheroleofJAK-STATsignalinginadiposetissuefunction[J].Biochim Biophys Acta，2013，1842(3)：431-439.

[10]SCHüBELERD.FunctionandinformationcontentofDNAmethylation[J].Nature，2015，517(7534):321-326.

[11]PANCY，JIAWC，WUXF，etal.DNAmethylationprofileofDNAmethyltransferase3b(Dnmt3b)anditsinfluenceongrowthtraitsingoat[J].JAPS，2013，23(2)：380-387.

[12]趙海諭． 奶山羊PITX1、PITX2基因的遺傳變異、mRNA表達(dá)及DNA甲基化修飾對(duì)泌乳性能的影響[D].楊凌：西北農(nóng)林科技大學(xué)，2014.

ZHAOHY.Geneticvariation，mRNAexpressionandDNAmethylationmodificationofPITX1andPITX2genesandtheirinfluenceonmilkproductiontraitsindairygoats[D].Yangling：NorthwestA&FUniversity，2014.(inChinese)

[13]LANXY，CRETNEYEC，KROPPJ，etal.MaternaldietduringpregnancyinducesgeneexpressionandDNAmethylationchangesinfetaltissuesinsheep[J].Front Genet，2013，4：49.

[14]李盛杰，杜曉華，羅玉柱，等.天祝白牦牛NGB基因的克隆及生物信息學(xué)分析[J].畜牧獸醫(yī)學(xué)報(bào)，2013，44(3)：395-398.

LISJ，DUXH，LUOYZ，etal.CloningandbioinformaticsanalysisofNGBgeneofTianzhuwhiteyak[J].Acta Veterinaria et Zootechnica Sinica，2013，44(3)：395-398.(inChinese)

[15]王萍，巨向紅，雍艷紅，等.豬SOCS3基因的克隆、生物信息分析及熱應(yīng)激條件下的表達(dá)[J].畜牧獸醫(yī)學(xué)報(bào)，2014，45(7)：1044-1052.

WANGP，JUXH，YONGYH，etal.CloningandbioinformaticsanalysisofporcineSOCS3cDNAanditsexpressionpatternduringheatstress[J].Acta Veterinaria et Zootechnica Sinica，2014，45(7)：1044-1052.(inChinese)

[16]張春香，郭云雁，白元生，等.綿羊GPx5基因克隆、生物信息學(xué)分析及表達(dá)特性的研究[J].畜牧獸醫(yī)學(xué)報(bào)，2014，45(7)：1075-1083.

ZHANGCX，GUOYY，BAIYS，etal.Clone，bioinformaticsanalysisandexpressioncharacteristicsofovineGPx5[J].Acta Veterinaria et Zootechnica Sinica，2014，45(7)：1075-1083.(inChinese)

[17]陳靜，翟飛，夏明秀，等.雞不同生精細(xì)胞中Piwil1基因啟動(dòng)子區(qū)DNA甲基化差異研究[J].畜牧獸醫(yī)學(xué)報(bào)，2014，45(9)：1404-1409.

CHENJ，ZHAIF，XIAMX，etal.DifferencesinDNAmethylationofthepromoterregionofPiwil1geneinchickenspermatogeniccells[J].Acta Veterinaria et Zootechnica Sinica，2014，45(9)：1404-1409.(inChinese)

[18]李芳芳，戴建軍，吳彩鳳，等.WT1和FGF9基因在克隆豬隱睪中mRNA表達(dá)和DNA甲基化狀態(tài)分析[J].畜牧獸醫(yī)學(xué)報(bào)，2014，45(6)：918-925.

LIFF，DAIJJ，WUCF，etal.mRNAexpressionandDNAmethylationstatusofWT1andFGF9geneincryptorchidismofclonedpig[J]. Acta Veterinaria et Zootechnica Sinica，2014，45(6)：918-925.(inChinese)

[19]王夢(mèng)楠，崔亞麗，吳國(guó)江，等.牛Ascl2基因印記及DNA甲基化狀態(tài)分析[J].畜牧獸醫(yī)學(xué)報(bào)，2014，45(12)：1949-1956.

WANGMN，CUIYL，WUGJ，etal.GenomicimprintingandDNAmethylationstatusofAscl2geneincattle[J]. Acta Veterinaria et Zootechnica Sinica，2014，45(12)：1949-1956.(inChinese)

[20]張文，王雅杰，薛春燕，等.STAT3基因在乳腺癌組織中的表達(dá)及意義[J].中國(guó)癌癥雜志，2007，17(7):546-551.

ZHANGW，WNAGYJ，XUECY，etal.ExpressionofSTAT3geneinbreastcarcinomaanditsclinicalsignificance[J].China Oncology，2007，17(7):546-551.(inChinese)

[21]DAVEYHW，WILKINSRJ，WAXMANDJ.STAT5signalinginsexuallydimorphicgeneexpressionandgrowthpatterns[J].Am J Hum Genet，1999，65(4):959-965.

[22]WUX，JIAW，ZHANGJ，etal.DeterminationofthenovelgeneticvariantsofgoatSTAT5AgeneandtheireffectsonbodymeasurementtraitsintwoChinesenativebreeds[J].Small Rumin Res，2014，121(2-3)：232-243.

[23]JIAW，WUX，LIX，etal.NovelgeneticvariantsassociatedwithmRNAexpressionofsignaltransducerandactivatoroftranscription3(STAT3)genesignificantlyaffectedgoatgrowthtraits[J].Small Rumin Res，2015，129：25-36.

[24]JONESPA.FunctionsofDNAmethylation：islands，startsites，genebodiesandbeyond[J].Nat Rev Genet，2012，13(7)：484-492.

[25]KUDOS.Methyl-CpG-bindingproteinMeCP2repressesSp1-activatedtranscriptionofthehumanleukosialingenewhenthepromoterismethylated[J].Mol Cell Biol，1998，18(9)：5492-5499.

[26]YANGJ，BAIWL，NIUPY，etal.AberranthypomethylatedSTAT3wasidentifiedasabiomarkerofchronicbenzenepoisoningthroughintegratingDNAmethylationandmRNAexpressiondata[J].Exp Mol Pathol，2014，96(3)：346-353.

[27]KIME，KIMM，WOODH，etal.PhosphorylationofEZH2activatesSTAT3signalingviaSTAT3methylationandpromotestumorigenicityofglioblastomastem-likecells[J].Cancer Cell，2013，23(6)：839-852.

(編輯郭云雁)

Molecular Cloning，Bioinformatics Analysis and Methylation Modification of GoatSTAT3 Gene

JIA Wen-chao1#，LIU Liang-liang2#，WU Xian-feng1，ZHAO Zhao-yan1，WANG Ke1，WANG Mao1，GAO Jing1，WANG Li-qiang1，CHEN Hong1，PAN Chuan-ying1，LAN Xian-yong1*

(1.ShaanxiKeyLaboratoryofMolecularBiologyforAgriculture，CollegeofAnimalScienceandTechnology/CollegeofLifeSciences，NorthwestA&FUniversity，Yangling712100，China；2.ComputerScienceandEngineeringSchoolofJust，JiangsuUniversityofScienceandTechnology，Zhenjiang212003，China)

Abstract:The aim of this study was to clone and analyze the entire coding sequence of goatSTAT3 gene as well as to detect theSTAT3 gene methylation modification and analyze its association with goat body weight using reverse transcription-polymerase chain reaction(RT-PCR)，TA cloning，bioinformatics and bisulfite sequencing(BSP) methods.The results showed that the entire coding sequence of goatSTAT3 gene was 2 313 bp in length，and encoded 770 amino acids.The shared amino acid sequence consistency between goat(caprahircus) and sheep(Ovisaries)，bovine(BosTaurus)，boar(Susscrofa)，mouse(Musmusculus)，rat(Rattusnorvegicus)，human(Homosapiens) were 99.7%，99.6%，99.4%，99.2%，99.4%，99.5%，respectively.The analysis of the methylation difference in the promoter region of goatSTAT3 gene between high body weight group and low body weight group demonstrated that the goats with high body weight had higher methylation level than those goats with low body weight(P=0.020).This finding indicated that the methylation modification ofSTAT3 gene had significant effect on body weight in goat，which could become the epigenetic marker in marker-assisted selection(MAS).These results would provide scientific data for studying genetic and epigenetic mechanism for meat production in goat.

Key words:goat；STAT3 gene；molecular cloning；sequence analysis；DNA methylation

doi:10.11843/j.issn.0366-6964.2016.03.006

收稿日期：2015-03-30

基金項(xiàng)目：國(guó)家自然科學(xué)基金(31172184)；陜西省留學(xué)人員科技活動(dòng)擇優(yōu)資助項(xiàng)目(2014-14)；西北農(nóng)林科技大學(xué)大學(xué)生創(chuàng)新性實(shí)驗(yàn)計(jì)劃(2013)

作者簡(jiǎn)介：賈文超(1990-)，男，土族，青海門源人，碩士生，主要從事動(dòng)物遺傳學(xué)研究，E-mail：jiawenchaozyq@163.com；劉亮亮(1979-)，男，江西萍鄉(xiāng)人，博士，講師，主要從事計(jì)算機(jī)軟件與信息分析，E-mail：lingyun79626@126.com。賈文超和劉亮亮為并列第一作者 *通信作者：藍(lán)賢勇，博士，副教授，碩士生導(dǎo)師，主要從事動(dòng)物遺傳育種與表觀遺傳修飾研究，E-mail：lanxianyong79@nwsuaf.edu.cn

中圖分類號(hào)：S827；S813.3

文獻(xiàn)標(biāo)志碼：A

文章編號(hào):0366-6964(2016)03-0457-10