桃GRAS家族全基因组鉴定与响应UV-B的表达模式分析

李晨，赵雪惠，王庆杰，王旭旭，肖伟，陈修德，付喜玲，李玲，李冬梅

桃家族全基因组鉴定与响应UV-B的表达模式分析

李晨，赵雪惠，王庆杰，王旭旭，肖伟，陈修德，付喜玲，李玲，李冬梅

（山东农业大学园艺科学与工程学院/作物生物学国家重点实验室/山东果蔬优质高效生产协同创新中心，山东泰安 271018）

【目的】基因家族成员在调节植物生长发育中发挥着关键作用。通过生物信息学分析在桃基因组中的分布、结构及进化，研究家族成员在不同组织中的表达特异性及其对UV-B的响应，解析基因家族的生物学功能。【方法】对设施油桃‘中油5号’（var.cv. Zhongyou5）补充适量UV-B（Ultraviolet-B）剂量，利用Plant TFDB数据库鉴定桃基因家族成员。采用Clustal W、MEGA6.0、ProtParam tool、MCScanX、Circos、SMART、NCBI-CDD、ExPASy、GSDS和MEME等软件构建系统进化树，绘制染色体定位图，预测蛋白的相对分子质量与等电点等理化性质等，分析基因家族成员在不同组织中的表达模式，利用qRT-PCR技术检测桃在UV-B处理下的表达情况。【结果】从桃全基因组中鉴定出48个GRAS转录因子家族基因，构建系统进化树将这48个成员分为9个亚家族，在桃的8条染色体上呈不均匀分布。对基因家族进行理化性质分析发现，其蛋白平均长度为590.52 aa，等电点在4.36—7.56。家族基因结构分析表明，有40个基因不含内含子，8个基因含有1个内含子。保守元件分析显示，家族包含20个保守元件，其中Motif 2和4在的家族中高度保守，同一个亚家族成员含有相同的保守元件，其可能具备相似的功能。然而，有些亚家族成员的表达模式不同，这可能与其保守基序之外的序列有关。在不同组织中具有不同的表达模式。叶片中经UV-B处理后上调表达最显著，而有多达15个基因表达量呈现下调。果实中经UV-B处理后上调表达，但有9个基因表达下调。韧皮部中，UV-B处理后有14个基因上调表达，而经处理后在韧皮部中表达下调最明显。【结论】从桃基因组中共鉴定出48个GRAS基因家族成员，分布于8条染色体上；多数能响应UV-B处理，但在不同组织中的表达不尽相同。

桃；基因家族；生物信息学分析；UV-B

0 引言

【研究意义】桃是中国北方的主栽果树之一，随着果树集约化栽培的推行，设施桃栽培得到了迅速发展。然而与露天栽培相比，设施栽培果树表现出果实甜度变低、风味偏淡、综合品质下降等问题，棚膜吸收、反射太阳光导致设施内UV-B弱化是重要的影响因素之一[1]。基因家族成员在植物生长发育的调节中发挥着关键作用，包括光信号转导、GA信号调节和腋芽分生组织形成等[2]。对桃基因家族进行生物信息学分析、组织特异性表达及响应UV-B处理的表达模式分析，有利于深入了解家族基因的潜在功能，并对桃响应UV-B光信号的基因挖掘具有重要意义。【前人研究进展】自20世纪80年代后期，人类活动造成了臭氧层空洞，而臭氧量的减少导致到达地表的UV-B辐射增强。UV-B是紫外线中的中波辐射（280—320 nm），约占太阳辐射总量的1.5%。其能量高，且能被核酸、蛋白质及脂类等物质吸收。研究表明UV-B具有损伤DNA、改变蛋白质结构；破坏光合作用；造成植株矮化、株型缩小，生长受抑制等负作用，对生态系统和动、植物生长产生深远影响[1,3-4]。事实上，UV-B对植物的影响效应具有双重性，高UV-B辐射对植物生长产生不良影响，低剂量的UV-B在基因表达、代谢活动以及生长发育等方面反而会起到正面的调节作用[5-6]。UV-B辐射是应激源还是调节器与其辐射强度和时间有关[7]。设施条件下的UV-B强度较露天弱，叶片中紫外吸收物质较露天少，补充UV-B后含量总体得到提升[8]。比较紫外线敏感品种‘旭日’和不敏感品种‘朝阳’叶片中的类黄酮含量发现，经较高强度UV-B辐射后，前者类黄酮含量先升高后降低，变化明显，后者则相对稳定，表明朝阳更能有效适应UV-B强度的增加[7]。光受体UV RESISTANCE LOCUS8（UVR8）被发现后，UV-B由环境胁迫因子变成了特殊的信号调控因子[1]，其在光感受中具有独特的调节机制[9-10]。在UV-B处理后，UVR8能与下游多种转录因子相互作用，从而启动光形态发生[11]。长时间以较低UV-B剂量照射作物，并未发现ROS积累及其编码基因的表达，反而会启动UVR8信号通路，引起良性应激，但UV-B剂量超出作物承受范围时会引起破坏性应激[12]。更有研究发现大多数植物在UV-B辐射下表现出植株形态、生理生化和激素等各个方面的变化[1]。转录因子在植物生长发育过程中发挥尤为重要的作用。基因家族通常被认为是植物特异性转录调节因子家族[13]，能够响应激素、光等多种信号来调节植物的生长发育。GRAS蛋白质由可变的N末端和高度保守的C末端组成，其被称为GRAS结构域。该结构域主要由两个亮氨酸七肽重复序列（LHRI和LHRII）组成的VHIID、PFYRE和SAW共3个基序组成[14-15]。在拟南芥中，将GRAS家族蛋白分为8个亚家族：DELLA、LS、SCR、SHR、PAT1、HAM、SCL9（LISCL）和SCL4/7[14]。而在水稻中[15]也认为可能会分为8个亚科，但包括一个新的亚家族SCL3，而不是SCL4/7。葡萄中GRAS家族蛋白序列也分成了8个亚家族，且其氨基酸序列主要为-螺旋和随机卷曲，其三维结构相似[16]。拟南芥、水稻和其他植物种类的GRAS蛋白序列构建系统进化树分析发现，其具有12个独立的分支[17]，这些研究表明该基因家族在植物中表现出多样化。其中，DELLA亚家族成员参与了拟南芥中的赤霉酸信号传导[18]，葡萄中/可能响应GA负调控参与无核果实的发育[19]，亚家族成员在番茄的腋芽分生组织中起作用[20]，亚家族成员、维持了拟南芥顶端分生组织的属性与更新[17,21]，分支中的和作为光敏色素信号通路的中间介质发挥作用[22]，亚家族成员对促进拟南芥的根细胞伸长起作用[23]。Ma等[24]证明在胡杨中过量表达（亚家族）增强了转基因拟南芥的抗旱性和耐盐性。【本研究切入点】目前在模式植物拟南芥中，对基因家族的结构、进化、相互作用机理以及抵御逆境等方面的研究较为深入。然而，在木本植物特别是果树中该家族基因的系统报道较少，桃基因家族成员响应UV-B光信号的功能有待进一步研究。【拟解决的关键问题】以设施栽培的‘中油5号’油桃的叶片、果实和韧皮部为研究材料，对家族中的48个成员进行生物信息学分析，旨在了解其基本理化性质、家族分类、保守结构域以及基因结构等；通过检测基因家族成员在UV-B处理后的叶片、果实和韧皮部中的表达水平，了解家族的基因性质，筛选该家族中可能响应UV-B光信号的重要基因。

1 材料与方法

试验于2018年在山东农业大学园艺科学与工程学院进行。

1.1 试验材料

试材为7年生设施油桃‘中油5号’（var.cv. Zhongyou5）。自开花至果实成熟期间，用1.44 kJ·m-2·d-1UV-B照射处理植株（此剂量由前期试验筛选得出）。UV-B由专用LED UV-B灯（40 W，297 nm，南京Kazhi）提供，悬挂在树体上方1.5 m处，通过电子自动控制装置控制紫外灯的开/关时间，每天上午9:30—10:30补充1 h，阴雨天和下雪天停止照射。使用配备有297 nm探针的UV-B型单通道UV辐照器（北京师范大学光电仪器厂）测定灯下80—120 cm范围内的UV-B辐射剂量为1.44 kJ·m-2·d-1，对该范围内的功能叶（以树体外围一年生枝的生长点为起始点，倒数第6—8片叶）、果实和一年生枝韧皮部进行取样，以相同设施条件下未进行UV-B照射的‘中油5号’为对照，3次重复。液氮速冻后-80℃保存备用。

1.2 数据来源与分析方法

1.2.1 桃的获取、理化性质预测与基因染色体定位分析 从Plant TFDB（http://planttfdb.cbi.pku.edu. cn/）下载桃家族数据，利用SMART（http://smart. embl-heidelberg.de/）和NCBI-CDD工具进行蛋白结构预测，检测候选蛋白的结构域，去除不含GRAS结构域的序列，筛选候选基因。利用ProtParam tool（http://web.expasy.org/protparam/）对所有GRAS家族氨基酸序列进行分子量和等电点预测。根据染色体的位置依次命名，利用MCScanX（http://chibba.pgml.uga. edu/mcscan2/）[25]和Circos软件（http://circos.ca/）[26]绘制桃基因家族成员在染色体上的定位分布。2个或多个同源基因在染色体上的物理位置不超过100 kb即视为串联重复基因[27]。

1.2.2 系统进化树构建、结构域分析及基因结构分析 参照Nakano等[28]的分类方法，将拟南芥与桃家族成员的蛋白序列用Clustal W进行多重序列比对后导入进化树分析软件MEGA 6.0，采用相邻连接法（neighbor-joining NJ；执行参数：Bootstrap method 1000，Poisson model，Pairwise deletion）构建系统进化树，对桃GRAS蛋白进行亚族分类。利用MEME（http://me-me-suite.org/index.html）对桃保守域结构进行分析。利用GSDS online（http://gsds.cbi. pku.edu.cn/）预测基因结构。

1.3 总RNA的提取和qRT-PCR

使用RNAprep Pure多糖多酚植物总RNA提取试剂盒（天根，Cat No.DP441）提取液氮冻样的RNA，利用HiScript®II Q RT SuperMix for qPCR（+gDNA wiper）试剂盒（Vazyme公司，中国）合成cDNA第一链。

采用实时荧光定量PCR的方法，用cDNA模板检测各基因表达水平，利用Primer 3.0设计荧光定量引物，引物序列见表1，内参为桃-Actin，由生工生物工程（上海）股份有限公司合成。使用ChamQTMUniversal SYBR®qPCR Master Mix（Vazyme）进行实时荧光定量PCR。qRT-PCR反应体系为：10.0 μL 2×ChamQ Universal SYBR qPCR Master Mix，上、下游引物分别为0.4 μL（10 μmol·L-1），1.0 μL模板，8.2 μL ddH2O。每个样品进行3次技术重复。反应条件为：95℃预变性30 s；95℃ 5 s，60℃ 30 s，40个循环，绘制溶解曲线，95℃ 15 s，60℃ 60 s，95℃ 15 s。用2-∆∆CT法对试验数据进行定量分析。

表1 48个PpGRAS的qRT-PCR引物序列

1.4 数据分析

采用Microsoft Excel 2003和SPSS 19.0软件进行数据处理和统计分析，包括单因素方差分析和差异显著性分析，应用GraphPad Prism 6.01软件绘图。

2 结果

2.1 桃GRAS转录因子家族成员与染色体分布分析

从桃基因组序列中鉴定得到48个转录因子家族成员，根据基因在染色体上的位置，对所有转录因子进行系统编号。开放式阅读框长度在1 218—2 550 bp，蛋白长度在406—850 aa，平均长度590.52 aa，分子质量在45 828.1—78 010.1 Da，等电点在4.36—7.56，保守结构域分析结果表明，鉴定出的48个PpGRAS蛋白均具有GRAS特征结构域。桃转录因子的详细信息见表2。

依据的染色体位置信息，利用Circos软件绘制其染色体定位图。由图1可知，中有2个基因（和）无匹配的染色体定位信息。染色体定位分析发现，在桃的8条染色体上呈现不均匀分布，其中第6号染色体上分布最多，为14个；其次是第3号染色体，为7个；第5和7号染色体分布最少，均为3个。由于基因重复能显著促进基因家族的扩展和蛋白质功能的多样化，本研究分析了该基因家族的串联重复。结果显示7对共线GRAS基因（/；；；；；；）和10个串联GRAS基因（；），占总数的20.83%。共线基因是能通过复制在另一基因组以相同连续顺序存在的旁系同源基因，除和外，所有共线都在同一亚家族中。

图1 桃GRAS家族基因染色体定位及其共线性预测

表2 桃GRAS及其相关信息

2.2 桃GRAS蛋白序列与拟南芥系统进化树的对比分析

通过与拟南芥蛋白聚类分析，将桃中48个GRAS蛋白分为9个亚家族（图2），这些亚家族是根据早期研究分类[15]或他们中的一名成员在新发现的亚家族中命名，分别为PAT1、SHR、HAM、LISCL、SCR、LS、DELLA、SCL3、U。每个亚家族的GRAS成员可能具有相似或相关的功能。一般而言，大多数亚家族都拥有来自2个物种的GRAS成员（图2）。然而，在桃中新鉴定出来的U亚家族不包括任何拟南芥成员，表明拟南芥中谱系特异性基因缺失。LISCL是最大的子家族，包含12个成员；而在拟南芥中LISCL也是最大的子家族，包含7个成员；其次是SCL3亚家族，有9个成员。PpGRAS家族成员共形成8个旁系同源基因对，其中3对的自展值为100，而且Prupe.6G073900.1.p/Prupe.6G074000.1.p的自展值为99，U家族仅有的两个成员Prupe.3G309800.1.p/Prupe. 8G054000.1.p的自展值达87。此外，PpGRAS与AtGRAS形成了14个直系同源基因对，其中12对的自展值为100（图2）。

图2 桃与拟南芥GRAS蛋白序列系统进化树比对

2.3 桃GRAS家族保守元件和基因结构分析

基因结构多样性作为进化的可能基础，可以为多基因家族的进化提供有价值的信息[15,28]。基于每个预测基因的编码序列（CDS）和基因组DNA序列分析基因结构（图3）。桃中家族基因的内含子数量最多为1，其中40个基因没有内含子，其余8个基因（Prupe.2G138400.1.p、Prupe.8G220900.1.p、Prupe. 6G300900.1.p、Prupe.6G073900.1.p、Prupe. 3G040400.1.p、Prupe.8G016200.1.p、Prupe.4G233000.1.p、Prupe. 3G309800.1.p）仅有1个内含子。通常，来自桃中相同亚家族（SCL3、DELLA和LS）的具有相同的结构。

图3 桃GRAS结构

为揭示桃中GRAS蛋白质的多样性，进一步确定桃家族成员间的关系，采用MEME预测桃家族基因中的元件。从候选的桃GRAS蛋白质中共鉴定出20个保守元件（图4），并列出了这些保守元件的结构域序列标识（即氨基酸组成）（图5）。本研究发现，所有的桃GRAS蛋白都含有Motif 4，表明其在的家族中高度保守；此外，高度保守的Motif 2同样存在于桃所有GRAS蛋白中。其中motif 2、11在SCL3亚家族成员的N端高度保守，而DELLA亚家族的3个成员的保守基序表现出较好的一致性，共有10个保守元件（motif 2、3、4、7、9、15、17、18、19和20），LS亚家族共有8个保守元件（motif 2、3、4、7、11、15、17和19），HAM亚家族中有5个共有保守元件（motif 2、4、5、9和19），SCR亚家族中共有7个保守元件（motif 2、3、4、5、7、9和17），LISCL家族中具有特异性保守元件Motif 13和在其N端高度保守的Motif 16，SHR亚家族中仅有3个共有保守元件（motif 4、5和17），PAT1亚家族中共有9个保守元件（motif 2、3、4、5、6、7、9、17和19），U亚家族共有7个保守元件（motif 2、3、4、5、9、17和19）。由此发现系统发育树中密切相关的含有相同的保守元件，而且有的亚家族中具备专有的保守元件，表明同一个亚家族的基因具备相似的功能。从桃的基因结构和保守元件位置分析，确定大多数成员在各个亚家族中是保守的。

图4 桃GRAS基因家族保守元件

图5 桃GRAS结构域的序列标识

2.4 桃GRAS家族基因的组织特异性表达分析及其对UV-B的响应

家族基因在叶片、果实和韧皮部中的组织特异性分析如图6-A所示。48个在叶片、果实和韧皮部3种组织中均有表达，表达模式存在一定的差异。33.33%的基因（16个）在叶片中的表达量高于其余两种组织。、和在桃韧皮部中高度表达；、、和在果实中的表达最高。UV-B处理后48个基因在3个组织中也均有表达，其表达模式与CK不尽相同。其中，和在叶片中的表达量高于其他组织；和在桃韧皮部中高表达；而有25个基因在果实中的表达量低于其他组织。并且，有8个基因（和）在3个组织中的表达水平基本一致。

响应UV-B信号的在不同组织内的表达量与CK的比值见图6-B。在叶片中，和在UV-B处理后上调表达，而31.25%的基因（15个）表达量下调。在果实中，经UV-B处理后上调表达，而有9个基因（和）下调表达。在韧皮部中，UV-B处理后有29.17%的基因（14个）上调表达，而和经处理后表达下调。另外，有7个基因（和）经UV-B处理后在各个组织内均未表现出明显变化。

3 讨论

3.1 PpGRAS转录因子家族的鉴定和分析

是涉及植物生长和发育等各方面重要的转录因子。基因家族在植物中广泛分布，且家族成员的数量与基因组的大小无关[21]。植物体在生物进化过程中主要通过基因组的广泛复制和多样化，产生基因家族中的多成员[29-30]。在拟南芥中，家族发生了2次串联重复事件[2]。而且，基因家族的串联重复现象在番茄、水稻、白杨和辣椒等物种中均被检测到[2,21,31]。本研究鉴定出48个基因家族成员，低于水稻（57个）[15]、辣椒（58个）[21]和番茄（53个）[32]等作物，可能是由于桃基因组中没有大规模的片段复制，且相对保守等原因造成的。通过染色体定位分析发现，在桃基因组中广泛分布，且存在明显的串联重复现象，串联重复基因占20.833%（10/48），这表明串联重复在基因家族扩增和进化中起着重要作用。

研究发现，水稻中的一些GRAS蛋白质并未找到合适的亚家族[2]，本研究发现，在桃家族中同样有两个基因（）在系统进化树分析中不属于已知的8个亚家族，而且其可靠性达87，将其划分为U亚家族，因而将桃中共划分为9个亚家族[15]。该U亚家族成员基因的表达有所差异，其中在果实中表达较低，而在果实和韧皮部中表达都较低；UV-B处理后的表达基本不受影响，但的表达在韧皮部中上调。而DELLA、LISCL和PAT1亚家族的基因结构具有明显的一致性，表明这些亚家族的基因具备更为相似的功能。而且还发现LISCL和PAT1亚家族的基因表达模式同样具有明显的特异性，主要在叶中表达。但是有些亚家族内基因的表达模式有所不同，可能参与了不同的生长发育过程。推测同一亚家族之间不同的表达差异可能与保守基序外的序列有关。此外，基因表达模式的研究可为家族基因的功能验证提供依据。基因家族的系统发育分析为进一步研究桃的功能基因组提供了理论基础。

3.2 PpGRAS响应UV-B处理的表达分析

GRAS蛋白在多种植物中被发现，功能呈多样性，并在不同品种、不同组织和不同生长阶段所发挥的作用不尽相同[2,32]。目前，模式植物拟南芥的基因家族中的多个亚家族成员的功能已被验证[15]。本研究发现，设施桃叶片、果实和韧皮部3个组织中均有多个响应UV-B处理，响应情况不尽相同。

3.2.1桃叶片中响应UV-B光信号的表达分析 DELLA蛋白作为负调控因子参与了GA信号的转导过程。当DELLA蛋白接收GA信号时，会迅速降解，植株表现正常的GA响应程序[23]，而突变体植株却表现出GA不敏感表型[33]。本研究发现，作为DELLA亚家族成员中同源基因，在UV-B处理后，在叶片中上调表达最为显著。有研究表明，在桃上进行UV-B处理后造成了叶片中内源赤霉素含量降低[34]。由此推测，UV-B处理后可能导致叶片中内源赤霉素含量降低，又经由GA信号调控了，从而使其上调表达。

矮牵牛中主要在侧生器官原基和茎原维管组织中表达，对维持茎尖分生组织活性发挥了重要作用。矮牵牛突变体植株的叶片数明显减少，茎尖分生组织丧失分化特性，阻止了器官的形成[35]。而本研究发现，桃HAM亚家族成员经UV-B处理后在叶片中表达下调最为明显，推测其可能响应UV-B处理后对叶片的叶原基分化起着重要作用，具体功能尚需进一步验证。

3.2.2桃果实中响应UV-B光信号的表达分析 远红光、红光和白光配合处理都降低了的转录水平，而且与SCL21相互作用的PAT1蛋白主要参与了phy A介导的光信号转导过程，在光信号转导中起着重要作用[22]。本研究发现，作为PAT1亚家族中的同源基因，经UV-B处理后在果实中显著下调表达，表明基因家族成员在响应UV-B光信号中也可能发挥了重要作用。

作为桃SCR亚家族的成员，在桃果实中经UV-B处理后上调表达最为显著。SHR和SCR蛋白的作用会导致细胞分裂过程的缺陷[36]，相似结果也在水稻[37]和玉米[38]中发现；而且不仅参与了根的发育调控，同时也参与了叶片保卫细胞的形成[39]。陈丽华等[40]研究发现，经GA处理后的‘无核白鸡心’葡萄的果形明显拉长，而且在GA处理后呈现整体性的上调表达，可能在响应GA处理后调控了葡萄果实细胞的分裂和分化。因而，推测可能具有响应UV-B光信号调控桃果实纵、横径的功能。

3.2.3桃韧皮部中响应UV-B光信号的表达分析 在拟南芥根部，SCL3能作为GAI和RGA的衰减子调控伸长区和分化区细胞的伸长，SCL3-DELLA互作调节的GA信号途径，能与SHR/SCR途径联合调控分裂区基本组织的分裂时间和程度[23]。相似的结论同样出现在毛竹中[39]。本研究发现，SCL3亚家族成员在韧皮部中经UV-B处理后下调表达最为明显，推测可能负响应UV-B光信号，减弱对桃韧皮部细胞分裂程度的抑制。

果实发育过程中，如果基因在幼果期的表达水平高于成熟期，则其可能在果实前期发育过程中发挥重要功能。黄伟[32]研究发现，番茄中的表达在果实转色期有显著上升，其可能对启动果实成熟具有重要作用。而本研究发现，作为的同源基因，经UV-B处理后在韧皮部中上调表达最为明显，说明其可能在“源-库”间有机物的运输中起着重要作用。

4 结论

通过生物信息学分析预测出48个桃基因家族成员，分为9个亚家族，基因结构、染色体分布和保守元件分布等多种生信分析表明这些基因在进化过程中比较保守。采用qRT-PCR技术对桃基因家族成员进行组织特异性分析，发现48个基因在叶片、果实和韧皮部中均有表达，但表达变化不尽相同。本研究还探索了基因家族成员对UV-B的响应模式，发现多数能响应UV-B处理，表明基因家族可能在光响应方面发挥着多种潜在功能。

[1] 李冬梅, 李少旋, 徐功勋, 李晨, 付喜玲, 陈修德, 张海森, 高东升. 设施作物响应UV-B辐射的研究进展. 植物生理学报, 2018, 54(1): 36-44.

Li D M, Li S X, Xu G X, Li C, FU X L, Chen X D, Zhang H S, Gao D S. Research advances of plant response to UV-B radiation in greenhouse., 2018, 54(1): 36-44. (in Chinese)

[2] Liu X Y, Widmer A. Genome-wide comparative analysis of the GRAS gene family in populus,and rice., 2014, 32: 1129-1145.

[3] Gill S S, Tuteja N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants., 2010, 48(12): 909-930.

[4] Lidon F C, Teixeira M, Ramalho J C. Decay of the chloroplast pool of ascorbate switches on the oxidative burst in UV-B-irradiated rice., 2012, 198(2): 130-144.

[5] Li F R, Peng S L, Chen B M, Hou Y P. A meta-analysis of the responses of woody and herbaceous plants to elevated ultraviolet-B radiation., 2010, 36(1): 1-9.

[6] Ballare C L, Caldwell M M, Flint S D, Robinson S A, Bornman J F. Effects of solar ultraviolet radiation on terrestrial ecosystems. patterns, mechanisms, and interactions with climate change., 2011, 10(2): 226-241.

[7] Frohnmeyer H, Staiger D. Ultraviolet-B radiation-mediated responses in plants, balancing damage and protection., 2003, 133(4): 1420-1428.

[8] 郭玉才. 不同棚膜及人工补充紫外线B对设施桃树形态建造和叶片生理的影响[D]. 泰安: 山东农业大学, 2009.

Guo Y C. Effects of film and UV-B on morphology and leaf physiology of peach tree in greenhouse [D]. Tai’an: Shandong Agricultural University, 2009. (in Chinese)

[9] Christie J M, Arvai A S, Baxter K J, Heilmann M, Pratt A J, O’hara A, Kelly S M, Smith B O, HITOMI K, JENKINS G I, GETZOFF E D. Plant UVR8 photoreceptor senses UVB by tryptophan-mediated disruption of cross-dimer salt bridges., 2012, 335(6075): 1492-1496.

[10] Wu D, Hu Q, Yan Z, Chen W, Yan C Y, HUANG X, ZHANG J, YANG P Y, DENG H T, WANG J W, DENG X W, SHI Y G. Structural basis of ultraviolet-B perception by UVR8., 2012, 484(7393): 214-219.

[11] Osterlund M T, Hardtke C S, Wei N, Deng X W. Targeted destabilization of HY5 during light-regulated development of., 2000, 405(6785): 462-466.

[12] Hideg E, Jansen M A K, Strid A. UV-B exposure, ROS, and stress: Inseparable companions or loosely linked associates?, 2013, 18(2): 107-115.

[13] Sun X l, Jones W T, Rikkerink E H a. GRAS proteins: The versatile roles of intrinsically disordered proteins in plant signalling., 2012, 442(1): 1-12.

[14] Bolle C. The role of GRAS proteins in plant signal transduction and development., 2004, 218(5): 683-692.

[15] Tian C g, Wan P, Sun S h, Li J y, Chen M s. Genome-wide analysis of the GRAS gene family in rice and., 2004, 54(4): 519-532.

[16] 孙欣, 王晨, 房经贵, 慕茜, 王西成. 葡萄GRAS基因家族生物信息学分析. 江西农业学报, 2011, 23(7):1-8.

Sun X, Wang C, Fang J G, Mu Q, Wang X C. Bioinformatics analysis of GRAS gene family in grapevine., 2011, 23(7): 1-8. (in Chinese)

[17] Engstrom E M. Phylogenetic analysis of GRAS proteins from moss, lycophyte and vascular plant lineages reveals that GRAS genes arose and underwent substantial diversification in the ancestral lineage common to bryophytes and vascular plants., 2011, 6(6): 850-854.

[18] Peng J R, Carol P, Richards D E, King K E, Cowling R J, Murphy G P, Harberd N P. TheGAI gene defines a signaling pathway that negatively regulates gibberellin responses., 2017, 11(23): 3194-3205.

[19] 张文颖, 朱旭东, 崔腾飞, 贾海峰, 房经贵, 王晨. 葡萄VvRGL基因应答GA调控无核果实发育的研究. 西北植物学报, 2019, 39(3): 381-390.

Zhang W Y, Zhu X D, Cui T F, Jia H F, Fang J G, Wang C. VvRGL genes regulate seedless fruit development by responding to GA in grapevine., 2019, 39(3):381-390. (in Chinese)

[20] Schumacher K, Schmitt T, Rossberg M, Schmitz G, Theres K. The Lateral suppressor (Ls) gene of tomato encodes a new member of the VHIID protein family., 1999, 96(1): 290-295.

[21] 张焕欣, 董春娟, 尚庆茂. 辣椒GRAS家族全基因组鉴定与表达分析. 园艺学报, 2017, 44(12): 2305-2317.

Zhang H X, Dong C J, Shang Q M. Genome-wide identification and expression analysis of GRAS gene family in pepper., 2017, 44(12): 2305-2317. (in Chinese)

[22] 周莲洁, 张富春, 王艳. GRAS家族基因在植物生长、代谢及逆境胁迫中的功能研究进展. 植物生理学报, 2013, 49(9): 855-860.

Zhou L J, Zhang F C, Wang Y. Research progress on the functional mechanism of GRAS family genes in plant growth, metabolism and stress., 2013, 49(9): 855-860. (in Chinese)

[23] Heo J O, Chang K S, Kim I A, Lee M H, Lee S A, Song S K, Lee M M, Lim J. Funneling of gibberellin signaling by the GRAS transcription regulator SCARECROW-LIKE 3 in theroot., 2011, 108(5): 2166-2171.

[24] Ma H S, Liang D, Shuai P, Xia X L, Yin W L. The salt- and drought-inducible poplar GRAS protein SCL7 confers salt and drought tolerance in., 2010, 61(14): 4011-4019.

[25] Wang Y P, Tang H B, Debarry J D, Tan X, Li J P, Wang X Y, Lee T H, Jin H Z, Marler B, Guo H, Kissinger J C, Paterson A H. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity., 2012, 40(7): e49.

[26] Krzywinski M, Schein J, BIROL I. Circos: An information aesthetic for comparative genomics., 2009, 19(9): 1639-1645.

[27] Huang S X, Gao Y F, Liu J K, Peng X L, Niu X L, Fei Z J, Cao S Q, Liu Y S. Genome-wide analysis of WRKY transcription factors in., 2012, 287(6): 495-513.

[28] Nakano T, Suzuki K, Fujimura T, Shinshi H. Genome-wide analysis of the ERF gene family inand rice., 2006, 140(2):411-432.

[29] 孙明岳, 周君, 谭秋平, 付喜玲, 陈修德, 李玲, 高东升. 苹果bZIP转录因子家族生物信息学分析及其在休眠芽中的表达. 中国农业科学, 2016, 49(7): 1325-1345.

Sun M Y, Zhou J, Tan Q P, Fu X L, Chen X D, Li L, Gao D S. Analysis of basic leucine zipper genes and their expression during bud dormancy in apple (×)., 2016, 49(7): 1325-1345. (in Chinese)

[30] Song J, Gao Z H, Huo X M, Sun H L, Xu Y S, Shi T, Ni Z J. Genome-wide identification of the auxin response factor (ARF) gene family and expression analysis of its role associated with pistil development in Japanese apricot (Sieb. et Zucc)., 2015, 37: 145.

[31] Huang W, Xian Z Q, Kang X, Tang N, Li Z G. Genome-wide identification, phylogeny and expression analysis of GRAS gene family in tomato., 2015, 15(1): 209.

[32] 黄伟. 番茄GRAS转录因子家族的鉴定及SlGRAS24基因的功能研究[D]. 重庆: 重庆大学, 2016.

Huang W. Comprehensive genome-wide analysis of the GRAS transcription factor family in tomato and functional identification of SlGRAS24[D]. Chongqing: Chongqing University, 2016. (in Chinese)

[33] Acharda P, Vriezen W H, van der Straeten D, Harberd N P. Ethylene regulates Arabidopsis development via the modulation of ELLA protein growth repressor function., 2003, 15(12): 2816-2825.

[34] 吴杏春. 水稻对UV-B辐射增强的生理生化响应及其适应机制研究[D]. 福州: 福建农林大学, 2001.

WU X C. Studies on physiochemical response and its adaptive mechanism of rice (L.) exposed to enhanced Ultraviolet- B radiation [D]. Fuzhou: Fujian Agriculture and Forestry University, 2001. (in Chinese)

[35] Huang W, Peng S Y, Xian Z Q, Lin D B, Hu G J, Yang L, Ren M Z, Li Z G. Overexpression of a tomato miR171 target gene, SlGRAS24 impacts multiple agronomical traits via regulating gibberellin and auxin homeostasis., 2016, 15(4): 472-488.

[36] Koizumi K, Wu S, MacRae-Crerar A, Gallagher K L. An essential protein that interacts with endosomes and promotes movement of the SHORT-ROOT transcription factor.2011, 21(18): 1559-1564.

[37] Kamiya N, Itoh J I, Morikami A, Nagato Y, Matsuoka M. The SCARECROW gene’s role in asymmetric cell divisions in rice plants., 2003, 36(1): 45-54.

[38] Lim J, Helariutta Y, Specht C D, Jung J, Sims L, Bruce W B, Diehn S, Benfey P N. Molecular analysis of the SCARECROW gene in maize reveals a common basis for radial patterning in diverse meristems., 2000, 12(8): 1307-1318.

[39] 董丽莉, 赵韩生, 李利超, 孙化雨, 王丽丽, 高志民. 毛竹PeSCL3基因表达特征及其启动子活性研究. 热带亚热带植物学报, 2016, 24(3): 252-258.

Dong L L, Zhao H S, Li L C, Sun H Y, Wang L L, Gao Z M. Expression characteristics and promoter activity analysis ofgene from phyllostachys edulis., 2016, 24(3): 252-258. (in Chinese)

[40] 陈丽华, 葛晖, 柴丽娟, 陈尚武, 马会勤. 葡萄SCARECROW基因家族的分析与表达. 中国农业大学学报, 2012, 17(1): 80-87.

Chen L H, Ge H, Chai L J, Chen S W, Ma H Q. Putative SCARECROW gene family ofL. and their differential expression under GA treatment., 2012, 17(1): 80-87. (in Chinese)

Genome Identification ofFamily and Expression Pattern Analysis of Responding to UV-B in Peach

LI Chen, ZHAO XueHui, WANG QingJie, WANG XuXu, XIAO Wei, CHEN XiuDe, Fu XiLing, LI Ling, LI DongMei

(College of Horticulture Science and Engineering, Shandong Agricultural University/State Key Laboratory of Crop Biology/Shandong Collaborative Innovation Center for Fruit & Vegetable Production with High Quality and Efficiency, Taian 271018, Shandong)

【Objective】transcription factor family genes play a key role in the regulation of plant growth and development. The objectives of this study were to analyze the distribution, structure and evolution ofin the peach genome by bioinformatics, to study the expression specificity of family members in different tissues and their responses to UV-B optical signal, and to investigate the biological function oftranscription factor family genes in peach. 【Method】The facility nectarine ‘var.cv. Zhongyou5 was supplemented with an appropriate dose of UV-B (Ultraviolet-B).The peachgene in the peach genome was identified by using the Plant TFDB database.Phylogenetic tree, chromosome localization, relative mass and isoelectric point and other physical and chemical properties ofmember were analyzed with Clustal W, MEGA6.0, ProtParam tool, MCScanX, Circos, SMART, NCBI-CDD, ExPASy, GSDS2.0, and MEME, respectively. The expression pattern ofgene family in different tissues was analyzed, and the expression of some members ofgene family in peach treated with UV-B was detected by real-time fluorescence quantitative PCR (qRT-PCR).【Result】48 members oftranscription factor family were identified from the whole genome of the peach, and they could be divided into 9 categories by constructing a phylogenetic tree. Thegene showed uneven distribution on 8 chromosomes of peach. The theoretical isoelectric point of the family protein was ranged from 4.36 to 7.56, and the average number of amino acids encoded was 590.52. Gene’s structure analysis showed that 40 genes contained no introns, and 8 genes contained 1 intron. Conservative elemental analysis revealed that thefamily contains 20 conserved elements, of which Motif 2 and Motif 4 were highly conserved in thefamily. Members of the same subfamily contained the same conserved elements, suggesting that members of the same subfamily might have similar functions. However, some subfamily members had different expression patterns, which might be related to sequences other than the conserved motif. Thegenes had different expression patterns in different tissues; and most ofgenes could respond to UV-B treatment, but the expression changes were different in different tissues. In leaves,was up-regulated after UV-B treatment, while up to 15 genes were down-regulated. In the fruit,was up-regulated by UV-B treatment, but 9 genes were down-regulated. In the phloem, 14 genes were up-regulated after UV-B treatment, whilewas most down-regulated in the phloem after treatment. 【Conclusion】A total of 48gene family members were identified from the peach genome and distributed on 8 chromosomes; mostgenes responded to UV-B treatment, but the expression changes were different in different tissues. This study laid the foundation for further analysis of thefamily of genes in response to UV-B light signals and other potential functions.

peach;gene family; bioinformatic analysis; UV-B

2019-06-11；

2019-08-03

国家自然科学基金（31601706）、山东省自然科学基金（ZR2016CM09）

李晨，Tel：15650453756；E-mail：15650453756@163.com。通信222作者李冬梅，Tel：0538-8246263；E-mail：dmli2002@sdau.edu.cn

（责任编辑 赵伶俐）