植物都有适宜的生长温度，温度过高或过低都会影响其生长发育，降低作物产量。转录因子(transcriptionfactor,TF)可以与不同基因启动子区域的顺式作用元件结合或者与其他转录因子家族的成员及调控蛋白相互作用，从而达到对生物体内某些特定生理生化过程的调控

东农冬麦1号是首例在黑龙江省安全越冬的强抗寒冬小麦品种。植物抗寒的生理生化反应是由低温改变基因的表达所引起，目前已经克隆了许多抗寒相关基因。了解植物安全越冬的分子机制对提高其抵抗低温的能力具有重要的作用，在生产实际中也具有广泛的应用价值。bZIP是植物中比较大的一类转录因子家族，几乎参与调控植物的所有生命现象，如抗逆。bZIP(basic-re-gionleucinezipper)的结构域由碱性氨基酸域和亮氨酸拉链域组成，总长度60-80个氨基酸CsC碱性氨基酸区域位于bZIP结构域的N端，由18个富含碱性氨基酸(Rrg}Lys)、不变的N-x7-R/K-x9基序组成;该区域负责与DNA顺式作用元件结合，决定了DNA结合的特异性，并具有核定位信号序列;该区域相比于亮氨酸拉链区域更加保守。亮氨酸拉链区域一般由7个或9个氨基酸组成一个重复单位，每个重复单位的第7或第9位一般是亮氨酸，但也有其他一些疏水性氨基酸(Ile}Val}Phe或Met)。该区域二聚化使富含亮氨酸区域形成亮氨酸拉链结构，然后结合DNA区。bZIP结构域的N端，也就是碱性区域负责与双链DNA的大沟结合，而C'端(亮氨酸拉链域)通过二聚化介导形成a一超螺旋结构，即亮氨酸拉链域。在拟南芥和水稻中，与抗逆相关的bZIP转录因子的功能研究相对比较深入。目前，在模式植物拟南芥和水稻中已分别鉴定出127个和140个bZIP转录因子，六倍体小麦发现了102个bZIP转录因子。小麦中与非生物胁迫相关的bZIP转录因子的研究相对较少，尤其是与低温相关的转录因子研究有限，所以在冬小麦抗寒过程中，转录因子起到一个什么样的作用，仍需深入研究。

小麦中与抗逆相关的bZIP转录因子基因(TabZIP)的转录水平表达受ABA、低温、高盐和机械伤害诱导，且参与条锈病菌的防卫反应。相关文献中的低温胁迫一般为4`C左右，而在0`C以下低温胁迫时TabZIPl的作用如何，尚未见报道。本研究对这方面进行了探讨，并对其进行生物信息学预测分析，同时通过定量方法测定其在低温下的相对表达水平，进而探索转录因子调控冬小麦抗寒性的机制。

1 材料与方法

1.1 材料

试验在东北农业大学试验田进行，参试冬小麦为强抗寒性品种东农冬麦1号和弱抗寒性品种济麦22，于2015年9月7日播种，完全区组设计，3次重复，行长2m，行距0.2m,10行区。播种量450粒·m-z，播深5cm，常规田间管理。

待麦苗长至三叶期(2015年9月23日)时，对冬小麦叶片喷施10-5mol"I,-1ABA(本课组前期试验得出大田试验的最适浓度)，喷施量为4mz"I,-i，对照组喷施等量的蒸馏水。于大田自然降温至5`C}0`C、一10`C和一25`C(连续tod最低温的平均值)时，分别取样分孽节和叶片，一80`C保存备用。

1.2 实验方法

1.2.1生物信息学分析

编码蛋白氨基酸序列的基本理化性质由Ex-PASy-ProtParamtool软件分析C}-8];蛋白的保守结构域由ConseveredDomains工具预测;二级结构用S()PMA软件进行预测与分析;利用NCBI中的Blastp在线工具进行其他物种的同源氨基酸序列搜索;利用ClustalX(version2.0)和DNAMAN(version8.0)进行氨基酸序列同源比对]9]，用MEUA6.0软件构建系统进化树;利用TargetP1.1Server和ChloroP1.1Server预测蛋白的亚细胞定位。

1.2.2 RNA提取及cDNA合成

Trizol法提取叶片及分孽节总RNAocDNA第一链的合成参照反转录试剂盒(购自TaKaRa公司)说明书进行。

1.2.3 RT-qPCR

根据已得到的TabZIPl基因(UenBank登录号FJ194457)序列设计实时定量引物，TabZIPl转录因子基因的正向引物和反向引物分别为5‘-TUUUAAUCAUTUAAUCAUAU-3‘和5'-CCACAUUAUAUUUAAAUAACC-3‘。

实时荧光定量PCR反应系统为Mx-3000pRe-al-TimePCRSystem。以小麦中的l1cti、为内参基因，其扩增引物序列为5'-CCTTAGTACCTTCCAACAUATUT-3‘和5'-CCAGACAACTCGCAACTTAUA-3'。荧光实时定量PCR(RT-qPCR)反应体系为20pI,:SYBR.PremixExTuyTMII(2X)10pL,cDNA模板2pI,，正向及反向引物(10pM)各0.8pL,ROXReferenceDyeII(50X)0.4pL,ddH2O6pLoPCR反应条件为95C预变性30s,95`C变性5s,60`C退火34s,40个循环。采取2-ooa法计算待测基因的相对表达量mo，C,，值即为达到检测阑值时的PCR循环数。3次生物学重复，数据显著性分析采用SPSS软件。

2 结果与分析

2.1 TabZIP1转录因子的生物学信息

2.1.1 TabZIPl的理化性质

NCBI中检索基因号(FJ194457)显示，TabZIPl基因全长1167by，编码388个氨基酸。ExPASy-ProtParamtool预测其表达蛋白的相对分子质量为94841.9Da，理论等电点(pI)为5.03，是弱酸性蛋白;其不稳定系数(Instabilityindex)为50.66，属于不稳定型蛋白质(X40，蛋白稳定);脂肪族系数(Aliphaticindex)为31.88，总平均亲水系数(Grandaverageofhydropathicity)为0.924(在2一一2范围内，正值表明此蛋白为疏水蛋白，负值表明为亲水蛋白)，表明该蛋白是疏水蛋白。该蛋白由4种氨基酸组成，分子式为}-sasiHs}s}Nms}}}iaosSz}s，其氨基酸所占比例分别为31.900(Ala)}23.600(Cys)}24.400(Gly)和20.10o(Thr)。保守结构域预测显示，该蛋白具有保守的BRLZ基序，保守区域从第291一第355个氨基酸，保守区域全长65个氨基酸，属于bZ-IP1家族。Signa1P4.1Server表明，该蛋白无信号肤或跨膜结构域，属于非分泌性蛋白。

2.1.2 TabZIPl蛋白的二级结构

利用S()PMA在线对TabZIPl蛋白进行二级结构预测，结果显示，二级结构中有a-螺旋(H)89个、延伸片段(E)37个担转角(T)15个和无规则卷曲(C)247个，分别占总蛋白的22.940o9.540o,3.87%和63.660o。当H}450o}E}S0o时该蛋白结构为all-alpha型;当H}50o}E}45%时为all-beta型;当H}300o}E}20%时为alpha-beta型;其他情况为mixed型。预测结果表明，TabZIPl蛋白二级结构为mixed型。

2.1.3 TabZIPl蛋白多序列比对结果

利用NCBI中的Blastp在线工具进行其他物种的同源氨基酸序列搜索，将TabZIPl与大麦(H}bZIPl)、水稻(OsZIP-la)、短药野生稻(()bH-BP-la-like)、玉米(ZmbZIPLOC100192709)、苹果(MdbZIPLOC103455650)的5个bZIP蛋白进行多序列比对，结果表明，TabZIPl蛋白N端具有完整的碱性氨基酸序列特征结构域，C端是典型的亮氨酸拉链域。

2.1.4 进化树

利用MEUA6.0软件将小麦与玉米(Zeamays)、拟南芥(llrabido}sisthaliana)、烟草(Nicotianatabacum)、大麦(Hordeu。vuln“二subsp.vuln“二)、葡萄CViti、二。fera)、二穗短柄草(Brachypodiumdistachyon)、水稻(Ury}asat:二“)、短药野生稻((}ry}aBrachyantha)、刚毛怪柳(Tarna:二hispida)、苹果(Malusdomesti-、)等其他物种的bZIP的碱基序列进行系统发育分析(图2)，结果表明，TabZIPl与HvbZIPlBdHBP-la-like、ObHBP-la-like、OsZIP-la和ZmLOC100192709蛋白比较相近，即这些单子叶植物的bZIP聚类关系较近，其中小麦与大麦聚到了一组，水稻与短药野生稻聚到一组，而二穗短柄草在进化上介于麦类和水稻之间。在双子叶植物中，拟南芥和烟草聚到了一组，而葡萄与其他几种植物的聚类关系最远。

2.1.5 TabZIPl蛋白的亚细胞定位预侧

蛋白质要参与正常的生命活动，就必须处于特定的细胞区域。胞质中合成的细胞器蛋白质与胞质蛋白质一样都由游离核糖体合成，随后必须转运到特定的细胞器中。转运进入线粒体及叶绿体中的许多蛋白质都含有特定的导向序列(tar-getingsequence)，负责细胞器外膜对蛋白质的识别。导向序列又称转运肤(transitpeptide)或导肤(targetingpeptide)，是蛋白质N一端一段约20}80个氨基酸的肤链，富含带正电荷的碱性氨基酸(Arg和Lys)和轻基氨基酸(Ser)，易形成。一螺旋结构。以是否具有转运肤序列为特征，利用TargetP和ChloroP在线程序对TabZIPl蛋白预测结果显示，TabZIPl定位于线粒体，具有mTP(mitochondrionchloroplasttransitpeptide，线粒体转运肤)。

2.2 TabZIP1在不同组织中的表达模式分析

利用RT-PCR分析4个温度下东农冬麦1号和济麦22分孽节及叶片中TabZIPl的表达量特点。结果显示，济麦22分孽节中TabZIPl表达量随着温度的降低表现出“升一降一升”的趋势，且在0`C时达到最高，叶片表达量呈现出先升后降的趋势，在一10`C时最高;东农冬麦1号分孽节中TabZIPl表达量先有所降低再增加，一10`C时达到峰值，叶片表达量表现为明显的先增后降趋势，峰值同样出现在一10`C时。济麦22分孽节中TabZIPl相对表达水平在5`C及0`C时显著高于叶片，率先对低温响应;东农冬麦1号叶片中TabZIPl相对表达水平在0`C时显著高于分孽节，对低温迅速响应，但在一10`C时却是分孽节的TabZIPl相对表达水平显著高于叶片。总体来说，济麦22分孽节及东农冬麦1号叶片中的TabZIPl表达相比东农冬麦1号分孽节均在0`C时率先对低温响应，而东农冬麦1号分孽节则在一10`C时积极响应低温，表明东农冬麦1号分孽节是重要的抗寒器官。

2.3 ABA对东农冬麦1号分集节TabZIP1表达模式的影响

由于低温下TabZIPl在东农冬麦1号分孽节中的表达量较高，所以利用RT-PLR分析低温胁迫下东农冬麦1号分孽节中TabZIPl对外源ABA的反应。RT-PLR结果显示，5`C和0`C时，对照的TabZIPl相对表达量高于ABA处理，但差异不显著;一10`C和一25`C时，ABA处理的TabZIPl表达量均显著高于对照(图5)。

随着低温的降低，ABA处理组分孽节的TabZIPl表达量线性升高，且一25`C时达到最高。可见，外源ABA对TabZIPl在低温(一10`C以下)下的表达有诱导作用，且表现为正调控。

3 讨论

随着分子生物学的快速发展，人类已经完成了拟南芥、水稻等多种植物的基因组测序工作，科研工作者也对基因组中转录因子做出了预测和统计，转录因子的数量在植物基因组中占有相当比重，可见其在植物生命过程中的重要性。近年来，越来越多的bZIP转录因子被揭示参与调节体内的多种生命现象。以往研究表明，bZIP类转录因子在植物组织中有表达差异，如许多在根中高量表达，在茎和叶片中却表达量很少或不表达。本研究结果显示，东农冬麦I号TabZIPl在分孽节和叶片中的表达量随着温度的降低而逐渐升高，且分孽节高于叶片，分孽节表达量在一I}时达到峰值，这与实验室前期证明一I}是冬小麦启动抗寒机制的关键温度点相一致;低温下济麦叶片中的TabZIPl表达量也升高，最高表达量出现在一I}时，而济麦分孽节TabZIPl表达量在低温下较低，且在oC时达到峰值，较东农冬麦I号提前应答低温胁迫。

ABA作为一个关键的调节因子，参与植物的各种生物和生理过程。通过ABA调控基因的启动子研究已经鉴定出很多ABA响应元件，如ABRE元件。很多ABRE元件包含1个ACUT核心序列，如小麦E。基因上游的Emla元件(UUACACUTUUC)，而其他的ABRE元件则不包含ACUT核心序列，如大麦HI几l1基因上游的CE3元件(ACUCUTUTCCTC)、水稻rab16B上游的“motifIII"(GCCGCGTGGC)和人工合成元件hex-3(GGACGCGTGGC)。研究表明，AtbZIPl转录因子基因通过与ABRE元件相结合，参与ABA信号传导通路。在极端环境条件下，有效的抗氧化系统是植物抵抗逆境的一个重要途径。研究发现，烟草在逆境胁迫条件下过量表达转基因ThbZIPl基因，可以提高过氧化物酶、超氧化物歧化酶的活性，增加可溶性糖和可溶性蛋白的含量，还可以加快活性氧的消除，促进可溶性渗透剂的积累，诱导或增加可溶性蛋白的生物合成。在正常情况下外源ABA可造成类似逆境的胁迫，而胁迫条件下外源ABA则可提高植物自身防御机制.I}IU等研究表明，在一I}及以上低温时外源ABA促进了分孽节中可溶性糖和淀粉的积累。外源ABA可通过增加可溶性糖含量提高小麦、甜菜、水稻等抗寒性。本实验结果表明，外源ABA促进低温(一10`C及以下温度)环境中东农冬麦1号分孽节中TabZ-IPI的表达，呈现正调控。但在5`C和0`C时ABA调控TabZIPl的机制并未启动，这可能是5℃和o `c时ABA处理过的东农冬麦1号分孽节Ta bZIPl还未对寒冷做出应答，这两个温度不足以启动ABA对Ta bZIPl的调控，而是在一10 `C及以下温度启动，这与实验室前期结果一10 `C是东农冬麦1号启动抗寒机制的关键温度点相一致。对于Ta bZIPl是否参与强抗寒冬小麦东农冬麦1号可溶性糖含量的提高、淀粉的积累及一些逆境相关酶活性的调节，有待进一步研究。