宋如昕，王蘭，張左兵

山西大學(xué) 生命科學(xué)學(xué)院，山西 太原 030006

中華鱉gp130基因全長(zhǎng)cDNA克隆及生物信息學(xué)分析

宋如昕，王蘭，張左兵

山西大學(xué) 生命科學(xué)學(xué)院，山西 太原 030006

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

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

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

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

1 材料與方法

1.1 材料

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

1.2 總RNA提取及反轉(zhuǎn)錄合成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℃）；取上清轉(zhuǎn)移到新的無RNA酶的離心管中，加入等體積的異丙醇，顛倒混勻，室溫靜置10 min，離心（12 000×g，10 min，4℃）得到RNA沉淀。沉淀用75%乙醇洗滌2次，室溫干燥5～10 min，待干燥后，加入DEPC處理的去離子水，55℃溶解10 min，冰浴20 min，分光光度計(jì)（Eppendorf公司）測(cè)量總RNA的濃度。1%瓊脂糖凝膠電泳檢測(cè)RNA的質(zhì)量。按照SMARTer RACE cDNA擴(kuò)增試劑盒（TaKaRa公司）操作手冊(cè)，將提取的脾臟、肝臟、腸道RNA混合后用于合成3'-與5'-RACE-ready cDNA文庫。

1.3 gp130全長(zhǎng)cDNA的克隆

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

表1 實(shí)驗(yàn)所用引物序列及用途

1.4 生物信息學(xué)分析

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

2 結(jié)果

2.1 中華鱉gp130基因全長(zhǎng)獲得及序列分析

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

2.2 中華鱉gp130基因結(jié)構(gòu)分析與共線性關(guān)系分析

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

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

2.3 中華鱉gp130基因編碼的蛋白特征、結(jié)構(gòu)預(yù)測(cè)及系統(tǒng)進(jìn)化樹構(gòu)建

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

圖1 中華鱉gp130基因全長(zhǎng)cDNA的PCR擴(kuò)增驗(yàn)證

圖2 中華鱉gp130基因全長(zhǎng)序列

3 討論

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

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

圖3 中華鱉gp130基因與基因組DNA結(jié)構(gòu)示意圖

圖4 幾種脊椎動(dòng)物gp130基因的共線性關(guān)系圖

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

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

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

圖6 NJ法構(gòu)建幾種脊椎動(dòng)物GP130氨基酸序列進(jìn)化樹

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

[1]張克英,唐君,方相春,等.JAK-STAT信號(hào)轉(zhuǎn)導(dǎo)通路在神經(jīng)系統(tǒng)疾病中作用的研究進(jìn)展[J].神經(jīng)解剖學(xué)雜志, 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及其受體與炎癥性疾病關(guān)系的新進(jìn)展[J].中國(guó)熱帶醫(yī)學(xué),2008,8(4):680-682.

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

[7]張臘娣,黃建安,於葛華.等.過敏性哮喘外周血細(xì)胞因子白細(xì)胞介素2、白細(xì)胞介素6及其受體的測(cè)定和臨床意義[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]黃麗英,何中央,丁詩華,等.中華鱉種質(zhì)資源的研究現(xiàn)狀及保護(hù)、利用對(duì)策[J].寧波大學(xué)學(xué)報(bào)(理工版), 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]朱錦芳,鄭仲承.細(xì)胞因子受體及其介導(dǎo)的信號(hào)轉(zhuǎn)導(dǎo)[J].生命的化學(xué),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]?生物學(xué)通報(bào),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

國(guó)家自然科學(xué)基金（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