桂花己糖激酶基因家族成员的序列及表达分析

庞天虹，钱婕妤，付建新，顾翠花，张 超

（1. 浙江农林大学 风景园林与建筑学院，浙江 杭州 311300；2. 浙江农林大学 浙江省园林植物种质创新与利用重点实验室，浙江 杭州 311300；3. 浙江农林大学 南方园林植物种质创新与利用国家林业和草原局重点实验室，浙江 杭州 311300）

桂花Osmanthus fragrans是中国十大名花之一，深受人们喜爱。前人研究报道桂花花香主要成分为α-紫罗兰酮、β-紫罗兰酮、芳樟醇及其氧化物、罗勒烯等挥发性物质[1]，而桂花花色主要成分是类胡萝卜素化合物[2]。桂花的花香和花色是其最主要的观赏性状，与花香和花色相关次生代谢物质的合成与积累直接影响到桂花的观赏品质。植物呼吸作用不仅提供植物所需的能量，一系列的中间产物还为多种次生代谢化合物合成提供原料，在植物体内有机物转变方面起到枢纽作用。已糖激酶(hexokinase，HXK)具有己糖磷酸化功能，经HXK磷酸化的己糖才能进行植物呼吸代谢的糖酵解途径，为植物生长和发育提供能量和中间代谢产物[3-6]。高等植物HXK基因以多基因家族形式存在[7-11]。拟南芥Arabidopsis thaliana的6个HXK家族成员中，AtHXK1～3能磷酸化葡萄糖[10]，其中AtHXK1还具有感知和传递己糖信号的功能[12-13]。AtHXK1在强光条件下对植物根、叶和花序的生长和发育有促进作用[13]。此外，拟南芥成花转变[14]和种子萌发[15]、植物花青素积累[16]和脂肪族硫甙生物合成[17]等植物生命活动都与HXK介导的己糖信号转导息息相关。目前，未见桂花HXK基因功能有关的研究报道。本研究基于桂花不同花色品种转录组数据筛选OfHXKs相关Unigene序列，分析不同OfHXKs基因家族成员核苷酸序列及其编码的氨基酸序列，并与其他物种HXK蛋白进行多重比对和进化树分析，同时分析不同桂花品种不同发育阶段花瓣以及不同组织中OfHXKs基因表达特征，为进一步揭示桂花不同OfHXKs基因成员功能奠定理论基础。

1 材料与方法

1.1 材料

选择浙江农林大学桂花资源圃生长状况良好的地栽桂花‘堰虹桂’Osmanthus fragrans‘Yanhong Gui’(YHG)、‘玉玲珑’O. fragrans‘Yulinglong’(YLL)、‘金球桂’O. fragrans‘Jinqiu Gui’(JQG)为材料，依据桂花花序不同发育阶段[18]，对这3个桂花品种顶壳期(S1)、铃梗期(S2)、初开期(S3)、盛开期(S4)花瓣分别进行取样，取样时间均为10:00。同时还对‘堰虹桂’1年生茎、2年生茎、嫩叶、成熟叶和盛开时花序进行取样，用于不同组织中OfHXKs基因表达分析。用液氮速冻暂存样品，后存于-80 ℃冰箱。

1.2 主要试剂和仪器

RNA提取试剂盒采用RNAprep Pure Plant Kit(天根，北京)；反转录试剂盒使用PrimeScript™ RT Master Mix(Perfect Real Time)(Takara，大连)；荧光定量试剂盒使用TB Green®Premix Ex Taq™ Ⅱ(Tli RNaseH Plus)(Takara，大连)。实时荧光定量PCR仪为 LightCycler®480 Ⅱ(Roche，德国)。

1.3 方法

1.3.1OfHXKs序列分析 从转录组中筛选得到HXK相关的Unigene序列，首先使用美国国家生物信息中心(NCBI)数据库中的BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi)对得到的基因序列进行对比；利用ORF FINDER (https://www.ncbi.nlm.nih.gov/orffinder/)寻找开放阅读框；用DNAMAN 6.0软件对不同品种相同HXK基因成员进行核苷酸相似性比较，使用‘堰虹桂’的4个HXK基因成员序列开展后续序列分析，包括应用Prot-Param在线软件(http://web.expasy.org/protparam/)预测所编码蛋白的分子量、理论等电点(PI)、不稳定系数、总平均亲水性等；使用YLoc (https://abi-services.informatik.uni-tuebingen.de/yloc/webloc.cgi)进行亚细胞定位预测；使用SOPMA分析软件(https://npsa-prabi.ibcp.fr/cgi-bin/secpred_sopma.pl)对OfHXKs家族的蛋白序列进行二级结构预测分析，判断目的基因的氨基酸残基是否处于于α-螺旋、β折叠(或其他状态)的二级结构；用SWISS-MODEL (https://swissmodel.expasy.org/)预测OfHXKs的三级结构；应用SignalP 5.0 (http://www.cbs.dtu.dk/services/SignalP/)对信号肽进行预测；利用TMHMM在线软件(http://www.cbs.dtu.dk/services/TMHMM/)预测分析桂花蛋白的跨膜结构；使用NetPhos 3.1 Server (http://www.cbs.dtu.dk/services/NetPhos/)对桂花花瓣OfHXKs蛋白质磷酸化位点进行预测；采用MEGAX软件中的邻位相邻法(neighbor-joining，NJ)进行同源聚类，建立系统发育树，并采用Bootstrap法(重复1 000次)评估检测系统进化树。

1.3.2 总RNA提取与cDNA的合成 采用RNAprep pure Plant Kit(天根，北京)多糖多酚植物总RNA提取试剂盒按照说明书步骤提取各样品的总RNA。使用紫外分光光度计和质量浓度为1%琼脂糖凝胶电泳检测总RNA浓度和质量。反转录体系为总RNA 1.0 μL，5×PrimeScript RT Master Mix(TaKara，大连)2.0 μL，无RNA酶双蒸水7.0 μL，PCR反应的程序设定为：37 ℃ 15 min，85 ℃ 5 s，4 ℃ 1 min。反转录合成的cDNA储存于-20 ℃备用。

1.3.3 荧光定量PCR 根据荧光定量PCR引物设计原则利用Primer Premier 5进行引物设计，以桂花OfACT为内参基因[19]，引物序列见表1。荧光定量PCR反应体系为：TB Green Premix Ex TaqⅡ(TaKaRa，大连) 5.0 μL，上下游引物 (10 μmol·L-1)各 0.4 μL，cDNA 模板(反转录 cDNA 稀释20 倍)2.0 μL，双蒸水2.2 μL，每个样品设置3个生物学重复。反应程序为两步法：95 ℃预变性30 s，95 ℃变性5 s，60 ℃复性30 s，重复40个循环；然后以95 ℃持续5 s，60 ℃持续1 min，95 ℃持续15 s作为溶解曲线程序。

表1 OfHXKs基因表达分析所用引物序列Table 1 Primer sequences of OfHXK genes of O. fragrans

2 结果与分析

2.1 桂花花瓣 OfHXKs基因序列分析

根据桂花不同花色品种转录组Unigene序列分析，筛选得到4个HXK基因的同源基因(OfHXK1～OfHXK4)，3个品种中均含有这4个HXK基因的同源基因。利用DNAMAN对不同品种OfHXK1～OfHXK4核苷酸序列进行多序列比对，核苷酸序列相似性比较结果见表2。不同品种桂花OfHXK1、OfHXK3和OfHXK4基因核苷酸序列相似性均高于99%，而不同品种桂花OfHXK2基因核苷酸序列相似性为94.42%～98.26%。选择‘堰虹桂’的OfHXK1～OfHXK4基因推测的氨基酸序列进行后续分析。

表2 桂花 OfHXKs核苷酸相似性比较表Table 2 Comparison of nucleotide sequences of OfHXKs in O. fragrans

利用DNAMAN对‘堰虹桂’OfHXK1～OfHXK4氨基酸序列比对。由图1可知：不同OfHXKs间具有高度保守的结构域，均包含2个保守的磷酸化作用位点phosphate 1和phosphate 2以及1个底物结合位点sugar binding。预测OfHXK1与OfHXK2序列中均包含12个疏水通道和11个甘氨酸残基，而OfHXK3包含11个疏水通道和10个甘氨酸残基，OfHXK4则包含12个疏水通道和10个甘氨酸残基。与其他序列相比，OfHXK4氨基酸序列在其C端的腺苷结合位点处多出11个氨基酸残基(KNEGAS RSKMR)。

图1 桂花OfHXKs氨基酸多序列比对图Figure 1 Multiple sequence alignment of OfHXKs amino acid sequences

2.2 桂花花瓣OfHXKs蛋白质基本理化性质分析及其跨膜区域和二级结构预测

利用ExPASy protparam预测分析OfHXKs家族蛋白质序列的基本理化性质，了解相对分子质量、氨基酸组成、等电点(PI)、不稳定系数、总平均亲水性等，结果见表3。以OfHXK1为例，该基因编码的氨基酸数目为498个，蛋白质相对分子量为53 881.09 Da，理论等电点为5.41，分子式为C2387H3857N653O719S21，氨基酸组成有20种，含有负电荷残基(Asp+Glu)66个，正电荷残基数(Arg+Lys)53个，总平均疏水性为0.055，不稳定系数为47.03，故推测OfHXK1为稳定蛋白质。结果显示：除OfHXK4为不稳定蛋白质外，其余均为稳定蛋白质。蛋白质信号肽预测使用SignalP在线软件默认条件下对OfHXKs氨基酸序列N-端开始前70位氨基酸进行信号肽预测分析，结果显示OfHXKs氨基酸序列均无信号肽存在。应用YLoc在线软件进行亚细胞定位预测，结果显示OfHXK1与OfHXK2均定位于细胞质中，而OfHXK3则定位于叶绿体，OfHXK4位于线粒体。

表3 桂花 OfHXKs 蛋白质基本理化信息Table 3 Information of OfHXK proteins in O. fragrans

使用TMHMM在线软件分析桂花OfHXKs蛋白质可能的跨膜区域，结果显示除OfHXK3外均具有跨膜区域(图2)，OfHXK1、OfHXK2与OfHXK4分别在第7～24位氨基酸残基、第9～28位氨基酸残基和第5～24位氨基酸残基处具有跨膜区域，因此判断OfHXK1、OfHXK2与OfHXK4为跨膜蛋白质，而OfHXK3则不属于跨膜蛋白质。

图2 桂花OfHXKs蛋白质跨膜区段预测Figure 2 Prediction of transmembrane protein segment of OfHXKs

蛋白质的二级结构预测是为了判断目的基因所编码的氨基酸序列残基是处于α-螺旋、β折叠还是其他状态的二级结构。通过SOPMA在线预测分析，结果显示桂花OfHXKs蛋白质二级结构主要以α-螺旋和无规则卷曲为主(表4)。使用SWISS-MODEL进行三级结构预测，结果显示为图3。

表4 桂花花瓣 OfHXKs 蛋白质二级结构预测Table 4 Secondary structures of OfHXK proteins in O. fragrans

图3 桂花OfHXKs蛋白质三级结构预测Figure 3 Prediction of tertiary structures of OfHXK proteins in O. fragrans

2.3 OfHXKs氨基酸序列同源性和系统进化分析

通过NCBI在线数据库比对，将4个桂花花瓣OfHXKs蛋白质序列分别进行同源性比对。Blastp比对结果显示：OfHXK1与油橄榄Olea europaeaXP_022882495.1同源性为96.39%，与芝麻Sesamum indicumXP_022882495.1同源性为88.15%；OfHXK2与油橄榄XP_022843013.1同源性高达97.80%；OfHXK3与油橄榄XP_022884724.1同源性高达96.16%，而与阿月浑子Pistacia veraXP_031283181.1以及粉红钟花Handroanthus impetiginosusPIM98716.1同源性分别为82.70%与85.45%；OfHXK4与油橄榄XP_022873061.1同源性达到95%以上，说明OfHXKs具有高度保守性。

利用MEGAX软件对4个桂花花瓣OfHXKs蛋白质和从NCBI在线数据库中找到的拟南芥HXKs蛋白质序列进行多序列比对，以100% bootstrap支持度构建系统进化树进行分析，预测桂花OfHXKs的功能。由图4所示：所选择蛋白质序列分别为拟南芥AtHXK1(U28214)、AtHXK2(U28215)、AtHXK3(BT030472)、AtHKL1(NM_103929)、AtHKL2(NM_112895)和AtHKL3(NM_119945)。进化树分析结果表明桂花花瓣OfHXK1与OfHXK2关系较近，而OfHXK3和OfHXK4分别聚集在不同小支，由此推断OfHXKs蛋白质之间的功能可能具有一定区别。其中OfHXK1、OfHXK2与拟南芥AtHXK1、AtHXK2聚集在一支，而OfHXK3则与拟南芥AtHXK3关系较近，OfHXK4则与拟南芥AtHKL1、AtHKL2关系较近。

图4 桂花花瓣OfHXKs和拟南芥AtHXKs蛋白质系统进化树分析图Figure 4 Phylogenetic tree of HXK proteins in O. fragrance and A. thaliana

2.4 桂花 OfHXKs基因表达分析

以桂花的顶壳期(S1)、铃梗期(S2)、初花期(S3)和盛花期(S4)花序以及1年生茎、2年生茎、嫩叶、成熟叶的cDNA为模板，以桂花ACT基因为内参，进行实时荧光定量PCR，检测不同品种桂花不同组织及不同发育阶段的OfHXKs基因表达情况。由图5可见：OfHXK1在‘堰虹桂’中的表达量呈现先上升后下降的趋势；在‘金球桂’中同样也在S3时期表达量最高，且显著高于‘堰虹桂’与‘玉玲珑’中相同时期的表达量(P＜0.05)；而OfHXK1基因在‘玉玲珑’中的表达呈现先上升后下降的趋势，于S2时期的表达量达到最高。OfHXK2基因在‘堰虹桂’中的表达量呈现逐步上升的趋势，于盛花期的表达量最高；在‘金球桂’中，OfHXK2的表达量随着开放过程未发生变化；OfHXK2的表达量在‘玉玲珑’中呈现先上升后下降的趋势。OfHXK3基因在3个品种的表达量随着花开放呈现先上升后下降的趋势，均在S2时期达到最高；S2时期‘金球桂’OfHXK3基因的表达量最高，‘玉玲珑’次之。OfHXK4基因在‘堰虹桂’和‘金球桂’中的表达量也是先上升后下降的趋势，在S2时期达到最高，而该基因在‘玉玲珑’中表达量逐渐下降。S2时期‘金球桂’OfHXK4基因的表达量最高，‘堰虹桂’次之。

图5 各品种桂花不同发育时期花序OfHXKs基因相对表达量变化Figure 5 Expression changes of OfHXK genes during different inflorescence of different O. fragrans cultivars

如图6所示：在不同组织中，OfHXK1在嫩叶中的表达量最高，在成熟叶和盛花期花序中表达量次之，在1年生茎中的表达量最低；OfHXK2在盛花期花序中的表达量最高，在2年生茎中的表达量最低；OfHXK3则在茎中表达量最高，在成熟叶中的表达量极低；OfHXK4在成熟叶中的表达量最高，嫩叶中表达量次之，显著高于1、2年生茎和盛开期花序中的表达量。

图6 ‘堰虹桂’不同组织中OfHXKs基因表达量变化Figure 6 Expression changes of OfHXK genes during different tissues of O. fragrans ‘Yanhong Gui’

3 讨论

HXK参与糖代谢的同时也参与己糖信号的感受和传导[20-22]，由于其重要的生物学功能，关于HXK蛋白的研究也逐渐成为热点[9,14,23-25]。本研究在桂花中发现HXK基因的4个同源基因，表明与其他植物一样[6]，HXK基因也是以多基因家族形式存在桂花中，同一HXK基因家族成员序列在不同桂花品种中保守。桂花OfHXK1～OfHXK4基因推导的氨基酸序列均包含葡萄糖和ATP的结合位点，推测OfHXK1～OfHXK4均具有催化己糖磷酸化的功能。其中，OfHXK4氨基酸序列在其N端的腺苷结合位点处多出11个氨基酸残基，这与拟南芥中HKL1和HKL2腺苷结合位点处多10个和6个氨基酸残基的情况[10]类似。

进化树表明OfHXK1、OfHXK2与拟南芥AtHXK1、AtHXK2聚集在一支，已知AtHXK1和AtHXK2具有催化己糖磷酸化的作用[10]，且具有糖感知和传导的功能[26]，推测OfHXK1与OfHXK2同样具有催化己糖磷酸化和感知并转导糖信号的功能。OLSSON等[27]将植物HXK分为A类和B类，A类HXK具有质体信号肽，如AtHXK3；B类HXK具有N端膜锚定结构，如AtHXK1与AtHXK2。根据进化树推测OfHXK1和OfHXK2可能属于B类HXK，此外，核心的糖结合区域(LGFTFSFP-Q-L/I)在OfHXK1和OfHXK2序列中极其保守。OfHXK3与AtHXK3聚集于同一小支，拟南芥AtHXK3具有催化作用但无糖信号传导功能[28]，推测无跨膜区域的OfHXK3属于A类HXK，具有催化作用但无糖信号传导功能。定位于线粒体中的拟南芥AtHKL1和AtHKL2属于酶催化活性不高的HXK亚型，其氨基酸序列C端的腺苷结合位点处多出若干个氨基酸残基[10]，而推测聚在一支的OfHXK4有相同特点。

桂花OfHXK1～OfHXK4基因在1年生茎、2年生茎、嫩叶、成熟叶和盛开时花序中均有表达。其他物种的大部分HXK基因，如拟南芥AtHXK1～AtHXK3、AtHKL1和AtHKL2[10]，番茄Lycopersicon esculentum LeHXK1～LeHXK4[11]， 水 稻Oryza sativa OsHXK2～OsHXK9[8]， 枸 杞Lycium barbarum LbHXK[29]和 苹果Malus×domestica MdHXK1[30]，在被检测的组织器官中均有表达。但是OfHXKs不同成员在不同组织中的表达水平呈现一定差异。OfHXK1和OfHXK2基因在叶和花中的表达量高于茎中，而OfHXK3基因在茎中的表达量最高，表明HXKs基因不同成员虽然磷酸己糖的功能存在冗余[10]，但是不同成员在桂花不同组织中发挥的主导作用存在差异。依据进化树分析，OfHXK1与OfHXK2基因被推测具有糖信号感知和转导功能，可能参与调控桂花花色和花香相关次生代谢物质的生物合成，这可能是OfHXK1与OfHXK2基因在花中的表达量较高的原因之一。

本研究分析了3个桂花品种不同发育阶段花瓣中OfHXKs基因的表达情况，整体上‘金球桂’中OfHXKs基因表达量最高；随着花开放的进程，OfHXK1、OfHXK3和OfHXK4基因的表达量呈现先上升后下降的趋势。枸杞中研究发现，随着枸杞果实的发育，HXK基因表达量在色变期达到峰值，成熟期表达量又下降[29]。牡丹Paeonia suffruticosa花瓣中的研究发现：随着发育过程PsHXK1和PsHXK2基因在花瓣中的表达量呈现先增加后下降的规律，该变化规律与己糖含量有密切的关联[24]。植物组织器官中的糖含量在一定范围内增加可以促进HXK基因转录本和蛋白质的积累，植物体内己糖的磷酸化水平随之增加；超过一定阈值后，植物体内糖含量增加则会抑制HXK基因转录和蛋白质的积累[5]。这也就解释了大部分OfHXKs基因表达量先上升后下降的原因。桂花不同品种色素物质和香气挥发物含量不同[1-2]，推测不同桂花品种对HXK己糖磷酸化和糖信号转导作用的需求不同，此外HXK不同成员还有功能冗余的特点，因此OfHXK2基因在3个桂花品种花瓣中的表达模式不同。

4 结论

本研究根据桂花转录组Unigene序列分析，筛选得到4个HXK基因的同源基因。生物信息学分析结果表明4个基因分别编码461～510个氨基酸残基，同时推测OfHXK1～OfHXK4均具有催化己糖磷酸化的功能。进化树分析表明OfHXK1和OfHXK2与拟南芥AtHXK1与AtHXK2亲缘关系较近；OfHXK3与AtHXK3亲缘关系较近；OfHXK4与AtHKL1和AtHKL2亲缘关系较近。实时荧光定量PCR分析表明桂花OfHXK1～OfHXK4基因在1年生茎、2年生茎、嫩叶、成熟叶和花序中均有表达；OfHXK1、OfHXK3和OfHXK4基因随着花开放的进程，整体上表达量呈现先上升后下降的趋势，而OfHXK2基因在3个品种花序发育过程中表达模式不同。