摘要

【目的】采用生物信息学研究方法,筛选分析果生刺盘孢效应子,并验证其功能,为其致病相关基因研究提供理论依据,同时为油茶病害管理提供参考。【方法】对果生刺盘孢非侵染期(分生孢子)和侵染时期(油茶叶片病斑组织)进行RNA-seq测序,并利用生物信息分析软件BUSCA、Target P、big-PI Predictor和GO、KEGG、PHI数据库分析转录组数据,选出符合条件的候选效应子,再利用qRT-PCR技术对效应子进行相对定量分析;克隆候选效应子全长基因,通过烟草瞬时表达系统和DAB染色验证相应基因功能。【结果】1)转录组结果数据显示,侵染时期相对于非侵染时期有7 850个差异表达基因。对差异表达基因的编码蛋白进行预测分析,345个为经典分泌蛋白,占总数的4.39%。经典分泌蛋白氨基酸长度在各区间段分布较为均匀,而符合效应子筛选条件(<300 aa)的蛋白约占经典分泌蛋白总数的36%,数量相对较少。信号肽长度集中在18~20 aa(占48.7%),信号肽切割位点以SPase Ⅰ型为主(占88.41%)。KEGG功能富集分析得到164条Pathway信息,多为真菌代谢途径、次生代谢产物、淀粉和蔗糖代谢相关的通路。PHI致病菌数据库同源比对结果表明,在经典分泌蛋白中含有80条与真菌致病力相关的基因,其中包括1个已知效应子。以氨基酸长度<300和半胱氨酸残基数≥4为筛选条件,并过滤已经注释过功能的蛋白后,在果生刺盘孢中筛选出17个候选效应子。2)随机抽样选取9个候选效应子做qRTPCR分析,结果显示9个基因均为上调表达,与转录组数据结果相符。3)克隆出4个候选效应子,并通过烟草瞬时表达系统和DAB染色,验证出4个候选效应子均能引起植物细胞坏死。【结论】果生刺盘孢在侵染油茶时期有多个效应子表达,而在分生孢子时期此类蛋白多为不表达或少量表达,这验证了效应子是植物与病原菌相互作用的关键因素。

全文