梁國平，李文芳，陳佰鴻，左存武，馬麗娟，何紅紅，萬鵬，安澤山，毛娟

?

不同糖源對葡萄試管苗蛋白激酶相關基因表達的影響

梁國平，李文芳，陳佰鴻，左存武，馬麗娟，何紅紅，萬鵬，安澤山，毛娟

（甘肅農業(yè)大學園藝學院，蘭州 730070）

【目的】探究不同外源糖對葡萄試管苗生長發(fā)育及蛋白激酶基因轉錄調控的影響，應用轉錄組測序挖掘蛋白質磷酸化過程中的基因,為葡萄蛋白激酶相關基因功能的驗證奠定一定基礎?！痉椒ā吭诨九囵B(yǎng)基中分別添加2%的蔗糖、葡萄糖和果糖，以無糖為對照,分別命名為S20、G20、F20和CK，經過37 d培養(yǎng)后，測定不同處理的地上和地下鮮重，并采用Illumina HiSeqTM2000對各處理葉片進行轉錄組測序，通過綜合生物信息學分析（參考基因組比對、差異基因（DEGs）篩選、COG（Cluster of Orthologous Groups of proteins）注釋、GO（Gene Ontology）注釋等）篩選出蛋白激酶相關基因，通過qRT-PCR分析該蛋白激酶相關基因的表達特性?！窘Y果】F20、G20和S20處理下的葡萄（‘紅地球’）試管苗與CK相比，地上鮮重具有明顯差異，且F20最高，而G20地下鮮重最高。SNP統(tǒng)計發(fā)現(xiàn)，轉換是主要的變異類型，顛換次之，且發(fā)生在基因間區(qū)的SNP數(shù)量最多，其次是下游；剪接位點供體和同義終止發(fā)生的基因數(shù)量最少且相等。4個樣品中共獲得了2 633個差異基因，3個處理與CK相比，共有差異基因180個且被聚類為3組，第一組中127個基因僅在CK中高表達，第二組19個基因僅在G20下高表達，而第三組34個基因在3個處理下表達模式不盡相同。這些共有的差異基因在COG中注釋到了26個基因并分在11個功能類別中，且主要注釋在一般的功能類別中。在GO分類中，共有的基因分別被注釋在分子功能、生物學過程和細胞組分的14、22和13個功能類別中。共篩選出7種蛋白激酶，分別為葡萄糖激酶（Glucokinase，GK）、絲裂原活化蛋白激酶（Mitogen-activated protein kinases，MAPKs）、鈣調蛋白激酶（Calcineurin protein kinase，CBL）、蛋白磷酸酶2（Protein phosphatase 2，PP2A）、己糖激酶（Hexokinase，HXK）、組氨酸蛋白激酶（Histidine protein kinase，HPK）和酪氨酸激酶（Tyrosine kinase，TK），其不同激酶的基因在不同處理中具有各自的表達模式，經qRT-PCR驗證，選擇的20個差異基因中有17個基因表達與轉錄組測序結果相一致?！窘Y論】在葡萄試管苗培養(yǎng)中，果糖較葡萄糖和蔗糖相比對生長較好。測序得出180個差異基因對3種不同糖均作出響應，這些基因在COG數(shù)據(jù)庫中主要富集在膜酯轉運和代謝、次級代謝物和碳水化合物的合成、轉運和分解；GO中大多注釋在蛋白激酶和氧化還原酶的活性中；篩選出了7種蛋白激酶，這些差異基因在數(shù)量、功能分類和代謝通路上對糖的響應各不相同。

RNA-Seq；外源糖；蛋白激酶；信號轉導

0 引言

【研究意義】在植物組織培養(yǎng)過程中，糖不僅調節(jié)細胞內的滲透勢，以維持水分平衡[1]，而且還是組培植物能量的供應者[2]和信號分子[3]。SOLAROVA等[4]指出，糖是初期組培苗幼芽生長和葉片光合作用所必須的。在葡萄良種繁育和脫毒苗培育過程中，糖對試管苗的形態(tài)建成具有決定性作用[5]。然而，在常規(guī)培養(yǎng)過程中，試管苗處在自養(yǎng)與異養(yǎng)的混合狀態(tài)[6]，且不同外源糖的供應對試管苗的生長也具有不同的作用。蛋白激酶可將ATP分子上γ位的磷酸基團轉移并共價結合到底物蛋白質氨基酸殘基（絲氨酸、蘇氨酸或絡氨酸）上，從而實現(xiàn)蛋白質磷酸化并引起蛋白質構象的改變，進而調控蛋白質的活性和功能，最終對物質轉運、能量代謝和信號轉導等過程起到開關作用[7]。因此，通過轉錄組測序深入分析不同外源糖對葡萄試管苗生長發(fā)育的調控，特別是對蛋白激酶基因表達的影響，為在組培條件下深入了解糖對植株的形態(tài)建成和進一步優(yōu)化碳源具有重要意義?！厩叭搜芯窟M展】在組織培養(yǎng)誘導植物產生不定根的過程中，培養(yǎng)基的初始糖濃度對其生長有很大影響，低濃度糖難以滿足植物生長的需要，而高濃度糖則導致細胞內滲透壓過高，最終對不定根的生長產生抑制[8]。鄒英寧[9]在培養(yǎng)基中分別添加蔗糖、葡萄糖、山梨醇、果糖4種糖源對中國李進行誘導增殖的研究中得出，葡萄糖濃度為30 g·L-1和15 g·L-1對促進增殖及生長和生根率及生根倍數(shù)最好。在誘導白樺愈傷組織的過程中發(fā)現(xiàn)，培養(yǎng)基中分別添加20—30 g·L-1的蔗糖、葡萄糖和果糖后，其鮮重積累量均顯著高于對照[10]。王思瑤等[11]研究了不同外源糖對檉柳叢生芽生長的影響，蔗糖濃度為10 g·L-1時有助于其叢生芽的生長。然而，KOZAI等[12]研究發(fā)現(xiàn)，將含有一片葉的單芽莖段葡萄試管苗繼代在無糖培養(yǎng)基中也可生長。這些結果表明，不同外源糖對組培植物的生長發(fā)育具有不同的作用，且目的不同需求糖的種類也不盡相同。植物蛋白激酶在信號感知、傳導、生長發(fā)育以及基因的表達調控中起重要作用。葡萄糖激酶（Glucokinase，GK）對葡萄糖具有特異性，在水稻[13]、豌豆[14]和馬鈴薯[15]中研究發(fā)現(xiàn)，該酶可優(yōu)先磷酸化葡萄糖，從而提前進入糖酵解途徑。促分裂原活化蛋白激酶（Mitogen-activated protein kinases，MAPK）是一種絲氨酸/蘇氨酸類蛋白激酶，它介導多種生物和非生物脅迫的相關信號通路，以磷酸化方式實現(xiàn)信號的級聯(lián)放大，最后引起相應的生物學效應，對植物的生長發(fā)育起重要的調控作用[16]。類鈣調磷酸酶B蛋白作為（Calcineurin B-like protein，CBL）Ca2+傳感蛋白之一[17]，以絲氨酸/蘇氨酸蛋白激酶CIPK（CBL-interacting protein kinase）作為靶蛋白[18-19]，并形成蛋白激酶復合物CBL-CIPK，調節(jié)下游基因的表達來維持細胞離子平衡[20]。在擬南芥突變體中研究發(fā)現(xiàn)，CBL-CIPK信號通路參與調控鹽脅迫下Na+穩(wěn)態(tài)分布和Na+、K+平衡[21]。蛋白磷酸酶2A（Protein phosphatase 2，PP2A）包含一系列絲氨酸/蘇氨酸磷酸酶[22]，該家族從酵母到哺乳動物在功能和結構上較為保守，并且在調節(jié)細胞的多種功能中起關鍵作用[23]。己糖激酶（Hexokinase，HXK）催化己糖磷酸化并使之參與糖酵解[24]。KARVE等[25]用RT-PCR法研究擬南芥HXK基因中不同家族成員在6種器官（幼葉、成熟葉、根、花、角果和莖）的表達情況發(fā)現(xiàn)，在所有器官中都具有較高的豐度；在葉中高表達，但在庫組織中表達略低。在根和角果中豐度高，但低于。和表達量低于，且這兩者在不同器官中表達豐度相近，僅在花中有表達。組氨酸蛋白激酶（Histidine protein kinase，HPK）是一個磷酸化組氨酸保守殘基的信號傳導酶家族[26]，與它們的下游靶蛋白一起構成了雙組分信號傳導系統(tǒng)，響應各種生物和非生物脅迫[27]。酪氨酸激酶（Tyrosinekinase，TK）是細胞信號轉導過程中最為重要的物質之一，對細胞調節(jié)、通訊和發(fā)育生物學方面具有重要的意義。在動物中研究發(fā)現(xiàn)，該酶的活性過高會導致其下游信號途徑激活，從而使細胞轉化、增殖、對抗細胞凋亡促進細胞生存，形成非正常組織[28]，而植物中報道甚少。【本研究切入點】雖然前人已在組培中對碳源影響植物生長發(fā)育方面做了大量報道，目的是確定最適碳源和濃度，以達到商品化生產和降低組培成本，然而在分子水平研究不同外源糖對其生長的影響，特別是對蛋白磷酸化中基因表達的研究鮮有報道?！緮M解決的關鍵問題】應用高通量測序手段，在分子層面深入分析在不同外源糖條件下‘紅地球’葡萄（L Red Globe）試管苗葉片中轉錄本的變化，特別是蛋白激酶相關基因如何響應不同的外源糖，為研究糖對植物生長發(fā)育的調控奠定一定的理論基礎。

1 材料與方法

1.1 試驗材料及處理

本試驗于2016—2017年在甘肅農業(yè)大學園藝學院果樹生理與生物技術實驗室進行，試驗材料為本實驗室保存的‘紅地球’試管苗。

選擇繼代健壯、長勢一致且無污染的試管苗，剪成單芽莖段（2—3 cm）并留一片葉，接種在不含糖以及含2%不同外源糖的50 mL MS培養(yǎng)基中，培養(yǎng)基分裝在150 mL的三角瓶中，每瓶接1個外植體，外源糖分別為蔗糖、葡萄糖和果糖，以無糖為對照，基本培養(yǎng)基組分為：4.43 g·L-1MS+6 g·L-1瓊脂+0.2 mg·L-1IAA，pH為5.8—6.0，分別命名為CK、S20、G20和F20，每個處理接種30瓶，培養(yǎng)條件為（25±1）℃，白熾光36 mmol·m-2·s-1光照16 h，黑暗8 h。經過37 d培養(yǎng)后，每個處理隨機挑選9瓶長勢一致的試管苗，測定其生物量，并迅速收集葉片液氮速凍，-80℃保存用于RNA的提取。每個處理進行3次生物學重復。

1.2 生物量測定

將培養(yǎng)37 d后的試管苗從培養(yǎng)基中取出，不損傷根系，從基部第一個芽處剪斷，測定其地上鮮重。用蒸餾水將根系表面的培養(yǎng)基洗干凈，后將根系在濾紙上輕輕滾動一圈，吸干表面多余水分，用于測定其地下鮮重，每個處理重復3次。

1.3 RNA提取和質檢

RNA提取采用植物提取試劑盒（中科泰瑞生物技術有限公司，北京）并按操作說明書進行。用Nanodrop和Agilent 2100（安捷倫有限責任公司，美國）檢測RNA的純度、濃度及完整性。

1.4 cDNA文庫的構建

cDNA文庫構建參照本課題組前期的試驗方案進行[29]。樣品檢測合格后，用帶有Oligo（dT）的磁珠，通過A-T互補配對與mRNA的ployA尾結合的方式富集mRNA。隨后加入Fragmentation Buffer將mRNA打斷成短的片段。以打斷的mRNA為模板，用六堿基隨機引物合成第一條cDNA鏈，然后加入緩沖液、dNTPs和DNA polymerase I合成第二條cDNA鏈，隨后利用AMPure XP beads純化雙鏈cDNA。純化的雙鏈cDNA再進行末端修復、加A尾并連接測序接頭，然后用AMPure XP beads進行片段大小選擇，最后通過PCR富集得到cDNA文庫。

1.5 cDNA庫檢及測序

文庫構建完成后，先使用Qubit 2.0進行初步定量，稀釋文庫至1 ng·μL-1，隨后使用Agilent 2100對文庫插入片段的長度進行檢測，插入片段大小符合預期后，使用qPCR方法對文庫的有效濃度進行準確定量（文庫有效濃度＞2 nmol·L-1）。

庫檢合格后，采用Illumina HiseqTM2000進行測序，測序讀長為PE100。

1.6 測序數(shù)據(jù)過濾

去除帶測序接頭的reads、無法確定堿基信息的比 例大于10%的reads及低質量reads（Qphred≤20的堿基數(shù)占整個read長度50%以上的reads）。

1.7 參考序列比對分析

利用TopHat將Clean Reads與指定的參考基因組進行序列比對[30]，獲取在參考基因組或基因上的位置信息以及測序樣品特有的序列特征信息。

1.8 差異基因篩選

通過設定FDR（False discovery rate）＜0.01和∣Log2（FC）∣≥2作為篩選差異基因的閾值進行差異基因的篩選[31]。

1.9 qPCR驗證

用差異基因的Gene ID對葡萄數(shù)據(jù)庫（http://www. genoscope.cns.fr/externe/GenomeBrowser/Vitis/）檢索并獲得相應的CDS區(qū)，通過在線引物設計軟件Primer 3 Input（http://primer3.ut.ee/）進行引物設計，引物見表1。對1.3的RAN使用Primer ScriptTMRT regent Kit with gDNA Eeaser（TaKaRa）試劑盒并按照說明書進行反轉錄獲得cDNA，并用SYBR Primer Ex TaqTMII（TaKaRa）試劑盒說明書進行定量，定量PCR儀為Light Cycler? 96 Real-Time PCR System（Roche，瑞士）。內參基因為（GenBank accession no. CB973647）。引物由生工生物工程股份有限公司合成。反應程序為：95℃預變性30 s，95℃變性10 s，60℃退火30 s，72℃延伸30 s，40個循環(huán)，重復3次。反應結束后分析熒光值變化曲線及熔解曲線?；虻南鄬Ρ磉_量采用2-ΔΔCT計算[32]。

2.0 統(tǒng)計分析

試驗結果數(shù)據(jù)經Excel 2010整理，統(tǒng)計分析采用單因素ANOVO的Duncan’s法，顯著性水平為＜0.05，用Origin 9.0作圖。

2 結果

2.1 不同外源糖對葡萄試管苗生長量的影響

以CK為對照，在2%的不同外源糖處理下，葡萄試管苗經37 d培養(yǎng)后，地上和地下生長具有明顯差異。其中地上鮮重表現(xiàn)出顯著差異，且在F20下生物量最大，而S20和G20處理下無顯著變化。地下的鮮重在S20和G20下沒有顯著差異，而F20與兩個處理均表現(xiàn)出顯著差異（圖1-A、B）。

2.2 測序數(shù)據(jù)統(tǒng)計

對4個樣品的轉錄組測序共獲得21.53 Gb Clean Data，各樣品Clean Data均達到4.69 Gb，不同處理獲得的Clean reads分別是：CK為19 840 139條，S20為19 619 720條，G20為24 923 975條和F20為20 579 635條。Q30堿基百分比均在90.26%以上（表2）。表明測序質量較好，可以用于后期的數(shù)據(jù)分析。

2.3 SNP/InDel分析

SNP是個體基因組DNA同一位置單個核苷酸變異（轉換、顛換、插入和缺失）所引起的多態(tài)性。對不同外源糖處理葡萄測序結果的單個核苷酸變異位點進行統(tǒng)計（表3），4個樣品的SNP數(shù)量分別為：403 332、489 805、358 014和331 527個，且主要發(fā)生的變異類型是轉換，其頻率在65.89%—67.18%，而顛換頻率在32.82%—34.11%。進一步對不同的SNP事件數(shù)量進行統(tǒng)計，結果表明發(fā)生在基因間區(qū)的SNP數(shù)量最多，為238 351個；其次是發(fā)生在下游的，為42 382個；而剪接位點供體和同義終止發(fā)生最少，且相等。（圖2）。

不同小寫字母表示不同糖處理水平上差異顯著（P＜0.05）。下同

表2 轉錄組測序數(shù)據(jù)

表3 SNP位點統(tǒng)計

圖2 SNP注釋分析

通過InDel注釋分類，結果如圖3所示。2%的不同外源糖處理‘紅地球’測序結果與參考基因組比對后，共獲得15種InDel類型，有21 479個差異表達基因的InDel類型為其他類型，而發(fā)生在基因間區(qū)的數(shù)量最多，為16 824個，其次是在下游，為3 343個；密碼子改變與缺失、終止獲得及剪接位點受體這3類事件的數(shù)量最少，都為1個。

2.4 可變剪接事件分析

對每個樣品存在的可變剪接類型統(tǒng)計，結果如圖4所示。不同外源糖處理下的葡萄試管苗中，其第一個外顯子可變剪接和最后一個外顯子可變剪接事件數(shù)量最多，其次是單內含子滯留，而多外顯子跳躍（模糊邊界）和單外顯子跳躍（模糊邊界）事件數(shù)量最少。

2.5 差異基因表達分析

通過設定FDR和Log2（FC）兩個閾值對不同外源糖處理的葡萄試管苗測序結果進行差異基因篩選，共獲得2 633個差異基因，其中共有的差異基因為180個（圖5-A）。對共有的差異基因以Log2（FPKM）值進行表達量分析（圖5-B），所有基因被聚類為3組，第一組包括127個基因，僅CK中高表達，而其余處理下均下調表達；第二組包括19個基因，僅在G20下高表達，其余處理中均下調；第三組的34個基因，與CK相比表達模式不盡一致，且表達在處理之間具有明顯差異。

圖3 InDel注釋分析

圖4 可變剪接分析

2.6 差異基因的COG分析

將共有的180個差異基因進行COG注釋，結果僅注釋到了26個，包括6個上調基因和20個下調基因（圖6）。其中，注釋到一般功能預測的差異基因數(shù)量最多，為10個，其次是次級代謝物的合成、轉運和分解，包括4個上調和3個下調差異基因。在轉錄，能量產生和轉運，氨基酸轉運和代謝、翻譯后修飾、蛋白質周轉、分子伴侶和防御機制4個功能分類中各注釋到1個基因。

圖5 差異基因的Venn圖和共有差異基因的表達分析

圖6 差異基因COG分析

2.7 差異基因的GO富集分析

通過對共有的差異基因進行GO功能注釋，發(fā)現(xiàn)這些差異基因主要被注釋在分子功能（14）、生物學過程（22）和細胞組分（13）功能類別中（圖7）。在分子功能中，注釋在ATP結合（ATP binding (GO: 0005524)）、絲氨酸/蘇氨酸蛋白激酶活性（protein serine/threonine kinase activity (GO: 0004674)）和氧化還原酶活性（oxidoreductase activity (GO: 0016491)）的基因數(shù)量最多；生物學過程中，注釋在蛋白質磷酸化（protein phosphorylation (GO: 0006468)）、氧化還原過程（oxidation-reduction process (GO: 0055114)）和碳水化合物代謝過程（carbohydrate metabolic process (GO: 0005975)）的基因個數(shù)分別為23、20和10個，而在細胞組分中，細胞外區(qū)域（extracellular region (GO: 0005576)）、質膜（plasma membrane (GO: 0005886)）占比最高。在分子功能中也注釋到了蛋白激酶活性（protein kinase activity (GO: 0004672)）、6-磷酸果糖激酶活性（6-phosphofructokinase activity (GO: 0003872)）和生長因子激活（growth factor activity (GO: 0008083)）等。生物學過程中也注釋到了負調控細胞程序性死亡（negative regulation of programmed cell death (GO: 0043069)）、果糖6-磷酸代謝過程（fructose 6-phosphate metabolic process (GO: 0006002)）及生長發(fā)育等（developmental growth (GO: 0048589)）功能類別。

2.8 7個蛋白激酶差異基因表達

蛋白激酶對植物的生長具有重要的作用。通過對本次測序結果進行篩選，獲得了7種與植物生長發(fā)育相關的蛋白激酶基因，包括GK的2個基因、MAPKs的34個基因、CBL的6個基因、PP2A的6個基因、HXK的6個基因、HPK的14個基因和TK的164個基因（圖8-A—G，G圖基因名稱見電子版附件）。與CK相比，2個GK基因在S20處理中上調表達（圖8-A）。MAPKs中的、、、、和在CK中上調表達，而在其余處理中均下調表達，S20下的4個基因包括、、和較CK、G20和F20表達量高（圖8-B）。CBL中的分別在CK和S20中出現(xiàn)下調，而在S20和F20中均上調表達（圖8-C）。有5個PP2A，包括、、、及僅在G20中上調表達，其余處理下均出現(xiàn)下調（圖8-D）。在HXK中，與CK相比在S20中僅有和為下調表達，其余的4個基因為上調表達（圖8-E）。在HPK中，僅在CK中高表達，而在S20、G20和F20中出現(xiàn)低表達，、、、和只在S20中上調，而、、、和只在G20中上調（圖8-F）。TK的164個基因在不同的處理中，表達具有明顯差異（圖8-G）。

2.9 qRT-PCR驗證

選擇了20個差異基因包括MAPK（5個）、HXK（1個）、PP2A（4個）、HPK（2個）和TK（8個）進行qRT-PCR分析，結果如圖9所示。有3個基因分別為GSVIVG01015297001（）GSVIVG01009192001（）和GSVIVG01014110001（）在不同的外源糖處理下，表達趨勢與轉錄組測序結不一致，其余的17個基因基本一致，表明測序結果可靠。

3 討論

糖對植物的形態(tài)建成和基因表達調控是其代謝產生的效應。近年來研究表明，糖作為信號分子在胞內形成不同信號轉導途徑，將信號轉導到細胞核或細胞器中目標基因的轉錄起始位點上，控制相關基因的轉錄，引起相應的生理生化變化，從而調控植物的生長和發(fā)育[33]。本次測序共獲得21.53 Gb Clean Data，各樣品Clean Data均達到4.69 Gb。對SNP統(tǒng)計發(fā)現(xiàn)，獲得轉換和顛換的SNP覆蓋率分別為65.89%—67.18%和32.82%—34.11%不等。在注釋到的15種InDel類型中，除21 479個差異基因的InDel類型未知外，最多的是發(fā)生在基因間區(qū)和下游中，獲得終止、密碼子改變與缺失和剪接位點受體這三類事件的數(shù)量最少，均為1個差異基因?？勺兗艚幽軌蚴挂粋€轉錄單位形成多個轉錄本和蛋白同系物，從而增加轉錄本和蛋白質的復雜性，這是導致真核生物基因和蛋白質數(shù)量較大差異的重要原因[34]。對本次測序結果進行可變剪接分析得出，12種可變剪接事件中TTS和TSS數(shù)量最多，其次是IR，而MSKIP和XMSKIP事件最少。

進行COG注釋，可從基因組水平上尋找直系同源體，預測未知的ORF生物學功能，并提高基因注釋的準確性和完整性[35]。對篩選出的180個共有的差異基因通過COG注釋，共獲得了26個基因，包括6個下調和20個上調，且注釋到最多的是一般功能預測，而轉錄相關和信號轉導機制功能類別次之，最少的則是核苷酸轉運與代謝和染色體結構與活力，在細胞核結構、細胞運動和胞外結構中均沒有注釋到差異基因，這一結果與張雪等[36]研究山葡萄果皮著色的轉錄組注釋結果基本一致，表明轉錄和信號轉導對植物生長發(fā)育和果實著色的調控具有重要作用。GO注釋分類發(fā)現(xiàn)這些差異基因主要被注釋在分子功能、生物學過程和細胞組分的12、22和13個功能類別中。在分子功能中，注釋到絲氨酸/蘇氨酸蛋白激酶活性的基因數(shù)量較多，其次是氧化還原酶活性；在生物學過程中，注釋到蛋白質磷酸化、氧化還原過程和碳水化合物代謝過程的基因個數(shù)分別為21、20和7個，而在細胞組分中，細胞外區(qū)域和質膜占比最高，達到11.2%和5.6%，這一結果與魏利斌等[37]在研究芝麻發(fā)育轉錄組分析結果一致。

大量研究表明蛋白質的磷酸化與去磷酸化過程在細胞的信號識別與轉導中起重要作用，它是生物體中普遍存在的一種調節(jié)機制，幾乎涉及所有的生理和病理過程，如糖代謝、細胞的生長發(fā)育、光合作用、基因表達等[38]。本研究通過設定差異基因閾值，篩選出了與植物生長發(fā)育相關的7種蛋白激酶，包括GKs、MAPKs、CBLs、PP2A、HXK、HPK和TK，這些蛋白激酶基因在不同外源糖處理下表達各不相同。

GKs能催化葡萄糖和其他己糖的磷酸化，包括果糖、甘露醇和半乳糖，也在大腸桿菌、酵母、哺乳動物和植物中充當信號感受器[39-40]。在擬南芥中研究發(fā)現(xiàn)，與線粒體胞膜相關的HXK1感知葡萄糖存在后，會特異性地下調編碼與光合作用過程中相關蛋白質基因的表達，例如Chl a/b結合蛋白（CAB）和Rubisco小亞基（RBCS）等[41-42]。在本研究中，GK中的兩個基因在S20處理中上調表達，在G20處理中下調，這一結果與CHO等[43]的研究一致。HXK作為雙功能酶存在于所有的有機體中，它不僅磷酸化己糖轉變?yōu)榧禾?磷酸，而且對糖信號的感知具有重要作用[44]。對測序結果進行HXK篩選得出，G20下分別有兩個和與其余3個處理相比均上調表達，表明這幾個基因較其余的處理更能快速感應葡萄糖的存在。

PP2A在植物中具有多重作用，例如對細胞周期性進程的調控，根皮質細胞伸長，組織發(fā)育和植物對生物和非生物脅迫的反應等[45]。RAZAVIZADEH等[46]研究發(fā)現(xiàn)，在不同程度的干旱脅迫和4 μmol·L-1的ABA處理下均表現(xiàn)出顯著上調。本研究中對PP2A進行差異基因篩選，獲得6個差異表達基因，其中、、、及僅在G20中上調表達，其余處理均出現(xiàn)下調，且與對照相比，在果糖中出現(xiàn)上調趨勢。擬南芥中研究發(fā)現(xiàn)，CIPK9與CBL3相互作用有助于在低K+脅迫下細胞K+的動態(tài)平衡，且CBL2和CBL3的過表達提高了對低K+的敏感性[47]，而本研究篩選到的6個CBL中，、和在G20下低表達，另一個、和卻表現(xiàn)出上調，表明這些基因在響應不同的外源糖時表達模式不盡相同。

在植物中，由環(huán)境條件或激素改變產生的信號均是通過MAPKs所介導的級聯(lián)反應進行傳遞的[48]。通過遺傳和生化方法研究發(fā)現(xiàn)，MAPKs不僅參與生物和非生物脅迫的信號轉導，而且還參與激素對生長發(fā)育信號的轉導過程[49]。如今，MAPKs已在大量植物物種中被發(fā)現(xiàn)，表明MAPKs級聯(lián)信號轉導機制在高等植物中也非常保守[50]。通過轉錄組測序，本研究篩選到了MAPKs的34個差異表達基因，其中、、、、和在CK中上調表達，而在其余處理下均下調表達。HPKs是一個大的信號轉導酶家族，能夠在保守的組氨酸殘基上自磷酸，并與其下游的靶蛋白形成雙組分信號系統(tǒng)響應生物和非生物脅迫[51]。本研究篩選到了10個組氨酸蛋白激酶基因，通過聚類分析發(fā)現(xiàn)，、、、和只在S20中上調，而、、、和只在G20中上調，表明這些基因對響應不同的糖具有明顯選擇性。酪氨酸（Tyr）磷酸化是動物體內受體酪氨酸激酶（RTKs）介導信號轉導的關鍵分子開關[52]，在植物中雖然缺乏典型的RTKs，但是擁有數(shù)百個與動物具有相似結構域結構的RTKs[53]。最近研究表明，在植物中酪氨酸磷酸化對于調控誘導免疫過程（PTI）的信號傳遞具有重要作用[54]。本研究篩選到酪氨酸蛋白激酶的164個差異基因，這些基因在響應不同的外源糖時具有各自的表達模式。

應用qRT-PCR對測序結果挑選了20個基因進行驗證，發(fā)現(xiàn)有17個基因的表達趨勢與測序結果一致，表明轉錄組測序結果可靠，這為通過基因表達研究蛋白激酶響應不同外源糖，并深入分析糖在植物生長發(fā)育中的作用提供了理論依據(jù)。

4 結論

本研究應用RNA-Seq技術對不同外源糖處理下的葡萄試管苗進行測序分析，共獲得了2 633個差異基因，有180個差異基因對3種不同外源糖都作出響應，且共有的這些基因表達對糖源具有選擇性。篩選到了7種蛋白質磷酸化相關的基因，這些差異基因在數(shù)量、功能分類和代謝通路上對糖的響應各不相同，且有各自的選擇表達特異性，為豐富糖在調節(jié)植物生長發(fā)育及進一步探究蛋白質磷酸化響應外源糖提供了參考。

[1] KADLECEK P, RANK B, TICHA I. Photosynthesis and photoprotection inL. in vitro-grown plantlets., 2003, 160(9): 1017.

[2] RUAN Y L. Sucrose metabolism: Gateway to diverse carbon use and sugar signaling., 2014, 65: 33-67.

[3] YUE C, CAO H L, WANG L, ZHOU Y H, HUANG Y T, HAO X Y, WANG Y C,WANG B, YANG Y J, WANG X C. Effects of cold acclimation on sugar metabolism and sugar-related gene expression in tea plant during the winter season., 2015, 88: 591-608.

[4] SOLAROVA J, POSPHILSILOVA J, CATSKY J, SANTRUCEK J. Photosynthesis and growth of tobacco plantlets independence on carbon supply., 1989, 23(4): 629-637.

[5] SMEEKENS S, MA J K, HANSON J, ROLLAND F. Sugar signals and molecular networks controlling plant growth., 2010, 13(3): 273-278.

[6] NGUYEN Q T, KOZAI T, NIU G, NGUYEN U V. Photosynthetic characteristics of coffee () plantletsin response to different CO2concentrations and light intensities., 1999, 55: 133-139.

[7] 張曦, 林金星, 單曉昳. 擬南芥無機氮素轉運蛋白及其磷酸化調控研究進展.植物學報, 2016, 51(1): 120-129.

ZHANG X, LIN J X, SHAN X Y. Progress in inorganic nitrogen transport proteins and their phosphorylation regulatory mechanism in, 2016, 51(1): 120-129. (in Chinese)

[8] 尹雙雙, 高文遠, 王娟, 劉輝, 左北梅. 藥用植物不定根培養(yǎng)的影響因素. 中國中藥雜志, 2012, 37(24): 3691-3694.

YIN S S, GAO W Y, WAN J, LIU H, ZUO B M. Influencing factors on culture of medicinal plants adventitious roots., 2012, 37(24): 3691-3694. (in Chinese)

[9] 鄒英寧. 不同碳源對中國李離體培養(yǎng)的影響. 安徽農業(yè)科學, 2010, 38(4): 1720, 1776.

ZOU Y N. Effect of different carbon sources on the culture of Chinese plum ().s, 2010, 38(4): 1720, 1776. (in Chinese)

[10] 王博, 范桂枝, 詹亞光, 李康. 不同碳源對白樺愈傷組織生長和三萜積累的影響. 植物生理學通訊, 2008, 44(1): 97-99.

WANG B, FAN G Z, ZHAN Y G, LI K. Effects of different carbon sources on callus growth and accumulation of triterpenoids in Birch ()., 2008, 44(1): 97-99. (in Chinese)

[11] 王思瑤, 崔瞳肸, 翟睿, 林香雨, 李欣, 孫璐, 詹亞光, 尹靜. 不同碳源對檉柳叢生芽生長、三萜及黃酮物質積累的影響. 植物生理學報, 2017, 53 (12): 2189-2196.

WANG S Y, CUI T X, ZHAI R, LIN X Y, LI X, SUN L, ZHAN Y G, YIN J. Effects of different carbon sources on the growth and accumulation of triterpenoids and flavonoids in tufted buds of., 2017, 53(12): 2189-2196. (in Chinese)

[11] 王愛民, 劉文, 傅中滇. 不同碳源對紅邊朱蕉組培苗生長的影響. 徐州師范大學學報 (自然科學版), 2003, 21(3): 76-78.

WANG A M, LIU W, FU Z D. A preliminary study on the reducing cost of tissue culture plantlet of cordyline terminalis.(), 2003, 21(3): 76-78. (in Chinese)

[12] KOZAI T, KUBOTA C, JEONGB R. Environmental control for the large-scale production of plants through in vitro techniques., 1997, 51: 49-56.

[13] GUGLIELMINETTI L, PERATA P, MORITA A, LORETI E, YAMAGUCHI J, ALPI A. Characterization of isoforms of hexose kinases in rice embryo., 2000, 53(2): 195-200.

[14] TURNER J F, CHENSEE, HARISSON D D. Glucokinase of pea seeds., 1997, 480: 367-375.

[15] MARTINEZ-BARAJAS E, RANGALL D D. Purification and characterization of a glucokinase from young tomato (. Mill.) fruct., 1998, 205: 567-573.

[16] 趙書平, 談宏斌, 鹿丹, 裴麗麗, 崔喜艷, 馬有志, 陳明, 徐兆師, 張小紅. 植物蛋白激酶介導的非生物脅迫和激素信號轉導途徑的研究進展. 植物遺傳資源學報, 2017, 18(2): 358-366.

ZHAO S P, TAN H B, LU D, PEI L L, CUI X Y, MA Y Z, CHEN M, XU Z S, ZHANG X H. Research progress of plant protein kinase mediated abiotic stress and hormone signal transduction pathway., 2017, 18(2): 358-366. (in Chinese)

[17] RANTY B, ALDON D, GALAUD J P. Plant calmodulins and calmodulin-related proteins: Multifaceted relays to decode calcium signals., 2006, 1(3): 96-104.

[18] BATISTIC O, KUDLA J. Integration and channeling of calcium signaling through the CBL calciumsensor/CIPK protein kinase network., 2004, 219: 915-924.

[19] WEINL S, KUDLA J. The CBL-CIPK Ca2+-decoding signaling network: Function and perspectives., 2009, 184(3): 517-528.

[20] LUAN S. The CBL-CIPK network in plant calcium signaling., 2009, 14: 37-42.

[21] YU Q Y, AN L J, LI W L. The CBL-CIPK network mediates different signaling pathways in plants., 2014, 33: 203-214.

[22] KENTARO F, TOMOYUKI F, SHUN-ICHI Y, TETSU S, YUSUKE K, YUTAKA Y, HIROMI K, HITOSHI N, TOMOTAKE K. The PP2A-like protein phosphatase ppg1 and the far complex cooperatively counteract CK2-mediated phosphorylation of Atg32 to inhibit mitophagy., 2018, 23: 3579-3590.

[23] 劉釗, 賈霖, 賈盟, 關明俐, 曹英豪, 劉麗娟, 曹振偉, 李莉云, 劉國振. 水稻PP2Ac類磷酸酶蛋白質在鹽脅迫下的表達. 中國農業(yè)科學, 2012, 45(12): 2339-2345.

LIU Z, JIA S, JIA M, GUAN M L, CAO Y H, LIU L J, CAO Z W, LI L Y, LIU G Z. Expression on profiling of rice PP2Ac type phosphatase proteins in seedlings under salt-stressed conditions., 2012, 45(12): 2339-2345. (in Chinese)

[24] 雍彬, 何兵, 徐攀, 虞傳洋, 王東. 甘薯基因的克隆、組織表達及生物信息分析. 四川大學學報(自然科學版), 2014, 51(2): 378-384.

YONG B, HE B, XU P, YU C Y, WANG D. Cloning, organ-specific expression pattern and sequence analysis ofgene in sweet potato.(), 2014, 51(2): 378-384. (in Chinese)

[25] KARVE A, RAUH B L, XIA X. Expression and evolutionary features of the hexokinase gene family in., 2008, 228(3): 411-425.

[26] HERICOURT F, CHEFDOR F, DJEGHDIR I, LARCHER M, LAFONTAINE F, COURDAVAULT V, AUGUIN D, COSTE F, DEPIERREUX C, TANIGAWA M, MAEDA T, GLEVAREC G, CARPIN S. Functional divergence of poplar histidine-aspartate kinase hk1 paralogs in response to osmotic stress., 2016, 17(12): 2061.

[27] WOLANIN P M, THOMASON P A, STOCK J B. Histidine protein kinases: Key signal transducers outside the animal kingdom., 2002, 3(10): reviews3013.1.

[28] HUPFELD T, CHAPUY B, SCHRADER V, BEUTLER M, VELTKAMP C, KOCH R, CAMERON S, AUNG T, HAASE D, LAROSEE P, TRUEMPER L, WULF G G. Tyrosinekinase inhibition facilitates cooperation of transcription factor SALL4 and ABC transporter A3 towards intrinsic CML cell drug resistance., 2013, 161(2): 204.

[29] MAO J, LI W F, MI B Q, DAWUDA M M, CALDERON-URREA A, MA Z H, ZHANG Y M, CHEN B H. Different exogenous sugars affect the hormone signal pathway and sugar metabolism in ‘Red Globe’ (L.) plantlets grown in vitro as shown by transcriptomic analysis., 2017, 246: 537-552.

[30] KIM D, PERTEA G, TRAPENLL C, PIMENTEL H, KELLEYR R, SALZBERG S L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions., 2013, 14(4): R36.

[31] ANDERS S, HUBER W. Differential expression analysis for sequence count data., 2010, 11(10): R106.

[32] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative and the 2-ΔΔCTmethod., 2001, 25: 402- 408.

[33] ROLLAND F, BAENA-GONZALEZ E, SHEEN J. Sugar sensing and signaling in plants: Conserved and novel mechanisms., 2006, 57: 675-709.

[34] FILICHKIN S A, HAMILYON M, DHARMAWARDHANA P D, SINGH S K, SULLIVAN C, BEN-HUR A, REDDY A S N, JAISWAL P. Abiotic stresses modulate landscape of poplar transcriptome via alternative splicing, differential intron retention, and isoform ratio switching., 2018. doi: 10.3389/fpls. 2018.00005.

[35] 賈新平, 孫曉波, 鄧衍明, 梁麗建, 葉曉青. 鳥巢蕨轉錄組高通測序及分析. 園藝學報, 2014, 41(11): 2329-2341.

JIA X P, SUN X B, DENG Y M, LIANG L J, YE X Q. Sequencing and analysis of the transcriptome of., 2014, 41(11): 2329-2341. (in Chinese)

[36] 張雪, 唐銘浩, 陳蒙, 李春艷, 劉海峰. 山葡萄不同著色時期果皮轉錄組測序分析. 果樹學報, 2017, 34(7): 781-789.

ZHANG X, TANG M H, CHEN M, LI C Y, LIU H F. Sequencing analysis of transcriptome ofduring different periods of coloration., 2017, 34(7): 781-789. (in Chinese)

[37] 魏利斌, 苗紅梅, 張海洋. 芝麻發(fā)育轉錄組分析. 中國農業(yè)科學, 2012, 45(7): 1246-1256.

WEI L B, MIAO H M, ZHANG H Y. Transcriptomic analysis of sesame development., 2012, 45(7): 1246-1256. (in Chinese)

[38] PROUD C G. Phosphorylation and signal transduction pathways in translational control., 2018. doi: 10.1101/cshperspect.a033050.

[39] IRWIN D M, TAN H. Evolution of glucose utilization: Glucokinase and glucokinase regulator protein., 2014, 70(1): 195-203.

[40] GUTIERREZNOGUES A, GARCIAHERRERO C M, ORIOLA J,VINCENT O, NAVAS M A. Functional characterization of MODY2 mutations in the nuclear export signal of glucokinase., 2018, 1864(7): 2385-2394.

[41] KAWAI S, MUKAI T, MORI S, MIKAMI B, MURATA K. Hypothesis: structures, evolution, and ancestor of glucose kinases in the hexokinase family., 2005, 99(4): 320-330.

[42] MOORE B, SHEEN J. Role of theglucose sensor, 2003, 300(5617): 332-336.

[43] CHO Y H, YOO S D, SHEEN J. Regulatory functions of nuclear hexokinase1 complex in glucose signaling., 2007, 127(2): 579-589.

[44] KELLY G, SADE N, DORON-FAIGENBOIM A, LERNER S, SHATIL-COHEN A, YESELSON Y, EGBARIA A, KOTTAPALLI J, SCHAFFER A A, MOSHELION M, GRANOT D. Sugar and hexokinase suppress expression of PIP aquaporins and reduce leaf hydraulics that preserves leaf water potential., 2017, 91: 325-339.

[45] HU R B, ZHU Y F, SHEN G X, ZHANG H. Overexpression of thegene confers increased salt tolerance in., 2017, 12(2): e1276687.

[46] RAZAVIZADEH R, SHOJAIE B, KOMATSU S. Characterization of, mRNA expression and growth patterns in, under drought stress and abscisic acid., 2018(1): 1-13.

[47] LIU L L, REN H M, CHEN L Q, WANG Y, WU W H. A protein kinase, calcineurin B-like protein-interacting protein Kinase9, interacts with calcium sensor calcineurin B-like Protein3 and regulates potassium homeostasis under low-potassium stress in., 2013, 161(1): 266-277.

[48] COLCOMBET J, HIRT H. Arabidopsis MAPKs: A complex signalling network involved in multiple biological processes., 2008, 413(2): 217-226.

[49] ANDREASSON E, ELLIS B. Convergence and specificity in theMAPK nexus., 2010, 15(2): 106-113.

[50] SINGHA H S, CHAKRABORTY S, DEKA H. Stress induced MAPK genes show distinct pattern of codon usage in,and., 2014, 10(7): 436-442.

[51] WOLANIN P M, THOMASON P A, STOCE J B. Histidine protein kinases: Key signal transducers outside the animal kingdom., 2002, 3(10): reviews3013.1.

[52] LIM W A, PAWSON T. Phosphotyrosine signaling: Evolving a new cellular communication system., 2010, 142: 661-667.

[53] LIN W W, LI B, LU D P, CHEN S X, ZHU N, HE P, SHAN L B. Tyrosine phosphorylation of protein kinase complex BAK1/BIK1 mediatesinnate immunity.2014, 111(9): 3632-3637.

[54] REDDY M M, RAJASEKHARAN R. Serine/threonine/tyrosine protein kinase fromis dependent on serine residues for its activity., 2007, 460(1): 122-128.

Effects of Different Sugar Sources on Protein Kinase Gene Expression in Grape Plantlets

LIANG GuoPing, LI WenFang, CHEN BaiHong, ZUO CunWu, MA LiJuan, HE HongHong, WAN Peng, AN ZeShan, MAO Juan

(College of Horticulture, Gansu Agricultural University, Lanzhou 730070）

【Objective】To explore the effects of different exogenous sugars on the growth and development of grape plantlets and the regulation of protein kinase gene transcription, the candidate genes were tapped in the process of protein phosphorylation by using transcription, which made a foundation for the verification of grape protein kinase-related gene function.【Method】Sucrose (2%), glucose (2%) and fructose (2%) were added to the basic medium, and the free-sugar treatment was as control, which were named as S20, G20, F20 and CK, respectively. After 37 days of culture, the fresh weight of the leaf-stem and root under different treatments was determined. Transcriptome sequencing of each treated foliages was performed by using Illumina HiSeqTM2000, and a series of protein kinases related genes were screened by integrated bioinformatics analysis, including reference genomic alignment, differentially expressed gene (DEGs) screening, COG (Cluster of Orthologous Groups of proteins) annotation, GO (Gene Ontology) annotation, etc., and the expression characteristic of these genes were further analyzed by qRT-PCR. 【Result】Compared with CK, ‘Red Globe’ grape plantlets under F20, G20 and S20 treatments exhibited significant differences in the fresh weight of leaf-stem, and the highest was obtained by F20 treatment, while the weight of fresh root under G20 was the highest. The SNP statistics found that the Transition was the main type of mutation, the second was Transversion. The highest number of SNPs that occurred in the Intergenic, and the next was the Upstream. Splice_Site_Donor and Synonymous_Stop events occurred with the least number of genes and equal. A total of 2 633 deferentially expressed genes were obtained in the 4 samples. The Venn diagram showed that there were a total of 180 differential genes under the 3 treatments compared with CK, and these genes were clustered into 3 groups. In the first group, 127 genes were only highly expressed under CK. The 19 genes of the second group were only highly expressed under G20, while the expression patterns of the 34 genes in third group were different under three treatments. The common 180 differential genes were annotated with 26 genes in the COG database to 11 functional categories, and these DEGs were mainly enriched in general functional categories. In the annotation of GO, the common genes were annotated in 14, 22 and 13 functional categories of molecular function, biological process and cellular component, respectively. Seven kinds of protein kinases were screened by this sequencing, including Glucokinase (GK), Mitogen-activated protein kinases (MAPKs), Calcineurin protein kinase (CBL), Protein phosphatase 2 (PP2A), Hexokinase (HXK), Histidine protein kinase (HPK) and Tyrosine kinase (TK), and these different protein kinases genes showed their own expression patterns among different treatments. By qRT-PCR analysis, 17 out of 20 screened genes expression were consistent with the transcriptome sequencing results. 【Conclusion】Compared with glucose and sucrose, fructose was the best sugar during grape culture process. The sequencing results showed that 180 DEGs all responded to three different sugars. In the COG annotation, these genes were mainly enriched in membrane ester transport and metabolism, the synthesis, transport and decomposition of secondary metabolites and carbohydrates. In the GO databases, the most of common DEGs were annotated in the activities of protein kinases and oxidoreductases. Seven protein kinases were identified, which were selectively in responses to different exogenous sugars in quantity, functional, category and metabolic pathways, and had their own choice of expression specificity.

RNA-seq; exogenous sugar; protein kinase; signaling transduction

10.3864/j.issn.0578-1752.2019.07.001

2018-07-25；

2018-10-15

國家自然科學基金（31460500）、甘肅省科技重大專項計劃（18ZD2NA006-4）、甘肅省現(xiàn)代水果產業(yè)體系崗位專家項目（GARS-SG-3）

梁國平，Tel：18298344227；E-mail：1143016341@qq.com。通信作者毛娟，E-mail：maojuan@gsau.edu.cn

（責任編輯 趙伶俐）