宋如昕，王蘭，張左兵

山西大學 生命科學學院，山西 太原 030006

中華鱉gp130基因全長cDNA克隆及生物信息學分析

宋如昕，王蘭，張左兵

山西大學 生命科學學院，山西 太原 030006

目的：從中華鱉脾臟、肝臟、腸道構建的SMART cDNA文庫中克隆中華鱉gp130基因cDNA全長。方法：采用RACE法克隆中華鱉gp130基因全長cDNA，并對其進行生物信息學分析。結果：中華鱉gp130基因cDNA全長3806 bp，對應基因組全長73 252 bp，包含16個內含子及17個外顯子。中華鱉gp130與鳥類、哺乳類相比均有保守的共線性關系。該基因編碼由927個氨基酸殘基構成的序列，二級結構主要包括α螺旋、延伸鏈、無規(guī)則卷曲、β轉角，三級結構包含配體結合、跨膜結構域、纖維鏈接蛋白結構域等。系統(tǒng)進化分析顯示與鳥類首先聚類，其次是哺乳類，最后是兩棲類與魚類。同時以人GP130蛋白為模型構建了中華鱉GP130蛋白的3D結構模型，一致性為58.32%。結論：獲得了中華鱉gp130基因全長并做了生物信息學分析，為深入研究GP130及相關信號通路提供了基礎。

中華鱉；糖蛋白130（GP130）；cDNA末端快速擴增（RACE）

JAK-STAT信號通路是近年研究較熱的一條由細胞因子刺激的信號轉導通路，參與細胞的增殖、分化、凋亡及免疫調節(jié)等許多重要的生物學過程[1]。在此信號通路中，糖蛋白130（glycoprotein 130，GP130）［又稱白細胞介素6信號轉導因子（IL-6 signal transducer，IL6st）］是白細胞介素6（interleukin-6，IL-6）、白血病抑制因子（leukemia inhibitory factor，LIF）、抑瘤素 M（oncostatin-M，OSM）、睫狀神經營養(yǎng)因子（ciliary neurotrophic factor，CNTF）、白細胞介素11（IL-11）等因子共用的受體和信號轉導因子[2]。GP130既可以膜型形式整合在機體的細胞膜上，又能以可溶性蛋白的形式分布于體液中，該蛋白在JAK-STAT信號通路的調控中起重要作用[3]。

有關GP130及相關信號通路的研究主要集中于哺乳動物。在腫瘤發(fā)生方面，Greten小組發(fā)現(xiàn)該基因的第757號酪氨酸至苯丙氨酸的點突變（gp130Y757F）會引起STAT3的超活化，這一位點突變的轉基因小鼠在氧化偶氮甲烷/葡聚糖硫酸鈉（AOM/DSS）的誘導下發(fā)生腸道腫瘤的比例明顯升高[4]。GP130蛋白可激活多個下游信號通路，如可磷酸化激活STAT轉錄因子，尤其是STAT3，激活的STAT3二聚體化，轉移至細胞核中，激活下游基因的轉錄；此外，還可引起酪氨酸磷酸酶SHP2聚集并激活下游基因Ras/Erk（胞外信號調節(jié)激酶）和PI3K（磷酸肌醇3激酶）信號通路[5-6]。免疫功能方面，可參與炎癥反應，如對過敏性哮喘患者的外周血檢測發(fā)現(xiàn)GP130異常增高[7]。IL-6/IL-6R/GP130激活STAT3的信號途徑在Th17細胞的分化過程中起重要作用，GP130信號通路在Th17的誘導中是必需的[8]。另外，在小鼠中的最新研究發(fā)現(xiàn)，肝細胞中IL-6-GP130-STAT3信號通路在T細胞介導的肝臟炎癥中起到一定的保護作用[9]。河豚中，IL-6可與IL-6R/GP130復合體結合，從而激活JAK-STAT3信號通路，促進IgM的產生[10]，在斑馬魚的研究中也發(fā)現(xiàn)gp130轉錄在魚體發(fā)育的早期就已產生，且在魚體前后軸線上均能檢測到，分布較廣泛[11]。但在進化上有重要過渡作用的爬行類中卻鮮有gp130的研究報道。中華鱉作為研究較多的爬行類物種，是我國重要的水產養(yǎng)殖品種[12]，廣泛分布于長江流域及華南地區(qū)[13]。故而中華鱉gp130基因的克隆及相關分析，對研究爬行動物免疫尤其是JAK-STAT3信號通路具有重要意義。

1 材料與方法

1.1 材料

實驗用中華鱉購自河北省玉田縣中華鱉良種場。于實驗室馴化2周，期間每日飼喂商品化飼料一次并換水，水溫維持在（25±1）℃。實驗時，斷頭處死中華鱉，快速分離脾臟、肝臟、腸道組織，液氮速凍，保存于-80℃?zhèn)溆谩?/p>

1.2 總RNA提取及反轉錄合成cDNA第一鏈

采用TRIzol法提取總RNA[14]。將-80℃保存的組織樣品取出，加入1 mL TRIzol（Invitrogen公司），組織勻漿，室溫靜置5 min，離心（12 000×g，5 min，4℃）；取上清，加入1/5體積的氯仿，渦旋15 s，室溫靜置5 min，離心（12 000×g，15 min，4℃）；取上清轉移到新的無RNA酶的離心管中，加入等體積的異丙醇，顛倒混勻，室溫靜置10 min，離心（12 000×g，10 min，4℃）得到RNA沉淀。沉淀用75%乙醇洗滌2次，室溫干燥5～10 min，待干燥后，加入DEPC處理的去離子水，55℃溶解10 min，冰浴20 min，分光光度計（Eppendorf公司）測量總RNA的濃度。1%瓊脂糖凝膠電泳檢測RNA的質量。按照SMARTer RACE cDNA擴增試劑盒（TaKaRa公司）操作手冊，將提取的脾臟、肝臟、腸道RNA混合后用于合成3'-與5'-RACE-ready cDNA文庫。

1.3 gp130全長cDNA的克隆

根據(jù)自NCBI數(shù)據(jù)庫中得到的中華鱉gp130基因預測序列設計引物，由上海生工公司合成（表1），按照SMARTer RACE cDNA擴增試劑盒操作說明，進行5'和3'RACE，從中華鱉cDNA文庫中擴增gp130的3'與5'端片段。其中，UPM與NUP為通用引物。按照普通PCR反應進行擴增。PCR產物經0.8%瓊脂糖凝膠電泳分離檢測后，切膠回收相應的目的片段，克隆到pMD19-T simple（TaKaRa公司）載體上，連接產物轉化DH5α大腸桿菌，挑單克隆培養(yǎng)，PCR篩選陽性克隆送上海生工公司測序，將測序結果進行BLAST比對并拼接全長，依此拼接序列設計gp130全長cDNA驗證引物，經PCR及測序驗證后獲得cDNA全長序列。

表1 實驗所用引物序列及用途

1.4 生物信息學分析

用BioEdit 7.2.3軟件對中華鱉gp130cDNA片段進行拼接、分析，并利用在線翻譯軟件Ex-PASy Translate（http://web.expasy.org/translate/）將gp130序列翻譯為氨基酸序列。各個脊椎動物基因結構及共線性關系通過搜索Ensembl數(shù)據(jù)庫（http://asia.ensembl.org/index.htm l）得到。中華鱉和其他脊椎動物的氨基酸序列多重比對采用ClustalX1.8 軟件結合 BoxShade（http://www.ch.embnet. org/software/BOX_form.htm l）在線軟件完成，并在MEGA 6.06軟件中用鄰接法進行進化樹構建（自舉檢驗1000次）。用到的其他物種的氨基酸序列均來自GenBank。采用SOPMA（https://npsa-prabi. ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html）在線軟件預測蛋白質二級結構，蛋白3D結構預測采用同源建模法，用SWISS-MODEL（http:// swissmodel.expasy.org）在線軟件完成。同時，結合SMART在線軟件（http://smart.embl-heidelberg.de/）預測蛋白三級結構。

2 結果

2.1 中華鱉gp130基因全長獲得及序列分析

經過3'及5'RACE，分別獲得中華鱉gp130基因的3'及5'cDNA片段，經全長PCR驗證，獲得了3.5 kb左右的惟一條帶（圖1）。將片段切膠回收，連接到pMD19T simple載體，送上海生工公司測序，將測序結果與NCBI數(shù)據(jù)庫進行BLAST比對確定為中華鱉gp130基因。序列分析表明基因全長為3806 bp，開放讀框長度為2784 bp，編碼927個氨基酸殘基，5'非翻譯區(qū)（UTR）長度為90 bp，3'UTR長度為936 bp（圖2）。分析顯示中華鱉GP130包含WSXWS（色氨酸-絲氨酸-X-色氨酸-絲氨酸）序列、血細胞生成素受體序列，氨基酸序列末端還具有STAT3結合位點等（圖2）。

2.2 中華鱉gp130基因結構分析與共線性關系分析

通過搜索Ensembl數(shù)據(jù)庫中中華鱉的基因組序列，并結合克隆得到的中華鱉gp130基因cDNA序列，我們獲得了該基因的結構示意圖。對應的基因組全長為73 252 bp，包含17個外顯子、16個內含子。圖3為中華鱉gp130基因的基因結構。

根據(jù)Ensembl數(shù)據(jù)庫獲得的基因相對位置信息，我們繪制了中華鱉及其他進化上具有代表性的脊椎動物（人、小鼠、原雞、爪蟾、斑馬魚）的gp130基因的共線性關系圖（圖4）。比較發(fā)現(xiàn)，gp130與鄰近的4～5個基因存在保守的共線性關系，除個別物種外（如斑馬魚），該基因附近均可發(fā)現(xiàn)DDX4、IL31Ra、ANKRD55及MAP3K1等基因；同時與哺乳類與鳥類相比，沒有發(fā)現(xiàn)基因錯位與顛倒現(xiàn)象，但與爪蟾相比，gp130基因的編碼方向發(fā)生了顛倒。

2.3 中華鱉gp130基因編碼的蛋白特征、結構預測及系統(tǒng)進化樹構建

圖5為幾種脊椎動物的GP130氨基酸全序列比對。構建的幾種脊椎動物的系統(tǒng)發(fā)育樹顯示（圖6），中華鱉GP130蛋白與西部錦龜首先聚類，其次是鳥類，然后與哺乳動物聚類，但與兩棲類、魚類等在進化上的距離較遠。采用SOPMA在線軟件預測結果表明，GP130蛋白的二級結構主要由α螺旋（17.80%）、延伸鏈（26.11%）、無規(guī)則卷曲（46.93%）、β轉角（9.17%）組成。用SMART在線蛋白結構分析軟件分析，發(fā)現(xiàn)該蛋白包含跨膜結構域、配體結合區(qū)、數(shù)個纖維鏈接蛋白（FN3）結構域。在SWISS-MODEL網(wǎng)站以人GP130為模型（PDB ID：3l5h.1.A），構建了中華鱉GP130蛋白的3D模型（圖7），兩者的一致性為58.32%。

圖1 中華鱉gp130基因全長cDNA的PCR擴增驗證

圖2 中華鱉gp130基因全長序列

3 討論

人類GP130首先獲得克隆，是一個包含918個氨基酸殘基的跨膜蛋白，相對分子質量為130× 103[15]。在鳥類及魚類中也有相應報道[11，16]。本研究以中華鱉為研究對象，采用SMART RACE技術，獲得了gp130基因的cDNA全長3806 bp，其中開放讀框2784 bp，編碼927個氨基酸殘基。分析顯示該序列含有許多功能結構：WSXWS序列，是造血細胞因子受體家族膜外部分受體與配體結合的關鍵序列[17]；血細胞生成素受體序列，是GP130家族特征結構；氨基酸序列末端還具有STAT3結合位點，該結構可促進與STAT3的結合，同時有助于激活STAT3[18]。在與中華鱉基因組比對時發(fā)現(xiàn)，gp130的基因組共跨越了16個內含子，總長近70 kb。同時分析比較了其他類型基因，發(fā)現(xiàn)比其他基因的內含子要長很多（結果未展示）。具體功能有待進一步探究。

圖2 中華鱉gp130基因全長序列（續(xù)）

圖3 中華鱉gp130基因與基因組DNA結構示意圖

圖4 幾種脊椎動物gp130基因的共線性關系圖

在比較了中華鱉與其他脊椎動物的gp130上下游4～5個基因的共線性關系后發(fā)現(xiàn)，在基因位置關系上，中華鱉與原雞、小鼠、人很相近，并沒有發(fā)現(xiàn)基因倒位或顛倒的現(xiàn)象；但在爪蟾中gp130的編碼方向發(fā)生顛倒，與斑馬魚相比附近基因種類差異較大。這表明，在位置保守性上中華鱉更加接近于鳥類與哺乳類。同樣，在進化樹中也可以看到，中華鱉及錦龜GP130首先與原雞聚類，其次是以小鼠和人為代表的哺乳類聚為一支，最后才是爪蟾與魚類。這也表明，在進化關系上，爬行類的GP130進化關系更為接近鳥類與哺乳類，而非兩棲類。這與傳統(tǒng)觀點一致[19]。

氨基酸序列分析表明我們獲得的GP130不含信號肽，說明該蛋白不是分泌型蛋白。包含的二級結構主要為α螺旋、延伸鏈、無規(guī)則卷曲和β轉角。通過蛋白三級結構分析，該蛋白含有配體結合區(qū)、跨膜結構域及FN3結構域。IL-6與其特異性受體IL-6R結合后使IL-6R發(fā)生構象變化，進而迅速與2分子信號轉導蛋白GP130結合，導致GP130同源二聚體的形成[20]。GP130分子的二聚化使得與之偶聯(lián)的JAK激酶相互接近，并通過交互的酪氨酸磷酸化作用而活化，進而磷酸化STAT3，激活下游通路[5，21]。采用同源建模法，以人GP130為模型預測了中華鱉GP130蛋白的3D結構，一致性為58.32%，兩者的結構較為相似。說明該蛋白較為保守。

綜上，我們克隆了中華鱉gp130基因cDNA全長，分析了其基因及蛋白結構，發(fā)現(xiàn)其在進化地位上更接近于哺乳類與鳥類。這為爬行動物的先天性免疫研究，特別是STAT信號通路相關研究奠定了基礎。

圖6 NJ法構建幾種脊椎動物GP130氨基酸序列進化樹

圖7 中華鱉GP130蛋白3D模型

[1]張克英,唐君,方相春,等.JAK-STAT信號轉導通路在神經系統(tǒng)疾病中作用的研究進展[J].神經解剖學雜志, 2014,30(2):237-240.

[2]Ip N Y,Nye S H,Boulton T G,et al.CNTF and LIF act on neuronal cells via shared signaling pathways that involve the IL-6 signal transducing receptor component gp130[J].Cell,1992,69(7):1121-1132.

[3]Scheller J,Kovaleva M,Rabe B,et al.Development of a monoclonal antibody-based enzyme-linked immunoabsorbent assay for the binding of gp130 to the IL-6/IL-6R complex and its competitive inhibition[J].J Immunol Methods,2004,291(1-2):93-100.

[4]Bollrath J,Phesse T J,von Burstin V A,et al. gp130-mediated Stat3 activation in enterocytes regulates cell survival and cell-cycle progression during colitis-associated tumorigenesis[J].Cancer Cell,2009,15 (2):91-102.

[5]林麗艷,張慧云,何韶衡.IL-6及其受體與炎癥性疾病關系的新進展[J].中國熱帶醫(yī)學,2008,8(4):680-682.

[6]劉徽,朱波,林治華.IL-6信號通路與腫瘤[J].細胞與分子免疫學雜志,2011,27(3):353-355.

[7]張臘娣,黃建安,於葛華.等.過敏性哮喘外周血細胞因子白細胞介素2、白細胞介素6及其受體的測定和臨床意義[J].臨床薈萃,2003,18(8):421-423.

[8]Kaneda M,Odaka T,Suetake H,et al.Teleost IL-6 promotes antibody production through STAT3 signaling via IL-6R and gp130[J].Dev Comp Immunol,2012,38 (2):224-231.

[9]Nishihara M,Ogura H,Ueda N,et al.IL-6-gp130-STAT3 in T cells directs the development of IL-17+Th with a minimum effect on that of Treg in the steady state[J].Int Immunol,2007,19(6):695-702.

[10]Klein C,Wüstefeld T,Assmus U,et al.The IL-6-gp130-STAT3 pathway in hepatocytes triggers liver protection in T cell-mediated liver injury[J].J Clin Invest,2005,115(4):860-869.

[11]Varela M,Dios S,Novoa B,et al.Characterisation,expression and ontogeny of interleukin-6 and its receptors in zebrafish(Danio rerio)[J].Dev Comp Immunol, 2012,37(1):97-106.

[12]Zhou F,Ding X Y,Feng H,et al.The dietary protein requirement of a new Japanese strain of juvenile Chinese soft-shelled turtle, Pelodiscus sinensis[J]. Aquaculture,2013,412-413(6):74-80.

[13]黃麗英,何中央,丁詩華,等.中華鱉種質資源的研究現(xiàn)狀及保護、利用對策[J].寧波大學學報(理工版), 2005,18(2):183-186.

[14]Zhang Z,Niu C,Storset A,et al.Comparison of Aeromonas salmonicida resistant and susceptible salmon families:a high immune response is beneficial for the survival against Aeromonas salmonicida challenge [J].Fish Shellfish Immunol,2011,31(1):1-9.

[15]朱錦芳,鄭仲承.細胞因子受體及其介導的信號轉導[J].生命的化學,1996,(3):14-18.

[16]Schmitz J,Dahmen H,Grimm C,et al.The cytoplasmic tyrosine motifs in full-length glycoprotein 130 have different roles in IL-6 signal transduction[J].J Immunol,2000,164(2):848-854.

[17]Hibi M,Murakami M,Saito M,et al.Molecular cloning and expression of an IL-6 signal transducer, gp130[J].Cell,1990,63(6):1149-1157.

[18]Geissen M,Heller S,Pennica D,et al.The specification of sympathetic neurotransmitter phenotype depends on gp130 cytokine receptor signaling[J].Development,1999,125(23):4791-4801.

[19]鄭光美.為什么近年有些教材將鳥類歸入“爬行綱”[J]?生物學通報,2010,45(11):17-18.

[20]Kishimoto T,Akira S,Narazaki M,et al.Interleukin-6 family of cytokines and gp130[J].Blood,1995,86(4): 1243-1254.

[21]Taga T,Kishimoto T.Gp130 and the interleukin-6 family of cytokines[J].Immunology,1997,15(15):797-819.

Molecular Cloning and Bioinformatics Analysis of gp130 in Chinese Soft-Shelled Turtle（Pelodiscus sinensis）

SONG Ru-Xin,WANG Lan*,ZHANG Zuo-Bing*
School of Life Sciences,Shanxi University,Taiyuan 030006,China

Objective:To clone thegp130full-length cDNA from Chinese soft-shelled turtle.Methods:RACE PCR and bioinformation were used to clone the full-length sequence and predict thegp130gene and protein.Result:The length ofgp130cDNA of Chinese soft-shelled turtle was 3806 bp,encoding a 927 amino acid sequence and its genomic DNA length was 73 252 bp,which spans 16 introns and 17 exons.The BLAST results suggested that thegp130is a homolog of those of birds and mammalians.Bioinformatics analysis showed that the secondary structure consists ofα-helix,extension strand,random coil andβ-turn.The result of tertiary structure analysis showed that the protein contains ligand binding domain,trans membrane domains and FN3 domains.Phylogenetic analysis demonstrated that the reptilegp130clustered firstly with birds,and secondly with mammals,and finally with amphibians and fish.Meantime,we use the human GP130 protein as a model to construct the turtle GP130,the identity was 58.32%.Conclusion:We cloned thegp130full-length cDNA from Chinese soft-shelled turtle.This study laid a basement for further study ofgp130and the related signal pathways.

Pelodiscus sinensis;glycoprotein 130（GP130）;rapid amplification of cDNA ends（RACE）

Q78

A

1009-0002（2017）03-0265-09

10.3969/j.issn.1009-0002.2017.03.006

2016-12-02

國家自然科學基金（31400343）

宋如昕（1992-），男，碩士研究生，（E-mail）srxsong@163.com

王蘭，（E-mail）lanwang@sxu.edu.cn；張左兵，（E-mail）zbzhang@sxu.edu.cn

*Co-corresponding authors,WANG Lan,E-mail:lanwang@sxu.edu.cn;ZHANG Zuo-Bin,E-mail:zbzhang@sxu.edu.cn