茶树CsDREB-A2转录因子基因的克隆与非生物胁迫响应分析
韩妙华, 滕瑞敏, 李辉, 刘昊, 林士佳, 庄静*
南京农业大学园艺学院/茶叶科学研究所,江苏 南京 210095
*通讯作者:庄静,女,教授,主要从事茶树分子生物学研究。E-mail:zhuangjing@njau.edu.cn

作者简介:韩妙华,女,主要从事茶树分子生物学研究。E-mail:2018104089@njau.edu.cn

摘要

干旱应答元件结合蛋白(DREB)类转录因子在植物逆境信号转导途径中具有重要的调控作用。为了解AP2/ERF转录因子在茶树逆境胁迫的分子调控机理,本研究从茶树龙井43叶片的cDNA中克隆得到一个编码CsDREB-A2转录因子的基因;对 CsDREB-A2基因及其编码蛋白序列特征进行分析,并利用实时荧光定量PCR法检测该基因在茶树不同非生物胁迫处理下的表达水平。结果表明, CsDREB-A2基因开放阅读框为1 056 bp,编码351个氨基酸,其编码氨基酸序列具有AP2保守结构域,包含典型的YRG元件和WLG基序。AP2结构域第14、第19位氨基酸分别为缬氨酸和谷氨酸。拟南芥AP2/ERF家族转录因子的进化分析表明,该转录因子属于DREB亚族的A2组。CsDREB-A2相对分子质量为39 080 Da,理论等电点为5.32,属于亲水性蛋白,主要由α-螺旋和随机卷曲组成,无序化特征明显且存在一个LM无序区域;可能定位于细胞核,不存在信号肽和跨膜结构,属于非分泌蛋白。 CsDREB-A2基因在高温(38℃)和干旱(200 g·L-1 PEG)胁迫下均能快速诱导表达,并显著高于对照,分别在处理4、2 h达到最大值,为对照的20.70和42.90倍,植物在渗透胁迫下,可能通过ABA信号途径调节该基因对干旱的耐受性。高盐(200 mmol·L-1 NaCl)胁迫下 CsDREB-A2基因的表达受抑制,推测其可能存在负调控结构域降低该基因在盐胁迫下的表达量。本试验结果为研究DREB类转录因子在茶树抗逆胁迫的分子调控机制提供了一定的理论参考。

关键词: 茶树; CsDREB-A2; 转录因子; 非生物胁迫; 表达分析
文章编号:1000-8551(2020)12-2647-11
Cloning of CsDREB-A2 Transcription Factor Gene and It's Response to Abiotic Stress in Camellia sinensis
HAN Miaohua, TENG Ruimin, LI Hui, LIU Hao, LIN Shijia, ZHUANG Jing*
College of Horticulture, Nanjing Agricultural University/Tea Science Research Institute, Nanjing, Jiangsu 210095
Abstract

DREB transcription factors play an important role in the regulation of stress signal transduction in plants. Studying the response of DREB transcription factors to abiotic stress is helpful to understand the molecular regulation mechanism of AP2/ERF transcription factors under stress in tea plants. A gene encoding AP2/ERF transcription factor, CsDREB-A2, was cloned from the cDNA of Longjing43. The sequence characteristics of CsDREB-A2 were analyzed. The expression profiles of CsDREB-A2 gene under different abiotic stresses were detected by real-time quantitative PCR. The results showed that the open reading frame of CsDREB-A2 gene was 1 056 bp, encoding 351 amino acids, containing AP2 DNA-binding domain, and with typical YRG elements and WLG motifs. The 14th and 19th amino acids of AP2 domain is valine and glutamic acid, respectively. Evolutionary analysis of transcription factors of AP2/ERF family in Arabidopsis thaliana indicated that CsDREB-A2 transcription factor belonged to A2 group of DREB subgroup. CsDREB-A2 was a hydrophilic protein with a molecular weight of 39 080 Da and theoretical isoelectric point of 5.32. It was mainly composed of alpha-helix and random coil. The disorder characteristic was obvious, and there was a LM disorder region. It may be located in the nucleus, and there was no signal peptide and transmembrane structure, which belonged to non-secretory protein. The results of real-time quantitative PCR showed that CsDREB-A2 gene could be induced rapidly under high temperature (38°C) and drought (200 g·L-1 PEG) treatments. The expression level was significantly higher than that of the control, reaching the maximum value at 4 h and 2 h, respectively, which was 20.70 and 42.90 times of the control. It is speculated that the ABA signaling pathway regulates the drought tolerance of plants under osmotic stress. The expression of CsDREB-A2 gene was mainly inhibited under high salt treatment (200 mmol·L-1 NaCl). It was speculated that CsDREB-A2 gene might have a negative regulatory domain to reduce the expression of CsDREB-A2 gene under salt stress. This study provides potential theoretical references for studying the molecular regulation mechanism of DREB transcription factors in stress resistance of tea plants.

Keyword: Camellia sinensis; CsDREB-A2; transcription factor; abiotic stress; expression analysis

茶树[Camellia sinensis (L.) O. Kuntze]作为一种重要的经济作物, 在我国已有几千年的栽培历史。我国茶区分布较广, 地形复杂, 且茶树具有喜湿耐荫、喜温惧寒的特性, 在引种和栽培过程中, 盐碱、高温、干旱、强光等各种非生物胁迫严重影响茶树的生长发育, 降低茶叶的产量和品质[1]。非生物胁迫下, 转录因子胁迫相关蛋白合成, 在植物逆境信号转导中起关键作用[2]

转录因子是一种能够与DNA启动子区域相互作用的蛋白质, 抑制或激活其他基因转录表达[3]。干旱应答元件结合蛋白(DRE-binding protein, DREB)属于APETALA2/ethylene responsive foutor(AP2/ERF)转录因子家族, 主要参与植物多种生理生化反应的信号转导, 如对逆境(低温、干旱、盐渍)、激素、病害等的应答[4]。DREB转录因子存在1个保守的DNA结合域AP2, 能够与脱水响应元件(dehyaration responsive element, DRE)或CRT(C-repeat)等顺式作用元件结合, 部分成员能特异性地结合其基因上游GCC-box元件或耦合元件(coupling element, CE1), 从而促进或抑制低温、高盐、干旱等相关基因表达。DRE核心序列元件TACCGACAT调节脱落酸(abscisic acid, ABA)独立途径, 可以不依赖ABA在逆境胁迫下完成抗性调节[2]。DREB蛋白碱性核定位信号区位于N端, 酸性转录激活区位于C端, 中间是AP2结构域, 由57~70个氨基酸残基组成, 含有RAYD和YRG保守序列块。RAYD元件存在一个由18个氨基酸残基组成的高度保守的核心序列, 形成1个双亲性α -螺旋结构, 可能参与AP2/ERF转录因子与其他转录因子或DNA的相互作用; YRG元件具有保守的YRG氨基酸基序。AP2结构域第2个β -折叠中第14和第19位氨基酸残基存在差异, 决定了AP2/ERF转录因子能够与不同顺式作用元件特异性结合。

植物中DREB转录因子最早在拟南芥中分离得到[5]。Liu等[6]从拟南芥的cDNA文库中克隆了5个与DRE元件结合的转录因子(DREB1A、DREB1B、DREB1C、DREB2A以及DREB2B)。随后, 在水稻[7]、棉花[8]、胡萝卜[9]、大白菜[10]、巴哈雀稗[11]等物种中均克隆得到DREB基因。在拟南芥中, DREB亚族被分为6个亚组(A1~A6), DREB1/CBF和DREB2蛋白分别归入A1和A2亚组[12]DREB-A1亚组基因响应低温诱导; DREB-A2亚组与植物在高盐、干旱下的逆境耐受性有关, 其他4个亚组基因响应其他非生物胁迫诱导。已有研究证明, 水稻、拟南芥、烟草DREB基因过表达可能使种子在发芽期间降低对高盐的敏感性, 并促进根的早期生长[13]。过表达菊花DmDREBaDmDREBb能够提高转基因烟草耐寒性, DmDREBb还可加强转基因植物在高盐下的耐受性[14]。齿肋赤藓ScDREB10过表达显著提高拟南芥种子在渗透和盐胁迫下的发芽率, 提高了活性氧的清除能力[15]

目前, 茶树中DREB转录因子的调控网络和作用机理尚不清楚, 因此, 研究茶树CsDREB-A2基因在非生物胁迫下的表达情况, 对于进一步了解CsDREB转录因子的功能具有重要的作用。本研究以茶树龙井43为材料, 克隆得到1个编码CsDREB-A2转录因子的基因, 对该基因编码的氨基酸序列进行了序列比对、亚细胞定位、理化性质、亲水性和疏水性等方面的分析, 并研究了其进化树、二级和三级结构。通过实时荧光定量PCR技术研究了茶树CsDREB-A2基因在低温、高温、干旱和高盐胁迫下的响应情况, 旨在为进一步研究茶树CsDREB基因对逆境胁迫的响应机理提供一定的参考。

1 材料与方法
1.1 试验材料及处理

供试材料龙井43(Camellia sinensis cv. Longjing 43)二年生茶树扦插幼苗, 种植于南京农业大学茶叶科学研究所。摘取龙井43无病虫害植株的幼嫩叶片为试验材料。

选取长势良好的龙井43盆栽幼苗, 分别进行高温(38℃)、低温(4℃)、高盐(200 mmol· L-1 NaCl)和干旱[200 g· L-1 聚乙二醇(polyethy lene glyco, PEG)]胁迫处理, 处理时间分别为2、4、8和24 h, 以未处理的植株作为对照(CK), 在各时间段采集不同处理的茶树叶片。取样时, 每处理从50株龙井43茶树群体中随机选取5株健康植株, 每株摘取一芽二叶初展叶片, 全部进行混样。每个处理设置3个生物学重复。

1.2 试验方法

1.2.1 总RNA提取和cDNA合成 茶树叶片总RNA提取方法参照北京华越洋生物科技有限公司Quick RNA Isolation Kit试剂盒说明书。RNA样品浓度采NanoDrop ND 1000微量紫外检测仪(上海谱元仪器有限公司)测定, 采用1.2%琼脂糖凝胶电泳检测RNA质量[16, 17]。按照大连TaKaRa公司Prime Script RT Reagent Kit试剂盒说明书将总RNA反转录成cDNA模板。

1.2.2 CsDREB-A2基因的克隆 基于南京农业大学茶树逆境生理课题组转录组数据[18], 设计1对特异引物, 上游引物CsDREB-A2-QF序列为5'-A T G A T A G C T C A G A A C A T T GA-3', 下游引物CsDREB-A2-QR序列为5'-T T A A A G A C C C A A A T C T G A G TA-3'。PCR扩增体系为20 μ L, 包含7 μ L ddH2O、10 μ L 2× Taq Plus Master Mix酶、1 μ L cDNA模板、上下游引物各1 μ L。PCR扩增程序:95℃预变性5 min; 94℃变性30 s, 55℃退火30 s, 72℃延伸1 min, 共35个循环; 最后72℃延伸10 min[19]。用1.2%琼脂糖凝胶电泳分离PCR产物, 参照杭州维特洁生化技术公司DNA回收试剂盒说明书将其回收, 与pMD 19-T载体连接, 转化至大肠杆菌DH 5α 中。挑取阳性克隆菌液送至南京金斯瑞公司测序。

1.2.3 生物信息学分析 在NCBI网站(http://www.ncbi.nlm.nih.gov)预测序列保守域, 并用BLAST进行序列同源性比较[20]。利用DNAMAN 8.0软件进行氨基酸序列多重比对和亲/疏水性分析。登陆SMS(http://www.bio-soft.net/sms)和ExPASy-ProtParam(http://web.expasy.org/protparam)分析氨基酸序列组成及理化性质[21]。利用FoldIndex网站(https://fold.weizmann.ac.il/fldbin/findex)分析氨基酸折叠无序化特性。蛋白亚细胞定位采用LocTree3(https://rostlab.org/services/loctree3)和SoftBerry ProtComp 9.0(http://linux1.softberry.com/berry.phtml)预测。拟南芥AP2/ERF转录因子家族蛋白全长序列登陆Tair(https://www.arabidopsis.org)下载, 并利用MEGA 7.0软件构建系统进化树。采用NetPhos 3.1(http://www.cbs.dtu.dk/services/NetPhos/)预测磷酸化位点。跨膜结构和信号肽分别利用TMpred(https://embnet.vital-it.ch/software/TMPRED_form.html)和SignalP4.1 (http://www.cbs.dtu.dk/services/SignalP-4.1)进行预测。蛋白质二级结构预测采用SOPMA在线网站(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)进行。登陆Swiss-Model(http://www.Swissmodel.expasy.org)预测蛋白质三级结构[19], 并利用PyMOL软件生成三级结构模型。

1.2.4 茶树CsDREB-A2基因的表达特性分析 提取非生物胁迫处理样品的RNA, 反转录成cDNA, 用于实时荧光定量PCR。按照大连TaKaRa公司SYBR Premix Ex Taq Kit试剂盒说明书以及CFX96TM real-time PCR system平台完成荧光定量, 分析CsDREB-A2基因在高温、低温、干旱、高盐处理下茶树中的表达情况。茶树GAPDH基因作为内参基因[22]。检测目标基因表达的引物为CsDREB-A2-qrtF(5'-T T A A A G A C C C A A A T C T G A G TA-3')和CsDREB-A2-qrtR(5'-T T A A A G A C C C A A A T C T G A G TA-3')。内参基因的引物为CsGAPDH-JF(5'-T T G G C A T C G T T G A G G G T CT-3')和CsGAPDH-JR(5'-C A G T G G G A A C A C G G A A A GC-3'), 采用2-Δ Δ CT[23]分析结果。采用IBM SPSS Statistic 25和Micosoft Office Excel 2010进行数据分析并制作图表。

2 结果与分析
2.1 茶树CsDREB-A2基因的克隆

以茶树叶片cDNA为模板, 通过引物CsDREB-A2-QF和CsDREB-A2-QR进行PCR扩增, 得到1 000 bp左右的扩增片段。对扩增产物进行序列测定与分析, 结果显示, 茶树CsDREB-A2基因开放阅读框为1 056 bp, 共编码351个氨基酸(图1), GenBank登录号为MN082027。

图1 茶树CsDREB-A2基因核苷酸序列与推测的氨基酸序列Fig.1 The nucleotide acid and deduced amino acid sequence of CsDREB-A2 gene in C. sinensis

2.2 茶树CsDREB-A2转录因子进化树分析

茶树CsDREB-A2转录因子与拟南芥AP2/ERF家族成员进化分析结果显示, 茶树CsDREB-A2属于DREB亚族的A2组(图2-A)。将其氨基酸序列进行BLAST比对发现, CsDREB-A2 转录因子与猕猴桃(Actinidia chinensis)、葡萄(Vitis vinifera)等物种有较近的亲缘关系, 与梅花(Prunus mume)、甜樱桃(Prunus avium)等物种亲缘关系较远(图2-B)。

图2 茶树CsDREB-A2转录因子进化树(A)及与部分物种氨基酸序列进化分析(B)Fig.2 Phylogenetic tree of CsDREB-A2 from C. sinensis (A) and phylogenetic analysis of amino acid sequences with some species (B)

2.3 茶树CsDREB-A2转录因子与其他植物相关氨基酸序列比对

茶树CsDREB-A2转录因子保守域预测结果显示(图3-A), AP2保守结构域位于该转录因子第9~第154氨基酸位点之间, N端含有典型的保守YRG元件和WLG基序, 其中AP2保守结构域第14、第19位氨基酸分别为缬氨酸(V)和谷氨酸(E)。第254~第263氨基酸位点之间存在1个LM无序化区域(FDxxELLxxL)。将茶树CsDREB-A2 氨基酸序列与其他物种DREB类转录因子氨基酸序列进行BLAST比对, 结果显示, 茶树与川桑(Morus notabilis)、木豆(Cajanus cajan)、马铃薯(Solanum tuberosum)等物种的DREB类转录因子相似性为52.56%(图3-B)。

图3 茶树CsDREB-A2转录因子保守域(A)与其他物种氨基酸序列的多重比对(B)Fig.3 Conserved domains of CsDREB-A2 transcription factor (A) and multiple alignment of relative amino acid sequences in C. sinensis and other species (B)

2.4 茶树CsDREB-A2转录因子氨基酸组成及理化性质分析

表1可知, DREB类转录因子蛋白相对分子量在33.78~44.25 kDa之间, 蛋白残基数为304~402, 理论等电点为4.63~5.32。不同植物中各氨基酸种类所占比例不同, 其中酸性氨基酸所占比例普遍高于碱性氨基酸, 脂肪族氨基酸高于芳香族氨基酸; 蛋白质总平均疏水性为-0.77左右, 表现出较强的亲水性。利用NetPhos 3.1对CsDREB-A2蛋白磷酸化位点进行预测, 结果表明, 该多肽链0.5以上分值的氨基酸位点为47个, 其中包含29个丝氨酸(S)、12个苏氨酸(T)和6个酪氨酸(Y)磷酸化位点。

表1 不同植物中DREB氨基酸组成成分及理化性质分析 Table 1 Analysis of composition and physicochemical properties of DREB amino acids in different plants
2.5 茶树CsDREB-A2蛋白亚细胞定位预测

CsDREB-A2蛋白亚细胞定位预测结果显示, CsDREB-A2蛋白主要定位于细胞核(表2)。在线预测网站LocTree3对该蛋白的亚细胞定位预测分析结果与SoftBerry ProtComp 9.0一致, 由此推断, 茶树CsDREB-A2蛋白主要在细胞核中发挥生物学作用。

表2 CsDREB-A2转录因子亚细胞定位预测 Table 2 The prediction of subcellular localization of CsDREB-A2 transcription factor
2.6 茶树CsDREB-A2蛋白氨基酸亲、疏水性和无序化分析

CsDREB-A2蛋白氨基酸序列的亲水性和疏水性分析结果显示, 该转录因子疏水性最强的位点是第351位亮氨酸(Leu), 其次是第261位甘氨酸(Gly); 亲水性最强的是第34位精氨酸(Arg), 其次第33位赖氨酸(Lys)、第35~第37位赖氨酸、丝氨酸(Ser)、精氨酸; 大部分氨基酸属于亲水性氨基酸(图4)。利用FoldIndex对CsDREB-A2转录因子进行氨基酸序列折叠无序化分析, 结果表明, 该转录因子存在5个氨基酸无序化区域, 共包含258个氨基酸, 无序化比例为73.50%, 无序化特征明显(图5)。

图4 茶树CsDREB-A2氨基酸序列疏水性和亲水性分析Fig.4 Hydrophobicity and hydrophilicity of amino acid sequence of CsDREB-A2 in C. sinensis

图5 茶树CsDREB-A2转录因子折叠状态的分析
注:正值表示蛋白的有序状态, 负值表示蛋白的无序状态; 下划线标注为有序状态氨基酸, 无下划线标注为无序状态氨基酸。
Fig.5 Analysis of folding state of CsDREB-A2 transcription factor in C. sinensis
Note:The positive value indicates the ordered state of the protein, while the negative value indicates the disordered state of the protein. The underline indicates the ordered state of the amino acid, and the ununderlined indicates the disordered state of the amino acid.

2.7 茶树CsDREB-A2转录因子信号肽、跨膜结构预测与分析

经信号肽预测分析可知(图6), CsDREB-A2转录因子mean S-score的值小于0.5, 表明CsDREB-A2蛋白不存在信号肽序列, 可能不属于分泌蛋白。利用Tmpred在线网站进行跨膜结构预测, 结果显示, CsDREB-A2蛋白在269~285氨基酸位点之间存在1处由内至外螺旋, 分值为88, 远低于500, 不存在由外至内螺旋, 由此推测CsDREB-A2转录因子属于非跨膜蛋白, 没有跨膜区域。

图6 茶树CsDREB-A2转录因子信号肽预测
注:C-score:原始切割位点分值; S-score:信号肽分值; Y-score:综合切割位点分值。
Fig.6 Signaling peptide prediction of CsDREB-A2 transcription factor in C. sinensis
Note:C-score:Raw cleavage site score. S-score:signal peptide score.Y-score:Combined cleavage site score.

2.8 茶树CsDREB-A2蛋白推导的二级和三级结构预测与分析

由图7可知, 该转录因子的组成为30.77% α -螺旋、3.99% β -折叠、8.83%延伸主链及56.42%随机卷曲。因此, CsDREB-A2主要由α -螺旋和随机卷曲组成, β -折叠散布于蛋白序列中。

图7 茶树CsDREB-A2转录因子二级结构预测
注:蓝色:α -螺旋; 红色:β -折叠; 绿色:β -转角; 粉色:随机卷曲。
Fig.7 Prediction of secondary structure of CsDREB-A2 transcription factor in C. sinensis
Note:Blue:Alpha helix. Red:Extended strand. Green:Beta turn. Pink:Random coil.

采用Swiss-Model对CsDREB-A2 转录因子进行三级结构建模分析, 结果表明, 其AP2结构域含有1个α -螺旋和3个β -折叠, 且该结构域第14位的缬氨酸和第19位的谷氨酸位于第2个β -折叠中, WLG基序存在于第3个β -折叠(图8)。该转录因子三级结构以随机卷曲和α -螺旋为主, 推测结果与二级结构相吻合。

图8 茶树CsDREB-A2的三级结构预测Fig.8 The three-dimension structures of CsDREB-A2 in C. sinensis

2.9 茶树CsDREB-A2基因在非生物胁迫下的表达情况

茶树CsDREB-A2基因在4种胁迫下均能诱导表达, 但是不同胁迫及不同处理时间下CsDREB-A2基因表达存在差异。在低温胁迫下, CsDREB-A2基因的相对表达量呈先降低后上升随后又降低的趋势, 在处理8 h时其相对表达量达到最大值, 为CK的5.61倍, 处理4、24 h时CsDREB-A2基因的相对表达量也明显高于CK, 分别为CK的1.74和4.39倍(图9-A)。高温和干旱胁迫下, CsDREB-A2基因相对表达量较CK表现出明显的上调, 分别在处理4 h和2 h时达到最大值, 为CK的20.70倍和42.90倍(图9-B、C), 其中高温处理2、8、24 h时CsDREB-A2基因的相对表达量分别为CK的18.23、12.65、13.12倍; 干旱处理8、24 h时CsDREB-A2基因的相对表达量显著高于CK, 分别为CK的24.05、39.98倍。高盐胁迫下茶树CsDREB-A2基因的相对表达量在处理2、4、8 h表现出下调趋势, 处理2 h时其相对表达量最低, 在处理24 h时其相对表达量显著高于CK, 为CK的2.05倍(图9-D)。上述结果表明, 茶树CsDREB-A2基因参与高温、低温、高盐、干旱等非生物胁迫应答过程, 对不同处理的响应模式存在差异。

图9 不同胁迫下茶树CsDREB-A2基因在叶片中的相对表达量
注:A:低温; B:高温; C:干旱; D:高盐。不同小写字母表示在0.05水平上差异显著。
Fig.9 Relative expression of CsDREB-A2 gene in leaves of C. sinensis under different stress treatments
Note:A:Low temperature. B:High temperature. C:Drought. D:High salinity. Different lowercase letters indicate significate differences at 0.05 level.

3 讨论

本研究对比了茶树CsDREB-A2基因在低温(4℃)、高温(38℃)、干旱(200 g· L-1 PEG)和高盐(200 mmol· L-1 NaCl)胁迫下的表达特性, 结果显示, 高温和干旱胁迫均能快速诱导CsDREB-A2基因表达。推测该转录因子参与茶树在高温、干旱胁迫下的应答过程。已有研究发现, 拟南芥叶片AtDREB2A不响应低温胁迫, 但受高温胁迫诱导表达[24]。大豆GmDREB2A基因在高温、干旱胁迫下显著上调[25]。在热激反应中, 热激引发拟南芥[26]和水稻[27]热激转录因子(heat shock transcription factor A1, HsfA1)激活, 诱导DREB2AHsfA2基因的转录; 从而使得DREB2A与其靶基因HsfA3的启动子结合, 激活HsfA3基因, 进一步诱导下游热激蛋白编码基因的表达[28]。以上结果进一步推测, CsDREB-A2可能参与高温胁迫下茶树转录水平的热激级联反应, 进而调控茶树应答高温和干旱胁迫。

拟南芥DREB2A蛋白中心区域包含一个负调控结构域(NRD), 影响DREB2A转录因子的结合和激活功能[29]。拟南芥AtDREB2A、大豆GsDREB2和向日葵HaDREB2A基因中的NRD结构域都含有一段调控蛋白质降解的PEST四肽序列(脯氨酸-谷氨酸-丝氨酸-苏氨酸)[30]。才华等[29]研究表明, 野生大豆GsDREB2基因的负调控结构域抑制GsDREB2的转录激活功能, 导致GsDREB2不能与DRE元件特异性结合。人工突变野生大豆GsDREB2基因的负调控结构域可获得GsDREB2-mNRD基因, 超量表达该基因的拟南芥显著提高对盐胁迫和渗透胁迫的耐受性[29]。本研究茶树CsDREB-A2转录因子AP2结构域和转录激活区存在丝氨酸(Ser)富集区域, 推测CsDREB-A2基因可能存在NRD结构域, 降低该基因对盐胁迫的调节功能, 使得CsDREB-A2基因在高盐胁迫下表现出下调趋势。ABA激活许多干旱胁迫基因的表达, 这些基因启动子区包含保守的ABA顺式反应元件ABRE[31]。植物中, DREB类转录因子基因启动子区域通常含有ABRE元件[30, 31]。但茶树CsDREB-A2转录因子对非生物胁迫响应过程中ABA的作用机制尚需要进一步深入研究。

本研究中, 茶树CsDREB-A2氨基酸序列折叠无序化分析显示, 该转录因子无序化比例为73.50%, N-末端和C-末端固有无序化程度较高, 这与Kragelund等[32]的研究一致。固有无序蛋白(intrinsic disordered proteins, IDPs)在天然状态下不存在稳定的三维结构, 但在行使生物学功能上具有重要作用。正常生长条件下, DREB2A基因表达不会激活下游基因, DREB2A需要翻译后修饰才能被激活, 但是其激活机制尚不明确。无序结构使IDPs比较灵活, 含有丰富的结合位点, 可对翻译后修饰, 如甲基化、磷酸化、乙酰化等进行调控, 同时在细胞生物信号多样性功能方面起作用[33]。DREB2A转录因子中的LM区域可能介导该转录因子与自由基诱导细胞死亡蛋白1(radical induced cell death protein 1, RCD1)之间的相互作用[33]。N-末端部分无序区域与DREB2A相互作用蛋白1(DREB2A-interacting protein1, DRIP1)和DRIP2在细胞核发生互作, 负调节干旱胁迫应答反应[32]。因此, 这些固有无序蛋白主要参与分子识别、调控以及信号转导等功能。但由于无序结构在茶树中的互作方式较为复杂, 其具体反应过程有待进一步研究。

4 结论

本研究克隆得到的茶树CsDREB基因属于AP2/ERF家族转录因子中DREB亚族的A2组, 与川桑、木豆、马铃薯等物种的DREB基因具有较高的同源性。CsDREB蛋白含有一个由60个氨基酸残基组成的高度保守的结构域AP2, YRG元件和WLG基序存在于该结构域中, 具有AP2/ERF家族转录因子的典型特征。该结构域第14和第19位氨基酸分别为缬氨酸和谷氨酸, 起到与DRE/CRT作用元件结合的重要作用。本研究证实高温和干旱胁迫可诱导茶树CsDREB-A2基因表达, 表明DREB类转录因子在茶树非生物胁迫过程中具有重要的调节作用。本试验结果为DREB类转录因子在茶树抗逆胁迫的分子调控机制研究奠定了一定的理论基础。

参考文献
[1] 曹丹, 马林龙, 金孝芳, 龚自明. 茶树抗非生物逆境的研究进展[J]. 湖南农业科学, 2015(10): 152-154 [本文引用:1]
[2] Xie Z L, Nolan T M, Jiang H, Yin Y H. AP2/ERF transcription factor regulatory networks in hormone and abiotic stress responses in Arabidopsis[J]. Frontiers in Plant Science, 2019, 10: 1-17 [本文引用:2]
[3] Shou H, Bordallo P, Fan J B, Yeakley J M, Bibikova M, Sheen J, Wang K. From The Cover: Expression of an active tobacco mitogen-activated protein kinase kinase kinase enhances freezing tolerance in transgenic maize[J]. Proceedings of the National Academy of Sciences, 2004, 101(9): 3298-3303 [本文引用:1]
[4] Agarwal P K, Gupta K, Lopato S, Agarwal P. Dehydration responsive element binding transcription factors and their applications for the engineering of stress tolerance[J]. Journal of Experimental Botany, 2017, 68(9): 2135 [本文引用:1]
[5] Stockinger E J, Gilmour S J, Thomashow M F. Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit[J]. Proceedings of the National Academy of Sciences, 1997, 94(3): 1035-1040 [本文引用:1]
[6] Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought-and low-temperature-responsive gene expression, respectively, in Arabidopsis[J]. Plant Cell, 1998, 10(8): 1391-406 [本文引用:1]
[7] Dubouzet J G, Sakuma Y, Ito Y, Kasuga M, Dubouzet E G, Miura S, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. OsDREB genes in rice, Oryza sativa L. encode transcription activators that function in drought-, high-salt-and cold-responsive gene expression[J]. The Plant Journal, 2003, 33(4): 751-763 [本文引用:1]
[8] 李付振, 邱新棉, 刘传亮. 棉花中一个新的类DREB转录因子(GhDREB2B)的克隆、序列特征及表达分析[J]. 浙江农业学报, 2010, 22(5): 564-569 [本文引用:1]
[9] 黄莹, 徐志胜, 王枫, 宋雄, 田畅, 熊爱生. 胡萝卜2个DcDREB-A1类转录因子基因的克隆与比较分析[J]. 核农学报, 2015, 29(1): 29-39 [本文引用:1]
[10] Li M Y, Wang F, Jiang Q, Li R, Ma J, Xiong A S. Genome-wide Analysis of the distribution of AP2/ERF transcription factors reveals duplication and elucidates their potential function in Chinese Cabbage ( Brassica rapa ssp. pekinensis)[J]. Plant Molecular Biology Reporter, 2013, 31(4): 1002-1011 [本文引用:1]
[11] 张宇君, 王普昶, 赵丽丽, 曾庆飞, 张雄, 陈超, 董瑞. 巴哈雀稗 PnDREB2基因克隆及时序表达图谱构建[J]. 核农学报, 2019, 33(1): 7-9 [本文引用:1]
[12] Eckardt N A. DREB duo defines distinct drought and cold response pathways[J]. Plant Cell, 2019, 31(6): 1196-1197 [本文引用:1]
[13] 毛彦芝. DREB转录因子与植物非生物胁迫抗性研究进展[J]. 农业科技通讯, 2014(3): 35-37 [本文引用:1]
[14] 武剑, 杨艳芳, 王慧, 朱凯, 刘黎卿, 陈发棣, 喻德跃. 菊花 DmDREBa/b基因提高转基因烟草耐逆性研究[J]. 分子植物育种, 2016, 14(11): 3063-3072 [本文引用:1]
[15] Liang Y Q, Li X S, Zhang D Y, Gao B, Yang H L, Wang Y C, Guan K Y, Wood A J. ScDREB8, a novel A-5 type of DREB gene in the desert moss Syntrichia caninervis, confers salt tolerance to Arabidopsis[J]. Plant Physiology and Biochemistry, 2017, 120: 242-251 [本文引用:1]
[16] 王文丽, 吴致君, 刘志薇, 王永鑫, 李辉, 崔新, 庄静. 茶树类黄酮3'-羟化酶基因的克隆与表达特性分析[J]. 茶叶科学, 2017, 37(1): 108-117 [本文引用:1]
[17] 刘志薇, 熊洋洋, 李彤, 严雅君, 韩洪润, 吴致君, 庄静. 茶树转录因子基因 CsDREB-A4的克隆与温度胁迫响应的分析[J]. 茶叶科学, 2015, 35(1): 24-26 [本文引用:1]
[18] Wu Z J, Li X H, Liu Z W, Xu Z S, Zhuang J. De novo assembly and transcriptome characterization: Novel insights into catechins biosynthesis in Camellia sinensis[J]. BMC Plant Biology, 2014, 14(1): 277 [本文引用:1]
[19] 刘昊, 王文丽, 滕瑞敏, 李辉, 王瑜, 王永鑫, 庄静. 茶树蛋白激酶基因 CsCIPK的克隆与表达特性分析[J]. 园艺学报, 2019, 46(1): 101-111 [本文引用:2]
[20] 李彤, 严雅君, 韩洪润, 刘志薇, 吴致君, 庄静. 茶树品种‘安吉白茶’ CsDREB-A4 b基因的克隆及其对非生物胁迫的响应[J]. 植物资源与环境学报, 2016, 25(2): 1-9 [本文引用:1]
[21] 刘志薇, 熊洋洋, 李彤, 严雅君, 韩洪润, 吴致君, 庄静. 茶树中两个ERF类转录因子的分离及不同茶树中温度胁迫的响应分析[J]. 植物生理学报, 2014(12): 1821-1832 [本文引用:1]
[22] Wu Z J, Tian C, Jiang Q, Li X Y, Zhuang J. Selection of suitable reference genes for RT-qPCR normalization during leaf development and hormonal stimuli in tea plant ( Camellia sinensis)[J]. Scientific Reports, 2016, 6(1): 19748 [本文引用:1]
[23] Livak K J, Schmittgen T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT Method[J]. Methods, 2001, 25(4): 402-408 [本文引用:1]
[24] 张志飞, 杨知建, 周倩, 赵志丽. DREB2 s转录因子基因的表达调控机制研究进展[J]. 植物科学学报, 2014, 32(3): 297-303 [本文引用:1]
[25] Mizoi J, Ohori T, Moriwaki T, Todaka D, Maruyama K, Kusakabe K, Osakabe Y, Shinozaki K, Yamaguchi-Shinoz aki K. GmDREB2A; 2, a Canonical DEHYDRATION-RESPONSIVE ELEM ENT-BINDING PROTEIN2-type transcription factor in soybean, is posttranslationally regulated and mediates dehydration-responsive element-dependent gene expression[J]. Plant Physiology, 2013, 161(1): 346-361 [本文引用:1]
[26] Yoshida T, Ohama N, Nakajima J, Kidokoro S, Mizoi J, Nakashima K, Maruyama K, Kim J M, Seki M, Todaka D, Osakabe Y, Sakuma Y, Schöffl F, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis HsfA1 transcription factors function as the main positive regulators in heat shock-responsive gene expression[J]. Molecular Genetics and Genomics, 2011, 286(5/6): 321-332 [本文引用:1]
[27] Matsukura S, Mizoi J, Yoshida T, Todaka D, Ito Y, Maruyama K, Shinozaki K, Yamaguchi-Shinozaki K. Comprehensive analysis of rice DREB2-type genes that encode transcription factors involved in the expression of abiotic stress-responsive genes[J]. Molecular Genetics and Genomics, 2010, 283(2): 185-196 [本文引用:1]
[28] Schramm F, Larkindale J, Kiehlmann E, Ganguli A, Englich G, Vierling E, von Koskull-Döring P. A cascade of transcription factor DREB2A and heat stress transcription factor HsfA3 regulates the heat stress response of Arabidopsis[J]. The Plant journal, 2008, 53(2): 264-274 [本文引用:1]
[29] 才华, 孙娜, 宋婷婷. 人工改造野生大豆 GsDREB2基因对植物耐盐和耐渗透胁迫能力的影响[J]. 草业学报, 2018, 27(6): 168-176 [本文引用:3]
[30] June-Sik K, Mizoi J, Yoshida T, Fujita Y, Nakajima J, Ohori T, Todaka D, Nakashima K, Hirayama T, Shinozaki K, Yamaguchi-Shinozaki K. An ABRE promoter sequence is involved in osmotic stress-responsive expression of the DREB2 A gene, which encodes a transcription factor regulating drought-inducible genes in Arabidopsis[J]. Plant and Cell Physiology, 2011, 52(12): 2136-2146 [本文引用:2]
[31] Yoshida T, Fujita Y, Sayama H, Kidokoro S, Maruyama K, Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE‐dependent ABA signaling involved in drought stress tolerance and require ABA for full activation[J]. The Plant Journal, 2010, 61(4): 672-685 [本文引用:2]
[32] Kragelund B B, Jensen M K, Skriver K. Order by disorder in plant signaling[J]. Trends in Plant Science, 2012, 17(11): 625-632 [本文引用:2]
[33] Xie H, Vucetic S, Iakoucheva L M, Oldfield C J, Dunker A K, Obradovic Z, Uversky Ⅴ N. Functional anthology of intrinsic disorder. 3. ligand s, post-translational modifications, and diseases associated with intrinsically disordered proteins[J]. Journal of Proteome Research, 2007, 6(5): 1917-1932 [本文引用:2]