我国苹果斑点落叶病菌携带dsRNA分析及一种dsRNA病毒的鉴定

曹钰晗，李紫腾，张静怡，张静娜，胡同乐，王树桐，王亚南，曹克强

我国苹果斑点落叶病菌携带dsRNA分析及一种dsRNA病毒的鉴定

曹钰晗，李紫腾，张静怡，张静娜，胡同乐，王树桐，王亚南，曹克强

河北农业大学植物保护学院，河北保定 071001

【目的】苹果斑点落叶病在我国各苹果主产区均有发生，给苹果产业造成了严重的经济损失。苹果斑点落叶病是由链格孢苹果专化型（f. sp.）侵染引起的一种气传病害。本研究旨在探明我国苹果斑点落叶病菌携带dsRNA的情况，为田间防治斑点落叶病提供新的生防资源，并为病毒的多样性以及进化研究提供新的认识。【方法】从全国8个省份采集有典型症状的组织样本，通过组织分离和单孢分离法获得纯培养；通过dsRNA的提取及凝胶电泳明确我国苹果斑点落叶病菌携带dsRNA的情况；通过菌株培养特性、菌丝生长速率、在果实及叶片上的致病力测定揭示携带dsRNA病原菌的生物学性状；通过高通量测序技术、分子克隆技术对QY-2菌株携带的dsRNA病毒进行全基因组序列、基因组结构和系统进化分析。【结果】自我国苹果产区获得102株苹果斑点落叶病菌菌株的纯培养，通过dsRNA的提取及凝胶电泳明确了5个菌株带有明显的dsRNA条带，携带的dsRNA分为4种类型，类型I（菌株YT-3-7）dsRNA在8、2.5和1.5 kb左右；类型II（菌株SJZ-4）dsRNA在8 kb左右；类型III（菌株QY-2）dsRNA在3和0.8 kb左右；类型IV（菌株CL-2-6和SQ-1-1）dsRNA在2.5和1.5 kb左右。含有dsRNA的菌株培养性状多样，菌落特征、生长速率与dsRNA携带与否及类型没有明显的关系，含有dsRNA的菌株大多属于弱致病力菌株。QY-2在叶片和果实上的致病力均最弱，确定其携带的dsRNA病毒基因组包括5个dsRNA片段，分别为dsRNA1—dsRNA5，片段大小依次为3 665、3 054、2 824、2 819、831 nt，提交至GenBank，登录号分别为MK672910、MK672913、MK672912、MK672911、MK836314。编码蛋白的分子量依次为124、83、84、83、13 kD。经系统进化关系分析，与chrysovirus 1遗传关系最近，确定其为产黄青霉病毒科（）、产黄青霉病毒属（）的chrysovirus，暂定命名为chrysovirus 2（AaCV2）。【结论】我国苹果斑点落叶病菌携带dsRNA群体多样，dsRNA的有无及类型与寄主病原菌培养性状没有明显相关性，携带dsRNA的菌株多为弱致病力菌株，致病力最低的QY-2菌株携带AaCV2。该菌株的弱致病力可能具有潜在的应用价值，可为苹果斑点落叶病提供新的生防资源。

苹果斑点落叶病；链格孢苹果专化型；dsRNA；生物学性状；致病力；生物防治

0 引言

【研究意义】苹果斑点落叶病（apple Alternaria blotch）又称褐纹病，是由链格孢苹果专化型（f. sp.）侵染而引起的一种世界性分布气传病害，该病危害苹果叶片、新梢和果实，严重影响苹果生产[1-2]，现已成为亚洲国家苹果产区的主要病害[3-4]。目前，苹果斑点落叶病的防治途径主要包括种植抗性品种、化学防治和生物防治，生产上一直以使用化学药剂防治为主[5-7]，但随着经济的发展，人们开始关注食品安全和环境问题，一些化学杀菌剂的使用已在发达国家和地区受到了严重限制[8-9]，生物防治日益受到人们的关注，利用携带病毒的弱毒株引起病原菌群体致病力衰退，将成为真菌病害生物防治的新途径[10-11]。【前人研究进展】真菌病毒是指能够侵染真菌并在其中完成自我复制的病毒。真菌病毒的分类依据主要有两个，首先依据真菌病毒的核酸类型将其分为DNA病毒和RNA病毒，就目前所发现的真菌病毒而言，大多数为RNA病毒，包括双链RNA（dsRNA）和单链RNA（ssRNA）[12-13]。从具有弱毒特性的病毒入手寻求新的生防因子，对深入了解病毒及其对寄主真菌的影响具有重要意义。链格孢属真菌中也发现了很多真菌病毒，1988年，Shepherd[14]从棉花种子中分离出链格孢菌株共21个，其中6株携带真菌病毒，且这6株均不产生色素；1997年，ZABALGOGEAZCOA等[15]从茄链格孢（）中检测到未经包被的两条大小分别为8.3和5.5 kb的dsRNA；1988年，Hayashi等[16]研究的12个链格孢日本梨致病型（Japanese pear pathotype）菌株中，有7个菌株含有dsRNA；2009年，Aoki等[17]从链格孢菌株EGS 35-193中分离到virus 1（AaV1），该病毒使菌株的菌丝生长变得异常，并且出现大量的囊泡结构；2015年，LIN等[18]从长柄链格孢（）HN28中发现dsRNA virus 1（AlRV1），其分类地位尚未明确；2016年，Komatsu等[19]从乔木链格孢（）中发现一个新的维多利亚病毒victorivirus 1（AaVV1）；ZHONG等[20]从芸薹生链格孢（）中发现病毒fusarivirus 1（AbFV1）；2017年，向均[21]从链格孢日本梨致病型中分离出病毒virus 1（ABRV1）；2018年，日本学者OKADA等[22]从链格孢日本梨致病型中分离出病毒chrysovirus 1（AaCV1）。AaCV1损害宿主真菌的生长，同时影响宿主AK毒素的产生。【本研究切入点】目前，已报道了多种作物链格孢属真菌病毒，但关于苹果斑点落叶病菌真菌病毒还未见报道，有待进一步解析。【拟解决的关键问题】在我国苹果产区8个省份采集苹果斑点落叶病样本，通过dsRNA提取，明确我国苹果斑点落叶病菌携带dsRNA的多样性；通过生物学测定揭示携带dsRNA菌株的生物学特征，分析dsRNA多样性与寄主病原菌生物学性状的关系；最后，鉴定弱毒力菌株携带dsRNA病毒的种类。为寻找对苹果斑点落叶病具有生防潜力的弱毒菌株打下基础。

1 材料与方法

1.1 供试材料

102株苹果斑点落叶病菌菌株：2018年5—8月自我国河北、河南、山东、甘肃、辽宁、黑龙江、陕西及云南8个省份采集有典型斑点落叶症状的叶片样本，采用组织分离法和单孢子稀释法进行病原菌的分离纯化[23]。在叶片和果实上利用柯赫氏法则进行验证。具体样本信息见表1。

表1 本研究中苹果斑点落叶病样本信息

1.2 主要试剂、仪器和引物

1.2.1 主要试剂和仪器 氯仿，北京酷莱博科技有限公司；纤维素粉，阿拉丁试剂（上海）有限公司；3M乙酸钠、氨苄，北京索莱宝科技有限公司；Taq Mix、载体PMD18-T，宝生物工程有限公司；感受态细胞Tans1-T、TRlzol、反转录试剂盒，北京全式金生物科技有限公司；DNA胶回收试剂盒，生工生物工程（上海）股份有限公司；RNeasy Mini Kit总RNA提取试剂盒，杭州沃森生物技术有限公司。组织研磨破碎仪，德国QIAGEN公司；台式冷冻高速离心机，默克生命科学（上海）有限公司；超低温冰箱、PCR仪，赛默飞世尔科技；DYY-12型电泳仪、DYY-31DN型琼脂糖水平电泳槽，北京六一仪器厂；BioRad凝胶成像系统，美国SYNGENE公司；超净工作台，江苏苏净集团公司。

1.2.2 引物 本研究所用引物具体信息见表2。

表2 本试验所用引物

1.3 dsRNA的提取

参考李波[24]的报道，采用纤维素粉的方法对102株苹果斑点落叶病菌纯化菌株进行dsRNA提取，最后加入无水乙醇和乙酸钠（3 mol·L-1）在-80℃沉淀，12 000 r/min离心后，用1 mL的70%乙醇洗涤沉淀，20 μL DEPC H2O溶解，取8 μL在1%（W/V）的琼脂糖凝胶电泳中检测（电压80 V，50 min），紫外灯下观察dsRNA电泳条带。

1.4 苹果斑点落叶病菌带毒菌株的生物学特性

1.4.1 培养特性及菌丝生长速率的测定 将5株带毒菌株及随机选取的20株不带毒菌株同时扩繁，25℃PDA培养7 d后，观察菌落的形状、颜色，每个菌株设置3个重复。用十字交叉法测菌落直径，计算生长速率。

1.4.2 在苹果叶片和果实上致病力的测定 采集海棠幼嫩无病叶片，清水洗净，放于铺有两层浸湿滤纸的瓷盘中，使用接种针将叶片的叶脉一侧刺伤；从25个苹果斑点落叶病菌菌株的PDA平板相同位置的边缘打取菌饼，将菌丝块倒置于叶片的伤口处，每个菌株设置3个重复；用保鲜膜将瓷盘包裹严实，置于25℃条件下培养7 d后，测量病斑大小。

取表观相近的同一批次健康富士苹果，先清水冲洗，再用75%的酒精擦洗干净，75%酒精淋湿，放于超净工作台中吹干。用接种针将果实表面刺一个针孔，从25个苹果斑点落叶病菌菌株的PDA平板相同位置边缘打取菌饼，有菌丝的一面贴于伤口处，保鲜膜包裹保湿，每个菌株至少3个重复。25℃培养7 d后测量病斑大小。

1.4.3 数据处理 数据采用SPSS Statistics 17.0单因素分析法中的均值比较法进行分析。

1.5 菌株QY-2携带病毒的高通量测序

1.5.1 RNA的抽提与质检 采用RNeasy Mini Kit进行菌株QY-2的total RNA抽提，然后通过Agilent Bioanalyzer 2100电泳分析仪进行RNA质量检测，合格后纯化total RNA。

1.5.2 文库构建与测序 委托上海伯豪生物技术有限公司采用宏转录组测序的方法，对纯化后的total RNA进行rRNA去除、片段化、第一链cDNA合成、第二链cDNA合成、末端修复、3′末端加A、连接接头、富集等步骤，完成测序样本文库构建，采用paired-end程序，进行双端测序。测序过程由Illumina提供的data collection software进行实时数据分析。

1.5.3 QY-2携带病毒种类分析及基因组序列分析 通过高通量测序结果分析QY-2菌株携带的病毒种类及占比。将样本测序获得的reads进行reference mapping，获得病毒reads覆盖度，用RPKM（reads per kilo bases per million reads）量化标准对reads覆盖度进行表达量计算，计算公式：transcription reads /（transcription length×total assembly reads in run）其中，transcription reads为覆盖病毒整个UniGene的reads数目，transcription length为UniGene长度，total assembly reads in run为该样本所有参与拼接的总reads数目。通过Vector NTI 11.0软件去除载体序列和序列比对；NCBI上预测开放阅读框（ORF）；运用MEGA 7.0中的邻接法进行序列系统发育分析，重复1 000次。

1.6 菌株QY-2携带病毒的RT-PCR鉴定

1.6.1 病毒特异性引物的设计 根据已报道的病毒序列，使用Premier5.0设计该病毒的特异性引物，委托生工生物工程（上海）股份有限公司进行合成。引物具体信息见表2。

1.6.2 总RNA的提取 取适量纯化好的新鲜菌丝，液氮研磨成粉末状，转移至2 mL离心管中，加入1 mL的Trizol试剂，颠倒混匀，室温下放置5 min；继续加入200 μL的氯仿，剧烈振荡15 s，室温下放置3 min；4℃，14 000 r/min离心15 min；取上清700 μL到1.5 mL的离心管中，加入700 μL的异丙醇，颠倒混匀，室温下放置10 min；4℃，14 000 r/min离心10 min，弃上清，吸除残留；加入1 mL 75%的乙醇洗涤沉淀；4℃，9 000 r/min离心5 min，弃上清，短暂离心后吸除残留的乙醇；室温干燥5—10 min，加入50 μL的DEPC-H2O。

1.6.3 反转录合成cDNA 使用全式金反转录试剂盒合成cDNA，操作方法参考说明书。具体步骤如下：dsRNA样品2 μL；Random Primer（N9）（0.1 μg·μL-1）1 μL；2×ES Reaction Mix 10 μL；Easy Script RT Mix 1 μL；DEPC-H2O 6 μL。25℃反应10 min，42℃反应30 min，然后85℃反应5 min，即为cDNA。

PCR反应体系：cDNA样品，1 μL；2×taq mix，10 μL；引物F 1 μL；引物R 1 μL；DEPC-H2O 7 μL。PCR反应程序：94℃预变性5 min，94℃变性30 s，58℃退火30 s，72℃延伸45 s，35个循环后，72℃完全延伸7 min。1%琼脂糖凝胶电泳检测目的片段。

1.6.4 目的片段的克隆测序 采用琼脂糖凝胶回收试剂盒回收DNA，按试剂盒说明书进行操作。按照pMD18-T载体试剂盒说明书进行载体连接，利用感受态细胞Tans1-T进行转化，挑取单菌落到含有氨苄抗生素的液体LB中，37℃，200 r/min振荡至液体浑浊。以菌液为模板，菌落PCR鉴定阳性克隆，阳性克隆委托华大基因测序公司进行测序。

2 结果

2.1 苹果斑点落叶病菌的dsRNA多态性分析

对102株苹果斑点落叶病菌菌株进行dsRNA的提取，使用DNase I及S1核酸酶消解，在102个菌株中有5个菌株携带dsRNA，如图1所示，菌株名称分别为CL-2-6、QY-2、SQ-1-1、SJZ-4、YT-3-7。dsRNA分为4种类型：菌株YT-3-7携带的dsRNA条带大小在8、2.5和1.5 kb左右；菌株SJZ-4携带的dsRNA条带在8 kb左右；菌株QY-2携带的dsRNA条带在3和0.8 kb左右；菌株CL-2-6和SQ-1-1携带的dsRNA条带在2.5和1.5 kb左右。

2.2 苹果斑点落叶病菌带毒菌株的培养性状

将苹果斑点落叶病菌带毒菌株和不带毒菌株在PDA培养基上25℃黑暗培养7 d，观察菌落颜色、气生菌丝形态等特征，将5个带毒菌株划分为两类。

A：酶处理前dsRNA Untreated dsRNA；B：酶处理后dsRNA Enzyme-treated dsRNA

第1类（图2）：菌丝灰白色或灰绿色，个别菌落后期有黑色素产生。菌丝分布均匀，多数菌株的菌丝层较薄，个别菌株的气生菌丝发达呈棉絮状。具有这一类特征的菌株有12个，包括两个带毒菌株SJZ-4和CL-2-6以及10个不带毒菌株，分别为HLJ-2-1、HLJ-2-3、CL-4-1、QY-5、QY-4、YT-3-9、YT-4-3、HLJ-1-6、HLJ-2-9、SQ-1-2。

第2类（图3）：菌丝前期白色，后期呈绿色，多数菌株内部的绿色菌丝和外围的白色菌丝间形成明显的环状界限。具有这一类特征的菌株有13个，包括3个带毒菌株，分别为QY-2、YT-3-7、SQ-1-1，以及10个不带毒菌株，分别为SQ-3-14、QY-8、CL-4-2、SQ-3-6、CL-2-9、CL-4-9、SMX-1-3、SQ-1-6、YA、LN-9。

2.3 苹果斑点落叶病菌带毒菌株菌丝的生长速率

苹果斑点落叶病菌不同菌株的菌丝生长速率测定结果如图4所示，带毒菌株与非带毒菌株相关性差异不明显，各菌株间生长速度的快慢因菌株而异。

2.4 苹果斑点落叶病菌带毒菌株的致病力

将带毒菌株和20个不带毒菌株接种海棠叶片，测定各菌株的致病力。测量接种第7天的病斑大小（图5），发现在5个带毒菌株中QY-2和YT-3-7的致病力较弱，其余菌株致病力中等。在20个不带毒的菌株中LN-9和CL-2-9的致病力较强，其他菌株的致病力中等（图6）。

将带毒菌株和20个不带毒菌株接种富士苹果，测量第7天的病斑大小，发现在5个带毒菌株中QY-2为弱致病力菌株，SJZ-4的致病力较强，其余菌株致病力中等。在20个不带毒的菌株中LN-9、CL-2-9、CL-4-9、YA、YT-3-9、HLJ-2-9、HLJ-1-6、SQ-3-14的致病力较强，其他菌株致病力中等（图7）。

带毒菌株QY-2在叶片和果实上的致病力均较弱，从寻找具有弱毒特性的菌株方面考虑，该菌株可能具有潜在的应用价值，选为进一步研究的对象。

2.5 苹果斑点落叶病菌菌株QY-2中dsRNA病毒的鉴定

2.5.1 病毒的全基因组序列 通过分析QY-2菌株中dsRNA高通量测序结果，获得了4条与真菌病毒相关的序列，经过Blast表明4个contig均与来源于链格孢日本梨致病型的AaCV1序列有较高的相似性，以AaCV1的4个dsRNA片段序列为参考序列，有1 451 570、2 331 542、1 167 609、1 268 863个reads与AaCV1 dsRNA片段序列相符，经比对和拼接获得了Contig 10、Contig 41、Contig 38、Contig 24，对reads覆盖度进行表达量计算，4个Contig的RPKM值在所有基因中是最高的，达到了4 820—9 292。E-value数值很小，说明了结果的可靠性（表3）。

A、B：带毒菌株The strains taking virus (SJZ-4, CL-2-6)；C—L：未带毒菌株The strains without virus (HLJ-2-1, HLJ-2-3, CL-4-1, QY-5, QY-4, YT-3-9, YT-4-3, HLJ-1-6, HLJ-2-9, SQ-1-2)

A—C：带毒菌株The strains taking virus (QY-2, YT-3-7, SQ-1-1)；D—M：未带毒菌株The strains without virus (SQ-3-14, QY-8, CL-4-2, SQ-3-6, CL-2-9, CL-4-9, SMX-1-3, SQ-1-6, YA, LN-9)

图4 苹果斑点落叶病菌不同菌株的菌丝生长速率

A：CK；B—F：接种带毒菌株Inoculation of strains taking virus (YT-3-7, SJZ-4, QY-2, SQ-1-1, CL-2-6)；G—Z：接种不带毒菌株Inoculation of strains without virus (LN-9, SQ-1-2, SQ-1-6, CL-2-9, CL-4-1, CL-4-2, CL-4-9, QY-5, QY-4, QY-8, SMX-1-3, YA, YT-3-9, YT-4-3, HLJ-2-3, HLJ-2-1, HLJ-2-9, HLJ-1-6, SQ-3-6, SQ-3-14)

图6 苹果斑点落叶病菌不同菌株在海棠叶片上的致病力

图7 苹果斑点落叶病菌不同菌株在苹果果实上的致病力

表3 高通量测序苹果斑点落叶病菌菌株QY-2携带病毒信息

E-value描述同源性的可靠性，值越小结果越可靠。E-value低于1.0E-5认为序列显著性同源

E-value describes the reliability of homology. The smaller the value, the more reliable the result is. E-value lower than 1.0E-5 is considered as significant homology

已报道的AaCV1基因组具有5条dsRNA片段，而高通量测序结果只获得了病毒的4条dsRNA序列，QY-2 dsRNA琼脂糖凝胶电泳（图1）显示在800 bp左右处也存在一条dsRNA，根据已报道的AaCV1设计其第5个片段的特异性引物CRV1-5进行分子克隆测序，最终获得了菌株QY-2 dsRNA病毒的第5个片段序列。综上，菌株QY-2 dsRNA病毒包括5个dsRNA片段，分别为dsRNA1—dsRNA5，片段大小依次为3 665、3 054、2 824、2 819、831 nt。5条dsRNA序列已全部提交到GenBank，登录号分别为MK672910、MK672913、MK672912、MK672911、MK836314。

2.5.2 病毒的基因组结构分析 将dsRNA1—dsRNA5的序列在NCBI网站上进行ORF分析，结果显示dsRNA1—dsRNA5分别编码一个推定的ORF，依次为ORF1—ORF5。 ORF1—ORF5依次含有1 117、774、774、769、115 aa，编码蛋白的分子量依次为124、83、84、83、13 kD。dsRNA1—dsRNA5的5′端非翻译区（UTR）大小分别为205、573、357、401、102 nt；3′ UTR大小分别为107、157、133、109、382 nt（图8）。

2.5.3 病毒开放阅读框序列相似性分析 将QY-2菌株dsRNA病毒ORF1—ORF5编码的氨基酸序列在NCBI上进行blastp搜索，结果显示与其相似性相对较高的病毒有4个，分别是AaCV1、chrysovirus1（BdCV1）、chrysovirus 1（AtCV1）、chrysovirus 2（PjCV2）（表4）。与QY-2菌株dsRNA病毒ORF1（RdRP）氨基酸序列相似性分别为93%、56%、55%、52%，核苷酸序列相似性为77%、59%、57%、54%；ORF2氨基酸序列相似性为91%、47%、44%、40%，核苷酸序列相似性为75%、54%、50%、49%；ORF3氨基酸序列相似性为80%、35%、32%、32%，核苷酸序列相似性为70%、40%、42%、42%；ORF4氨基酸序列相似性为90%、47%、43%、43%，核苷酸序列相似性为76%、51%、50%、56%；与ORF5氨基酸序列相似性较高的只有AaCV1编码的假定蛋白，相似性为87%，核苷酸序列相似性为95%。根据国际病毒分类委员会（ICTV）2020年发布的内种的分类标准（https://talk.ictvonline.org/taxonomy/），从寄主种类、RdRP、CP氨基酸序列相似性（≤70%、≤53%）、dsRNA片段数目、5′-UTR长度几个方面分析，确定菌株QY-2所携带的dsRNA病毒为产黄青霉病毒科（）、产黄青霉病毒属（）的chrysovirus。

图8 苹果斑点落叶病菌菌株QY-2 dsRNA病毒基因组模式图

表4 苹果斑点落叶病菌菌株QY-2的dsRNA与其他病毒的核苷酸及氨基酸序列的相似性

2.5.4 病毒系统进化关系分析 为分析病毒系统进化关系，将苹果斑点落叶病菌菌株QY-2携带的AaCV2 RdRP氨基酸序列与、以及一些未分类病毒的RdRP氨基酸序列构建系统发育树。结果表明，菌株QY-2携带的病毒与的聚在一个分支，并且与病毒AaCV1遗传关系最近，其自展支撑值为100（图9）。因此将QY-2菌株携带的病毒暂定命名为chrysovirus 2（AaCV2）。

2.6 菌株QY-2携带病毒的RT-PCR鉴定

根据获得的QY-2 dsRNA病毒基因序列，设计特异性引物CRV1-1对其进行RT-PCR检测。琼脂糖凝胶电泳显示，菌株QY-2出现了453 bp的病毒特异性目的条带，对照菌株及清水对照均没有条带出现，进一步说明了菌株QY-2中AaCV2的存在（图10）。

3 讨论

目前，苹果斑点落叶病防治的主要方式为化学防治。为避免化学防治对人体和环境产生的不利影响，寻求安全有效的绿色生防措施颇为重要。部分真菌病毒可导致真菌致病力的衰退，因此可作为一种生防因子用于真菌病害的防治[25-26]。

图9 以真菌病毒的RdRP氨基酸序列为基础构建的系统发育树

M：DL2000；1：QY-2；2：对照菌株Control strain；3：清水对照Water control

本研究对采自我国8个省份苹果主产区的102株苹果斑点落叶病菌菌株进行dsRNA检测，根据琼脂糖凝胶电泳结果，明确了其中5个菌株具有明显的dsRNA条带，可能携带真菌病毒，约占5%的菌比例，在苹果斑点落叶病菌中dsRNA具有一定的普遍性，且根据条带大小来看，病毒种类具有多样性。据报道链格孢属真菌目前发现的病毒包括：从棉花种子中分离出6个携带真菌病毒的链格孢菌株，dsRNA片段在1.0—5.1 kb[14]；从茄链格孢中检测到两条dsRNA大小分别为8.3和5.5 kb[15]；另外还有真菌病毒AaV1（dsRNA片段：3 617、2 794、2 576、1 420 bp）[17]、AlRV1（dsRNA片段：3 145 bp）[18]、AaVV1（dsRNA片段：5 203 bp）[19]、AbFV1（ssRNA片段：6 639 bp）[20]、ABRV1（dsRNA片段：6 188、5 903 bp）[21]、AaCV1（dsRNA片段：3 647、2 857、2 785、2 772、836 bp）[22]。不同的链格孢菌株携带dsRNA条带差异较大，且与本研究所分离到的dsRNA病毒条带有较大差异。因此，本研究发现的苹果斑点落叶病菌真菌病毒有可能是未揭示的新病毒。

随着真菌病毒种类及数量不断增多，人们开始研究病毒对真菌生物学功能方面的影响，发现大多数病毒在菌株中能够稳定遗传[27-28]，但因病毒种类不同，对真菌致病力及其他方面会有不同的影响。真菌病毒与寄主间的相互作用成为现代病毒学研究的热点领域。本研究主要针对携带dsRNA病毒的5个菌株进行了形态学观察、生长速率测定及致病力鉴定，同时与不携带dsRNA病毒的其他菌株进行比较，以明确携带病毒是否引起了菌株的特异性或其是否具有生防意义。通过测定发现，带毒菌株间菌落形态、菌丝的生长速率和菌株的致病能力存在显著差异，不带毒菌株间这3个方面同样也存在显著差异。菌株的菌落特征、生长速率与dsRNA携带与否及类型没有明显的关系，致病力测定表明含有dsRNA病毒的供试菌株大多属于低毒力类型，其中的线性关系或多因素关系，还有待进一步研究。

在含有dsRNA病毒的低致病力菌株中，QY-2致病力最低，与清水对照无显著性差异。经鉴定该菌株中携带大量的产黄青霉病毒科、产黄青霉病毒属的chrysovirus，这是苹果斑点落叶病菌中该病毒的首次报道。王利华等[29]明确了产黄青霉病毒科成员Chrysovirus 1（BdCV1）对梨轮纹病菌（）的菌落形态、生长速度有明显抑制作用，是引起梨轮纹病菌致病力衰退的主要因子；Ma等[30]报道细链格孢（）携带AaCV1-AT1，AaCV1-AT1的感染降低了宿主真菌的菌落生长速度和产孢能力；Okada等[22]报道了AaCV1对链格孢日本梨致病型具有两个方面的作用，一是影响真菌宿主的生长速度，二是影响真菌宿主毒素的产生。因此，AaCV2也可能是引起QY-2菌株致病力衰退的主要原因。下一步将对QY-2携带的AaCV2进行深入研究，以开发对苹果斑点落叶病具有生防潜力的dsRNA病毒。

4 结论

我国苹果斑点落叶病菌菌株携带dsRNA的概率在5%左右，dsRNA群体多样，有4种类型：类型I的dsRNA在8、2.5和1.5 kb左右；类型II的dsRNA在8 kb左右；类型III的dsRNA在3和0.8 kb左右；类型IV的dsRNA在2.5和1.5 kb左右。dsRNA的有无及类型与寄主病原菌培养性状没有明显相关性，携带dsRNA的菌株多为弱致病力，致病力最低的QY-2菌株携带产黄青霉病毒科、产黄青霉病毒属的chrysovirus 2。

[1] JUNG K H. Growth inhibition effect of pyroligneous acid on pathogenic fungus,, the agent of alternaria blotch of apple. Biotechnology and Bioprocess Engineering, 2007, 12: 318-322.

[2] 胡同乐, 曹克强, 王树桐, 甄文超. 生长季苹果斑点落叶病流行主导因素的确定. 植物病理学报, 2005, 35(4): 374-377.

HU T L, CAO K Q, WANG S T, ZHEN W C. Study on the main factor for the epidemic of alternaria blotch on apple. Acta Phytopathologica Sinica, 2005, 35(4): 374-377. (in Chinese)

[3] MIYAMOTO Y, ISHII Y, HONDA A, MASUNAKA A, TSUGE T, YAMAMOTO M, OHTANI K, FUKUMOTO T, GOMI K, PEEVER T L, AKIMITSU K. Function of genes encoding Acyl-CoA synthetase and enoyl-CoA hydratase for host-selective ACT-toxin biosynthesis in the tangerine pathotype of. Phytopathology, 2009, 99(4): 369-377.

[4] LI Y, ZHANG L Y, ZHANG Z, CONG P H, CHENG Z M. A simple sequence repeat marker linked to the susceptibility of apple to alternaria blotch caused byapple pathotype. Journal of the American Society for Horticultural Science, 2011, 136(2): 109-115.

[5] 吴玉星, 李美娜, 周宗山, 仇贵生, 金纪艳, 乔壮. 8种杀菌剂防治苹果斑点落叶病试验. 中国果树, 2007(2): 28-30.

WU Y X, LI M N, ZHOU Z S, QIU G S, JIN J Y, QIAO Z. Experiment of eight fungicides against apple Alternaria blotch. China Fruits, 2007(2): 28-30. (in Chinese)

[6] 任璐, 史晓晶, 姚众, 韩巨才, 赵晓军. 苹果斑点落叶病菌对戊唑醇敏感基线建立及抗性突变体适合度. 植物病理学报, 2017, 47(3): 380-388.

REN L, SHI X J, YAO Z, HAN J C, ZHAO X J. Baseline sensitivity off. sp.to tebuconazole and fitness of resistant mutants. Acta Phytopathologica Sinica, 2017, 47(3): 380-388. (in Chinese)

[7] 赵国康, 李焰, 张树武, 徐秉良, 刘佳. 5种植物源农药对苹果斑点落叶病的防效评价. 中国果树, 2020(5): 46-49.

ZHAO G K, LI Y, ZHANG S W, XU B L, LIU J. Evaluation the control effects of five botanical fungicides on apple Alternaria leaf spot (). China Fruits, 2020(5): 46-49. (in Chinese)

[8] FERNANDO W G D, NAKKEERAN S, ZHANG Y, SAVCHUK S. Biological control of(Lib.) de Bary byandspecies on canola petals.Crop Protection, 2007, 26(2): 100-107.

[9] ANAGNOSTAKIS S L. Biological control of chestnut blight. Science, 1982, 215(4532): 466-471.

[10] NUSS D L. Biological control of chestnut blight: an example of virus-mediated attenuation of fungal pathogenesis. Microbiological Reviews, 1992, 56(4): 561-576.

[11] HOLLINGS M. Viruses associated with a die-back disease of cultivated mushroom. Nature, 1962, 196(4858): 962-965.

[12] GHABRIAL S A, NIBERT M L., a new genus of fungal viruses in the family. Archives of Virology, 2009, 154(2): 373-379.

[13] YU J, KWON S J, LEE K M, SON M, KIM K H. Complete nucleotide sequence of double-stranded RNA viruses fromstrain DK3. Archives of Virology, 2009, 154(11): 1855-1858.

[14] SHEPHERD H S. Viruslike particles in tentoxin-producing strains of. Journal of Virology, 1988, 62(10): 3888-3891.

[15] ZABALGOGEAZCOA I, PETRUNAK D, CHRIST B J, GILDOW F E. Unencapsidated double-stranded RNA associated with membrane vesicles in isolates of. Mycological Research, 1997, 101(5): 604-608.

[16] HAYASHI N, TSUGE T,KOBAYASHI H, NISHIMURA S. The presence of double-stranded RNAs inJapanese pear pathotype and their participation in AK-toxin productivity. Japanese Journal of Phytopathology, 1988, 54(2): 250-252.

[17] AOKI N, MORIYAMA H, KODAMA M, ARIE T, TERAOKA T, FUKUHARA T. A novel mycovirus associated with four double- stranded RNAs affects host fungal growth in. Virus Research, 2009, 140(1/2): 185-187.

[18] LIN Y H, ZHANG H L, ZHAO C J, LIU S X, GUO L H. The complete genome sequence of a novel mycovirus fromstrain HN28. Archives of Virology, 2015, 160(2): 577-580.

[19] KOMATSU K, KATAYAMA Y, OMATSU T, MIZUTANI T, FUKUHARA T, KODAMA M, ARIE T, TERAOKA T, MORIYAMA H. Genome sequence of a novel victorivirus identified in the phytopathogenic fungus. Archives of Virology, 2016, 161(6): 1701-1704.

[20] ZHONG J, SHANG H H, ZHU C X, ZHU J Z, ZHU H J, HU Y, GAO B D. Characterization of a novel single-stranded RNA virus, closely related to fusariviruses, infecting the plant pathogenic fungus. Virus Research, 2016, 217: 1-7.

[21] 向均. 我国梨黑斑病病原菌种类及一种dsRNA病毒的鉴定[D]. 武汉: 华中农业大学, 2017.

XIANG J. Identification of the pathogenic fungi of pear black spot disease from china and a dsRNA virus in them[D]. Wuhan: Huazhong Agricultural University, 2017. (in Chinese)

[22] OKADA R, ICHINOSE S, TAKESHITA K, URAYAMA S I, FUKUHARA T, KOMATSU K, ARIE T, ISHIHARA A, EGUSA M, KODAMA M, MORIYAMA H. Molecular characterization of a novel mycovirus inmanifesting two-sided effects: Down-regulation of host growth and up-regulation of host plant pathogenicity. Virology, 2018, 519: 23-32.

[23] 李朋华. 部分链格孢属真菌的形态学及多基因鉴定[D]. 郑州: 河南农业大学, 2016.

LI P H. Identification of partialby morphology and multiple-gene analysis[D]. Zhengzhou: Henan Agricultural University, 2016. (in Chinese)

[24] 李波. 核盘菌弱毒菌株SX276生物学特性及其携带RNA病毒的研究[D]. 武汉: 华中农业大学, 2016.

LI B. Research on biological characteristics and RNA viruses from the hypovirulent strain SX276 of[D]. Wuhan: Huazhong Agricultural University, 2016. (in Chinese)

[25] CHIBA S, SALAIPETH L, LIN Y H, SASAKI A, KANEMASTU S, SUZUKI N. A novel bipartite double-stranded RNA mycovirus from the white root rot fungus: Molecular and biological characterization, taxonomic considerations, and potential for biological control. Journal of Virology, 2009, 83(24): 12801-12812.

[26] YU X, LI B, FU Y P, XIE J T, CHENG J S, GHABRIAL S A, LI G, YI X H, JIANG D H. Extracellular transmission of a DNA mycovirus and its use as a natural fungicide. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(4): 1452-1457.

[27] MELZER M S, DUNN M, ZHOU T, BOLAND G J. Assessment of hypovirulent isolates offor potential in biological control of chestnut blight. Canadian Journal of Plant Pathology, 1997, 19(1): 69-77.

[28] 赵迎彤.病毒dsRNA的多样性研究[D]. 泰安: 山东农业大学, 2014.

ZHAO Y T. Diversity of dsRNA fromfungal viruses[D]. Taian: Shandong Agricultural University, 2014. (in Chinese)

[29] 王利华, 罗慧, 王国平, 王利平. 来源于产黄青霉病毒科成员的分离物对梨轮纹病菌生长及其致病力的影响. 果树学报, 2017, 34(10): 1330-1339.

WANG L H, LUO H, WANG G P, WANG L P. Effect ofchrysovirus 1 isolate belonging to thefamily on growth and pathogenicity of thestrain infection in pears. Journal of Fruit Science, 2017, 34(10): 1330-1339. (in Chinese)

[30] MA G P, ZHANG X F, HUA H H, ZHOU T, WU X H. Molecular and biological characterization of a novel strain ofchrysovirus 1 identified from the pathogencausing watermelon leaf blight. Virus Research, 2020, 280: 197904.

Analysis of dsRNA carried byf. sp.in China and identification of a dsRNA virus

Cao YuHan, Li ZiTeng, Zhang JingYi, Zhang JingNa, Hu TongLe, Wang ShuTong, Wang YaNan, Cao KeQiang

College of Plant Protection, Hebei Agricultural University, Baoding 071001, Hebei

【Objective】Apple Alternaria blotch occurs in all major apple producing areas in China, causing serious economic loss. Apple Alternaria blotch is a worldwide airborne disease caused byf. sp. The objective of this study is to investigate the diversity of dsRNA off. spin China, and to provide new biological control resources for the disease and new understanding of virus diversity and evolutionary.【Method】Tissue samples with typical symptom were collected from eight provinces in China. Pure culture was obtained by tissue isolation and single spore separation. The diversity of dsRNA carried byf. spin China was determined by dsRNA extraction and gel electrophoresis. The biological characteristics of dsRNA-carrying pathogens were determined by evaluation of culture characteristics, mycelial growth rate, and pathogenicity on fruits and leaves. By high throughput sequencing and molecular cloning techniques, the whole genome sequence, genome structure and phylogenetic for the mycovirus carried by QY-2 strain were analyzed.【Result】The pure culture of 102 strains off. spwas obtained and five strains with obvious dsRNA bands were found. The dsRNA virus off. spin China could be divided into four types, there were about 8, 2.5 and 1.5 kb dsRNA in type I (strain yt-3-7), about 8 kb dsRNA in type II (strain sjz-4), about 3 and 0.8 kb dsRNA in type III (strain qy-2), and about 2.5 and 1.5 kb dsRNA in type IV (strains cl-2-6 and sq-1-1). The strains carrying dsRNA had diverse culture characteristics, and there was no significant relationship between the colony characteristics, growth rate and whether dsRNA carrying or not and its type. The strains carrying dsRNA were mostly of low pathogenicity type. It was found that QY-2 strain showed weak pathogenicity on leaves and fruits, so QY-2 was selected for virus identification. The genome of the dsRNA virus carried by the strain QY-2 consists of five dsRNA fragments, dsRNA1-dsRNA5, with the sequence sizes of 3 665, 3 054, 2 824, 2 819 and 831 nt, which were submitted to GenBank. The accession numbers are MK672910, MK672913, MK672912, MK672911, and MK836314. The molecular weights of coding proteins are 124, 83, 84, 83 and 13 kD, respectively. Phylogenetic analysis showed that it was closely related to AaCV1, so it was identified aschrysovirus belonging toe ands. The tentative designation ischrysovirus 2 (AaCV2). 【Conclusion】There are diverse dsRNAs inf. spin China, and there is no significant correlation between dsRNA diversity and host pathogen culture traits. Most of the strains carrying dsRNA showed low pathogenicity. QY-2 strain with the least pathogenicity carries AaCV2. The low pathogenicity characteristic of this strain may have potential application value, which provides a new resource for biological control of apple Alternaria blotch.

apple Alternaria blotch;f. sp; dsRNA; biological characteristic; pathogenicity; biological control

2021-04-06；

2021-05-10

国家现代农业产业技术体系建设专项（CARS-27）、河北省自然科学基金（C2019204327）、河北省重点研发计划（21326506D，19226508D，20327405D）、河北省引进留学人员资助项目（C201839）

曹钰晗，E-mail：c897836047@163.com。李紫腾，E-mail：lzt304838566@163.com。曹钰晗和李紫腾为同等贡献作者。通信作者王亚南，E-mail：wyn3215347@163.com。通信作者曹克强，E-mail：cao_keqiang@163.com

（责任编辑 岳梅）