谷子WRKY转录因子基因家族分析_邢国芳
（Ｊ．ＳＮＳＥＨＡＮＸＩＡＧＲＩＣ，ＵＮＩＶ．ａｔｕｒａｌｃｉｅｎｃｅｄｉｔｉｏｎ
）
（）：自然科学版）学报（０１６，３６１２８３７２
３４８３０
谷子ＷＲＫＹ转录因子基因家族分析
２
，张莉１，冯万军１，谢圣杰１，贾亚涛１，苑乂川１，陈小雨１邢国芳１，
（）山西农业大学农学院，杂粮种质资源发掘与遗传改良山西省重点实验室，山西太谷０山西太谷０１．３０８０１；２．３０８０１
［谷子是抗旱研究的模式植物，谷子基因组测序的完成，为分离鉴定谷子抗旱基因提供了便利，摘 要：目的］ＫＹＷＲ基因参与植物非生物胁迫应答，在植物抗旱方面也发挥重要作用。为分离鉴定谷子ＷＲ从ＫＹ类抗旱转录因子基因，而为谷子抗旱机理研究和抗旱育种提供参考。［方法］采用生物信息学方法，对谷子转录因子ＳｉＷＲＫＹ基因家族进行了分离鉴定和表达分析。［结果］我们发现，谷子基因组中存在１分布在谷子的９条染色体上，根据０３个ＷＲＫＹ基因，在染色体上的位置命名为Ｓ结论］表明其ＳｉＷＲＫＹ１－ｉＷＲＫＹ１０３。［８０％的谷子ＷＲＫＹ基因在干旱胁迫后上调表达，参与谷子抗旱机制，为植物抗旱研究提供了候选基因。关键词：谷子；干旱胁迫；基因表达ＷＲＫＹ转录因子；
（）中图分类号：Ｓ５１５；Ｑ９４ 文献标识码：６７１１５１２０１６１２８３７９Ａ 文章编号：１８００－－－
DOI:10.13842/j.cnki.issn1671-8151.2016.12.014
ＡＷＲｎａｌｓｉｓｏｆＫＹｔｒａｎｓｃｒｉｔｉｏｎｅｎｅｆａｍｉｌｉｎｆｏｘｔａｉｌｍｉｌｌｅｔ ｙｐｇｙ
１，２１１１１１
，Ｘ，Ｊ，Ｙ，ＣＸｉｎＧｕｏｆａｎｈａｎＬｉｅｎＷａｎｕｎｉｅＳｈｅｎｉｅｉａＹａｔａｏｕａｎＹｉｃｈｕａｎｈｅｎ ，Ｆ ｇｇ，Ｚｇｇｊｇｊ
１
Ｘｉａｏｕｙ
（．Ｃｏｌｌｅｅｉｃｕｌｔｕｒｅ，ｈａｎｘｉｉｃｕｌｔｕｒａｌｎｉｖｅｒｓｉｔａｉｕ３０８０１，ｈｉｎａ，．Ｓｈａｎｘｉ ｅａｂｏｒａｔｏｒ１ｏＡｒＳＡｒＵＴ０Ｃ２ＫＬｏ ｇｙ，ｇｙ ｙｆ ｇｇｆ
ＲａＧＩｏＭＣＴ０Ｃｅｎｅｔｉｃｅｓｏｕｒｃｅｓｎｄ ｅｎｅｔｉｃｍｒｏｖｅｍｅｎｔｉｎｏｒｒｏｓ，ａｉｕ３０８０１，ｈｉｎａ）Ｇ ｆ ｐｐｇ
：［］ＡｂｓｔｒａｃｔＯｂｅｃｔｉｖｅＷＲＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｉｓｏｎｅｏｆｔｈｅｂｉｅｓｔｅｎｅｆａｍｉｌｉｅｓｔｈａｔａｒｅｉｍｏｒｔａｎｔｒｅｕｌａｔｏｒｓｏｆｍａ －ｊｐｇｇｇｐｇ，，，，ｏｒｌａｎｔｒｏｃｅｓｓｅｓｉｎｃｌｕｄｉｎｒｏｗｔｈｄｅｖｅｌｏｍｅｎｔａｎｄｒｅｓｏｎｓｅｓｔｏｓｔｒｅｓｓ．Ｈｏｗｅｖｅｒｔｈｅｒｅｉｓｌｉｍｉｔｅｄｉｎｆｏｒｍａｔｉｏｎａ－ ｊｐｐｇｇｐｐ ，）ｏｕｔｔｈｅｎｕｍｂｅｒｓｓｔｒｕｃｔｕｒｅｓａｎｄｆｕｎｃｔｉｏｎｓｏｆｔｈｏｓｅＷＲＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｓｉｎｆｏｘｔａｉｌｍｉｌｌｅｔ（ｅｔａｒｉａｔａｌｉｃａ（ｂＳｉＬ． ｐ）［］，Ｂｅａｕｖ．．ＭｅｔｈｏｄｓＩｎｔｈｉｓａｅｒｔｈｅＷＲＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｅｎｅｆａｍｉｌｉｎｆｏｘｔａｉｌｍｉｌｌｅｔｗｅｒｅｓｅａｒａｔｅｄａｎｄｉＰ． －ｐｐｐｐｇｙ ，，［］ｄｂｉｏｉｎｆｏｒｍａｔｉｃｓａｎｄｔｈｅｅｘｒｅｓｓｉｏｎｗｅｒｅａｎａｌｚｅｄｂＲＴ－Ｐ１０３ＷＲＫＹｔｒａｎｅｎｔｉｆｉｅｄｂＣＲ．ＲｅｓｕｌｔｓＩｎｏｕｒｓｔｕｄ－ ｐｙｙｑｙｙ ｃｒｉｔｉｏｎｆａｃｔｏｒｎｃｏｄｉｎｅｎｅｓｗｅｒｅｉｄｅｎｔｉｆｉｅｄｉｎｔｈｅｔａｌｉｃａｎｏｍｅ．ＴｈｅｅｎｅｓｗｅｒｅｎａｍｅｄｉＷＲＫＹ１ｏｉｓｅＳ．ｉｅＳｔＳ － －ｐｇｇｇｇ
ＷＲＫＹａｆｏｕｒｄｉｆｆｅｒｅｎｔｔｉｓｓｕｅｓｗａｓａｎａｌｚｅｄｕｔｏｔｈｅｉｒｏｒｄｅｒｏｎｔｈｅｃｈｒｏｍｏｓｏｍｅｓ．Ｔｈｅｅｘｒｅｓｓｉｏｎａｍｏｎ１０３ｃｃｏｒｄｉｎ － ｐｇｙｇ ａｖａｉｌａｂｌｅＲＮＡｓｅｕｅｎｃｉｎｄａｔａ．ＴｈｅｅｘｒｅｓｓｉｏｎｏｆｉｎｕｂｌｉｃｌｉＷＲＫＹｇｎｅｓｉｎｒｅｓｏｎｓｅｔｏｄｅｈｄｒａｔｉｏｎｗｅｒｅａｎａｅｓＳ －ｐｑｇｇｐｙｐｙ ［］，，ｅｌｚＳｒｄｒｏｕｈｔｓｔｒｅｓｓｉｔｅｄｂＲＴ－ＰＣＲ．ＣｏｎｃｌｕｓｉｏｎＩｎｃｏｎｃｌｕｓｉｏｎｍｏｓｔｉＷＲＫＹｇｎｅｓ（０％）ｗｅｒｅｕｅｕｌａｔｅｄｂ－ ＞８ ｙｙｑｇｙｇｐ ，ｉｍｏｒｔａｎｔｒｏｌｅｓｆｏｒｄｒｏｕｈｔｖｅｒｍａｉｎｄｉｃａｔｅｄｔｈａｔｔｈｅｓｅｅｎｅｓｌａｅｓｉｓｔａｎｔｉｎｍｉｌｌｅｔａｎｄｏｕｒｃｏｎｃｌｕｓｉｏｎｗｉｌｌｒｏｖｉｄｅｒ －ｐｇｙｇｐｙｙｐ
’ａｓｅｔｏｆｃａｎｄｉｄａｔｅｓｅｎｅｓｆｏｒｔｈｅｒｅｓｅａｒｃｈｏｆｌａｎｔｓｒｅｓｏｎｓｅｔｏｄｒｏｕｈｔｓｔｒｅｓｓ． ｇｐｐｇ：，，ＷＲ，ＫｅｗｏｒｄｓＦｏｘｔａｉｌｍｉｌｌｅｔＤｒｏｕｈｔｓｔｒｅｓｓＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒＧｅｎｅｅｘｒｅｓｓｉｏｎ ｇｐｐｙ
人们在对谷子抗旱机制的相关研究中 近年来，
取得了一定进展。特别是谷子全基因组序列测序的完成，为谷子基因组中的转录因子基因家族的全面生物信息学分析和鉴定提供了良好的契机。通过生物信息学和基因表达的方法筛选出不少响应干旱胁迫的基因。对谷子抗旱的机制做了很多补
充。然而，谷子的抗旱性是一个极其复杂多基因调
１］
。到目前为止，控网络［对该调控网络的了解可能
仅仅是冰山一角，还有许多方面的研究并未开展，例如近年来研究比较多的与抗／耐旱有关的关键基因ＷＲＫＹ转录因子的生物学功能和抗旱机制的
研究。
收稿日期：修回日期：００００１６９９０１６９－３０２２－－－
），（），，，作者简介：女汉山西夏县人邢国芳（副教授博士，研究方向：作物功能基因组学９８１１－
；））山西省高等学校创新基金（基金项目：国家自然科学基金（２０１５１４５３１５０１３２３
８８３
山西农业大学学报（自然科学版）
０１６２
ＷＲＫＹ转录因子是近年来在植物中发现的一类重要转录因子，因其Ｎ－末端含有高度保守的ＫＹＧＱＫ氨基酸序列而得名，ＫＹ域能特ＷＲＷＲ异地与靶基因启动子Ｗ－从而调控ｂｏｘ序列结合，
２］
。根据ＷＲ靶基因的表达［ＫＹ域的数量及锌指
结构的组成，ＫＹ蛋白家族可以分为３个亚家ＷＲ［３］
。从２族（组）０多年前首次报道ＷＲＫＹ以来，
物学重复。
．２ 保守域预测１
基因组数据和ｃ谷子（Ｓｉｅｔａｒｉａｔａｌｉｃａ）ＤＮＡ 数据均下载于谷子最新数据库版本Ｐｈｔｏｚｏｍｅｙ
（／／／。利ｈｔｏｚｏｍｅ．ｎｅｔｓｅａｒｃｈ．ｈｈｔｔｗｗｗ．ｐｙｐｐ）ｐ：／／／）用ＳｈＭＡＲＴ（ｈｔｔｍａｒｔ．ｅｍｂｌｅｉｄｅｌｂｅｒ．ｄｅｓ－ｐ：ｇ
在线工具分析ＷＲＫＹ蛋白结构域，利用ＯＧ２．０绘制９４个ＷＲＫＹ家族蛋白的功能结Ｄ
构域示意图。
．３ 进化树与基因结构分析１
应用ＭＥＧＡ５．０对９４个谷子ＷＲＫＹ基因和 ２个拟南芥的ＷＲＫＹ转录因子基因进行多序列７
比对分析，将分析结果用邻接法生成系统进化树，ｏｏｔｓｔｒａ００。Ｂ０ ｐ值设置为１１．４ 染色体定位和基因组复制的预测
通过ＭａＤｒａｗ程序绘制出每个基因在染色ｐ体上的位置，并根据其在染色体上的顺序将其命名。串联重复被认为是位于同一条染色体上３０ｋｂ以内的以基因簇形式存在的染色体重复，采用将ｌａｓｔＰ对谷子染色体重复片段进行同源搜索，Ｂ
最初的评价值Ｅ＜１ｅ－０５的５条序列作为靶标序）作列，采用ＭＣＳｃａｎｖ０．８和相似性值（ｅ－５Ｅ＜１
１６，１７］
。为共线性关系评价的重要依据［
ＫＹ转录因子已被发现在植物的防御应答反应ＷＲ
和生长发育过程中起着至关重要的作用。如拟南
［］
水稻的芥中的ＡｔＷＲＫＹ７５与侧根的生长有关４；［］
ＯｓＷＲＫＹ３１也与侧根的形成和延伸有关５。Ａｔ－
［］
而ＡＷＲＫＹ６可以引发衰老的进程６，ｔＷＲＫＹ７０７］
。然而目前对于ＷＲ可以负调控衰老反应［ＫＹ
转录因子报道最多的是在植物应对外界非生物胁
８，９］
。迄今为止，迫时所发挥的作用［已报道的与非
生物胁迫相关的ＷＲＫＹ转录因子主要有１３个。
其中ＯＧｓＷＲＫＹ４５、ｓＷＲＫＹ２０和ＡｔＷＲＫＹ５７在植物应对干旱胁迫过程中发挥重要作用。
谷子具有耐旱性强、水分利用率高、蒸腾系数低等生物学特性，是发掘抗旱基因的优异种质。谷比水稻的子基因组小（谷子单倍体只有４７０Ｍｂ，，所以谷子远远小于玉米、珍珠粟等）３０Ｍｂ略大，４
被当作黍亚科的模式植物来研究。自从第一个中克ＩｏｂＫＹ基因ＳＰＦ１从白薯（ｍｏｅａａｔａｔａｓ）ＷＲ ｐ
１０］
，许多植物的ＷＲ隆出后［ＫＹ基因相继被克隆和
［１１］
、水稻鉴定，包括拟南芥（Ａｔｒａｂｉｄｏｓｉｓｈａｌｉａｎａ） ｐ［［２］３］１１
、、（大豆（玉米ＯｓＧｍｒｚａａｔｉｖａ）ｌｃｉｎｅ ａｘ） ｙｙ
［１４］
（二穗短柄草（ＺｍＢｄｅａ ａｒａｃｈｏｄｉｕｍ ｉｓ－ｙ）、ｙｐ
］［５１
等２有关谷子ｔａｃｈｏｎ）０多种高等植物。然而，ｙ
１．５ 电子表达谱分析
来自谷子穗子、茎秆、叶片和根系四种组织的ＮＡｓｅｕｅｎｃｅ数据全部来自于欧洲核酸数据库Ｒ ｑ
／／／（ｔｔｅｂｉ．ａｃ．ｕｋｅｎａＳＲＸ１２８２２６（ｓｉｗｗｗ．ｈ－ ｐ：ｐ），Ｓ，Ｓ，ｃａＲＸ１２８２２５Ｚ（ｓｔｅｍ）ＲＸ１２８２２４（ｌｅａｆ）
（））。所有比对上的来自于四种组ＲＸ１２８２２３ｒｏｏｔＳ
织的Ｉｌｌｕｍｉｎａｒｅａｄｓ序列均采用Ｂｏｗｔｉｅ２软件定位
１８］
，并采用Ｒ到谷子的基因序列上［ＰＫＭ（ｒｅａｄｓｅｒ ｐ）方法进行均一化处理，基于每ｋｉｌｏｂａｓｅｅｒｍｉｌｌｉｏｎ ｐ
采用个ＳｉＷＲＫＹ基因在特定组织中的ＲＰＫＭ值，软ＩＧＲ ＭｕｌｔｉＥｘｅｒｉｍｅｎｔＶｉｅｗｅｒ（ＭｅＶ４．９．０）Ｔ ｐ
件获得了ＷＲＫＹ基因的组织特异表达图。从谷／／子数据库Ｐｈｔｏｚｏｍｅ（ｈｔｔｈｔｏｚｏｍｅ．ｗｗｗ．ｙｐ：ｐｙ
）／中选取谷子ＷＲｎｅｔｓｅａｒｃｈ．ｈＫＹ基因起始密ｐｐ
利用植物码子上游２００ｂ０ ｐ序列作为启动子区，／／顺式作用元件数据库ＰＬＡＣＥ２６．０（ｈｔｔｗｗｗ．ｐ：／分析基因的顺式作用ｄｎａ．ａｆｆｒｃ．ｏ．ＰＬＡＣＥ）ｇｊｐ
元件。
ｔ．６ ＲＮＡ提取和ＲｅａｌｉｍｅＰＣＲ分析１－
对生长３周的谷子幼苗分别采用２０％ＰＥＧ模
数量、结构和功能等还未ＫＹ转录因子的种类、ＷＲ
见相关报道。因此，我们拟采取生物信息学方法，对对谷子转录因子ＳｉＷＲＫＹ基因家族进行分析，其在干旱胁迫前后的表达进行分析。
１ 材料和方法
１．１ 实验材料
）由山西省农业科学院品质谷子种子（晋谷２１资源研究所温其汾研究员提供。谷子幼苗在温度／相对湿度６光照周期１２℃、５％、６ｈ８ｈ的条件下２
培养。生长５周的幼苗采用２０％ＰＥＧ６０００和
（·盐）处理，在０ｈ，ＮａＣｌ１ｈ，６ｈ，２５０ｍｍｏｌＬ－１
将根和叶片分开收集，于１２ｈ和２４ｈ后分别取样，
液氮中速冻后，置于－８设置３次生０℃冰箱备用，
（）３６１２邢国芳等：谷子ＷＲＫＹ转录因子基因家族分析
９８３
拟干旱处理，分２ｈ、４ｈ、６ｈ、１２ｈ和２４ｈ后，０ｈ、北京）别取样，利用Ｔ提ＲＩｚｏｌ试剂盒（ＩＡＮＧＥＮ，Ｔ采用Ｒ取总ＲｆＮＡ，Ｎａｓｅｒｅｅ的ＤＮａｓｅＩ（ＲＱ１，－
）处理，并采用Ｓ祛除ＤＰｒｏｍｅａＮＡ，ｕｅｒＳｃｒｉｔＩＩＩ ｇｐｐ）和Ｏ反转录酶（ｌｉｏ（ｄＴ）１８（Ｐｒｏｍｅａｎｖｉｔｒｏｅｎ）Ｉｇｇｇ
按照说明书操作，以Ｓ反转录生成ｃＹＢＲＤＮＡ，ｔＧｒｅｅｎ染料法进行荧光定量ＲｅａｌｉｍｅＰＣＲ。所采－ 用实时荧光定量ＰＣＲ仪为ＡＢＩ公司的Ｐｒｉｓｍ７０００。反应体系包括：ａＰＣＲ ＭａｓｔｅｒＭｉｘ２×Ｔ ｑ （含荧光染料）０μｍｏｌＬ－１引物各０．５μＬ、９μＬ、１ ９μＬ和ｃＤＮＡ模板１μｄｄＨ２ＯＬ。反应条件为
；，，，５℃５ｍｉｎ９５℃１５ｓ５６℃２０ｓ７２℃３０ｓ９
，６０℃３０ｓ３５个循环。ＰＣＲ引物采用ＧｅｎＳｃｒｉｔｐ
／／ｔＲｅａｌｉｍｅＰＣＲＰｒｉｍｅｒＤｅｓｉｎ在线软件（ｔｔｈ－ ｇｐ：／／／）设计获ｂｅｎｓｃｒｉｔ．ｃｏｍｓｓｌｉｎａｒｉｍｅｒｗｗｗ．－ｇｐｐｐｐ
）为内标，引物序得。以谷子Ａｃｔｉｎ基因（ｉ００１８７３Ｓ列为ＡＦ：５ｃｔｉｎ′ＧＣＡＡＡＣＡＧＧＧＡ ＧＡＡ－－Ｇ
和Ａ３５ＧＡＴＧＡ′ｃｔｉｎ′ＧＴＴＧＴＣＧＧＴ－－Ｒ：－ＧＡＧ
ｔ
。用２－ΔΔＣ法计算各基因表３ＡＡＧＧＴＣＡＣＧ′－
达量。
２ 结果分析
２．１ 谷子ＳｉＷＲＫＹ基因的克隆
为了得到所有的谷子ＷＲＫＹ基因的潜在序应用Ｂ列，以ＷＲＫＹ基因序列作为靶标，ＬＡＳＴ软件对谷子全基因组进行了扫描，共鉴定出１０３个ＳｉＷＲＫＹ序列。应用在线分析软件ＳＭＡＲＴ（ＳｉｍｌｅＭｏｄｕｌａｒＡｒｃｈｉｔｅｃｔｕｒｅＲｅｓｅａｒｃｈＴｏｏｌ） ｐ
／（／／）对每一个候ｈｍａｒｔ．ｅｍｂｌｅｉｄｅｉｂｅｒ．ｄｅｈｔｔｓ－ｇｐ：选基因进行了鉴定。其中４个基因（Ｓｉ０３２７４８ｍ；
Ｓｉ００３８７６ｍ；Ｓｉ０１３３２６ｍ和Ｓｉ０２７９８２ｍ）缺少
另５个基因（ＳＫＹ基因的保守域，ｉ００２９５７ｍ；ＷＲ
仅ＳＳＳｉ０１５４３３ｍ；ｉ０１０８３２ｍ；ｉ０２７８７６ｍ和Ｓｉ０３９５１９ｍ）
有部分Ｃ其－末端的锌指结构域而均被排除在外，余的９４个转录本根据其在染色体上的位置而被分。Ｓ图１）别命名为ＳＳｉＷＲＫＹ－ｉＷＲＫＹ９４（ｉＷＲＫＹ蛋白编码平均大小为３５８个氨基酸的残基。但其编码氨基酸数目大小仍然存在很大差异，其中最大的Ｓ而最小的ｉ０２８７１０ｍ预测为１９１个氨基酸，２ ｉ０３８９５５ｍ则仅有９４个氨基酸
。Ｓ
图１ １０３个ＳｉＷＲＫＹ基因在９条谷子染色体上的分布
ＦｅＳｉ．１ Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆ１０３ｉＷＲＫＹｇｎｅｓｏｎｔｏｎｉｎｅｆｏｘｔａｉｌｍｉｌｌｅｔｃｈｒｏｍｏｓｏｍｅｓ ｇ
串联重复基因用红色方框标示。染色体物理距离用Ｍｂ来表示 注：
：ｏｔｅＴａｎｄｅｍｄｕｌｉｃａｔｅｄｅｎｅｓｏｎａａｒｔｉｃｕｌａｒｃｈｒｏｍｏｓｏｍｅａｒｅｉｎｄｉｃａｔｅｄｉｎｔｈｅｒｅｄｂｏｘ．ＣｈｒｏｍｏｓｏｍａｌｄｉｓｔａｎｃｅｓＮ ｐｇｐａｒｅｉｖｅｎｉｎＭｂ ｇ
２．２ ＳｉＷＲＫＹ基因在基因组中的分布
根据Ｐｈｔｏｚｏｍｅｖ９．０电子定位结果显示，４９ ｙ个谷子ＳｉＷＲＫＹ基因分别定位于谷子的９条染色，该图谱显示Ｓ体上（图１）ｉＷＲＫＹ基因在基因组中的分布是非随机的。其中，在５号染色体上分布有１每个遗传距离分布有０７个基因，．３６个基因；在４号染色体上分布的基因有４个，每个遗传距离分布有０．１１个基因。
串联重复是同一条染色体上以基因簇形式存在的３我们发现许多Ｓ０ｋｂ以内的染色体重复，ｉ－图１中ＷＲＫＹ基因以２～５个基因簇的形式存在（
。最大的以５个串联重复存在的基红色框子所示）
因簇位于第７号染色体上，总共有３０个ＳｉＷＲＫＹ基因被发现存在串联重复现象，表明这种串联重复现象对于谷子ＷＲＫＹ基因家族的扩增具有重要
作用。
０８４
山西农业大学学报（自然科学版）
０１６２
２．３ ＳｉＷＲＫＹ基因的同源进化分析为了分析谷子ＷＲＫＹ基因与拟南芥中的同源基因的关系，我们采用了多序列比对的方法对谷子中的９４个ＳｉＷＲＫＹ基因和拟南芥的７２个Ａｔ－
２．４ ９４个ＳｉＷＲＫＹ蛋白序列保守性分析
ＷＲＫＹ转录因子是由位于Ｎ末端带有ＫＹＧＱＫ基序的６０个高度保守的氨基酸和位ＷＲ
［，，９］
。于Ｃ末端的Ｃ２ＨＣ２Ｈ２的锌指结构域２３１－或Ｃ为了分析９４个谷子ＳｉＷＲＫＹ基因成员在ＷＲＫＹ基序和锌指结构域的保守和差异，我们采用在线分析软件ＭＥＭＥ４．８．０对９４个ＳｉＷＲＫＹ蛋白的
结果发现ＫＹ和锌指核心结构域进行了分析，ＷＲ
ＫＹＧＱＫ结构域在不同的谷子ＷＲＫＹ家族成ＷＲ
员间是高度保守的，仅有个别序列差异（ＷＫＫＹ－。而与之相反，锌指结构域ＧＥＫ和ＷＲＫＹＧＫＫ）
。）在不同的成员之间则存在较大的变异（图３
序列分析发现谷子ＷＲＫＹ基因的氨基酸序每个结构域列保守域中存在１５个保守的结构域，包含的氨基酸残基不同，在１其中保１～５０个之间；守域１、而其他保守域２和３是每个成员都具有的，。）则在一些成员中缺失，可能与其功能有关（表１．５ ＳｉＷＲＫＹ基因的组织特异性电子表达谱２
基因的组织特异性表达常常被认为是基因在该组织中具有特定功能的一个标志，由于ＷＲＫＹ基因与植物的抗逆性有关，而谷子又是一个耐旱性较高的作物，因此，我们首先关注了ＳｉＷＲＫＹ基因在根系中的表达，随后也分析了其在穗子、叶片和茎秆中的表达，采用ＴＩＧＲ ＭｕｌｔｉＥｘｅｒｉｍｅｎｔ ｐ
）绘制了ＳＶｉｅｗｅｒ软件（ｅＶ４．９．０ｉＷＲＫＹ基因在Ｍ。结果发现，不同组织器官中的表达模式图（图４）０个ＳｉＷＲＫＹ基因在根系中都有很高的表达。８
很有意思的是没有ＳｉＷＲＫＹ基因在四种组织中组成型表达。Ｓ５３ｉＷＲＫＹ１和ＳｉＷＲＫＹ１在叶片 组织中高表达；ＳｉＷＲＫＹ８６和ＳｉＷＲＫＹ８７在茎中高表达；Ｓ４，６，３ｉＷＲＫＹ４，５４，５７和６７在穗中特异。）表达（图４
ＷＲＫＹ基因进行了同源进化分析，并采用
所ＧＡ５．０软件构建了系统进化树。结果显示，ＭＥ
有的拟南芥和谷子的ＷＲＫＹ基因可以被聚类为）。其中，３个同源群（ｒｏｕ１，２和３ｒｏｕ１和２ＧＧｐｐ 中的ＷＲＫＹ蛋白具有高度保守的ＷＲＫＹ结构域，在同一分支中存在谷子和拟南芥家族的成员，。与此相反，说明他们可能具有相同的功能（图２）但是，ｒｏｕ３中的成员虽然被分在同一个类群，Ｇｐ
在小的分支上仍然存在差异，谷子和谷子的聚为一类，拟南芥和拟南芥的聚在一个分支中，表明ｒｏｕ３中的ＷＲＫＹ基因在单子叶和双子叶植物Ｇｐ
间存在很大差异，说明Ｇｒｏｕ３中的成员在不同的ｐ。）
植物间也可能存在功能差别（图２
图２ 拟南芥（ＡｔＷＲＫＹｓ）和谷子（ＳｉＷＲＫＹｓ）
ＷＲＫＹ蛋白的同源进化树
Ｆｉ．２ ＰｈｌｏｅｎｅｔｉｃｔｒｅｅｏｆＷＲＫＹｒｏｔｅｉｎｓｆｒｏｍ ｙｇｐｇ
）ａｆｏｘｔａｉｌｍｉｌｌｅｔ（ＳｉＷＲＫＹｓｎｄｒａｂｉｄｏｓｉｓＡ ｐ
（）ＡｔＷＲＫＹｓ．
序列分析采用ＭＥＧＡ５．０软件的ＣｌｕｓｔａｌＷ程序， 注： 进化树采用邻接法构建，蛋白被分为３个类群，不同的颜色），，，，和亚群（代表不同的类群。群（Ｉ和ＩＩＩＩａＩＩｂ，ＩＩｃＩＩｄＩＩＩ）在圈外显示ＩＩｅ
：ＴｏｔｅｈｅｓｅｕｅｎｃｅｓｗｅｒｅａｌｉｎｅｄｂＣｌｕｓｔａｌＷａｔＮ ｑｇｙ ＭＥＧＡ５ａｎｄｔｈｅｈｌｏｅｎｅｔｉｃｔｒｅｅｗａｓｃｏｎｓｔｒｕｃｔｅｄｂ ｐｙｇｙｎｅｉｈｂｏｒｏｉｎｉｎｍｅｔｈｏｄ．Ｔｈｅｒｏｔｅｉｎｓｗｅｒｅｃｌａｓｓｉｆｉｅｄｉｎｔｏ－ ｇｊｇｐ ｔｈｒｅｅｄｉｓｔｉｎｃｔｃｌｕｓｔｅｒｓ．Ｅａｃｈｒｏｕｗａｓａｓｓｉｎｅｄａｄｉｆｆｅｒｅｎｔ ｇｐｇ ｃｏｌｏｒ．Ｔｈｅｎａｍｅｏｆｒｏｕｓ（ｎｄⅢ）ａｎｄｓｕｂｒｏｕ Ⅰ，Ⅱａ ｇｐｇｐ，Ⅱｅ）（，Ⅱｂ，Ⅱｃ，Ⅱｄｗｅｒｅｓｈｏｗｎａｔｔｈｅｏｕｔｓｉｄｅｏｆｔｈｅ Ⅱａｃｉｒｃｌｅ
Ｓ．６ｉＷＲＫＹ家族基因在干旱和盐胁迫下的表达２
前人对ＷＲＫＹ基因在干旱胁迫下的功能已
２０，２１］
。为了分析Ｓ经进行了深入的研究［ｉＷＲＫＹ
基因对于干旱胁迫的响应，我们随机选取了１０个其中包括ＧＳＳ１和ｉＷＲＫＹ基因，ｒｏｕ１（ｉＷＲＫＹｐ
，Ｇ，２Ｓ４）１ｉＷＲＫＹｒｏｕ２ａ（ＳｉＷＲＫＹ８）ｂｐ
（），，Ｓ４ＳＳ４ｉＷＲＫＹ１ｃ（ｉＷＲＫＹ８９）ｄ（ｉＷＲＫＹ３２２，，和Ｓ以及Ｇ５Ｓ１ｉＷＲＫＹ０）ｅ（ｉＷＲＫＹ３）ｒｏｕ３２ｐ（。我们将苗期的谷子材料Ｓ２ｉＷＲＫＹ７，８７和９２）经２０％ＰＥＧ处理后，１ｈ、４ｈ、６ｈ、１２ｈ和２４ｈ０ｈ、后，分别取样，采用ｑＰＴ－ＣＲ检测了上述１０个Ｒ
（）３６１２邢国芳等：谷子ＷＲＫＹ转录因子基因家族分析
１８４
图３ 不同类群的ＳｉＷＲＫＹ蛋白的氨基酸基序保守域分析
Ｆｉ．３ ＳｃｈｅｍａｔｉｃｄｉａｒａｍｏｆａｍｉｎｏａｃｉｄｍｏｔｉｆｓｏｆＳｉＷＲＫＹｒｏｔｅｉｎｓｆｒｏｍｄｉｆｆｅｒｅｎｔｒｏｕｓ ｇｐｇｐｇ
黑线代表相关蛋白的长度，不同颜色的框子代表ＷＲ选择ＷＲＫＹ蛋白信息在左侧标示，ＫＹ蛋白的不同结构域和位置 注：
：ＷＲＫＹｒｏｔｅｉｎＮｏｔｅＴｈｅｓｅｌｅｃｔｅｄＷＲＫＹｒｏｔｅｉｎｓａｒｅｌｉｓｔｅｄｏｎｔｈｅｌｅｆｔ．Ｔｈｅｂｌａｃｋｓｏｌｉｄｌｉｎｅｒｅｒｅｓｅｎｔｓｔｈｅｃｏｒｒｅｓｏｎｄｉｎ ｐｐｐｐｇ ａｎｄｉｔｓｌｅｎｔｈ．Ｔｈｅｄｉｆｆｅｒｅｎｔｃｏｌｏｒｅｄｂｏｘｅｓｒｅｒｅｓｅｎｔｄｉｆｆｅｒｅｎｔｍｏｔｉｆｓａｎｄｔｈｅｉｒｏｓｉｔｉｏｎｉｎｅａｃｈ ＷＲＫＹｓｅｕｅｎｃｅ － ｇｐｐｑ
基因的表达，结果发现，ＳｉＷＲＫＹ１和ＳｉＷＲＫＹ８６的表达量在２开始增加，然后０％ＰＥＧ处理１ｈ后，逐渐减少。而Ｓ４、Ｓ１ｉＷＲＫＹｉＷＲＫＹ３、Ｓｉ－ＷＲＫＹ２４５７、ＳｉＷＲＫＹ１、ＳｉＷＲＫＹ０和Ｓｉ－
活。另外，Ｓ５ｉＷＲＫＹ０在盐胁迫和干旱胁迫下表。）现出相同的表达模式（图５
３ 讨论
ＷＲＫＹ转录因子作为一种重要的植物生长调
控因子，在植物的生长、发育和对逆境胁迫的方面起着重要的调控作用。因此，人们在不同作物中对其开展了深入研究。本文对谷子中的ＷＲＫＹ家族成员进行了生物信息学分析，总共分离鉴定了０３个ＳｉＷＲＫＹ家族成员。并对拟南芥和谷子的１
构建了系统进化树。ＫＹ蛋白序列进行了比对，ＷＲ
，）。在将其分为３个大的类群（ｒｏｕＩＩＩ和ＩＩＩＧｐ
占拟南芥ｒｏｕ４个拟南芥基因，ＩＩＩ中仅有１Ｇｐ ＫＹ家族的１５％。而有４０个谷子ＷＲＫＹ家ＷＲ
族成员，占４２％。ＧｒｏｕＩＩＩ中的这种扩增现象可ｐ 能与进化过程中基因的复制有关。再者，这种扩增
ＷＲＫＹ８９的表达量在２０％ＰＥＧ处理后持续增加，
在１Ｓ８２ｈ的表达量达到最高。其中，ｉＷＲＫＹ９的。）表达量增加最为明显（图５
为了分析这些基因对于盐胁迫的响应，我们同样采用ｑＰＴ－ＣＲ检测了上述基因在盐胁迫下的Ｒ表达，结果发现除Ｓ４和Ｓ９ｉＷＲＫＹｉＷＲＫＹ２外，其他基因在盐胁迫下的表达水平均低于干旱胁迫下的表达，而Ｓ８ｉＷＲＫＹ７在盐胁迫后下调表达，
Ｓ４在盐胁迫１ｉＷＲＫＹ２ｈ后表达量上升了６倍，
在盐胁Ｓ１作为干旱胁迫的早期响应基因，ｉＷＲＫＹ迫１Ｓ８２ｈ后上调表达。与之相反，ｉＷＲＫＹ９作为干旱胁迫的后期响应基因，在盐胁迫１ｈ后就被激
２８４
山西农业大学学报（自然科学版）表１ ＳｉＷＲＫＹ蛋白的保守域分析
Ｔａｂｌｅ１ ＴｈｅｃｏｎｓｅｒｖｅｄｍｏｔｉｆｓｏｆｔｈｅＳｉＷＲＫＹｒｏｔｅｉｎｓ ｐ
０１６２
保守域Ｎｏ．ＣｏｎｓｅｒｖｅｄｍｏｔｉｆＭｏｔｉｆ１ Ｍｏｔｉｆ２ Ｍｏｔｉｆ３
Ｅｖａｌｕｅ－
最佳匹配氨基酸类型
ＢｅｓｔＭａｔｃｈｉｎＡｍｉｎｏＡｃｉｄＴｅｓ ｇｙｐ
［］［Ｎ］ＰＩＬＤＤＧＹ［ＲＱ］ＷＲＫＹＧＱＫｘ［ＩＶ］Ｋ［ＧＮ］Ｓ］［ＰＹＰＲＳＹ［ＹＦ］Ｋ］Ｃ［ＴＳＲ
［［［ＳＳＤＤＰＳＭ［ＬＶＦ］Ｖ］ＴＹＥＧＥＨ［ＴＮ］ｘｘＰＩＶＴ］ＨＣ］ ［［［ＱＧＣＰＡＶ］ＫＴ］ＫＱＶ［ＱＥ］ＲＲＱ］
氨基酸数目／ａａ
Ｎｏｍｂｅｒｏｆ ｔｈｅｒｏｔｅｉｎ ｐ
２９２１１１
ＳｉＷＲＫＹｓ
２７．３ｅ１７８－７５９３．３ｅ－９．０ｅ５２３－
ａｌｌＳｉＷＲＫＹｓｍｅｍｂｅｒｓ ａｌｌＳｉＷＲＫＹｓｍｅｍｂｅｒｓ ａｌｌＳｉＷＲＫＹｓｍｅｍｂｅｒｓ
ＳｉＷＲＫＹ１０，１１，１７，１８，３２，１９，
Ｍｏｔｉｆ４ ５０１４．９０Ｅ－
［［［ＶＲ［ＥＱ］ＰＲ［ＦＩＶ］Ｖ］Ｔ［ＲＫＴ］Ｓ［ＥＤ］ＡＶ］ＦＱＭ］ＶＤ
１５
３３，３６，５０，５１，５２，５４，７０，５，８，９４，９５，１０１，２２，２４，２７，２８，２９，４５，４７，６１，６６，８２，８９
［］］［［［［［［ＬＩＫ］ＩＱ］Ｈ［ＱＰ］Ｌ］Ｌ］Ｌ［ＱＨ］Ｙ］ＤＩＫＳＰＲＦ
Ｍｏｔｉｆ５
８．２０Ｅ１３３－
］［］［［［［［［ＳＱ］Ｐ］ＴＴ［ＮＱ］ＬＡＤＳＶＱＨ］ＱＨ］ＬＨ］ＳＡ［［［ＱＨ］Ｎ］ＬＬ［ＮＳ］Ｄ］Ａ［ＬＭ］ＲＡＬ［ＨＮ］ＬＡＬＳＮ［］［［［［［ＣＳＶＭ［ＮＫ］Ｑ］Ｐ］Ａ］Ａ［ＨＡ］Ｇ］ＰＱＨ］ＴＳＣ］［［［［［［［Ｑ［ＲＩＳＳＹ］Ｃ］Ｎ］Ｎ［ＮＧ］Ｇ］Ａ］ＲＤＧＨ］ＴＡ］ＩＳ
Ｍｏｔｉｆ６
７．９０Ｅ１０９－
］［［［［［［［Ｎ］ＩＱ］Ｅ］Ｑ［ＥＰ］Ｐ］ＱＡ］Ｍ［ＧＮ］ＤＡＭＫ］ＩＲＱ［］［［［ＧＶ］ＬＬＰＰＲ［ＡＬ］Ｅ［ＶＧＩＳ［ＴＡ］Ｌ］Ｐ］ＶＦＭＦ］［］［］［ＤＴＭＳＷＣ］ＥＥ［ＰＴ］ＳＶＬＳＮＳＦ
ＳｉＷＲＫＹ１８，３２，３３，３６，７０，９０， ９４，１０１，６１，６４，６６，８２
ＳｉＷＲＫＹ１８，３２，３３，３６，７０，９４， １０１，６１，６４，６６，８２
５０
ＳｉＷＲＫＹ７３，８１，８３，８４，８５，８６， ８８
５０
ＳｉＷＲＫＹ７３，８１，８３，８４，８５，８６， ８８
Ｍｏｔｉｆ７ ６．５０Ｅ１０７－
［］］］［［ＰＳＳＥ］ＤＧＹ［ＮＳＷＲＫＹＧＱＫＱＶＫＧＳＡＤ ［］ＳＬＤＧ［ＱＨ］ＩＴＥＩＶＩＹＫＧ［ＡＴ］ＨＮＨＰ［ＫＲ］Ｐ［ＱＰ］ＮＴＲＲ
［［［［ＬＴ］ＨＬＦＳＰＲ］Ｓ［ＡＤ］ＲＭ］ＰＧＡＥ［ＳＧ］ＳＳＲＳＮＨ］
１９
Ｍｏｔｉｆ８ ９１．９０Ｅ８－２５
Ｍｏｔｉｆ９ １．２０Ｅ９２－
［］［［［［ＡＡＧＤＶＧ［ＧＡ］ＩＴＭ］Ｒ］Ｑ［ＱＰＲ］Ｒ］Ｒ］ＳＳＫＫＧ］［ＫＲＲＲ［ＡＳＮ［ＥＡ］ＥＴＳＷＱＲ］ＶＩＬＴ［ＥＤ］Ａ
５０
ＳｉＷＲＫＹ８３，８４，８５，８６，８８
Ｍｏｔｉｆ１０ ７６２．３０Ｅ－
］Ｅ［ＮＨ］ＰＲＳＹＹＫＣＴ［ＨＦ］Ｐ［ＮＳＣＰ［ＶＴ］ＫＫＫＶＥ ］［［［［ＨＣＳＫ［ＫＲ］ＲＫＮ［ＲＱ］Ｖ［ＫＲ］Ｋ］Ｓ］ＩＲ］ＲＶＶＣ［］ＶＩＰ［ＡＶ］ＩＳ
［［ＬＳＰＮＨＡＧＹ［ＨＬ］ＴＳＴＤＤＧＡＳＤＦＨ［ＹＣ］Ｇ］ＭＴ］Ｅ［ＳＴＤＧ［ＹＴ］Ｙ］Ｉ
［］］［［ＱＧ］Ｃ［ＥＫ］ＬＶ［ＡＱＳＱ［ＬＭ］ＦＱ［ＤＮＳＩＬＦ］ＲＱ］［［［［Ｆ］Ｌ［ＱＥ］Ｌ［ＧＨ］ＤＣＳＳＫ［ＡＶ］ＩＳＴ］ＥＤＨ］ＮＱ］［［［［［Ａ［ＬＶ］ＫＶ］ＥＳＴ］ＥＬ［ＲＳ］Ｅ］ＩＭ］ＥＥＥＡＲＶＧＫ］［］［［［ＮＶ］Ｒ［ＲＱ］Ｌ［ＡＥ］Ｇ］ＩＥＨＭ］Ｌ
［］［［［［［ＳＰ］Ｙ］ＬＶ］ＳＶ］Ｐ［ＰＡ］ＧＬＳＰ［ＡＴ］ＥＴ］ＰＦＴＩＡ［ＬＬ［ＤＥ］ＳＰ［ＶＧ］Ｆ］ＬＬ
２１
ＳｉＷＲＫＹ１８，３３，３６，１０１，６１，６６， ８２
ＳｉＷＲＫＹ１３，１４，３７，１，７１，７２，４， ６，９２，９３，９６，９９，１００，４８，６３，６８
Ｍｏｔｉｆ１１ ３．２０Ｅ６０－２０
Ｍｏｔｉｆ１２ １．１０Ｅ３８－２９
ＳｉＷＲＫＹ８３，８４，８５，８６
Ｍｏｔｉｆ１３ ８．２０Ｅ３８－２９
ＳｉＷＲＫＹ２，７５，７６，７９，２１，８７
Ｍｏｔｉｆ１４ ２．４０Ｅ３１－２０
ＳｉＷＲＫＹ１３，１４，１５，１６，３８，３， ７５，７，７９，２５，２６，４２，４３，４４，４８ＳｉＷＲＫＹ１７，１８，３３，３６，７０，１０１， ６１，６６
Ｍｏｔｉｆ１５ １．８０Ｅ２７－２１
（）３６１２邢国芳等：谷子ＷＲＫＹ转录因子基因家族分析
３８４
图４ Ｓ叶片、根系、茎和穗）的表达ｉＷＲＫＹ基因在４种谷子组织中（
，，ｓｈｏｗｉｎｔｈｅｅｘｒｅｓｓｉｏｎａｔｔｅｒｎｏｆｌｅａｆｒｏｏｔｓｔｅｍａｎｄｓｉｃａｉ．４ｅａｔａｉＷＲＫＹｇｎｅｓｉｎｆｏｕｒｔｉｓｓｕｅｓｎａｍｅｌＦｅＳ Ｈ －ｍｐｐｇｐｐｙｇ
右边颜色梯度表示表达倍数的变化。蓝色－红色表达量升高 注：
：Ｎ））ｏｔｅｏｔｅｔｈａｔｅｘｒｅｓｓｉｏｎｖａｌｕｅｓｍａｅｄｔｏａｃｏｌｏｒｒａｄｉｅｎｔｆｒｏｍｌｏｗ（ｂｌｕｅｔｏｈｉｈｅｘｒｅｓｓｉｏｎ（ｏｒａｎｅＮ ｐｐｐｇｇｐｇ
现象或多或少的与ＳｉＷＲＫＹ基因的串联重复有关。我们总共检测到有３０个基因为串联重复序。列，并且最大的串联重复基因来自于ＧｒｏｕＩＩＩｐ 这种串联重复现象在玉米和二穗短柄草中同样存
２２］
，表明ＷＲＫＹ基因的连续复制在禾本科作物在［
中可能是一种普遍事件。
对于解析抗旱机制，进行遗传改良具有重要意义。本研究发现，ＳｉＷＲＫＹ基因与谷子的抗旱性密切相关，大部分家族成员对于干旱胁迫能够快速响应，同样印证了前人在其他作物中的研究结果。例如，超表达Ｏ４ｓＷＲＫＹ５基因提高
２３］
，同样，水稻抗旱耐盐的能力［超表达大豆
［４］
，ＷＲＫＹ基因提高了转基因拟南芥抗旱能力２
前人已经在不同的物种中分析了ＷＲＫＹ
转录因子对于干旱胁迫的响应，但这些作物本身并不是耐旱作物。因此，作为抗旱耐瘠薄的模式植物，发掘谷子自身的ＷＲＫＹ基因资源，
本研究发现与拟南芥ＡｔＷＲＫＹ５７直系同源的
谷子Ｓ虽９基因在盐胁迫后差异表达，ｉＷＲＫＹ然间接证明了该基因参与谷子对逆境胁迫的响
４８４
山西农业大学学报（自然科学版）
０１６２
图５ １０个ＳｉＷＲＫＹ基因在干旱和盐胁迫后的表达
，，，ＦｅＳｉ．５ｈｅｒｅｌａｔｉｖｅｅｘｒｅｓｓｉｏｎｒａｔｉｏｏｆ１０ｉＷＲＫＹｇｎｅｓａｎａｌｚｅｄｗｉｔｈＲＴＣＲａｔ０１６１２ａｎｄ２４ｈｏｆｄｅｈｄｒａｔｉｏｎ－Ｐ Ｔ ｐｙｑｙｇ
应，但是要证实其对抗逆境胁迫的能力仍然需要进一步的功能分析数据。
电子表达谱的数据结果进一步暗示谷子
是，超过８０％的ＳｉＷＲＫＹ基因在根中高度表达，其中９个基因在根中特异表达。而根系在抵抗非
２５］
。生物胁迫中具有重要功能［
ＷＲＫＹ基因对于干旱胁迫具有重要作用。特别
（）３６１２邢国芳等：谷子ＷＲＫＹ转录因子基因家族分析
５８４
参 考 文 献
［］：［］１ＬａｔａＣ，ＧｕｔａＳ，ＰｒａｓａｄＭ．ＦｏｘｔａｉｌｍｉｌｌｅｔａｍｏｄｅｌｃｒｏｆｏｒｅｎｅｔｉｃａｎｄｅｎｏｍｉｃｓｔｕｄｉｅｓｉｎｂｉｏｅｎｅｒｒａｓｓｅｓＪ．ＣｒｉｔｉｃａｌＲｅｖｉｅｗｓｉｎＢｉｏ－ ｐｐｇｇｇｙｇ ，ｔ３ｅｃｈｎｏｌｏ２０１３，３３：３２８４３．－ｇｙ
［］，［］，２ＲｕｓｈｔｏｎＰＪＳｏｍｓｓｉｃｈＩＥ，ＲｉｎｌｅｒＰ，ｅｔａｌ．ＷＲＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｓＪ．ＴｒｅｎｄｓｉｎＰｌａｎｔＳｃｉｅｎｃｅ２０１０，１５：２４７２５８． －ｇｐ
［］，［］，２０６．３Ｅｕｌｅｍ Ｔ，ＲｕｓｈｔｏｎＰＪＲｏｂａｔｚｅｋＳ，ｅｔａｌＴｈｅＷＲＫＹｓｕｅｒｆａｍｉｌｏｆｌａｎｔｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｓＪ．ＴｒｅｎｄｓＰｌａｎｔＳｃｉ２０００，５：１９９－ ｇｐｙｐｐ ［］４ＤｅｖａｉａｈＢＮ，ＫａｒｔｈｉｋｅａｎＡＳ，ＲａｈｏｔｈａｍａＫ Ｇ．ＷＲＫＹ７５ｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｉｓａｍｏｄｕｌａｔｏｒｏｆｈｏｓｈａｔｅａｃｕｉｓｉｔｉｏｎａｎｄｒｏｏｔｄｅｖｅｌｏ－ ｙｇｐｐｐｑｐ
［］，Ｊ．ｍＡ１ｅｎｔｉｎｒａｂｉｄｏｓｉｓＰｌａｎｔＰｈｓｉｏｌｏ２００７，１４３：１７８９８０１．－ ｙｇｙｐ
［］，ｉ５ＺｈａｎＪＰｅｎＹＬ，ＧｕｏＺＪ．ＣｏｎｓｔｉｔｕｔｉｖｅｅｘｒｅｓｓｉｏｎｏｆａｔｈｏｅｎｎｄｕｃｉｂｌｅＯｓＷＲＫＹ３１ｅｎｈａｎｃｅｓｄｉｓｅａｓｅｒｅｓｉｓｔａｎｃｅａｎｄａｆｆｅｃｔｓｒｏｏｔ－ ｇｇｐｐｇ ［］，５２１．ｒｏｗｔｈａｎｄａｕｘｉｎｒｅｓｏｎｓｅｉｎｔｒａｎｓｅｎｉｃｒｉｃｅｌａｎｔｓＪ．ＣｅｌｌＲｅｓｅａｒｃｈ２００８，１８：５０８－ ｇｐｇｐ
［］，６ＲｏｂａｔｚｅｋＳ，ＳｏｍｓｓｉｃｈＩＥ．Ａｎｅｗ ｍｅｍｂｅｒｏｆｔｈｅＡＷＲｒａｂｉｄｏｓｉｓＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｆａｍｉｌＡｔＷＲＫＹ６，ｉｓａｓｓｏｃｉａｔｅｄｗｉｔｈｂｏｔｈｓｅ － ｐｙｐ［］，ｎｒ１ｅｓｃｅｎｃｅｎｄｄｅｆｅｎｃｅｅｌａｔｅｄｒｏｃｅｓｓｅｓＪ．ＴｈｅＰｌａｎｔＪｏｕｒｎａｌ２００１，２８：１２３３３．－ａ－－ ｐ
［］７üｌｋｅｒＢ，ＭｕｋｈｔａｒＭ Ｓ，ＳｏｍｓｓｉｃｈＩＥ．ＴｈｅＷＲＫＹ７０ｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｏｆｒａｂｉｄｏｓｉｓｎｆｌｕｅｎｃｅｓｂｏｔｈｔｈｅｌａｎｔｓｅｎｅｓｃｅｎｃｅａｎｄｄｅｆｅｎｓｅＡｉ ｐｐｐ［］，ｓｉｎａｌｉｎａｔｈｗａｓＪ．Ｐｌａｎｔａ２００７，２２６，１２５１３７．－ｇｇｐｙ
［］（）］：李摇冉，娄永根．植物中逆境反应相关的ＷＲ３８ＫＹ转录因子研究进展［０１１，３１１１２２３２３１Ｊ．生态学报，３２－
［］９ＴｒｉａｔｈｉＰ，ＲａｂａｒａＲＣ，ＲｕｓｈｔｏｎＰＪ．ＡｓｓｔｅｍｓｂｉｏｌｏｅｒｓｅｃｔｉｖｅｏｎｔｈｅｒｏｌｅｏｆＷＲＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｓｉｎｄｒｏｕｈｔｒｅｓｏｎｓｅｓｉｎ ｐｙｇｙｐｐｐｇｐ ［］，（）：２５５６６．ｌａｎｔｓＪ．Ｐｌａｎｔａ２０１４，２３９２２－ｐ
［］，ｂ１０ＩｓｈｉｕｒｏＳ，ＮａｋａｍｕｒａＫ．ＣｈａｒａｃｔｅｒｉｚａｔｉｏｎｏｆａｃＤＮＡｅｎｃｏｄｉｎａｎｏｖｅｌＤＮＡｉｎｄｉｎｒｏｔｅｉｎＳＰＦ１，ｔｈａｔｒｅｃｏｎｉｚｅｓＳＰ８ｓｅｕｅｎｃｅｓｉｎ－ ｇｇｇｐｇｑ
’［］，ｔｈｅ５ｕｓｔｒｅａｍｒｅｉｏｎｓｏｆｅｎｅｓｃｏｄｉｎｆｏｒｓｏｒａｍｉｎａｎｄｂｅｔａｍｌａｓｅｆｒｏｍｓｗｅｅｔｏｔａｔｏＪ．ＭｏｌｅｃｕｌａｒａｎｄＧｅｎｅｒａｌＧｅｎｅｔｉｃｓ１９９４，ａ －ｐｇｇｇｐｙｐ （）：２４４６５６３７１．５－
［］，：［］１１ＲｉｅｃｈｍａｎｎＪＬ，ＨｅａｒｄＪＭａｒｔｉｎＧ，ｅｔａｌ．ＡｒａｂｉｄｏｓｉｓｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｓｅｎｏｍｅｉｄｅｃｏｍａｒａｔｉｖｅａｎａｌｓｉｓａｍｏｎｅｕｋａｒｏｔｅｓＪ．－ｗ ｐｐｇｐｙｇｙ
，（）：Ｓｃｉｅｎｃｅ２０００，２９０５４９９１０５１１０．２２－
［］ｉｎｒｉｃｅｕｎｄｅｒｖａｒｉｏｕｓｏｆｔｈｅＷＲＫＹｅｎｅｆａｍｉｌ１２ＲａｍａｍｏｏｒｔｈＲ，ＪｉａｎＳＹ，ＫｕｍａｒＮ，ｅｔａｌ．Ａｃｏｍｒｅｈｅｎｓｉｖｅｔｒａｎｓｃｒｉｔｉｏｎａｌｒｏｆｉｌｉｎ ｇｙｇｐｐｐｇｙ
，（）［］：ａｂｉｏｔｉｃａｎｄｈｔｏｈｏｒｍｏｎｅｔｒｅａｔｍｅｎｔｓＪ．ＰｌａｎｔＣｅｌｌＰｈｓｉｏｌｏ２００８，４９６８６５７９．８ －ｐｙｙｇｙ
［］，，［］，（）：１３ＳｃｈｍｕｔｚＪＣａｎｎｏｎＳＢ，ＳｃｈｌｕｅｔｅｒＪｅｔａｌ．ＧｅｎｏｍｅｓｅｕｅｎｃｅｏｆｔｈｅａｌａｅｏｏｌｌｏｉｄｓｏｂｅａｎＪ．Ｎａｔｕｒｅ２０１０，４６３７２７８７８８３．１１ －ｑｐｐｙｐｙ［］，ｉｎ１４ＷｅｉＫＦ，ＣｈｅｎＪＣｈｅｎＹＦ，ｅｔａｌ．ＭｏｌｅｃｕｌａｒｈｌｏｅｎｅｔｉｃａｎｄｅｘｒｅｓｓｉｏｎａｎａｌｓｉｓｏｆｔｈｅｃｏｍｌｅｔｅＷＲＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｆａｍｉｌ ｐｙｇｐｙｐｐｙ
［］，（）：ｍａｉｚｅＪ．ＤＮＡ Ｒｅｓｅｒｃｈ２０１２，１９２６５３４．１－
［］１５ＷｅｎＦ，ＺｈｕＨ，ＬｉＰ，ｅｔａｌ．ＧｅｎｏｍｅｃｈａｒａｃｔｅｒｉｚａｔｉｏｎａｎｄｅｘｒｅｓｓｉｏｎａｎａｌｓｅｓｏｆＷＲＫＹｆａｍｉｌｅｎｅｓｉｎｂｒａｃｈｏｄｉｕｍｉｄｅｅｖｏｌｕｔｉｏｎａｒ －ｗｐｙｙｇｙｐｙ
［］（）：ｄｉｓｔａｃｈｏｎＪ．ＤＮＡ Ｒｅｓｅａｒｃｈ，２０１４，２１３３２７９３－ｙ
［］，［］４ｉｎｌａｎｔｅｎｏｍｅｓＪ．Ｓｃｉｅｎｃｅ２００８，３２０：４８６ａｎｄｃｏｌｌｉｎｅａｒｉｔＸ，ｅｔａｌ．Ｓｎｔｅｎ１６ＴａｎＨ，ＢｏｗｅｒｓＪＥ，Ｗａｎ８８．－ ｐｇｙｙｙｇｇ
［］，１７ＰｕｒａｎｉｋＳ，ＳａｈｕＰＰ，ＭａｎｄａｌＳＮ，ｅｔａｌ．ＣｏｍｒｅｈｅｎｓｉｖｅｅｎｏｍｅｏｆｔｈｅＮＡＣｉｄｅｓｕｒｖｅｅｎｏｍｉｃｃｏｎｓｔｉｔｕｔｉｏｎａｎｄｅｘｒｅｓｓｉｏｎｒｏｆｉｌｉｎ －ｗｐｇｙｇｐｐｇ
）［］ｉｎｆｏｘｔａｉｌｍｉｌｌｅｔ（ｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｆａｍｉｌＳｉＬ．ｅｔａｒｉａｔａｌｉｃａＰＬｏＳＯＮＥ，２０１３，８：ｅ６４５９４．Ｊ． ｐｙ
［］［］，ｒ３１８ＬａｎｍｅａｄＢ，ＳａｌｚｂｅｒＳＬ．ＦａｓｔａｅｄｅａｄａｌｉｎｍｅｎｔｗｉｔｈＢｏｗｔｉｅ２Ｊ．ＮａｔｕｒｅＭｅｔｈｏｄｓ２０１２，９：３５７５９．－－ ｇｐｐｇｇｇ
［］，１９ＲｕｓｈｔｏｎＰＪＴｏｒｒｅｓＪＴ，ＰａｒｎｉｓｋｅＭ，ｅｔａｌ．ＩｎｔｅｒａｃｔｉｏｎｏｆｅｌｉｃｉｔｏｒｎｄｕｃｅｄＤＮＡ－ｉｎｄｉｎｒｏｔｅｉｎｓｗｉｔｈｅｌｉｃｉｔｏｒｒｅｓｏｎｓｅｅｌｅｍｅｎｔｓｉｎｔｈｅｉｂ －ｇｐｐ
［］，５ＰＲ１ｇｅｎｅｓＪ．ＴｈｅＥＭＢＯＪｏｕｒｎａｌ１９９６，１５：５６９０ｒｏｍｏｔｅｒｓｏｆａｒｓｌｅ７００．－ ｐｐｙ
［］［］ｉＯ０ＡＪ．Ｓ，ＹｕＤ．ＯｖｅｒｅｘｒｅｓｓｉｏｎｏｆｔｈｅｓｔｒｅｓｓＹ，Ｊｉｎ２０ＳｏｎｎｄｕｃｅｄｓＷＲＫＹ８ｉｍｒｏｖｅｓｏｓｍｏｔｉｃｓｔｒｅｓｓｔｏｌｅｒａｎｃｅｉｎｒａｂｉｄｏｓｉｓＣｈｉｎｅｓｅ－ ｐｇｇｐｐ
，ＳｃｉｅｎｃｅＢｕｌｌｅｔｉｎ２０１０，５４：４６７１６７８．４ －
［］ａｇｒｏｕＩＷＲＫＹｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｏｆＹ Ｎ，ＤａｎＦＦ，ＬｉｕＺＱ，ｅｔａｌ．ＣａＷＲＫＹ５８，ｅｎｃｏｄｉｎ２１ＷａｎｎＣａａｓｉｃｕｍ ｎｎｕｕｍ，ｅａｔｉｖｅｌ ｐｇｐｇｇｇｙｐ
，［］Ｒｓ１ｒｅｕｌａｔｅｓｒｅｓｉｓｔａｎｃｅｔｏ２０１３，１４：１３１ａｌｓｔｏｎｉａｏｌａｎａｃｅａｒｕｍｉｎｆｅｃｔｉｏｎＪ．ＭｏｌｅｃｕｌａｒＰｌａｎｔＰａｔｈｏｌｏ４４． － ｇｇｙ
［］ａｎｄｆｕｎｃｔｉｏｎａｌｃｈａｒａｃｔｅｒ２２ＢｅｎｃｋｅａｌａｔｏＭ，ＣａｂｒｅｉｒａＣ，ＷｉｅｂｋｅｔｒｏｈｍＢ，ｅｔａｌ．Ｇｅｎｏｍｅｉｄｅａｎｎｏｔａｔｉｏｎｏｆｔｈｅｓｏｂｅａｎ ＷＲＫＹｆａｍｉｌＳ －Ｍ－－ｗ－ｙｙ
［］，ａｉＰｉｚａｔｉｏｎｏｆｅｎｅｓｉｎｖｏｌｖｅｄｉｎｒｅｓｏｎｓｅｔｏｈａｋｏｓｏｒａｃｈｒｈｉｚｉｎｆｅｃｔｉｏｎＪ．ＢＭＣＰｌａｎｔＢｉｏｌｏ２０１４，１４：２３６． ｇｐｇｙｐｐｙ
［］ｅｉＡｓ２３ＱｉｕＹＰ，ＹｕＤ Ｑ．ＯｖｅｒｘｒｅｓｓｉｏｎｏｆｔｈｅｓｔｒｅｓｓｎｄｕｃｅｄＯｓＷＫＲＹ４５ｅｎｈａｎｃｅｓｄｉｓｅａｓｅｒｅｓｉｓｔａｎｃｅａｎｄｄｒｏｕｈｔｔｏｌｅｒａｎｃｅｉｎｒａｂｉｄｏｓｉ－－ ｐｇｐ
，［］ＥｎｖｉｒｏｎｍｅｎｔａｌａｎｄＥｘｅｒｉｍｅｎｔａｌＢｏｔａｎ２００９，６５：３５７．４Ｊ． －ｐｙ
［］２４ＪｉａｎＹ Ｑ，ＤｅｈｏｌｏｓＭ Ｋ．ＦｕｎｃｔｉｏｎａｌｃｈａｒａｃｔｅｒｉｚａｔｉｏｎｏｆＡＮｒａｂｉｄｏｓｉｓａＣｌｉｎｄｕｃｉｂｌｅＷＲＫＹ２５ａｎｄＷＲＫＹ３３ｔｒａｎｓｃｒｉｔｉｏｎｆａｃｔｏｒｓｉｎ ｇｙｐｐ
［］，ＭＰ１ａｂｉｏｔｉｃｓｔｒｅｓｓｅｓＪ．ｌａｎｔｏｌｅｃｕｌａｒＢｉｏｌｏ２００９，６９：９１０５． － ｇｙ
［］：［］，２５ＧｈｏｓｈＤ，ＸｕＪ．ＡｂｉｏｔｉｃｓｔｒｅｓｓｒｅｓｏｎｓｅｓｉｎｌａｎｔｒｏｏｔｓａｐｒｏｔｅｏｍｉｃｓｅｒｓｅｃｔｉｖｅＪ．ＦｒｏｎｔｉｅｒｓｉｎＰｌａｎｔＳｃｉｅｎｃｅ２０１４，５：６． ｐｐｐｐ
（编辑：武英耀）