小麦4个WAKs/WAKLs基因的分子特征及其响应条锈菌侵染的表达
鉴定4个小麦WAKs/WAKLs基因，研究其响应条锈菌(Puccinia striiformis f. sp. tritici)侵染的表达模式。
从小麦条锈病抗病材料云337和感病材料云402转录组差异表达的基因中筛选出4个WAKs/WAKLs基因，对其基因结构、蛋白质理化性质、保守结构域和进化进行分析，通过RT-qPCR验证4个基因在小麦接种条锈菌后的表达模式。
所鉴定的TraesCS1B02G004100、TraesCS6B02G459100、TraesCSU02G079300和TraesCS3A02G122300基因长度分别为3 194、5 117、3 208和2 934 bp，分别编码714、960、737和652个氨基酸，均包括胞外GUB结构域和胞内激酶结构域，均位于细胞膜。系统进化分析表明：TraesCS6B02G459100与水稻Os05g04460、TraesCSU02G079300与水稻Os02g02120、TraesCS3A02G122300与拟南芥AtWAKL15、TraesCS1B02G004100与大麦HORVU.MOREX.r3.5HG0503880的亲缘关系分别最近。激酶部位氨基酸序列分析表明：TraesCS1B02G004100和TraesCSU02G079300为与植物先天免疫机制相关的非RD激酶，TraesCS6B02G459100和TraesCS3A02G122300为RD激酶。基因表达分析表明：4个WAKs/WAKLs基因在小麦接种条锈菌后均差异表达，其中TraesCS1B02G004100在抗病材料中的表达被抑制、在感病材料中不表达。
TraesCS3A02G122300和TraesCSU02G079300在小麦与条锈菌亲和及非亲和互作体系中均差异表达；TraesCS6B02G459100和TraesCS1B02G004100仅在非亲和互作中特异表达。非RD激酶TraesCS1B02G004100表达受到抑制，其调控抗病性的功能和机制需进一步研究。
-
关键词:
- 小麦 /
- WAKs/WAKLs基因 /
- 表达分析 /
- 条锈菌
Molecular Characteristics of Four WAKs/WAKLs Genes in Wheat and Their Expression in Response to Puccinia striiformis f. sp. tritici Infection
To identify four wheat WAKs/WAKLs genes and analyze their expression patterns responding to Puccinia striiformis f.sp. tritici (Pst) infection.
Four WAKs/WAKLs genes were screened from transcriptome differentially expressed genes of wheat, which inoculated Pst in resistant variety Yun337 and susceptible variety Yun402. The gene structure, physicochemical property of protein, conserved domains and phylogenetic relationship were analyzed, and the expression patterns of four genes in wheat inoculated Pst were verified by RT-qPCR.
The full lengths of TraesCS1B02G004100, TraesCS6B02G459100, TraesCSU02G079300, and TraesCS3A02G122300 were 3194, 5117, 3208, and 2934 bp, respectively; they encoded 714, 960, 737 and 652 amino acids, respectively; they all included extracellular GUB domain and intracellular Ser/Thr kinase domain; and they were all located on the cell membrane. Phylogenetic analysis showed that TraesCS6B02G459100 with Os05g04460, TraesCSU02G079300 with Os02g02120, TraesCS3A02G122300 with AtWAKL15, TraesCS1B02G004100 with HORVU.MOREX.r3.5HG0503880, had a close relationship, respectively. The kinase site amino acid analysis showed that TraesCS1B02G004100 and TraesCSU02G079300 were non-RD kinases related to plant innate immune mechanism, while TraesCS6B02G459100 and TraesCS3A02G122300 were RD kinases. Gene expression analysis showed that four WAKs/WAKLs genes in wheat were expressed differently upon infection of Pst, and TraesCS1B02G004100 was down-regulated in resistant variety but scarcely expressed in susceptible variety.
TraesCS3A02G122300 and TraesCSU02G079300 are expressed differently in both compatible and incompatible interactions between wheat and Pst; TraesCS6B02G459100 and TraesCS1B02G004100 are only respectively expressed in incompatible interaction. The expression of non-RD receptor kinase gene TraesCS1B02G004100 is inhibited, and its function and mechanism in disease resistance deserve further investigation.
-
类受体激酶(receptor-like kinase，RLK)是植物受体最主要的组成成分，在细胞间的跨膜转导信号中起重要作用，控制植物的繁殖、生长、发育及对各种环境的适应[1]。RLK家族一般包含胞外结构域(extracellular domain，ED)、跨膜区(transmembrane region，TM)和胞内激酶结构域(Ser/Thr kinase domain，STK)[2]。在植物病原菌互作的信号转导通路中，RLK作为模式识别受体 (pattern-recognition receptors，PRRs)可以识别病原菌分子模式 (pathogen-associated molecular patterns，PAMPs)，从而在细胞膜上响应植物的天然免疫系统(pattern-triggered immunity，PTI)，激活植物的抗病能力[3-5]。当外界信号PAMPs分子被受体类激酶的ED识别后，通过STK将信号传导给下游的信号通路，引发一系列免疫反应，包括ROS爆发、MAPK激活、胼胝质沉淀、钙离子流入、水杨酸积累等[3-5]。如油菜素内酯信号转导通路中，共受体BAK1和受体FLS2共同识别病原菌的鞭毛蛋白，BRI1受体类激酶和BAK1相互作用并经交互磷酸化形成1个活跃的受体复合物，从而激活油菜素内酯抵御病原菌。越来越多的RLKs已被证明是PRRs[3-5]。
细胞壁受体激酶(wall-associated kinases，WAK)及其类基因WAKLs (WAK like genes，WAKLs)属于RLK家族的亚家族[6]，合称为WAK/WAKL家族。WAKs和WAKLs的结构非常相似，主要是胞外表皮生长因子结构域(epidermal growth factor domain，EGF)、跨膜结构域和胞内STK，但部分WAKLs缺少一些结构域，如WAKL7缺少跨膜域和激酶结构域[7]。在植物防御反应中，拥有胞外结构域的WAKs/WAKLs基因可以作为细胞外界信号接收员，而胞内激酶部位参与胞质信号级联反应[7]，并且在生长发育和对非生物压力的耐受性中起作用。WAKs/WAKLs基因在植物防御病原体方面发挥着关键作用，如：AtWAK1编码WAKs蛋白，与病原菌互作，在拟南芥的防御反应中起重要作用[8]；从玉米(Zea mays L.)中克隆的主效数量性状基因ZmWAK可抗玉米丝黑穗病[9]；在棉花(Gossypium hirsutum L.)中鉴定到的GhWAK7A参与调节对镰刀菌和黄萎病的反应[10]；在山新杨(Populus davidiana×P. bolleana)中过表达克隆到PdPapWAK基因后，其抵抗非生物和生物胁迫的能力进一步增强[11]。
植物WAK/WAKL家族是响应病原菌信号的重要成员。长久以来，基因对基因假说提供了植物抗病模式PTI与ETI系统、植物PRRs与抗病基因(resistance gene，R)的明确界限，但这一界限面临着大量研究的挑战。WAK/WAKL家族作为信号转导途径中的重要受体，可以激活植物的PTI免疫反应模式，检测到病原菌的效应蛋白，可能原因之一是其在植物R基因与病原菌无毒基因互作过程中发挥“监测”作用[12-13]。在小麦对禾谷丝核菌的防御反应中，鉴定到WAK7的转录丰度升高[14]，积极响应了小麦禾谷丝核菌的防御反应。小麦基因组中定位到1个WAKL4/stb6基因，该基因识别病原菌效应物Avr Stb6，且绕过过敏性反应赋予小麦对枯叶病的抗性，符合基因对基因假说[15]。当过表达OsWAK91时，水稻激活ETI系统，产生H2O2反应并增强相关R基因的表达[16]。SlWAK1通过番茄中的FLS2/FLS3复合物参与调节PRRs介导的免疫应答[17]。总之，WAK/WAKL家族在PTI和ETI系统中的调控机制很复杂，具体情况尚不清楚。
本研究从课题组前期测定的小麦响应条锈病侵染的转录组数据中筛选出4个在抗感材料中差异表达的WAKs/WAKLs基因，对其进行基因结构、蛋白基本理化性质、保守结构域、系统进化等生物信息学分析，结合实时荧光定量 PCR (RT-qPCR)在抗病和感病材料中分析其响应条锈菌的表达模式，为进一步研究小麦WAKs/WAKLs基因响应条锈菌的侵染作用机制及功能提供参考。
1. 材料与方法
1.1 试验材料
供试小麦品种为抗条锈病品种云337和感条锈病品种云402；小麦条锈菌(Puccinia striiformis f.sp.tritici)为CYR32。将小麦播种在人工培养室中，温度控制在昼16~18 ℃、夜12~15 ℃，每天光照16 h；当小麦生长到二叶期时，将CYR32与滑石粉按体积比1∶20混合后进行撒粉接种；分别采集条锈菌侵染后0、24、48、72、120和168 h的叶片样本，于−80 ℃保存备用。
1.2 转录组测序及差异基因表达分析
转录组测序由北京百迈客生物科技有限公司完成。文库质检合格的测序数据采用DESeq2进行差异基因分析，以差异倍数(fold change，FC)≥1.5且错误发现率(false discovery rate，FDR)<0.01作为差异基因筛选标准，从中挖掘在抗感材料中差异表达的小麦WAKs/WAKLs基因，并对其进行COG分类、GO功能富集、KEGG注释和KEGG通路富集分析。
1.3 小麦WAK/WAKL家族基因鉴定
下载拟南芥、水稻、大麦和小麦基因组数据库文件，采用TBtools v1.108 从中比对出差异表达的小麦WAKs/WAKLs基因的同源基因组和氨基酸序列，并筛选同源性较高的基因；利用NCBI的Batch CD-Search在线软件(https//www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi)对筛选得到的候选基因编码蛋白进行蛋白结构域分析。
1.4 小麦WAK/WAKL家族成员生物信息学分析
利用GSDS在线软件(http//gsds.cbi.pku.edu.cn/)分析和绘制小麦WAKs/WAKLs基因结构；利用ExPASyProtParam Tool (http//web.expasy.org/protparam)预测蛋白质的等电点和分子质量；利用ExPASy (http//web.expasy.org/protscale)预测蛋白质的亲水性和疏水性；利用MEME 在线软件(https://meme-suite.org/meme/tools/meme)预测蛋白基序。为了研究小麦WAKs/WAKLs蛋白与其他物种间的系统进化关系，对小麦、大麦、拟南芥和水稻4个物种中同源性较高的WAKs/WAKLs蛋白全长序列进行比对，采用DNAMAN V6比对蛋白质多序列，采用MEGA 7.0的邻接法构建系统进化树，采用 TBtools 绘制系统进化树与基序预测图。
1.5 小麦WAKs/WAKLs基因的表达模式分析
根据从转录组筛选出的小麦WAKs/WAKLs基因序列，使用Beacon Designer 8设计RT-qPCR引物序列(表1)，引物由北京擎科生物股份有限公司合成。采用RNA Easy Fast 植物组织 RNA快速提取试剂盒(天根生化科技有限公司)提取小麦叶片的总RNA，再采用Aid lab公司反转录试剂盒(TURE script 1st Stand cDNA SYNTHESIS Kit)进行反转录。以小麦GAPDH基因为内参，使用2×SYBR® Green预混液(中国DF)试剂盒在Analytikjena-qTOWER2.2型荧光定量PCR仪(德国)上进行RT-qPCR反应。反应体系总体积10.0 μL，包括cDNA 模板1.0 μL，SYBR Green 1 Master 5.0 μL，上、下游引物各 0.5 μL，ddH2O 3.0 μL。反应程序为：95 ℃ 3 min；95 ℃ 10 s，60 ℃ 30 s，39 个循环 ，溶解曲线(60~95 ℃，每个循环增加1 ℃，每次循环时间为4 s)设置2个生物学重复，并使用2-ΔΔCt[18]计算各个样品及各组中的基因相对表达量。采用Excel 2010 处理数据；采用GraphPad Prism 8.4.2 作图。
表 1 供试引物Table 1. Primers for test引物名称
primer names引物序列 (5′→3′)
primer sequencesGAPDH-F GCCTTGCTCCTCTTGCTA GAPDH-R TCCACCTCTCCAGTCCTT TraesCS1B02G004100-F CAATCAGCGGACGGTCTA TraesCS1B02G004100-R ATGTCGCCGATGACTTGT TraesCS6B02G459100-F TTCATCGGAATGTGGTCAAG TraesCS6B02G459100-R GCTGAGTGGAGATAAGAAAGG TraesCSU02G079300-F GAGATGTCGCTGAGAGGTT TraesCSU02G079300-R CCATGCTGTCCAGTGCTA TraesCS3A02G122300-F AGCTCGACACAATGAAGG TraesCS3A02G122300-R CCCAGCCTCAATGTTGAT 2. 结果与分析
2.1 4个WAKs/WAKLs基因的结构
经转录组筛选出4个小麦WAKs/WAKLs基因，分别为TraesCS1B02G004100、TraesCS6B02G459100、TraesCSU02G079300和TraesCS3A02G122300基因(图1)。其中，基因组序列最长的是TraesCS6B02G459100，长度为5 117 kb，包含4个外显子和3个内含子；序列最短的是TraesCS3A02G122300，长度为2 934 kb，包含3个外显子和2个内含子；TraesCSU02G079300和TraesCS1B02G004100的长度相近，分别为3208和3194 bp，均包含3个外显子和2个内含子。
2.2 WAKs/WAKLs蛋白的理化性质和亚细胞定位
由表2可知：TraesCS1B02G004100、TraesCS6B02G459100、TraesCSU02G079300和TraesCS3A02G122300编码的氨基酸序列长度分别为 714、960、737和652 aa；蛋白分子质量分别为78.68、106.29、81.86和69.67 ku；理论等电点分别为8.50、5.98、5.65和6.25；脂肪族氨基酸指数分别为89.59、85.33、81.38和90.49；4个基因的蛋白亲疏水性值均为负值，为亲水性蛋白。亚细胞定位预测显示：4个基因均定位于细胞膜上，为膜蛋白。
表 2 4个WAKs/WAKLs 基因理化性质Table 2. Physicochemical properties of four WAKs/WAKLs genes基因
gene蛋白质序列号
accession number of protein多肽长度/aa
peptide length蛋白分子质量/ku
protein molecular weight等电点
pI氨基酸指数
amino acid
index蛋白亲水性
protein
hydrophilicityTraesCS1B02G004100 XP_044361219.1 714 78.68 8.50 89.59 −0.119 TraesCS6B02G459100 XP_044410226.1 960 106.29 5.98 85.33 −0.120 TraesCSU02G079300 XP_044407301.1 737 81.86 5.65 81.38 −0.214 TraesCS3A02G122300 XP_044339162.1 652 69.67 6.25 90.49 −0.024 2.3 WAKs/WAKLs的系统进化以及基序和功能结构域
小麦、大麦、拟南芥和水稻中同源性较高的43个WAKs/WAKLs蛋白质序列按照亲缘关系被分为4个分支，TraesCS6B02G459100和来自水稻、大麦、拟南芥的成员聚为一类，处于分支Ⅰ，且与水稻Os05g04460亲缘关系最近；TraesCSU02G079300和TraesCS3A02G122300处于分支Ⅱ，分别与水稻Os02g02120、拟南芥AtWAKL15亲缘关系最近；TraesCS1B02G004100处于分支Ⅲ，与大麦HORVU.MOREX.r3.5HG0503880亲缘关系最近；分支Ⅳ仅包括小麦成员(图2a)。
![]()
图 2 小麦与拟南芥、水稻、大麦 WAK/WAKLs基因中保守蛋白质的系统进化关系(a)、蛋白结构域(b)和基序结构(c)注：红色框为从转录组中筛选的4个WAKs/WAKLs基因；下同。Figure 2. Phylogenetic relationship (a), protein domain(b) and motifs structure (c) based on the amino acid sequences of WAKs/WAKLs from Triticum aestivum and Arabidopsis thaliana, Oryza sativa, and Hordeum vulgareNote: In the red box are four WAKs/WAKLs genes screened from the transcriptome; the same as below.蛋白结构域预测结果显示：小麦、大麦、拟南芥和水稻WAKs/WAKLs蛋白序列均包括激酶部位(图2b)。在进化分支Ⅰ中，小麦WAKs/WAKLs蛋白都包括2个GUB_WAK_bind结构域，1个EGF和1个STK，除水稻Os01g26300和拟南芥AtWAKL16编码的蛋白质不具有胞外结构域外，其他成员均具有胞外结构域。在进化分支Ⅱ中，大部分水稻成员不具备胞外结构域，而其他成员均具备；部分小麦成员还具有跨膜TM_FGFR-like结构域。在进化分支Ⅲ中，除小麦TraesCS1B02G004100和大麦HORVU.MOREX.r2.1HG0004690、HORVU.MOREX.r3.5HG0503880编码的蛋白质具有完整的胞内和胞外结构域外，其他成员均不具备。在进化分支Ⅳ中，只有TraesCS3A02G049500编码的蛋白质具有胞外和胞内结构域。说明胞内STK为所有成员共有，预测WAKs/WAKLs家族基因发挥功能与激酶部位有较大关联。
基序预测结果显示：小麦、大麦、拟南芥和水稻的WAKs/WAKLs蛋白序列共包括12种基序(图2c)。进化分支Ⅰ包括11种基序，TraesCS6B02G459100与大麦HORVU.MOREX.r3.3HG0499530编码的蛋白质基序相似；进化分支Ⅱ包括10种基序，TraesCSU02G079300与水稻Os02g02120、大麦HORVU.MOREX.r3.7HG0649430编码的蛋白基序高度相似，TraesCS3A02G122300与拟南芥AtWAKL15编码的蛋白基序高度相似；进化分支Ⅲ包括10种基序，TraesCS1B02G004100与大麦HORVU.MOREX.r3.5HG0503880编码的蛋白基序高度相似；进化分支Ⅳ包括9种基序，均是来自小麦WAKs/WAKLs家族的成员，且编码的蛋白基序基本相似。此外，基序1、2、3、4、5、7、8、11为所有成员共有，均包括STK，预测胞内激酶部位保守性较高，是决定蛋白质功能的关键基序；基序9、10、12包括胞外EGF和半乳糖醛酸结合结构域(galacturonic acid binding domain，GUB)，说明胞外结构域多样，但具有一定的保守性。
2.4 WAKs/WAKLs蛋白激酶部位的同源性比对
为了进一步研究激酶部位的保守性，预测其发挥功能的部位，对涉及的成员进行激酶部位蛋白序列同源性比对分析。该部位通常有2个与激酶发挥功能有关的重要区域，包括蛋白激酶ATP结合位点、丝氨酸/苏氨酸蛋白激酶活性位点和酪氨酸蛋白激酶活性位点(图3)，分析结果显示：胞内结构域同源性达73.34%，蛋白激酶ATP结合位点、丝氨酸/苏氨酸蛋白激酶活性位点和酪氨酸蛋白激酶活性位点保守性强，均符合这2个结合位点序列构成特征。根据蛋白激酶ATP结合位点中不变的天冬氨酸(D)附近是否存在保守精氨酸(R)，将激酶分为RD型激酶和非RD型激酶，由图4可知：TraesCS6B02G459100和TraesCS3A02G122300为RD型激酶，TraesCS1B02G004100和TraesCSU02G079300为非RD型激酶。
![]()
图 3 WAKs/WAKLs蛋白激酶催化域中的保守结构域基序注：a) 蛋白激酶ATP结合位点；b) 丝氨酸/苏氨酸蛋白激酶活性位点或酪氨酸蛋白激酶活性位点。Figure 3. Conserved domain motifs in the catalytic domain of WAKs/WAKLs protein kinasesNote: a) ATP binding; b) serine/threonine protein kinase active site, or tyrosine protein kinases specific active site.![]()
图 4 小麦与拟南芥、水稻、大麦 WAK/WAKLs基因氨基酸多序列同源比对Figure 4. Homology alignment of amino acid sequences of WAKs/WAKLs genes in T. aestivum and A. thaliana, O. sativa, and H. vulgare2.5 小麦 WAKs/WAKLs基因响应条锈病的表达
由图5可知：4个基因在小麦接种条锈菌亲和互作和非亲和互作过程中差异表达。侵染0 h，TraesCS1B02G004100和TraesCS6B02G459100在小麦与条锈菌非亲和互作中被短暂诱导表达后下调，随后在侵染早期(24~48 h)上调表达，而TraesCS3A02G122300则表达相反；TraesCSU02G079300在早期表达量无显著差异。接种48 h后，TraesCS1B02G004100和TraesCSU02G079300表达受到抑制而下调直至接种168 h；TraesCS6B02G459100表达则呈上调直至接种120 h。综上所述，亲和互作过程中，TraesCS1B02G004100和TraesCS6B02G459100几乎不表达，TraesCS3A02G122300早期表达趋势和非亲和互作一致，TraesCSU02G079300在早期表达量低于非亲和互作。
![]()
图 5 4个WAKs/WAKLs基因在抗病品种云337 (Y337)和感病品种云402 (Y402)接种CYR32后的转录组表达量注：“**”表示与上一个时间相比，在0.01水平上差异极显著；下同。Figure 5. Transcriptome expression level of four WAKs/WAKLs genes in resistant variety Yun337 (Y337) and susceptible variety Yun402 (Y402) post inoculation with CYR32Note: “**” indicate extremely significant differences at the 0.01 level compared to the previous time; the same as below.由图6可知：非亲和互作过中，在条锈菌侵染小麦的24~48 h，TraesCS1B02G004100上调表达，TraesCS6B02G459100显著上调表达，TraesCS3A02G122300和TraesCSU02G079300表达受到抑制而下调；在接种条锈菌48 h后，TraesCS1B02G004100和TraesCS6B02G459100先下调后上调表达；TraesCS6B02G459100、TraesCS3A02G122300和TraesCSU02G079300的表达量分别在接种条锈菌120、24和168 h后达到峰值。亲和互作过程中，TraesCS1B02G004100几乎不表达；转录组中几乎未表达的TraesCS6B02G459100，在接种0~24 h被抑制表达，24 h后上调表达，168 h表达量达到峰值，但由于该基因转录丰度较低，可能导致RT-qPCR检测结果不能有效地反映其实际表达情况；TraesCS3A02G122300和TraesCSU02G079300在接种早期下调表达，接种120 h后表达上调。综上所述，4个基因在小麦和条锈菌互作过程中利用转录组和RT-qPCR分析的表达模式基本一致。
![]()
图 6 RT-qPCR验证4个WAK/WAKLs基因在抗病品种云337 (Y337)和感病品种云402 (Y402)接种CYR32后的表达模式Figure 6. Expression patterns of four WAKs/WAKLs in resistant variety Yun337 (Y337) and susceptible variety Yun402 (Y402) post inoculation with CYR32 by RT-qPCR3. 讨论
WAKs/WAKLs基因在植物中属于大家族，目前从拟南芥中鉴定了26个WAKs/WAKLs基因，水稻中125个，大麦中91个，小麦中320个[19-21]，作为病原菌的模式识别受体，它在植物信号转导过程中发挥着重要作用。本课题组从转录组中筛选的4个WAKs/WAKLs基因，其结构域存在一定差异，暗示这4个基因功能上存在差异；通过RT-qPCR验证，阐明TraesCS1B02G004100基因参与调控小麦对条锈病的抗性，研究结果为揭示WAKs/WAKLs基因调控抗病性的功能和机制提供了参考。
本研究的4个WAKs/WAKLs基因具有典型的胞外富含半胱氨酸的GUB_WAK_bind和胞内丝氨酸/苏氨酸激酶结构域。从进化分支上看，小麦TraesCS6B02G459100与水稻OsWAK91的亲缘关系最近，OsWAK91正向调控水稻抵御稻瘟病菌侵染[16]。已有研究报道：2个胞外结构域变异的OsWAK10基因以不同的亲和力结合果胶寡糖，以控制水稻次生壁的合成[22]，由此可以说明：GUB作为首先感知到外界信号的受体，可以与果胶结合激活WAKs/WAKLs基因，以调节植物细胞形态。TraesCS6B02G459100与上述2个水稻WAKs基因同为RD激酶，由此可以推测，该基因可能编码与调节植物生长发育、抗病相关的激酶，相较于其他3个基因，其具有的2个GUB结构域在信号转导过程中的作用机制值得探索。TraesCS3A02G122300与未知功能的激酶聚为同一分支，在抗感材料中都被抑制表达，但在后期感病材料的表达量高于抗病材料，尽管属于RD激酶，但已有报道中也出现能够调控植物抗病反应的RD激酶，如OsWAK91、OsWAK92、OsWAK94和OsWAK14基因是水稻对稻瘟病菌数量抗性的调节因子[16]，因此TraesCS3A02G122300的具体功能机制值得进一步研究。
WAK/WAKL家族胞内的ATP结合位点及丝氨酸/苏氨酸结合位点是调节植物先天免疫的重要结构[23]，尤其是丝氨酸/苏氨酸蛋白激酶活性位点中保守的天冬氨酸附近是否存在精氨酸很大程度上决定了激酶的功能。已知在先天免疫中作为PRRs的15个激酶中有12个是非RD型[24]，表明所有非RD 激酶可能都参与先天免疫。TraesCS1B02G004100和TraesCSU02G079300为非RD型激酶。其中，TraesCS1B02G004100在亲和互作体系中几乎不表达；而在非亲和体系中，该基因早期被诱导，24 h以后被抑制，这与已报道的OsWAK112d基因表达模式一致，且过表达OsWAK112d基因及其突变体分别增加了水稻的易感性和抗性，表明其在植物感染期间是抑制基础抗性的负调控因子[16]，故推测TraesCS1B02G004100可能是小麦响应条锈菌侵染的负调控基因，并且早期可能被PAMPs诱导。TraesCSU02G079300与OsWAK11处于同一进化分支，OsWAK11通过GUB结合果胶并监测果胶甲基化的变化以解毒过量的铜，并且可以调节细胞的伸长率以适应外界环境的变化[25]，其转录丰度的趋势与TraesCS1B02G004100相似，并且在2种互作体系中差异表达，该基因可能也在防御反应中发挥作用。
本研究筛选的小麦WAK1基因TraesCS1B02G004100和已经在小麦中鉴定到的Snn1基因同源性最高，达99.86%。在腐生病原菌(Parastagonospora nodorum)和小麦互作过程中，SnTox1作为病原菌分泌的效应蛋白与Snn1基因相互作用，引起小麦细胞程序性死亡并赋予其对叶斑病的易感性，Snn1-SnTox1是小麦和叶斑病互作的多个重要模式之一[26-27]。本研究发现：在小麦与条锈菌非亲和互作过程中，TraesCS1B02G004100基因响应条锈菌侵染且表达受到抑制，而在亲和互作过程中不受条锈菌诱导表达，说明该基因极大可能参与调控小麦对条锈病的抗性。
4. 结论
4个小麦WAKs/WAKLs基因中，TraesCS3A02G122300和TraesCSU02G079300在小麦与条锈菌亲和及非亲和互作体系中差异表达；TraesCS6B02G459100和TraesCS1B02G004100在小麦非亲和互作中特异表达。非RD型激酶TraesCS1B02G004100在表达过程中受到抑制，作为病原菌效应蛋白的受体，其在参与调控小麦对条锈病的抗性功能值得进一步研究；其他3个WAKs/WAKLs基因在植物生长发育过程中的作用也值得探索。
-
![]()
图 2 小麦与拟南芥、水稻、大麦 WAK/WAKLs基因中保守蛋白质的系统进化关系(a)、蛋白结构域(b)和基序结构(c)
注：红色框为从转录组中筛选的4个WAKs/WAKLs基因；下同。
Figure 2. Phylogenetic relationship (a), protein domain(b) and motifs structure (c) based on the amino acid sequences of WAKs/WAKLs from Triticum aestivum and Arabidopsis thaliana, Oryza sativa, and Hordeum vulgare
Note: In the red box are four WAKs/WAKLs genes screened from the transcriptome; the same as below.
![]()
图 3 WAKs/WAKLs蛋白激酶催化域中的保守结构域基序
注：a) 蛋白激酶ATP结合位点；b) 丝氨酸/苏氨酸蛋白激酶活性位点或酪氨酸蛋白激酶活性位点。
Figure 3. Conserved domain motifs in the catalytic domain of WAKs/WAKLs protein kinases
Note: a) ATP binding; b) serine/threonine protein kinase active site, or tyrosine protein kinases specific active site.
![]()
图 4 小麦与拟南芥、水稻、大麦 WAK/WAKLs基因氨基酸多序列同源比对
Figure 4. Homology alignment of amino acid sequences of WAKs/WAKLs genes in T. aestivum and A. thaliana, O. sativa, and H. vulgare
![]()
图 5 4个WAKs/WAKLs基因在抗病品种云337 (Y337)和感病品种云402 (Y402)接种CYR32后的转录组表达量
注：“**”表示与上一个时间相比，在0.01水平上差异极显著；下同。
Figure 5. Transcriptome expression level of four WAKs/WAKLs genes in resistant variety Yun337 (Y337) and susceptible variety Yun402 (Y402) post inoculation with CYR32
Note: “**” indicate extremely significant differences at the 0.01 level compared to the previous time; the same as below.
![]()
图 6 RT-qPCR验证4个WAK/WAKLs基因在抗病品种云337 (Y337)和感病品种云402 (Y402)接种CYR32后的表达模式
Figure 6. Expression patterns of four WAKs/WAKLs in resistant variety Yun337 (Y337) and susceptible variety Yun402 (Y402) post inoculation with CYR32 by RT-qPCR
表 1 供试引物
Table 1 Primers for test
引物名称
primer names引物序列 (5′→3′)
primer sequencesGAPDH-F GCCTTGCTCCTCTTGCTA GAPDH-R TCCACCTCTCCAGTCCTT TraesCS1B02G004100-F CAATCAGCGGACGGTCTA TraesCS1B02G004100-R ATGTCGCCGATGACTTGT TraesCS6B02G459100-F TTCATCGGAATGTGGTCAAG TraesCS6B02G459100-R GCTGAGTGGAGATAAGAAAGG TraesCSU02G079300-F GAGATGTCGCTGAGAGGTT TraesCSU02G079300-R CCATGCTGTCCAGTGCTA TraesCS3A02G122300-F AGCTCGACACAATGAAGG TraesCS3A02G122300-R CCCAGCCTCAATGTTGAT 
下载: 导出CSV
表 2 4个WAKs/WAKLs 基因理化性质
Table 2 Physicochemical properties of four WAKs/WAKLs genes
基因
gene蛋白质序列号
accession number of protein多肽长度/aa
peptide length蛋白分子质量/ku
protein molecular weight等电点
pI氨基酸指数
amino acid
index蛋白亲水性
protein
hydrophilicityTraesCS1B02G004100 XP_044361219.1 714 78.68 8.50 89.59 −0.119 TraesCS6B02G459100 XP_044410226.1 960 106.29 5.98 85.33 −0.120 TraesCSU02G079300 XP_044407301.1 737 81.86 5.65 81.38 −0.214 TraesCS3A02G122300 XP_044339162.1 652 69.67 6.25 90.49 −0.024
下载: 导出CSV
-
[1] SHIU S H, BLEECKER A B. Plant receptor-like kinase gene family: diversity, function, and signaling[J]. Science’s STKE, 2001, 2001(113): re22. DOI: 10.1126/stke.2001.113.re22.
[2] VAN DER GEER P, HUNTER T, LINDBERG R A. Receptor protein-tyrosine kinases and their signal transduction pathways[J]. Annual Review of Cell Biology, 1994, 10(1): 251. DOI: 10.1146/annurev.cb.10.110194.001343.
[3] LIU C, HAN L H, CHU H L. Research advances in molecular mechanism between plant and pathogen interaction[J]. Microbiology China, 2018, 45(10): 2271. DOI: 10.13344/j.microbiol.china.170969.
[4] JONES J D G, DANGL J L. The plant immune system[J]. Nature, 2006, 444(7117): 323. DOI: 10.1038/nat- ure05286.
[5] SANABRIA N M, HUANG J C, DUBERY I A. Self/non-self-perception in plants in innate immunity and defense[J]. Self/Nonself, 2010, 1(1): 40. DOI: 10.4161/self.1.1.10442.
[6] 张素巧, 孙颖, 郭毅, 等. 细胞壁连接的类受体激酶(WAK): 植物细胞壁与细胞质联系的重要蛋白[J]. 科学通报, 2006, 51(1): 1. DOI: 10.1360/csb2006-51-1-1. [7] 师莹莹, 李大勇, 张慧娟, 等. 植物细胞壁介导的抗病性及其分子机制[J]. 植物生理学报, 2011, 47(7): 661. DOI: 10.13592/j.cnki.ppj.2011.07.009. [8] BRUTUS A, SICILIA F, MACONE A, et al. A domain swap approach reveals a role of the plant wall-associated kinase 1 (WAK1) as a receptor of oligogalacturonides[J]. Proceedings of the National Academy of Sciences, 2010, 107(20): 9452. DOI: 10.1073/pnas.1000675107.
[9] 左为亮. 玉米抗丝黑穗病主效QTL的克隆与抗病机理研究[D]. 北京: 中国农业大学, 2014. [10] WANG P, ZHOU L, JAMIESON P A, et al. The cotton wall-associated kinase GhWAK7A mediates responses to fungal wilt pathogens by complexing with the chitin sensory receptors[J]. The Plant Cell, 2020, 32(12): 3978. DOI: 10.1105/tpc.19.00950.
[11] 邓俊杰. 山新杨PdPapWAK基因抗逆功能研究[D]. 哈尔滨: 东北林业大学, 2021. [12] 马占鸿. 中国小麦条锈病研究与防控[J]. 植物保护学报, 2018, 45(1): 1. DOI: 10.13802/j.cnki.zwbhxb.2018.2018900. [13] 何锋, 王长春, 王锋青, 等. 植物抗病基因 ( R) 与病原物无毒基因 ( Avr) 相互作用机制的研究进展[J]. 中国细胞生物学学报, 2011, 33(9): 1037. [14] QI H J, ZHU X L, GUO F, L, et al. The wall-associated receptor-like kinase TaWAK7D is required for defense responses to Rhizoctonia cerealis in wheat[J]. International Journal of Molecular Sciences, 2021, 22(11): 5629. DOI: 10.3390/ijms22115629.
[15] SAINTENAC C, LEE W S, CAMBON F, et al. Wheat receptor-kinase-like protein Stb6 controls gene-for-gene resistance to fungal pathogen Zymoseptoria tritici[J]. Nature Genetics, 2018, 50(3): 368. DOI: 10.1038/s41588-018-0051-x.
[16] DELTEIL A, GOBBATO E, CAYROL B, et al. Several wall-associated kinases participate positively and negatively in basal defense against rice blast fungus[J]. BMC Plant Biology, 2016, 16(1): 17. DOI: 10.1186/s12870- 016-0711-x.
[17] ROSLI H G, ZHENG Y, POMBO M A, et al. Transcriptomics-based screen for genes induced by flagellin and repressed by pathogen effectors identifies a cell wall-associated kinase involved in plant immunity[J]. Genome Biology, 2013, 14(12): R139. DOI: 10.1186/gb-2013- 14-12-r139.
[18] 梁子英, 刘芳. 实时荧光定量PCR技术及其应用研究进展[J]. 现代农业科技, 2020(6): 1. DOI: 10.3969/j.issn.1007-5739.2020.06.001. [19] ZHANG S B, CHEN Z J, LI L, et al. Evolutionary expansion, gene structure, and expression of the rice wall-associated kinase gene family[J]. Plant Physiology, 2005, 139(3): 1107. DOI: 10.1104/pp.105.069005.
[20] TRIPATHI R K, AGUIRRE J A, SINGH J. Genome-wide analysis of wall associated kinase ( WAK) gene family in barley[J]. Genomics, 2021, 113(1): 523. DOI: 10.1016/j.ygeno.2020.09.045.
[21] XIA X B, ZHANG X, ZHANG Y C, et al. Characterization of the WAK gene family reveals genes for FHB resistance in bread wheat (Triticum aestivum L.)[J]. International Journal of Molecular Sciences, 2022, 23(13): 7157. DOI: 10.3390/ijms23137157.
[22] CAI W G, HONG J, LIU Z Y, et al. A receptor-like kinase controls the amplitude of secondary cell wall synthesis in rice[J]. Current Biology, 2023, 33(3): 498. DOI: 10.1016/j.cub.2022.12.035.
[23] LIU J, LI W J, WU G, et al. An update on evolutionary, structural, and functional studies of receptor-like kinases in plants[J]. Frontiers in Plant Science, 2024, 15(1): 1305599. DOI: 10.3389/fpls.2024.1305599.
[24] DARDICK C, RONALD P. Plant and animal pathogen recognition receptors signal through non-RD kinases[J]. PLoS Pathogens, 2006, 2(1): e2. DOI: 10.1371/journal.ppat.0020002.
[25] YUE Z L, LIU N, DENG Z P, et al. The receptor kinase OsWAK11 monitors cell wall pectin changes to fine-tune brassinosteroid signaling and regulate cell elongation in rice[J]. Current Biology, 2022, 32(11): 2454. DOI: 10.1016/j.cub.2022.04.028.
[26] LI D, WALKER E, FRANCKI M. Genes associated with foliar resistance to Septoria nodorum blotch of hexaploid wheat (Triticum aestivum L.)[J]. International Journal of Molecular Sciences, 2021, 22(11): 5580. DOI: 10.3390/ijms22115580.
[27] SHI G J, ZHANG Z C, FRIESEN T L, et al. The hijacking of a receptor kinase-driven pathway by a wheat fungal pathogen leads to disease[J]. Science Advances, 2016, 2(10): e1600822. DOI: 10.1126/sciadv.1600822.
计量
- 文章访问数: 967
- PDF下载量: 32
