陳錦玲 徐媛 陳玉梅 李璐璐 李惠敏 秦新民

摘要 [目的] 研究沙田柚類蛋白激酶基因編碼的蛋白質(zhì)序列所包含的生物信息學。[方法] 利用高通量測序?qū)ι程镨肿越缓彤惤换ㄖM行轉(zhuǎn)錄組測序，通過差異分析得到沙田柚類蛋白激酶基因序列。采用生物信息學方法，對其編碼的蛋白質(zhì)從序列特征、理化性質(zhì)、跨膜結(jié)構(gòu)域、高級結(jié)構(gòu)以及功能域等方面進行預測和分析。[結(jié)果] 該基因全長為2 323 bp，開放閱讀框（ORF）全長為1 560 bp（Genbank登錄號：MG925820），共編碼519個氨基酸，分子質(zhì)量為57.39 kD，等電點PI為8.55。該蛋白質(zhì)含有1個植物STKc IRAK蛋白家族的保守結(jié)構(gòu)域，為親水性穩(wěn)定蛋白。氨基酸序列分析表明，其編碼的氨基酸與克萊門柚和甜橙的同源性分別為 100%、99%。系統(tǒng)進化樹表明沙田柚類蛋白激酶基因與克萊門柚和甜橙親緣關(guān)系很近，屬于同一進化分支。[結(jié)論] 該研究結(jié)果可為今后深入研究沙田柚自交不親和機理提供參考。

關(guān)鍵詞 沙田柚；類蛋白激酶基因；生物信息學分析

中圖分類號 S188文獻標識碼 A文章編號 0517-6611（2018）16-0088-05

Abstract [Objective]To study the biological information contained in the protein sequence encoded receptorlike protein kinases gene.[Method]The transcriptome of the selfpollinated style and crosspollinated style of Citrus grandis var.Shatinyu Hort were sequenced by highthroughput sequencing technology.Receptorlike protein kinases gene sequence of C.grandis var.Shatinyu Hort was obtained though differential analysis method， and some characters of the receptorlike protein kinases protein were analyzed and predicted by bioinformatics method， including the composition of amino acid sequence， physicochemical parameters， hydrophobicity， transmembrane domain， secondary structure and function of protein，etc. [Result] Receptorlike protein kinases gene was 2 323 bp in length with an open reading frame （ORF） of 1 560 bp， encoding 519 amino acids with deduced molecular weight of 57.39 kD， and theoretical pI value of 8.55， and contained a STKc IRAK domain.Bioinformatics analysis showed that receptorlike protein kinases protein was a hydrophobic and stable protein.The homology analysis of amino acid sequence indicated that the receptorlike protein kinases protein shared high homology with that of Citrus clementina （100%） and Citrus sinensis （99%）.Phylogentic analysis revealed that receptorlike protein kinases gene showed closer kinship with that of C.clementina and C.sinensis，indicating that they belong to the same evolutionary branch.[Conclution]The results could provide a theoretical reference for further study of selfincompatibility molecular mechanism in C.grandis var.Shatinyu Hort.

Key words Citrus grandis var.Shatinyu Hort；Receptorlike protein kinases gene；Bioinformatics

沙田柚（Citrus grandis var.Shatinyu Hort）隸屬于蕓香科柑橘屬，具有配子體自交不親和的特性。植物自交不親和（selfincompatibility， SI）是指具有完全花且能夠產(chǎn)生雌雄配子，雌蕊的柱頭或花柱能特異性識別與自體基因型相同的花粉，導致花粉不萌發(fā)或花粉管停止生長的一種現(xiàn)象[1]。這種現(xiàn)象在避免植物近親繁殖方面具有重要的意義。目前，關(guān)于蕓香科植物的自交不親和性機理尚不明確，但在沙田柚自交不親和的形態(tài)學、細胞學和蛋白質(zhì)化學等方面的研究已取得了一定進展[2]。秦新民等[3]分離并鑒定了沙田柚花粉管特異蛋白；楊繼華等[4-5]利用IEF/SDS-PAGESHUA雙向電泳的方法測定了沙田柚自交花柱特異蛋白的分子量、等電點以及N-末端氨基酸序列；薛妙男等[6-7]確定了沙田柚自交花粉管在花柱中生長的受阻部位；秦新民等[8-9]還確定了花柱通道細胞中特異蛋白和花粉管中特異蛋白的產(chǎn)生部位及分布。

植物中的類蛋白激酶通過磷酸化和去磷酸化作用調(diào)節(jié)植物的代謝，這些過程幾乎涉及所有的生理和病理過程[10]，主要作用于植物的細胞抗逆反應[11] 、生長發(fā)育[12-13] 、抗病防衛(wèi)[14-15] 、自交不親和[16]以及胞外信號轉(zhuǎn)導[17]等方面。筆者通過對沙田柚自交和異交花柱進行轉(zhuǎn)錄組測序，獲得了沙田柚類蛋白激酶基因，并對該基因編碼的蛋白質(zhì)理化特征和基因表達進行了分析，旨在為深入研究沙田柚自交不親和的分子機理提供分子理論依據(jù)。

1 材料與方法

1.1 材料

所用材料為廣西靈川縣潮田鄉(xiāng)大山口村果園內(nèi)10年生的沙田柚果樹花柱，分別采集使用人工自交（沙田柚×沙田柚）和異交（酸柚×沙田柚）方式授粉后1～3 d的花柱，立即放置于液氮中，隨后存放到超低溫冰箱（-80 ℃）中備用。

1.2 方法

1.2.1 RNA的提取。根據(jù)改良的 Trizol 法對植物總 RNA 進行提取[18]。

1.2.2 測序。通過檢測選擇合格RNA送至深圳華大基因科技服務有限公司進行建庫和測序。 用 Illumina HiSeqTM 2000 對制備好的文庫進行測序。

1.3 序列分析和系統(tǒng)樹構(gòu)建

將得到的測序結(jié)果序列在NCBI的Blast程序上進行基因序列的同源性比較分析；并利用DNAman、Finder、TMpred、Swissmodel、NetPhos2.0等軟件進行序列和蛋白質(zhì)理化分析；在GenBank上下載與該序列同源性相近的其他物種的序列，用DNAman軟件對類蛋白激酶基因編碼的氨基酸序列構(gòu)建系統(tǒng)進化樹。

2 結(jié)果與分析

2.1 基因的生物信息學分析

沙田柚類蛋白激酶基因（Unigene10208_All）序列全長2 323 bp（GenBank登錄號：MG925820），在NCBI在線數(shù)據(jù)庫（https：//blast.ncbi.nlm.nih.gov/Blast.cgi）中進行比對，尋找該基因的開放閱讀框（ORF），用DNAman軟件對該基因進行翻譯，結(jié)果顯示：該序列有1個1 560 bp的開放閱讀框，編碼的蛋白質(zhì)含519個氨基酸（圖1）。

2.2 編碼蛋白質(zhì)的分析和疏水性預測

將類蛋白激酶基因編碼的蛋白質(zhì)序列放入理化性質(zhì)在線預測網(wǎng)站（http：//web.expasy.org/protparam/）中進行蛋白質(zhì)理化性質(zhì)的預測分析，結(jié)果顯示：該基因編碼的蛋白質(zhì)分子式為C2569H4044N698O745S24，分子質(zhì)量為57.39 kD，等電點PI為8.55。構(gòu)成該蛋白質(zhì)的氨基酸中，亮氨酸（Leu）所占比例最大，達13.7%，組氨酸（His）所占比例最少，僅為1.5%。蛋白質(zhì)不穩(wěn)定系數(shù)（instability index）為33.94（<40），表明該蛋白質(zhì)為穩(wěn)定蛋白。其中，該蛋白質(zhì)攜帶的負電荷氨基酸（Asp + Glu）總數(shù)為 52，正電荷氨基酸（Arg + Lys）總數(shù)為58。利用DNAman軟件對該蛋白質(zhì)進行疏水性分析，結(jié)果見圖2。從圖2可以看出，該基因編碼的肽鏈中疏水性最大值約為3.48，位于第141位氨基酸，最小值約為-3.71，位于第172位氨基酸。該蛋白質(zhì)疏水性平均值為-0.12，鑒定該蛋白質(zhì)為親水性蛋白。

2.3 跨膜預測 運用跨膜蛋白數(shù)據(jù)庫 TMHMM（http：//www.cbs.dtu.dk/services/TMHMM/）對沙田柚類蛋白激酶跨膜區(qū)域進行預測，結(jié)果表明該蛋白質(zhì)屬于跨膜蛋白，其中1～141位部分位于膜外，142～164位部分跨膜，165～519位部分位于膜內(nèi)（圖3）。

2.4 蛋白質(zhì)磷酸化位點預測

蛋白質(zhì)磷酸化位點采用在線預測軟件NetPhos 3.1 Server（http：//www.cbs.dtu.dk/services/NetPhos/）進行預測，結(jié)果顯示，類蛋白激酶基因編碼的蛋白質(zhì)共有39個可能的磷酸化位點，其中絲氨酸（Ser）磷酸化位點共有26個，分別位于蛋白質(zhì)的第4、23、53、77、85、94、107、135、183、191、200、220、244、297、298、303、370、371、374、375、398、442、443、446、483、493位；蘇氨酸（Thr）磷酸化位點共有9個，分別位于蛋白質(zhì)的第68、103、230、232、234、235、253、300和415位；酪氨酸（Tyr）磷酸化位點僅4個，分別位于蛋白質(zhì)的第236、386、391、404位。

可知此蛋白質(zhì)的磷酸化以絲氨酸磷酸化為主，其次為蘇氨酸和酪氨酸。

2.5 功能結(jié)構(gòu)域分析

將編碼的氨基酸序列用NCBI上的Structure/cdd/wrpsb.cgi）進行分析，發(fā)現(xiàn)該蛋白質(zhì)具有1個與STKc IRAK蛋白相同的保守結(jié)構(gòu)域，屬于PKc like超家族，預測結(jié)果如圖4所示。

2.6 蛋白質(zhì)二級以及三級結(jié)構(gòu)的預測

運用在線分析軟件SOPMA（https：//npsa-prabi.ibcp.fr/cgi-bin/secpred_sopma.pl）對類蛋白酶蛋白二級結(jié)構(gòu)進行預測，結(jié)果表明，在該蛋白質(zhì)的二級結(jié)構(gòu)中， α-螺旋所占比例為38.54%，延伸鏈13.29%，無規(guī)則卷曲48.17%，二級結(jié)構(gòu)預測結(jié)果如圖5所示。

利用在線軟件Phyre2 [19]對RING Finger蛋白氨基酸序列三級結(jié)構(gòu)進行預測，結(jié)構(gòu)見圖6。

2.7 同源性分析

從NCBI上下載其他9種植物的類蛋白激酶基因編碼的氨基酸序列，利用DNAman構(gòu)建系統(tǒng)發(fā)育樹，結(jié)果表明沙田柚類蛋白激酶基因編碼的蛋白質(zhì)與蕓香科的克萊門柚（Citrus clementina）和甜橙（Citrus sinensis）的類蛋白激酶基因編碼的蛋白質(zhì)相似性很高，相似度約為100%和99%。該同源樹顯示該物種與克萊門柚和甜橙親緣關(guān)系很近，屬于同一分支（圖7）。

2.8 類蛋白激酶基因在沙田柚自交和異交花柱中的表達

該研究克隆的類蛋白激酶在沙田柚自交花柱和異交花柱中的表達有較為明顯的差異：在未授粉花柱中其表達量（reads per kb per million reads，RPKM） [20]為 0.01，自交1 d花柱中的表達量上升至1.61，自交2 d花柱中的表達量繼續(xù)上升至5.97，自交3 d的其表達量上升至6.23。而在異交1 d花柱中的表達量上升至4.06，異交2 d花柱中的表達量繼續(xù)上升至6.41，但異交3 d花柱中的表達量迅速下降至1.81。從上述結(jié)果可以看出無論自交授粉還是異交授粉，第1天花柱中類蛋白激酶基因的表達都表現(xiàn)為迅速升高，利用 RPKM對該基因在自交和異交花柱中的表達水平進行估算，設(shè)定錯誤檢測率（false discovery rate，F(xiàn)DR）≤0.001 且差異倍數(shù) ≥2 倍的基因視作顯著差異表達基因， 以 RPKM 值取 2 的對數(shù)值（log2）對自交1 d/異交 1 d、自交2 d/異交2 d，以及自交3 d/異交3 d花柱中類蛋白激酶基因的表達水平進行統(tǒng)計，其數(shù)值分別為-1.33、-0.10、1.78。

3 討論

蛋白激酶是一類催化蛋白質(zhì)磷酸化反應的酶，能將腺苷三磷酸（ATP）上的γ-磷酸轉(zhuǎn)移到蛋白質(zhì)分子的氨基酸殘基上[21]。類蛋白激酶蛋白的磷酸化以絲氨酸磷酸化為主，其次為蘇氨酸和酪氨酸，且與STKc IRAK蛋白有相同的保守結(jié)構(gòu)域，屬于PKc like超家族。該研究通過序列同源性分析，沙田柚類蛋白激酶基因編碼的氨基酸與克萊門柚（Citrus clementina）和甜橙（Citrus sinensis） 的同源性分別為 100%、99%，表明Unigene10208_All為類蛋白激酶基因。

植物蛋白激酶在植物體內(nèi)起著信號轉(zhuǎn)導的作用[17，22-23]。杜秀貞等[10]發(fā)現(xiàn)小麥類蛋白激酶基因片段TA50-10可以誘導植物抗病，且具有譜性。蛋白激酶在生長發(fā)育方面也有顯著的作用。李濤等[24]發(fā)現(xiàn)類蛋白激酶基因OsAGWI可以通過控制水稻穎殼外表皮細胞以及維管束的發(fā)育來實現(xiàn)對籽粒大小的控制。該研究獲得的類蛋白激酶基因自交1 d花柱中的表達量上升至1.61，自交2 d花柱中的表達量繼續(xù)上升至5.97，自交3 d的表達量仍上升至6.23。而在異交1 d花柱中的表達量上升至4.06，異交2 d花柱中的表達量也繼續(xù)上升至6.41，但異交3 d花柱中的表達量卻下降至1.81。自交 1 d/異交 1 d、自交3 d/異交3 d花柱的 log2（RPKM ratio）分別為-1.33，1.78，達到差異表達的水平。此外，該基因與沙田柚授粉相關(guān)，在未授粉花柱中其表達量非常低（0.01），自交授粉1 d時達1.61，異交授粉1 d也達4.06。但該基因與沙田柚自交不親和性的關(guān)系尚待進一步研究。

參考文獻

[1] DE NETTANCOURT D.Incompatibility in angiosperms[J].Sexual plant reproduction，1997，10（4）：185-199.

[2] 秦新民，張渝，劉玉潔，等.沙田柚S-RNase基因的克隆及序列分析[J].廣西師范大學學報（自然科學版），2015，33（1）：139-145.

[3] 秦新民，李惠敏，薛妙男，等.沙田柚自交、異交花粉管蛋白的雙向電泳分析[J].廣西植物，2004，24（6）：566-569.

[4] 楊繼華，李紅艷，薛妙男.沙田柚花柱S- 糖蛋白的分離及鑒定[J].廣西師范大學學報（自然科學版），2000，18（4）：66-70.

[5] 楊繼華，饒桂榮，薛妙男.沙田柚花柱S- 糖蛋白的純化和 N-端序列測定[J].廣西師范大學學報（自然科學版），2001，19（1）：72-79.

[6] 薛妙男，陳騰土，楊繼華.沙田柚自交和異交親和性觀察[J].園藝學報，1995，22（2）：127-132.

[7] 薛妙男，李義平，張杏輝，等.沙田柚自交花柱中識別蛋白的免疫金定位[J].園藝學報，2001，28（1）：59 -61.

[8] 秦新民，莫花濃，萬珊，等.沙田柚花粉管特異蛋白的免疫細胞化學研究[J].廣西師范大學學報（自然科學版），2008，26（4）：112-115.

[9] 秦新民，莫花濃，石菁萍，等.沙田柚花粉管 S1-RNase 免疫膠體金定位研究[J].廣西農(nóng)業(yè)科學，2009，40（5）：483-485.

[10] 杜秀貞，李培培，潘小玫，等.小麥類蛋白激酶基因片段TA50-10的克隆及其TMV抗性誘導功能[J].江蘇農(nóng)業(yè)學報，2011，27（6）：181-185.

[11] KOMATSU S，YANG G X，KHAN M，et al.Overexpression of calcium-dependent protein kinase 13 and calreticulin interacting protein 1 confers cold tolerance on rice plants[J].Molecular genetics & genomics，2007，277（6）：713-723.

[12] MIZUNO S，OSAKABE Y，MARUYAMA K，et al.Receptorlike protein kinase 2 （RPK 2） is a novel factor controlling anther development in Arabidopsis thaliana[J].Plant journal，2007，50（5）： 751-766.

[13] BECRAFT P W.Receptor kinases in plant development [J].Trends in plant science，1998，3（10）： 384-388.

[14] YOSHIMURA S，YAMANOUCHI U，KATAYOSE Y，et.al.Expression of Xa1，a bacterial blightresistance gene in rice，is induced by bacterial inoculation[J] Proc Natl Acad Sci USA，1998，95（4）：1663-1668.

[15] CHEN X W，SHANG J J，CHEN D X，et.al.A Blectin receptor kinase gene conferring rice blast resistance[J].Plant journal，2006，46（5）： 794-804.

[16] MORILLO S A，TAX F E.Functional analysis of receptorlike kinases in monocots and dicots[J].Current opinion in plant biology，2006，9（5）： 460-469.

[17] TORII K U.Leucinerich repeat receptor kinases in plants：Structure， function， and signal transduction pathways [J].International review of cytology，2004，234：1-46.

[18] 孫德權(quán)，郭啟高，胡玉林，等.改良Trizol法提取香蕉葉片總RNA[J].廣東農(nóng)業(yè)科學，2009（5）：162-164.

[19] KELLEY L A， MEZULIS S， YATES C M， et al.The Phyre2 web portal for protein modeling， prediction and analysis [J].Nature protocols， 2015， 10（6）：845-858.

[20] MORTAZAVI A，WILLIAMS B A，MCCUE K，et al.Mapping and quantifying mammalian transcriptomes by RNASeq[J].Nature methods，2008，5（7）：621-628.

[21] 張春寶，趙麗梅，趙洪錕，等.植物蛋白激酶研究進展[J].生物技術(shù)通報 ，2011（10）：17-23.

[22] NAKAGAMI H，PITZSCHKE A，HIRT H.Emerging MAP kinase pathways in plant stress signaling[J].Trends in plant science，2005，10（7）：339-346.

[23] STONE J M，WALKER J C.Plant protein kinase families and signal transduction[J].Plant physiology，1995，108（2）：451-457.

[24] 李濤，江潔明，張盛春，等.類蛋白激酶基因OsAGSW1調(diào)控水稻籽粒發(fā)育的研究[C]//廣東省遺傳學會.廣東省遺傳學會第九屆代表大會暨學術(shù)研討會論文及摘要匯編.廣州：廣東省遺傳學會，2014：1.