Identification and Pathogenic Response of Tomato WRKY Transcription Factors
-
摘要:目的 深入了解番茄WRKY转录因子的组学特征及其生物胁迫响应。方法 基于最新公共数据,利用生物信息学和比较基因组学方法对番茄WRKY进行系统鉴定,结合抗、感2个番茄自交系在青枯菌侵染前后的RNA-seq数据,挖掘青枯病抗性相关WRKY。结果 85个番茄WRKY转录因子被鉴定,可分为I、IIa+b、IIc、IId+e和III等类别,IIe基因最多。其中,9个基因七肽基序发生了单一氨基酸变异,WRKYGKK为优势突变型。这些WRKY主要分布在5号染色体,且具有端部和成簇分布现象,尤其是IIe亚类。45.88%的番茄WRKY具有共线性。58.82%的番茄WRKY(主要是I和IIc类)与拟南芥和辣椒WRKY形成73对直系同源基因,其选择压力(Ka/Ks)均小于1。16个番茄WRKY(主要是IIa+b和IIc类)对几种生物胁迫反应强烈,且主要在根中表达。12个差异表达WRKY(主要是III和IIb类)被鉴定,其中Solyc03g095770.3(III)与Solyc09g014990.4(I)互作在番茄青枯病响应中发挥重要作用。结论 综合鉴定了番茄WRKY转录因子,筛选到12个青枯病响应基因。Abstract:Objective Characteristics and biotic stress response of WRKY transcription factors (TFs) in tomato plants were investigated.Method The latest available bioinformatics and genomics methods were employed to identify the tomato WRKY TFs. RNA-seq of disease-resistant and susceptible tomato inbred lines before and after artificial Ralstonia solanacearum infection were obtained to identify the TFs associated with the pathogenic resistance of the plants.Result Eighty-five tomato WRKY TFs were identified and divided into I, IIa+b, IIc, IId+e, and III categories. The IIe group had the highest number of the TFs. The conserved motif of 9 TFs had one single amino acid variation, and WRKYGKK was the dominant mutant. The TFs, especially those in the IIe group, were mainly found on chromosome 5, at the ends, and in clusters. In them, 45.88% showed collinearity and 58.82% (mainly in I and IIc groups) formed 73 pairs of orthologs with those in Arabidopsis and chili pepper at a Ka/Ks ratio below 1. Sixteen of them, mainly belonging to IIa+b and IIc, responded significantly to the biotic stress with expressions largely in the roots. There were 12 differentially expressed WRKY TFs identified mainly in III and IIb. Of which, the interaction between Solyc03g095770.3 (III) and Solyc09g014990.4 (I) played a significant role in the response of the tomato plant to bacterial wilt.Conclusion The WRKY TFs were identified in tomato plants. Twelve genes responded to the bacterial wilt were isolated.
-
0. 引言
【研究意义】多花黄精(Polygonatum cyrtonema Hua)为百合科(Liliaceae)黄精属(Polygonatum)多年生草本植物。其根状茎肥厚,富含多糖、甾体皂苷、黄酮类等多种对人体有益物质,是我国传统药食同源的珍稀药材,具补气养阴、健润脾肺、益肾之功效,开发利用前景广阔[1-3]。多花黄精根状茎作为主要药用部位,生长缓慢,极易受到外界不良生长环境的影响,部分适宜生长地区土壤盐渍化日益严重逐年扩散,对多花黄精的生长和品质都有较大的影响,目前相关研究多局限于药理作用[4]、生理指标[5]、组织培养[6]等方面,对分子生物学方面的研究少之又少。因此对多花黄精根状茎生长发育及抗盐胁迫相关基因的深入研究具有非常重要的意义。【前人研究进展】转录因子(Transcription factor,TF)是能与真核基因启动子区域中顺式作用元件发生特异性结合,和下游靶基因顺式作用元件相互作用,进而调控靶基因表达的一类蛋白[7],在植物的生长发育、形态建成、抗逆及次生代谢物等方面具有重要的调控作用。Plethora(PLT)属于AP2/ERF型转录因子,含有两个 AP2 重复结构。PLT基因在胚胎发育、干细胞龛(静止中心,QC)、分生组织维护、器官生长发育等过程中具有核心作用[8];在盐胁迫、干旱胁迫、高温、低温胁迫等非生物胁迫条件下,调节相关功能基因的表达以提高抗逆性[9]。QC是干细胞的组织中心,植物干细胞具有自我更新和持续分裂子细胞的能力,是植物根、茎、叶和花等器官发生的源泉。PLT1 和 PLT2是干细胞维持静止中心所必需的两个基因,静止中心的正常形成要求PLT1、PLT2 、SHORT-ROOT(SHR)和SCARECROW(SCR )共同参与,PLT、SCR、植物特异性转录因子TCP的3个基因互作形成 PLT-TCP-SCR 复合体,与WOX5 启动子中的 PLT 结合位点结合进行诱导表达,在胚胎发生过程中特化干细胞龛,并在胚胎后维持干细胞属性[8,10]。PLT基因与生长素关联密切,参与较多的基因调控网络,转录响应生长素积累并依赖于生长素响应转录因子进行转录,通过生长素的局部合成[11]及PIN蛋白[12]作为介导定向运输转运生长素在QC 中形成了浓度梯度[13-14],作用于下游的PLT基因,形成一个调控网络来维持稳定的生长素浓度和根尖干细胞龛的位置[15],对根分生组织的大小和维持至关重要[16]。PLT3、PLT5 和 PLT7调节侧根原基的启动及建立正确的对称生长轴中基因表达程序[17-18]。Kareem等[19]揭示了PLT 控制不定芽再生的机制:第一步,PLT3、PLT5 和 PLT7 通过 PLT1 和 PLT2激活根干细胞建立多能性;第二步,PLT3、PLT5 和 PLT7通过增强因子 Cupshaped cotyledon2(CUC2)调节来完成芽形成。PLT活性高能维持干细胞活性,PLT活性低有利于干细胞进行有丝分裂,PLT活性更低则促进细胞分化[8]。PLT基因在根干细胞区域和不定芽的生成过程中都表现出特异性表达[19]。PLT基因不仅能响应种子萌发、生根等器官发育,响应激素信号,还能调控非生物胁迫[9,20]。【本研究切入点】PLT基因家族的功能研究大多集中在模式植物拟南芥(Arabidopsis thaliana L.)和水稻(Oryza sativa L.)中,且主要集中在根系中,目前国内有关PLT基因家族的研究较少,在多花黄精中更是未见报道,PLT基因家族是否调控多花黄精生长发育及行使相关功能还未曾揭示。【拟解决的关键问题】以实验室前期获得的多花黄精不同组织部位转录组数据为基础,对多花黄精PLT家族进行生物信息学分析并对关键成员进行克隆,测定其在不同组织部位和调控盐胁迫方面的表达情况,以期为多花黄精PLT家族基因的生物学功能研究、多花黄精遗传改良奠定基础。
1. 材料与方法
1.1 材料
野生植株不同组织部位(果实、叶片、根、茎秆、根状茎)转录组数据材料由采自福建省邵武市大竹镇(北纬27°11′6.80″,东经117°35′4.60″)野生多花黄精测序获得(另文发表),盐胁迫转录组数据由多花黄精组培苗经不同浓度NaCl处理24 h后采样测序获得。
1.2 方法
1.2.1 总RNA的提取及cDNA的合成
参考RNAprep Pure总RNA提取试剂盒(TIANGEN)说明书提取多花黄精样品总RNA,用Thermo超微量核酸检测仪测定RNA浓度,通过1%琼脂糖凝胶电泳检测RNA纯度,然后采用All supeiMix(TransGen Biotech公司)和GeneRacerTM试剂盒(TaKaRa)逆转录成cDNA用于PLT基因定量和克隆,用SYBR Green Master Mix(No Rox)进行qRT-PCR反应。
1.2.2 多花黄精 PLT 基因家族的鉴定
从 Pfam 数据库(http://pfam.xfam.org/)获取 PLETHORA(PF00847)结构域的隐马尔可夫模型文件,使用 TBtools软件对多花黄精转录组数据库进行搜索。使用 CDD(https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi)和SMART(http://smart.embl-heidelberg.de/)对候选序列进行保守结构域鉴定和筛选。采用在线网站 ExPASy(https://web.expasy.org/protparam/)分析多花黄精PLT成员的氨基酸数量、蛋白分子量、等电点、不稳定指数和亲水性。利用SignalP 4.1 Server(https://services.healthtech.dtu.dk/service.php?SignalP)、TMHMM Server 2.0(https://services.healthtech.dtu.dk/service.php?TMHMM-2.0)、 Cell-PLoc 2.0(http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/)、SOPMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)在线网站进行 PcPLT 蛋白的信号肽、跨膜结构、亚细胞定位和二级结构的分析与预测。
1.2.3 多花黄精 PLT 家族基因结构、保守基序和保守结构域分析
运用 GSDS(http://gsds.gao-lab.org/)在线网站分析 PcPLT基因结构。使用 MEME 在线软件(https://meme-suite.org/meme/tools/meme)分析该家族成员编码蛋白的保守基序。同时结合结构域分析结果,使用 TBtools 绘制基因结构、保守基序和保守结构域图。
1.2.4 PcPLT 家族系统进化关系分析
模式植物拟南芥(Arabidopsis thaliana)的 PLT 家族蛋白序列下载自 TAIR 数据库(https://www.arabidopsis.org),水稻(Oryza sative L.)的PLT 家族蛋白序列下载自 RiceData 数据库(http://www.ricedata.cn/gene),从NCBI [National Center for Biotechnology Information (nih.gov)]下载石刁柏(Asparagus officinalis L.)、油棕(Elaeis guineensis)、葡萄(Vitis vinifera)蛋白序列。使用 MEGA 11.0 软件中的最大似然法(Maximum likelihood,ML)构建系统发育树,并使用在线软件 iTOL(https://itol.embl.de/upload.cgi)对进化树进行美化。
1.2.5 PcPLT 家族GO功能注释
为了解到PcPLT 全部基因信息,进行GO数据库功能注释搜索,分析注释信息。
1.2.6 多花黄精PLT基因特异表达分析
为进一步了解PLT 基因在多花黄精不同组织器官中的表达模式,通过qRT-PCR检测PLT 基因在不同组织部位以及NaCl不同浓度处理下的表达情况。根据HieffTM qPCR SYBR Green Master Mix(YEASEN)说明书进行qRT-PCR反应。 反应体系为20 μL,包括Mix 10 μL,ddH2O 7.2 μL,cDNA 模版2 μL,上、下引物各0.4 μL。内参基因选用β-tubulin基因,依据2−△△Ct法分析计算PcPLT基因的相对表达量,并使用GraphPad Prism进行图表绘制。分析多花黄精PLT基因在不同组织器官中及不同浓度NaCl处理下的表达情况(引物序列见表1)。
表 1 引物序列Table 1. Sequences of primers used引物名称
Primer name引物序列
Primer sequence (5′-3′)退火温度Tm/℃ 用途 Usage PLT2-2-Q-F
PLT2-2-Q-RCGTCATCTTATCGTGGAGTTAC
CATCTTTCCGCCTTGTTACT55 qRT-PCR引物 PLT2-7-Q-F
PLT2-7-Q-RTCAAGTGTAATGGGAGGGA
AGAGAATAGCGGCAATGTTG55 qRT-PCR引物 β-tubulin-Q-F
β-tubulin-Q-RTTGTCGAAAATGCTGACGAG
CAAGCTTCCGGAGATCAGAG60 内参基因引物 PLT2-2-GFP-F
PLT2-2-GFP-RCATGCCATGGATGAAGAAGTTTGCGGTGCC
GGACTAGTTAACTCAAGTTTGCAGGGTT60 载体构建、亚细胞定位 PLT2-7-GFP-F
PLT2-7-GFP-RCATGCCATGGATGGGGAACTTGACCAAG
GGACTAGTGCTCCTGGGGTAGAAATAGG64 载体构建、亚细胞定位 PLT1-3-Q-F
PLT1-3-Q-RTACCTCCTCCACCACTTCTC
TCTTCAAATCAGAGTCCGC60 qRT-PCR引物 PLT1-4-Q-F
PLT1-4-Q-RATGGTGAGTGCTTGGCTT
TGCCTCCCTTTCCTTGTT50 qRT-PCR引物 1.2.7 过表达载体的构建及亚细胞定位
将PcPLT2-2、PcPLT2-7的CDS序列导入DNAMAN 7.0 软件中,对其进行pCAM-BIA1302的酶切位点预测,本试验选择Nco Ⅰ和Spe I 酶切位点进行后续试验。经扩增、胶回收、连接后将连接产物导入DH5α大肠杆菌感受态细胞,挑取单菌落摇菌后送测,提取重组质粒pCAMBIA1302-目的片段GFP。将带目的片段的重组质粒与pCAMBIA1302植物表达载体分别进行酶切,并回收。采用冻融法转化根癌农杆菌3101,然后采用农杆菌介导法侵染洋葱内表皮PLT进行蛋白亚细胞定位,洋葱37 ℃暗培养3 d后取表皮进行观察(引物序列见表1)。
2. 结果与分析
2.1 多花黄精PLT家族成员鉴定
利用PLT结构域对多花黄精不同组织部位转录组进行搜索筛选,最终鉴定出15个PLT家族成员,依次命名为PcPLT1-1–PcPLT1-4、PcPLT2-1–PcPLT2-11(表2)。PcPLT编码长度159~601 aa、分子量17827.93 Da(PcPLT2-4)~66506.19(PcPLT2-1)、等电点(pI)为6.47~9.64的蛋白质。所有蛋白的不稳定系数都大于 40,表明均为不稳定蛋白;亲水性指数都为负值,表明均为亲水蛋白。信号肽预测结果显示均不含信号肽且15个成员中都不具有跨膜结构。亚细胞定位显示, PcPLT1-4、PcPLT2-1、PcPLT2-2、PcPLT2-5、PcPLT2-8、PcPLT2-9、PcPLT2-10定位在细胞质、细胞核中,PcPLT1-1、PcPLT1-2、PcPLT1-3、PcPLT2-3、PcPLT2-4、PcPLT2-6、PcPLT2-7、PcPLT2-11定位在细胞核中。通过在线软件 SOPMA 分析预测多花黄精 PLT 蛋白序列的二级结构,发现多花黄精 PLT 家族蛋白中大部分成员无规则卷曲结构所占比例高,β-转角结构所占比例最低,所有 α-螺旋结构所占比例都高于延伸链结构。
表 2 多花黄精PLT 家族成员信息Table 2. Information on PcPLT family members基因名称
Gene name基因ID
Gene ID注释到
拟南芥成员
Annotated
Arabidopsis氨基酸数量
Amino
acids/aa分子量
Molecular
weight/Da等电点
pI脂肪系数
Aliphatic
index二级结构
Secondary structure/%亚细胞定位
Subcellular
localizationα-螺旋
Alpha
helixβ-转角
Beta
turn延伸链
Extended
strand无规卷曲
Random
coilPcPLT1-1 >Cluster-11683.110576.p1 AtPLT1 317 35195.85 8.96 54.64 32.18 5.36 10.73 51.74 细胞核 PcPLT1-2 >Cluster-11683.110575.p1 AtPLT1 469 51308.34 6.47 53.41 25.59 5.33 12.15 56.93 细胞核 PcPLT1-3 >Cluster-11683.173509.p1 AtPLT1 468 51236.27 6.47 52.88 25.21 5.13 12.61 57.05 细胞核 PcPLT1-4 >Cluster-11683.233671.p1 AtPLT1 421 47120.43 7.71 58.24 24.47 7.84 13.78 53.92 细胞质、细胞核 PcPLT2-1 >Cluster-11683.222146.p1 AtPLT2 601 66506.19 6.51 59.63 18.47 3.99 12.48 65.06 细胞质、细胞核 PcPLT2-2 >Cluster-11683.100274.p1 AtPLT2 277 30714.02 7.25 61.66 34.30 7.22 15.52 42.96 细胞质、细胞核 PcPLT2-3 >Cluster-11683.77849.p1 AtPLT2 180 20995.90 9.64 70.94 33.89 7.78 26.11 32.22 细胞核 PcPLT2-4 >Cluster-11683.163200.p1 AtPLT2 159 17827.93 6.48 67.61 42.77 5.66 16.35 35.22 细胞核 PcPLT2-5 >Cluster-11683.105355.p1 AtPLT2 312 34191.20 6.94 73.91 26.92 6.73 13.14 53.21 细胞质、细胞核 PcPLT2-6 >Cluster-11683.151415.p2 AtPLT2 199 22672.77 9.11 74.02 39.70 7.54 22.61 30.15 细胞核 PcPLT2-7 >Cluster-11683.204809.p1 AtPLT2 332 36898.44 9.11 61.84 25.00 5.12 12.65 57.23 细胞核 PcPLT2-8 >Cluster-11683.147730.p1 AtPLT2 395 44135.24 9.07 63.01 27.59 6.84 15.95 49.62 细胞质、细胞核 PcPLT2-9 >Cluster-11683.151416.p1 AtPLT2 391 43773.69 8.89 59.64 24.55 7.42 13.81 54.22 细胞质、细胞核 PcPLT2-10 >Cluster-11683.133871.p1 AtPLT2 380 42322.23 9.27 61.89 24.47 7.37 13.68 54.47 细胞质、细胞核 PcPLT2-11 >Cluster-11683.182808.p2 AtPLT2 184 20988.80 9.26 69.46 39.67 7.07 18.48 34.78 细胞核 2.2 多花黄精PLT基因结构、保守结构域及保守基序分析
为进一步揭示 PcPLT的结构和功能特性,采用在线网站 MEME 对多花黄精 PLT 蛋白保守基序进行分析,共鉴定到 10 个保守基序,分别命名为 motif1–motif10(图1),每个蛋白序列上分布的基序种类和顺序有一定差异性,揭示了功能的多样性,其中除PcPLT2-11外,其他成员均含有motif1、motif2、motif3、motif5,具有高度保守性;所有成员均含有2个AP2结构域且无内含子。
2.3 多花黄精 PLT 家族系统进化分析
为分析多花黄精PLT家族的进化关系,参考水稻、拟南芥、石刁柏、油棕和葡萄PLT成员氨基酸序列构建系统进化树(图2)。参考水稻PLT家族的进化分类,将多花黄精PLT家族分为2个大类(Ⅰ和Ⅱ)8个亚家族[Ⅰ(a)、Ⅰ(b)、Ⅰ(c)、Ⅰ(d)、Ⅱ(a)、Ⅱ(b)、Ⅱ(c)、Ⅱ(d)],两个大类各包含6和9个多花黄精成员,进化树显示多花黄精与石刁柏亲缘关系较近,油棕其次,推测是因为多花黄精与石刁柏都属于百合科植物,具有一定相同的序列,而与拟南芥亲缘关系最远,可能是拟南芥是双子叶植物,多花黄精是单子叶植物,因此同源性较其他单子叶植物更远,这有利于进一步探索PcPLT的生物学功能。
2.4 多花黄精PLT基因家族的不同组织部位表达分析
为了探究PLT家族基因不同组织部位的表达差异,使用多花黄精不同组织部位和盐胁迫转录组中PLT家族基因的表达量(FPKM)绘制热图(图3)。不同组织部位转录组结果(图3A)显示,15个PcPLT基因在不同组织部位中均有不同程度的表达,PcPLT2-7、PcPLT2-3、PcPLT2-2、PcPLT2-8在根状茎和根中有较高的表达量,PcPLT2-2和PcPLT2-8在所有部位表达量都较高,PcPLT2-3和PcPLT2-7在根状茎里的表达量最高。根据不同部位转录组筛选到的PcPLT成员在盐胁迫转录组中进行比对,共比对到8个成员(图3B),在500 mmol·L−1 NaCl处理下PcPLT1-3有明显上调表达,PcPLT2-2、PcPLT2-8、PcPLT1-4、PcPLT2-5和PcPLT2-9相比较于CK都有不同程度的下调表达。选取PcPLT2-2、PcPLT2-3、PcPLT2-7进行不同组织部位的表达量分析,选取PcPLT1-3、PcPLT1-4、PcPLT2-7进行不同浓度盐胁迫的qRT-PCR表达分析。结果表明(图4):PcPLT2-2、PcPLT2-3、PcPLT2-7在多花黄精各个组织部位均有表达,以茎秆为CK进行表达量计算,其中PcPLT2-2在根的表达量最高,在叶的表达量最低;PcPLT2-3在根状茎的表达量最高,在叶的表达量最低,根状茎的表达量是叶的45倍;PcPLT2-7在根状茎的表达量最高,在果实的表达量最低,根状茎的表达量是果实的135倍,推测PcPLT基因在根状茎和根中行使重要功能。PcPLT1-3在不同浓度NaCl处理下呈现出上调表达,在500 mmol·L−1 NaCl时表达量最高;PcPLT1-4和PcPLT2-7随着NaCl浓度的提高相对表达量表现出下降的趋势,在500 mmol·L−1 NaCl表达量最低,表明PcPLT1-3正调控盐胁迫,PcPLT1-4、PcPLT2-7负调控盐胁迫。
2.5 多花黄精PcPLT2-2、PcPLT2-7的载体构建及亚细胞定位
为进一步研究PcPLT基因家族亚细胞定位,将构建好的pCAMBIA1302-PcPLT-GFP载体进行洋葱注射,洋葱亚细胞定位结果显示(图5):1302空载GFP充斥在洋葱整个细胞中,1302-PcPLT2-2-GFP定位在细胞质、细胞核,1302-PcPLT2-7-GFP定位在细胞核,均与亚细胞定位预测结果一致。
图 5 PcPLT蛋白在洋葱表皮中的亚细胞定位A:pCAMBIA1302-GFP荧光信号;B:pCAMBIA1302-PcPLT2-2-GFP荧光信号;C:pCAMBIA1302-PcPLT2-7-GFP荧光信号。Figure 5. Subcellular localization of PcPLT in onion epidermal cellsA: pCAMBIA1302-GFP fluorescence signal; B: pCAMBIA1302-PcPLT2-2-GFP fluorescence signal; C: pCAMBIA1302-PcPLT2-7-GFP fluorescence signal.3. 讨论与结论
3.1
多花黄精PcPLT家族成员鉴定与进化特性 PLT 作为AP2/ERF型转录因子,在植物体各个组织中广泛表达,尤其在幼嫩的组织中表达量高。已被证明具有促进根部分生组织生长、芽从头再生、激活根部干细胞建立全能性、茎尖分生组织定位、启动花序分生组织原基等功能,通过与其他转录因子及生长素进行互作影响植物器官的发育与形态建成,对多花黄精遗传改良和增产具有重要意义。
本研究通过生物信息学方法从多花黄精不同组织部位转录组信息鉴定到15个PLT成员,该家族成员均为不稳定亲水蛋白,不含信号肽且不具有跨膜结构,无内含子。系统进化关系表明,PcPLT与石刁柏和油棕的亲缘关系最近,通过在线软件 SOPMA 分析预测多花黄精 PLT 蛋白序列的二级结构,发现多花黄精 PLT 家族蛋白中大部分成员都是无α-螺旋结构和延伸链结构占比较多。蛋白功能的发挥与所在位置关联密切,只有被运输到合适位置才能参与生命活动行使功能,蛋白的亚细胞定位与其蛋白的功能密不可分,因此研究亚细胞定位是研究蛋白功能的一种手段[21]。细胞核是细胞遗传和代谢的调控中心,控制着细胞的分裂和死亡,在细胞的代谢、生长、分化中起着重要作用,PcPLT2-2、PcPLT2-7亚细胞定位结果表明,PcPLT2-2定位在细胞质、细胞核中,PcPLT2-7定位在细胞核中,这与预测结果一致,由此推测 PLT基因在多花黄精生长发育中有重要作用。
3.2 PcPLT可能促进多花黄精的生长发育及参与生长素的调控
Nole-Wilson等[22]研究表明,PLT基因在新生幼嫩组织中表达量更高,而多花黄精不同组织部位qRT-PCR表明,PcPLT2-3、PcPLT2-7基因在根状茎中表达量最高,因此推测PLT基因可能在根状茎新生组织中高表达,更有可能促进多花黄精根状茎的生长发育。生长素是一类重要的植物激素,参与调控植物生长发育等众多生物学过程,尤其在调控根系生长、侧根发育等过程中发挥着不可替代的作用,参与根干细胞龛的形态建成。在拟南芥中,生长素合成基因YUC3和YUC8是PLT2激活靶基因[23],芽中的PLT5-GR在DEX存在的情况下能调控YUC5 表达,PLT与生长素之间相互激活,PLT蛋白通过激活YUC促进生长素积累[24]。不同PLT基因成员与生长素早期反应的共表达模式不同,其中与PLT2和PLT5共表达的生长素早期反应基因明显多于PLT1和PLT7,而与PLT4/BBM共表达的生长素早期反应基因仅为少数。PLT在调节器官发育的细胞分裂和分化阶段中发挥作用的分子机制可能涉及根和地上部的不同靶基因,这可能是由于这些器官系统的不同结构和生长策略。需要对拟南芥和其他物种的根和芽中的AIL/PLT基因调控网络进行更详细的分析,以阐明生长调节剂是如何在不同的发育环境中进化定位。
3.3 PcPLT可能参与多花黄精非生物胁迫的调控
植物在受到非生物胁迫时,会通过升高体内脱落酸(ABA )水平进行响应,这一途径已被充分报道[25]。宋婷婷等[9]通过对水曲柳进行低温(4 ℃)、NaCl、NaHCO3 和 PEG6000 处理,利用qRT-PCR检测PLT基因的表达情况,结果表明均有明显响应,证明PLT基因积极参与植物非生物胁迫耐受性的调控;本研究通过不同浓度NaCl处理表明, PcPLT1-3 的表达量随着NaCl浓度的提高出现上升的趋势, PcPLT1-4 和 PcPLT2-7 则表现出下调,初步说明 PcPLT 调控盐胁迫,因此推测PLT基因在多花黄精中同样存在以上功能,但如何响应或参与非生物胁迫体调控的相关机理尚未明晰。
-
图 3 番茄、辣椒和拟南芥WRKY转录因子系统发生树
氨基酸多序列比对后,用IQ-Tree基于最大似然法构建系统发生树。红色标注指图位克隆到的首个抗青枯病基因RRS1。
Figure 3. Phylogenetic trees of WRKY TFs in tomato, chili pepper and Arabidopsis
IQ-Tree based on maximum likelihood (ML) method after multiple amino acid sequence alignment was used to construct phylogenetic tree. Red indicates 1st bacterial wilt resistance gene RRS1 from map-based cloning.
图 4 番茄WRKY在不同组织部位(a)和病菌侵染(b)下的相对表达
a: 栽培番茄Heinz1706,1. 芽,2. 花,3. 叶,4. 根,5~7. 果(1~3 cm),8. 绿熟期,9. 破色期,10. 破色后10 d;醋栗番茄(S. pimpinellifolium),11. 绿果期,12. 破色期,13. 破色后5 d,14. 叶。b:1. 细菌鞭毛蛋白/对照病原相关分子模式 6 h,2~4. 不同丁香假单胞菌/对照细菌 6 h,5. 荧光假单胞菌/对照细菌 6 h,6. 恶臭假单胞菌/对照细菌 6 h,7. 根癌农杆菌/对照细菌 6 h。
Figure 4. Expression profiles of WRKY TFs in different tissues (a) and with pathogen infection (b) in tomato
a: cultivated tomato Heinz1706, 1. bud, 2. flower, 3. leaf, 4. root, 5–7. fruit (1-3 cm), 8. mature green (MG), 9. breaker (B), 10. breaker+10 (B10); wild S. pimpinellifolium, 11. immature green (IM), 12. breaker (B), 13. breaker+5 (B5), 14. leaf. b: 1. flagellin (flgII-28)/mock pathogen-associated molecular pattern (PAMP) 6 h, 2-4. different Pseudomonas syringae strains (DC3000, DC3000ΔhrcQ-UΔfliC, and DC3000ΔAvrPtoΔAvrPtoB)/bacterial mock 6 h, 5. Pseudomonas fluorescens/bacterial mock 6 h, 6. Pseudomonas putida/bacterial mock 6 h, and 7. Agrobacterium tumefaciens/bacterial mock 6 h.
表 1 番茄WRKY转录因子基因家族
Table 1 WRKY TF gene family in tomato
基因 ID
Gene IDAA MW/kD pI II SL 基因 ID
Gene IDAA MW/kD pI II SL Solyc01g058540.3 324 35.49 4.74 72.92 nucl Solyc05g050060.1 322 36.76 9.23 56.65 nucl Solyc01g079260.4 347 39.15 6.45 40.95 extr Solyc05g050065.1 126 14.62 9.11 58.95 nucl Solyc01g079360.4 240 26.88 5.03 49.87 nucl Solyc05g050300.3 195 22.57 6.06 53.10 cyto Solyc01g089960.3 288 32.57 8.60 46.70 nucl Solyc05g050330.3 244 27.54 6.02 54.36 nucl Solyc01g095100.4 315 34.95 6.50 63.97 nucl Solyc05g050340.4 218 24.58 5.05 46.20 nucl Solyc01g095630.3 336 37.73 5.55 58.66 nucl Solyc05g053380.4 304 34.36 7.75 56.40 nucl Solyc01g104550.3 475 52.98 6.09 50.54 nucl Solyc05g055750.3 459 51.65 8.77 45.35 nucl Solyc02g021680.3 381 42.63 5.37 47.45 nucl Solyc06g008610.3 348 39.05 9.64 53.70 nucl Solyc02g032950.3 504 56.31 8.17 49.81 nucl Solyc06g048870.3 243 26.85 6.15 38.10 nucl Solyc02g067430.3 440 48.40 6.99 60.86 nucl Solyc06g066370.4 549 61.30 6.72 60.62 vacu Solyc02g071130.4 329 37.85 6.55 55.76 nucl Solyc06g068460.3 360 39.74 8.37 43.78 nucl Solyc02g072190.4 304 33.99 6.29 70.76 nucl Solyc06g070990.3 649 71.66 6.07 52.36 nucl Solyc02g080890.3 550 59.66 7.25 45.54 nucl Solyc07g005650.4 513 56.29 6.41 56.26 nucl Solyc02g088340.4 460 50.95 6.20 59.45 nucl Solyc07g047960.3 420 45.96 6.93 69.06 nucl Solyc02g093050.3 326 36.21 9.65 48.66 nucl Solyc07g051840.4 660 72.08 6.32 50.50 nucl Solyc02g094270.2 131 15.47 9.51 35.84 cyto Solyc07g055280.4 275 29.90 5.48 55.13 nucl Solyc03g007380.2 353 40.06 6.09 51.86 nucl Solyc07g056280.3 322 36.65 6.00 61.72 nucl Solyc03g007640.1 359 41.14 6.98 64.88 nucl Solyc07g065260.4 601 65.55 6.58 49.22 nucl Solyc03g082750.1 178 20.73 9.64 38.54 chlo Solyc07g066220.3 739 79.85 6.04 51.97 nucl Solyc03g082810.1 219 25.05 6.55 43.49 nucl Solyc08g006320.4 335 36.44 9.70 37.80 nucl Solyc03g095770.3 273 31.53 6.02 55.04 nucl Solyc08g008280.3 360 40.71 5.51 47.00 nucl Solyc03g104810.3 486 53.90 6.83 66.50 nucl Solyc08g067340.4 279 31.95 9.31 50.58 nucl Solyc03g113120.4 530 58.87 6.83 43.04 nucl Solyc08g067360.3 258 29.35 5.60 51.10 nucl Solyc03g116890.3 350 39.18 8.73 48.24 nucl Solyc08g081610.4 303 34.69 5.49 56.18 nucl Solyc04g050205.1 117 13.73 9.54 44.50 cyto Solyc08g081630.2 231 26.41 8.97 56.81 pero Solyc04g050210.1 366 41.61 5.25 51.49 nucl Solyc08g082110.4 380 43.45 6.82 45.73 chlo Solyc04g051540.3 237 27.14 8.74 57.58 nucl Solyc09g010960.3 290 32.37 5.27 48.26 nucl Solyc04g051690.4 174 19.96 7.05 40.62 nucl Solyc09g014990.4 529 58.67 7.66 56.35 nucl Solyc04g056360.4 413 47.29 5.83 57.52 nucl Solyc09g015770.3 291 33.25 5.51 60.87 nucl Solyc04g072070.3 255 29.59 8.32 37.35 nucl Solyc09g066010.3 331 37.30 9.66 57.28 nucl Solyc04g078550.3 351 38.49 9.63 54.29 nucl Solyc10g005680.2 703 76.13 5.18 47.47 nucl Solyc05g007110.2 392 44.41 8.89 48.16 nucl Solyc10g007970.2 255 28.60 7.07 47.03 nucl Solyc05g012500.3 327 35.43 5.26 57.84 nucl Solyc10g009550.3 290 33.65 5.45 50.93 nucl Solyc05g012770.3 508 55.40 7.65 65.68 nucl Solyc10g011910.4 348 39.17 5.97 50.91 nucl Solyc05g014040.1 290 32.49 9.41 59.32 nucl Solyc10g084380.1 422 46.99 9.49 70.75 nucl Solyc05g015850.4 176 20.31 9.22 40.11 nucl Solyc12g006170.2 549 60.45 6.24 49.52 nucl Solyc05g045710.3 238 26.99 9.00 63.69 cyto Solyc12g011200.3 335 37.80 5.87 63.14 nucl Solyc05g045800.1 253 28.49 9.29 58.72 nucl Solyc12g014610.2 611 65.94 6.23 51.69 nucl Solyc05g045880.1 321 36.33 8.74 50.56 nucl Solyc12g042590.2 252 28.56 8.82 38.23 nucl Solyc05g045927.1 241 27.83 9.25 48.99 nucl Solyc12g056745.1 324 37.80 6.33 38.64 nucl Solyc05g050040.3 278 31.44 9.37 57.09 nucl Solyc12g056750.3 188 22.46 8.78 24.31 nucl Solyc05g050050.1 322 36.75 9.26 58.14 nucl Solyc12g096350.2 338 36.94 9.75 34.56 nucl Solyc05g050057.1 129 14.92 9.71 62.07 nucl 1)AA:氨基酸数;MW:分子量;pI:理论等电点;II:不稳定系数;SL:亚细胞定位。2)nucl:细胞核;cyto:细胞质;chlo:叶绿体;pero:过氧化物酶体;extr:胞外间隙;vacu:液泡。
1) AA: Number of amino acids; MW: molecular weight; pI: Theoretical isoelectric point; II: Instability index; SL: Subcellular localization. 2) nucl: Nucleus; cyto: Cytoplasm; chlo: Chloroplast; pero: Peroxysome; extr: Extracellular; vacu: Vacuole.表 2 番茄WRKY转录因子染色体分布
Table 2 Distribution of WRKY TFs in tomato chromosomes
染色体
Chromosome基因数
Genes基因簇
Gene clusters成簇基因
Genes in clusters单簇基因
Genes per cluster成簇基因占比
Percentage/%串联重复基因
Tandem repeat genes重复基因占比
Percentage/%Chr.1 7 0 0 — 0.0 0 0.0 Chr.2 9 0 0 — 0.0 0 0.0 Chr.3 8 1 2 2 25.0 0 0.0 Chr.4 7 0 0 — 0.0 0 0.0 Chr.5 19 3 12 4 63.2 9 47.4 Chr.6 5 0 0 — 0.0 0 0.0 Chr.7 7 0 0 — 0.0 0 0.0 Chr.8 7 2 4 2 57.1 2 28.6 Chr.9 4 0 0 — 0.0 0 0.0 Chr.10 5 0 0 — 0.0 0 0.0 Chr.11 0 0 0 — — 0 — Chr.12 7 1 2 2 28.6 2 28.6 —:表示不适用,下同。
—: Not applicable. The same below.表 3 番茄与辣椒、拟南芥之间的WRKY同源基因对
Table 3 WRKY homolog between tomato and chili pepper or Arabidopsis
基因1
Gene 1基因2
Gene 2类别
GroupKa Ks Ka/Ks 基因1
Gene 1基因2
Gene 2类别
GroupKa Ks Ka/Ks Solyc01g079260.4 Capana01g002803 IIc 0.114 0.654 0.174 Solyc08g082110.4 Capana01g004471 III 0.268 0.970 0.276 Solyc01g089960.3 Capana01g003441 IIc 0.088 0.355 0.249 Solyc10g005680.2 Capana10g000205 I 0.236 0.283 0.834 Solyc01g095100.4 Capana08g001012 IIe 0.042 0.301 0.140 Solyc10g009550.3 Capana08g001044 III 0.331 1.184 0.279 Solyc01g095630.3 Capana08g001044 III 0.103 0.439 0.235 Solyc10g084380.1 Capana10g001791 I 0.072 0.231 0.314 Solyc01g104550.3 Capana08g001961 IIb 0.052 0.262 0.199 Solyc12g006170.2 Capana00g004057 I 0.088 0.188 0.470 Solyc02g067430.3 Capana02g000918 IIb 0.070 0.316 0.222 Solyc12g011200.3 Capana09g000676 IIc 0.071 0.228 0.310 Solyc02g072190.4 Capana02g001642 IIe 0.066 0.389 0.169 Solyc01g079260.4 AT2G47260.1 IIc 0.488 1.512 0.323 Solyc02g080890.3 Capana02g002230 IIb 0.054 0.373 0.144 Solyc01g089960.3 AT2G44745.1 IIc 0.289 2.165 0.134 Solyc02g088340.4 Capana02g003339 I 0.080 0.222 0.359 Solyc01g095100.4 AT4G01250.1 IIe 0.382 1.931 0.198 Solyc02g093050.3 Capana02g003053 IId 0.116 0.532 0.218 Solyc01g095630.3 AT2G46400.1 III 0.590 — — Solyc02g094270.2 Capana02g003661 IIc 0.099 0.278 0.356 Solyc02g021680.3 AT2G34830.2 IIe 0.446 — — Solyc03g007380.2 Capana03g002072 III 0.103 0.486 0.213 Solyc02g071130.4 AT4G18170.1 IIc 0.446 1.921 0.232 Solyc03g007640.1 Capana03g002134 IIe 0.388 0.708 0.547 Solyc02g071130.4 AT5G46350.1 IIc 0.424 2.659 0.159 Solyc03g082810.1 Capana03g001962 IIe 0.142 0.308 0.460 Solyc02g071130.4 AT1G29860.1 IIc 0.366 2.306 0.159 Solyc03g095770.3 Capana03g002635 III 0.148 0.309 0.480 Solyc02g088340.4 AT1G13960.1 I 0.382 1.643 0.233 Solyc03g113120.4 Capana03g001099 IIb 0.083 0.399 0.209 Solyc02g088340.4 AT3G01080.2 I 0.512 — — Solyc03g116890.3 Capana03g000473 IIa 0.102 0.455 0.224 Solyc02g094270.2 AT3G01970.1 IIc 0.384 — — Solyc04g051540.3 Capana12g001851 IIc 0.042 0.194 0.215 Solyc03g116890.3 AT1G80840.1 IIa 0.340 2.918 0.117 Solyc04g051690.4 Capana12g001826 IIc 0.080 0.300 0.266 Solyc05g007110.2 AT1G69810.1 IIb 0.491 — — Solyc04g056360.4 Capana05g002502 I 0.394 1.071 0.368 Solyc05g012500.3 AT1G69310.1 IIc 0.376 3.568 0.106 Solyc04g078550.3 Capana04g000568 IId 0.066 0.660 0.099 Solyc05g012770.3 AT1G13960.1 I 0.285 2.483 0.115 Solyc05g012500.3 Capana11g001905 IIc 0.132 0.357 0.369 Solyc05g012770.3 AT2G03340.1 I 0.290 2.177 0.133 Solyc05g012770.3 Capana11g001882 I 0.045 0.346 0.130 Solyc06g048870.3 AT2G40740.3 NG 0.389 2.044 0.190 Solyc05g055750.3 Capana05g002502 I 0.118 0.324 0.363 Solyc06g070990.3 AT1G18860.1 IIb 0.464 2.330 0.199 Solyc06g008610.3 Capana06g003072 IId 0.042 0.251 0.168 Solyc06g070990.3 AT5G15130.1 IIb 0.524 3.356 0.156 Solyc06g066370.4 Capana06g001506 I 0.032 0.261 0.123 Solyc07g047960.3 AT2G04880.1 I 0.513 2.372 0.216 Solyc06g068460.3 Capana06g001110 IIa 0.053 0.282 0.187 Solyc07g065260.4 AT4G26640.2 I 0.370 2.286 0.162 Solyc06g070990.3 Capana06g001008 IIb 0.058 0.396 0.146 Solyc08g006320.4 AT2G24570.1 IId 0.343 1.396 0.245 Solyc07g005650.4 Capana07g000181 I 0.131 0.307 0.427 Solyc08g006320.4 AT4G31550.1 IId 0.401 2.820 0.142 Solyc07g047960.3 Capana07g001256 I 0.089 0.162 0.553 Solyc08g081610.4 AT4G23550.1 IIe 0.552 — — Solyc07g051840.4 Capana07g001387 IIb 0.061 0.259 0.235 Solyc08g082110.4 AT4G11070.1 III 0.439 2.534 0.173 Solyc07g055280.4 Capana07g001809 IIe 0.087 0.445 0.195 Solyc08g082110.4 AT4G23810.1 III 0.436 — — Solyc07g056280.3 Capana07g001968 IIc 0.051 0.251 0.202 Solyc09g014990.4 AT2G38470.1 I 0.368 — — Solyc07g065260.4 Capana07g002350 I 0.054 0.176 0.308 Solyc10g005680.2 AT4G26440.1 I 0.524 1.613 0.325 Solyc08g008280.3 Capana01g004471 III 0.131 0.401 0.327 Solyc10g084380.1 AT2G37260.1 I 0.496 2.115 0.235 Solyc08g081610.4 Capana01g000167 IIe 0.076 0.228 0.334 Solyc12g014610.2 AT4G26640.2 I 0.394 2.070 0.190 Solyc08g081630.2 Capana01g000165 IIc 0.187 0.522 0.358 表 4 番茄青枯菌诱导转录差异表达WRKY
Table 4 Transcriptionally DE WRKY TFs with RsI in tomato
基因ID
Gene ID类别
GroupRC_vs_RT
FDRRC_vs_RT
log2FCSC_vs_ST
FDRSC_vs_ST
log2FCGO/KEGG注释
GO/KEGG annotationSolyc01g089960.3 IIc 0.00 −0.99 0.90 0.05 GO:0003700、GO:0006355、GO:0043565 Solyc01g104550.3 IIb 0.03 0.65 0.62 0.22 GO:0003700、GO:0006355、GO:0043565、GO:0044212 Solyc02g032950.3 IIb 0.01 1.05 0.19 0.49 GO:0003700、GO:0006355、GO:0043565 Solyc03g095770.3 III 0.80 0.12 0.02 0.78 GO:0003700、GO:0006355、GO:0043565 Solyc04g051540.3 IIc 0.00 −1.08 0.54 −0.24 GO:0003700、GO:0006355、GO:0043565 Solyc05g050340.4 III 0.46 0.26 0.05 0.71 GO:0003700、GO:0006355、GO:0043565 Solyc05g055750.3 I 0.91 0.06 0.02 −0.67 GO:0003700、GO:0006355、GO:0043565 Solyc07g051840.4 IIb 0.00 1.17 0.00 0.76 GO:0003700、GO:0006355、GO:0043565 Solyc09g014990.4 I 0.01 1.06 0.01 0.88 GO:0003700、GO:0006355、GO:0043565、ko04626 Solyc09g015770.3 III 0.20 0.61 0.00 1.44 GO:0003700、GO:0006355、GO:0043565 Solyc09g066010.3 IId 0.37 −0.24 0.00 −0.66 GO:0003700、GO:0006355、GO:0043565 Solyc10g009550.3 III 0.00 1.18 0.00 0.90 GO:0003700、GO:0006355、GO:0043565 R和S分别代表抗病和感病番茄自交系,C和T分别代表对照和青枯菌侵染。GO:0003700表示转录因子活性和序列特异性DNA结合,GO:0006355表示转录调控,GO:0043565表示序列特异性DNA结合,GO:0044212表示转录调控区DNA结合,ko04626表示植物-病原互作。
R and S: resistant and susceptible tomato inbred lines, respectively; C: control; T: RsI. GO:0003700: transcription factor activity and sequence-specific DNA binding; GO:0006355: regulation of transcription; GO:0043565: sequence-specific DNA binding; GO:0044212: transcription regulatory region DNA binding; ko04626: plant-pathogen interaction. -
[1] 刘强, 张贵友, 陈受宜. 植物转录因子的结构与调控作用 [J]. 科学通报, 2000, 45(14):1465−1474. DOI: 10.3321/j.issn:0023-074X.2000.14.002 LIU Q, ZHANG G Y, CHEN S Y. Structure and regulation of plant transcription factors [J]. Chinese Science Bulletin, 2000, 45(14): 1465−1474.(in Chinese) DOI: 10.3321/j.issn:0023-074X.2000.14.002
[2] RUSHTON P J, SOMSSICH I E, RINGLER P, et al. WRKY transcription factors [J]. Trends in Plant Science, 2010, 15(5): 247−258. DOI: 10.1016/j.tplants.2010.02.006
[3] EULGEM T, RUSHTON P J, ROBATZEK S, et al. The WRKY superfamily of plant transcription factors [J]. Trends in Plant Science, 2000, 5(5): 199−206. DOI: 10.1016/S1360-1385(00)01600-9
[4] AGARWAL P, REDDY M P, CHIKARA J. WRKY: Its structure, evolutionary relationship, DNA-binding selectivity, role in stress tolerance and development of plants [J]. Molecular Biology Reports, 2011, 38(6): 3883−3896. DOI: 10.1007/s11033-010-0504-5
[5] 赵楠楠, 刘立峰. 植物WRKY转录因子及其生物学功能 [J]. 分子植物育种, 2019, 17(21):7040−7046. DOI: 10.13271/j.mpb.017.007040 ZHAO N N, LIU L F. WRKY transcription factors and their biological functions in plants [J]. Molecular Plant Breeding, 2019, 17(21): 7040−7046.(in Chinese) DOI: 10.13271/j.mpb.017.007040
[6] XIE Z, ZHANG Z L, ZOU X L, et al. Annotations and functional analyses of the rice WRKY gene superfamily reveal positive and negative regulators of abscisic acid signaling in aleurone cells [J]. Plant Physiology, 2005, 137(1): 176−189. DOI: 10.1104/pp.104.054312
[7] BAKSHI M, OELMÜLLER R. WRKY transcription factors: Jack of many trades in plants [J]. Plant Signaling & Behavior, 2014, 9(2): e27700.
[8] WU K L, GUO Z J, WANG H H, et al. The WRKY family of transcription factors in rice and Arabidopsis and their origins [J]. DNA Research, 2005, 12(1): 9−26. DOI: 10.1093/dnares/12.1.9
[9] ZHANG Y J, WANG L J. The WRKY transcription factor superfamily: Its origin in eukaryotes and expansion in plants [J]. BMC Evolutionary Biology, 2005, 5(1): 1. DOI: 10.1186/1471-2148-5-1
[10] CHEN L G, SONG Y, LI S J, et al. The role of WRKY transcription factors in plant abiotic stresses [J]. Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms, 2012, 1819(2): 120−128. DOI: 10.1016/j.bbagrm.2011.09.002
[11] ROSS C A, LIU Y, SHEN Q J. The WRKY gene family in rice (Oryza sativa) [J]. Journal of Integrative Plant Biology, 2007, 49(6): 827−842. DOI: 10.1111/j.1744-7909.2007.00504.x
[12] WEI K F, CHEN J, CHEN Y F, et al. Molecular phylogenetic and expression analysis of the complete WRKY transcription factor family in maize [J]. DNA Research, 2012, 19(2): 153−164. DOI: 10.1093/dnares/dsr048
[13] BENCKE-MALATO M, CABREIRA C, WIEBKE-STROHM B, et al. Genome-wide annotation of the soybean WRKY family and functional characterization of genes involved in response to Phakopsora pachyrhizi infection [J]. BMC Plant Biology, 2014, 14: 236. DOI: 10.1186/s12870-014-0236-0
[14] ZHANG C, WANG D D, YANG C H, et al. Genome-wide identification of the potato WRKY transcription factor family [J]. PLoS One, 2017, 12(7): e0181573. DOI: 10.1371/journal.pone.0181573
[15] HUANG S X, GAO Y F, LIU J K, et al. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum [J]. Molecular Genetics and Genomics, 2012, 287(6): 495−513. DOI: 10.1007/s00438-012-0696-6
[16] YANG Y, LIU J, ZHOU X H, et al. Identification of WRKY gene family and characterization of cold stress-responsive WRKY genes in eggplant [J]. PeerJ, 2020, 8: e8777. DOI: 10.7717/peerj.8777
[17] SCIENCE F I P. Retraction: Genome-wide identification and expression analysis of WRKY gene family in Capsicum annuum L [J]. Frontiers in Plant Science, 2016, 7: 1727.
[18] LING J, JIANG W J, ZHANG Y, et al. Genome-wide analysis of WRKY gene family in Cucumis sativus [J]. BMC Genomics, 2011, 12: 471. DOI: 10.1186/1471-2164-12-471
[19] PHUKAN U J, JEENA G S, SHUKLA R K. WRKY transcription factors: Molecular regulation and stress responses in plants [J]. Frontiers in Plant Science, 2016, 7: 760.
[20] 史建磊, 熊自立, 苏世闻, 等. 基于RNA-seq的番茄青枯病响应基因鉴定与表达分析 [J]. 华北农学报, 2022, 37(2):171−182. DOI: 10.7668/hbnxb.20192621 SHI J L, XIONG Z L, SU S W, et al. Identification and expression analysis of bacterial wilt response genes based on RNA-seq in tomato [J]. Acta Agriculturae Boreali-Sinica, 2022, 37(2): 171−182.(in Chinese) DOI: 10.7668/hbnxb.20192621
[21] CHEN C J, CHEN H, ZHANG Y, et al. TBtools: An integrative toolkit developed for interactive analyses of big biological data [J]. Molecular Plant, 2020, 13(8): 1194−1202. DOI: 10.1016/j.molp.2020.06.009
[22] ASHBURNER M, BALL C A, BLAKE J A, et al. Gene Ontology: Tool for the unification of biology [J]. Nature Genetics, 2000, 25(1): 25−29. DOI: 10.1038/75556
[23] KANEHISA M, GOTO S, KAWASHIMA S, et al. The KEGG resource for deciphering the genome [J]. Nucleic Acids Research, 2004, 32(Suppl_1): D277−D280.
[24] LEE S W, HAN S W, SRIRIYANUM M, et al. A type I–secreted, sulfated peptide triggers XA21-mediated innate immunity [J]. Science, 2009, 326(5954): 850−853. DOI: 10.1126/science.1173438
[25] 张红, 姜景彬, 许向阳, 等. 番茄WRKY基因家族的生物信息学分析 [J]. 分子植物育种, 2016, 14(8):1965−1976. DOI: 10.13271/j.mpb.014.001965 ZHANG H, JIANG J B, XU X Y, et al. Bioinformatics analysis of WRKY gene family in tomato [J]. Molecular Plant Breeding, 2016, 14(8): 1965−1976.(in Chinese) DOI: 10.13271/j.mpb.014.001965
[26] SONG H, SUN W H, YANG G F, et al. WRKY transcription factors in legumes [J]. BMC Plant Biology, 2018, 18(1): 243. DOI: 10.1186/s12870-018-1467-2
[27] MOHANTA T K, PARK Y H, BAE H H. Novel genomic and evolutionary insight of WRKY transcription factors in plant lineage [J]. Scientific Reports, 2016, 6(1): 1−22. DOI: 10.1038/s41598-016-0001-8
[28] CHEN F, HU Y, VANNOZZI A, et al. The WRKY transcription factor family in model plants and crops [J]. Critical Reviews in Plant Sciences, 2017, 36(5/6): 311−335.
[29] 刁卫平, 王述彬, 刘金兵, 等. 辣椒全基因组WRKY转录因子的分析 [J]. 园艺学报, 2015, 42(11):2183−2196. DIAO W P, WANG S B, LIU J B, et al. Genome-wide analysis of the WRKY transcription factor family in pepper [J]. Acta Horticulturae Sinica, 2015, 42(11): 2183−2196.(in Chinese)
[30] LIU J L, LIU X L, DAI L Y, et al. Recent progress in elucidating the structure, function and evolution of disease resistance genes in plants [J]. Journal of Genetics and Genomics, 2007, 34(9): 765−776. DOI: 10.1016/S1673-8527(07)60087-3
[31] ÜLKER B, SOMSSICH I E. WRKY transcription factors: From DNA binding towards biological function [J]. Current Opinion in Plant Biology, 2004, 7(5): 491−498. DOI: 10.1016/j.pbi.2004.07.012
[32] RINERSON C I, RABARA R C, TRIPATHI P, et al. The evolution of WRKY transcription factors [J]. BMC Plant Biology, 2015, 15: 66. DOI: 10.1186/s12870-015-0456-y
[33] DESLANDES L, OLIVIER J, THEULIERES F, et al. Resistance to Ralstonia solanacearum in Arabidopsis thaliana is conferred by the recessive RRS1-R gene, a member of a novel family of resistance genes [J]. Proceedings of the National Academy of Sciences of the United States of America, 2002, 99(4): 2404−2409. DOI: 10.1073/pnas.032485099
[34] MUKHTAR M S, DESLANDES L, AURIAC M C, et al. The Arabidopsis transcription factor WRKY27 influences wilt disease symptom development caused by Ralstonia solanacearum [J]. The Plant Journal, 2008, 56(6): 935−947. DOI: 10.1111/j.1365-313X.2008.03651.x
[35] DANG F F, WANG Y N, YU L, et al. CaWRKY40, a WRKY protein of pepper, plays an important role in the regulation of tolerance to heat stress and resistance to Ralstonia solanacearum infection [J]. Plant, Cell & Environment, 2013, 36(4): 757−774.
[36] CAI H Y, YANG S, YAN Y, et al. CaWRKY6 transcriptionally activates CaWRKY40, regulates Ralstonia solanacearum resistance, and confers high-temperature and high-humidity tolerance in pepper [J]. Journal of Experimental Botany, 2015, 66(11): 3163−3174. DOI: 10.1093/jxb/erv125
[37] HUSSAIN A, LI X, WENG Y H, et al. CaWRKY22 acts as a positive regulator in pepper response to RalstoniaSolanacearum by constituting networks with CaWRKY6, CaWRKY27, CaWRKY40, and CaWRKY58 [J]. International Journal of Molecular Sciences, 2018, 19(5): 1426. DOI: 10.3390/ijms19051426
[38] DANG F F, WANG Y N, SHE J J, et al. Overexpression of CaWRKY27, a subgroup IIe WRKY transcription factor of Capsicum annuum, positively regulates tobacco resistance to Ralstonia solanacearum infection [J]. Physiologia Plantarum, 2014, 150(3): 397−411. DOI: 10.1111/ppl.12093
[39] YANG S, ZHANG Y W, CAI W W, et al. CaWRKY28 Cys249 is required for interaction with CaWRKY40 in the regulation of pepper immunity to Ralstonia solanacearum [J]. Molecular Plant-Microbe Interactions:MPMI, 2021, 34(7): 733−745. DOI: 10.1094/MPMI-12-20-0361-R
[40] HUSSAIN A, KHAN M I, ALBAQAMI M, et al. CaWRKY30 positively regulates pepper immunity by targeting CaWRKY40 against Ralstonia solanacearum inoculation through modulating defense-related genes [J]. International Journal of Molecular Sciences, 2021, 22(21): 12091. DOI: 10.3390/ijms222112091
[41] IFNAN KHAN M, ZHANG Y W, LIU Z Q, et al. CaWRKY40b in pepper acts as a negative regulator in response to Ralstonia solanacearum by directly modulating defense genes including CaWRKY40 [J]. International Journal of Molecular Sciences, 2018, 19(5): 1403. DOI: 10.3390/ijms19051403
[42] WANG Y N, DANG F F, LIU Z Q, et al. CaWRKY58, encoding a group WRKY transcription factor of Capsicum annuum, negatively regulates resistance to Ralstonia solanacearum infection [J]. Molecular Plant Physiology, 2013, 14(2): 131−144.
[43] 谢政文, 王连军, 陈锦洋, 等. 植物WRKY转录因子及其生物学功能研究进展 [J]. 中国农业科技导报, 2016, 18(3):46−54. DOI: 10.13304/j.nykjdb.2015.605 XIE Z W, WANG L J, CHEN J Y, et al. Studies on WRKY transcription factors and their biological functions in plants [J]. Journal of Agricultural Science and Technology, 2016, 18(3): 46−54.(in Chinese) DOI: 10.13304/j.nykjdb.2015.605
[44] CHENG Y, ZHOU Y, YANG Y, et al. Structural and functional analysis of VQ motif-containing proteins in Arabidopsis as interacting proteins of WRKY transcription factors [J]. Plant Physiology, 2012, 159(2): 810−825. DOI: 10.1104/pp.112.196816
[45] PARK C Y, LEE J H, YOO J H, et al. WRKY group IId transcription factors interact with calmodulin [J]. FEBS Letters, 2005, 579(6): 1545−1550. DOI: 10.1016/j.febslet.2005.01.057
[46] 黄幸, 丁峰, 彭宏祥, 等. 植物WRKY转录因子家族研究进展 [J]. 生物技术通报, 2019, 35(12):129−143. DOI: 10.13560/j.cnki.biotech.bull.1985.2019-0626 HUANG X, DING F, PENG H X, et al. Research progress on family of plant WRKY transcription factors [J]. Biotechnology Bulletin, 2019, 35(12): 129−143.(in Chinese) DOI: 10.13560/j.cnki.biotech.bull.1985.2019-0626