-
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进行研究。
-
欧美杨被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)。
图 5 欧美杨R2R3-MYBs靶基因网络
Figure 5. Populus × euramericana regulatory networks of R2R3-MYBs to target genes
表 1 靶基因GO富集分析
Table 1. GO enrichment analysis of target genes
序号
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
MicroRNA转录后调控欧美杨R2R3-MYBs抗锈菌表达
MicroRNA-regulated Compatible Poplar R2R3-MYB Gene Expression Following the Infection of Virulent Melampsora larici-populina
-
摘要:
目的 强致病锈菌亲和型树种欧美杨适合做寄主感病分子机理研究。研究表明转录因子及其调控miRNA在杨树抗/感病过程中起重要的信号传导和调控作用,决定杨树与锈菌的亲和性。为深入研究杨树感病性,对欧美杨R2R3-MYB转录因子及其调控miRNA进行研究。 方法 构建锈菌侵染和未侵染miRNA组、降解组文库和基因表达谱文库,进行高通量测序。对基因表达谱结果进行生物信息学分析鉴定R2R3-MYBs,根据miRNA组和降解组数据确定R2R3-MYBs调控miRNA。应用定量PCR技术鉴定候选R2R3-MYBs和其调控miRNA的表达量。 结果 共鉴定到230个欧美杨R2R3-MYBs,其中55个R2R3-MYBs在锈菌侵染和未侵染叶片间表达量差异显著,其余表达量变化不显著。基于降解组,共鉴定到由22个miRNA调控18个R2R3-MYBs的86个转录后调控关系。定量PCR结果表明PC-3p-2521022_1与PndMYB173存在转录后负调控关系。预测的R2R3-MYB靶基因涉及水杨酸、茉莉酸、乙烯和脱落酸等多个植物激素抗性路径。 结论 E4强致病锈菌侵染导致亲和型树种欧美杨23.9%的R2R3-MYBs表达量发生变化,锈菌侵染不但可以影响R2R3-MYBs的表达量,还可以启动或关闭R2R3-MYBs的表达。欧美杨在E4强致病锈菌胁迫条件下,R2R3-MYBs转录后主要受miR159和miR858家族的miRNA调控。 Abstract:Objective The poplar 'Robusta', susceptible to virulent rust E4 Melampsora larici-populina, provides suitable host material for the study of the plant host resistance/susceptibility to pathogen. Previously study indicates that the transcription factor (TF) plays an important role in determining the poplar susceptibility. In this study, the authors tried to determine the R2R3-MYB TF genes driving basal disease resistance and focus on microRNAs (miRNAs) regulation. Method Two small RNA (sRNA) libraries, a degradome cDNA library and 2 digital gene expression libraries were constructed for rust-inoculated (Rust+) and rust-free (Rust-) susceptible poplar 'Robusta' leaves through high-throughput sequencing. The expression of selected miRNA and target R2R3-MYBs were identified by RTq-PCR. Result Altogether, 230 R2R3-MYBs were identified. Among those R2R3-MYBs, the expression of 55 genes were suppressed/promoted following the drive straight of the rust, the others were unchanged. Based on the miRNA and degradome sequencing, 86 regulating interactions (22 miRNA to 18 R2R3-MYB) were identified. RTq-PCR of selected PC-3p-2521022_1 and PndMYB173 showed a negative post-transcriptional regulation. The R2R3-MYBs target genes predication indicated R2R3-MYBs related to many resistance and phytohormone pathways, such as response to salicylic acid stimulus, response to stress, response to jasmonic acid stimulus and response to abscisic acid stimulus. Conclusion The infection of virulent rust changed the expression of 23.9% R2R3-MYBs. The differential expression of R2R3-MYBs indicated that the rust fungi infection could not only switch on but also switch off R2R3-MYB TF genes. Under the stress of rust infection, the post-transcriptional regulation of R2R3-MYBs was regulated by the families of miR159 and miR858. -
Key words:
- microRNA
- / R2R3-MYB
- / Populus× euramericana
- / Melampsora larici-populina
- / compatible interaction
-
表 1 靶基因GO富集分析
Table 1. GO enrichment analysis of target genes
序号
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 -
[1] Bartel D P. MicroRNAs: target recognition and regulatory functions[J]. Cell, 2009, 136: 215-233. doi: 10.1016/j.cell.2009.01.002 [2] Schwab R, Palatnik JF, Riester M, et al. Specific effects of microRNAs on the plant transcriptome[J]. Developmental Cell, 2005, 8: 517-527. doi: 10.1016/j.devcel.2005.01.018 [3] Jover-Gil S, Candela H, Ponce MR. Plant microRNAs and development[J]. The International Journal of Developmental Biology, 2005, 49: 733-744. doi: 10.1387/ijdb.052015sj [4] Willmann M R, Poethig RS. Conservation and evolution of miRNA regulatory programs in plant development[J]. Current Opinion in Plant Biology, 2007, 10: 503-511. doi: 10.1016/j.pbi.2007.07.004 [5] 刘蕾, 杜海, 唐晓凤, 等. MYB转录因子在植物抗逆胁迫中的作用及其分子机理[J]. 遗传, 2008, 30(10): 1265-1271. doi: 10.3321/j.issn:0253-9772.2008.10.005 [6] Li L. Arabidopsis MYB30 is a direct target of BES1 and cooperates with BES1 to regulate brassinosteroid-induced gene expression[J]. Plant Journal, 2009, 58: 275-286. doi: 10.1111/j.1365-313X.2008.03778.x [7] Shim J S, Jung C, Lee S, et al. AtMYB44 regulates WRKY70 expression and modulates antagonistic interaction between salicylic acid and jasmonic acid signaling[J]. Plant Journal, 2012, 73(3): 483-495. [8] Johnson E T, Dowd P F. Differentially enhanced insect resistance, at a cost, in Arabidopsis thaliana constitutively expressing a transcription factor of defensive metabolites[J]. J Agric Food Chem, 2004, 52: 5135-5138. doi: 10.1021/jf0308049 [9] 牛春阳, 王峰, 李丹蕾, 等. C14族R2R3-MYB基因调控杨树抗锈菌过敏性反应[J]. 北京林业大学学报, 2016, 38(7): 16-23. [10] Reyes J L, Chua N H. ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination[J]. Plant Journal, 2007, 49(4): 592-606. doi: 10.1111/j.1365-313X.2006.02980.x [11] Li Danlei, Feng Wang, Chao Wang, et al. MicroRNA-mediated susceptible poplar gene expression regulation associated with the infection of virulent Melampsora larici-populina[J]. Bmc Genomics, 2016, 17(1): 1-16. [12] Min C, Cao Z. Genome-wide expression profiling of microRNAs in poplar upon infection with the foliar rust fungus Melampsora larici-populina[J]. Bmc Genomics, 2015, 16(1): 173-184. doi: 10.1186/s12864-015-1362-2 [13] 陈俏丽, 李丹蕾, 王峰, 等. 欧美杂交杨Pnd-LRR3基因克隆及其抗锈菌侵染表达[J]. 东北林业大学学报, 2015, 3: 95-98. [14] By Yoav B, Daniel Y. The control of the false discovery rate in multiple testing under dependency[J]. Annals Statistics. 2001; 29: 1165-88. [15] Li Q, Lin Y C, Sun Y H, et al. Splice variant of the SND1 transcription factor is a dominant negative of SND1 members and their regulation in Populus trichocarpa[J]. Proceedings of the National Academy of Sciences, 2012, 109(36): 14699-14704. doi: 10.1073/pnas.1212977109 [16] Duan Y, Jiang Y, Ye S, et al. PtrWRKY73, a salicylic acid-inducible poplar WRKY transcription factor, is involved in disease resistance in Arabidopsis thaliana[J]. Plant Cell Reports, 2015, 34(5): 831-841. doi: 10.1007/s00299-015-1745-5 [17] Nawrath C, Metraux J P. Salicylic acid induction-deficient mutants of Arabidopsis express PR-2 and PR-5 and accumulate high levels of camalexin. Plant Cell. 1999; 11: 1393-1404. [18] Navarro L, Dunoyer P, Jay F, et al. A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science. 2006; 312: 436-439. doi: 10.1126/science.1126088 [19] Liu P T, Wheelwright M, Teles R, et al. MicroRNA-21 targets the vitamin D-dependent antimicrobial pathway in leprosy. Nat Med. 2012; 18: 267-273. doi: 10.1038/nm.2584