0 引言
【研究意义】苹果是世界上栽培最广泛的水果之一,因其味道鲜美,且富含营养价值,深受人们喜爱。作为多年生木本植物,其产量和果实品质易受生物和非生物胁迫的影响,特别是病原菌侵害,给苹果产业造成巨大的经济损失[1]。水杨酸(salicylic acid,SA)作为一种植物激素,不仅能够调控植物的生长发育和代谢过程,而且在生物和非生物胁迫反应中起到积极的调控作用[2]。研究表明SA在鼠李糖脂诱导的植物抗性反应中处于中心位置,SANCHEZ等[3]通过对细菌中鼠李糖脂诱导的植物对活体营养型病菌Hyaloperonospora arabidopsidis、腐生营养型病菌Botrytis cinerea及兼性寄生菌Pseudomonas syringae pv. tomato的抗性研究发现,SA参与了这3种病菌的抗性,而乙烯仅仅参与了Hyaloperonospora arabidopsidis和Pseudomonas syringae pv. tomato的抗性反应,茉莉酸仅参与了对Botrytis cinerea的抗性。外源SA处理能够增强苹果叶片对炭疽病的抗性[4],在富士苹果中过表达SA信号途径的关键调控因子NPR1能够提高对白粉病的抗性[5]。外源SA处理能够诱导辣椒的抗疫病[6],以及水稻幼苗对白叶枯病的抗性[7]。除此之外,SA还能够提高植物对非生物胁迫的抗性,如SA能够诱导黄瓜中冷害响应基因的表达,并减弱冷害胁迫对黄瓜的伤害,增强植物对冷害的抗性[8]。SA信号途径对植物抗病抗逆性调控具有重要作用。虽然根据前人的研究结果,显示SA在植物的抗病抗逆信号转导中具有举足轻重的作用,但是,在苹果这种多年生的高等植物中,SA的信号转导途径一直不清楚。利用RNA-seq探究SA对苹果叶片转录调控的影响,对解析SA信号转导途径具有重要意义。【前人研究进展】尽管SA在模式植物拟南芥中取得了很大的进展,但在苹果中,SA如何进行信号转导,依然知之甚少,苹果SA 信号通路目前尚不明确。张计育等[9]构建了SA处理湖北海棠的全长cDNA文库,并克隆了与病菌抗性相关的多聚半乳糖醛酸酶抑制蛋白基因PGIP全长。不同浓度SA处理对PGIP表达产生不同影响,高浓度SA促进苹果叶片中PGIP的表达,低浓度SA则抑制其表达[10]。外源SA提高了苹果叶片中过氧化氢酶(CAT)、多酚氧化酶(PPO)、过氧化物酶(POD)及苯丙氨酸解氨酶(PAL)的活性并诱导PR5和PR8的上调表达[11]。罗昌国等[12]克隆了调控富士苹果白粉病抗性的关键基因WRKY40b,该基因表达受SA的强烈诱导。SA通过与靶蛋白互作发挥其功能,目前其调控途径的靶蛋白在很大程度上还未知[13]。SA信号通路主要围绕NPR1、TGA转录因子(TGACA元件绑定因子)及病程相关蛋白基因PR等进行研究。NPR1只在SA诱导条件下表达,是调控SA介导的系统获得抗性(SAR)的关键基因。拟南芥npr1突变体无法诱导PR的表达,增强了对病原菌的感病性[14]。研究表明TGA1-TGA7均可与NPR1互作,调控SA响应基因[15-16]。NPR1增强TGA转录因子与防卫基因启动子as-1元件的结合[17]。目前SA诱导植物抗病性研究主要利用不同植物和病原菌的互作。例如外源SA诱导显著减弱了番茄维管束褐变及叶片黄萎现象,增强对镰刀菌的抗性[18]。【本研究切入点】苹果SA信号途径研究较少,其信号通路和分子机制有待明确。【拟解决的关键问题】RNA-seq是利用高通量测序研究细胞内基因表达的方法,本研究以SA处理的苹果叶片为研究试材,采用RNA-seq的研究方法,以获得在苹果中响应水杨酸的基因,构建能够对SA进行特异性响应的启动子鉴定体系,为深入研究SA信号转导途径奠定基础。1 材料与方法
1.1 植物材料及处理
供试材料为‘嘎啦’组培苗叶片,组培苗于24℃,12 h/12 h光暗环境下培养30 d,选取长势一致的叶片、将其浸入2 mmol∙L-1 SA溶液中,24℃下处理12 h,滤纸吸干表面水分,液氮迅速冷冻,-80℃保存,以备RNA提取;以0.2%乙醇水溶液处理的叶片为对照。重复3次。1.2 RNA提取及测序
利用天根公司的植物总RNA提取试剂盒(DP424)提取RNA(操作步骤见说明书)。检测合格后,经DNase I处理,用带有Oligo(dT)的磁珠富集mRNA,高温条件下加入适量打断试剂使其片断化,再以片段mRNA为模板合成cDNA,经过磁珠纯化、末端修复、3′末端加碱基A、加测序接头后,进行PCR扩增,完成文库制备。Agilent 2100 Bioanalyzer和ABI StepOnePlusRealTime PCR System进行文库质量和产量检测。测序得到的数据经碱基识别分析后得到原始测序序列,去除含接头、含N比例大于10%、低质量的序列后得到适合分析的数据。使用比对软件BWA[19]和Bowtie[20]将数据比对‘金冠’苹果基因组。1.3 差异基因筛选、功能分析及条件特异性表达分析
差异基因筛选:利用RSEM工具进行基因表达定量,表达定量结果以FPKM(fragments per kilobase of transcript per million fragments mapped)为单位,FDR(false discovery rate)≤0.001且倍数差异在2倍以上的基因为差异基因。将差异基因与Nr数据库比对后进行注释。条件特异性表达分析:条件特异性表达分析通过几种不同条件下差异基因的比较,鉴定在某种特定条件下才表达的基因。利用3种植物激素处理苹果叶片进行转录组学研究,对CTRL1(无菌水)、CTRL2(0.2%乙醇)、水杨酸(2 mmol∙L-1 SA)、1-氨基环丙烷羧酸(10 μmol∙L-1 ACC)、茉莉酸甲酯(100 μmol∙L-1 MeJA)这5个处理样品的测序数据进行分析,根据其他两类激素转录组学的研究结果,应用如下统计学研究方法,得到SA特异性诱导表达的基因。其计算公式为:
\[P_i(g)=\sum
{E(g)}{X=e_i(g)}(
{E(g)_X})(1-P_i)
{E(g)-X}
其中,g为基因,i为样品,x=ei(g)为基因g在样品i中的数据数目,E(g)为基因g在所有样品中数据数目。si为样品i中所有数据数目,Pi=si/∑isi。P≤0.05的基因定义为条件特异表达基因。
差异基因功能分析:所有SA处理得到的差异基因提交Gene Ontology数据库(http://www.geneontology.org/),用WEGO软件[21]对其进行GO功能分类分析。数据提交KEGG数据库(kyoto encyclopedia of gene and genomes),以KEGG途径为单位,Qvalue≤0.05的途径定义为在差异基因中显著富集的途径,分析差异基因参与的最主要生化代谢途径和信号途径。
1.4 差异基因的qRT-PCR荧光定量分析
采用RevertAid First Strand cDNA Synthesis Kit试剂盒(Thermo公司)对1.2中RNA进行反转录,以SYBRGreen染料进行荧光定量PCR(qRT-PCR)分析。苹果肌动蛋白基因ACTIN(序列号为GQ339778)作为内参,基因表达倍数通过2-△△Ct法计算。3次重复。差异基因的特异性引物由生工生物工程有限公司合成(详细引物序列见表1)。Table 1
表1
表1用于qRT-PCR分析的基因及引物
Table 1The qRT-PCR primers for 24 genes
基因号 Transcript-id | 引物序列 Primer sequence | |
---|---|---|
MDP0000830129 | F: CTTCTTCCACTTCGTCAT | R: TCATAATCAAGGAGGTCTTC |
MDP0000449615 | F: TGGACACCTGAATTGATG | R: CCACCAAGAACCAACTAA |
MDP0000383328 | F: ATGTAAGAAGGAAGGAAGAC | R: TATATGCTGCTCCGCTAT |
MDP0000289278 | F: ATGAAGAAGTGCGATCTG | R: GAGTTCTTGAATGCTTGG |
MDP0000270446 | F: GTTGGCTTCCTTCTGATTA | R: TGTCCTTCCGATCTTAGT |
MDP0000226971 | F: GATGTGGAATTGCTTGTC | R: GAGGTCCTTCATAAGAGTAG |
MDP0000168650 | F: CGAGACTTACATCTTTGCT | R: GTCCAGTTGGTGTAGAAG |
MDP0000263349 | F: CGAGACTTACATCTTTGCT | R: GTCCAGTTGGTGTAGAAG |
MDP0000247896 | F: TCGGAACAGGATGAAATTAG | R: TGAAAATTCATCCGAAAACC |
MDP0000710349 | F: TGCTGCTCGGTCTTTTAATG | R: CCATGCATATGGACCATCTG |
MDP0000450565 | F: CTGATTTCGAAGCCCATTTATC | R: GGTCAGTATCTCCATCAATACA |
MDP0000785886 | F: GAGGATCTGAAGAGCTACAT | R: CAATTTTATCATGGCCTTGA |
MDP0000774112 | F: GGAAGAGTTGGAAGTTGAAG | R: CTTGCAGCCTATCTATGAATC |
MDP0000300920 | F: CTTTCGCTGCTATGATTGAT | R: GAGAAGTGGTTGTAGTTGTC |
MDP0000254260 | F: AAGAAGACGATGGAGATCAA | R: CTTCAACATCCGGTGAATTC |
MDP0000767097 | F: CTATTTCTCAATGTCTCAGCAG | R: GACCGTCAAATGAGAAGTTG |
MDP0000293114 | F: AATCATGCCATGCCTTTC | R: GCTCATCTCCACTGTCTT |
MDP0000256492 | F: AAATAGCGACTTTGGGAAAC | R: GTTGTTCTCGGGGAGTATTC |
MDP0000138686 | F: GGAACCATATCCAAGGATTG | R: TGTGTGACCTGATCTCTT |
MDP0000280322 | F: ATGAGTGAGGATGATGATGA | R: AGCAGTAGTAGCAGTATCAT |
MDP0000277999 | F: TATCACAGGGAATGGACAAA | R: GTGAGGATTCGAGACATTTG |
MDP0000868782 | F: CCTCCCTCATTTGTTTTACC | R: GTTCCCTAACCTTCTGATTTG |
MDP0000232344 | F: TACATCACTCAAGACGAACT | R: AATCTGCATTGCCTTTCTTC |
MDP0000711379 | F: ATGACAATGTAGCAGGCTAT | R: ATAACTGTAGTCGGCTTTCT |
新窗口打开
1.5 SA信号途径特异性启动子的鉴定及其生物信息学分析
根据qRT-PCR分析结果,选定SA诱导上调最显著的基因以及受SA诱导调控不显著的基因,克隆其启动子序列,将启动子克隆到含有编码荧光素酶基因的表达载体中,并置于荧光素酶基因之前,用以驱动荧光素酶基因。采用PEG法转化苹果愈伤原生质体,常温表达3 h后,利用2 mmol∙L-1 SA、100 μmol∙L-1 MeJA、10 μmol∙L-1 ACC、10 μmol∙L-1 ABA、5 mmol∙L-1 EGTA(Ca2+螯合剂)分别处理3 h,收集原生质体,加入裂解液,用酶标仪(PerkinElmer公司,型号为VICTOR X4)测定LUC(荧光素酶,luciferase)和GUS值,分析LUC/GUS比值,鉴定SA途径的特异性响应启动子序列。表达载体构建及原生质体分离和转化的详细流程参见HE等[22]的方法。构建该基因的MBP融合表达载体pMAL-gene- 2HA,诱导蛋白表达并进行纯化,用于后续研究。利用PlantCARE数据库(http://bioinformatics. psb.ugent. be/webtools/plantcare/html/)进行启动子顺式作用元件分析。
2 结果
2.1 数据测序质量评估及差异基因、条件特异表达基因的筛选
对照和SA处理分别得到750 439 459 bp和751 596 153 bp测序数据,其中,71.72%和71.29%定位于苹果基因组,44.77%和43.88%的序列与‘金冠’苹果基因组数据[23]中的基因完全匹配(表2)。结果表明:共筛选到3 329个差异表达基因,其中上调基因1 085个,下调基因2 244个(图1)。在本转录组学的研究中,除对SA进行转录组学研究外,还对ACC(乙烯合成前体物质)、MeJA(茉莉酸甲酯)的转录组学进行对比与统计学计算,筛选到苹果基因组中33个只被SA诱导的特异表达基因(图2),包括NAC转录因子(MDP0000252435)、WRKY转录因子(MDP0000263349)、NIMIN1(MDP0000280322)、锌指结构基因(MDP0000289278)、ERF转录因子(MDP0000830129)以及未知功能基因等(表3)。
显示原图|下载原图ZIP|生成PPT
图1差异基因的筛选
左图:横坐标表示CTRL样品表达量的对数值,纵坐标表示SA样品表达量的对数值。橙色三角代表上调基因,蓝色方块代表下调基因,褐色圆点代表非显著差异基因。右图:1代表上调基因,2代表非显著差异基因,3代表下调基因
-->Fig. 1The screening of differential expressed genes
Left chart: X-axis represents log10 (gene expression level of CTRL), Y-axis represents log10 (gene expression level of SA). The orange triangle represents up-regulated genes, the blue block represents down-regulated genes, the brown dot represents no significant difference genes. Right chart: 1 represents up-regulated genes, 2 represents no significant difference genes, 3 represents down-regulated genes
-->
显示原图|下载原图ZIP|生成PPT
图2条件特异表达基因的筛选
-->Fig. 2The screening of condition specificity genes
-->
Table 2
表2
表2样品与‘金冠’苹果基因组序列的比对结果
Table 2The alignment result of genome sequences between sample and ‘golden delicious’ apple
样品 Sample | 总数据 Total Reads | 总碱基对 Total Base Pairs | 总定位序列 Total Mapped Reads | 完全匹配 Perfect Match | 错配 Mismatch | 唯一匹配 Unique Match | 多重位置匹配 Muti-position Match | 总的未定位数据 Total Unmapped Reads |
---|---|---|---|---|---|---|---|---|
对照 CTRL | 15315091 (100.00%) | 750439459 (100.00%) | 10983398 (71.72%) | 6855951 (44.77%) | 4127447 (26.95%) | 9011351 (58.84%) | 1972047 (12.88%) | 4331693 (28.28%) |
水杨酸SA | 15338697 (100.00%) | 751596153 (100.00%) | 10934376 (71.29%) | 6731351 (43.88%) | 4203025 (27.40%) | 9134003 (59.55%) | 1800373 (11.74%) | 4404321 (28.71%) |
新窗口打开
Table 3
表3
表3SA诱导的条件特异表达基因
Table 3Condition specificity genes induced by SA
基因序列号 Gene ID | Nr注释 Nr annotation | 表达比率 Log2 (SA/CTRL) |
---|---|---|
MDP0000830129 | 乙烯响应转录因子 Ethylene-responsive transcription factor | 18.21 |
MDP0000144930 | 赖氨酸组氨酸转运体 Lysine histidine transporter | 15.05 |
MDP0000281098 | 热激蛋白 Heat shock protein | 14.40 |
MDP0000289278 | 锌指结构蛋白 Zinc finger protein | 14.08 |
MDP0000798443 | 油质蛋白 Oleosin | 13.24 |
MDP0000319640 | 糖酯类转运蛋白GLTP1 Glycolipid transfer protein 1 | 13.17 |
MDP0000267022 | NRT1蛋白 NRT1 protein | 12.70 |
MDP0000263349 | WRKY40转录因子 WRKY40 transcription factor | 8.89 |
MDP0000130244 | 黄酮醇合成酶 Flavonol synthase | 7.49 |
MDP0000448602 | 肉桂醇脱氢酶 Cinnamyl alcohol dehydrogenase | 6.71 |
MDP0000280322 | NIMIN1 NIM1-Interacting 1 | 6.56 |
MDP0000592869 | 转录起始因子 Transcription initiation factor | 6.10 |
MDP0000138157 | NADP+还原酶 NADP+ reductase | 6.07 |
MDP0000252435 | NAC蛋白 NAC protein | 5.36 |
MDP0000286915 | 乙烯响应转录因子 Ethylene-responsive transcription factor | 5.12 |
MDP0000715898 | 葡萄糖基转移酶 Glucosyltransferase | 4.88 |
MDP0000190316 | 谷胱甘肽转移酶蛋白 Glutathione transferase | 4.70 |
MDP0000454217 | 外皮分泌蛋白 Skin secretory protein | 4.21 |
MDP0000157971 | 未知蛋白 Unknown protein | 14.92 |
MDP0000159774 | 三角状五肽重复包含蛋白 Pentatricopeptide repeat-containing protein | 3.26 |
MDP0000163041 | 未知蛋白 Unknown protein | 15.01 |
MDP0000253428 | 醣脂类转移相关蛋白 Glycolipid transfer protein-related | 12.77 |
MDP0000262678 | 未知蛋白 Unknown protein | 3.93 |
MDP0000280686 | 短链还原酶家族蛋白 Short-chain reductase (SDR) family protein | 4.72 |
MDP0000386815 | 未知蛋白 Unknown protein | 4.89 |
MDP0000512219 | 未知蛋白 Unknown protein | 14.97 |
MDP0000549101 | 未知蛋白 Unknown protein | 12.99 |
MDP0000682771 | NADH还原酶 NADH dehydrogenase | 13.29 |
MDP0000697527 | 未知蛋白 Unknown protein | 15.63 |
MDP0000729348 | 未知蛋白 Unknown protein | 5.76 |
MDP0000794504 | 未知蛋白 Unknown protein | 12.52 |
MDP0000809400 | 线粒体ATP合酶 Mitochondrial ATP synthase complex. | 13.09 |
MDP0000890154 | DNA绑定/转录因子 DNA binding/transcription factor | 12.80 |
新窗口打开
2.2 差异基因的GO分类
GO分析将差异基因分为参与生物过程、细胞成分和分子功能3部分(图3)。参与生物过程差异基因分为16类,其中细胞过程(cellular process,54.2%)、代谢过程(metabolic process,78.4%)所占比例最多,其次为单细胞过程(single-organism process,33.6%)。细胞成分中,细胞(cell,64.1%)、细胞组分(cell part, 64.1%)所占比例最多,其次为膜结构(membrane,39.6%)以及细胞器(organelle,33%)。分子功能的绑定(binding,55.4%)、催化活性(catalytic activity,77%)差异基因所占比例最多,其次为转运活性(transporter activity,7.2%)、分子转导活性(molecular transducer activity,3.6%)以及核酸绑定转录因子活性(nucleic acid binding transcription factor activity,2.2%)等。显示原图|下载原图ZIP|生成PPT
图3差异基因的GO功能分类
1:生物相;2:生物调节;3:细胞组分组织或生物形成;4:细胞过程;5:发育过程;6:免疫系统过程:;7:定位;8:代谢过程;9:多细胞有机体过程;10:生物过程积极调控;11:生物过程调控;12:复制;13:复制过程;14:刺激反应;15:信号;16:单一生物过程;17:细胞;18:细胞成分;19:胞外区:;20:胞外区部分;21:大分子复合物;22:膜;23:膜附着腔;24:膜组分;25:细胞器;26:细胞器部分;27:抗氧化活性;28:绑定;29:催化活性;30:酶调节器活性;31:分子转导活性;32:核酸绑定转录因子活性;33:蛋白绑定转录因子活性;34:受体活性;35:结构分子活性;36:转运蛋白活性
-->Fig. 3GO function classification of differentially expressed genes
1: Biological phase; 2: Biological regulation; 3: Cellular component organization or biogenesis; 4: Cellular process; 5: Developmental process; 6: Immune system process; 7: Localization; 8: Metabolic process; 9: Multi-cellular organismal process; 10: Positive regulation of biological process; 11: Regulation of biological process; 12: Reproduction; 13: Reproduction process; 14: Response to stimulus; 15: Signaling; 16: Single-organism process; 17: Cell; 18: Cell part; 19: Extracellular region; 20: Extracellular region part; 21: Macromolecular complex; 22: Membrane; 23: Membrane-enclosed lumen; 24: Membrane part; 25: Organelle; 26: Organelle part; 27: Antioxidant activity; 28: Binding; 29: Catalytic activity; 30: Enzyme regulator activity; 31: Molecular transducer activity; 32: Nucleic acid binding transcription factor activity; 33: Protein binding transcription factor activity; 34: Receptor activity; 35: Structural molecular activity; 36: Transporter activity。
-->
2.3 差异基因KEGG富集分析及受SA影响的次生代谢基因
通过KEGG富集分析,确定差异基因参与的最主要生化代谢途径和信号途径。根据P≤0.01,找到30个与整个基因组背景相比显著富集的途径(图4)。KEGG分析的结果发现差异基因显著富集途径主要有4个,其中代谢途径差异基因608个,占29.49%;次生代谢物质合成途径411个,占19.93%;植物-病原菌互作途径235个,占11.4%;植物信号转导途径182个,占8.83%。植物次生代谢途径在植物抵抗生物与非生物胁迫反应中起重要的调控作用。本研究发现在次生代谢合成途径中,SA上调基因124个,下调基因287个(图5)。在24个上调比率≥5的基因中(表4),5个苯丙烷生物合成途径基因,8个类黄酮合成途径基因,3个花青苷生物合成途径基因,4个维生素代谢途径涉及基因,说明SA主要调控苯丙烷和类黄酮生物合成途径。涉及苯丙烷生物合成途径的基因包括羟基脯氨酸富集糖蛋白(MDP0000180381),甘露醇脱氢酶(MDP0000164361、MDP0000609114)和肉桂醇脱氢酶CAD(MDP0000448602、MDP0000858930),其中甘露醇脱氢酶通过调节甘露醇的合成增强植物的
显示原图|下载原图ZIP|生成PPT
图4差异基因的前30个KEGG富集分析途径
1:次生代谢的生物合成;2:代谢途径;3:果糖和甘露醇代谢;4:植物激素信号转导;5:维生素B6代谢;6:昼夜节律;7:光合作用-天线蛋白;8:脂肪酸代谢;9:类黄酮生物合成;10:卟啉和叶绿素代谢;11:甘油酯代谢;12:类胡萝卜素生物合成;13:丙酮酸代谢;14: 缬氨酸、亮氨酸和异亮氨酸退化;15:油菜素内酯生物合成;16:谷胱甘肽代谢:17:半乳糖代谢:18:糖酵解和糖质新生;19:苯丙素的生物合成;20:亚麻酸代谢;21:维生素C代谢;22:植物-病原菌互作;23:其他多糖降解;24:萜类化合物生物合成;25:ABC转运蛋白;26:糖胺聚糖降解;27:二羧酸代谢;28:类单萜生物合成;29:半胱氨酸和蛋氨酸代谢;30:氨基糖和核苷酸糖代谢
-->Fig. 4The top 30 pathways of KEGG enrichment result for differentially expressed genes
1: Biosynthesis of secondary metabolites; 2: Metabolic pathways; 3: Fructose and mannose metabolism; 4: Plant hormone signal transduction; 5: Vitamin B6 metabolism; 6: Circadian rhythm; 7: Photosynthesis - antenna proteins; 8: Fatty acid metabolism; 9: Flavonoid biosynthesis; 10: Porphyrin and chlorophyll metabolism; 11: Glycerolipid metabolism; 12: Carotenoid biosynthesis; 13: Pyruvate metabolism; 14: Valine, leucine and isoleucine degradation; 15: Brassinosteroid biosynthesis; 16: Glutathione metabolism; 17: Galactose metabolism; 18: Glycolysis / Gluconeogenesis; 19: Phenylpropanoid biosynthesis; 20: alpha-Linolenic acid metabolism; 21: Ascorbate and aldarate metabolism; 22: Plant-pathogen interaction; 23: Other glycan degradation; 24: Terpenoid backbone biosynthesis; 25: ABC transporters; 26: Glycosaminoglycan degradation; 27: Glyoxylate and dicarboxylate metabolism; 28: Monoterpenoid biosynthesis; 29: Cysteine and methionine metabolism; 30: Amino sugar and nucleotide sugar metabolism
-->
显示原图|下载原图ZIP|生成PPT
图5前4个显著富集途径中差异基因分析
1:代谢途径;2:次生代谢生物合成;3:植物病原菌互作;4:植物激素信号转导
-->Fig. 5Analysis of DEGs in the top 4 enrichment pathway
1: metabolic pathway; 2: biosynthesis of secondary metabolites; 3: plant-pathogen interaction; 4: plant hormone signal transduction
-->
抗性[24],其上调表达为研究SA增强抗性的分子机制提供理论依据。木质素是植物细胞壁的一部分,增强了细胞壁的韧性,同时增强植物对胁迫反应的抵抗能力。肉桂醇脱氢酶CAD作用于整个木质素合成途径的最后一步,在木质素合成中起到关键作用[25]。本研究中SA诱导CAD的上调表达,间接表明SA增强植物抵抗胁迫反应的能力,为研究SA诱导植物抗性的分子机制奠定了基础。苯丙氨酸解氨酶PAL是苯丙烷代谢途径的关键酶和限速酶,参与植物木质素、类黄酮等次生代谢物质的生物合成,其活性在SA诱导下呈现先上升后下降的趋势[26]。本研究中SA未诱导PAL(MDP0000787168,-3.03)的上调,与前人研究结果不一致[11],很可能是因为SA处理时间过长,导致PAL活性下降。在类黄酮生物合成途径中,莨菪碱双加氧酶(MDP0000472462、MDP0000183313、MDP0000229093),细胞色素P450家族(MDP0000221432、MDP0000196375),异黄酮羟化酶(MDP0000138677、MDP0000165503),黄酮醇合成酶(MDP0000130244)被SA诱导显著上调。其中细胞色素P450一方面参与植物次生代谢的生物合成,一方面参与代谢解毒功能。维生素代谢途径中的醛酮类还原酶(MDP0000804510、MDP0000046192、MDP0000158194、MDP0000246999)在SA诱导下表达量显著提高,它们也是植物中重要的解毒酶,在植物的非生物胁迫反应中发挥重要作用[27]。SA调控花青苷的生物合成,其中3-氧-葡萄糖苷5-氧-葡萄糖基转移酶(MDP0000229280、MDP0000545122、MDP0000506825)在SA诱导下显著上调。β-1,3-葡聚糖酶降解真菌细胞壁,有直接攻击病原真菌的潜在能力[28],在本研究中,SA诱导β-1,3-葡聚糖酶(MDP0000254363)的显著上调,同时SA诱导上调细胞色素P450(MDP0000196375),该基因调控油菜素内酯的生物合成。
Table 4
表4
表4次生代谢生物合成途径中的上调基因(表达比率≥5)
Table 4The up-regulated genes involved in biosynthesis of secondary metabolites pathway (Log2 (SA/CTRL)≥5)
基因号 Gene ID | 注释 Annotation | 表达比率 Log2 (SA/CTRL) | |
---|---|---|---|
MDP0000180381 | 羟基脯氨酸富集糖蛋白 Hydroxyproline-rich glycoprotein family protein | 1 | 12.68 |
MDP0000164361 | 甘露醇脱氢酶 Mannitol dehydrogenase | 11.85 | |
MDP0000609114 | 甘露醇脱氢酶 Mannitol dehydrogenase | 7.53 | |
MDP0000448602 | 肉桂醇脱氢酶 Cinnamyl alcohol dehydrogenase 1 | 6.71 | |
MDP0000858930 | 肉桂醇脱氢酶 Cinnamyl alcohol dehydrogenase 1 | 5.63 | |
MDP0000472462 | 莨菪碱双加氧酶 Hyoscyamine 6-dioxygenase-like | 2 | 12.08 |
MDP0000183313 | 莨菪碱双加氧酶 Hyoscyamine 6-dioxygenase-like | 11.9 | |
MDP0000229093 | 莨菪碱双加氧酶 Hyoscyamine 6-dioxygenase-like | 8.42 | |
MDP0000196375 | 细胞色素P450 Cytochrome P450 CYP72A219-like | 6.69 | |
MDP0000221432 | 细胞色素P450 Cytochrome P450 82C4 | 11.11 | |
MDP0000138677 | 异黄酮羟化酶 Isoflavone hydroxylase | 9.41 | |
MDP0000165503 | 异黄酮羟化酶 Isoflavone hydroxylase | 7.22 | |
MDP0000130244 | 黄酮醇合成酶 Flavonol synthase | 7.49 | |
MDP0000196375 | 细胞色素P450 Cytochrome P450 CYP72A219-like | →3 | 6.69 |
MDP0000229280 | 3-氧-葡萄糖苷5-氧-葡萄糖基转移酶 Anthocyanidin 3-O-glucoside 5-O-glucosyltransferase-like | 4 | 16.82 |
MDP0000545122 | 3-氧-葡萄糖苷5-氧-葡萄糖基转移酶 Anthocyanidin 3-O-glucoside 5-O-glucosyltransferase-like | 6.43 | |
MDP0000506825 | 3-氧-葡萄糖苷5-氧-葡萄糖基转移酶 Anthocyanidin 3-O-glucoside 5-O-glucosyltransferase-like | 6.16 | |
MDP0000233325 | 薄荷醇脱氢酶 (+)-neomenthol dehydrogenase | 5 | 5.88 |
MDP0000702120 | α-蒎烯合酶 Alpha-pinene synthase | 12.6 | |
MDP0000254363 | β-葡聚糖酶 Beta-glucanase | 7.09 | |
MDP0000804510 | 醛酮类还原酶 Aldo-keto reductase | 6 | 6.5 |
MDP0000046192 | 醛酮类还原酶 Aldo-keto reductase | 5.85 | |
MDP0000158194 | 醛酮类还原酶 Aldo-keto reductase | 5.81 | |
MDP0000246999 | 醛酮类还原酶 Aldo-keto reductase | 5.31 |
新窗口打开
2.4 植物病原菌互作途径中受SA影响的抗病相关基因
在进化过程中,植物与病原菌的较量也在自然选择过程中不断改进,使植物形成了复杂的防御机制。植物的抗性分为基础抗性和主动抗性。基础抗性在病原菌入侵前就已存在,是植物抵抗病原菌的第一道防线,如植物本身的气孔、皮孔、角质层等。主动抗性是病原菌突破基础抗性后激发的第二道防御系统,即免疫反应,激活一系列防卫相关基因的表达,此过程需要转录因子等重要抗性功能基因的调控。其中植物与病原菌互作途径在早期抵抗病原菌入侵中发挥重要作用。本研究植物-病原菌互作途径中被SA诱导61个基因上调,174个基因下调(图5)。分析得到了该途径中调控幅度较大的抗病相关基因(表5)。其中MAPK级联反应在植物与病原菌互作过程中具有重要意义。MPK4能够负调控植物的系统获得抗性(SAR)。研究表明拟南芥mpk4突变体SA合成量提高,PR表达增强,表现出组成性的系统获得抗性(SAR),提高了植物对病原菌的抗性[29]。利用外源SA对丹参处理2 h后能够诱导MKK5及MAPK级联途径上游基因EDR1上调表达[30],拟南芥EDR1被SA诱导且负调控SA介导的植物抗病反应[31]。本研究中SA诱导MAPK级联途径上游基因MEKK1(MDP0000431417)上调,但对下游基因MKK4/5的转录水平无影响。MPK3/4/6调控下游的转录因子,促进或抑制转录因子与下游靶蛋白启动子结合,调控防卫基因表达以及植物抗毒素的积累[32]。本研究SA对MPK3/4/6无影响,但下游WRKY转录因子被SA显著诱导表达,包括与调控防卫基因相关的转录因子MDP0000754989、MDP0000228304、MDP0000175240、MDP0000767097等。与MPK3/6相比,MPK4下游的WRKY转录因子MDP0000263349、MDP0000253189、MDP0000300712调控比率较大(表5)。由此推测,SA更侧重于调控MPK4下游WRKY转录因子的表达。
植物病原菌互作途径中上调趋势较为明显的还有抗病蛋白受体基因RPM1(MDP0000774112),PBS1(MDP0000300920),热激蛋白基因HSP90(MDP0000254260),此外,钙调蛋白CaM(MDP0000785886)同样被SA诱导上调表达。
热激蛋白HSP90是一种高度保守且在胁迫反应中大量积累的蛋白,通常以分子伴侣的形式与靶蛋白形成复合物调控生物抗性反应[33]。核酸绑定和寡聚化域类受体(NLRs)感知不同病原菌的受体蛋白,病原菌受体的识别触发植物防卫反应,伴随一系列细胞程序性死亡,即过敏反应[34]。RPM1属于NLRs家族成员,拟南芥中HSP90对于NLRs的积累是必需的[35]。
HSP90与抗病蛋白RPM1结合,可通过改变分子伴侣活性改变RPM1的分子功能。本研究SA处理诱导HSP90(MDP0000254260)和RPM1(MDP0000774112)的上调,为研究SA增强植物抗病性的调控机制提供理论依据。
Table 5
表5
表5植物病原菌互作途径调控基因(表达比率≥3)
Table 5The up/down-regulated genes involved in plant-pathogen interaction pathway (Log2(SA/CTRL) ≥3)
基因序列号 Gene ID | Nr注释 Nr annotation | 表达比率 Log2(SA/CTRL) |
---|---|---|
MDP0000321613 | 核苷酸门控离子通道 Cyclic nucleotide-gated ion channel 1-like isoform | 4.57 |
MDP0000545337 | 钙离子结合蛋白 Calcium-binding protein CML19 | 4.78 |
MDP0000785886 | 钙调蛋白 Calmodulin | 4.29 |
MDP0000431417 | MEKK1 mitogen-activated protein kinase kinase kinase A-like | 5.26 |
MDP0000217310 | MEKK1 mitogen-activated protein kinase kinase kinase ANP1-like | 3.68 |
MDP0000263349 | WRKY 40转录因子 WRKY 40 transcription factor | 8.89 |
MDP0000253189 | WRKY 51转录因子 WRKY 51 transcription factor | 4.67 |
MDP0000754989 | WRKY 70转录因子 WRKY 70 transcription factor | 4.01 |
MDP0000774112 | 抗病蛋白RPM1 Disease resistance protein RPP8-like | 3.82 |
MDP0000414786 | 抗病蛋白 Disease resistance protein | 3.71 |
MDP0000288873 | 丝氨酸/苏氨酸蛋白激酶 Serine/threonine-protein kinase RLK1 | 10.67 |
MDP0000300920 | 类受体蛋白激酶PBS1 Receptor-like protein kinase ANXUR2 | 4.61 |
MDP0000254260 | 热激蛋白HSP90 Heat shock protein 83-like | 3.63 |
MDP0000300712 | WRKY 61转录因子 WRKY 61 transcription factor | -13.18 |
MDP0000767097 | WRKY 53转录因子 WRKY 53 transcription factor | -3.29 |
MDP0000293114 | 抗病蛋白 Putative disease resistance protein RGA3 | -10.82 |
MDP0000223219 | 丝氨酸/苏氨酸蛋白激酶 Serine/threonine-protein kinase CES101 | -11.91 |
MDP0000270500 | 白介素受体相关激酶 Interleukin-1 receptor-associated kinase 4 | -5.41 |
MDP0000256492 | 转录因子MYC2 Transcription factor MYC2 | -12.09 |
新窗口打开
2.5 差异基因qRT-PCR荧光定量分析
选定24个差异表达基因(15个上调基因,9个下调基因)进行qRT-PCR分析。如图6所示,15个上调基因中,除NPR3(MDP0000868782,1.7倍)诱导表达倍数较弱,其他基因均被SA诱导显著上调。WRKY40(MDP0000263349,4 615倍)、UPG0481(MDP0000226971,3 241倍)、钙调蛋白CaM(MDP0000785886,3 219倍)和类固醇还原酶(MDP0000383328,1 293倍)上调倍数最显著,其次为NIMIN1(MDP0000280322,572倍)、HSP90(MDP0000254260,359倍)等,抗病蛋白基因PRM1(MDP0000774112,128倍)、PBS1(MDP0000300920,112倍)也被SA诱导上调。9个下调基因中,TGA(MDP0000277999)、PR1(MDP0000711379)、Rboh(MDP0000138686)、MYC2(MDP0000256492)下调趋势较明显。ChiB-1(MDP0000710349)、WRKY29(MDP0000247896)和CDPK(MDP0000232344)虽然调控趋势不明显,但转录组分析其调控比率均≤1.5,表明SA对这3个基因的调控较弱。综上所述,qRT-PCR结果表明各基因表达与转录组分析结果一致,只是在差异倍数上与转录组数据有所差别。显示原图|下载原图ZIP|生成PPT
图6差异基因的qRT-PCR分析
-->Fig. 6qRT-PCR analysis of DEGs
-->
2.6 苹果中SA特异性响应启动子的鉴定
为进一步研究苹果对SA的响应机制,本研究利用启动子驱动的荧光素酶活性进行分析。选定SA特异性响应基因MdWRKY40作为研究对象,克隆了该基因的启动子序列:从该基因开放阅读框起始位点ATG,向上游2 500 bp的核苷酸序列(图7)。以苹果中对SA不敏感的基因ChiB-1(MDP0000710349)的启动子(图8)及拟南芥中对细菌鞭毛蛋白进行响应的FRK1(At2g19190)的启动子(图9)[36]为参考对照,检测MdWRKY40启动子驱动的荧光素酶活性。荧光素酶活性分析表明,MdWRKY40能够对SA进行响应(图10-A)。结果表明,经SA处理,其荧光素酶活性为未经SA处理的20.6倍;FRK1也能够被SA诱导激活,但是,诱导倍数显著低于MdWRKY40;而另一基因MdChiB-1(MDP0000710349)启动子驱动的荧光素酶,经SA处理后其活性仅为未处理的2.29倍,说明该片段对SA不具有或仅具有微弱的响应能力。
对MdWRKY40基因的启动子序列进行生物信息学分析,发现序列中存在ABA/MeJA/GA/SA响应顺式作用元件、病菌响应元件、防卫胁迫响应元件等(图7)。但是,进一步应用该启动子序列的研究表明,除SA外,其他激素不能够引起该启动子片段的响应(图10-B),说明该片段是对SA进行特异性响应的启动子序列。同时,对FRK1和ChiB-1的启动子序列进行分析发现,FRK1启动子序列含有4个SA响应元件位点(图9),而ChiB-1启动子序列中不含SA响应元件(图8)。
为进一步精细区分界定MdWRKY40启动子中对SA响应的顺式作用元件,分别克隆了从MdWRKY40开放阅读框的翻译起始位点ATG向上游500 bp、1 000 bp、1 500 bp、2 000 bp的核苷酸序列,并构建到表达载体中置于荧光素酶编码基因的上游,转化原生质体,检测对SA的响应能力。研究表明,500 bp、1 000 bp的核苷酸序列,不具有对低浓度SA(0.2 mmol∙L-1)的响应能力,而在1 500 bp时,对SA表现出微弱响应,但与1 000 bp相比,二者在统计学上差异不显著;而长度达到2 000 bp时,显著提高了荧光素酶活性,为未用SA处理的4.1倍,分别是500 bp、1 000 bp时SA处理条件下荧光素酶活性的2.3倍、3.3倍。说明,在1 000 bp—2 000 bp间具有显著提高SA响应能力的核苷酸序列。而2 000 bp、2 500 bp的核苷酸序列在响应SA的水平上,不具有显著差异(图10-C,左)。
应用高浓度的SA处理,同样证明了这一结论。用0.5 mmol∙L-1 SA处理原生质体细胞,发现500 bp、1 000 bp的核苷酸序列即具有对SA的响应能力,荧光素酶活性分别为未用SA处理的4.2倍、3.3倍,500 bp、1 000 bp二者响应能力处于同一水平。而长度为1 500 bp时,荧光素酶活性进一步升高,为未用SA处理的7.4倍,与1 000 bp、500 bp核苷酸序列相比,具有显著差异。核苷酸序列长度增加到2 000 bp时,对SA响应能力进一步增强,荧光素酶活性升高倍数为11.5倍,与1 500 bp核苷酸序列相比,具有显著差异。该研究结果说明,在1 000 bp—1 500 bp及1 500—2 000 bp,具有增强SA响应能力的顺式作用元件。而2 000—2 500 bp,在对SA响应的荧光素酶活性变化上,没有进一步提高(图10-C,右)。
显示原图|下载原图ZIP|生成PPT
图7MdWRKY40启动子序列的生物信息学分析
-->Fig. 7Bioinformatic analysis of MdWRKY40 promotor sequence
-->
结合低浓度(0.2 mmol∙L-1)与高浓度(0.5 mmol∙L-1)SA的研究结果可以看出,从ATG开始向上游500 bp的核苷酸序列已具有对SA的响应能力,但是,该区段不能响应低浓度SA,而在1 000—1 500 bp及1 500—2 000 bp,尤其是后者,具有显著增强苹果细胞对SA敏感性的未知核苷酸序列。
显示原图|下载原图ZIP|生成PPT
图8MdChiB-1启动子序列的生物信息学分析
-->Fig. 8Bioinformatic analysis of MdChiB-1 promotor sequence
-->
显示原图|下载原图ZIP|生成PPT
图9FRK1启动子序列的生物信息学分析
-->Fig. 9Bioinformatic analysis of FRK1 promotor sequence
-->
2.7 MdWRKY40转录调控存在反馈抑制机制
研究表明,一些重要的转录因子存在自我催化或反馈抑制的转录调控机制,以达到快速响应外界信号或及时阻抑过高表达的目的,前者如调控苹果果实着色的重要转录因子MdMYB10[37]与欧芹PcWRKY1[38],后者如拟南芥中调控植物衰老进程的重要转录因子WRKY6[39]。为探究MdWRKY40是否存在自我调控机制,将CaMV 35S启动子驱动的MdWRKY40蛋白与Promoter- WRKY40-Luciferase在原生质体中进行共转化。共转后的荧光素酶活性仅为对照的35%,Western blot显示MdWRKY40蛋白在原生质体中成功表达,其分子量为37 kD(图10-D),说明MdWRKY40蛋白的表达对其自身RNA水平的转录具有反馈抑制作用。
在本研究中,还克隆了MdWRKY40的开放阅读框序列,构建到含有MBP融合蛋白表达载体中,转化大肠杆菌BL21,成功得到WRKY40诱导蛋白(MBP蛋白分子量为40 kD,加上WRKY40后,总的蛋白分子量约为77 kD,图10-E),为解析该蛋白在SA信号途径的功能奠定了生化基础。
显示原图|下载原图ZIP|生成PPT
图10苹果原生质体鉴定SA途径中的特异性启动子
A:SA处理各启动子的LUC/GUS比值;B:各种激素处理MdWRKY40启动子的LUC/GUS比值;C:MdWRKY40启动子的不同区段对不同浓度SA处理的LUC/GUS比值,0.2 mmol∙L-1 SA(左)、0.5 mmol∙L-1 SA(右);D:MdWRKY40蛋白超表达对LUC/GUS比值的影响;E:MdWRKY40诱导蛋白表达。*表示P<0.05,**表示P<0.01,***表示P<0.001
-->Fig. 10The identification of special promotor response to SA in apple protoplasts
A: LUC/GUS expression of the identified promotors with SA treatment; B: LUC/GUS expression of WRKY40 promotor with plant hormone treatment; C: LUC/GUS expression of different regions of WRKY40 promotor with different concentrations of SA treatment, 0.2 mmol∙L-1 SA (left), 0.5 mmol∙L-1 SA (right); D: Identification of WRKY40 self-regulation; E: The induce expression of WRKY40 protein. * represents P<0.05,** represents P<0.01,*** represents P<0.001
-->
3 讨论
3.1 转录组数据分析筛选到SA介导抗病信号途径的重要功能基因
由于SA调控苹果抗病的分子机制尚不明确,SA信号途径的关键基因尚未找到,本研究利用外源SA处理苹果叶片,构建SA诱导的转录组文库,并进行转录组测序分析,探讨SA处理后差异表达基因、基因富集通路的变化,寻找SA调控抗病反应的关键基因。对照和SA处理测序数据可靠,覆盖率高,适合用于数据分析。NPR1(nonexpressor of PR genes 1)是SA信号途径的关键调控因子[40]。拟南芥和小麦过量表达NPR1能够增强二者对病原菌的抗性[41]。CHAI等[42]的研究表明,利用0.1 mmol∙L-1 SA处理拟南芥离体叶片2 d,NPR1转录水平表达倍数为对照的3倍,处理4 d,NPR1的表达提高到9倍,而在本研究中,未检测到SA处理后苹果NPR1(MDP0000292425)基因表达的变化,一方面可能与SA处理浓度(2 mmol·L-1)与取样时间(12 h)有关;另一方面,在SA信号转导中,NPR1作为SAR(系统获得抗性)的分子开关,其泛素化降解启动SAR信号转导过程,而NPR1旁系同源基因NPR3和NPR4作为SA的受体,介导了NPR1泛素化降解[43]。在本研究中NPR3(MDP0000868782)被SA处理诱导表达(表达比率为1.7)。因此,推测NPR1在表达水平上不受调控,而主要发生在蛋白水平上。NPR1表达与SA处理浓度及时间的关系,以及与SA信号转导的关系,需要进一步研究。NIMIN1(MDP0000280322)是NPR1的互作蛋白,负调控PR表达[44-45]。本研究NIMIN1(MDP0000280322)被SA诱导显著上调(其表达比率为6.56),是SA诱导的条件特异表达基因之一,且PR1(MDP0000711379,-1.14)的转录水平同时降低,因此推测苹果与拟南芥NIMIN1负调控PR1的分子机制一致。其在苹果中的具体功能尚不清楚,需要进行后续研究。
在本研究中,一些调控次生代谢的基因也能够被SA诱导上调,其中包括葡萄糖基转移酶,该基因在烟草中同样被SA所诱导表达[46-47]。研究表明,病菌通过分泌甘露醇到质外体来抑制活性氧介导的的宿主防卫反应,植物中的甘露醇脱氢酶通过将病菌分泌的甘露醇转变为甘露糖,从而增强植物的抗病反应,本研究中甘露醇脱氢酶基因被SA诱导上调表达,这与CHENG等[48]的研究一致。肉桂醇脱氢酶(CAD)是木质素合成的关键酶,对转录组数据进行分析后发现,SA上调CAD的表达,在甘薯和茶树中同样证实了这一观点[49-50]。表明不同物种之间SA响应次生代谢途径的保守性很强。SA对次生代谢的调控,是调控植物抗性的一个非常重要的方面。
研究中还发现了参与活性氧途径代谢的基因,在植物防卫反应中,还有一种重要的次生代谢产物活性氧(ROS),对触发过敏反应至关重要[51]。呼吸爆发氧化酶(Respiratory burst oxidase homologues,Rboh)是植物活性氧(ROS)的主要生产者[52],PTI(pattern triggered immunity)和ETI(effector triggered immunity)均能调控植物防卫反应,但ETI能够更强更快速的引起局部过敏反应[53]。SA处理下调呼吸爆发氧化酶(Rboh)的转录水平,表明PTI途径中由ROS引发的过敏反应被抑制,同时SA诱导苹果中HSP90(MDP0000254260)和RPM1(MDP0000774112)的上调,激发ETI途径中抗病基因介导的过敏反应,对于SA诱导苹果抗病性的调控机制指明了研究方向。SA对非生物胁迫起到正面调控作用。研究表明SA处理显著降低苹果幼苗内ROS水平,提高抗氧化酶活性,降低膜脂过氧化水平,从而减轻低氧胁迫对植株的伤害,其调控机制尚不明确[54]。本研究中SA诱导呼吸爆发氧化酶Rboh下调,为SA能够增强苹果对低氧胁迫的抗性提供了分子证据。
研究表明持续的MAPK级联途径的激活有助于激发SA诱导的大部分响应基因的表达[55]。本研究SA诱导了苹果MAPK级联途径顶端基因MEKK1的上调,同时诱导下游WRKY转录因子的表达,但对MKK4/5及MPK3/4/6则无转录水平的影响。对差异基因进行条件特异表达分析发现SA特异性诱导了33个差异基因的表达,包括NAC转录因子、NIMIN1、WRKY40、ERF转录因子等,这在之前苹果转录组数据分析中未见报道。33个特异表达基因为揭示SA途径的特异性生物过程提供了理论依据,其中的WRKY转录因子及植物病原菌互作相关蛋白,在SA诱导的植物抗病反应中发挥重要作用。PANDEY等[56]的研究表明水稻中多种WRKY转录因子参与抗病反应,例如 WRKY45、WRKY53和WRKY89在稻瘟病菌的侵染过程中表现为差异性表达,白叶枯菌能够诱导WRKY62和WRKY76的特异性表达,并与WRKY28和WRKY71互作从而负调控XA21介导的抗病反应[57-58],本研究诱导的WRKY转录因子在苹果抗病反应的具体功能尚不清楚,仍需进行深入研究。
从本研究转录组的数据可以看出,SA显著调控NIMIN,次生代谢合成基因,活性氧代谢基因以及WRKY转录因子家族基因表达,SA广泛参与了对于苹果系统抗性中各个过程的调控,其中,对每一方面的调控方式和作用位点,需要在本研究基础上,进一步应用分子以及生化手段进行揭示。
3.2 原生质体转化技术与SA信号途径
启动子负责对外界诱导条件的响应,在本研究中只克隆了2 500 bp的核苷酸序列,有可能在上游序列中,仍然含有进一步提高对SA进行响应的特异性元件。FRK1为拟南芥中对细菌鞭毛蛋白flg22的22个氨基酸序列进行强烈诱导表达的基因,是细菌鞭毛蛋白响应的指示性基因[36],在苹果基因组中并未找到FRK1的同源基因,但是,将FRK1-LUC转化于苹果细胞中,经SA处理,发现FRK1升高10.1倍,说明在苹果细胞中FRK1启动子序列仍然能够被识别,这不仅表明苹果细胞与拟南芥细胞中,信号转导途径的保守性,同时也说明,在苹果细胞中,SA信号途径具有影响PTI信号转导的能力,说明SA在植物抗病信号转导中,具有重要的作用。在FRK1的启动子中,发现了对SA进行响应的4个序列元件,暗示FRK1的启动子除了能够对细菌鞭毛蛋白flg22进行响应外,还能够对SA进行响应,而试验结果也表明FRK1的启动子可以响应SA,与生物信息学分析结果相一致。启动子序列通常包含多种顺式作用元件,确定对信号应激反应起关键作用的顺式作用元件对理解细胞的信号转导极为重要。在本研究中,通过不同区段克隆的方式,发现在1 000—2 000 bp,存在增强SA信号响应能力的作用元件。生物信息学分析表明,在1 000 bp以内,存在3个响应SA的顺式作用元件,其中500 bp以内含有2个,500—1 000 bp含有1个(图7),这与500 bp的核苷酸序列已具有对SA的响应能力的试验结果相一致,但是,500 bp、1 000 bp对SA的响应能力没有区别,说明这3个顺式元件可能功能上具有重叠。同时,该区段不能响应低浓度SA,而当核苷酸序列长度增长到2 000 bp时,具有了对低浓度SA的响应能力,说明,在1 000—2 000 bp,含有进一步增强SA信号转导能力的作用元件。将这一区段一分为二,即1 000—1 500 bp、1 500—2 000 bp,发现这两个区段在高浓度SA时,都能够显著增强SA信号转导的功能,说明在这两个区段中,分别具有提高对SA敏感性的作用元件,而在低浓度SA时,只有从1 000 bp增加到2 000 bp时,才能显著提高苹果细胞对SA的敏感性,说明这两个区段的SA顺式作用元件具有累积与加和效应。在启动子的研究中,经常发现,多个元件具有累积效应,元件拷贝数目的不同会对启动子响应的强度和诱导性产生不同的影响[59]。生物信息学分析表明,在1 000—1 500 bp存在2个茉莉酸甲酯响应元件(TGACG),1 500—2 000 bp存在3个病菌响应元件(ACGTG)、2个生长素响应元件(AACGAC)和1个W-box(TTGAC),并未预测到SA的响应元件,说明在这一区段中,可能存在增强细胞对SA敏感性的未知核苷酸序列。
在植物中一些具有重要功能的转录因子具有自我催化或自我抑制机制。研究发现,MdWRKY40蛋白能够抑制自身基因的转录,这种自我抑制模式与拟南芥中调控衰老的关键基因WRKY6的表达模式相同。ROBATZEK等[39]发现,将WRKY6蛋白在拟南芥细胞中超表达,其自身的转录活性降低12倍。这种反馈抑制机制可能有效阻止了由自身基因的过量表达而导致的植物早衰。在本研究中发现MdWRKY40能够快速高效响应SA,说明MdWRKY40是抗性信号物质SA信号转导中的一个重要基因,该基因的反馈抑制作用,可能及时阻止了自身过量表达对植物带来的危害。笔者实验室已经构建植物稳定超表达载体,进行转化苹果愈伤和苹果植株,对其基因功能进行进一步研究。
原生质体转化技术广泛应用于植物的基因功能、亚细胞定位、蛋白互作及瞬时表达等方面,可以高效准确的解析植物细胞信号转导途径,是分子生物学和生物化学研究的有效手段[60]。不同植物原生质体技术的成功应用,促进了植物基因组学及蛋白质组学功能的研究。例如,水稻中利用原生质体转化技术成功得到了MYBS3的亚细胞定位性质[61]。拟南芥中利用原生质体研究了不同基因的表达、信号转导途径[62]及蛋白互作情况[63]等。本研究利用苹果愈伤的原生质体转化技术首次鉴定了在苹果中对SA进行特异性响应的启动子序列,即MdWRKY40(MDP0000263349)的启动子序列,为苹果这种多年生木本植物中开展相关激素信号转导途径的研究提供了稳定的鉴定体系,也为以后获得更多特异途径的靶标基因与调控基因提供了可能。大肠杆菌中成功诱导WRKY40蛋白的表达,为研究WRKY结构域的结合元件以及鉴定下游的靶基因奠定了生化基础。WRKY40的启动子序列中包含各种激素响应元件、防卫胁迫响应元件等,表明WRKY40可能调控多种激素和胁迫反应过程。作为SA信号途径的特异性启动子,SA响应元件在诱导WRKY40表达过程中占主导地位,图8中3个位点的SA响应元件对于WRKY40是否存在功能冗余,可对3个SA响应元件进行缺失突变,利用原生质体转化技术进一步鉴定其关键调控元件。
SA调控的途径复杂而精细,在植物抗性反应中起到显著作用[64-66]。目前苹果中关于SA诱导抗病的分子机制尚不完善。与模式植物相比,苹果中SA诱导的途径更为复杂,是苹果抗病的重要机制之一。分析SA诱导苹果叶片转录水平的基因表达差异,利于更好的理解SA信号途径与植物病原菌互作的分子机制。本研究的转录组学数据及启动子的克隆,为SA在苹果中信号转导途径的深入研究奠定了基础。
4 结论
利用RNA-Seq技术获得了水杨酸(SA)诱导苹果叶片转录水平的差异基因,包括次生代谢、植物病原菌互作途径相关基因以及33个条件特异表达基因。利用原生质体转化技术鉴定了SA特异性响应基因WRKY40(MDP0000263349)的启动子序列,在1 000—1 500 bp及1 500—2 000 bp具有显著提高启动子对SA敏感性的未知核苷酸序列,另外,MdWRKY40转录调控存在反馈抑制机制,为苹果中SA信号途径的研究奠定了基础。The authors have declared that no competing interests exist.