杜仲CONSTANS-like全基因组鉴定、系统进化及表达模式分析

  目的  揭示CONSTANS-like在杜仲Eucommia ulmoides基因组中的分布、结构特征及表达模式。  方法  利用生物信息学方法，对杜仲CONSTANS-like基因家族进行鉴定及理化性质、进化关系、基因结构、启动子元件和表达模式分析。  结果  杜仲基因组中共鉴定到8个EuCOLs基因，分别命名为EuCOL1~EuCOL8，氨基酸数目为315~469，理论等电点分布范围为5.10~6.47，分子量为35.21~52.65 kDa。亚细胞定位预测均定位在细胞核中，为亲水性蛋白，分布于8条染色体。系统进化分为2个亚家族(群组Ⅰ和群组 Ⅲ)，分别包含2和6个EuCOLs蛋白，同一亚家族基序具有相似性。EuCOLs基因结构简单，启动子中含有多个光周期响应元件。表达模式分析显示:EuCOLs在杜仲叶片发育中表达水平相对较低，EuCOL7在杜仲胶形成中表达量最高，各家族成员表达特征存在差异。蛋白互作预测显示:EuCOL7可与多个光周期响应蛋白互作。  结论  杜仲CONSTANS-like基因家族含有典型的CCT和B-box结构域，可能参与叶片发育及杜仲胶的形成。图8表1参56     

毛竹GRF基因家族全基因组鉴定与表达分析

  目的  探究毛竹Phyllostachys edulis GRF基因家族的性质、结构特点以及在不同组织中的表达水平,为进一步研究GRF在毛竹生长发育中的分子作用机制奠定基础。  方法  采用生物信息学方法,对毛竹全基因组和转录组信息进行分析,筛选出13个毛竹GRF基因家族成员,并对GRF基因家族的理化性质、进化、基因结构、保守结构域、启动子、基因家族表达模式及三级结构进行分析。  结果  毛竹GRF基因家族成员按照其在scaffold上分布的位置,分别被命名为PeGRF01~PeGRF13;将含有3个内含子的PeGRF04和PeGRF12归为非ε组,其余为ε组。毛竹各GRF基因家族成员理化性质存在一定的差异,但是其结构域相对保守,均含有14/3/3结构域。毛竹GRF家族启动子区含有大量与光响应、低温响应、激素调控等相关的顺式作用元件。毛竹GRF存在基因复制扩增的现象,与水稻Oryza sativa的共线性关系明显高于拟南芥Arabidopsis thaliana。毛竹GRF在不同组织器官中均有表达,各家族成员的表达量存在差异;在毛竹花和根组织中的表达量略高于叶和鞭,同时家族成员间的表达量也存在一定差异。GRF蛋白质由2个单体构成,每个单体由9个α-螺旋构成,整体结构呈“W”型。  结论  毛竹GRF基因家族具有典型的14/3/3结构域,可能参与根、鞭、叶、花序及笋芽的生长发育过程。图6表1参39     

番茄LBD基因家族的全基因组序列鉴定及其进化和表达分析

【目的】从番茄全基因组中鉴定LBD基因,并进行基因进化、基因结构、染色体定位以及组织表达和诱导表达分析,为番茄LBD基因的功能研究与利用奠定基础。【方法】利用番茄基因组数据库,通过生物信息学手段,鉴定番茄LBD家族成员;采用MEGA5软件进行系统进化树分析;通过perl程序、MapDRAW及GSDS工具进行基因结构及染色体定位分析;利用已有的番茄芯片数据进行组织表达谱和基因表达响应分析。【结果】系统分析鉴定了 46 个番茄LBD 家族基因,根据基因结构及系统进化分析将其分成class I与class II两类,细分为5个亚家族（Ia、Ib、Ic、Id与II）。基因定位表明,12 条染色体中的 10 条均有LBD 基因,该基因家族的分布具有广泛性。各个发育阶段表达谱的分析结果表明,LBD 家族基因具有不同的组织表达模式,在各个发育时期均有LBD 基因的表达。基因响应分析发现,不同LBD基因响应不同的外界信号。【结论】通过全基因组分析,番茄LBD家族基因包括 46 个成员,在进化上分为两大类,5 个亚家族,分布于 10 条染色体上,组织表达模式及基因响应具有多样性。这些信息为番茄LBD基因家族的功能分析奠定了基础。     

毛竹磷转运蛋白Ⅰ家族基因鉴定及表达模式

  目的  鉴定毛竹Phyllostachys edulis磷转运蛋白Ⅰ (phosphate transporter 1, PHTⅠ)家族基因,分析其表达模式。  方法  利用生物信息学方法,鉴定毛竹PHTⅠ家族成员,分析基因启动子调控元件、编码蛋白的理化性质、基因结构、氨基酸保守基序、基因在染色体的位置、组织表达特异性、基因适应性进化及系统进化等。  结果  毛竹中共鉴定出20个PHTⅠ家族基因 (PePHTs),分布在10条染色体上,均定位于细胞膜。每个基因都含有1~2个内含子,PePHTs启动子序列中包含干旱、低温等非生物胁迫以及赤霉素等激素类响应元件。毛竹PHTⅠ大部分为碱性蛋白质,分子量为48.61~76.37 kDa,理论等电点为6.84~9.30,疏水性值均大于0,都属于疏水蛋白。基因适应性进化分析显示:大多数PePHTs的选择压力值小于0,说明多数基因受到负选择压力。转录组表达图谱表明:PePHTs在不同组织中的表达存在差异性,说明该基因家族在毛竹生长发育过程中发挥着不同的作用。系统进化树表明:PePHTs都聚类在第Ⅰ亚家族,并且优先和水稻Oryza sativa聚类在同一支上。  结论  PHTⅠ家族在植物吸收和转运磷的过程中扮演了重要作用。本研究结果为深入研究毛竹PHTⅠ家族基因的功能奠定了理论基础。图5表1参45     

黄瓜AQP基因家族的鉴定与生物信息学分析

  目的  深入研究黄瓜Cucumis sativus水通道蛋白(aquaporin, AQP)基因家族(CsAQP)的相关功能。  方法  通过全基因组分析技术鉴定其家族成员,对其蛋白质理化性质、系统进化关系、选择压力、基因结构、保守基序、顺式作用元件、蛋白质互作进行分析。  结果  黄瓜基因组共有33个AQP基因,含有2~5个数量不等的外显子,在染色体上不均匀分布;根据物种进化关系将黄瓜AQP基因家族划分为5个亚家族;基因组重复序列分析表明:5号和6号染色体上各有2~3对基因串联重复;计算这些基因的同义替换(synonymous, K_s)和非同义替换(nonSynonymous, K_a)的比率,结果显示均小于1,表明其进化受纯化选择作用;顺式作用调控元件分析发现,大部分基因启动子区所含元件与激素调节、光响应、胁迫密切相关。  结论  通过黄瓜全基因组扫描,获得黄瓜基因组的33个AQP家族成员,分属于5个亚族,映射于7条染色体上。上游启动子区含逆境相关作用元件,且部分基因参与串联复制,历经纯化选择。图6表4参26     

外施赤霉素对薄壳山核桃幼苗生长及相关代谢基因表达的影响

  目的  研究赤霉素(gibberellin)对植物株高的影响,探讨了GA20ox、GA3ox和GA2ox在赤霉素合成过程中的反馈调节机制。  方法  以薄壳山核桃Carya illinoensis幼苗为材料,进行100 mg·L?1 赤霉素叶面喷施处理。对外施赤霉素后薄壳山核桃的株高、节间长度和主根长以及薄壳山核桃GA20ox、GA3ox和GA2ox的时空变化进行了研究。  结果  外施赤霉素28 d后,薄壳山核桃节间长度和主根长的伸长量与对照相比都存在显著差异(P＜0.05),节间数量没有发生变化;植株平均生长量达2.9 cm,为对照的1.93倍,差异极显著(P＜0.01)。实时荧光定量PCR结果表明:GA20ox、GA3ox和GA2ox在薄壳山核桃生长期存在时空表达差异。外施赤霉素能够使CiGA20ox的表达量持续下降,28 d后的表达量仅为初始值的38.6%;而CiGA3ox表达量则是在7 d时跌落低谷,下降至初始值的55.4%,28 d后回升,表达量为初始值的350%;CiGA2ox表达量总体呈波浪形变化,在7 d后达到顶峰,之后有所回落,至21 d时恢复到了初始水平,28 d后上升为初始值的220%。CiGA3ox在薄壳山核桃植株茎杆内发生了累积,相较于CiGA20ox和CiGA2ox在转录水平上发生了较大的变化。  结论  外施赤霉素促进了薄壳山核桃茎杆伸长,同时引起了薄壳山核桃体内赤霉素代谢关键基因表达模式的时空变化。图6表1参29     

油松JAZ基因家族特征及其与DELLA蛋白互作的功能域鉴定

  目的  对油松JAZ基因家族特征及其与赤霉素负调控因子DELLA蛋白互作的功能域深入分析，旨在为解析针叶树以JAZ-DELLA为核心模块、茉莉酸（JA）-赤霉素（GA）介导的生长/防御平衡策略奠定基础。  方法  以油松全基因组数据为基础，筛选鉴定油松全部的JAZ家族基因成员，并分析其基本特征；构建多物种JAZ基因家族系统发育进化树，解析油松JAZ基因家族在系统进化过程中的特点；利用酵母双杂交技术，明确油松中JAZ与DELLA蛋白互作的功能域。  结果  油松中共鉴定得到53个JAZ基因家族成员，均具有TIFY和Jas保守结构，而且在degron序列中进化出了更丰富的变异。多个油松JAZ基因家族成员启动子区域包含响应JA和GA的顺式作用元件，并与模式植物中多个JAZ蛋白有着较近的进化距离。进一步实验发现，油松中5个JAZ蛋白（TIFY4、TIFY11、TIFY16、TIFY25、TIFY59）的Jas结构域与油松DPL（DELLA-like）蛋白存在相互作用，明确了油松中Jas结构域是JAZ蛋白与DELLA蛋白互作的功能域。  结论  明确了油松JAZ基因家族的基本序列特征，确定了油松中JAZ与DELLA蛋白互作的Jas功能域，为针叶树JAZ基因家族及JA-GA信号转导通路的深入研究提供了重要依据。      

基于MaxEnt模型预测薄壳山核桃在中国的种植区

  目的  预测当前和未来中国适宜种植薄壳山核桃Carya illinoinensis的地区,对薄壳山核桃的引种栽培提供参考。  方法  基于北美薄壳山核桃栽培区的分布数据和世界环境气候数据,利用MaxEnt模型和ArcGIS预测薄壳山核桃当前以及未来不同气候情景下在中国的适生区范围。  结果  ①训练集和检验集的AUC值分别为0.987和0.985,表明预测结果精度高。②刀切法检验结果和贡献率结果表明:年均气温、最暖季平均气温、年降水量、最干月降水量、最干季平均气温、昼夜温差月均值和气温的季节性是薄壳山核桃生长的关键气候变量。③预测薄壳山核桃当前在中国的适生区主要集中在中东部省份,包括江西、湖南、湖北、安徽、浙江、江苏和河南;在未来碳排放增强的气候情景下,薄壳山核桃在中国中东部省份的潜在高适生区的面积大幅增加,适生范围有向北移动的趋势,少量延伸到了辽宁东部和吉林南部。  结论  利用MaxEnt模型预测的薄壳山核桃适生区主要集中在中国的中东部省份;随着未来气温的升高,薄壳山核桃的潜在高适生区在中国相应的中东部地区内扩张。图3表3参36     

板栗Hsp90基因家族的鉴定及在胚珠发育过程中的表达分析

  目的  热激蛋白（Hsp90）是一种在细菌、植物、动物中均广泛存在的高度保守的蛋白质分子，在植物应激反应和生长发育等过程中发挥着重要的作用。板栗作为我国重要的坚果树种，雌雄花比例不同、授粉受精困难及胚败育等现象致使板栗单产低，严重制约板栗的产量。近期有研究发现Hsp90基因家族参与植物的胚形成，本研究首次对板栗中Hsp90基因家族进行鉴定及表达分析，旨在揭示Hsp90基因家族的功能及其在授粉受精和胚珠发育过程中的作用，进而为提高板栗产量提供理论依据。  方法  利用生物信息学方法对板栗Hsp90基因家族成员进行了鉴定，分析家族成员的理化性质、基因结构、进化关系、启动子元件。基于转录组数据分析CmHsp90基因在板栗胚珠不同发育时期的表达情况。  结果  本研究在板栗中共鉴定到10个CmHsp90基因，该基因家族成员的编码区（CDS）长度范围是2 112 ~ 2 481 bp，含有4 ~ 22个外显子，编码703 ~ 826个氨基酸，相对分子量范围为80.67 ~ 94.21 kDa，等电点pI范围介于4.81 ~ 5.26；不稳定系数介于36.15 ~ 46.75之间，除CmHsp90-6、CmHsp90-4、CmHsp90-9外，其他基因编码的蛋白质相对比较稳定。对CmHsp90基因家族蛋白进行亚细胞定位预测发现:5个成员（CmHsp90-1、CmHsp90-4、CmHsp90-5、CmHsp90-9、CmHsp90-10）定位到了细胞质上，2个成员（CmHsp90-7、CmHsp90-8）定位到了内质网上，其余3个成员分别定位到细胞核（CmHsp90-3）、线粒体（CmHsp90-2）和叶绿体上（CmHsp90-6）上。对CmHsp90基因家族蛋白信号肽预测发现，仅CmHsp90-7和CmHsp90-8存在潜在的信号肽位点。CmHsp90基因家族蛋白的motif、保守结构域基本一致，均包含Hsp90结构域（PF00183）和HATPase_C结构域（PF02518）。基于系统进化树分析，10个CmHsp90家族蛋白被分为5个组。对CmHsp90基因家族启动子区域进行顺式作用元件预测，发现光响应元件是CmHsp90基因家族启动子上的主要调控元件，其次为ABA响应元件，此外，有4个基因（CmHsp90-1、CmHsp90-5、CmHsp90-7、CmHsp90-8）存在温度响应元件。通过对板栗胚珠不同发育阶段的RNA-seq数据进行分析发现:除CmHsp90-3基因外，其他基因在受精前胚珠中的表达量均高于受精后胚珠；CmHsp90-2、CmHsp90-4、CmHsp90-5、CmHsp90-7、CmHsp90-8、CmHsp90-9、CmHsp90-10基因在发育胚珠中的表达量高于败育胚珠。  结论  CmHsp90基因在胚珠不同发育过程中的表达存在差异，在同一子房内发育胚珠和败育胚珠中的表达情况也不尽相同，研究结果为下一步研究CmHsp90基因在板栗授粉受精以及胚发育过程中的具体作用奠定理论基础。      

毛白杨PtYABBY基因家族3个成员的克隆与表达分析

  目的  YABBY是高等植物中特有的转录因子家族,在侧生器官发育和开花过程中发挥重要作用。本研究对毛白杨PtYABBY基因进行克隆和生物信息学分析,并对雌雄花芽8个关键发育时期和根茎叶等组织中PtYABBY基因表达模式进行探究,以期为阐明PtYABBY基因在毛白杨侧生器官发育过程中的功能奠定前期基础。  方法  以毛白杨为试材,采用同源基因克隆法从毛白杨中分离FILAMENTOUS FLOWER (FIL)/YABBY3 (YAB3)同源基因PtYABBY3、PtYABBY4和PtYABBY11的编码区序列,并开展生物信息学分析。通过qRT-PCR研究这3个基因在雌雄花芽8个发育时期及根、茎、幼叶、成熟叶片中的表达规律。  结果  PtYABBY3、PtYABBY4和PtYABBY11编码区长度分别为633、642、639 bp,分别编码210、213、212个氨基酸。所推测的氨基酸序列包含C₂C₂锌指结构域和YABBY结构域。系统进化分析进一步表明这3个基因属于FIL/YAB3亚家族成员。qRT-PCR结果显示这3个基因在根、茎、叶及8个时期雌雄花芽组织中均有表达,但PtYABBY基因间的表达水平存在明显差异。PtYABBY3表达水平在雌雄花芽发育初期下调,发育后期上调;花原基形成到休眠期内PtYABBY3在雌雄花芽中呈现相反的表达变化趋势。PtYABBY4和PtYABBY11基因在雌雄株成花诱导期表达水平最高,成花诱导到伸长期内呈现下降趋势,在雄花芽休眠期和小孢子发生期基因表达水平无明显变化。营养组织表达量测定结果显示这3个基因在幼叶和成熟叶片中表达水平相对较高,根系中表达量最低。  结论  YABBY基因家族成员PtYABBY3、PtYABBY4和PtYABBY11同属FIL/YAB3亚类,在毛白杨雌雄花芽发育的8个时期中表达水平存在较大差异,且在叶片和花芽中表达量较高,表明其可能与叶片和花芽发育相关。本研究为今后深入研究YABBY家族成员在杨树器官生长发育过程中的作用奠定基础。      

光皮桦AP2/ERF基因家族鉴定与表达分析

  目的  深入研究AP2/ERF基因家族在光皮桦Betula luminifera生长发育及环境胁迫响应中的生物学功能。  方法  利用光皮桦基因组数据,通过生物信息学方法开展AP2/ERF基因家族鉴定、基因特征、系统进化、基因结构、保守基序、顺式作用元件、蛋白互作和表达分析。  结果  在光皮桦基因组中共鉴定到77个AP2/ERF基因,其编码的蛋白理化性质存在差异,大多数蛋白(60个)的理论等电点小于7.0。系统进化分析显示:这77个AP2/ERF转录因子属于5个亚家族,其中ERF亚家族最大,包含34个成员。光皮桦AP2/ERF各亚家族间的基因结构存在较大差异,其中AP2亚家族成员均具有6~9个内含子,而DREB亚家族基因则没有内含子;但AP2/ERF各亚家族内的不同成员具有相似的保守基序类型和分布。同时,AP2/ERF基因启动子上都存在大量与激素、调节、胁迫响应及生长发育相关的顺式作用元件。此外,互作网络分析预测不同的光皮桦AP2/ERF亚家族蛋白间存在广泛的互作关系。进一步的表达分析显示:绝大多数光皮桦AP2/ERF基因(71个)的表达存在较强组织特异性,且在高温胁迫下,多数ERF或DREB基因的表达发生显著变化,表明ERF和DREB基因在高温胁迫应答中可能发挥着重要作用。  结论  通过生物信息学分析获得光皮桦77个AP2/ERF基因,分属于5个亚家族。不同亚家族基因具有相似的基因结构、保守基序等特征。基因启动子区含激素、胁迫响应等相关的作用元件。基因表达具有较强的组织特异性,且多数ERF或DREB基因对高温胁迫有明显响应。图6表1参39     

2023 年 1 期目次

  目的  Wuschel (WUS)相关的同源异型盒(Wuschel-related homeobox,WOX)转录因子家族在植物生长发育中发挥重要作用。本研究旨在探索WOX转录因子在杜仲Eucommia ulmoides中的分布及表达特征。  方法  以杜仲基因组数据库为基础,利用生物信息学方法对杜仲WOX家族进行全基因组鉴定;基于转录组数据分析EuWOXs在叶片发育及杜仲胶形成中的表达特征,通过实时荧光定量PCR (RT-qPCR)检测EuWOXs在‘紫叶’杜仲‘Ziye’叶片不同发育时期的表达模式。  结果  杜仲基因组中共鉴定出8条EuWOXs,分布于8条染色体;EuWOXs蛋白质长度为182~352个氨基酸,理论等电点为5.10~6.47,分子量为20.7~40.4 kDa;亚细胞定位预测EuWOXs均定位在细胞核中,均为亲水性蛋白。根据系统进化关系,杜仲WOX家族包括3个亚家族,分别含2、1和5个EuWOXs基因。EuWOXs均含有内含子,并包含多个基序,启动子中富含激素、胁迫和光周期响应元件。大部分EuWOXs在杜仲叶片中表达量较低,EuWOX13-1随叶片发育表达量逐渐降低,EuWOX13-2在生长叶中表达量最高。  结论  杜仲中有8个EuWOXs基因,EuWOX13-1和EuWOX13-2可能在杜仲叶片发育中发挥重要作用。图10表2参50     

毛果杨ZHD家族全基因组水平鉴定及在干旱胁迫下的表达分析

  目的  对毛果杨Populus trichocarpa ZHD (PtrZHD)家族进行生物信息学以及干旱胁迫下表达特性分析,为研究PtrZHD在干旱胁迫中的功能提供参考。  方法  利用生物信息学方法从全基因组水平鉴定出毛果杨ZHD家族全部成员,并对其进化、理化性质、基因结构、保守基序、启动子顺式作用元件和表达特性进行分析。  结果  毛果杨ZHD家族包括21个成员,可分为7个亚家族;有8对同源基因,且非同义替换率(K_a)/同义替换率(K_s)值远小于1。该家族成员理化性质存在差异,但其结构较为保守,均含有Motif 1;启动子区含有数量不等的激素和非生物胁迫响应元件,不同基因之间响应元件的种类存在差异。在毛果杨PtrZHDs中,分别有1、7和13个基因在根、茎和叶组织中具有偏好性表达特征;PtrZHD家族成员对干旱胁迫的响应具有组织和时间表达特异性,在根、茎和叶部组织中各成员的表达量不同,但随着干旱胁迫时间的增加均呈先上升后下降的趋势。  结论  PtrZHD家族基因对干旱胁迫有不同程度的响应,可调控毛果杨对干旱胁迫的应答。图6表2参27     

陆地棉磷高效基因GhMGD3的克隆与表达分析

  目的  基于前期陆地棉Gossypium hirsutum根部低磷胁迫基因表达谱芯片差异表达序列数据分析，挖掘相关基因，并对其克隆与表达分析。  方法  克隆GhMGD3基因并进行基因组DNA与cDNA测序分析，借助生物信息学方法分析GhMGD3的基因结构和进化关系；采用半定量RT-PCR技术与实时荧光定量PCR (RT-qPCR)的方法检测该基因在根、茎、叶、花4个组织中的基因表达量的变化以及低磷胁迫下的表达模式。  结果  成功克隆了陆地棉GhMGD3基因。GhMGD3基因的编码序列全长为681 bp，共编码226个氨基酸，存在3个内含子，分子量是26 610.54 Da，等电点为8.74，是稳定的亲水碱性蛋白。二级结构以α-螺旋和无规则卷曲为主。该蛋白不存在信号肽、跨膜结构域、N-糖基化位点，但含有多个磷酸化位点。亚细胞定位结果显示:该基因编码的蛋白定位于叶绿体。GhMGD3蛋白与木槿Hibiscus syriacus氨基酸序列相似性较高，其亲缘关系最近。半定量RT-PCR和RT-qPCR试验结果均表示:GhMGD3基因主要表达于根，中量表达于茎，微量表达于叶和花，在低磷胁迫72 h时其相对表达量达到最高值。  结论  首次成功克隆到了陆地棉GhMGD3基因，获得了GhMGD3基因的组织表达以及低磷胁迫下的表达模式，GhMGD3基因在棉花磷高效利用信号调控中具有重要作用。图7表1参27     

荆芥HD-Zip基因家族的全基因组鉴定及分析

【目的】鉴定荆芥Schizonepeta tenuifolia的HD-Zip基因家族,利用生物信息学方法分析其在全基因组中的分布和相关特征以及在不同时期中的表达规律,为该家族基因的进一步研究奠定基础。【方法】根据已经表征的HD-Zip基因,筛选荆芥基因组内的HD-Zip基因序列,利用MEME、PlantCARE、NCBI、MEGA X、MCScanX、Circos等在线网站及软件对蛋白序列进行基本理化性质分析、进化树构建、染色体定位、基因结构分析、共线性基因分析等。【结果】在荆芥全基因组中共鉴定到42条HD-Zip基因序列,它们可被分为4个亚家族,分别含有16、7、5、14个基因,亚家族之间的基因长度、结构及保守基序差异显著,但在亚家族内部保守,荆芥基因组与拟南芥Arabidopsis thaliana基因组共线性分析发现有37对基因,可能具有相似的生物学功能。荆芥的4个亚家族基因的顺式元件中均高频出现了光响应、脱落酸响应、MeJA响应等元件,在不同生长时期的叶片及部位的转录组数据中具有不同的表达趋势,Ⅰ亚家族主要在幼叶中表达,Ⅱ和Ⅲ亚家族主要在根中在表达,Ⅳ亚家族主要在叶中表达。【结...     

牡丹响应高温胁迫的转录组分析及PsHSP基因表达

  目的  夏季高温可导致牡丹Paeonia suffruticosa生长受限,花期缩短,观赏品质降低。研究牡丹耐热基因对牡丹热胁迫的分子机制。  方法  以牡丹品种‘羽红’‘Yuhong’为材料,对高温(40 ℃)和对照(25 ℃)处理后的叶片进行转录组测序,分析响应高温表达的差异基因;同时利用实时荧光定量PCR(qRT-PCR)验证差异基因,并对热激蛋白基因HSPs进行时空表达分析。  结果  测序共获得45.97 Gb数据,与对照组(25 ℃)相比,高温处理后4 220个基因上调表达,3 453个基因下调表达。通过基因本体(GO)分析发现,这些差异基因主要富集于代谢等生物学过程、细胞和膜结构等细胞成分以及结合和催化活性等分子功能条目;另外,通过全基因组及代谢途径数据库(KEGG)分析发现,碳代谢途径的差异基因数量最多。对牡丹热激蛋白PsHSP基因定量分析发现,PsHSP基因随着高温处理时间呈上升表达趋势,24 h达到最高,之后便呈现下降趋势。  结论  高温显著影响牡丹的代谢和生物合成等过程,进而影响牡丹的生长和发育。PsHSP基因可在短期内迅速响应高温胁迫并参与牡丹的耐热调控。图 4表6参34