红麻非生物逆境胁迫响应基因HCWRKY71表达分析及转化拟南芥

李 辉 李德芳 邓 勇 潘 根 陈安国 赵立宁 唐慧娟

专题

红麻非生物逆境胁迫响应基因表达分析及转化拟南芥

李 辉1,2李德芳2,*邓 勇2潘 根2陈安国2赵立宁2唐慧娟2

1湖南文理学院生命与环境科学学院, 湖南常德 415000;2中国农业科学院麻类研究所, 湖南长沙 410205

WRKY转录因子在植物响应非生物逆境胁迫过程中具有重要的调控作用。本研究根据红麻转录组unigene序列(CL3883.Contig4), 设计引物进行PCR扩增, 经sanger测序获得全长为957 bp的基因cDNA序列。该基因开放读码框为957 bp, 编码1个含有318个氨基酸的蛋白, 具有1个WRKY功能保守结构域, 属于II类WRKY转录因子。在盐胁迫下,基因表达量随NaCl溶液浓度升高而增加; 在干旱胁迫下,基因表达量随干旱胁迫时间的增加呈现先下降再上升然后再下降的趋势; 在重金属镉胁迫下,基因表达量随CdCl2溶液浓度的增加而降低。表明该基因表达受盐、干旱和重金属镉胁迫的诱导。利用农杆菌介导花序浸染法将该基因转化拟南芥发现,基因提高了转基因拟南芥幼苗的耐盐性, 这为进一步研究基因的耐逆机制奠定了坚实的基础。

红麻; 耐盐; 干旱; 重金属镉;

红麻(L.)是一年速生草本纤维作物。随着人们环保意识的增强, 全球对于天然纤维的需求日益增加, 大力发展红麻产业是满足天然纤维需求的重要途径之一[1]。然而随着全球人口的增加, 人口、粮食、耕地之间的矛盾日益加剧, 红麻的种植必须转向盐碱地、干旱坡地、不能生产粮食的重金属污染土地。因此, 研究红麻的耐逆机制, 选育红麻耐逆品种对维持红麻产业的可持续发展具有重要意义。

高等植物在进化过程中形成了精细的调控网络来应对不同的逆境胁迫。转录因子如MYB[2]、WRKY[3]、NAC[4]在调控网络中发挥极其重要的作用, 它们通过激活/抑制下游基因表达来提高植物本身对逆境的适应能力。其中WRKY转录因子是植物体内特有的转录因子家族, 其特点是DNA结合域的N-端至少包含一段保守的WRKYGQK序列, C-端是一个类锌指结构域。DNA结合域通过与启动子中的W-box序列结合来调节靶标基因的表达[5]。

WRKY转录因子家族在植物响应干旱、盐、重金属、温度、病害等多种非生物/生物胁迫过程中具有重要的调控功能, 且不断有新的家族成员被鉴定[6]。Wu等[7]在菜豆基因组中发现了19个WRKY响应基因, 在干旱胁迫条件下, 其中11个表现为下调表达, 8个表现为上调表达。王岩岩等[8]研究表明, 野生大豆基因表达受干旱胁迫的诱导, 且随干旱胁迫时间的增加,基因表达呈现先上升后下降的趋势。柯丹霞等[9]研究表明, 大豆基因表达受盐胁迫诱导, 并且能够显著提高转基因百脉根的耐盐性。甘蔗基因在NaCl胁迫下显著上调表达, 推测该基因在甘蔗响应NaCl胁迫过程中具有一定的调控作用[10]。在低温胁迫条件下, 转基因的拟南芥与野生型相比, 具有较高的存活率。在低温胁迫和外源ABA诱导下,基因出现上调表达[11]。王影等[12]研究表明, 东南景天基因表达受镉胁迫的诱导, 且随镉胁迫时间的增加, 其表达量逐渐降低, 推测其在镉胁迫响应早期具有重要作用。以上研究证实, WRKY转录因子在植物抵御各种非生物胁迫过程中具有正向或负向的调控功能。但目前关于WRKY基因家族成员在红麻非生物逆境胁迫响应过程中的作用鲜见报道。

本课题组前期通过转录组测序获得unigene序列(CL3883.Contig4)[13], 通过NCBI的在线工具Blast分析发现, 其编码产物为转录因子WRKY71, 且该unigene受盐胁迫的诱导。因此, 本研究克隆了基因序列, 利用实时荧光定量PCR的方法分析了其在盐、干旱和重金属不同逆境胁迫下的时空表达特征, 并转化拟南芥进行初步功能验证, 旨在为进一步研究该基因耐盐的调控机制奠定基础。

1 材料与方法

1.1 试验材料

红麻16号由中国农业科学院麻类研究所提供。挑选籽粒饱满、大小均匀的红麻种子, 用0.1%的氯化汞消毒10 min, 蒸馏水冲洗3次, 然后播种于发芽盒中, 置光照培养箱中培养。培养10 d后, 选取植株健壮、长势均匀的红麻幼苗, 移植到1/2 Hoagland营养液中培养。红麻幼苗植株长到四叶一心时, 选取20株幼苗进行NaCl胁迫处理, 每个处理5株, 1/2 Hoagland营养液中NaCl的最终浓度分别为0、50、100、200 mmol L-1; 同样选取20株幼苗进行干旱处理, 将选取的幼苗根部暴露在空气中, 不进行供水处理, 处理时间为0、0.5、1.0、2.0、4.0 h, 每个处理5株; 同样选取20株幼苗进行重金属镉胁迫处理, 每个处理5株, 1/2 Hoagland营养液中氯化镉(CdCl2)的最终浓度分别为0、50、100、150 μmol L-1。分别取每个处理组红麻植株的根、茎、叶, 液氮速冻, 放置于-80℃冰箱备用。

1.2 红麻HcWRKY71基因的克隆及测序

利用植物总RNA提取试剂盒(DP432) (天根生化科技北京有限公司)提取200 mmol L-1NaCl处理的红麻植株叶片总RNA, 然后逆转录合成cDNA, 进行PCR扩增。根据课题组转录组unigene序列(编号: CL3883.Contig4), 利用PrimerPremier 5.0设计基因全长扩增引物HcWRKY71-F/R (表1)。PCR反应的体系包括2.5 μL 10× PCR buffer (Mg2+plus)、0.25 μL酶 (TaKaRa, 中国大连)、1 μL 2.5 mmol L-1dNTP、0.25 μL 10 μmol L-1HcWRKY71-F、0.25 μL 10 μmol L-1HcWRKY71-R、20 ng cDNA, ddH2O补足25 μL。PCR扩增的反应程序为94℃预变性2 min; 94℃变性30 s, 55℃退火30 s,72℃延伸60 s, 30个循环; 72℃延伸3 min, 4℃终止反应。反应产物经1.2%琼脂糖凝胶电泳检测。参照琼脂糖凝胶DNA回收试剂盒(天根生化科技北京有限公司)的说明书切胶回收, 然后将回收片段连接到pMD19-T载体上, 转化大肠杆菌DH5α感受态细胞, 挑取阳性菌落送生工生物工程(上海)股份有限公司测序。

表1 本研究所用引物

1.3 红麻HcWRKY71基因的生物信息学分析

利用NCBI在线程序ORF finder (http://www. ncbi.nlm.nih.gov/gorf/gorf.html)和NCBI Conservered Domain (https//www.ncbi.nlm.nih.gov/cdd)预测基因序列的最大开放读码框及其保守功能结构域; 运用在线工具(http://www.expasy.ch/tools/pi_tool.html)和NetPhos2.0 (http://www.cbs.dtu.dk/services/NetPhos/)分别预测基因编码蛋白质的分子量和等电点及其蛋白磷酸化位点; 利用在线工具BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi)对HcWRKY71的氨基酸序列进行在线序列比对; 利用DNAMAN软件进行基因氨基酸序列一致性比对; 利用MEGA6.0软件构建基因系统进化树。

1.4 红麻HcWRKY71基因在盐、干旱、重金属镉胁迫下的表达分析

红麻幼苗植株分别在NaCl和CdCl2胁迫下生长2周后, 剪取每个处理的红麻植株叶片, 进行浓度梯度的表达分析, 同时剪取红麻植株的根、茎、叶进行组织特异性表达分析。干旱胁迫条件下, 分别剪取不同处理时间的红麻植株叶片, 进行不同胁迫时间的表达分析; 红麻干旱处理8 h时, 剪取红麻植株的根、茎、叶, 进行组织特异性表达分析。每个处理取样后, 分别提取总RNA, 反转录合成cDNA, 利用荧光定量PCR检测基因表达量。荧光定量PCR引物为HcWRKY71-QF/QR (表1), 内参基因为, 内参引物为Actin-QF/QR (表1)。荧光定量PCR反应体系10 µL, 包括Power 2×SYBR PCR Premixture 5 µL、上下游引物各0.3 µL、cDNA 2 µL, ddH2O补足10 µL。荧光定量PCR反应程序为95℃预变性2 min, 95℃变性15 s, 61℃退火15 s, 72℃延伸20 s, 40个循环。每个反应设置3个生物学重复, 每个生物学重复设置3个技术重复。

1.5 红麻HcWRKY71基因通过花序浸染法转化拟南芥

设计含有H I和I 2个限制性内切酶位点的表达引物PC1301-35S--GFP-F/R (表1), 以阳性克隆载体pMD19-T为模板进行PCR扩增, 并回收PCR扩增产物。利用限制性内切酶H I,I分别对表达载体pC1301-35S-GFP和PCR产物进行双酶切, 利用T4DNA聚合酶将目的基因片段和载体连接。将连接好的重组产物转化到DH5α感受态细胞中, LB培养基培养1 h, 涂含有卡那霉素的平板, 37℃倒置培养18 h, 经菌落PCR测序验证获得重组表达质粒pC1301-35S-- GFP。将重组质粒通过电击法转化到农杆菌EHA105中, 用含有卡那霉素(50 mg L-1)和利福平(50 mg L-1)双抗的LB固体培养基划线培养、筛选, 获得含有功能质粒的农杆菌。挑取单菌落接种到含有卡那霉素(50 mg L-1)和利福平(50 mg L-1)的LB液体培养基, 培养至OD600为1.2。离心收集菌液, 然后加入50 mL悬浮液(1/2 MS培养基, 50 g L-1蔗糖, 0.5 g L-12-吗啉乙磺酸, 用NaOH (1 mol L-1)溶液调pH 5.7, 加表面活性剂Silwet至终浓度为0.02%)悬浮菌液。将培养4周左右刚开花的拟南芥植株上的果荚去掉, 将整个花序浸入悬浮好的菌液10 s, 暗培养24 h, 正常培养1周后再浸染1次。果荚变黄后收获T0代转基因种子。

1.6 转基因拟南芥纯合株系的筛选及幼苗耐盐性鉴定

用次氯酸钠消毒T0代转基因种子6~7 min, 然后用蒸馏水洗涤3~5次。用移液枪将消毒后的种子点种在含有40 mg L-1潮霉素的MS培养基上, 对转基因种子进行抗性筛选。将健壮绿色的转基因幼苗移栽到灭菌的营养土培养, 然后单株收获T1代种子。根据孟德尔遗传定律T1代转基因种子含有杂合种子需要进一步进行筛选, 筛选方法同上, 然后单株收获T2代种子。T2代种子继续筛选, 如果T2代单株收获的种子没有发生分离, 则该单株为纯合植株, 此时筛选单株收获的种子为T3代, 则该单株收获的T3代种子为纯合株系。

将筛选获得纯合的转基因拟南芥种子和野生型拟南芥种子消毒后分别播种于MS培养基上, 待拟南芥幼苗长到四叶一心时移栽到装有灭菌营养土(营养土与蛭石按3︰1的比例混合)的塑料花盆。光照培养室培养5 d后用100 mmol L-1NaCl溶液胁迫处理, 对转基因幼苗耐盐性进行鉴定。

1.7 转基因拟南芥目的基因PCR鉴定及荧光定量分析

将筛选获得纯合的转基因拟南芥种子和野生型拟南芥种子消毒后分别播种于MS培养基上, 待拟南芥幼苗长到四叶一心时移栽到装有灭菌营养土(营养土与蛭石按3︰1的比例混合)的塑料花盆。光照培养室培养7 d后用100 mmol L-1NaCl溶液胁迫处理, 24 h后分别提取对照条件下野生型拟南芥植株和转基因拟南芥植株的总RNA以及盐胁迫下野生型拟南芥植株和转基因拟南芥植株的总RNA, 进行目的基因检测和荧光定量分析。用于目的基因PCR检测和荧光定量分析的引物HcWRKY71-PF/R (表1)。

2 结果与分析

2.1 红麻HcWRKY71基因的克隆和测序

根据红麻转录组unigene序列(编号: CL3883.Contig4), 设计特异性引物, 进行PCR扩增, 获得长度为957 bp的cDNA序列, 扩增产物经1.2%琼脂糖凝胶电泳检测, 结果如图1所示。对PCR产物进行切胶回收, 然后连接到pMD19-T载体上, 送生工生物工程(上海)股份有限公司测序。

M: 2K Plus II; 1: PCR产物。

M: DNA marker 2K Plus II; 1: PCR product.

2.2 红麻HcWRKY71基因的生物信息学分析

PCR产物经sanger测序, 获得长度为957 bp cDNA序列, 经NCBI ORF Finder在线预测,基因的开放读码框为957 bp, 编码1个含有318个氨基酸的蛋白(图2)。NCBI在线预测该蛋白的保守功能域发现, 该蛋白含有1个WRKY功能保守结构域(172~229位) (附图1), 属于第二类WRKY转录因子。利用在线工具预测该蛋白的理论分子量为35.5 kD, 等电点为7.24。蛋白质磷酸化位点预测结果(附图2)显示, 该蛋白质含有32个丝氨酸、7个酪氨酸、18个苏氨酸, 其中位于第189位、194位、212位的丝氨酸和位于第199位、200位、204位、218位和221位的苏氨酸位的磷酸化位点, 位于WRKY的保守功能域内, 推测这8个磷酸化位点可能对HcWRKY71转录因子活性的调控具有重要作用。

HcWRKY71转录因子氨基酸序列相似性通过NCBI BlastP比对发现, 红麻与雷蒙德氏棉(XP_012445471.1)、可可(XP_017978107.1)、榴莲(XP_022720946.1)、哥伦比亚锦葵(XP_021293950.1)、葡萄(XP_002272089.1)、枣(XP_015876832.1)、梨(XP_009353824.1)、苹果(XP_008383508.1)、橄榄(XP_022879638.1)的氨基酸序列相似性分别为84%、80%、79%、79%、60%、60%、57%、57%、53%。DNAMAN的分析结果(图3)表明, 上述物种氨基酸序列的一致性为67.92%。WRKY的功能结构域高度保守, 其他位置氨基酸序列, 由于物种的不同表现出不同的变化。WRKY71转录因子可能在不同的物种中具有相似的生物学功能。

利用MEGA6.0软件分析HcWRKY71转录因子氨基酸序列与具有代表性的雷蒙德氏棉(XP_012445471.1)、可可(XP_017978107.1)、榴莲(XP_022720946.1)、哥伦比亚锦葵(XP_ 021293950.1)、葡萄(XP_002272089.1)、枣(XP_ 015876832.1)、梨(XP_009353824.1)、苹果(XP_ 008383508.1)、橄榄(XP_022879638.1)氨基酸序列的进化关系(图4)发现, 红麻与雷蒙德氏棉的关系最近, 这可能是两者同属于锦葵科的缘故。而红麻与橄榄的关系最远。

红色框为WRKYGQK基序; 黑色框为C2H2基序(CX4CX23HXH)。

The sequence of WRKYGQK motif is highlighted in red box, that of C2H2motif (CX4CX23HXH) in black box.

2.3 盐胁迫下红麻HcWRKY71基因的表达分析

分别提取0、50、100、200 mmol L-1NaCl胁迫下红麻叶片总RNA, 并反转录合成cDNA作为荧光定量PCR反应的模板, 进行定量分析。基因表达量随NaCl浓度的增加而增加。当NaCl浓度为100 mmol L-1和200 mmol L-1时基因表达量与对照表达量相比, 达到极显著差异水平(< 0.01) (图5-A)。

分别提取对照和200 mmol L-1NaCl胁迫下红麻根、茎、叶的总RNA, 反转录合成cDNA进行荧光定量PCR分析。基因在根、茎、叶中均有表达, 其表达量没有显著差异。在NaCl胁迫下, 根、茎、叶中基因表达量均有增加, 其中根部的表达量上升最多, 与对照表达量相比, 具有极显著差异(< 0.01) (图6-A)。可见基因的表达是受NaCl胁迫诱导。

2.4 干旱胁迫下红麻HcWRKY71基因的表达分析

分别提取0、0.5、1.0、2.0、4.0 h干旱胁迫下红麻叶片总RNA, 并反转录合成cDNA作为荧光定量PCR反应的模板, 进行定量分析。基因表达量随干旱胁迫时间的增加先下降再上升然后再下降。干旱胁迫0.5 h时,基因的表达量下调, 其表达量与对照的表达量相比, 达到显著差异水平(< 0.05); 干旱胁迫2 h时,基因的表达量达到最高, 其表达量与对照的表达量相比, 达到极显著差异水平(< 0.01) (图5-B)。

不同颜色代表不同氨基酸残基的保守性。蓝色表示氨基酸完全保守; 粉红色、青色、黄色分别表示氨基酸的保守性为75%以上、50%以上及33%以上; 白色表示氨基酸的保守性不足33%。

The different color indicates that the different levels of the conservation amino acid residues. The blue means amino acids are complete conservation, the pink, cyan, yellow means that the conservation of amino acids is no less than 75%, no less than 50%, no less than 33% and white means no more than 33%.

分别提取对照和干旱胁迫8 h时红麻植株根、茎、叶的总RNA, 反转录合成cDNA进行荧光定量PCR分析。基因在根、茎、叶中均有表达, 其表达量没有显著差异。在干旱胁迫条件下, 根、茎、叶中基因表达量均有增加, 且根中表达量增加最多, 与对照相比差异不显著(图6-B)。这可能是干旱胁迫后期基因表达量出现下降的原因。综上所述,基因的表达受干旱胁迫诱导。

2.5 CdCl2胁迫下红麻HcWRKY71基因的表达分析

分别提取0、50、100、150 μmol L-1CdCl2胁迫下红麻叶片总RNA, 并反转录合成cDNA作为荧光定量PCR反应的模板, 进行定量分析。基因的表达量随着CdCl2浓度的增加而减少, 其表达量与对照的表达量相比, 达到极显著差异水平(< 0.01) (图5-C)。

分别提取对照和CdCl2(100 μmol L-1)胁迫下红麻根、茎、叶的总RNA, 反转录合成cDNA进行荧光定量PCR分析。基因在根、茎、叶中均有表达, 其表达量没有显著差异; 在CdCl2胁迫条件下, 根、茎、叶中基因表达量均有下降, 与对照相比根和叶中的表达量下降达显著水平, 其中根部的表达量下降最多(图6-C)。由此可见,基因的表达受CdCl2胁迫诱导。

A: 盐胁迫; B: 干旱胁迫; C: 镉胁迫。*和*分别表示在0.05和0.01水平上差异显著。误差线为每组处理的标准误差(= 3)。

A: salt stress; B: drought stress; C: cadmium stress. *, ** mean significant difference at the 0.05 and 0.01 probability levels, respectively. Error bars represent the standard error of each treating group (= 3).

A: 盐胁迫; B: 干旱胁迫; C: 镉胁迫。*和*分别表示在0.05和0.01水平上差异显著。误差线为每组处理的标准误差(= 3)。

A: salt stress; B: drought stress; C: cadmium stress. *, ** mean significant difference at the 0.05 and 0.01 probability levels, respectively. Error bars represent the standard error of each treating group (= 3).

2.6 转基因拟南芥纯合株系的获得及幼苗耐盐性鉴定

通过农杆菌浸染花序法将红麻基因转化拟南芥, 收获转基因种子。利用含有40 mg L-1潮霉素的MS培养基对收获的转基因拟南芥种子进行筛选, 结果如图7-A所示。绿色健壮的幼苗为转基因苗, 黄色弱小的幼苗为非转基因苗。转基因种子T0代、T1代、T2代经过潮霉素筛选最终获得纯合的T3代转基因株系。纯合转基因植株的获得为进一步研究基因的功能奠定坚实的基础。野生型和转基因拟南芥幼苗在100 mmol L-1NaCl溶液胁迫下生长10 d后发现, 转基因拟南芥幼苗的耐盐性显著高于野生型拟南芥幼苗的耐盐性。

2.7 转基因拟南芥HcWRKY71基因的PCR鉴定和实时定量PCR分析

在转基因拟南芥植株的泳道能够观察到明显的条带, 而在野生型拟南芥植株的泳道没有看到明显的电泳条带(图7-B)。这表明基因已经整合到转基因拟南芥植株的基因组内。基因在对照野生型拟南芥植株和盐胁迫下拟南芥植株中的表达量没有显著差异。基因在转基因拟南芥植株中的表达量显著高于其在野生型拟南芥植株中的表达量(< 0.05), 在盐胁迫下转基因拟南芥植株中的表达量极显著高于其在野生型拟南芥植株中的表达量(< 0.01) (图7-C)。由此可见,基因在转基因拟南芥中同样受到盐胁迫的诱导, 在转基因拟南芥响应盐胁迫的过程中同样具有重要的作用。

A: 潮霉素鉴定; B: PCR鉴定; C: 实时定量PCR鉴定。*和*分别表示在0.05和0.01水平上差异显著。误差线为每组处理的标准误差(= 3)。

A: identification of transgenicby hygromycin; B: PCR identification of transgenicplant; C: real time quantitative PCR identification. *, ** mean significant difference at the 0.05 and 0.01 probability levels, respectively. Error bars represent the standard error of each treating group (= 3).

3 讨论

WRKY转录因子家族是植物体内最大的转录因子家族之一, 在多种植物中都发现了WRKY转录因子家族成员。本课题组通过对转录组序列CL3883.Contig4进行序列分析发现, CL3883.Contig4中间含有1个完整的阅读框, 即包含具有完整的起始密码子和终止密码子序列。根据完整的阅读框序列设计引物、扩增、测序获得基因完整的cDNA序列。根据生物信息学分析发现,基因的编码产物为红麻HcWRKY71转录因子。

转录因子WRKY71在植物生长发育[14]、逆境胁迫[15]、防御[16]等方面具有重要的调控功能。目前有关转录因子WRKY71的调控功能主要集中在植物的生长发育和生物逆境胁迫响应方面。在拟南芥中, AtWRKY71通过调控基因提高腋芽分生组织的活性, 从而增加拟南芥的分枝数; 通过直接激活开花基因和加速拟南芥开花[17-19]。Qin等[20]研究表明, 在转基因拟南芥中, 通过小麦基因的过表达调控拟南芥植株内IAA-MeIAA的激素水平来调控拟南芥下部叶片的发育。邹克琴等[21]研究表明, 水稻基因的表达受水稻稻曲菌的侵染诱导, 表现出先升高后下降的趋势, 推测其可能在稻曲菌防御过程中发挥重要作用。

在植物响应非生物胁迫过程中WRKY71转录因子同样具有重要的调控作用。在拟南芥中,基因表达受盐胁迫的诱导, 其过程快速而短暂。在转基因的拟南芥超表达株系中, AtWRKY71通过抑制基因和的表达促进拟南芥开花[22]。Guo和Qin[18]的研究表明,基因可能参与多种非生物胁迫的生理响应过程。在水稻中, OsWRKY71通过调控下游基因和的表达来提高转基因水稻植株的耐寒性[23]。本研究通过荧光定量PCR分析发现, 在NaCl溶液胁迫下,基因的表达量随着NaCl溶液浓度的提高而增加, 这与耐盐转录组的测序结果[13]相一致; 在干旱胁迫下,基因的表达量随着时间的延长呈现出先下降在上升然后在下降的规律; 在CdCl2溶液胁迫下,基因的表达量随着CdCl2溶液浓度的增加而降低。在转基因拟南芥植株中,基因的表达同样受到盐胁迫的诱导, 其表达量与野生型拟南芥的表达量相比显著上调。表明,基因同时参与红麻对盐胁迫、干旱胁迫和重金属镉胁迫的响应过程, 且在不同的逆境胁迫下具有不同的表达方式, 推测该基因在红麻非生物逆境胁迫调控网络中具有关键枢纽作用。

在NaCl溶液、干旱和CdCl2溶液不同逆境胁迫下,基因在根中的表达量变化最大(图6), 推测该基因在逆境胁迫响应过程中主要通过调控根部的基因表达来提高红麻自身的抗逆性。

4 结论

从红麻16号中成功克隆出1个基因。该基因编码1个含有318个氨基酸的II类WRKY转录因子。基因是一个组成型且受盐、干旱和重金属镉胁迫诱导表达的基因。通过农杆菌介导的花序浸染法和潮霉素筛选最终获得了转基因拟南芥的纯合植株。苗期耐盐性鉴定发现,基因提高了转基因拟南芥幼苗的耐盐性。这为进一步研究该基因耐盐的调控机制奠定了基础。

[1] Zhang L W, Xu Y, Zhang X T, Ma X K, Zhang L L, Liao Z Y, Zhang Q, Wan X B, Cheng Y, Zhang J S, Li D X, Zhang L M, Xu J T, Tao A F, Lin L H, Fang P P, Chen S, Qi R, Xu X M, Qi J M, Ming R. The genome of kenaf (L.) provides insights into bast fibre and leaf shape biogenesis., 2020, 18: 1796–1809.

[2] Yang A, Dai X, Zhang W H. A R2R3-type MYBgene, is involved in salt, cold, and dehydration tolerance in rice., 2012, 63: 2541–2556.

[3] Duan Y J, Jiang Y Z, Ye S L, Karim A, Ling Z Y, He Y Q, Yang S Q, Luo K M. PtrWRKY73, a salicylic acid-inducible poplar WRKY transcription factor, is involved in disease resistance in., 2015, 34: 831–841.

[4] Wang G, Zhang S, Ma X, Wang Y, Kong F, Meng Q. A stress-associated NAC transcription factor (SlNAC35) from tomato plays a positive role in biotic and abiotic stresses., 2016, 158: 45–64.

[5] 田云, 卢向阳, 彭丽莎, 方俊. 植物WRKY转录因子结构特点及其生物学功能. 遗传, 2006, 28: 1607–1612. Tian Y, Lu X Y, Peng L S, Fang J. The structure and function of plant WRKY transcription factors., 2006, 28: 1607–1612 (in Chinese with English abstract).

[6] 史书婷, 龚小庆, 邹养军. WRKY转录因子与植物逆境响应. 中国生物化学与分子生物学报, 2017, 33: 674–680. Shi S T, Gong X Q, Zou Y J. The role of WRKY transcription factors in biotic and abiotic stress responses in plants., 2017, 33: 674–680 (in Chinese with English abstract).

[7] Wu J, Chen J B, Wang L F, Wang S M. Genome-wide investigation of WRKY transcription factors involved in terminal drought stress response in common bean., 2017, 8: 380.

[8] 王岩岩, 张永兴, 郭葳, 代文君, 周新安, 矫永庆, 沈欣杰. 野生大豆转录因子基因的克隆与抗旱性功能分析. 中国油料作物学报, 2019, 41: 524–530. Wang Y Y, Zhang Y X, Guo W, Dai W J, Zhou X A, Jiao Y Q, Shen X J. Cloning and function oftranscription factors gene response to drought stress., 2019, 41: 524–530 (in Chinese with English abstract).

[9] 柯丹霞, 彭昆鹏, 夏远君, 朱玉莹, 张丹丹. 盐胁迫应答基因的克隆及转基因百脉根的抗盐分析. 草业学报, 2018, 27(8): 95–106. Ke D X, Peng K P, Xia Y J, Zhu Y Y, Zhang D D. Clonging of salt-stressed responsive geneand salt resistance analysis of transgenic., 2018, 27(8): 95–106 (in Chinese with English abstract).

[10] 王玲, 刘峰, 戴明剑, 孙婷婷, 苏炜华, 王春风, 张旭, 毛花英, 苏亚春, 阙友雄. 甘蔗基因的克隆与表达特征性分析. 作物学报, 2018, 44: 1367–1379.Wang L, Liu F, Dai M J, Sun T T, Su W H, Wang C F, Zhang X, Mao H Y, Su Y C, Que Y X. Cloning and expression characteristic analysis ofgene in sugarcane., 2018, 44: 1367–1379 (in Chinese with English abstract).

[11] Zhang Y, Yu H, Yang X, Li Q, Ling J, Wang H, Gu X, Huang S, Jiang W. CsWRKY46, a WRKY transcription factor from cucumber, confers cold resistance in transgenic-plant by regulating a set of cold-stress responsive genes in an ABA dependent manner., 2016, 108: 478–487.

[12] 王影, 邱文敏, 李鹤, 贺雪莲, 刘明英, 韩小娇, 曲同宝, 卓仁英. 东南景天基因对镉胁迫的响应研究. 南京林业大学学报(自然科学版), 2019, 43(3): 59–66.Wang Y, Qiu W M, Li H, He X L, Liu M Y, Han X J, Qu T B, Zhuo R Y. Research on the response ofgene to cadmium stress inhance.(Nat Sci Edn), 2019, 43(3): 59–66 (in Chinese with English abstract).

[13] Li H, Li D F, Chen A G, et al.RNA-seq for comparative transcript profiling of kenaf under salinity stress., 2017, 130: 365–372.

[14] Tan X L, Fan Z Q, Li L L, Wu Y, Kuang J F, Lu W J, Chen J Y. Molecular characterization of a leaf senescence-related transcription factor BrWRKY75 of Chinese flowering cabbage., 2017, 2: 272–278.

[15] Wang C, Deng P, Chen L, Wang X, Ma H, Hu W, Yao N, Feng Y, Chai R, Yang G, He G. A wheat WRKY transcription factor TaWRKY10 confers tolerance to multiple abiotic stresses in transgenic tobacco., 2013, 8: e65120.

[16] 禹阳, 贾赵东, 马佩勇,郭小丁, 谢一芝, 边小峰. WRKY 转录因子在植物抗病反应中的功能研究进展. 分子植物育种, 2018, 16: 7009–7020. Yu Y, Jia Z D, Ma P Y, Guo X D, Xie Y Z, Bian X F. Research progress on the role of WRKY transcription factors in plant defense., 2018, 16: 7009–7020 (in Chinese with English abstract).

[17] Guo D S, Zhang J Z, Wang X L, Han X, Wei B Y, Wang J Q, Li B X, Yu H, Huang Q P, Gu H Y, Qu L J, Qin G J. The WRKY transcription factor WRKY71/EXB1 controls shoot branching by transcriptionally regulationgenes in., 2015, 27: 3112–3127.

[18] Guo D S, Qin G J. EXB1/WRKY71 transcription factor regulates both shoot branching and responses to abiotic stresses., 2016, 11: e1150404.

[19] Yu Y C, Liu Z H, Wang L, Kim S G, Seo P J, Qiao M, Wang N, Li S, Cao X, Park C M, Xiang F. WRKY71 accelerates flowering via the direct activation of FLOWERING LOCUS T and LEAFY in., 2016, 85: 96–106.

[20] Qin Z, Lyu H J, Zhu X L, Chen M, Quan T Y, Wang M C, Xia G M. Ectopic expression of a wheat WRKY transcription factor generesults in hyponastic leaves in., 2013, 8: e63033.

[21] 邹克琴, 汪得凯, 胡东维, 李素芳, 陈振煌. 水稻WRK转录因子OsWRKY71的cDNA分离和表达分析. 浙江农业学报, 2014, 26: 1405–1411. Zou K Q, Wang D K, Hu D W, Li S F, Chen Z H. Isolation and expression analysis oftranscription factors in rice., 2014, 26: 1405–1411 (in Chinese with English abstract).

[22] Yu Y C, Wang L, Chen J C, Liu Z H, Park C M, Xiang F X. WRKY71 acts antagonistically against salt-delayed flowering in., 2018, 59: 414–422.

[23] Kim C Y, Vo K T X, Nguyen C D, Jeong D H, Lee S K, Kumar M, Kim S R, Park S H, Kim J K, Jeon J S. Functional analysis of a cold-responsive rice WRKY gene,., 2016, 10: 13–23.

附图1 红麻HcWRKY71蛋白的保守功能结构域预测

Fig. S1 Conserved domain prediction of HcWRKY71 protein in kenaf

附图2 HcWRKY71转录因子磷酸化位点的预测

Fig. S2 Prediction of phosphorylation sites of transcription factor HcWRKY71

Expression analysis of abiotic stress response genein kenaf and transformation of

LI Hui1,2, LI De-Fang2,*, DENG Yong2, PAN Gen, CHEN An-Guo2, ZHAO Li-Ning2, and TANG Hui-Juan2

1College of Life and Environment Science, Hunan University of Arts and Science, Changde 415000, Hunan, China;2Institute of Bast Fiber Crops, Chinese Academy of Agricultural Sciences, Changsha 410205, Hunan, China

WRKY transcription factor plays an important role in plant responded to abiotic stress. In this study, the unigene sequence (CL3883.Contig4) of the transcriptome in kenaf was used as a reference. Primers were designed for PCR amplification. The full length ofgene was 957 bp obtained by the sanger sequencing.gene had an open reading frame length of 957 bp, and encoded a protein containing 318 amino acids, with a conserved functional domain of WRKY, which belonged to WRKY transcription factor II. Under the salt stress, its relative expression level increased with the increase of NaCl concentration; under the drought stress, with the extension of drought time, the relative expression of thegene first decreased, then increased and finally decreased; under the stress of heavy metal cadmium, its expression decreased with the increase of CdCl2concentration, indicating that the expression of the gene was induced by salt, drought and heavy metal cadmium stress. The gene was transformed intoby-mediated inflorescence impregnation. It was found that thegene improved the salt tolerance of transgenicseedlings. This laid a solid foundation for further study of stress tolerance mechanism ofgene.

kenaf; salt tolerance; drought; heavy metal cadmium;

10.3724/SP.J.1006.2021.04201

本研究由国家现代农业产业技术体系建设专项(CARS-19-E07), 中国农业科学院科技创新工程一年生麻类育种项目(ASTIP-IBFC03), 湖南省教育厅项目(18C0737)和湖南文理学院博士科研启动项目(17BSQD13)资助。

This study was supported by the China Agriculture Research System (CARS-19-E07), the Agricultural Science and Technology Innovation Program at the Chinese Academy of Agricultural Science (ASTIP-IBFC03), the Hunan Education Department Project (18C0737), and the Doctoral Research Start-up Project of Hunan University of Arts and Sciences (17BSQD13).

李德芳, E-mail:18973612200@163.com

E-mail: guangjunmuzi@126.com

2019-12-19;

2020-10-14;

2020-09-08.

URL: https://kns.cnki.net/kcms/detail/11.1809.S.20200905.1226.002.html