免费文献传递   相关文献

The comparative analysis of secreted CAZymes in phytopathogenic fungi with different lifestyle

植物病原真菌寄生性与分泌蛋白组CAZymes的比较分析



全 文 :植物病理学报
ACTA PHYTOPATHOLOGICA SINICA  44(2): 163 ̄172(2014)
收稿日期: 2013 ̄06 ̄04ꎻ 修回日期: 2013 ̄10 ̄09
基金项目: 国家重点基础研究发展计划(2013CB127803)ꎻ国家自然科学基金(31171799ꎻ 31200113)
通讯作者: 戴小枫ꎬ研究员ꎬTel: 010 ̄62815976ꎻ E ̄mail: daixiaofeng@caas.cn
第一作者: 陈相永ꎬ男ꎬ安徽萧县人ꎬ博士研究生ꎬ研究方向为植物病原菌基因组学ꎻE ̄mail: xy_ch@163.comꎮ
植物病原真菌寄生性与分泌蛋白组
CAZymes的比较分析
陈相永ꎬ 陈捷胤ꎬ 肖红利ꎬ 桂月靖ꎬ 李 蕾ꎬ 戴小枫∗
(中国农业科学院农产品加工研究所ꎬ北京 100193)
摘要:通过信息学方法对植物病原真菌分泌蛋白组进行预测以及对碳水化合物酶类 CAZymes注释和比较ꎬ分析病原真菌寄
生性与 CAZymes家族功能的关系ꎮ 结果显示非专性寄生(半活体营养型和死体营养型)真菌编码的分泌蛋白占基因组编码
基因的比例较专性寄生真菌(活体营养型)的偏高ꎻCAZymes家族中ꎬ非专性寄生真菌的糖基水解酶家族 GH和多糖裂解酶家
族 PL较专性寄生真菌显著扩增ꎻ功能聚类分析发现ꎬ非专性寄生真菌参与纤维素、果胶、木聚糖等植物细胞壁组分降解相关
的基因较专性寄生真菌明显增多ꎮ 结果表明了 CAZymes与寄生类型的高度关联性ꎬ揭示了其在非专性寄生真菌致病过程中
可能的关键作用ꎮ
关键词:植物病原真菌ꎻ寄生性ꎻ分泌蛋白组ꎻ碳水化合物酶类ꎻ植物细胞壁降解
The comparative analysis of secreted CAZymes in phytopathogenic fungi with diffe ̄
rent lifestyle  CHEN Xiang ̄yongꎬ CHEN Jie ̄yinꎬ XIAO Hong ̄liꎬ GUI Yue ̄jingꎬ LI Leiꎬ DAI Xiao ̄feng
(Institute of Agro ̄Products Processing Science and Technologyꎬ Chinese Academy of Agricultural Sciencesꎬ Beijing 100193ꎬ
China)
Abstract: To understand the correlation between lifestyle and extracellular CAZymesꎬ secretomes and
CAZymes were predicted and annotated in phytopathogenic fungiꎬ followed by comparative analysis. The results
indicated that the ratio of secreted proteins to total putative proteins in hemibiotrophic and necrotrophic fungi was
higher than that in biotrophic fungi. The CAZymes especially glycoside hydrolases (GH) and polysaccharide
lyases ( PL) subfamilies were expanded significantly in hemibiotrophic and necrotrophic fungi. Functional
analysis of secreted CAZymes indicated that the plant cell wall ( PCW) degrading enzymes involved in
degradation of celluloseꎬ pectin and xylanꎬ were more abundant in hemibiotrophic and necrotrophic fungi. The
distinct expansion of PCW degrading enzymes indicated the high correlation between fungal lifestyle and
CAZymesꎬ meanwhile revealing their potential roles in pathogenesis in hemibiotrophic and necrotrophic fungi.
Key words: phytopathogenic fungiꎻ pathogenic lifestyleꎻ secretomeꎻ CAZymesꎻ plant cell wall degradation
中图分类号: S432.4          文献标识码: A          文章编号: 0412 ̄0914(2014)02 ̄0163 ̄10
    植物病原真菌依据寄生的程度可以分为专性寄
生物和非专性寄生物ꎮ 专性寄生物的营养方式为活
体营养型(biotroph)ꎻ非专性寄生物除了寄生生活
外ꎬ还可以在死的植物组织上生活ꎬ营养方式包括半
活体营养型(hemibiotroph)和死体营养型(necrotro ̄
ph) [1]ꎮ 病原菌不同寄生性反映了其对寄主植物侵
染策略或致病机制的多样性ꎬ而对寄主植物成功侵
染的前提是侵入的病原菌必须适应寄主植物并从中
获取生长、发育以及繁殖所必须的各种营养物质ꎬ在
这个过程中ꎬ与寄主植物直接互作的胞外蛋白具有
 
植物病理学报 44卷
重要的作用[2~4]ꎮ
    胞外分泌蛋白是病原菌在细胞内合成后通过各
种分泌途径转运到细胞表面和质膜空间感知周围环
境并做出各种适应性反应的蛋白ꎬ它们可与寄主细胞
中的相关蛋白互作引发寄主一系列与病害发展相关
的生理生化反应并最终形成病害ꎮ 分泌的胞外蛋白
除了作为效应物参与侵染寄主过程外ꎬ还可以作为毒
力因子或者毒素改变寄主细胞的结构或者功能ꎬ或者
作为无毒因子或者激发子引发寄主防卫反应[2ꎬ 5]ꎮ 还
存在大量的碳水化合物酶类(carbohydrate active en ̄
ZYmesꎬ CAZymes)参与了寄主植物细胞壁的降解作
用ꎬ包括酯酶、纤维素酶、果胶酶、葡聚糖酶、木葡聚糖
酶等ꎬ促进病原菌在寄主植物细胞中的附着、侵入、定
殖以及营养吸收等过程ꎮ 因此ꎬ胞外分泌蛋白中的碳
水化合物酶类尤其是植物细胞壁降解酶对病原菌的
侵染能力和特性具有重要作用[6]ꎮ
    随着大量植物病原菌基因组测序的完成ꎬ通过
比较基因组学分析病原菌的寄生性与基因组变异
的关系已经成为研究热点[7]ꎬ尤其是胞外分泌蛋
白 CAZymesꎮ CAZymes 参与碳水化合物的合成
与降解ꎬ包括 5 大类:(1)糖基水解酶( glycoside
hydrolasesꎬGH) [8]ꎻ(2)糖基转移酶(glycosyltrans ̄
ferasesꎬGT) [9]ꎻ ( 3) 多糖裂解酶 ( polysaccharide
lyasesꎬPL) [10]ꎻ(4)碳水化合物酯酶( carbohydrate
esterasesꎬCE) [10]ꎻ(5)碳水化合物绑定结构( car ̄
bohydrate binding modulesꎬCBM) [11]ꎮ 不同种类
CAZymes都包含了较多的亚家族ꎬ针对不同结构
的碳水化合物行使相关催化功能ꎬ其中部分亚家族
直接参与了寄主植物细胞壁组分(纤维素、果胶、
木聚糖等) 的降解ꎬ如 GH28、 GH43、 PL1、 PL3、
PL9、CE8、CE12等ꎬ以利于病原菌对寄主植物的成
功侵染[12ꎬ 13]ꎮ 已有研究表明非活体营养型病原菌
基因组较活体营养型编码更多的碳水化合物水解
酶[14~16]ꎮ 通过黄萎病菌 ( Verticillium dahliae 和
V. albo-atrum)和稻瘟病菌(Magnaporthe oryzae)
分泌蛋白组的 CAZymes 比较分析发现ꎬ黄萎病菌
编码的碳水化合物结合模块 CBM1(carbohydrate-
binding modules 1)亚家族明显扩增ꎬ推测其与黄萎
病菌死体营养型的广谱寄生相关[13]ꎮ 因此ꎬ在基
因组层面上分析比较分泌到胞外的 CAZymes对理
解病原菌寄生性差异和侵染特性有重要意义ꎮ
    本研究以不同寄生性的 12个植物病原真菌为
研究对象ꎬ包括 3 个活体营养型、5 个死体营养型
和 4个半活体营养型ꎬ通过全基因组编码蛋白的分
泌特性预测获得分泌蛋白组ꎬ对分泌蛋白组的
CAZymes进行注释和比较分析ꎬ进而阐明病原菌
编码的植物细胞壁降解酶类差异与其寄生性(营
养型)的关系ꎮ
1  材料与方法
1.1  植物病原菌及其基因组数据 
    本研究对 12种不同寄生性的植物病原真菌分泌
蛋白组的细胞壁降解酶类进行了比较分析ꎬ包括 3种
活体营养型真菌:玉米黑粉病菌(Ustilago maydis)、玉
米丝黑穗病菌(Sporisorium reilianum)和杨树叶锈病
菌(Melampsora laricis ̄populina)ꎻ4种半活体营养型
真菌:稻瘟病菌(Magnaporthe grisea)、小麦白斑病菌
(Mycosphaerella graminicola)、禾谷镰孢菌(Fusarium
graminearum)和禾生炭疽菌(Colletotrichum gramini ̄
cola)ꎻ5 种死体营养型真菌:玉米小斑病菌(Coch ̄
liobolus heterotrophus)、小麦叶斑病菌(Stagonospora
nodorum)、核盘菌(Sclerotinia sclerotiorum)、灰霉菌
(Botrytis cinerea)和大丽轮枝菌(Verticillium dahli ̄
ae)ꎮ 其中玉米丝黑穗病菌基因组编码蛋白数据下载
自 NCBI 数据库(http: / / www􀆰 ncbi􀆰 nlm􀆰 nih􀆰 gov / ge ̄
nome)ꎻ灰霉菌、禾生炭疽菌、玉米黑粉病菌、稻瘟病
菌、核盘菌、禾谷廉孢菌和大丽轮枝菌基因组编码蛋
白序列下载自美国布罗德研究所真菌基因组数据库
(Broad Instituteꎬhttp: / / www􀆰 broadinstitute􀆰 org / / sci ̄
entific ̄community / data)ꎻ其他真菌基因组编码蛋白下
载自联合基因组研究所真菌基因组数据库(Joint Ge ̄
nome Instituteꎬhttp: / / genome􀆰 jgi􀆰 doe􀆰 gov / programs /
fungi / index.jsfꎮ
1.2  分泌蛋白预测
    分泌蛋白预测按如下流程进行:利用 WoLF
PSORT程序对 12 个真菌基因组的编码蛋白进行
细胞亚定位预测ꎬ获得定位于胞外的蛋白基因集ꎬ
extr≥18[17]ꎻ通过 Signal P4.1 预测蛋白集的信号
肽ꎬD-score≥0.5[18]ꎻ利用跨膜结构预测软件 TM ̄
HMM2.0分析编码蛋白的跨膜结构域[19]ꎬ应用跨
膜拓扑学和信号肽预测软件 Phobius 预测蛋白的
跨膜结构域和信号肽[20]ꎮ 分泌蛋白要求满足下列
条件:亚细胞定位于胞外ꎬ具有信号肽ꎬ无跨膜结构
461
 
  2期 陈相永ꎬ等:植物病原真菌寄生性与分泌蛋白组 CAZymes的比较分析
域ꎮ 获得的分泌蛋白组用于 CAZymes的注释ꎮ
1.3  CAZymes注释与功能聚类 
    从 CAZy 网站( http: / / www. cazy. org / )下载
CAZymes数据库[21]ꎬ通过 BLASTP 方法分别比较
12个基因组分泌蛋白组与数据库 CAZymes 序列
的同源性[22]ꎬE 值≤10-5ꎬ序列一致性≥30%ꎬ取
最佳匹配结果ꎬ得到 12 个真菌分泌蛋白组的
CAZymesꎮ 根据 CAZyme家族的功能进行统计分
析ꎬ并对参与植物细胞壁组分降解的 CAZymes 进
行分类统计[21ꎬ 23]ꎮ
1.4  CAZymes 亚家族进化分析 
    利用 MEGA5.2 软件中的 MUSCLE 程序将建
树需要的同源序列进行联配( alignment)ꎬ用邻近
法(Neighbor-Joining)构建进化树ꎬ选择泊松修正
方法(Poisson correction method)计算进化距离ꎬ
Bootstrap= 1 000[24]ꎮ
2  结果与分析
2.1  不同寄生性植物病原菌分泌蛋白预测和比较
分析
    本研究对 12个不同寄生性的植物病原真菌的
编码基因统计发现ꎬ总编码基因数量在 6 522 ~ 16
694之间ꎬ玉米黑粉病菌最少ꎬ为 6 522 个ꎮ 类似于
玉米黑粉病菌ꎬ另一活体营养型玉米丝黑穗病菌编
码的基因也相对较少ꎬ为 6 673个ꎬ但同样为活体营
养型真菌杨树叶锈病菌则编码 16 694个基因ꎻ半活
体营养型和死体营养型真菌基因组编码的基因均在
10 000个以上ꎬ说明植物病原菌寄生特性与其全基
因组编码基因数目之间没有必然的联系(图 1A)ꎮ
    综合亚细胞定位预测程序 WoLF PSORT、信
号肽预测程序 Signal P4.1、跨膜结构域及其拓扑结
构预测程序 TMHMM2.0 和 Phobius 分析结果ꎬ获
得了 12 个植物病原真菌的分泌蛋白集ꎮ 结果表
明ꎬ不同寄生性病原真菌编码的分泌蛋白数量在
303~1 059 之间ꎬ且编码分泌蛋白的数量和基因组
总编码基因无显著相关性ꎬ如稻瘟病菌共编码
11 054个基因ꎬ含有潜在分泌蛋白 1 059个ꎬ死体营
养型病原菌核盘菌编码 14 522 个基因ꎬ仅含有潜
在分泌蛋白 530 个(图 1A)ꎮ 进一步研究发现ꎬ半
活体营养型和死体营养型真菌中ꎬ除了核盘菌外ꎬ
预测分泌蛋白占总编码蛋白的比例为 4. 8% ~
9.6%ꎬ均高于活体营养型真菌 4. 5% ~ 4. 6% (图
1B)ꎬ一定程度上表现为非专性寄生菌编码分泌蛋
白的比例高于专性寄生物的特征ꎮ
2.2  分泌蛋白组的 CAZymes 
    结果表明ꎬ兼性寄生真菌编码的分泌的
CAZymes的数量显著高于专性寄生真菌ꎬ半活体
营养型真菌编码 CAZyme 的基因为 128 ~ 438 个ꎬ
死体营养型真菌编码的 CAZymes 的数量为 176 ~
289个ꎬ而活体营养型真菌编码的 CAZymes 的数
量仅为64 ~ 119个(图2A) ꎮ另一方面ꎬ不同寄生
Fig. 1  The genomic gene sets and secretome of phytopathogen
561
 
植物病理学报 44卷
性病原真菌编码 CAZymes基因占分泌蛋白基因的
比例也存在较大差异ꎬ除小麦白斑病菌外ꎬ半活体
营养型和死体营养型真菌 CAZymes占分泌蛋白基
因的比例均达到 26.0%以上ꎬ对禾本科植物危害严
重的禾生炭疽菌的 CAZymes 所占的比例高达
41.01%ꎻ而活体营养型真菌中 CAZymes 占分泌蛋
白的比例均偏低ꎬ尤其是杨树叶锈病菌编码的
CAZymes 仅占分泌蛋白集的 15.0%(图 2B)ꎮ 结
果说明ꎬ真菌编码的分泌型 CAZymes 和其寄生性
存在关联ꎬ相对于活体营养型真菌ꎬ非专性寄生真
菌的 CAZymes无论是数量还是占分泌蛋白组的比
例均得到显著增强ꎮ
    对分泌型 CAZymes 的 5 个家族(CBM、CE、
GH、GT和 PL)的比较分析发现ꎬ非专性寄生真菌
的糖基水解酶家族 GH和多糖裂解酶家族 PL相对
于专性寄生真菌发生了扩增ꎮ 活体营养型真菌中ꎬ
杨树叶锈病菌编码的 GH 家族最多为 61 个ꎻ但在
半活体营养型和死体营养型真菌中ꎬ编码 GH家族
基因最少的真菌小麦白斑病菌也达到了 81 个ꎬ禾
生炭疽菌更是高达 212个ꎮ 类似地ꎬ活体营养型真
菌几乎不编码 PL家族基因ꎬ但在非专性寄生类型
真菌中则普遍存在ꎬ多数病原编码 5 个以上的 PL
家族基因ꎬ尤其是大丽轮枝菌ꎬ编码的分泌型 PL
家族基因高达 30个(图 2A)ꎮ GH和 PL家族中多
个亚家族参与了植物细胞组分(纤维素、果胶、木
聚糖等)的降解作用ꎬ而病原真菌 GH 和 PL 家族
的扩增现象呈现出和寄生性相关的特性ꎬ暗示着植
物细胞壁降解酶类可能参与影响了病原菌对寄主
的寄生性ꎮ
2.3   纤维素降解酶与寄生性 
    纤维素降解酶类主要降解植物细胞壁组分纤
维素ꎬ通常具有 β-D-1ꎬ4-葡萄糖酶和纤维二糖酶
功能ꎬ主要包括糖基水解酶的 8 个亚家族(GH1、
GH3、GH5、GH6、GH7、GH9、GH45 和 GH61) [21ꎬ 23]ꎮ
本研究发现ꎬ不同植物寄生性真菌编码的纤维素降
解酶存在较大差异ꎬ活体营养型真菌编码的纤维素
降解酶类相对较少(8~ 16 个)ꎬ而非专性寄生类型
真菌中则显著增多ꎬ共有 6个真菌编码的纤维素降
解酶类超过了 30 个(32 ~ 62 个) (表 1)ꎮ 尤其是
GH3 ( β ̄glucosidase)和 GH61 ( copper - dependent
lytic polysaccharide monooxygenases)亚家族基因
在非活体营养型真菌中发生了显著的扩增ꎮ 3 个
活体营养型真菌中各编码 1个 GH3亚家族基因ꎬ4
个半活体营养型真菌编码 9~ 12 个ꎬ死体营养型真
菌编码 3~8个ꎮ 活体营养型真菌没有预测出编码
GH61亚家族的基因ꎬ4 个半活体营养型真菌编码
GH61亚家族基因差异较大ꎬ小麦白斑病菌编码 1
个ꎬ其它 3种病原真菌分别编码 12 ~ 28 个ꎬ5 种死
体营养型真菌编码的 GH61亚家族基因扩增明显ꎬ
分别为 7~27个(表 1)ꎮ 结果说明ꎬ参与纤维素降
解的 GH3家族和 GH61家族的扩增可能增强了病
Fig. 2  Annotation the CAZymes of secretome in phytopathogens
A: Classification of CAzymes predicted in fungal genomesꎻB: The ratios of secreted CAZymes to fungal secretomes respectively.
661
 
  2期 陈相永ꎬ等:植物病原真菌寄生性与分泌蛋白组 CAZymes的比较分析
Table 1  Comparison of the subfamilies of cellulose degrading enzymes
Pathogenic lifestyle Species GH1 GH3 GH5 GH6 GH7 GH9 GH45 GH61 Total
Biotroph Sporisorium reilianum 0 1 5 0 0 0 2 0 8
Biotroph Ustilago maydis 0 1 6 0 0 0 3 0 10
Biotrophy Melampsora laricis ̄populina 0 1 7 0 7 1 0 0 16
Hemibiotroph Magnaporthe grisea 1 12 7 3 3 0 1 15 42
Hemibiotroph Mycosphaerella graminicola 1 9 4 0 1 0 1 1 17
Hemibiotroph Colletotrichum graminicola 0 10 13 3 6 0 2 28 62
Hemibiotroph Fusarium graminearum 1 9 7 1 1 0 1 12 32
Necrotroph Botrytis cinerea 0 6 5 1 3 0 1 10 26
Necrotroph Cochliobolus heterotrophus 2 8 8 3 5 0 3 18 47
Necrotroph Stagonospora nodorum 0 3 13 4 2 0 2 27 51
Necrotroph Sclerotinia sclerotiorum 2 4 8 1 3 0 2 7 27
Necrotroph Verticillium dahliae 1 6 8 4 6 0 2 19 46
 
原菌由活体营养转向为半活体营养或者死体营养ꎮ
2.4   果胶降解酶与寄生性 
    GH和 PL家族中多个亚家族参与了植物细胞
壁果胶组分的降解作用ꎬ包括 GH2、GH28、GH35、
GH43、GH51、GH53、GH54、GH78、GH93、GH105、
PL1、PL3、PL4、PL9、PL11等ꎬ以及碳水化合物酯酶
CE1、CE8和 CE12[21ꎬ 23]ꎮ 对不同寄生性植物病原
菌编码的果胶降解酶比较分析同样发现ꎬ非专性寄
生类型真菌编码的果胶降解酶类相对于活体营养
型病原真菌发生了显著扩增ꎬ小麦白斑病菌相对较
少ꎬ共编码 18个果胶降解酶类ꎬ其它非专性寄生物
编码 36 ~ 84 个ꎬ而 3 个活体营养型真菌则只编码
10 ~ 14 个(表 2)ꎮ 对个别亚家族比较分析发现ꎬ
GH28( polygalacturonase)、GH43 ( β -xylosidase)、
PL1( pectatelyase)和 CE1 ( acetyl xylan esterase)
4个亚家族的扩增最为明显ꎬ部分非专性寄生菌编
码的基因数量远高于活体营养型真菌ꎬ如玉米黑粉
病菌仅编码 1个 GH28亚家族基因ꎬ而死体营养型
的灰霉菌编码的 GH28 亚家族基因高达 18 个ꎻ此
外ꎬ部分果胶降解酶亚家族(GH53、GH78、GH93、
GH105、PL3、PL4 和 CE12)在活体营养型真菌中
未发现同源基因ꎬ但在多数非专性寄生菌中均有发
现(表 2)ꎮ 结果同样说明病原真菌的寄生性和果
胶降解酶相关ꎬ其通过增强对植物细胞壁组分果胶
降解的能力增强了对寄主的兼性寄生能力ꎮ
    通过在非专性寄生真菌中扩增明显的 PL1 家
族的系统发育树可以看出ꎬ除了与活体营养病原真
菌编码的 2 个 PL1 基因 Spre _ CBQ69665. 1 和
Mellp_72436近缘 Clade1 分枝ꎬ非活体营养类型的
植物致病真菌的 PL1基因出现扩增外(图 3)ꎬ还在
其他进化分枝出现大量的扩增ꎮ 如黄萎病菌在
Clade1具有 4 个 PL1 基因ꎬ在其他进化分支出现
具有多糖裂解功能的 12 个 PL1 基因ꎮ 这意味着
黄萎病菌降解植物细胞壁果胶成分能力的增强ꎬ以
及适应(降解)更多类型的植物细胞壁结构ꎬ进一
步适应更多的植物寄主范围ꎮ 可以推测果胶降解
相关 CAZymes 亚家族 GH28、GH43、CE1 等在进
化上也出现类似情况ꎬ以增强病原真菌适应其非专
性寄生的能力ꎮ
2.5  木聚糖降解酶与寄生性
    除了纤维素和果胶之外ꎬ由 β-D-1ꎬ4 木糖苷
键链接的木聚糖也是植物细胞壁的重要组分ꎮ 参
与木聚糖降解酶的亚家族包括 GH3ꎬGH10ꎬGH11ꎬ
GH43ꎬGH62ꎬGH67ꎬGH115 和 CE1ꎬ其中 GH3 也
同时参与了纤维素的降解ꎬGH43和 CE1 则同时参
与果胶的降解作用[21ꎬ 23]ꎮ 相对于活体营养型真
菌ꎬ非专性寄生真菌中 GH3、GH43 和 CE1 亚家族
的扩增除了增强对纤维素和果胶的降解之外(表
1ꎬ表 2)ꎬ同样也对木聚糖的降解有正向作用(表
3) ꎻ除此之外ꎬ部分非专性寄生真菌的GH10和
761
 
植物病理学报 44卷
Ta
bl

2 
Co

pa
ris
on
of
th

su
bf
am
ili
es
of
pe
ct
in
de
gr
ad
in

en
zy

es
Pa
th
og
en
ic
lif
es
ty
le
Sp
ec
ie






28


35


43


51


53


54


78


93


10

PL

PL

PL

PL

PL
11
CE

CE

CE
12
To
ta

Bi
ot
ro
ph
Sp
or
iso
riu

re
ili
an
um


















11
Bi
ot
ro
ph

sti
la
go

ay
di



















10
Bi
ot
ro
ph

el
am
ps
or

la
ric
is
- p
op
ul
in



















14

em
ib
io
tro
ph

ag
na
po
rth

gr
ise




12














39

em
ib
io
tro
ph

yc
os
ph
ae
re
lla
gr
am
in
ic
ol



















18

em
ib
io
tro
ph
Co
lle
to
tri
ch
um
gr
am
in
ic
ol




11














54

em
ib
io
tro
ph
Fu
sa
riu

gr
am
in
ea
ru




14














59

ec
ro
tro
ph
Bo
try
tis
ci
ne
re


19
















44

ec
ro
tro
ph
Co
ch
lio
bo
lu

he
te
ro
tro
ph
us


















51

ec
ro
tro
ph
St
ag
on
os
po
ra
no
do
ru
















10


42

ec
ro
tro
ph
Sc
le
ro
tin
ia
sc
le
ro
tio
ru


18
















36

ec
ro
tro
ph
Ve
rti
ci
lli
um
da
hl
ia


11

15






16







84
861
 
  2期 陈相永ꎬ等:植物病原真菌寄生性与分泌蛋白组 CAZymes的比较分析
Fig. 3  Phylogenetic analysis of polysaccharide lyases family 1
GH11 亚家族相对于活体营养型真菌也发生了一
定程度的扩增ꎮ 因此ꎬ非专性寄生寄生真菌可能通
过增强木聚糖降解酶类家族基因来影响其对植物
的寄生性ꎮ
3  讨论
    病原菌从寄主植物体内夺取养分和水分等生
活物质的能力存在差异ꎬ形成了活体营养、半活体
营养和死体营养等不同的寄生特性[1]ꎬ然而对寄
生性遗传基础的研究尚不清楚ꎮ 近年来通过比较
基因组学研究发现ꎬ病原菌的胞外分泌蛋白尤其是
碳水化合物酶类 CAZymes 和寄生性有相关
性[14~16]ꎬ因此ꎬ本研究聚焦于不同寄生性病原真菌
分泌蛋白组的 CAZymes 的比较分析ꎬ有助于解析
胞外碳水化合物酶类与植物病原菌寄生性的关系ꎮ
    分泌蛋白组预测发现ꎬ尽管 12 个病原菌基因
组编码基因的数量存在较大差异ꎬ但多数半活体营
养型和死体营养型真菌编码的分泌蛋白占总基因
的比例高于活体营养型真菌(图 1)ꎮ 非专性寄生
物需要更多的胞外酶类或者毒素物质与寄主互作ꎬ
进而杀死寄主的细胞和组织ꎬ而后以死亡的植物组
织作为营养来源[4ꎬ 25]ꎬ而专性寄生物侵入后需要
保持寄主植物的活性[26]ꎬ因此ꎬ非活体病原菌可能
需要编码更多的分泌蛋白来完成半活体营养或死
体营养的特性ꎮ 病原菌编码的胞外分泌蛋白中含
有大量的碳水化合物活性酶ꎬ参与对寄主植物细胞
961
 
植物病理学报 44卷
Table 3  Comparison of the subfamilies of xylan degrading enzymes
Pathogenic lifestyle Species GH3 GH10 GH11 GH43 GH62 GH67 GH115 CE1 Total
Biotroph Sporisorium reilianum 1 1 0 2 1 0 1 1 7
Biotroph Ustilago maydis 1 1 1 4 1 0 1 1 10
Biotroph Melampsora laricis-populina 1 1 0 4 0 0 0 0 6
Hemibiotroph Magnaporthe grisea 12 5 4 12 4 1 1 8 47
Hemibiotroph Mycosphaerella graminicola 9 1 1 5 1 0 1 2 20
Hemibiotroph Colletotrichum graminicola 10 10 6 11 2 1 1 6 47
Hemibiotroph Fusarium graminearum 9 4 3 14 1 1 1 3 36
Necrotroph Botrytis cinerea 6 2 2 4 1 0 1 1 17
Necrotroph Cochliobolus heterotrophus 8 5 5 6 2 1 0 7 34
Necrotroph Stagonospora nodorum 3 6 7 6 3 1 1 10 37
Necrotroph Sclerotinia sclerotiorum 4 2 3 1 0 0 1 0 11
Necrotroph Verticillium dahliae 6 3 5 15 0 1 0 5 35
 
壁的降解作用ꎬ促进病原菌在寄主植物细胞中的附
着、侵入、定殖以及营养吸收等过程[4ꎬ 6]ꎮ 本研究
发现ꎬ病原真菌编码的 CAZymes与寄生性相关ꎬ非
专性寄生真菌编码的 CAZymes的数量和占分泌蛋
白组的比例显著高于专性寄生真菌ꎬ尤其是参与植
物细胞壁降解的糖基水解酶家族 GH 和多糖裂解
酶家族 PL发生了显著扩增(图 2)ꎮ 因此ꎬ非专性
寄生真菌编码植物细胞壁降解酶能力的强化可能
与其通过降解更复杂的植物细胞壁组份来完成死
体营养的生活史密切相关ꎬ而专性寄生真菌胞外分
泌蛋白的功能则更集中于营养的摄取ꎬ而非强制性
的降解植物细胞壁组织[27]ꎮ
    病原菌侵染寄主植物的过程首先需要突破植
物的表皮和细胞壁障碍ꎬ分泌到胞外的 CAZymes
扮演着重要的角色ꎬ尤其是参与植物细胞壁组份降
解的亚家族[7ꎬ 14ꎬ 23ꎬ 26]ꎬ特别是 GH、PL 和 CE 家族
参与了植物多糖组分的降解作用ꎬ是病原菌成功侵
染寄主的重要因素[28ꎬ 29]ꎮ 本研究发现ꎬ参与细胞
壁组份纤维素、果胶和木聚糖降解的亚家族在不同
寄生真菌中发生了显著分化ꎬ相对于活体营养型真
菌仅编码少量的细胞壁组份降解酶类ꎬ半活体营养
和死体营养编码的细胞壁降解酶得到强化ꎬ如纤维
素降解酶类 GH3、 GH61 等 (表 1 )ꎬ果胶酶类
GH28、PL1 等 (表 2 )ꎮ 不同营养类型真菌间
CAZymes亚家族的进化分析也表明ꎬ非活体寄生
真菌不仅具有活体寄生近缘的编码基因出现扩增ꎬ
同时衍生出多拷贝的 CAZymes 同源基因(图 3)ꎬ
以适应非专性寄生菌扩大寄主范围的特点(图 3)ꎮ
一般认为ꎬ死体营养型真菌较活体营养型真菌编码
更多的胞外碳水化合物酶类ꎬ通过植物细胞系的降
解来实现广谱侵染的特性[30]ꎬ而活体营养型真菌
编码植物细胞壁降解酶的能力退化[31]ꎬ如玉米黑
粉病菌通过弱化编码 CAZymes的能力来适应在寄
主植物体内进行活体营养的生活方式[32]ꎮ 因此ꎬ
推测病原菌的寄生性和其编码的植物细胞壁降解
酶类能力相关ꎬ相对于专性寄生菌ꎬ非专性寄生菌
需要编码更多的植物细胞壁降解酶来适应其摄取
营养的生活方式ꎮ
    本研究对不同寄生性的 12种植物病原真菌分
泌蛋白组的 CAZymes 进行了比较分析ꎬ发现兼性
寄生物编码的 CAZymes相对于专性寄生菌发生了
扩增ꎬ特别是参与植物细胞壁组分(纤维素、果胶
和木聚糖)降解的亚家族在兼性寄生菌中得到明
显强化ꎬ而专性寄生菌则弱化编码植物细胞壁降解
酶类ꎬ研究结果为解释植物病原菌寄生性的遗传基
础提供了理论依据ꎮ
参考文献:
[1]   Xu Z G. General plant pathology ( in Chinese) [M] .
Beijing: China Agriculture Press (北京:中国农业出
版社)ꎬ 2006.
071
 
  2期 陈相永ꎬ等:植物病原真菌寄生性与分泌蛋白组 CAZymes的比较分析
[2]   Kamoun S. A catalogue of the effector secretome of
plant pathogenic oomycetes [ J] . Annual Review of
Phytopathologyꎬ 2006ꎬ 44: 41-60.
[3]   Goker Mꎬ Voglmayr Hꎬ Riethmuller Aꎬ et al. How do
obligate parasites evolve? A multi ̄gene phylogenetic
analysis of downy mildews [ J] . Fungal Genetic and
Biologyꎬ 2007ꎬ 44(2): 105-122.
[4]   Knogge W. Fungal infection of plants [ J] . The Plant
Cellꎬ 1996ꎬ 8(10): 1711-1722.
[5]   Rep M. Small proteins of plant ̄pathogenic fungi secre ̄
ted during host colonization [ J] . FEMS Microbiology
Lettersꎬ 2005ꎬ 253(1):19-27.
[6]   Schafer W. Molecular mechanisms of fungal pathoge ̄
nicity to plants [ J] . Annual Review of Phytopatholo ̄
gyꎬ 1994ꎬ 32: 461-477.
[7]   Schmidt S Mꎬ Panstruga R. Pathogenomics of fungal
plant parasites: what have we learnt about pathogenesis
[J] . Current Opinion of Plant Biologyꎬ 2011ꎬ 14(4):
392-399.
[8]   Henrissat Bꎬ Davies G J. Glycoside hydrolases and
glycosyltransferases. Familiesꎬ modulesꎬ and implica ̄
tions for genomics [J] . Plant Physiologyꎬ 2000ꎬ 124:
1515-1519.
[9]   Campbell J Aꎬ Davies G Jꎬ Bulone Vꎬ et al. A
classification of nucleotide ̄diphospho ̄sugar glycosyl ̄
transferases based on amino acid sequence similarities
[J] . Biochemistry Journalꎬ 1997ꎬ 326: 929-942.
[10] Lombard Vꎬ Bernard Tꎬ Rancurel Cꎬ et al. A hierar ̄
chical classification of polysaccharide lyases for
glycogenomics. Biochemistry Journalꎬ 2010ꎬ 432: 437
-444.
[11] Boraston A Bꎬ Bolam D Nꎬ Gilbert H Jꎬ et al. Carbo ̄
hydrate ̄binding modules: fine ̄tuning polysaccharide
recognition [ J] . Biochemistry Journalꎬ 2004ꎬ 382:
769-781.
[12] van den Brink Jꎬ de Vries R P. Fungal enzyme sets for
plant polysaccharide degradation. Applied Microbiology
and Biotechnologyꎬ 2011ꎬ 91(6): 1477-1492.
[13] Klosterman1 S Jꎬ Subbarao K Vꎬ Kang Sꎬ et al. Com ̄
parative genomics yields insights into niche adaptation
of plant vascular wilt pathogens [ J] . PLoS Pathogenꎬ
2011ꎬ 7(7): e1002137.
[14] Lévesque C Aꎬ Brouwer Hꎬ Cano Lꎬ et al. Genome
sequence of the necrotrophic plant pathogen Pythium
ultimum reveals original pathogenicity mechanisms and
effector repertoire [ J] . Genome Biologyꎬ 2010ꎬ 11:
R73.
[15] Hane J Kꎬ Lowe R Gꎬ Solomon P Sꎬ et al. Dothideo ̄
mycete - plant interactions illuminated by genome
sequencing and EST analysis of the wheat pathogen
Stagonospora nodorum [ J] . The Plant Cellꎬ 2007ꎬ
19: 3347-3368.
[16] Ellwood S Rꎬ Liu Z Hꎬ Syme R Aꎬ et al. A first
genome assembly of the barley fungal pathogen
Pyrenophora teres f. teres [ J ] . Genome Biologyꎬ
2010ꎬ 11: R109.
[17] Horton Pꎬ Park K Jꎬ Obayashi Tꎬ et al. WoLF
PSORT: protein localization predictor [ J] . Nucleic
Acids Researchꎬ 2007ꎬ 35: W585-587.
[18] Emanuelsson Oꎬ Brunak Sꎬ von Heijne Gꎬ et al.
Locating proteins in the cell using TargetPꎬ SignalP
and related tools [J] . Nature Protocolꎬ 2007ꎬ 2: 953-
971.
[19] E Lꎬ von Heijne Gꎬ Krogh A. A hidden Markov model
for predicting transmembrane helices in protein
sequences [ J ] . Proceedings of the International
Conference on Intelligent Systems for Molecular Biolo ̄
gyꎬ 1998ꎬ 6: 175-182.
[20] Käll Lꎬ Krogh Aꎬ Sonnhammer E L L. Advantages of
combined transmembrane topology and signal peptide
prediction—the Phobius web server [J] . Nucleic Acids
Researchꎬ 2007ꎬ 35: W429-432.
[21] Cantarel B Lꎬ Coutinho P Mꎬ Rancurel Cꎬ et al. The
carbohydrate ̄active enzymes database (CAZy): an ex ̄
pert resource for glycogenomics [ J] . Nucleic Acids
Researchꎬ 2009ꎬ 37(Database issue): 233-238.
[22] Altschul S Fꎬ Madden T Lꎬ Schäffer A Aꎬ et al.
Gapped BLAST and PSI ̄BLAST: a new generation of
protein database search programs [ J] . Nucleic Acids
Researchꎬ 1997ꎬ 25(17): 3389-3402.
[23] Battaglia Eꎬ Benoit Iꎬ van den Brink Jꎬ et al. Carbohy ̄
drate ̄active enzymes from the zygomycete fungus
Rhizopus oryzae: a highly specialized approach to car ̄
bohydrate degradation depicted at genome level [ J] .
171
 
植物病理学报 44卷
BMC Genomicsꎬ 2011ꎬ 12: 38.
[24] Tamura Kꎬ Peterson Dꎬ Peterson Nꎬ et al. MEGA5:
Molecular evolutionary genetics analysis using
maximum likelihoodꎬ evolutionary distanceꎬ and maxi ̄
mum parsimony methods [ J] . Molecular Biology and
Evolutionꎬ 2011ꎬ 28: 2731-2739.
[25] Pegg G Fꎬ Brady B L. Verticillium wilts. New York:
CABI Publishingꎬ 2002.
[26] Schirawski Jꎬ Mannhaupt Gꎬ Munch Kꎬ et al. Patho ̄
genicity determinants in smut fungi revealed by genome
comparison [J] . Scienceꎬ 2010ꎬ 330(6010): 1546-
1548.
[27] Mueller Oꎬ Kahmann Rꎬ Aguilar Gꎬ et al. The secre ̄
tome of the maize pathogen Ustilago maydis [J] . Fun ̄
gal Genetic and Biologyꎬ 2008ꎬ 45 Suppl. 1: S63-70.
[28] Ospina-Giraldo M Dꎬ Griffith J Gꎬ Laird E Wꎬ et al.
The CAZyome of Phytophthora spp.: a comprehensive
analysis of the gene complement coding for carbohy ̄
drate ̄active enzymes in species of the genus Phytoph ̄
thora [J] . BMC Genomicsꎬ 2010ꎬ 11: 525.
[29] Amselem Jꎬ Cuomo C Aꎬ van Kan J Aꎬ et al. Genomic
analysis of the necrotrophic fungal pathogens Sclerotinia
sclerotiorum and Botrytis cinerea [ J] . PLoS Geneticsꎬ
2011ꎬ 7(8): e1002230.
[30] Oliver R Pꎬ Ipcho S V S. Arabidopsis pathology
breathes new life into the necrotrophs ̄vs. ̄biotrophs
classification of fungal pathogens [J] . Molecular Plant
Pathologyꎬ 2004ꎬ 5: 347-352.
[31] Kemen Eꎬ Jones J D. Obligate biotroph parasitism: can
we link genomes to lifestyles [J] . Trends in Plant Sci ̄
enceꎬ 2012ꎬ 17 (8): 448-57.
[32] Kämper Jꎬ Kahmann Rꎬ Bölker Mꎬ et al. Insights
from the genome of the biotrophic fungal plant patho ̄
gen Ustilago maydis [ J] . Natureꎬ 2006ꎬ 444(7115):
97-101.
责任编辑:李晖
271