欢迎访问林业科学,今天是

林业科学 ›› 2021, Vol. 57 ›› Issue (9): 110-120.doi: 10.11707/j.1001-7488.20210911

• 论文与研究报告 • 上一篇    下一篇

油茶炭疽病菌果生刺盘孢效应子的筛选

陈星州,周国英,陈行钢,江玲玉,包安华,刘君昂*   

  1. 中南林业科技大学 南方人工林病虫害防治国家林业局重点实验室 森林有害生物防控湖南省重点实验室经济林培育与保护教育部重点实验室 长沙 410004
  • 收稿日期:2020-09-01 出版日期:2021-09-25 发布日期:2021-11-29
  • 通讯作者: 刘君昂
  • 基金资助:
    国家自然科学基金项目(31971661);湖南省自然科学基金项目(2021JJ31145);湖南省研究生科研创新项目(CX20200712);中南林业科技大学研究生创新基金(X20201008)

Screening of Effectors of Colletotrichum fructicola in Camellia oleifera

Xingzhou Chen,Guoying Zhou,Xinggang Chen,Lingyu Jiang,Anhua Bao,Jun Liu*   

  1. Key Laboratory of Southern Plantation Pests and Diseases Prevention and Control of National Forestry and Grassland Administration Hunan Provincial Key Laboratory of Control of Forest Disease and Pests Key Laboratory for Non-Wood Forest Cultivation and Conservation of Ministry of Education Central South University of Forestry and Technology Changsha 410004
  • Received:2020-09-01 Online:2021-09-25 Published:2021-11-29
  • Contact: Jun Liu

摘要:

目的: 采用生物信息学研究方法,筛选分析果生刺盘孢效应子,并验证其功能,为其致病相关基因研究提供理论依据,同时为油茶病害管理提供参考。方法: 对果生刺盘孢非侵染期(分生孢子)和侵染时期(油茶叶片病斑组织)进行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 I型为主(占88.41%)。KEGG功能富集分析得到164条Pathway信息,多为真菌代谢途径、次生代谢产物、淀粉和蔗糖代谢相关的通路。PHI致病菌数据库同源比对结果表明,在经典分泌蛋白中含有80条与真菌致病力相关的基因,其中包括1个已知效应子。以氨基酸长度 < 300和半胱氨酸残基数≥4为筛选条件,并过滤已经注释过功能的蛋白后,在果生刺盘孢中筛选出17个候选效应子。2)随机抽样选取9个候选效应子做qRT-PCR分析,结果显示9个基因均为上调表达,与转录组数据结果相符。3)克隆出4个候选效应子,并通过烟草瞬时表达系统和DAB染色,验证出4个候选效应子均能引起植物细胞坏死。结论: 果生刺盘孢在侵染油茶时期有多个效应子表达,而在分生孢子时期此类蛋白多为不表达或少量表达,这验证了效应子是植物与病原菌相互作用的关键因素。

关键词: 油茶炭疽病, 果生刺盘孢, 转录组, 生物信息学, 效应子, 烟草瞬时表达, 分泌蛋白

Abstract:

Objective: Colletotrichum fructicola is the pathogen of anthracnose in Camellia oleifera. In this study, bioinformatics research methods were used to screen and analyze the effectors of Colletotrichum fructicola, and their functions were verified, so as to provide a theoretical basis for the research of disease-related genes of Colletotrichum fructicola, and to provide a reference for disease management in Camellia oleifera cultivation. Method: RNA-seq sequencing was performed at the non-infection stage (conidia) and infection stage (lesion tissue of Camellia oleifera leaves) of C. fruiticola. A series of bioinformatics software, such as BUSCA, target P, big PI predictor, and databases of GO, KEGG and PHI, were used to screen out the candidate effectors and then the effectors were quantitatively analyzed with qRT-PCR. The full-length gene of the candidate effector was cloned, and the corresponding gene function was verified by tobacco transient expression system and DAB staining. Result: The transcriptome data showed that there were 7 850 differentially expressed genes at the infection stage compared to the conidia stage. Predictive analysis of the encoded proteins of differentially expressed genes revealed that 345 proteins were classical secretory proteins, accounting for 4.39% of the total. The amino acid length of the classical secretory protein was evenly distributed in each segment, and the protein that meets the effector screening conditions (< 300 aa) accounted for about 36% of the total protein, which was relatively small. The length of the signal peptide was concentrated in 18-20 aa (48.7%), and the signal peptide cleavage site was SPase I type(88.41%). A total of 164 pieces of pathway information were gained by KEGG functional enrichment analysis, most of which were fungal metabolic pathways, secondary metabolites, and pathways related to starch and sucrose metabolism. The homology comparison result of PHI pathogen database showed that there were 80 genes related to fungal pathogenicity in the classical secretory protein, including a known effector. With the length of amino acid < 300 and cysteine residue number ≥ 4 as screening conditions, 17 candidate effector were screened out in C. fructicola after filtering the proteins with annotated functions. Nine candidate effectors were randomly selected for qRT-PCR analysis, and the results showed that the expression of all nine genes were up-regulated, which was in accordance with the results of transcriptome data. Four candidate effectors were cloned and verified by tobacco transient expression system and DAB staining, indicating that all candidate effectors was able to lead to plant cell necrosis. Conclusion: Multiple effectors are expressed during the infection period of Camellia oleifera by C. fructicola, while most of these proteins are not expressed or slightly expressed during the conidia period. This proves that the effectors is the key factor in the interaction between plant and pathogen.

Key words: Camellia oleifera anthracnose, Colletotrichum fructicola, transcriptome, bioinformatics, effector, Tobacco transient expression, secretory protein

中图分类号: