论文部分内容阅读
甘蔗(Saccharum spp.)是我国最重要的糖料作物和能源作物,蔗糖占我国食糖总产的92%。由黑粉菌(Sporisorium scitamineum, S. scitamineum)引起的甘蔗黑穗病(sugarcane smut),严重危害甘蔗生产,导致产量和蔗糖分降低,给我国甘蔗产业造成了巨大的经济损失。培育抗黑穗病甘蔗品种是控制甘蔗黑穗病危害的最经济有效的手段。植物与病原互作的分子机制包括转录水平、转录后水平、翻译水平、翻译后水平等多方面调控机制。随着生物技术的发展,有关甘蔗与黑穗病菌互作分子机制的研究有了长足进步,但是尚未见利用高通量测序技术鉴定甘蔗受黑穗病菌胁迫后microRNA (miRNA)差异表达及其靶标基因预测和功能分析的报道。本研究首次利用高通量测序技术鉴定甘蔗受黑穗病菌胁迫下miRNA的差异表达,从甘蔗小RNA高通量测序与niRNA鉴定注释、差异miRNA筛选和实时荧光定量PCR技术(qRT-PCR)表达验证与表达模式分析、差异miRNA的功能分析及其靶基因预测和靶基因表达分析等三个方面入手,系统比较了甘蔗与黑穗病菌互作系统的亲和与不亲和互作间在转录后水平上miRNA的调控机制,以期为甘蔗抗黑穗病分子机理的研究提供niRNA水平的佐证,丰富和深化甘蔗抗黑穗病的分子机理研究,进而为培育抗黑穗病甘蔗品种提供miRNA水平的理论依据。miRNA具有重要的生物学功能,在植物受到生物和非生物胁迫的过程中发挥重要调控作用,其介导的转录后水平的基因调控是生物体内一条重要的小RNA调控路径。本研究以甘蔗抗黑穗病无性系崖城05-179 (YA05-179)和感黑穗病甘蔗品种新台糖22 (ROC22)为研究材料,以接种黑穗病菌蔗芽作为处理组(YAT、RT),接种清水作为对照组(YACK、RCK),应用高通量HiSeq测序技术鉴定甘蔗受黑穗病胁迫下48h的miRNA差异表达,对所获得的小RNA序列注释并筛选出差异表达的miRNA,利用生物信息学技术与数据库结合对其靶向基因进行预测。实验对筛选出的部分差异miRNA进行实时定量验证,并进一步分析了甘蔗与黑穗病菌互作Oh、12h、48h和96h不同时间点miRNA与其部分靶基因的表达模式。研究旨在找到甘蔗在黑穗病菌胁迫下发挥关键调控作用的miRNA及其靶基因,为甘蔗抗黑穗病分子育种提供miRNA和基因资源。主要研究结果如下:1.四个甘蔗样品的送样质量良好,分别测序得到符合要求的干净序列为RCK 36396588条、RT27812972条、YACK 27464468条、YAT 28290231条,且杂质序列含量都少于1%,故高通量测序获得的小RNA质量高。结合多个数据库对测得的sRNAs进行筛选注释,其中RCK测序获得264个已知miRNA,RT有263个,YACK有260个,YAT有262个。此外,在RCK、RT、YACK和YAT中分别预测出137、140、111和119个新miRNA。有碱基分布统计分析结果显示,四个小RNA文库中在植物中含量较多的长度为21 nt和22 nt的新niRNA序列首位点碱基U所占比例较大;在序列碱基11位点,四个样品的新miRNA序列的碱基都偏向于A,符合miRNA序列特征,说明新miRNA序列的预测结果可行性高。2.在ROC22中,筛选出231个已知miRNA,其中35个在黑穗病菌胁迫处理样品和清水接种对照样品间发生显著差异表达(|log2一ratio|>1,P<0.05);涉及10个上调表达和25个下调表达的miRNA:在对照组与处理组样品中表达量差异倍数最高达4.6倍。在YA05-179中,两个样品中共有的已知miRNA为208个,对照组与处理组中发生显著差异表达niRNA有11个;涉及2个上调表达和9个下调表达的miRNA:对照组与处理组样品中表达量差异倍数最高达2.1倍。在ROC22和YA05-179中,筛选出的差异表达的已知miRNA,下调表达的占多数,且大部分miRNA的种类和表达趋势是相同的。在ROC22中,筛选出的显著差异表达(|log2-ratio|>1,P<0.05)的新miRNA有16个,而YA05-179中有6个;其在单一品种特有的比较多,且在两个品种中没有发现共同显著差异表达的新miRNA。本研究针对筛选出的16个差异表达的已知miRNA和新miRNA,利用qRT-PCR对其在ROC22与YA05-179两个品种中27个表达结果进行测序结果验证,相符率达到了85.2%,可见本研究的测序结果可靠性较高。3.甘蔗是个高度杂合的异源多倍非整倍体作物,其全基因组测序尚未完成,因此本研究注释的miRNA未知序列占小RNA序列的绝大多数,随后对所注释的miRNA在甘蔗受黑穗病菌胁迫时对应靶基因的预测及分析数据也不完整。差异表达的已知miRNA预测出的靶基因数,YA05-179少于ROC22,但差异表达的新miRNA对应的靶基因数目在YA05-179中明显多于R022。靶基因生物功能分析中,GO归类为生物过程(biological process)、细胞组分和元件(cellular components)分子功能(molecular function)。两个品种中,差异表达的已知miRNA靶基因的分类情况与差异表达的新miRNA所匹配的靶基因序列统计情况相似,即在生物过程分类中靶基因主要分布在细胞过程与代谢过程,细胞组分和元件分类中靶基因主要分布在细胞、细胞组分和细胞器中,分子功能分类中靶基因主要分布在结合与催化活性。KEGG通路富集分析显示,差异miRNA调控的靶基因参与了一系列miRNA调控靶基因生理生化途径或抗病相关生理代谢和信号转导途径,表明甘蔗与黑穗病菌互作的分子机制是一个复杂网络调控系统。研究发现,miR5077、 miR5671、miR5044、miR5261、miR5783、miR5221、miR6478和miR948的靶向基因在MAPK信号通路(MAPK signaling pathway)、植物激素信号转导(Plant hormone signal transduction)和植物-病原菌互作 (Plant-pathogen interaction)中发挥重要作用并构建了相关蛋白网络调控图。4.从miRNA靶基因中筛选了在甘蔗EST文库中有生物功能注释的序列,并针对所筛选出的差异表达的已知miRNA和新miRNA所对应的部分靶基因进行了qRT-PCR表达分析。结果显示,在已知miRNA与其靶基因的表达分析中,大部分miRNA与其靶基因呈负调控模式,即niRNA过表达时,靶基因的表达量呈下降趋势;其中8个靶基因中有7个靶基因的表达在12h后基本相符,且在48h时匹配度最高。这也从另一方面说明甘蔗受黑穗病菌胁迫后miRNA对靶基因的调控作用在48h达到最高。研究还揭示,本研究中已知miRNA与靶基因表达丰度预测与功能分析的可靠性高,但在候选的新miRNA及其靶基因定量表达结果之间的负调控趋势的符合率就没有那么高,故还需对新miRN A的预测与分析作进一步的挖掘。