蘭州大尾羊MAPK13基因cDNA克隆及生物信息學分析

金方園徐紅偉達小強臧榮鑫柏家林曹忻蔡勇馮玉蘭楊具田

（1. 西北民族大學生命科學與工程學院，蘭州 730030；2. 西北民族大學實驗中心，蘭州 730030）

蘭州大尾羊MAPK13基因cDNA克隆及生物信息學分析

金方園1徐紅偉2達小強1臧榮鑫1柏家林1曹忻2蔡勇2馮玉蘭1楊具田2

（1. 西北民族大學生命科學與工程學院，蘭州 730030；2. 西北民族大學實驗中心，蘭州 730030）

克隆蘭州大尾羊促分裂素原活化蛋白激酶MAPK13基因，并分析其序列及其編碼蛋白的生物學特性，為研究綿羊MAPK13基因的功能和生產(chǎn)應用提供參考。根據(jù)綿羊MAPK13基因CDS序列設計特異引物，利用RACE和RT-PCR技術(shù)克隆獲得蘭州大尾羊MAPK13基因全長序列，并結(jié)合生物信息學方法分析其生物學特性?？寺～@得蘭州大尾羊MAPK13基因cDNA 序列全長1 397 bp，其CDS區(qū)片段長1 102 bp，編碼367個氨基酸。預測蘭州大尾羊MAPK13蛋白分子量為42.29 kD，理論等電點為8.82，為非跨膜的疏水性蛋白，亞細胞定位主要在細胞質(zhì)中，無信號肽，不屬于分泌蛋白。預測其氨基酸序列有19個磷酸化位點，3個糖基化位點，3個磷酸化功能結(jié)構(gòu)域，4個其他結(jié)構(gòu)域，1個LCR片段，二級結(jié)構(gòu)以α-螺旋為主。同源性分析顯示蘭州大尾羊MAPK13基因序列與已發(fā)布的綿羊MAPK13 mRNA序列（登錄號：NM_001139455.1）相比，其第852位發(fā)生堿基轉(zhuǎn)換（C?A），導致編碼蛋白第265位氨基酸發(fā)生堿基轉(zhuǎn)換（S?T），同時，第951位發(fā)生堿基轉(zhuǎn)換（T?G），但其所編碼氨基酸不變。構(gòu)建的基因進化樹分析結(jié)果顯示蘭州大尾羊與牛親緣關系最近。蘭州大尾羊與其他物種MAPK13基因在結(jié)構(gòu)上相似性較高，說明該基因具有高度的保守性，其序列包含的S_TKc結(jié)構(gòu)域可將ATP的γ磷?；D(zhuǎn)移到蛋白質(zhì)絲氨酸/蘇氨酸殘基上，導致一系列肥胖相關基因的功能失調(diào)和表達變化，為進一步研究MAPK13基因與成脂分化過程的相關性提供了參考。

蘭州大尾羊 MAPK13基因 cDNA 末端快速擴增 生物信息學

促分裂素原活化蛋白激酶（Mitogen-activated protein kinases，MAPK）鏈是真核生物信號傳遞網(wǎng)絡中的重要途徑之一［1］。MAPK通路主要包括ERK、p38MARK和JNK三條途徑，主要作用于脂肪細胞組織中，其逐級磷酸化過程［2］涉及成脂分化的調(diào)控，與胰島素抵抗、脂肪細胞分化密切相關［3］。p38 delta MAPK也稱為MAPK13，是脂肪細胞分化的主要候選基因。蘭州大尾羊（Lanzhou fat-tailed sheep）是清朝同治年間（1 862-1 875年）由陜西大荔一帶引進的同羊與蘭州當?shù)孛晒叛螂s交選育而成的地方綿羊品種，因其尾部脂肪過度沉積，導致尾大、脂肪充實而得名［4，5］。綿羊MAPK13基因克隆及其相關功能研究報道較少。1997年，Cohern等［6］用Ser/Thr磷酸酶處理使MAPKK失活，提出了存在MAPKK激酶的假說。Emanuelli［7］研究證實，胰島素抵抗，二型糖尿病和肥胖均與MAPK有關。2005年，Bost等［8］研究發(fā)現(xiàn)，ERK小鼠脂肪細胞含量少于野生型小鼠，與對照組相比，高脂飼料喂養(yǎng)的ERK小鼠有較高的餐后代謝率，不易患胰島素抵抗和肥胖。2009年，趙芳等［3］研究表明，ERK1/2，JNK和p38MAPK三條MAPK信號通路通過獨立作用及復雜的協(xié)調(diào)作用，參與調(diào)控脂肪細胞分化過程的各個環(huán)節(jié)，其調(diào)節(jié)作用具有多樣性。2013年，呂雯［9］脂肪細胞等研究表明，ERK和p38MAPKK兩種MAPK信號通路對脂分化的調(diào)節(jié)在不同的試驗模型中表現(xiàn)為正調(diào)控和負調(diào)控兩種不同形式；而另一成員JNK能使胰島素受體底物1的絲氨酸發(fā)生磷酸化，進而干擾胰島素信號，從而抑制骨髓間充質(zhì)干細胞（BMSCs）的成脂分化，即對脂肪細胞分化發(fā)揮負調(diào)控作用。關于MAPK通路對脂肪細胞分化發(fā)揮正調(diào)控作用，Haridas Nidhina［10］研究表明，香草精通過激活ERK通路，上調(diào)了前體脂肪細胞中PPARg的表達，增加了細胞葡萄糖的攝取量，促進了成脂分化。關于MAPK通路對脂肪細胞分化發(fā)揮負調(diào)控作用，Chuang［11］研究表明，ERF信號通路可以誘導白細胞介素表達，使其下游PPARg被磷酸化，導致PPARg下游生脂相關的靶基因表達水平降低，從而抑制了脂肪細胞分化。此外，Bordicchia等［12-15］的研究，也表明MAPK基因與脂肪代謝關系密切。基因調(diào)控是細胞分化的核心問題。目前，雖然有關MAPK13基因的研究在報道中較多，但蘭州大尾羊是具有耐粗飼、抗病能力強、飼料利用率高、適應性強、肉質(zhì)鮮美等優(yōu)良地方綿羊品種，其尾部脂肪細胞分化與其它綿羊差異顯著。迄今為止，MAPK13等脂肪細胞分化基因克隆及表達的相關報道較少［16］。脂肪細胞分化過程與脂肪細胞特異性功能基因的表達調(diào)控密切相關，涉及多種脂肪特異性轉(zhuǎn)錄因子的激活。本研究以蘭州大尾羊為研究對象，應用RT-PCR技術(shù)獲得cDNA，利用RACE法對MAPK13基因進行克隆和序列分析，為進一步研究MAPK13基因?qū)μm州大尾羊尾部脂肪細胞分化的影響提供參考。

1 材料與方法

1.1 材料

1.1.1 試驗動物 2012年3月于永登明鑫良種肉羊繁育場選購6月齡蘭州大尾羊（羯羊），頸靜脈放血后，快速切取尾部脂肪組織，切成小塊，用錫箔紙包裹后液氮速凍，-80℃保存?zhèn)溆?。蘭州大尾羊脂尾肥大，方圓平展，自然下垂飛節(jié)上下，尾有中溝將尾部分為左右對稱兩瓣，尾尖外翻，并緊貼中溝，尾面著生被毛，內(nèi)面光滑無毛，呈淡紅色。由于該品種瀕臨滅絕，純種存欄不足百只，利用正常繁育羊只開展研究工作受到了限制，因而屠宰了1只羯羊進行了研究。該羊尾脂重1.71 kg，占體重的3%，具有理想羊只的特征。以下實驗于西北民族大學實驗中心進行。

1.1.2 主要試劑 RNAiso Plus（TaKaRa）、SMART-erTMRACE cDNA Amplification Kit（Clontech）、DH5α感受態(tài)細胞（TaKaRa）、普通瓊脂糖凝膠DNA回收試劑盒、pGM-T克隆試劑盒、質(zhì)粒小提試劑盒、X-gal、IPTG、DNA Marker，6×loading buffer、胰蛋白胨、酵母提取物、瓊脂粉、氯仿、異丙醇、0.1%DEPC、75%乙醇、RNase-free水、5×TBE 緩沖液、瓊脂糖、氯化鈉等均購自天根生物公司。

1.2 方法

1.2.1 脂肪組織總RNA的提取 通過TRIzol法提取總 RNA。

1.2.2 引物設計與合成 根據(jù)綿羊MAPK13mRNA序列（登錄號：NM_001139455.1），按照SMARTerTMRACE cDNA Amplification Kit要求［17］，利用Primer 5.0軟件設計引物，由上海生工合成。公司測序。

1.2.5 數(shù)據(jù)分析 利用DNASTAR軟件對已獲得的MAPK13序列片段進行拼接，對拼接的cDNA序列進行分析并推導其氨基酸序列；利用ProtParam（http：//us.expasy.org/tools/protparam.htmL）分析氨基酸序列的理化性［18］；利用IBCP（http：//npsa-pbil. ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_gor4.html）預測蛋白質(zhì)的二級結(jié)構(gòu)［19］；利用PHYRE2（http：// www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index）預測蛋白質(zhì)的三級結(jié)構(gòu)［20］；利用 TMHMM（http：// genome.cbs.dtu.dk/services /TMHMM/）分析蛋白質(zhì)的跨膜結(jié)構(gòu)［21］；利用Protscale（http：//web.expasy. org/protscale/）分析蛋白質(zhì)親疏水性；利用NetPhos（http：//www.cbs.dtu.dk/services/NetPhos/）分析蛋白質(zhì)磷酸化位點；利用NetNGlyc（http：//www.cbs.dtu. dk/services/NetNGlyc/）分析蛋白質(zhì)糖基化位點；利用 TargetP（http：//www.cbs.dtu.dk/services/TargetP/）分析蛋白潛在信號肽；利用SignalP4.0（http：//www. cbs.dtu.dk/services/SignalP/）分析蛋白分泌途徑；利用 PSORT II Prediction（http：//psort.hgc.jp/form.html）預測蛋白質(zhì)的亞細胞定位；利用 Smart（http：//smart. embl-heidelberg.de/）分析預測氨基酸序列的保守結(jié)構(gòu)域；利用Blast和DNAman軟件分析其同源性及構(gòu)建系統(tǒng)發(fā)育樹［22］。

1.2.3 cDNA第一鏈合成 cDNA 第一條鏈的合成根據(jù)SMARTerTMRACE cDNA Amplification Kit說明書進行操作。

1.2.4 RACE反應及產(chǎn)物回收 根據(jù)SMARTerTMRACE cDNA Amplification Kit中 3'-Full RACE Core Set 和 5'-Full RACE Core Set 2 5'-full race試劑盒進行蘭州大尾羊MAPK13的 3'端與 5'端序列的擴增，具體操作步驟參考產(chǎn)品說明書，用CDS區(qū)引物CDS1和 CDS2按常規(guī)方法進行 RT-PCR 擴增。反應結(jié)束后各取 3'-RACE PCR 產(chǎn)物、5'-RACE PCR 產(chǎn)物和CDS區(qū)PCR 產(chǎn)物5.0 μL，用 12 g/L瓊脂糖凝膠電泳。目標3'-RACE cDNA、5'-RACE cDNA和CDS區(qū)PCR產(chǎn)物經(jīng)DNA凝膠回收試劑盒進行回收和純化，連接至pDM19-T載體，提取質(zhì)粒并送上海生工生物工程

2 結(jié)果

2.1 蘭州大尾羊MAPK13基因cDNA末端快速擴增（RACE）

根據(jù)SMARTerTMRACE cDNA擴增試劑盒中5'-full race試劑盒和3'-full race試劑盒說明書進行的操作步驟。NGSP1擴增得到90 bp 5'-RACE cDNA，NGSP2擴增得到741 bp 3'-RACE cDNA；CDS區(qū)引物CDS1和CDS2擴增獲得1 102 bp中間片段（圖1）。

2.2 蘭州大尾羊MAPK13基因cDNA 序列分析

將中間片段、5'-RACE 和 3'-RACE序列測序拼接后獲得蘭州大尾羊MAPK13基因1 397 bp的全長cDNA序列。MAPK13基因核苷酸序列及其編碼蛋白的氨基酸序列如圖2所示，蘭州大尾羊MAPK13基因的開放閱讀框為1 101 bp，起始密碼子為ATG位于58 nt，終止密碼子為TAG位于1 158 nt，編碼367個氨基酸。編碼區(qū)兩側(cè)翼分別具有5'-UTR 57 bp（1-57 nt）和3'-UTR 214 bp（1 158-1 397 nt）。其中C+G%=57.55%，分子量42.29 kD。

2.3 蘭州大尾羊MAPK13蛋白質(zhì)結(jié)構(gòu)分析

2.3.1 基因蛋白理化性質(zhì)分析 DNAstar軟件分析表明，蘭州大尾羊MAPK13基因cDNA編碼蛋白質(zhì)367個氨基酸殘基中，有52個堿性氨基酸（K，R），占氨基酸總數(shù)的14.2%；47個酸性氨基酸（D，E），占氨基酸總數(shù)的12.8%；124個疏水氨基酸（A，I，F(xiàn)，W，V），占氨基酸總數(shù)的33.8%；87個極性氨基酸（N，C，Q，S，T，Y），占氨基酸總數(shù)的23.7%。

2.3.2 基因蛋白二級結(jié)構(gòu)與三級結(jié)構(gòu)預測 ProtParam程序分析表明蘭州大尾羊MAPK13蛋白的的分子量為分子量42.29 kD，理論等電點為8.82。蘭州大尾羊MAPK13蛋白中，α-螺旋（Alpha helix）占41.26%；隨機卷曲（Random coil）占39.89%；延伸鏈（Extended strand）占18.85%（圖2）。通過PHYRE2預測其三級結(jié)構(gòu)（圖3）。

2.3.3 MAPK13蛋白的保守結(jié)構(gòu)域 根據(jù) Smart軟件推測，MAPK13蛋白質(zhì)其第306-319個氨基酸殘疾有一段LCR序列（表2，圖4）。

2.3.4 蛋白基因跨膜結(jié)構(gòu)與親疏水性分析 TMHMM分析也表明蘭州大尾羊MAPK13蛋白無跨膜螺旋，為非跨膜蛋白。利用ExPASy的Protscale分析蛋白疏水性（圖5）發(fā)現(xiàn)，該基因編碼蛋白疏水性最大值為2.078，最小值為-3.044。該氨基酸序列內(nèi)絕大多數(shù)為疏水性殘基，表明綿羊MAPK13基因編碼蛋白為不溶性蛋白。

2.3.5 亞細胞定位 根據(jù)TargetP 和SignalP 軟件預測，推測該蛋白無信號肽，不是分泌蛋白，大部分在游離核糖體上合成（表3）。將MAPK13蛋白序列輸入 PSORT II Prediction，可以基本確認該基因蛋白不屬于分泌蛋白，主要在細胞質(zhì)（45%）和過氧化物酶體（36.5%）中發(fā)揮生物學作用，少量作用于溶酶體（10%）和線粒體基質(zhì)（10%）（表4）。

2.3.6 蛋白基因磷酸化與糖基化位點分析 根據(jù)NetPhos 2.0 Server軟件分析，結(jié)果（圖6）表明MAPK13蛋白含有12個絲氨酸Ser，3個蘇氨酸Thr，4個酪氨酸Tyr可以成為蛋白質(zhì)激酶磷酸化的位點。根據(jù)NetNGlyc 1.0 Server軟件分析，結(jié)果表明其含有3個糖基化位點：N15、N201、N347。

2.4 蘭州大尾羊MAPK13基因及其編碼蛋白與其他動物同源性比較

用Blast進行核苷酸同源性在線分析，結(jié)果表明蘭州大尾羊MAPK13基因核苷酸序列與牛、人、海象、犬、猴、虎鯨、海豚和鼠核苷酸序列間的同源性分別為97%、92%、91%、91%、91%、90%、90%和88%，說明蘭州大尾羊與綿羊和牛的親緣關系較近，與鼠類親緣關系最遠。用DNAman軟件進行氨基酸同源性分析，結(jié)果（圖7）表明，蘭州大尾羊MAPK13氨基酸序列與牛、猴、犬、海象、海豚、虎鯨、鼠和人核苷酸序列間的同源性分別為99.2%、87.8%、94.8%、93.7%、91.8%、92.1%、92.1%和81.1%，其中與牛同源性較高，與鼠和人類同源性較差。系統(tǒng)發(fā)育樹（圖7）也表明，蘭州大尾羊與牛的進化水平更為接近，與人類較遠。蘭州大尾羊MAPK13基因序列與已發(fā)布的綿羊MAPK13mRNA序列（登錄號：NM_001139455.1）相比第852位發(fā)生堿基轉(zhuǎn)換（C←→A），導致編碼蛋白第265位氨基酸發(fā)生堿基轉(zhuǎn)換（S←→T），同時，第951位發(fā)生堿基轉(zhuǎn)換（T←→G），但其所編碼氨基酸不變（圖8）。

3 討論

本研究通過RT-PCR和RACE技術(shù)首次克隆了蘭州大尾羊MAPK13基因cDNA 全長序列1 397 bp，其CDS區(qū)片段長1 102 bp，編碼367個氨基酸。蘭州大尾羊MAPK13基因核苷酸序列與牛、人、海象、犬、猴、虎鯨、海豚和鼠核苷酸序列間的同源性分別為97%、92%、91%、91%、91%、90%、90%和88%。蘭州大尾羊MAPK13基因序列與已發(fā)布的綿羊MAPK13mRNA序列（登錄號：NM_001139455.1）相比，其第852位發(fā)生堿基轉(zhuǎn)換（C←→A），導致編碼蛋白第265位氨基酸發(fā)生堿基轉(zhuǎn)換（S←→T），同時，第951位發(fā)生堿基轉(zhuǎn)換（T←→G），但其所編碼氨基酸不變。系統(tǒng)發(fā)育樹符合進化規(guī)律。由此可見，蘭州大尾羊與各物種MAPK13基因序列間具有高度的保守性，從而預測生物體內(nèi)由MAPK13基因所指導的生物學功能非常穩(wěn)定。

預測MAPK13分子量為42.49 kD，為非跨膜的疏水性蛋白，無信號肽，不屬于分泌蛋白。亞細胞定位主要在細胞質(zhì)（45%）和過氧化物酶體（36.5%）中發(fā)揮生物學作用，少量作用于溶酶體（10%）和線粒體基質(zhì)（10%）。綜上推測MAPK13蛋白大部分在游離核糖體上合成，其功能與過氧化物酶體有關。該預測結(jié)果符合其促分裂素原活化蛋白激酶的身份。

MAP級聯(lián)反應主要通過絲/蘇/酪氨酸的磷酸化作用推進。S_TKc結(jié)構(gòu)域由十幾個保守半胱氨酸殘基組成，具有絲/蘇氨酸激酶催化功能，經(jīng)上游MAPK激酶MKK3和MKK6激活，可將ATP的γ磷?；D(zhuǎn)移到蛋白質(zhì)絲/蘇氨酸殘基上，是細胞外信號反應的重要介質(zhì)［24］。其催化核心的N端存在一個LLLL富集區(qū)，這個區(qū)域可能與ATP的結(jié)合有關。TyrKc酪氨酸激酶結(jié)構(gòu)域與S_TKc功能基本相同，包含活性位點、ATP結(jié)合位點和底物結(jié)合位點，可催化ATP的γ磷酸基團轉(zhuǎn)移到專一性底物如蛋白絲氨酸、蘇氨酸或酪氨酸殘基的羥基基團［24］。MAPK級聯(lián)放大過程中存在一個Ras-Raf-MAPKKMAPK信號轉(zhuǎn)導通路，該通路則需要RasGap結(jié)構(gòu)域的參與，Ras在PTKs介導的信號通路中也是一種關鍵組分。Ras蛋白是ras基因表達產(chǎn)物，具有GTPase活性，分布于質(zhì)膜胞質(zhì)一側(cè)，結(jié)合GTP時為活化狀態(tài)，而結(jié)合GDP時為失活態(tài)，所以Ras蛋白也是GTPase開關蛋白［25］。Yamakawa［26］通過對逆轉(zhuǎn)錄病毒法處理的小鼠B細胞進行篩選，發(fā)現(xiàn)RasGAP結(jié)合蛋白Dok-1作為一個酪氨酸磷酸化激活分子作用于其途徑的下游，但其本身并沒有被磷酸化，從而證明Dok-1是連接激活素受體與Smad蛋白的重要分子。Smads家族蛋白可以將TGF-β信號從細胞表面受體傳導至細胞核，因此，RasGAP結(jié)構(gòu)域間接地激活或抑制了靶基因的轉(zhuǎn)錄［27］。預測MAPK13氨基酸序列有19個可以成為蛋白質(zhì)激酶磷酸化的位點，某些位點被磷酸化后，可導致一系列肥胖相關基因的功能失調(diào)和表達變化，如胰島素增敏激素，脂聯(lián)素（adiponectin）的表達減少［28，29］，從而推測其磷酸化反應過程可能與脂肪代謝相關聯(lián)。保守結(jié)構(gòu)域預測結(jié)果分析表明，蘭州大尾羊MAPK13蛋白序列中含有3個磷酸化功能結(jié)構(gòu)域：S_TKc絲/蘇氨酸激酶結(jié)構(gòu)域、TyrKc酪氨酸激酶結(jié)構(gòu)域和Rio絲氨酸激酶結(jié)構(gòu)域，1個RasGAP信號傳導結(jié)構(gòu)域，1個CPSase_sm谷氨酰胺氨基轉(zhuǎn)移酶結(jié)構(gòu)域，1個IENR1核酸內(nèi)切酶結(jié)構(gòu)域，1個MgtE_N Mg2+轉(zhuǎn)運結(jié)構(gòu)域。說明MAPK13基因可能通過磷酸化調(diào)控蘭州大尾羊脂肪代謝和局部沉積，是研究蘭州大尾羊尾部脂肪細胞分化的重要候選基因。

4 結(jié)論

應用RACE技術(shù)首次克隆了蘭州大尾羊MAPK13基因1 397 bp全長cDNA序列（登錄號：KF770837），并對其核苷酸與編碼氨基酸序列特性進行了研究，為進一步研究綿羊MAPK13基因功能提供參考。

［1］范衡宇, 佟超, 孫青原. 絲裂原活化蛋白激酶（MAPK）信號通路的研究進展［J］. 動物學雜志, 2002, 37（5）：98.

［2］劉玉鯤, 李德全. MAPK級聯(lián)途徑調(diào)控植物細胞胞質(zhì)分裂［J］.中國生物化學與分子生物學報, 2008, 24（12）：1085-1091.

［3］趙芳, 孫瑩璞. MAPK信號轉(zhuǎn)到通路及其在脂肪細胞分化中的作用［J］. 國際生殖健康/計劃生育雜志, 2009, 28（5）：301-304.

［4］ 夏吉鵬, 張永升, 孟圣杰, 等. 瀕危物種蘭州大尾羊調(diào)查分析［J］. 中國科技論文在線http：//www.paper.edu.cn, 2006, 5.

［5］丁紹殿, 吳筱華, 李性明. 蘭州大尾羊調(diào)查報告［J］. 中國畜牧雜志, 1986, 1：20-21.

［6］ Cohen P. The search for physiological substrates of MAP and SAP kinases in mammalian cells［J］. Trands Cell Biol, 1997, 7（2）：353-361.

［7］ Emanuelli B, Eberle D, Suzuki R, et al. Overexpression of the dualspecificity phosphatase MKP-4/DUSP-9 protects against stressinduced insulin resistance［J］. Proc Natl Acad Sci USA, 2008, 105（9）：3545-3550.

［8］ Bost F, Aouadi M, Caron L, et al. The extracellular signal-regulated kinase isoform ERK1 is specifically required for in vitro and in vivo adipogenesis［J］. Diabetes, 2005, 54（2）：402-411.

［9］ 呂雯, 盧榮勝, 李捷. MAPK信號通路對脂肪細胞分化的調(diào)控［J］. 動物醫(yī)學進展, 2012, 33（3）：110-114.

［10］ Haridas Nidhina PA, Poulose N, Gopalakrishnapillai A. Vanillin induces adipocyte differentiation in 3T3-L1 cells by activating extracellular signal regulated kinase 42/44［J］. Life sci, 2011, 88（15-16）：675-680.

［11］ Chung S, Brown JM, Provo JN, et al. Conjugated linoleic acid promotes human adipocyte insulin resistance through NFκB-dependent cytokine production［J］. J Biol Chem, 2005, 280（46）：38445-38456.

［12］ Bordicchia M, Liu DX, Amri EZ, et al. Cardiac natriuretic peptides act via p38 MAPK to induce the brown fat thermogenic program in mouse and human adipocytes［J］. J Clin Invest, 2012, 122（3）：1022-1036.

［13］ Moule SK, Denton RM, et al. The activation of p38 MAPK by the β-adrenergic agonist isoproterenol in rat epiddymal fat cells［J］. University Bristol, 1998, 439（3）：287-290.

［14］ Sun C, Wang L, Yan J, et al. Calcium ameliorates obesity induced by high-fat diet and its potential correlation with p38 MAPK pathway［J］. Molecular Biology Reports, 2012, 39（2）：1755-1763.

［15］ Yoon MJ, Lee GY, Chung JJ, et al. Adiponectin increases fatty acid oxidation in skeletal muscle cells by sequential activation of AMP-activated protein kinase, p38 mitogen-activated protein kinase, and peroxisome proliferator-Activated Receptor α［J］. Diabetes, 2006, 55（9）：2562-2570.

［16］徐紅偉, 柏家林, 馮玉蘭, 等. 蘭州大尾羊心臟型脂肪酸結(jié)合蛋白（H-FABP）基因克隆及其同源性比較［J］. 中國農(nóng)業(yè)科學, 2013, 46（3）：639-646.

［17］薛霖莉, 董常生, 赫曉燕, 等. 羊駝垂體催乳素（PRL）基因全長cDNA的克隆及序列分析［J］. 中國農(nóng)學通報, 2011, 27（1）：356-361.

［18］Wilkins MR, Gasteiger E, Bairoch A, et al. Protein identification and analysis tools in the ExPASy server［J］. Methods Mol Biol, 1999, 112：531-552.

［19］Geourjon C, Deleage G. SOPM：a self-optimized method for protein secondary structure prediction［J］. Protein Eng, 1994, 7：157-164.

［20］Tyagi N, Srinivasan N, et al. Recognition of nontrivial remote homology relationships involving proteins of Helicobacter pylori：implications for function recognition［J］. Methods in Molecular Biology, 2013, 993：155-175.

［21］Moller S, Croning MD, Apweiler AR. Evaluation of methods for the prediction of membrane spanning regions［J］. Bioinformatics, 2001, 17（7）：646-653.

［22］Zhang B, Huang K, Guo J, et al. Origin and the recombinant model of H7N9 virus prevailing in China［J］. Nan Fang YiKe Da Xue Xue Bao, 2013, 33（7）：1017-1021.

［23］鄒先瓊. 人類基因C4orf13、CDK5RAP1_v3、CDK5RAP1_v4和PXK的克隆和功能研究［D］. 廣州：中山大學, 2005.

［24］ 梁聰, 彭景琨, 曹廣力, 等. 家蠶hop基因的克隆及序列分析［J］. 江蘇農(nóng)業(yè)科學, 2012, 40（1）：10-13.

［25］Rudack T, Xia F, Schlitter J. Ras and GTPase-activating protein（GAP）drive GTP into a precatalytic state as revealed by combining FTIR and biomolecular simulations［J］. Proc Natl Sci USA, 2012, 109（38）：15295-15300.

［26］Yamakawa N, Tsuchida K, Sugino H. The rasGAP-binding protein, Dok-1, mediates activin signaling via serine/threonine kinase receptors［J］. EMBO J, 2002, 21（7）：1684-1694.

［27］ 孫琰, 葉棋濃. Smad4蛋白的研究進展［J］. 生物技術(shù)通訊, 2005, 16（3）：299-302.

［28］ 羅森. 抗糖尿病藥通過Cdk5抑制肥胖相關的PPARγ磷酸化［J］. 中國病理生理雜志, 2010, 446（7305）：451-456.

［29］ 李詠, 楊俊卿. PPARγ——中樞神經(jīng)系統(tǒng)損傷治療的新靶點［J］.中國老年學雜志, 2009（18）：2405-2407.

（責任編輯 李楠）

Cloning and Bioinformatics Analysis of MAPK13 Gene in Lanzhou Fat-tailed Sheep

Jin Fangyuan1Xu Hongwei2Da Xiaoqiang1Zang Rongxin1Bai Jialin1Gao Xin2Cai Yong2Feng Yulan1Yang Jutian2
（1. College of Life Science and Engineering，Northwest University for Nationalities，Lanzhou 730030；2. Science Experimental Center，Northwest University for Nationalities，Lanzhou 730030）

The objective of this study is to reveal the biological function of mitogen-activated protein kinases（MAPK）gene in sheep and provide theoretical data for application by cloning MAPK and analyzing its sequence as well as genetic features. The specific primers on the CDS template of MAPK13 gene were designed and then the MAPK13 gene sequence of Lanzhou fat-tail sheep was cloned using RACE and RTPCR technology. After that, its genetic characteristics were analyzed by bioinformatics methods. The 1 397 bp full-length cDNA of MAPK13 of Lanzhou fat-tailed sheep was obtained, in which the length of CDS is 1 102 bp. It encodes 367 amino acids in total. It was predicted that the molecular weight of MAPK13 protein of Lanzhou fat-tail sheep is 42.29 kD and its theoretical isoelectric point is 8.82. It is a non-transmembrane hydrophobin with no signal peptide, which locates mostly in cytoplasm by the result of subcellular level prediction and does not belong to the secretory protein. Its amino acid sequence contains 19 phosphorylation sites, 3 glycosylation sites, 3 phosphorylation domains, 4 other domains and 1 LCR domain. Also, its secondary structure is mainly randomly curly. Nucleotide sequence analysis revealed that the gene nucleotidesequence of Lanzhou fat-tail sheep has some different from the known sheep MAPK13 mRNA（Genbank No.：NM_001139455.1）. Base transition occur in the 852th site is（C vs. A）, respectively, the corresponding amino acids in the 265th is（S vs. T）. Base transition occur in the 951th site is（T vs. G）, but the corresponding amino acids do not change.The phylogenetic tree indicates that Lanzhou fat-tail sheep is close to bos. The fact that the structure of MAPK13 of Lanzhou fat-tail sheep highly similar with other species indicates that the gene is highly conservative. The S_TKc domain in the sequence can transfer the gamma phosphoryl of ATP to the serine / threonine residues of protein, leading to function imbalance and expression change in a series of obesity related gene.

Lanzhou fat-tailed sheep MAPK13 gene RACE Bioinformatics analysis

2014-01-23

國家自然科學基金項目（31160440，31260533，31360529）,甘肅省科技支撐計劃項目（1011NKCA051）,蘭州市科技計劃項目（2011-1-113），西北民族大學研究生科研創(chuàng)新項目（ycx13175）

金方園，E-mail：373995197@qq.com

楊具田，教授，碩士生導師，E-mail：jutianyang988@163.com；臧榮鑫，教授，碩士生導師，E-mail：rongxinzang2000@163.com