欢迎访问《热带作物学报》，
热带作物学报

热带作物学报 >

2019 , Vol. 40 >Issue 12: 2456 - 2465

DOI: https://doi.org/10.3969/j.issn.1000-2561.2019.12.020

植物病理、作物病虫害及其防治

我国热带作物多主棒孢种群多样性及致病力分化分析

  • 李博勋 ,
  • 冯艳丽 ,
  • 刘先宝 ,
  • 蔡吉苗 ,
  • 陆翠梅 ,
  • 郑肖兰 ,
  • 黄贵修 , *
展开
  • 中国热带农业科学院环境与植物保护研究所/农业农村部热带作物有害生物综合治理重点实验室/海南省热带农业有害生物监测与控制重点实验室,海南海口 571101
黄贵修,E-mail：

李博勋（1986—）,女,硕士,助理研究员,研究方向：植物病理学。

收稿日期: 2019-03-15

  要求修回日期: 2019-05-13

  网络出版日期: 2019-12-20

基金资助

国家重点研发计划项目(No.2017YFC1200600)

海南省自然科学基金面上项目(No.317233)

农业部现代农业人才支撑计划项目(No.0316001)

版权

版权所有，未经授权，不得转载、摘编本刊文章，不得使用本刊的版式设计。
收起

Genetic Diversity and Pathogenic Variability Among Isolates of Corynespora cassiicola from Tropical Crops in China

  • LI Boxun ,
  • FENG Yanli ,
  • LIU Xianbao ,
  • CAI Jimiao ,
  • LU Cuimei ,
  • ZHENG Xiaolan ,
  • HUANG Guixiu , *
Expand
  • Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences / Key Laboratory of Integrated Pest Management on Tropical Crops, Ministry of Agriculture and Rural Affairs / Hainan Key Laboratory for Monitoring and Control of Tropical Agricultural Pests, Haikou, Hainan 571101, China
HUANG Guixiu,E-mail：

Received date: 2019-03-15

  Request revised date: 2019-05-13

  Online published: 2019-12-20

Copyright

Copyright reserved © 2019. Office of Acta Agronomica Sinica All articles published represent the opinions of the authors, and do not reflect the official policy of the Chinese Medical Association or the Editorial Board, unless this is clearly specified.
Fold

摘要

多主棒孢（Corynespora cassiicola）是一种全球分布的重要植物病原真菌,寄主范围十分广泛。近年来,该病原菌侵染危害多种热带作物,并在局部地区大面积暴发流行,造成巨大的经济损失。本研究利用核糖体转录间隔区（Internal transcribed spacer, ITS）、翻译延伸因子（translation elongation factor 1 alpha, EF-1α）、β微管蛋白（β-Tubulin, TUB）,对来自橡胶树、木薯、番木瓜和瓜菜等多种热带作物的70株多主棒孢进行种群多样性分析。结果发现,利用最大似然法,在相似系数为0.97时可将供试菌株划分为两大遗传类群：类群Ⅰ为国内和一些国外橡胶树的多主棒孢;类群Ⅱ为其他寄主多主棒孢,遗传类群与寄主来源具有显著的相关性。致病力测定表明,多主棒孢种内致病力分化明显,类群Ⅰ的橡胶树多主棒孢仅能侵染橡胶树,不能侵染其他寄主,具有明显的寄主专化性,且不同地理来源的菌株致病力差异不显著。类群Ⅱ的多主棒孢在不同寄主间可以相互侵染,但菌株在致病力和发病症状上存在一定差异,表现出对其原寄主的高度致病力。多基因序列聚类分析与致病力分化分析结果具有较好的一致性,这为研究我国热带作物多主棒孢的种群结构、主要致病型以及为多主棒孢病害的发生、流行和防治提供了理论依据。

关键词： 热带作物; 多主棒孢病菌; 种群多样性; 致病力分化

本文引用格式

李博勋 , 冯艳丽 , 刘先宝 , 蔡吉苗 , 陆翠梅 , 郑肖兰 , 黄贵修 . 我国热带作物多主棒孢种群多样性及致病力分化分析[J]. 热带作物学报, 2019 , 40(12) : 2456 -2465 . DOI: 10.3969/j.issn.1000-2561.2019.12.020

Abstract

Corynespora cassiicola is an important plant pathogenic fungus distributed globally and has a wide host range. In recent years, C. cassiicola had infected a variety of tropical crops. The diseases caused by C. cassiicola has outbroken in some areas, causing huge economic losses. In this study, we used ITS, EF-1α and β-tubulin genes to analyze the population diversity of isolates of C. cassiicola collected from rubber trees, cassava, papaya, tropical fruits and vegetables. When the similarity coefficient was 0.97, the tested isolate could be divided into two divergent clusters by the maximum likelihood method. Cluster I was the isolates of rubber trees, which came from China and some foreign countries, and cluster II was the isolates of other hosts. There was significant correlation between genetic clusters and host source. Pathogenicity test showed that the pathogenic differentiation of C. cassiicola was obvious. The isolates of cluster I could only infect rubber trees, but could not infect other hosts. They had obvious host specificity, and there was no significant difference in pathogenicity among isolates from different geographical sources. C. cassiicola in cluster II could infect each other among different hosts, but the pathogenicity and symptoms of the isolates were different, showing a high pathogenicity to the original hosts. The results of genetic diversity and pathogenic variability showed a good consistency, which would provide a theoretical basis for understanding the population structure and main pathogenic types among the isolates of C. cassiicola from tropical crops in China, as well as the occurrence, epidemic and control method of the diseases caused by C. cassiicola.

Key words： tropical crops; Corynespora cassiicola; population diversity; pathogenic variability

多主棒孢病菌（Corynespora cassiicola）是为害我国热带作物的主要病原真菌,其寄主范围广、症状类型多样、生理小种复杂。由该病原菌侵染引起的橡胶树棒孢霉落叶病是橡胶树最为严重的病害之一,在全球主要产胶国普遍发生,造成的经济损失巨大[1]。同时该菌也能为害木薯、番木瓜、香蕉等热带果蔬,并在局部地区暴发流行。李宝聚等[2]对我国蔬菜多主棒孢病害进行了调查,发现海南、云南、广东、广西、山东、河北、北京、辽宁、吉林、河南等地的黄瓜、茄子、番茄、菜豆、辣椒受多主棒孢的危害十分严重,每年发生面积超过66.7万hm2,损失超过50亿元。再加上该病原菌的侵染破坏力极强,给防治带来了巨大困难,近年来已成为我国主要热带作物的主要病害之一。
多主棒孢在分子水平上存在明显的差异,具有丰富的遗传多样性。Silva利用RFLP和RAPD标记将采自斯里兰卡和澳大利亚地区的多个寄主植物上的多主棒孢划分为7个遗传类群,发现遗传类群和样本采集地、致病力、菌株的生长速率之间具有显著相关性[3]。但也有学者指出,多主棒孢的遗传类群与致病型、地理或寄主来源不具有相关性[4]。Smith等[5]利用高变异核位点GA4,分析了4个国家58个寄主上多主棒孢进行系统发育分析,发现其种内遗传多样性丰富,而且不同类群与致病力和寄主来源密切相关。Dixon等[6]利用Cassiicolin毒素基因、ITS、ga4、caa5以及act1基因序列分析了不同寄主植物上的143株多主棒孢,进一步证实了遗传类群与寄主来源、致病性和生长速率有关,但是与地理来源无关。Shimomoto等[7]利用多基因（β-微管蛋白、EF-1α、钙调蛋白和肌动蛋白）分析了64株多主棒孢,发现遗传类群与致病型一致。随着分子生物学技术的发展和对多主棒孢病害研究的深入,学者们开始利用不同毒素类群的Cassiicolin毒素基因对多主棒孢的致病型进行了详细的划分,这也是现阶段研究多主棒孢分类的主要依据[8]。Déon等[8]分析了橡胶树上6个Cassiicolin毒素类群（Cas1~ Cas6型）基因的11个菌株对两个橡胶树品种（RRIM600、PB260）的致病力,发现不同毒素类型的菌株在致病力强弱上存在显著差异,Cas1型的致病力最强。而对于不同寄主来源多主棒孢的遗传多样性与致病力分化之间的研究进展很少。为此,本研究利用多基因序列和致病力分化测定等方法,分析我国主要热带作物和果蔬多主棒孢病菌的遗传种群结构,进一步明确我国多主棒病菌的优势种群及其致病力分化情况,这为该病原菌的发生和流行规律、生理小种的确定以及防控策略的制定提供理论依据。

1 材料与方法

1.1 材料

1.1.1 供试菌株和材料 70株多主棒孢（表1）由中国热带农业科学院环境与植物保护研究所热带工业原料作物病害研究组分离和保存。致病力测定选用的橡胶树和木薯种苗均由中国热带农业科学院环境与植物保护研究所演丰基地保存,瓜菜种苗为市场购买种子培育。
表1 多主棒孢信息和多基因序列登录号

Tab. 1 Information and GenBank accession numbers of the C. cassiicola isolates

菌株编号
Stains		 采集年
Year		 地理来源
Geographical origin		 寄主
Host		 ITS登录号
GenBank		 EF-1α登录号
GenBank		 TUB登录号
GenBank
HCcGX01 2007 中国广西省东兴市马路镇 Hevea brasiliensis MK571285 MK589895 MK589805
CpCcHNDZ02 2007 中国海南省儋州市两院 Carica papaya MK571291 MK589901 MK589807
CpCcYN02 2007 中国云南省西双版纳景洪农场 Carica papaya MK571295 MK589905 MK589809
CpCcHN07 2007 中国海南省澄迈县红光农场 Carica papaya MK571297 MK589907 MK589810
HCcYNMX1504 2013 中国云南省景洪市勐醒农场 Hevea brasiliensis MK571313 MK589923 MK589811
HCcYNMX1501 2013 中国云南省景洪市勐醒农场 Hevea brasiliensis MK571314 MK589924 MK589812
HCcYNRZS1506 2013 中国云南省景洪市热带作物研究所 Hevea brasiliensis MK571315 MK589925 MK589813
CpCcYNMX1501 2013 中国云南省景洪市勐醒农场 Carica papaya MK571316 MK589926 MK589814
HCcGDMM1402 2014 中国广东省化州市农垦热带作物科学研究所 Hevea brasiliensis MK571317 MK589927 MK589815
VsCcHN10 2010 中国海南省三亚南繁科学技术研究院 Vigna unguiculata MK571324 MK589934 MK589816
VsCcHN09 2010 中国海南省三亚南繁科学技术研究院 Vigna unguiculata MK571325 MK589935 MK589817
CmCcHN03 2010 中国海南省三亚南繁科学技术研究院 Cucumis melo MK571328 MK589938 MK589818
SmCcHN01 2010 中国海南省三亚南繁科学技术研究院 Solanum melongena MK571330 MK589940 MK589819
LJCcHN01 2016 中国海南省三亚南繁科学技术研究院 Capsicu mannuum MK571331 MK589941 MK589820
LJCcHN02 2016 中国海南省三亚南繁科学技术研究院 Capsicu mannuum MK571332 MK589942 MK589821
HCcYNRL04 2013 中国云南省瑞丽市遮放农场 Hevea brasiliensis MK571333 MK589943 MK589822
HCcYNRL05 2013 中国云南省瑞丽市遮放农场 Hevea brasiliensis MK571334 MK589944 MK589823
HCcYNBN06 2013 中国云南省景洪市景洪农场 Hevea brasiliensis MK571335 MK589945 MK589824
HCcHNWZS01 2013 中国海南省五指山市潘阳镇 Hevea brasiliensis MK571338 MK589948 MK589826
VsCcHN07 2010 中国海南省三亚市三亚南繁科学技术研究院 Hevea brasiliensis MK571342 MK589952 MK589828
HCcYNBD01 2013 中国云南省保定县市 Hevea brasiliensis MK571343 MK589953 MK589829
HCcGDYC06 2013 中国广东省阳春市红十月农场 Hevea brasiliensis MK571344 MK589954 MK589830
HCcHNCJ02 2013 中国海南省昌江县十月田镇 Hevea brasiliensis MK571349 MK589959 MK589833
HCcHNQZ01 2013 中国海南省琼中县阳江农场 Hevea brasiliensis MK571351 MK589961 MK589834
HCcHNBS03 2013 中国海南省白沙县邦溪大岭农场 Hevea brasiliensis MK571352 MK589962 MK589835
HCcHNBS01 2013 中国海南省白沙县邦溪大岭农场 Hevea brasiliensis MK571353 MK589963 MK589836
HCcHNTC01 2013 中国海南省屯昌县中建二队 Hevea brasiliensis MK571356 MK589966 MK589837
HCcYNHK02 2013 中国云南省河口市蚂蝗堡农场 Hevea brasiliensis MK571357 MK589967 MK589838
HCcHNWC01 2013 中国海南省文昌市海垦橡胶所 Hevea brasiliensis MK571358 MK589968 MK589839
HCcHNWZS03 2013 中国海南省五指山市潘阳镇 Hevea brasiliensis MK571359 MK589969 MK589840
HCcGDZJ04 2013 中国广东省湛江市热科院南亚所 Hevea brasiliensis MK571360 MK589970 MK589841
HCcGDYJ1603 2013 中国广东省阳江市鸡山农场 Hevea brasiliensis MK571361 MK589971 MK589842
HCcYNJP02 2013 中国云南省金平县 Hevea brasiliensis MK571362 MK589972 MK589843
PaCcHB03 2007 中国河北省廊坊市 Cucumis sativus MK571366 MK589976 MK589846
PaCcSD02 2007 中国山东省寿光市 Cucumis sativus MK571367 MK589977 MK589847
PaCcSD03 2007 中国山东省寿光市 Cucumis sativus MK571368 MK589978 MK589848
PaCcSD04 2007 中国山东省寿光市 Cucumis sativus MK571369 MK589979 MK589849
HCcHNDZ1608 2016 中国海南省儋州国家种质圃 Hevea brasiliensis MK571371 MK589981 MK589850
HCcHNDZ1615 2016 中国海南省儋州国家种质圃 Hevea brasiliensis MK571372 MK589982 MK589851
HCcHNHK1630 2016 中国海南省海口市演丰镇 Hevea brasiliensis MK571374 MK589984 MK589852
HCcHNHK1631 2016 中国海南省海口市演丰镇 Hevea brasiliensis MK571375 MK589985 MK589853
HCcHNWN1633 2016 中国海南省万宁市东兴农场 Hevea brasiliensis MK571376 MK589986 MK589854
HCcHNWN1639 2016 中国海南省万宁市东兴农场 Hevea brasiliensis MK571377 MK589987 MK589855
HCcHNLD1648 2016 中国海南省乐东县志仲镇 Hevea brasiliensis MK571379 MK589989 MK589856
菌株编号
Strains		 采集年
Year		 地理来源
Geographical origin		 寄主
Host		 ITS登录号
GenBank		 EF-1α登录号
GenBank		 TUB登录号
GenBank
HCcHNLD1654 2016 中国海南省乐东县保显农场 Hevea brasiliensis MK571383 MK589993 MK589857
HCcHNBS1679 2016 中国海南省白沙县邦溪大岭农场 Hevea brasiliensis MK571389 MK589999 MK589858
HCcHNQZ1689 2016 中国海南省琼中县阳江农场 Hevea brasiliensis MK571390 MK590000 MK589859
CpCcHNQZ1603 2016 中国海南省琼中县阳江农场 Carica papaya MK571391 MK590001 MK589860
MeCcHN01 2007 中国海南省乐东县保显农场 Hevea brasiliensis MK571392 MK590002 MK589861
LsCcHNSY1101 2011 中国海南省三亚南繁科学技术研究院 Lactuca sativa MK571393 MK590003 MK589862
LsCcHNSY1102 2011 中国海南省三亚南繁科学技术研究院 Lactuca sativa MK571394 MK590004 MK589863
LeCcHNSY1702 2011 中国海南省三亚南繁科学技术研究院 Hevea brasiliensis MK571395 MK590005 MK589864
PaCcHNSY1708 2011 中国海南省三亚市天涯区 Cucumis sativus MK571396 MK590006 MK589865
PaCcHNSY1705 2011 中国海南省三亚市涯城区 Cucumis MK571397 MK590007 MK589866
PaCcHNSY1703 2011 中国海南省三亚市天涯区 Cucumis MK571398 MK590008 MK589867
MeCcGD02 2007 中国广东省湛江市 Manihot esculenta MK571399 MK590009 MK589868
HCcYN49 2007 中国云南省河口市蚂蝗堡农场 Hevea brasiliensis MK571401 MK590011 MK589870
CpCcHNQZ1605 2016 中国海南省琼中县阳江农场 Carica papaya MK571402 MK590012 MK589871
PaCcHNSY1711 2017 中国海南省三亚市涯城区 Cucumis sativus MK571404 MK590014 MK589872
CmCcHN01 2010 中国海南省三亚市南繁院 Cucumis melo MK571326 MK589936 MK589874
HCcYNHZM01-2 2015 越南 Hevea brasiliensis MK571405 MK590015 MK589873
HCcYNHZM01-1 2015 越南 Hevea brasiliensis MK571363 MK589973 MK589844
HCcYNHZM01-3 2015 越南 Hevea brasiliensis MK571364 MK589974 MK589845
HCcJPZ03 2013 柬埔寨 Hevea brasiliensis MK571347 MK589957 MK589831
HCcJPZ05 2013 柬埔寨 Hevea brasiliensis MK571348 MK589958 MK589832
HCcJPZ04 2013 柬埔寨 Hevea brasiliensis MK571341 MK589951 MK589827
CpCcJPZ01 2013 柬埔寨 Carica papaya MK571337 MK589947 MK589825
HCcYD1 2007 印度 Hevea brasiliensis MK571293 MK589903 MK589808
CYDB2 2007 印度 Hevea brasiliensis MK571286 MK589896 MK589806
CYDB3 2007 印度 Hevea brasiliensis MK571370 MK589980 MK589869
1.1.2 培养基和试剂 PDA培养基参照《植病研究方法》配制[9],用于多主棒孢菌株的培养。引物（见表1）均由华大基因股份有限公司合成;克隆载体pMD18-T、PCR扩增反应试剂购于TAKARA公司;pMD-18T载体、Taq DNA聚合酶、dNTP等购自宝生物工程（大连）有限公司。DNA Marker、Trans5α大肠杆菌E. coli感受态购自北京全式金生物技术有限公司。其他化学试剂均为国产分析纯。
1.1.3 仪器与设备 Bio-Rad T100型梯度PCR仪,美国伯乐公司; UVI FireReader凝胶成像系统,英国UVItec公司。

1.2 方法

1.2.1 DNA提取、PCR扩增和测序 采用CTAB方法[10]提取病原菌DNA,用超微量分光光度计进行DNA质量测定,用ITS、EF-1α、TUB、Cassiicolin 等多基因的引物序列（表2）对病原菌目的基因进行PCR扩增,反应条件参照Déon等[8]的方法,扩增产物用1%的琼脂糖凝胶电泳进行检测,经Omega Bio-Tek试剂盒进行目的条带回收,之后克隆到pMD18-T载体,转化Trans5α大肠杆菌感受态细胞,筛选阳性克隆转化子,送其中的3个克隆至华大基因技术有限公司进行序列测定。
表2 多基因序列引物

Tab. 2 Polygenic sequence primers

基因
Genes		 引物
Primers		 序列(5'-3')
Sequence (5'-3')		 退火温度
Annealing temperature/℃		 文献来源
Reference
ITS ITS1
ITS4		 TCCGTAGGTGAACCTGCGG
TCCTCCGCTTATTGATATGC		 53 Shimomoto等[7]
EF-1α EF1-728F
EF1-986R		 CATCGAGAAGTTCGAGAAGG
TACTTGAAGGAACCCTTACC		 53 Shimomoto等[7]
TUB Tub-F
Tub-R		 CCTCCAAACCGGTCAATG
CTGGGTCAACTCGGGGAC		 56 Shimomoto等[7]
Cas5 Cas5-F
Cas5-R		 CGGGGAGGTATCAGGTGTGAGATA
CAGAACAAGCCAAAAGAGAACTAC		 57 Déon等[8]
Cas2 Cas2-F
Cas2-R		 GGATTTGCCTGAGATCCTA
CAAACAATGCTAACCAAACAAAC		 53 Déon等[8]
1.2.2 多基因序列分析、遗传距离计算和系统发育树构建 将序列在NCBI中进行BLAST分析,利用CodonCode Aligner V7.0.1去除测序峰图的低质量区域、完成峰图校对拼接和引物序列切除,应用PAUP 4.0软件计算种内遗传距离;利用MEGA 7.0软件最大似然法中的Tamura-Nei model构建系统发育树。
1.2.3 致病力分化分析 根据多基因序列分析结果,从每个类群里挑选代表性的27个菌株,用无菌水配制孢子悬浮液浓度为105个/mL,于盆栽黄瓜长出3片真叶,番茄、辣椒、豇豆长出4片真叶时进行喷雾接种,每株接种15 mL孢子悬浮液,每个菌株设3次重复,以清水接种为对照,每个重复15株苗,接种后置于25 ℃光照培养箱中保湿培养。在接种后7~14 d,按瓜菜棒孢叶斑病情分级标准[11]进行统计。橡胶树叶片接种法参照农业行业标准《热带作物种质资源抗病虫鉴定技术规程橡胶树棒孢霉落叶病》（NY/T 3195-2018）[12]。致病性分析所获数据均采用SAS 9.0软件进行统计分析,多重比较采用Duncan’s新复极差检验（LSR test）。

2 结果与分析

2.1 ITS、EF-1α、TUB单基因序列分析

用ITS1/ITS4、EF1-728F/EF1-986R和TUB-F/ R的3种引物对供试的70株多主棒孢菌株进行ITS区、EF-1α和TUB序列的扩增,分别得到大小约为531、304、954 bp的条带。所有菌株序列在NCBI经BLAST比对分析,均为多主棒孢（C. cassiicola）,并将各单基因的序列提交至GenBank（表1）。通过序列间的比对发现,在ITS1区的第87位点存在单核苷酸多态性（C/T）,其他序列完全一致。
系统发育树的分析结果显示,ITS、EF-1α、TUB 3个单基因序列均能将供试的71个多主棒孢菌株划分为2大遗传类群（图1,图2,图3）,第Ⅰ类群均为橡胶树多主棒孢,来源于国内各橡胶主产区以及境外植胶国,不同地理来源的菌株之间没有显著的差异;第Ⅱ类群为其他寄主多主棒孢,包括木薯、番木瓜、黄瓜、番茄、豇豆、茄子、西甜瓜等,但也有8株橡胶多主棒孢菌株（占供试橡胶多主棒孢的17%）,其中3株（HCcYD1、CYDB3、CYDB2）来自印度,3株（HCcYNHZM01-1、HCcYNHZM01-2、HCcYNHZM01-3）来自越南,2株（HCcGDMM1402、HCcGDMM1406）来自广东茂名,这8个菌株与第Ⅰ类群的橡胶多主棒孢遗传距离较远,存在显著的差异。从聚类图的分支上看,印度橡胶多主棒孢（HCcYD1、CYDB3、CYDB2）与番木瓜多主棒孢同在一个小分支,亲缘关系最近;越南橡胶多主棒孢（HCcYNHZM01-2）与莴笋、木薯多主棒孢亲缘关系最近,而HCcYNHZM01-1、HCcYNHZM01-3与辣椒多主棒孢同在一个小分支;来自广东茂名的2株橡胶多主棒孢（HCcGDMM1402、HCcGDMM1406）与黄瓜多株棒孢同在一个小分支,亲缘关系最近。
图1 采用Maximum Likelihood构建多主棒孢ITS基因系统发育树

Fig. 1 Phylogenetic trees of C. cassiicola isolates based on ITS sequence

图2 采用Maximum Likelihood构建多主棒孢EF-1α基因系统发育树

Fig. 2 Phylogenetic trees of C. cassiicola isolates based on EF-1α sequence

图3 采用Maximum Likelihood构建多主棒孢TUB基因系统发育树

Fig. 3 Phylogenetic trees of C. cassiicola isolates based on TUB sequence

从上述结果可以看出,3个单基因的序列的聚类图基本上能将橡胶多主棒孢和其他寄主来源的多主棒孢进行有效区分,而且不同寄主之间的遗传差异显著高于地理来源之间的差异。说明寄主因素与多主棒孢的遗传变异间的相关性更为显著,但也有少部分橡胶多主棒孢与其他寄主来源多主棒孢的遗传距离较为接近。从遗传进化的角度分析,这少部分橡胶多主棒孢可能来源于其他寄主,存在交互侵染的现象,也有可能为橡胶多主棒孢的一个新的类群。

2.2 ITS- EF-1α-TUB多基因序列分析

采用3种多基因的部分序列进行拼接,共同构建系统发育树,对来自柬埔寨、越南、印度以及我国的6个省份,10个寄主（橡胶树、木薯、番木瓜、黄瓜、茄子、辣椒、豇豆、莴笋、西甜瓜、番茄）上的70株多主棒孢进行遗传进化树的构建（图4）,发现不同来源的多主棒孢菌株间存在丰富的遗传多样性。根据多基因序列聚类分析,利用最大似然法,在相似系数为0.97时可将供试菌株被划分为两大遗传类群。第Ⅰ类群,为橡胶多主棒孢,来源于我国云南、海南、广东和广西等橡胶主产区以及柬埔寨、越南等国外的橡胶多主棒孢,这一类群的菌株与年份、地理来源之间不具有相关性。但越南（HCcYNHZM01-1）和柬埔寨（HCcJPZ03、HCcJPZ04、HCcJPZ05）的橡胶多主棒孢菌株虽然与国内的橡胶多主棒孢同属于一个类群,但处于不同的小分支上,推测国内和国外的橡胶多主棒孢菌株间存在略微的差异,而国内各植胶区的橡胶多主棒孢间无明显差异,所有菌株都聚类在一个分支上,支持率为98%（图4）。第Ⅱ类群,为不同寄主来源的多主棒孢菌株,与单基因序列聚类分析结果一致。其中供试的所有番木瓜多主棒孢CpCcHNQZ1603、CpCcYNMX1501、CpCcHNDZ02、CpCcHNQZ1605、CpCcYN02、CpCcJPZ01（柬埔寨番木瓜）和所有黄瓜多主棒孢PaCcHB03、PaCcSD02、PaCcSD04、PaCcHNSY1703、PaCcHNSY1705、PaCcHNSY1708、PaCcHNSY1711各自成为一个小分支,单独聚为一类,支持率为99%和63%。西甜瓜多主棒孢（CmCcHN01、CmCcHN03）、豇豆多主棒孢（VsCcHN07、VsCcHN09、VsCcHN10）和茄子多主棒孢（SmCcHN01）分属于同一个小分支,支持率为65%。木薯多主棒孢（MeCcHN01、MeCcGD02）和莴笋多主棒孢（LsCcHNSY1101、LsCcHNSY1102）的亲缘关系最近。
从多基因序列分析的结果可以看出,不同寄主来源的多主棒孢菌株之间存在显著的差异,能将国内、外的橡胶多主棒孢菌株与其他寄主来源的多主棒孢菌株严格区分开,不同地理来源的菌株之间没有显著的差异,而且种内相似度较高,推测这类菌株可能存在寄主专化性,需要通过致病力测定进一步验证。除橡胶之外的其他寄主多主棒孢的种群分化关系不明显,尽管有些寄主的多主棒孢聚类在不同的小分支上,但遗传距离还是比较近,推测第 Ⅱ 类群的多主棒孢之间存在交互侵染的可能。

2.3 不同寄主来源多主棒孢的致病力分化分析

根据多基因序列分析结果,从每个类群中选取出27株不同寄主来源的多主棒孢进行致病力测定。在此之前,参照法国学者报道的多主棒孢Cassiicolin毒素类型的划分[8],明确了这27株多主棒孢的毒素类型（表3）。采用SAS软件的Ward最小方差进行聚类分析,将27株多主棒孢分为Clade A和Clade B 2种主要的致病型（图4,图5,表3）,其中Clade A以Cas5毒素类型为主,Clade B型均为Cas2毒素类型,目前供试的多主棒孢菌株仅存在2种毒素类型,与刘先宝等[13]2016年报道的结果一致。
图4 采用Maximum Likelihood构建多主棒孢ITS-EF-1α-TUB的多基因系统发育树及对应的致病型

Fig. 4 Phylogenetic trees of C. cassiicola isolates based on ITS-EF-1α-TUB sequence and Phylogenetic clade

图5 27株不同来源多主棒孢致病力分化聚类分析图

Fig. 5 Pathogenicity differentiation of 27 isolates with Ward cluster

表3 不同寄主来源多主棒孢病菌的致病分化分析

Tab. 3 Pathogenicity of C. cassiicola isolates derived from different hosts

菌株编号
Isolate code		 寄主
Host		 采集地点
Geographic location		 分离年份
Collection year		 致病型
Phylogenetic clade		 毒素类型
Cassiicolin isoforms		 病情指数Disease index
橡胶
Rubber		 辣椒
Pepper		 黄瓜
Cucumber		 豇豆
Cowpea		 番茄
Tomato
HCcYN15 橡胶 中国云南版纳 2007 A1 Cas5 27.7 0 0 0 0
HCcGX01 橡胶 中国广西东兴 2007 A1 Cas5 19.4 0 0 0 0
HCcYNBN06 橡胶 中国云南版纳 2013 A1 Cas5 34.6 0 0 0 0
HCcHNBS03 橡胶 中国海南白沙 2013 A1 Cas5 25.1 0 0 0 0
HCcYNHK02 橡胶 中国云南河口 2013 A1 Cas5 31.0 0 0 0 0
HCcYNRL04 橡胶 中国云南瑞丽 2013 A1 Cas5 18.3 0 0 0 0
HCcGDYC06 橡胶 中国广东阳春 2013 A1 Cas5 17.2 0 0 0 0
菌株编号
Isolate code		 寄主
Host		 采集地点
Geographic location		 分离年份
Collection year		 致病型
Phylogenetic clade		 毒素类型
Cassiicolin isoforms		 病情指数Disease index
橡胶
Rubber		 辣椒
Pepper		 黄瓜
Cucumber		 豇豆
Cowpea		 番茄
Tomato
HCcGDZJ04 橡胶 中国广东湛江 2013 A1 Cas5 20.1 0 0 0 0
HCcJPZ03 橡胶 柬埔寨 2013 A1 Cas5 30.2 0 0 0 0
HCcHNDZ1619 橡胶 中国海南儋州 2016 A1 Cas5 40.5 0 0 0 0
HCcHNLD1654 橡胶 中国海南乐东 2016 A1 Cas5 42.7 0 0 0 0
HCcHNQZ1689 橡胶 中国海南琼中 2016 A1 Cas5 13.2 0 0 0 0
CYDB3 橡胶 印度 2007 B1 Cas2 20.3 11.2 26.8 0 0
HCcGDMM1402 橡胶 中国广东茂名 2014 B4 Cas2 33.5 36.8 39.0 27.1 9.2
HCcYNHZM01-3 橡胶 越南 2014 A2 Cas2 19.0 23.7 15.5 7.2 16.3
CmCcHN03 西甜瓜 中国海南三亚 2010 B2 Cas2 23.3 15.3 17.2 21.9 18.6
CmCcHN01 西甜瓜 中国海南三亚 2010 B2 Cas2 16.8 27.1 34.6 38.9 19.5
VsCcHN07 豇豆 中国海南三亚 2010 B2 Cas2 6.5 12.4 30.7 39.2 20.1
VsCcHN09 豇豆 中国海南三亚 2010 B2 Cas2 16.9 28.6 34.1 40.3 19.8
SmCcHN01 茄子 中国海南三亚 2010 B2 Cas2 14.2 7.0 0 6.9 0
MeCcGD02 木薯 中国广东湛江 2007 B3 Cas2 15.8 27.4 32.8 37.1 19.2
CpCcHNQZ1603 番木瓜 中国海南琼中 2016 B1 Cas5 19.9 9.1 25.7 0 0
CpCcYNMX1501 番木瓜 中国云南版纳 2015 B1 Cas5 21.8 11.4 26.1 0 0
LeCcHNSY1702 番茄 中国海南海口 2017 A2 Cas2 21.5 37.9 14.3 6.3 18.2
LJCcHN01 辣椒 中国海南海口 2016 A2 Cas2 18.7 38.8 16.2 6.9 17.7
PaCcSD04 黄瓜 中国山东寿光 2007 B4 Cas2 34.2 36.6 44.3 25.5 8.8
PaCcHB03 黄瓜 中国河北廊坊 2007 B4 Cas2 33.6 37.7 43.9 26.8 9.3
A1致病型的菌株为多基因序列聚类分析第Ⅰ类群的12株橡胶树多主棒孢,来自国内橡胶主产区（云南、海南、广东、广西）以及柬埔寨,这一致病型的菌株只能侵染橡胶树,不能侵染其他寄主,表现出明显的寄主专化性（表3）,进一步验证了多基因聚类分析的结果。A2致病型的橡胶多主棒孢HCcYNHZM01-3（越南）、辣椒多主棒孢LJCcHN01和番茄多主棒孢LeCcHNSY02既能侵染橡胶树,也能侵染其他瓜菜,但对豇豆表现出弱致病力。分属于第 Ⅱ 遗传类群的其他寄主多主棒孢在橡胶、辣椒、黄瓜、豇豆、番茄上能进行相互侵染,其中B1致病型的番木瓜多主棒孢对黄瓜致病力最强,但不侵染豇豆和番茄;B3致病型的木薯多主棒孢不能侵染黄瓜和番茄,对辣椒和豇豆表现出弱致病力;另外,B4致病型的黄瓜多主棒孢（PaCcSD04、PaCcHB03）和橡胶多主棒孢（HCcGDMM1402）除了对番茄表现出弱致病力外,对其他寄主的致病力都很强。
根据致病力测定结果发现,多主棒孢的遗传类群与致病力分化具有明显的相关性,致病力的结果能较好地和多基因序列的聚类结果对应。第 Ⅰ 类群的多主棒孢,在致病型上属于Clade A,毒素类型以Cas5型为主,致病力测定发现,A1致病型只能侵染橡胶树,不能侵染其他寄主,表现出明显的寄主专化性,A2致病型的橡胶多主棒孢既能侵染橡胶树也能侵染其他寄主。第 Ⅱ 类群,主要以其他寄主多主棒孢为主,也有部分橡胶多主棒孢,致病型上属于Clade B,毒素类型均为Cas2型,致病力测定不同寄主来源的多主棒孢之间能进行相互侵染,但致病力和发病症状上存在一定差异。

3 讨论

本研究采用分子系统学的分析方法,利用ITS、EF-1α、TUB等多基因序列分析了我国主要热带作物和果蔬上的70株多主棒孢的遗传种群结构和优势种群,发现不论是单基因还是多基因构建的系统发育树,均能较好的将供试菌株划分为2大遗传类群,第Ⅰ类群有40株菌,均为国内和国外橡胶多主棒孢,国内（云南、海南、广东、广西）的菌株与柬埔寨和越南的菌株高度同源,从地理来源上没有显著的差异。第Ⅱ类群有30株菌,为其他寄主多株棒孢和少数几株橡胶多主棒孢,而且相同寄主的多主棒孢能聚到一起,遗传进化上具有高度的相似性。说明供试的多主棒孢与寄主来源关系密切,但与地理来源不具有显著的相关性,这与Qi等[4]利用RFLP和RAPD以及Dixon 等[6]利用肌动蛋白、ITS区以及2个多态性位点的分析结果一致。
致病力分化分析结果显示,致病型A1的橡胶多主棒孢只能侵染橡胶树不能侵染其他寄主,具有明显的寄主专化性,但不同地理来源的橡胶多主棒孢在致病力上没有显著的差异。致病型A2和B1的菌株既能侵染橡胶也能侵染其他瓜菜寄主,不如A1型的致病力强。致病型B3和B4的多主棒孢菌株之间均能进行交互侵染,但菌株间的致病力存在一定的差异,表现为部分多主棒孢菌株对其分离寄主的致病力往往高于其他寄主,而有些菌株对非自身寄主之外的其他寄主表现出弱致病力或不致病。说明多主棒孢菌尽管寄主范围十分广泛,但对其原寄主的致病力是最强的,这与多主棒孢与其原寄主植物可能存在协同进化的关系,在寄主长期作用的影响下,也导致多主棒孢具有寄主专化性[14,15]。另外,发现采集自相同地点、相同寄主的瓜菜（如西甜瓜和豇豆）多主棒孢在致病力上也存在较大的差异,推测在同一个生态系统中这些菌株在个体之间可能存在生理生化和遗传分化上的差异,这导致病原菌对寄主的侵染力、定殖、症状以及危害程度上的差异。因此,对多主棒孢主要优势种群和主要致病型的确定,能为多主棒孢的流行特点分析、抗病品种的选育以及防控策略的制定奠定重要的理论依据。
总体来看,结合多基因序列聚类分析和致病力测定的结果,能较好地反应多主棒孢的遗传结构、优势种群以及主要致病型等,而且多基因序列聚类分析与致病力分化分析结果能相互印证,具有较好的一致性。就橡胶树多主棒孢而言,第Ⅰ类群是我国橡胶树多主棒孢的优势种群具有明显的寄主专化性,致病型A1是其主要的致病型,结合橡胶树棒孢霉落叶病在我国的发生历史[16],不排除橡胶多主棒孢是从我国边境的植胶国通过种苗调运、贸易往来等途径传入国内的可能性,但仍需要更多地理来源的菌株和大量的证据进行证实。就其他热带果蔬多主棒孢而言,虽然不像橡胶多主棒孢具有明显的寄主专化性,但不同寄主间多主棒孢的致病力存在一定的差异,表现出对其原寄主的高度致病力,加上寄主范围十分广泛,在不同寄主上表现的症状也具有显著的差异,本研究基本查明了这些热带果蔬多主棒孢的种群结构和主要致病型,为多主棒孢病害的发生、流行和防治提供了理论依据。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
Rahman A Y A, Usharraj A O, Misra B B , et al. Draft genome sequence of the rubber tree Hevea brasiliensis[J]. BMC Genomics, 2013,14(1):75.

DOI PMID

[2]
李宝聚, 赵彦杰, 于淑晶 , 等. 2008 年秋季河北青县黄瓜棒孢叶斑病大发生[J]. 中国蔬菜, 2008(11):51-52, 66.

[3]
Silva W P K, Deverall B J, Lyon B R . Molecular, physiological and pathological characterization of Corynespora leaf spot fungi from rubber plantations in Sri Lanka[J]. Plant Pathology, 1998,47(3):267-277.

DOI

[4]
Qi Y X, Zhang X, Pu J J , et al. Morphological and molecular analysis of genetic variability within isolates of Corynespora cassiicola from different hosts[J]. European Journal of Plant Pathology, 2011,130(1):83-95.

DOI

[5]
Smith L J, Datnoff L E, Pernezny K , et al. Phylogenetic and pathogenic characterization of Corynespora cassiicola isolates[J]. Acta Horticulturae, 2009(808):51-56.

DOI PMID

[6]
Dixon L J, Schlub R L, Pernezny K , et al. Host specialization and phylogenetic diversity of Corynespora cassiicola[J]. Phytopathology, 2009,99(9):1015-1027.

DOI PMID

[7]
Shimomoto Y, Sato T, Hojo H , et al. Pathogenic and genetic variation among isolates of Corynespora cassiicola in Japan[J]. Plant Pathology, 2010,60(2):253-260.

DOI PMID

[8]
Déon M, Fumanal B, Gimenez S , et al. Diversity of the cassiicolin gene in Corynespora cassiicola and relation with the pathogenicity in Hevea brasiliensis[J]. Fungal Biology, 2014,118(1):32-47.

DOI

[9]
Fang Z D. Research method for plant pathology [M]. 2nd ed. Beijing: China Agriculture Press, 1998: 63.

[10]
刘先宝, 蔡吉苗, 潘羡心 , 等. 橡胶树多主棒孢菌rDNA-ITS区的分子鉴定及检测[J]. 热带作物学报, 2008,29(4):489-493.

[11]
阚琳娜, 李宝聚, 纪明山 , 等. 黄瓜褐斑病防治药剂的活体筛选[J]. 中国蔬菜, 2007,1(4):22-24.

[12]
李博勋, 刘先宝, 林春花 , 等. 国内橡胶树主要种质对棒孢霉落叶病的抗性评价[J]. 植物保护, 2014,40(5):86-92.

[13]
刘先宝, 李博勋, 陈珊 , 等. 国内橡胶树多主棒孢Cassiicolin毒素多样性及致病性分析[J]. 热带作物学报, 2016,37(10):1969-1973.

[14]
高苇 . 我国黄瓜棒孢叶斑病的病原学及诊断技术的研究[D]. 北京: 中国农业科学院, 2011.

[15]
Cutrim F A, Silva G S . Pathogenicity of Corynespora cassiicola to different plant species[J]. Fitopatologia Brasileira, 2003,28(2):193-194.

DOI PMID

[16]
Pu J J, Zhang X, Qi Y X , et al. First record of Corynespora leaf fall disease of Hevea rubber tree in China[J]. Australasian Plant Disease Notes, 2007,2(1):35-36.

DOI

Options
文章导航

/

〈 〉

    联系我们

  • 地址：海南省海口市龙华区学院路4号
  • 邮编：571101
  • 电话：0898-66989102
    　　    0898-66890382
               0898-66986832
  • E-mail：rdzx136@163.com
网站版权 © 《热带作物学报》编辑部 　琼ICP备19001573号