-
MicroRNAs(miRNA)是一类由20~22个核苷酸构成的非编码小RNA,可通过结合靶基因mRNA的3′非翻译区(3′-untranslated region,3′-UTR)使其降解或者阻碍其转录,从而在转录后水平控制基因表达[1]。作为基因表达负调控子,miRNA主要在转录后水平调节植物基因表达[2],在植物新陈代谢、组织生长、器官发育分化以及细胞程序性凋亡中起重要作用[3-4]。此外,miRNA在植物感受逆境胁迫而做出适应性调整的过程中发挥着重要的作用。特别是在病原侵染过程中,miRNA能够调控植物抗病途径中的关键靶基因来决定病原与寄主的亲和关系(compatible/incompatible interaction)。
MYB(v-myb avian myeloblastosis viral oncogene homolog)是一类在其N端含有MYB结构域的转录因子(transcription factor,TF)[5]。MYB是植物中最大的一类转录因子家族,涉及水杨酸(salicylic acid,SA)、脱落酸(abscisic acid,ABA)及茉莉酸(jasmonic acid,JA)等植物抗病信号传导通路。拟南芥(Arabidopsis thalian)MYB30编码通过调控超长链(very-long-chain)脂肪酸的合成应答病原入侵[6]。MYB44经由WRKY70激活SA介导的防卫反应,抑制JA介导的防卫反应,通过拮抗调节方式调控拟南芥对生物胁迫应答[7]。水稻OsLTR1编码MYB转录因子调控JA防御反应[8]。根据蛋白含有结构域的数量,MYB转录因子分为4种类型,其中R2R3-MYB转录因子参与植物响应生物和非生物胁迫等过程[9]。目前发现的植物MYB在非生物胁迫表达过程主要受miR159、miR828、miR172和miR319等家族的miRNA调控。例如拟南芥ABA刺激促进miR159积累,miR159过表达则抑制MYB101和MYB33转录[10]。
对落叶松-杨栅锈菌(Melampsora larici-populina)E4生理小种侵染欧美杨(Populus × euramericana)研究过程中[11],我们发现欧美杨抗病相关基因的表达并未发生预想的显著上调。以NBS-LRR基因为例,E4生理小种侵染后仅8个TIR-NBS-LRR表达上调,另有11个TIR-NBS-LRR和3个CC-NBS-LRR基因表达下调,其余21个NBS-LRR基因在E4生理小种侵染后表达量未发生显著变化。但,与抗病相关的NAC、WRKY和MYB等转录因子基因表达量却发生显著变化。大量研究表明miRNA参与植物生长发育和胁迫应答中的各个方面,可以在不同层面对靶基因进行调控,是基因表达精细调控必不可少的重要手段,处于复杂调控网络中的关键位置。关于木本植物miRNA在病原菌胁迫中的调控机制的研究相对较少,特别是“杨树-锈菌”体系的miRNA调控机制研究仍处于探索阶段[11, 12]。为揭示欧美杨对E4锈菌感病的分子调控机理,本研究通过miRNA组、降解组(degradome)和基因表达谱(digital gene expression profiling,DGE)等高通量测序技术对锈菌侵染前后欧美杨R2R3-MYB基因及其调控miRNA进行研究。
HTML
-
欧美杨被E4强致病性锈菌侵染至7 dpi时叶组织坏死,无过敏反应(hypersensitive response,HR)症状产生,且已产生大量夏孢子堆(图 1 Rust+),ddH2O对照处理组无明显变化(图 1 Rust-)。Rust+组叶片在2 hpi~1 dpi未表现明显症状,2 dpi叶片局部出现细小褪绿斑点,4 dpi叶片可见明显黄色斑点,7 dpi病斑面积扩大,锈菌完成一个生长周期开始形成夏孢子堆。发病过程中未见HR和细胞程序化死亡(programmed cell death,PCD)等抗性症状出现,表明欧美杨为高感病树种,与E4强致病性锈菌构成亲和型互做关系。
-
DGE测序共鉴定到欧美杨230个R2R3-MYB基因,其中log2(rust+/rust-)表达量有差异的基因55个(36个锈菌接种后表达量上调,19个表达量下调),175个表达量未发生显著变化。与毛果杨基因组相比较,有36个毛果杨R2R3-MYBs未在DGE测序中找到同源基因。此外,DGE测序检测到208个MYB相关基因(myb-like gene),其中21个表达量下调、3个表达量下调、其余表达量未发生显著变化。进一步分析表达发生变化的R2R3-MYB基因,有10个锈菌侵染前不表达的R2R3-MYBs在锈菌侵染后表达,另有9个锈菌侵染前表达的R2R3-MYBs在锈菌侵染后停止表达(图 2)。表明,锈菌侵染不但可以调控R2R3-MYBs的表达量,还可以启动或关闭R2R3-MYBs的表达。
-
综合分析miRNA组和降解组数据,鉴定到18个R2R3-MYBs与22个miRNA共组成86对转录后调控关系(图 3)。仅有7对转录后调控是单一miRNA对单一R2R3-MYB,其余均为1个miRNA调控多个R2R3-MYBs且每个R2R3-MYB受2个以上(包括2个)miRNA调控。调控miRNA主要集中在miR159(6个)和miR858(5个)2个家族,同时检测到PC-3p-2521022_1和PC-5p-301019_3(图 3黄色标识)2条新miRNA。
-
鉴于鉴定到的miRNA与R2R3-MYB转录后调控关系比较复杂,本研究暂选择一对一调控关系的新miRNA调控关系进行荧光定量PCR验证。PC-3p-2521022_1在锈菌接种后7个时间点(2 hpi、6 hpi、12 hpi、24 hpi、48 hpi、96 hpi和168 hpi)表达量变化处于-0.5~0.5之间,变化幅度较小,在12~96 hpi出现往复式波动(图 4柱状图)。PndMYB173表达量变化处于-2~3之间,变化幅度较大,在12~96 hpi同样出现往复式波动,且与其调控miRNA PC-3p-2521022_1呈负调控关系(图 4折线图)。两配对样本t检验结果表明有PC-3p-2521022_1表达量有5个相邻时间点差异极显著(n=3),PndMYB173表达量有4个相邻时间点差异极显著(n=3)。
-
根据结合域,55个差异表达R2R3-MYBs中的48个成功预测到共计336个靶基因。依据STRING数据库预测48个差异表达R2R3-MYBs与靶基因构成400个互做关系(图 5)。其中PndMYBx30具有最多的139个靶基因,构成139个互做关系。根据亚细胞定位,这些靶基因主要定位于细胞核(nucleus)、质膜(plasma membrane)、高尔基体(golgi)、液泡(vacuole)和叶绿体(chloroplast)等细胞器。GO富集分析表明139个靶基因涉及多个与免疫、植物激素相应和抗逆相关的生物过程(biological process,表 1)。
序号
No.GO术语
GO annotation功能描述
Description基因数
Gen numberp值
p value伪发现率
FDR1 GO: 0045087 先天免疫应答innate immune response 5 0.00035 0.0021 2 GO: 0002376 免疫系统进程immune system process 5 0.00046 0.0027 3 GO: 0006955 免疫应答immune response 5 0.00046 0.0027 4 GO: 0050896 响应刺激response to stimulus 34 1.40E-16 1.40E-15 5 GO: 0042221 响应化学刺激response to chemical stimulus 31 6.40E-22 7.90E-21 6 GO: 0010033 响应有机物response to organic substance 29 8.80E-25 1.20E-23 7 GO: 0009719 响应内源性刺激response to endogenous stimulus 28 4.30E-26 6.70E-25 8 GO: 0009725 响应激素刺激response to hormone stimulus 26 3.20E-24 4.30E-23 9 GO: 0005634 细胞核nucleus 21 1.50E-09 3.80E-08 10 GO: 0009751 响应水杨酸刺激response to salicylic acid stimulus 19 1.40E-27 2.20E-26 11 GO: 0006950 应激响应response to stress 19 7.90E-09 6.90E-08 12 GO: 0009753 响应茉莉酸刺激response to jasmonic acid stimulus 18 3.20E-25 4.60E-24 13 GO: 0009737 响应脱落酸刺激response to abscisic acid stimulus 18 4.90E-21 5.90E-20 14 GO: 0009628 响应非生物刺激response to abiotic stimulus 17 3.70E-10 3.30E-09 15 GO: 0006970 响应渗透压刺激response to osmotic stress 15 5.20E-16 4.80E-15 16 GO: 0009739 响应赤霉素刺激response to gibberellin stimulus 14 5.10E-20 5.70E-19 17 GO: 0009723 响应乙烯刺激response to ethylene stimulus 14 9.90E-19 1.10E-17 18 GO: 0010035 响应无机物response to inorganic substance 14 8.80E-17 8.90E-16 19 GO: 0009733 响应生长素刺激response to auxin stimulus 13 7.00E-14 6.30E-13 Table 1. GO enrichment analysis of target genes