茶树BES1转录因子全基因组鉴定与分析

植物特异性转录因子(BRI1-EMS-SUPPRESSOR1,BES1)家族参与调节油菜素内酯(brassinosteroids,BR)信号通路,在植物生长和应对非生物胁迫的反应中发挥重要作用。该研究利用生物信息学方法对茶树BES1转录因子家族进行鉴定,分析茶树CsBES1基因在8个茶树组织中的表达模式,并对CsBES1转录因子家族在低温、干旱及ABA激素处理下的表达进行分析,以揭示茶树CsBES1转录因子家族的功能及其对环境胁迫的响应。结果显示:(1)从茶树全基因组中鉴定获得10个茶树BES1转录因子家族成员,均具有完整的BES1_N结构域;根据系统进化关系将10个CsBES1转录因子家族分成3组,分别包含2个、3个和5个成员;该家族成员每个基因含有2～11个外显子。(2)组织表达模式分析显示,茶树BES1转录因子家族在生长活跃的茶树嫩梢和成熟叶中表达量较高,不同成员之间表达存在差异性。(3)上游启动子区域分析发现大量与植物激素和非生物胁迫响应相关的顺式作用元件。(4)荧光定量检测发现,BES1基因在茶树不同胁迫处理下的表达模式不同,在低温和脱落酸处理下,CsBES1-1、CsBES1-7、CsBES1-8和CsBES1-9基因的表达量显著提高,而在干旱处理下,CsBES1-2、CsBES1-4、CsBES1-5、CsBES1-6、CsBES1-7、CsBES1-8和CsBES1-9基因的表达量显著提高。     

唐菖蒲赤霉素受体基因克隆及表达分析

以唐菖蒲子球为材料,采用RT PCR技术,克隆了1个赤霉素（GA）受体基因,命名为GhGID1a（GenBank登录号为KU525107）。其开放阅读框为1 032 bp,编码343个氨基酸。序列比对结果表明,GhGID1a推导氨基酸序列与百子莲、油棕和海枣等植物的GID1的氨基酸序列相似性较高,分别为82%、78%和77%。50 mg/L GA₃对子球萌发具有促进作用,150 mg/L GA₃会抑制其萌发。实时荧光定量PCR结果显示,GA3处理对GhGID1a基因表达具有反馈抑制作用,GhGID1a表达量随着子球休眠解除而逐渐降低,推测唐菖蒲子球休眠解除可能与GA及其受体基因有关,GA及其受体基因GhGID1a可能参与了调控唐菖蒲的子球休眠与萌发过程。     

茶树CsLEA5基因的克隆及表达分析

该研究以茶树‘龙井长叶’为材料,克隆获得了茶树胚胎发育晚期丰富蛋白基因CsLEA5的cDNA序列,该序列全长515 bp,包含一个375 bp的开放阅读框,编码124个氨基酸,预测蛋白分子量为13.5 kD,理论等电点为5.92。蛋白序列分析结果显示,CsLEA5为高亲水性和稳定性蛋白,且含有一个典型的LEA_3保守结构域,属于LEA蛋白中LEA_3亚家族成员。CsLEA5基因启动子区域包含多种与逆境响应相关的顺式作用元件,如乙烯响应元件（ERE）、胁迫响应元件（STRE）、创伤响应元件（WUN motif）及MYB、MYC转录因子识别位点等。qRT PCR分析显示,CsLEA5基因表达具有明显的组织特异性,在叶片中的表达量最高,其次是嫩茎,而在其他组织器官中的表达量较低,且CsLEA5基因表达受低温和干旱胁迫的诱导。研究表明,CsLEA5基因可能在茶树响应低温和干旱胁迫过程中发挥重要作用。该研究对了解茶树抗逆分子机制,筛选抗性候选基因资源提供了重要理论依据。     

香菜CsEF 1α基因克隆与表达分析

翻译延伸因子EF 1α(elongation factor 1 alpha)是细胞中最丰富的蛋白质之一,其在确保mRNA正确解码以产生细胞蛋白质方面发挥重要作用。该研究采用RT PCR扩增方法克隆香菜CsEF 1α基因序列,利用生物信息学对CsEF 1α基因结构、序列特征及系统进化等进行分析,并采用qPCR探究CsEF 1α基因在香菜不同生长时期和非生物胁迫下的表达模式,为进一步揭示EF 1α基因调控机制的研究奠定基础。结果显示：(1)成功克隆获得香菜CsEF 1α基因序列;CsEF 1α基因包含1个1 344 bp的开放阅读框,编码447个氨基酸,分子式为C₂₂₀₂H₃₅₄₄N₅₉₄O₆₄₄S₂₀,蛋白质分子量为49.29 kD,等电点为9.12;氨基酸序列组成中赖氨酸数量最多(49个,占11.0%),色氨酸数量最少(3个,占0.7%);属碱性蛋白。(2)CsEF 1α蛋白主要由无规则卷曲(36.91%)和α 螺旋(30.43%)构成,定位于细胞质;系统进化树分析显示,CsEF 1α与胡萝卜、野生番茄、青蒿素和非洲菊的亲缘关系较接近;启动子分析包括4种植物生长发育元件、3种激素响应元件和3种胁迫响应元件。(3)qRT PCR结果显示,CsEF 1α基因的表达量随着香菜生长发育时间的延长而升高,并且与转录丰度的变化一致;CsEF 1α基因对4种不同非生物胁迫的响应表达模式有所差异;随着胁迫时间的延长,在盐胁迫下CsEF 1α基因表现出先升高后降低趋势,而在低温、高温和干旱胁迫下表现出先降低再升高的趋势。研究表明,CsEF 1α基因参与了香菜对非生物胁迫的应答,在香菜生长发育和非生物胁迫中具有重要调控作用。     

牡丹PsAGL6基因克隆与表达分析

该研究在转录组数据分析基础上,通过RT-PCR方法从牡丹(Paeonia suffruticosa.L)品种‘洛阳红’中克隆AGL6基因,并分析了其在不同组织和花型中的表达模式。结果显示:(1)AGL6开放阅读框长度732bp,编码244个氨基酸;结构域分析显示,该基因具有高度保守的MADS MEF2区、K区、AGL6I与AGL6II基序,归类于MADS-box基因家族,命名为PsAGL6,GenBank登录号为MF563611。(2)同源比对分析显示,牡丹PsAGL6与葡萄VvAGL6氨基酸序列相似度最高,达79%。(3)qRT-PCR分析结果显示,PsAGL6在牡丹各器官中均有表达,但营养器官中表达量较低,花器官中表达量较高,其中以萼片最高,花瓣与雌蕊次之;对4种牡丹花型花瓣的表达分析显示,PsAGL6基因在不同花型花瓣中的表达量差异明显,且蔷薇型牡丹花瓣相对表达量最高。研究表明,AGL6基因参与牡丹花器官的形成,为深入研究花器官发育的分子机制提供了帮助。     

白菜型油菜RbohC和RbohF基因克隆与表达分析

该研究以白菜型油菜（Brassica rapa L.）‘陇油6号’为实验材料,采用RT PCR方法克隆油菜RbohC和RbohF基因,并采用实时荧光定量PCR技术对RbohC、RbohF基因在不同组织及非生物胁迫下的表达进行分析,为深入研究油菜RbohC、RbohF基因的生物学功能提供依据。结果显示：(1) 成功克隆得到2个全长分别为3 050 bp和2 995 bp的油菜RbohC (GenBank登录号：XM_009134386) 和RbohF (GenBank登录号：XM_009114548) 基因序列。(2) 生物信息学分析显示,油菜RbohC和RbohF基因开放阅读框（ORF）分别为2 733 bp和2 847 bp,编码910和948个氨基酸,推测二者的蛋白质分子量分别为103 kDa和108 kDa,理论等电点分别为9.47和9.21; 油菜RbohC、RbohF编码的氨基酸序列与萝卜等多种植物相应蛋白氨基酸序列具有较高的同源性,且这些序列高度保守并含有NADPH氧化酶的典型保守结构域,包括2个可以与Ca²⁺结合的EF手性模体结构、6个跨膜结构域、黄素腺嘌呤二核苷酸结合结构域、NAD焦磷酸结合结构域和C末端区域中的NADP核糖保守结合位点。(3) 油菜RbohC和RbohF基因在根、茎、叶和下胚轴中均表达,无组织特异性,但RbohC基因在根中表达量最高, RbohF基因在下胚轴中表达量最高。(4) 低温、干旱、盐、ABA、H₂O₂处理都能够诱导油菜RbohC和RbohF基因的表达,但抗寒性强的 ‘陇油6号’的RbohC和RbohF基因对胁迫的响应更敏感,且RbohC基因的表达量均高于RbohF基因。(5) 用H₂O₂清除剂DMTU、NADPH氧化酶抑制剂DPI和IMD、MAPKK抑制剂U0126处理后,油菜RbohC和RbohF基因的表达均较对照下降,说明U0126和DMTU对油菜RbohC和RbohF基因的表达有抑制作用。研究认为,油菜RbohC和RbohF基因在油菜适应逆境胁迫中具有重要作用,两基因的表达均受MAPK激酶信号途径的调节,并受到H₂O₂的反馈调节,而且抗寒性强的‘陇油6号’品种中RbohC和RbohF基因对H2O2和MAPK激酶信号途径的响应更敏感。     

菜心SOD基因的克隆及其在低温条件下的表达分析

采用同源克隆法从菜心中获得3个SOD基因,并进行生物信息学分析,采用qRT-PCR分析3个基因在不同组织器官和低温胁迫下的表达模式。结果表明:(1)获得Cu/Zn-SOD、Fe-SOD、Mn-SOD基因的ORF,分别命名为BclCZSD、BclFSD、BclMSD,序列长分别为459、639、696bp,分别编码152、212、231个氨基酸。(2)生物信息学分析显示,3种蛋白均为稳定的亲水性蛋白,均不存在跨膜结构和信号肽,BclCZSD的二级结构以无规则卷曲为主,含有2个Cu/Zn-SOD结构域,BclFSD和BclMSD的二级结构以α-螺旋为主,含有一个相同的Mn/Fe-SOD结构域;进化分析显示BclCZSD和BclMSD与油菜最先聚在一个分支,BclFSD与甘蓝、萝卜、油菜、芜菁聚在一个分支。(3)qRT-PCR结果显示,BclCZSD、BclFSD和BclMSD基因在菜心根、茎、叶和叶柄中均有表达,且在根、茎、叶和叶柄中的表达模式不完全相同;低温条件下,3个基因的表达量随着胁迫时间的延长均呈先升高后降低的变化趋势。研究结果为进一步探讨菜心SOD基因在低温胁迫下的响应机制奠定了基础。     

藜麦CqSAP8基因克隆及其在非生物胁迫下的表达分析

该研究根据NCBI公布的藜麦胁迫相关蛋白(SAP)基因CqSAP8序列进行克隆,利用生物信息学分析CqSAP8蛋白序列、理化性质和结构特点,并采用qRT PCR方法检测CqSAP8基因表达的组织特异性及在非生物胁迫下的相对表达。结果表明：(1)藜麦CqSAP8基因CDS全长528 bp,编码175个氨基酸;预测CqSAP8蛋白分子量为18.73 kD,理论等电点为7.46,属于稳定的亲水性蛋白质;CqSAP8蛋白在N端和C端分别含有A20和AN1保守域,是SAP蛋白最典型的类型。(2)序列比对与进化分析显示,藜麦CqSAP8与甜菜BvSAP8、菠菜SoSAP8亲缘关系最近,序列相似性分别为89.66%和89.47%。(3)qRT PCR分析表明,藜麦CqSAP8基因在根、茎、叶、花和种子中均有表达,且在种子中表达量最高;CqSAP8基因在干旱和高温胁迫12 h表达量达到最大值,分别是对照的13.09和17.47倍,高盐和低温胁迫下的最大表达量均为对照的3.91倍,说明藜麦CqSAP8基因响应多种非生物胁迫应答;另外,藜麦CqSAP8的表达量在ABA胁迫下24 h急剧升高,推测CqSAP8基因在非生物胁迫前期的响应不依赖ABA。该研究为进一步研究CqSAP8基因功能及抗逆分子机制奠定了基础。     

茶树CsCML24的克隆及表达分析

类钙调素(calmodulin-like protein,CML)是植物体内一类钙受体蛋白,介导Ca²⁺与下游靶蛋白的相互作用,在植物抗逆反应中发挥重要作用。探究茶树中CML蛋白在逆境胁迫中的功能,为进一步研究茶树CsCML24对逆境胁迫的响应机理提供理论依据。以龙井43一年生茶树扦插苗为材料,克隆得到类钙调蛋白基因CsCML24(GenBank登录号为MZ325391),并进行了生物信息学分析。同时,利用实时荧光定量PCR技术分析茶树CsCML24组织表达特异性以及不同非生物胁迫处理下的表达模式。结果显示,CsCML24的CDS长为480 bp,编码159个氨基酸。CsCML24为无信号肽及跨膜结构的稳定性亲水蛋白,含有EF-hand保守结构域,能与Ca²⁺结合。qRT-PCR表明,CsCML24在茶树各组织中均表达,在成熟叶的表达量显著高于其他组织,茎中表达量较低;该基因在低温(10℃)、干旱(20% PEG 6000)、高盐胁迫(200 mmol/L NaCl)和ABA(100 μmol/L)处理均能诱导表达,且在不同胁迫下表达差异显著...     

茶树CsGME1基因的克隆与逆境胁迫响应分析

该研究基于茶树转录组和基因组信息,以茶树‘龙井43’为实验材料,从其cDNA中克隆获得茶树CsGME1基因,并对其蛋白序列特征、基因表达模式及其在不同非生物胁迫下的表达水平进行实时荧光定量分析。结果显示:(1)茶树CsGME1开放阅读框长度为1131 bp,编码376个氨基酸;该序列与多个相关物种的GME氨基酸序列一致性为94.25%,均含有NAD结合域。(2)进化树分析表明,茶树CsGME1基因与番茄SlGME1亲缘关系较近,与水稻OsGME2亲缘关系最远。(3)CsGME1蛋白属于亲水性蛋白,理论相对分子量为42046.84 Da,理论等电点为5.73,具有4个无序化区域,无序化程度较低;CsGME1蛋白二级结构由39.25%α-螺旋,13.26%延伸主链,5.84%β-转角和41.38%随机卷曲组成;三级结构分析结果显示,CsGME1包含螺旋和随机卷曲,与二级结构吻合。(4)荧光定量分析结果显示,在高温(38℃)、低温(4℃)、干旱(20%PEG)和高盐(200 mmol·L-1 NaCl)4种非生物胁迫处理下,茶树CsGME1基因均有响应,且表达存在差异,推测CsGME1参与了茶树的逆境胁迫响应过程。     

茶树CsGIF1基因的克隆及其对非生物胁迫的响应分析

GRF-INTERACTING FACTOR(GIF)基因是植物叶发育相关的一类重要调控因子,调节植物叶器官的发育。该研究采用RT-PCR方法从茶树‘龙井43’的叶片cDNA中克隆得到CsGIF1基因,并利用荧光定量PCR分析了高温(38℃)、低温(4℃)、干旱(200g·L-1 PEG)、盐胁迫(200mmol·L-1 NaCl)下CsGIF1基因的表达水平,以明确CsGIF1基因对非生物胁迫的应答特性,为茶树CsGIF1基因的逆境调控以及功能研究奠定基础。结果表明:(1)CsGIF1基因长666bp,编码221个氨基酸,具有高度保守的SNH结构域;CsGIF1蛋白为亲水性蛋白,理论相对分子质量为23 380,理论等电点为6.30;酸性氨基酸、碱性氨基酸、芳香族氨基酸和脂肪族氨基酸所占比例分别为7%、9%、5%和13%;蛋白质二级结构显示,茶树CsGIF1蛋白由38.01%的α-螺旋、10.41%的β-折叠、12.22%的延伸主链和39.37%的随机卷曲组成。(2)qRT-PCR分析显示,茶树CsGIF1基因对于4种非生物胁迫均有响应,但不同处理响应不同。其中在胁迫2h时,高温和低温胁迫下,CsGIF1基因相对表达量显著增加,分别为对照的3.54和5.69倍,且高低温胁迫下,CsGIF1基因响应明显大于干旱和盐胁迫。在高温、低温、干旱、盐等不同胁迫处理下,茶树‘龙井43’中CsGIF1基因的表达水平与对照相比均差异明显。     

茶树谷丙转氨酶基因的克隆及其表达分析

该研究基于茶树的转录组数据,采用RT PCR方法从茶树‘黄金芽’cDNA中克隆获得茶树谷丙转氨酶基因（CsAlaAT）,利用荧光定量PCR方法,对CsAlaAT在茶树材料‘迎霜’和‘黄金芽’不同组织、温度胁迫与激素处理的表达进行分析。结果显示：CsAlaAT基因开放阅读框长度为1 662 bp,编码553个氨基酸,含有天冬氨酸转氨酶家族（Aspartate aminotransferase family）典型的AAT like保守结构域。多序列比对显示,该序列与多个相关物种的序列一致性达78.83%,与磷酸吡哆醛（PLP）结合的10个氨基酸残基以及第358位赖氨酸催化位点在物种间高度保守。茶树CsAlaAT蛋白属亲水性蛋白,相对分子质量为60 877.5 D,等电点为6.11,碱性、酸性、脂肪族和芳香族氨基酸比例分别为12%、11%、22%和8%,无序化特征不明显,与大麦HvAlaAT具有相似的三维结构。实时定量PCR分析表明,CsAlaAT在茶树‘迎霜’和‘黄金芽’中的表达均具有组织特异性,且均在根部的表达量最高;CsAlaAT响应高温（38 ℃）和低温（4 ℃）胁迫的表达上调;外源施用脱落酸（ABA）和赤霉素（GA）能够抑制茶树中CsAlaAT基因的表达。     

茶树CsBAP1基因的克隆与非生物胁迫和激素响应分析

该研究基于茶树转录组数据,从茶树‘龙井43’cDNA中克隆获得茶树CsBAP1基因,对其蛋白质序列特征和基因表达模式以及CsBAP1在茶树不同组织、不同激素处理及不同逆境胁迫下的表达水平进行荧光定量检测。结果显示：（1）克隆得到的茶树CsBAP1基因开放阅读框长度为543 bp,编码180个氨基酸;CsBAP1蛋白为亲水性蛋白,理论相对分子质量为19 398 Da,理论等电点为9.30,含有保守的Ca²⁺依赖性的C2结构域;茶树CsBAP1蛋白二级结构由8.89％的α 螺旋、7.78％的β 转角、35.56％的β 折叠和47.78％的随机卷曲组成;三级结构分析结果显示,CsBAP1包括螺旋和折叠结构,与二级结构吻合。（2）荧光定量分析结果显示, CsBAP1基因在茶树的花、花蕾、幼叶、老叶和成熟叶中均有表达,且在老叶中的表达量最高;外源施用ABA、IAA、MeJA、GA³以及SA均能够抑制茶树CsBAP1基因的表达;在高温（38 ℃）、低温（4 ℃）、干旱（200 g·L^-1 PEG）、高盐（200 mmol·L^-1 NaCl）胁迫下,CsBAP1的表达量均有所上调,且不同时间段的表达存在显著差异。该研究为进一步研究茶树CsBAP1基因的功能奠定了基础。     

茶树硝酸还原酶基因克隆及表达分析

利用茶树全器官转录组文库中硝酸还原酶(NR)的EST,通过RACE技术扩增出NR基因的cDNA,并利用实时荧光定量PCR检测了NR基因在不同茶树品种中的表达。结果表明:NR基因cDNA全长2 927bp,开放阅读框2 652bp,编码一个有884个氨基酸蛋白质,GenBank登录号为JX987133。经BlastX比对,与GenBank中登录的烟草NR相似性达到74%。茶树NR蛋白属于亲水性蛋白,可能为胞质蛋白。25个茶树品种叶片中NR表达水平差异明显,最高值是最低值的22.75倍。因NR是植物氮代谢过程中的关键限速酶,推测25个茶树品种间氮吸收利用能力存在差异。     

茶树CsMYB123转录因子的克隆及表达特性研究

该实验以茶树品种‘紫娟’为试验材料,利用RT-PCR方法,从茶树cDNA中克隆得到一个R2R3-MYB型基因(CsMYB123)。生物信息学分析显示,CsMYB123基因的开放阅读框为915bp,编码304个氨基酸,蛋白分子量约34.07kD,理论等电点为8.69,含有2个保守的MYB结构域,编码1个R2R3-MYB蛋白;CsMYB123蛋白与拟南芥(Arabidopsis thaliana)MYB转录因子家族第五亚组的AtMYB123亲缘关系最近;CsMYB123属于亲水性蛋白,无N端信号肽,可能定位于细胞核上。荧光定量PCR分析表明,CsMYB123基因在茶树各组织的表达量大小依次为:一芽一叶第二叶第三叶第四叶老茎嫩茎,且在一芽一叶中的表达量是嫩茎的15.68倍;但CsMYB123的表达受IAA、ABA、ETH和GA3的抑制。花青素含量检测显示,茶树‘紫娟’各组织中花青素的含量高低依次为:第二叶一芽一叶第三叶第四叶嫩茎老茎,且第二叶和一芽一叶的含量分别为老茎的15倍和11倍。研究发现,CsMYB123基因在茶树‘紫娟’的新稍中高水平表达,且其表达模式与不同组织中的花青素含量呈较好的正相关关系,推测CsMYB123基因与茶树花青素合成的调控相关。     

茶树GRF基因家族的全基因组鉴定及表达分析

生长调控因子(growth regulating factor,GRF)是植物特有的转录因子家族,在植物生长发育过程中起重要调控作用。该研究利用生物信息学的方法,从茶树基因组中鉴定获得11个CsGRF转录因子家族成员,且均具有完整的特征结构域QLQ和WRC。CsGRF家族成员含有3~6个外显子,基于系统进化关系将CsGRF家族分为6组,且与葡萄及猕猴桃GRF家族的亲缘关系更为接近。不同组织转录组数据分析结果表明,该家族在生长活跃的茶树嫩梢部位高表达。上游启动子区域分析发现大量与植物发育、激素及胁迫响应密切相关的顺式作用元件。荧光定量检测发现,分别有10个和2个CsGRF成员在低温和干旱胁迫处理后呈现显著上调表达,其中CsGRF8和CsGRF11对2种非生物胁迫均有响应;并发现受ABA、MeJA和GA激素处理诱导分别有9个、3个和6个CsGRF基因的表达水平差异显著。研究表明,CsGRF基因家族在茶树生长发育过程及应激反应中起作用,推测CsGRF基因可能通过各种激素信号转导途径参与胁迫响应过程。     

茶树中2个NAC转录因子基因的克隆及温度胁迫的响应

利用茶树转录组数据库,检索得到2个NAC家族转录因子基因CsNAC1和CsNAC2。通过RT-PCR方法,将其从茶树‘迎霜’中分离克隆,利用荧光定量PCR方法,对CsNAC1和CsNAC2基因在‘迎霜’和‘安吉白茶’2个茶树品种不同组织以及温度胁迫处理下的表达进行分析,以探讨NAC家族转录因子在温度胁迫下的响应特征。结果表明:(1)CsNAC1和CsNAC2基因开放阅读框长度分别为1 044和1 047bp,分别编码347个和348个氨基酸;蛋白功能域预测和多重对比显示,CsNAC1和CsNAC2蛋白N端均含有典型NAC家族成员所具有的NAM保守结构域。(2)进化分析表明,CsNAC1和CsNAC2分别属于NAC家族的NAP和AtNAC3亚家族。(3)三维分子模型建模显示,CsNAC1和CsNAC2蛋白分别含有3个和2个α-螺旋,6个和7个β-折叠。(4)荧光定量PCR结果显示,CsNAC1在2个茶树品种中具有较相似的组织特异性,均在茶树成熟叶中表达量最高;CsNAC2则分别在‘安吉白茶’的幼叶中,‘迎霜’的根中表达量最高;高温(38℃)和低温(4℃)处理下,CsNAC1和CsNAC2基因的表达均受不同温度胁迫影响,不同茶树品种、不同时间段的表达存在差异。     

茶树CsCBF1基因克隆和转录活性分析

采用RACE技术克隆了一个受冷诱导的茶树CBF基因全长cDNA,命名为CsCBF1(GenBank登录号为EU563238)。CsCBF1cDNA全长序列为1 211bp,开放阅读框编码259个氨基酸。氨基酸序列分析表明,CsCBF1具有CBF家族典型的保守结构域,与其他植物的CBF具有较高的相似性;与拟南芥、辣椒和橡胶树编码的CBF相似性分别为56%、63%和56%。亚细胞定位结果表明,CsCBF1位于细胞核内。分别将10个CsCBF1缺失突变体与GAL4DNA结合域融合的结果显示,CsCBF1的羧基末端酸性结构域(第137位氨基酸至259位氨基酸)在酵母中具有转录激活活性。实时定量RT-PCR分析表明,CsCBF1基因受低温的快速诱导表达。