段亞飛, 張 喆, 李吉濤, 李 健, 劉 萍

( 1.中國水產(chǎn)科學(xué)研究院 南海水產(chǎn)研究所，農(nóng)業(yè)部南海漁業(yè)資源開發(fā)利用重點實驗室，廣東 廣州 510300；2. 中國水產(chǎn)科學(xué)研究院 黃海水產(chǎn)研究所，農(nóng)業(yè)部海洋漁業(yè)可持續(xù)發(fā)展重點實驗室，山東 青島 266071；3. 青島海洋科學(xué)與技術(shù)國家實驗室，海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實驗室，山東 青島 266071 )

脊尾白蝦血細(xì)胞ESTs的生物信息學(xué)與微衛(wèi)星序列特征分析

段亞飛1, 張 喆1, 李吉濤2, 3, 李 健2, 3, 劉 萍2, 3

( 1.中國水產(chǎn)科學(xué)研究院 南海水產(chǎn)研究所，農(nóng)業(yè)部南海漁業(yè)資源開發(fā)利用重點實驗室，廣東 廣州 510300；2. 中國水產(chǎn)科學(xué)研究院 黃海水產(chǎn)研究所，農(nóng)業(yè)部海洋漁業(yè)可持續(xù)發(fā)展重點實驗室，山東 青島 266071；3. 青島海洋科學(xué)與技術(shù)國家實驗室，海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實驗室，山東 青島 266071 )

對前期測序得到的脊尾白蝦血細(xì)胞2853條EST序列進(jìn)行了生物信息學(xué)和微衛(wèi)星序列特征分析。EST序列拼接得到1053條Unigenes，包括329條Contigs和724條Singlets。BLAST分析表明，593 (56.3%)條Uingenes與數(shù)據(jù)庫中已知基因具有相似性。KEGG代謝途徑分析表明，181條Unigenes映射到120條代謝途徑。通過EST-SSR分析，共得到416條微衛(wèi)星序列，檢出率為14.58%。其中，兩堿基重復(fù)序列374條，占89.90%，AG重復(fù)類型最多；三堿基35條，占8.41%，AAT重復(fù)類型最多；四堿基7條，占1.68%，AAGT重復(fù)類型最多。本研究可為脊尾白蝦功能基因資源挖掘及分子標(biāo)記篩選提供有效數(shù)據(jù)。

脊尾白蝦；血細(xì)胞；EST；生物信息學(xué)；微衛(wèi)星

脊尾白蝦(Exopalaemoncarinicauda)是我國黃、渤海地區(qū)重要的經(jīng)濟蝦類，具有繁殖能力強、生長速度快和環(huán)境適應(yīng)性強等優(yōu)點[1-2]。近年來，脊尾白蝦養(yǎng)殖產(chǎn)業(yè)發(fā)展迅速，其遺傳育種研究日益受到重視。脊尾白蝦基因信息和分子遺傳標(biāo)記的發(fā)掘，可以為其功能基因篩選、遺傳多樣性分析、遺傳連鎖圖譜構(gòu)建及標(biāo)記輔助育種等奠定基礎(chǔ)。

表達(dá)序列標(biāo)簽(EST)是指從特定組織來源的cDNA文庫中隨機挑取單克隆進(jìn)行單向測序而得到的部分cDNA序列，代表了生物體特定組織在某一時期的表達(dá)基因[3]。表達(dá)序列標(biāo)簽分析可以用于新基因發(fā)現(xiàn)、功能基因篩選、分子標(biāo)記開發(fā)和基因時空表達(dá)狀況分析等[3]。微衛(wèi)星標(biāo)記(SSR)是水產(chǎn)動物遺傳育種重要的分子標(biāo)記，具有多態(tài)性高、共顯性、通用性好、數(shù)量豐富和在基因組上分布均勻等特點[4-5]。EST-SSR可以避免傳統(tǒng)方法構(gòu)建、篩選文庫的繁瑣步驟和人力物力的耗費，廣泛應(yīng)用于水生動物新基因的分離鑒定、比較基因組學(xué)、遺傳圖譜構(gòu)建和基因芯片制備等研究[6-7]。

近年來，許多甲殼動物已獲得EST-SSRs標(biāo)記，如中國明對蝦(Fenneropenaeuschinensis)[8]、凡納濱對蝦(Litopenaeusvannamei)[9]、斑節(jié)對蝦(Penaeusmonodon)[10]、三疣梭子蟹(Portunustrituberculatus)[11]和擬穴青蟹(Scyllaparamamosain)[12]等。而目前脊尾白蝦微衛(wèi)星標(biāo)記僅見通過構(gòu)建基因組微衛(wèi)星富集文庫篩選獲得[13]。本研究利用前期構(gòu)建的脊尾白蝦血細(xì)胞全長cDNA文庫，對測序獲得的EST序列進(jìn)行生物信息學(xué)和微衛(wèi)星序列特征分析，以期為脊尾白蝦功能基因資源挖掘及分子標(biāo)記篩選提供有效數(shù)據(jù)。

1 材料與方法

1.1脊尾白蝦血細(xì)胞EST序列測定和分析

利用前期構(gòu)建的脊尾白蝦血細(xì)胞全長cDNA文庫[3]，對其進(jìn)行批量測序，引物為5′pTriPLEx2 (5′-CTCCGAGATCTGGACGAGC-3′)。測得的EST序列，利用Cross-match軟件進(jìn)行引物序列、載體序列及低質(zhì)量序列(<100 bp)去除；利用Phrap軟件對高質(zhì)量ESTs進(jìn)行序列拼接，拼接結(jié)果為由Contigs和Singlets組成的Unigenes。

1.2 ESTs的相似性比對、分類及代謝途徑分析

利用BLAST軟件對拼接到的Unigenes與NCBI核苷酸和蛋白質(zhì)數(shù)據(jù)庫進(jìn)行相似性比對(E-value < 1e-5)。將所有的Unigenes與日本京都大學(xué)的《京都基因和基因組百科全書》(KEGG)的PATHWAY子數(shù)據(jù)庫進(jìn)行比較和代謝途徑分析，獲得脊尾白蝦血細(xì)胞的ESTs的部分代謝網(wǎng)絡(luò)信息。

1.3 ESTs的微衛(wèi)星序列特征分析

利用RepeatMasker軟件從拼接的Unigenes中查找微衛(wèi)星序列，最小重復(fù)參數(shù)設(shè)置為重復(fù)7次以上的兩堿基序列，重復(fù)5次以上的三堿基序列，重復(fù)5次以上的四堿基序列。單核苷酸重復(fù)序列未進(jìn)行篩選，因為一般認(rèn)為其作為多態(tài)性標(biāo)記的意義不大[18]。

2 結(jié) 果

2.1 ESTs的生物信息學(xué)分析

2.1.1 ESTs序列分析

對脊尾白蝦血細(xì)胞全長cDNA文庫批量測序獲得的EST序列，去除載體序列、引物序列和低質(zhì)量載體序列(<100 bp)后，共得到2853條高質(zhì)量EST序列，平均長度為493 bp(表1)。利用Phrap軟件進(jìn)行序列拼接，共獲得1053條Unigenes，包括329條(31.24%)Contigs和724條(68.76%)Singlets。對Unigenes中的ESTs數(shù)量分析可知，Contigs由2～284條ESTs拼接而成，其中大部分是由2～5條ESTs拼接而成，10條Contigs是由多于50條的ESTs拼接而成，只有1條Contigs是由284條ESTs拼接而成(圖1)。

表1 脊尾白蝦血細(xì)胞ESTs基本信息

BLAST分析表明，593(56.3%)條Uingenes與數(shù)據(jù)庫中已知基因具有相似性，22(2.1%)條Uingenes顯示為未知基因，438(41.6%)條Uingenes沒有任何相似性。ESTs相似比對的E-value分布表明，276(44.9%)條Uingenes具有高度相似性(E-value < 1e-100)，256(41.6%)條Uingenes具有中度相似性(1e-99

圖1 脊尾白蝦血細(xì)胞ESTs與其所拼接形成的Unigene之間的關(guān)系

圖2 脊尾白蝦血細(xì)胞ESTs相似比對的E-value分布 圖3 脊尾白蝦血細(xì)胞ESTs匹配的相似物種分布

2.1.2 ESTs的代謝途徑分析

脊尾白蝦血細(xì)胞ESTs的代謝途徑分析表明，共有181條Unigenes得到注釋，相關(guān)酶類參與120條代謝途徑，主要包括代謝、細(xì)胞過程、環(huán)境信息的加工、遺傳信息的加工和人類疾病5類，其相關(guān)的基因序列分別占33.1%、21.6%、12.2%、16.0%和17.1%。其中，代謝途徑最多的類別是氨基酸代謝(11%)，其次為傳染病(10%)，第三為信號轉(zhuǎn)導(dǎo)(9%)(表2)。

表2 脊尾白蝦血細(xì)胞ESTs的代謝途徑分析

2.2 ESTs的微衛(wèi)星序列分析

從血細(xì)胞2853條ESTs序列中，共篩選到微衛(wèi)星序列416個，占全部ESTs總數(shù)的14.58%?？紤]到堿基互補配對和計數(shù)拷貝數(shù)起始堿基順序的差異，本研究將同類重復(fù)兼并為一種重復(fù)類型代表，如AAG代表六種重復(fù)類型，分別為AAG、AGA、GAA、TTC、TCT和CTT。統(tǒng)計所有重復(fù)類型，以兩堿基重復(fù)序列數(shù)量最多，為374條，占重復(fù)序列總數(shù)的89.90%；其次為三堿基，35個，占8.41%；再次是四堿基，7個，占1.68%(表3)。

表3 不同重復(fù)類型的重復(fù)序列數(shù)目及其百分比

同類型的SSR重復(fù)序列中，各重復(fù)拷貝類別所占的比例也各不相同。在兩堿基重復(fù)類型中，AG重復(fù)拷貝類別最多，為355個，占94.92%；其次為AC和AT，分別為12(3.21%)和7(1.87%)；未發(fā)現(xiàn)GC重復(fù)拷貝類別。三堿基重復(fù)類型中，共有6種重復(fù)拷貝類別，分別為AAT、AAG、ATC、ACT、AGC和AGG。其中，以AAT最多，為21個，占60%；其次為ACT，為5個，占14.29%；再次為AAG和ATC，這兩種重復(fù)拷貝類別數(shù)量相同，為3個，占8.57%。四堿基重復(fù)類型中，共有2種重復(fù)拷貝類別，分別為AAGT和GCCG，其中以AAGT最多，為7個，占85.71%；其次為GCCG，為1個，占14.29%(表4)。

表4脊尾白蝦血細(xì)胞ESTs篩選到的微衛(wèi)星序列

重復(fù)類型重復(fù)單元數(shù)量比例/%二堿基AT71.87AG35594.92AC123.21三堿基AAT2160.00AAG38.57ATC38.57ACT514.29AGC12.86AGG25.71四堿基AAGT685.71GCCG114.29

3 討 論

3.1 ESTs生物信息學(xué)分析

通過對前期測序得到的脊尾白蝦血細(xì)胞2853條ESTs進(jìn)行序列拼接，共得到1053條Uingenes。BLAST分析表明，41.6%的Uingenes與美國國立生物技術(shù)信息中心數(shù)據(jù)庫中的已知基因序列無相似性，與三疣梭子蟹[14]、中華絨螯蟹(Eriocheirsinensis)[15]等的研究一致。此類無相似性的ESTs可能為新基因序列，表明甲殼類尚有大量基因未被發(fā)掘；又可能為已知基因，由于蝦類與其他脊椎動物的親緣關(guān)系太遠(yuǎn)，而無法通過BLAST比對確定其相似關(guān)系。此外，Xu等[14]認(rèn)為，ESTs是通過一次測序得到的序列，未必包含編碼區(qū)，可能是ESTs未能注釋的原因之一。因此，此類未得到注釋的ESTs，有待進(jìn)一步的深入研究。

KEGG數(shù)據(jù)庫涵蓋了物種的基因與基因組、酶促反應(yīng)及其途徑和各種生化物質(zhì)的信息，在基因與基因組、代謝途徑與代謝網(wǎng)絡(luò)的相關(guān)研究中具有重要作用[16]。KEGG PATHWAY數(shù)據(jù)庫不僅能夠提供生化物質(zhì)相互轉(zhuǎn)化所有可能的代謝途徑，還可以對催化各步反應(yīng)的酶進(jìn)行全面注解[16]。通過代謝途徑分析，本研究共有181條Unigenes(17.2%)獲得注釋，發(fā)現(xiàn)酶類參與120條代謝途徑。代謝途徑分析在多個水產(chǎn)物種中已有相關(guān)報道。Meng等[17]對蝦夷扇貝(Patinopectenyessoensis)肌肉EST序列進(jìn)行KEGG代謝途徑分析，共有11.3%的Unigenes得到注釋，涉及103個代謝途徑；其中以代謝方面的Unigenes比例最大，占30.5%，與本研究結(jié)果一致。

3.2 ESTs的微衛(wèi)星序列特征分析

本研究通過對脊尾白蝦血細(xì)胞ESTs進(jìn)行SSR篩選，共得到416條(14.58%)SSR序列，表明脊尾白蝦ESTs-SSR序列比例較高。研究表明，動物的ESTs-SSR以兩堿基重復(fù)為主，大部分植物則以三堿基重復(fù)為主[18]。脊尾白蝦血細(xì)胞ESTs-SSR共發(fā)現(xiàn)3種堿基重復(fù)類型，以兩堿基重復(fù)類型最多，占89.90%；其次為三堿基重復(fù)類型(8.41%)和四堿基重復(fù)類型(1.68%)。本研究未發(fā)現(xiàn)五堿基和六堿基重復(fù)類型，可能與所分析的ESTs序列的組織差異有關(guān)。

兩堿基重復(fù)類型中，共發(fā)現(xiàn)3種重復(fù)拷貝類別，其中以AG最多，其次是AC，與斑節(jié)對蝦[19]、羅氏沼蝦(Macrobrachiumrosenbergii)[20]和三疣梭子蟹[21]等一致。賈舒雯等[13]通過構(gòu)建脊尾白蝦微衛(wèi)星富集文庫對多態(tài)性標(biāo)記進(jìn)行了篩選，發(fā)現(xiàn)脊尾白蝦基因組微衛(wèi)星核心序列以AG重復(fù)最多，占88.94%；其次是AC重復(fù)，占5.03%。這與本研究兩堿基重復(fù)類型結(jié)果(AG 94.92%，AC 3.21%)一致，由此說明AG重復(fù)拷貝類別在脊尾白蝦基因組中含量豐富。但本研究并未發(fā)現(xiàn)GC重復(fù)，而GC兩堿基重復(fù)在其他物種基因組中的含量均較少[22-23]。

三堿基重復(fù)類型中，共發(fā)現(xiàn)6種重復(fù)拷貝類別，其中以AAT最多，占60%。研究表明，中國明對蝦[18]和紅鰭東方鲀(Takifugurubripes)[24]中，以AAT重復(fù)最多。高煥等[18]通過分析中國明對蝦基因組微衛(wèi)星特征，發(fā)現(xiàn)三堿基重復(fù)類型中以AAT重復(fù)最多，占50.34%。四堿基重復(fù)類型中，共發(fā)現(xiàn)7種重復(fù)拷貝類別，總體數(shù)量較少。而本研究中未發(fā)現(xiàn)五堿基和六堿基重復(fù)類型，可能與所分析的ESTs序列相對較少、覆蓋面較小有關(guān)。

目前脊尾白蝦SSR的研究，主要是利用SSR核心序列進(jìn)行探針雜交或以核心重復(fù)序列為引物進(jìn)行PCR篩選，而從ESTs中篩選SSR序列的方法則更為簡單，不僅能夠節(jié)省人力和物力，而且由于ESTs-SSR標(biāo)記位于基因編碼區(qū)中，因此具有高度的物種間轉(zhuǎn)移性[25]。本研究所篩選獲得的SSR序列，可以用于脊尾白蝦SSR標(biāo)記篩選，對脊尾白蝦遺傳圖譜、基因定位和遺傳多樣性分析等奠定基礎(chǔ)。

[1] 梁俊平, 李健, 劉萍, 等. 脊尾白蝦生物學(xué)特性與人工繁育的研究進(jìn)展[J]. 中國農(nóng)學(xué)通報, 2012,28(17):109-116.

[2] 梁俊平, 李健, 李吉濤, 等. 氨氮對脊尾白蝦幼蝦和成蝦的毒性試驗[J]. 水產(chǎn)科學(xué), 2012,31(9):526-529.

[3] 段亞飛, 劉萍, 李吉濤, 等. 脊尾白蝦血細(xì)胞全長cDNA文庫的構(gòu)建及EST序列分析[J]. 中國水產(chǎn)科學(xué), 2013, 20(2):243-249.

[4] 孫立元, 郭華陽, 朱彩艷, 等. 卵形鯧鲹育種群體遺傳多樣性分析[J]. 南方水產(chǎn)科學(xué), 2014,10(2):67-71.

[5] 楊兵, 林琳, 李純厚, 等. 基于高通量測序的二長棘鯛微衛(wèi)星標(biāo)記開發(fā)與評價[J]. 南方水產(chǎn)科學(xué), 2015,11(4):116-120.

[6] Ji X S, Chen S L, Ma H Y, et al. Isolation and characterization of 19 EST-linked ploymorphic microsatellite loci for olive flounder (Paralichthysolivaceus)[J]. Aquaculture Research, 2009,40(8): 980-983.

[7] Yu H, Li Q. Exploiting EST databases for the development and characterization of EST-SSRs in the Pacific oyster (Crassostreagigas)[J]. J Heredity,2008,99(2):208-214.

[8] Wang H X, Li F H, Xiang J H. Polymorphic EST-SSR markers and their mode of inheritance inFenneropenaeuschinensis[J]. Aquaculture, 2005,249(1/4):107-114.

[9] Perez F, Ortiz J,Zhinaula M, et al. Development of EST-SSR markers by data mining in three species of shrimp:Litopenaeusvannamei,Litopenaeusstylirostris, andTrachypenaeusbirdy[J]. Marine Biotechnol, 2005,7(5):554-569.

[10] Maneeruttanarungroj C, Pongsomboon S, Wuthisuthimethavee S, et al. Development of polymorphic expressed sequence tag-derived microsatellites for the extension of the genetic linkage map of the black tiger shrimp (Penaeusmonodon) [J]. Animal Genetics, 2006,37(4):363-368.

[11] 呂建建, 王渝, 高保全, 等. 三疣梭子蟹Ⅰ型微衛(wèi)星標(biāo)記的發(fā)掘及多態(tài)性分析[J]. 水產(chǎn)學(xué)報, 2013,37(6):816-822.

[12] Zhou Z H, Zhang Z P, Wang Y L, et al. EST analysis on the gonad development related organs and microarray screen for differentially expressed genes in mature ovary and testis ofScyllaparamamosain[J]. Comp Biochem and Physiol D, 2011,6(2):150-157.

[13] 賈舒雯, 劉萍, 韓智科, 等.脊尾白蝦微衛(wèi)星富集文庫的構(gòu)建與多態(tài)性標(biāo)記的篩選[J]. 水產(chǎn)學(xué)報, 2011,35(12):1787-1794.

[14] Xu Q H, Liu Y, Liu R L. Expressed sequence tags from cDNA library prepared from gills of the swimming crab,Portunustrituberculatus[J]. J Experimental Marine Biol Ecol, 2010,394(1/2):105-115.

[15] Gai Y C, Wang L L, Zhao J M, et al. The construction of a cDNA library enriched for immune genes and the analysis of 7535 ESTs from Chinese mitten crabEriocheirsinensis[J]. Fish and Shellfish Immunol, 2009,27(6):684-694.

[16] Priebe S, Linde J, Albrecht D, et al. FungiFun:a web-based application for functional categorization of fungal genes and proteins [J]. Fungal Genetics Biol, 2011,48(4):353-358.

[17] Meng X Y, Chang Y Q, Qiu X M, et al. Generation and analysis of expressed sequence tags from adductor muscle of Japanese scallopMizuhopectenyessoensis[J]. Comp Biochem Physiol D, 2010,5(4):288-294.

[18] 高煥, 劉萍, 孟憲紅, 等. 中國對蝦基因組微衛(wèi)星特征分析[J]. 海洋與湖沼, 2004,35(5):424-431.

[19] Xu Z,Dhar A K, Wyrzykowski J, et al. Identification of abundant and informative microsatellites from shrimp (Penaeusmonodon) genome [J]. Animal Genetics, 1999,30(2):150-156．

[20] Charoentawee K, Poompuang S, Nakom U. Isolation and characterization of microsatellites in giant freshwater prawnMacrobrachiumrosenbergii[J]. Mol Ecol Resources, 2006,6(3):823-825.

[21] 宋來鵬， 劉萍, 李健, 等. 三疣梭子蟹基因組微衛(wèi)星特征分析[J]. 中國水產(chǎn)科學(xué), 2008,15(5):738-743.

[22] Zhang L,Bao Z, Cheng J, et al. Fosmid library construction and initial analysis of end sequences in Zhikong scallop (Chlamysfarreri) [J]. Marine Biotechnol, 2007,9(5):606-612.

[23] Wang Y, Ren R, Yu Z, et al.Bioinformatic mining of EST-SSR loci in the Pacific oyster,Crassostreagigas[J]. Animal Genetics, 2008,39(3):287-289.

[24] 崔建洲, 申雪艷, 楊官品, 等. 紅鰭東方鲀基因組微衛(wèi)星特征分析[J]. 中國海洋大學(xué)學(xué)報：自然科學(xué)版, 2006,36(2):249-254.

[25] 吳文婷, 張磊, 宋志民, 等. 5個紫菜屬物種絲狀體的微衛(wèi)星初步遺傳分析[J]. 南方水產(chǎn)科學(xué), 2012,8(4):29-36.

BioinformaticsandMicrosatelliteSequencesAnalysisofESTSequenceinRidgeTailShrimpExopalaemoncarinicauda

DUAN Yafei1, ZHANG Zhe1, LI Jitao2, 3, LI Jian2, 3, LIU Ping2, 3

( 1. Key Laboratory of South China Sea Fishery Resources Exploitation & Utilization, Ministry of Agriculture, South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou 510300, China; 2. Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, China; 3. Function Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China )

The analysis of bioinformatics and microsatellite sequences characteristics of expressed sequence tag (EST) sequence were studied, basing on the hemocyte cDNA library of ridge tail shrimpExopalaemoncarinicaudaconstructed by our laboratory. The 1053 unique sequences were yielded from 2853 high quality ESTs, including 329 contigs and 724 singletons. BLAST analysis revealed 593 (56.3%) of the unique sequences as orthologs of genes from other organisms, 22 (2.1%) of which were hypothetical protein, and the remaining 438 (41.6%) had no homology to any sequences. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that 181 of the unique sequences mapped to 120 of metabolic pathways. A total of 416 (14.58%) of the microsatellite repeat sequences were found using EST-SSR analysis. In the 416 repeat sequences, most was dinucleotide with repeats of 374 (89.90%), followed by trinucleotide with repeats of 35 (8.41%), and the third tetranucleotide with repeats of 7 (1.68%). The most common repeat types were AG, AAT and AAGT repeat. The findings provide important information for screening the functional gene resource and molecular markers in ridge tail shrimp.

Exopalaemoncarinicauda; hemocyte; EST; bioinformatics ; microsatellite

10.16378/j.cnki.1003-1111.2016.05.018

S917

A

1003-1111(2016)05-0562-06

2015-11-03；

2015-12-30.

國家蝦產(chǎn)業(yè)技術(shù)體系項目(CARS-47)；國家自然科學(xué)基金資助項目(31472275)；山東省泰山產(chǎn)業(yè)領(lǐng)軍人才工程項目(LNJY2015002)；青島海洋科學(xué)與技術(shù)國家實驗室鰲山科技創(chuàng)新計劃項目(2015ASKJ02).

段亞飛(1989—)，男，助理研究員,碩士；研究方向: 蝦類免疫學(xué). E-mail: duanyafei89@163.com. 通訊作者：劉萍(1962—)，女，研究員；研究方向: 海水養(yǎng)殖生物種質(zhì)資源與遺傳育種. E-mail: liuping@ysfri.ac.cn.