唐建民 齐力旺 贺小婷 蒋雨欣 李哲馨

摘 要：生姜（zingiber officinale roscoe）是药食兼用的高附加值特色经济作物，姜辣素是生姜的核心药用成分，6-姜酚则是姜辣素的主要成分。以川渝特色生姜品种竹根姜根茎为研究材料，克隆zowrky1基因并分析其在竹根姜不同发育阶段的表达及其与6-姜酚含量增量的相关性，探讨不同浓度nacl胁迫下zowrky1与6-姜酚合成酶基因zopal、zoc4h和zo4cl的表达及其相关性。结果表明：克隆的zowrky1基因cdna全长1 026 bp， 登录号为mt414710;氨基酸序列分析显示zowrky1属于wrky家族的第二类成员;zowrky1的表达水平在竹根姜发育的第二阶段即播种后1个月左右高度上调，且与6-姜酚含量增量存在显著的相关性;在25 g/l nacl胁迫下，zowrky1与6-姜酚合成酶基因zopal、zoc4h和zo4cl高度上调表达，且zowrky1与zoc4h和zo4cl的表达存在显著相关性，推测zowrky1基因与zoc4h及zo4cl基因存在调控关系并进一步影响6-姜酚的生物合成。

关键词：生姜;wrky转录因子;6-姜酚合成;相关性

中图分类号：q786 文献标识码：a文章编号：1006-060x（2020）10-0001-05

abstract： ginger （zingiber officinale rosc.） is a special economic crop with high added value for both medicine and food. pungent prinsiple is the core medicinal component of ginger， and 6-gingerol is the main component of pungent prinsiple. zowrky1 gene was cloned from rhizomes of zhugen ginger in sichuan and chongqing. the of zowrky1 gene in different developmental stages and its correlation with 6-gingerol content increment were analyzed. the and correlation of zowrky1 and 6-gingerol synthase genes zopal， zoc4h and zo4cl under different concentrations of nacl stress were investigated. the results showed that the full length cdna of zowrky1 gene was 1026 bp and the accession number was mt414710; amino acid sequence analysis showed that zowrky1 belonged to the second class member of wrky family; the level of zowrky1 was highly up-regulated in the second stage of zhugen ginger development about one month after sowing， and there was a significant correlation between zowrky1 and 6-gingerol content increment; under 25 g/l nacl stress， the of zowrky1 and 6-gingerol synthetase genes zopal， zoc4h and zo4cl?were highly up-regulated， and the of zowrky1was significantly correlated with the of zoc4h and zo4cl. it is speculated?that there is a regulatory relationship between zowrky1and zoc4h and zo4cl genes， and further affect the biosynthesis of 6-gingerol.

key words： ginger （zingiber officinale rosc.）; wrky tran ion factor; 6-gingerol synthesis; correlation

生姜（zingiber officinale rosc.）是姜科（zingiber-aceae）姜属（zingiber）多年生宿根草本单子叶植物，也是药食两用的经济作物。姜辣素（pungent prinsiple）是生姜中的辛辣味物质和药用核心成分[1-4]。姜酚类（gingerols）物质是姜辣素的主要成分，尤以6-姜酚含量最高，可达姜辣素含量的75% 以上。含有苯丙烷结构的物质都是通过苯丙烷代谢途径直接或者间接生成的。目前在已知的包括6-姜酚在内的含有苯丙烷結构的物质生物合成通路中，苯丙氨酸氨裂解酶（phenylalanine ammonia lyase，pal）、肉桂酸-4-羟化酶（cinnamate 4-hydroxylase，c4h）、4-香豆酸辅酶a连接酶（4-coumarate：coa ligase，4cl）等酶在6-姜酚的合成过程中扮演着重要角色[5-6]。大量试验数据表明这些酶基因控制着苯丙烷类代谢物的合成。在阿月浑子和银杏树中，pal基因的表达量与黄酮含量均表现出正相关关系[7-8]。通过胁迫处理降低菊花c4h基因表达量，也会导致黄酮类物质量的减少[9]。桉树中下调c4h的表达会导致木质素含量减少18.9%[10]。在短莛飞蓬和凹唇姜中，4cl基因的表达水平与黄酮类物质的积累呈现显著的正相关性[11-12]。在烟草或拟南芥中抑制pal、c4h、4cl的表达，均会导致木质素含量的降低，并造成木质素构成的改变[13]。将c4h反义基因转入紫花苜蓿，发现与野生型相比，基因沉默植株木质素含量均有不同程度的降低[14]。

植物中诱导型基因的表达主要受转录水平的调控，而在转录水平上起主要作用的是各类转录因子。wrky转录因子在植物木质素、黄酮类等苯丙烷类物质合成过程中起着重要的调控作用。欧洲葡萄vvwrky2基因可激活木质素生物合成途径中vvc4h的表达，进而调控木质素的合成与积累，这一功能在烟草中获得验证[15]。在毛白杨ptrwrky25超量表达的转基因株系中木质素含量相较于野生型减少约1倍，表明ptrwrky25可能主要在木质素合成途径中发挥作用[16]。其他植物如苜蓿、云杉等也有wrky转录因子参与调控木质素合成的相关报道[17-18]。拟南芥中，atwrky23能通过上调编码黄酮生物合成途径相关酶基因的表达来刺激黄酮类生物合成，且该调控机制与生长素信号途径密切相关[19]。葡萄中筛选的vvwrky26转录因子在花青素合成中具有重要的调控作用[20]。过表达苜蓿mdwrky11基因可促进f3h、fls、dfr、ans和ufgt的表达，进而导致苹果愈伤组织中黄酮和花青素的积累量增加[21]。

植物在与环境相互作用时获得了一系列的防卫反应机制，在应对外界各种生物和非生物的胁迫时，植物能通过体内抗性基因的表达，合成并积累次生代谢产物来抵御逆境胁迫以减轻伤害[22]。根据前期已经获得的竹根姜肉质根、地下根茎、地上茎及叶片的rna-seq測序结果[23]，筛选到23个wrky基因（zowrky），按顺序编号为1～23，其中有5个zowrky（命名为zowrky1～zowrky5）在地下根茎中高度表达。其中zowrky1 基因被盐胁迫诱导后表达量变化较大。该研究重点对zowrky1进行克隆与生物信息学分析，探索其在竹根姜发育不同阶段以及不同浓度盐胁迫下的表达，并结合6-姜酚含量及其生物合成关键酶基因的表达量分析，阐释 zowrky1与6-姜酚合成可能存在的调控关系，以期为zowrky1在生姜6-姜酚生物合成调控中的功能研究奠定基础，并为实现植物化学物质合成与积累的人工调控提供依据。

1 材料与方法

1.1 试验材料及样品采集

选取川渝特色生姜品种竹根姜[24] “渝姜1号”作为研究材料，每年4月开始种植母姜（催芽1个月），10月采收。 根据不同的研究内容分别采样。（1）不同发育阶段竹根姜样品采集：每月15日采集第一分支子姜根茎材料，4月15日—9月15日共6个阶段，记为stage 1～stage 6。（2）nacl胁迫下竹根姜样品采集：对播种一个月的竹根姜植株叶片喷施不同浓度nacl溶液（浓度分别为0、25、50、75、100和125 g/l）[25]，2 d喷施一次，第10天收获地下根茎。上述样品采集后均置于液氮中快速冷冻，于-80℃冰箱保存备用。

1.2 总rna提取与cdna合成

样品总rna提取采用rnaiso plus total rna 提取试剂盒（takara，大连），琼脂糖凝胶电泳检测结果显示各阶段样品28s和 18s rna 清晰可见且无降解，利用trans all-in-one first-strand cdna synthesis supermix for qpcr（one-step gdna removal）试剂盒（transgen，北京）反转录合成cdna。

1.3 基因克隆与序列分析

根据转录组测序结果设计zowrky1基因序列的特异引物（表 1），以不同发育阶段竹根姜混合cdna为模板进行 pcr 扩增，反应体系（20 μl）：0×buffer?10 μl、1.5 mmol/ l mg2+ 1.5 μl、0.1 mmol/ l dntps 4 μl、0.4 μmol/ l f、r引物各0.5 μl、0.5 u rtaq dna 聚合酶0.2 μl、200 ng/μl cdna 模板1 μl、ddh2o 2.3 μl。pcr 扩增反应程序：95℃ 5 min;95℃ 10 s，56℃ 30 s，72 ℃ 1 min，35个循环;72℃ 7 min。pcr扩增产物于1.0%的琼脂糖凝胶电泳检测。扩增产物条带单一，切胶回收后连接至克隆载体，克隆试剂盒选用peasy-blunt simple cloning kit，转化大肠杆菌top10，筛选阳性克隆，对菌液进行pcr检测和测序验证。采用dnaman软件分析zowrky1序列组分，并采用在线工具http：//smart.embl-heidelberg.de/预测其蛋白结构。

1.4 qrt-pcr分析

将不同发育阶段及nacl胁迫处理的竹根姜根茎组织cdna分别稀释10倍作为反应模板，采用tb greentm premix ex taqtm （tli rnaseh plus）试剂盒（takara，大连）在qtower2.0实时荧光定量pcr系统（analytik jena，德国）进行qrt-pcr检测与分析。引物序列如表1所示。反应体系（20 μl）：cdna 2 μl、tb green Ⅱ 10 μl、引物各0.8 μl、ddh2o 6.4 μl。反应条件：95℃ 30 s;95℃ 5 s，55℃ 15 s，72℃ 10 s，40个循环。

1.5 相关性分析

前期采用lc-ms技术测得不同发育阶段竹根姜中6-姜酚的含量。6-姜酚各阶段增量为各阶段6-姜酚总量减去其前一阶段6-姜酚总量。采用ibm spss statistics18.0统计软件进行相关性分析。

2 结果与分析

2.1 生姜zowrky1基因的克隆与序列分析

从竹根姜转录组测序文库中获得zowrky1 基因的 cdna 序列，在ncbi中进行blastx比对及 orf 分析，显示该序列具有完整的orf，推测其为全长基因。设计其特异性引物，以不同发育阶段竹根姜根茎组织混合cdna 为模板，进行 pcr 扩增，产物与预期大小相符（图1a）。将 pcr 产物连接到pmd19-t 载体上， 测序后获得竹根姜zowrky1基因的cdna序列，全长1 026 bp， 利用dnaman软件对zowrky1基因序列进行分析，结果显示其编码的蛋白质由341个氨基酸组成。zowrky1基因 cdna 的序列信息已提交到 ncbi genbank，登录号为mt414710。经ncbi比对，该基因与姜目芭蕉科小果野蕉mdwrky21（xm_018829735）的同源性高达92%。氨基酸序列分析结果显示，zowrky1含有1个 wrky保守结构域以及2个lcr（low complexity region）区域，并包含有c2h2 型锌指结构（图1b和 c），说明zowrky1属于wrky家族的第二类成员。

2.2 zowrky1在生姜不同发育阶段的表达

采用qrt-pcr技术检测了竹根姜不同发育阶段zowrky1基因的表达水平，结果如图2所示，zowrky1基因的表达在竹根姜发育过程中变化明显，竹根姜早期发育阶段，即催芽播种后1～2个月（对应stage 2—3）时表达量骤然上升，较播种时增加了5倍之多;随后几个月的表达水平则大幅下降，与刚播种时差异不明显。

2.3 zowrky1表达量与6-姜酚含量增量的相关性分析

经测定，stage 1～stage 6竹根姜样品中6-姜酚含量分别为182.37、498.47、528.56、486.70、478.14、489.34 μg/g。采用ibm spss statistics18.0统计软件对zowrky1基因的表达量与6-姜酚含量增量进行了相关性分析，结果显示，zowrky1基因的表达与6-姜酚含量相关系数为0.831 （p=0.040