0 引言
【研究意义】衔接蛋白（adaptor protein,AP）在哺乳动物细胞内信号的传递、蛋白质的分拣及向囊泡转运过程中发挥着重要作用[1],但关于AP蛋白在真菌中的功能鲜有报道。以GcAP1复合体β亚基基因为出发点,揭示其在苹果炭疽叶枯病菌生长发育和致病过程中的功能,可为植物病害防治的新途径、新方法提供线索,并对植物病原真菌致病机理研究具有重要意义。【前人研究进展】近几年,苹果一种重要的新病害——苹果炭疽叶枯病（Glomerella leaf spot of apple）,在中国苹果主产区连年大发生,给苹果产业带来巨额损失[2]。现已发现,围小丛壳菌（Glomerella cingulata）（无性态：胶胞炭疽菌Colletotrichum gloeosporioides）[3-5]、尖孢炭疽菌（C. acutatum）[6]、喀斯特炭疽菌（C. karstii）[7]、果生刺盘孢（C. fructicola）和隐秘刺盘孢（C. aenigma）[8]均可引起苹果炭疽叶枯病。炭疽菌通过分生孢子萌发,产生芽管和附着胞,附着胞下形成的侵染钉可直接侵入寄主[4]。像其他植物病原真菌一样,炭疽菌（Colletotrichum）侵染过程受cAMP和MAPK路径调控[9]。胶胞炭疽菌通过分泌果胶酶,例如聚半乳糖醛酸酶（polygalacturonase,PG）、果胶裂解酶（pectin lyase,PNL）和果胶酸酯裂解酶（pectate lyase,PEL）,克服寄主细胞壁中的果胶,使其成功在寄主体内定殖。研究表明,PG1、PG2、pnl-1、pnl-2、pelA和pelB单个基因突变或双基因突变均明显地降低病原菌的致病力[10-14]。最近的研究表明,在碱性条件下,转录因子pacC调控pelB的表达[15],但关于果胶酶基因的调控机制尚需进一步研究。AP蛋白以复合体形式存在于生物体中,具有重要的功能[16-17]。在哺乳动物中,现已发现5种类型的AP复合体,即AP-1、AP-2、AP-3、AP-4和AP-5[18-19]。AP蛋白由4个不同的亚基构成（δ、β、μ和σ）,2个大亚基（δ和β）具有识别和结合网格蛋白的功能,1个中亚基（μ）负责蛋白的分拣和跨膜装配,小亚基（σ）具有稳定复合体的作用[18,20-21]。目前,胶胞炭疽菌菌株23的基因组测序已完成,其Contig序列公布在“Comparative Fungal Genomics Platform”。此外,博德研究所（http://www.broadinstitute.org）还公布了玉米炭疽菌（C. graminicola）和希金斯刺盘孢（C. higginsiamim）全基因序列。这些基因组信息为炭疽菌的致病机理研究创造了条件。【本研究切入点】为开展苹果炭疽叶枯病菌致病机制研究,笔者课题组之前利用农杆菌介导的遗传转化（ATMT）方法,构建了苹果炭疽叶枯病菌强致病力菌株W16的T-DNA插入突变体库[5]。在对一些T-DNA插入突变体进行致病性测验中,发现一株突变体A7346对苹果叶片及果实致病力明显降低。进一步揭示突变体A7346致病力降低的分子机制,将是研究苹果炭疽叶枯病菌致病机制的一个很好契机。【拟解决的关键问题】阐明GcAP1复合体β亚基在苹果炭疽叶枯病菌致病过程中的功能,明确GcAP1复合体是否调控CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB的表达,为深入开展AP蛋白在苹果炭疽叶枯病菌致病信号传导途径中的分子机制研究打下基础。1 材料与方法
试验于2015—2016年在中国农业科学院果树研究所完成。1.1 供试菌株及培养条件
苹果炭疽叶枯病菌（G. cingulata）强致病力菌株W16,由笔者课题组2014年分离自辽宁省绥中市一个病害重发生果园[5],W16和由其衍生的转化子采用PDA培养基进行培养。农杆菌（Agrobacterium tumefacien）菌株LBA4404（由云南农业大学何月秋教授惠赠）采用LB培养基进行培养,用于ATMT。大肠杆菌（Escherichia coli）菌株TG1采用LB培养基进行培养,用作质粒的宿主。为检测GcAP1β的表达模式,RNA提取样品分别取自菌丝（PDA培养3 d的菌丝体）、分生孢子（无菌蒸馏水刮洗PDA平板10 d的培养物,2层擦镜纸过滤除去菌丝后,离心收集分生孢子）、芽管（用毛笔将106个/mL分生孢子悬浮液涂在叶片正面,28℃条件下保湿培养6 h后,用灭菌刀片轻刮叶子叶面,收集萌发的分生孢子）、附着胞（方法同芽管,保湿培养10 h,收集叶片表面形成的附着胞）和叶片接种48 h（方法同芽管,保湿培养48 h,将整个叶片液氮速冻后研磨）。为检测CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB在GcAP1β突变体中的表达量,RNA提取样品分别取自PDA培养6 d的菌丝体。1.2 核酸操作
基因组DNA的提取采用HE[22]的方法进行。PCR扩增体系和程序按ZHANG等[23]的方法进行。PCR产物的纯化采用生工生物工程（上海）股份有限公司（简称：生工生物）的PCR产物纯化试剂盒进行。质粒的小量提取采用生工生物的质粒提取试剂盒进行。Southern blot分析采用ZHANG等[24]的方法进行。总RNA的提取采用北京天根生化科技有限公司的植物总RNA提取试剂盒（HL缓冲液）进行。RT-PCR分析采用大连宝生物工程有限公司的反转录试剂盒进行。qRT-PCR分析采用北京全式金生物技术有限公司的荧光定量PCR试剂盒进行。qRT-PCR每个反应包含2xTransStart Tip Green qPCR SuperMix 10 μL,引物（10 μmol·L-1）各0.4 μL,cDNA（RT-PCR产物100倍稀释液）1 μL,RNase-free H2O 8.2 μL。qRT-PCR采用美国BIO-RAD公司的荧光定量PCR分析系统（Real-Time PCR System）进行检测,扩增程序为1个循环（94℃ 30 s）,40个循环（94℃ 5 s,60℃ 15 s,72℃ 10 s）,并以β-微观蛋白基因（β-tubulin）作为内参基因,用2-ΔΔCt方法[25]计算各被检测基因的相对表达量。PCR引物（表1）由生工生物合成。Table 1
表 1
表 1PCR引物序列
Table 1The primers used in this study
| 引物Primer | 序列Sequence (5′-3′) |
|---|---|
| S2G2 | AAAGAATTCGAAGGATCCACCCAGCTCTG# |
| S2Df2 | AAAGTCGACGCAATGGGGCAACAGGTAGTAGTAG |
| S2Dr1 | ATAGTCGACAAACATGGGCCAGAAGGAGAGG |
| S3S1 | AAAACTAGTGATGATGTAGGCGGAAATTAGC |
| HB1 | AAAGTCGACCAGCACGAGCAGGTCCCC |
| HB2 | TTTGGATCCTTGACCTTCAACGTGTCCATC |
| S2F | TCACAGCACGAGCAGGTCC |
| S2R | ATGTTCATGGCCTATGTGCCC |
| S2R1 | ATCCTCTGCCGACTTTCGTAC |
| S2Q1 | CGCATCCTGGCATAGTTGAC |
| S2Q2 | CGTTGTTCCCCGACATCATCG |
| TubF | CTTCCGGCAACAAGTACGT |
| TubR | GCGTCCTGGTATTGCTGGT |
| H-F1 | CGTTATGTTTATCGGCACT |
| H-R1 | TTGGCGACCTCGTATTGG |
| P-f | TACGAAAGTCGGCAGAGGAT |
| P-r | CTTAGCAGCAACTGGAGGTCA |
| CgPG1F | TACCACCCTCGACATGACCG |
| CgPG1R | CACTCCTTGTAGCCGAAGGT |
| CgPG2F | TCCTGCACCGACATCGTG |
| CgPG2R | GGTGCCGTCTTTCAGGCC |
| Pnl-1F | AAGGGCGCTATCAAGGGCA |
| Pnl-1R | GGGTTCAAGTCGGAGACCTGAAT |
| Pnl-2F | GCATCACCAACTCAGCGAAG |
| Pnl-2R | ATCACCTCCCCAGACGTACT |
| PelAF | GGCACTGGTCACATCTACAACTC |
| PelAR | TGCTCTGAATGAGGACCTGC |
| PelBF | TGTCCGCAACATGAAGATCTCC |
| PelBR | AGTGGTCAACCCAGACCTTGGA |
| TubQF | TGATGGCTGCTTCTGACTTCC |
| TubQR | TCCTCGACATCCTTCATAGCG |
新窗口打开
1.3 突变体A7346 T-DNA右翼序列分析
以质粒pCamhybgfp1（GenBank accession number：KX223837）为模板,用引物H-F1和H-F2进行PCR扩增,产物纯化后作为探针（P1）进行A7346 T-DNA插入拷贝数的Southern blot分析。A7346 T-DNA右翼序列采用hiTAIL-PCR方法[26]进行PCR扩增,产物连结到pUCm-T载体（生工生物）,由北京三博远志生物技术有限责任公司测序。获得的T-DNA右翼序列参考围小从壳菌菌株23基因组（http://genome.jgi.doe. gov/Gloci1/Gloci1.home.html）,进行本地blast分析。目的GcAP1复合体β亚基基因（GcAP1β）利用比较真菌基因组平台[27]进一步进行分析,并提取基因结构及其上下游序列和蛋白序列信息。1.4 GcAP1复合体β亚基基因的敲除与恢复
以野生型菌株W16基因组为模板,用引物S2Dr1和S3S1扩增GcAP1β基因上游序列,SalⅠ/SpeⅠ酶切后,连结到SalⅠ/SpeⅠ酶切的pCambiaMX9（GenBank accession number：KX755248）,生成质粒pCamS1。再以野生型菌株W16基因组为模板,用引物S2G2和S2DF2扩增GcAP1β基因下游序列,EcoRⅠ/SalⅠ酶切后,连结到EcoRⅠ/SalⅠ酶切的pCambiaMX9（GenBank accession number：KX755248）,生成质粒pCamS2。约2.0 kb的潮霉素抗性基因（hph）序列从质粒pTFCM[28]（由华中农业大学姜道宏教授惠赠）SalⅠ酶切得到后,连结到SalⅠ酶切的pCamS2,形成GcAP1β基因敲除载体pCamS2KN1（图1）。pCamS2KN1被电转到农杆菌LBA4404中后,利用ATMT方法[29]将GcAP1β从W16基因组中敲除。筛选潮霉素B（生工生物,终浓度100 μg·mL-1）抗性转化子,对假定的GcAP1β基因突变体用引物HB1和S2R1进行特定位点的PCR扩增分析（图1）,并以探针P2（以野生型菌株W16基因组为模板,由引物P-f和P-r进行PCR扩增及纯化获得）进行Southern blot鉴定。以W16基因组为模板,用引物HB1和HB2扩增GcAP1β及其假定的启动子序列,BamHⅠ/SalⅠ酶切后,连结到载体pGapneoR12（GenBank accession number：KY363244）,生成的GcAP1β基因恢复载体pCamGcAP1β（形成GcAP1β-gfp融合表达结构）被电转到农杆菌LBA4404中,利用ATMT[29]将GcAP1β-gfp融合基因转移到Δgcap1β突变体基因组中,筛选G418（生工生物,终浓度500 μg·mL-1）抗性转化子,对假定的GcAP1β基因恢复菌株用引物S2F和S2R进行PCR扩增分析和GFP信号检测（DM5000B型荧光显微镜,德国Leica）进行确定。
显示原图|下载原图ZIP|生成PPT图 1GcAP1β的克隆和鉴定
-->Fig. 1The cloning and identification of GcAP1β
-->
1.5 表型分析
真菌的生长速度和产孢能力测验采用WU等[30]的方法进行。挑取气生菌丝接种于PDA平板上（直径9 cm）,培养6 d后,测量菌落直径;平板继续培养至10 d,用10 mL无菌水刮洗平板,用擦镜纸滤去菌丝后,采用血球计数板进行测定,计算每毫升孢子悬浮液中的产孢量。分生孢子的萌发、附着胞的形成能力测验采用XU等[31]的方法。将无菌的疏水载玻片至于水琼脂（8%）平板上,将10 μL分生孢子悬浮液（104个/mL）滴于载玻片,28℃黑暗条件下,保湿培养至8 h,观测生孢子萌发情况,保湿培养至20 h,观测附着胞形成情况。每个重复观测100个分生孢子。计算公式：孢子萌发率（%）=萌发的孢子数/检测的孢子总数×100;附着胞形成率（%）=形成附着胞的萌发孢子数/检测的萌发孢子总数×100。叶片致病性测定：摘取成熟度一致的金冠（感病品种）叶片,无菌水冲洗后,用2 μL分生孢子悬浮液（1×105个/mL）滴在叶片正面,每个叶片接种6个悬滴,接种后,置于塑料盒中冷凝的10%水琼脂上。28℃条件下,黑暗培养10 h后,转入光照14 h,以此循环。每天观察发病情况。苹果致病性测定：苹果果实被75%酒精消毒后,用牙签刺约3 mm深的伤口,2 μL分生孢子悬浮液（1×105个/mL）接种于伤口处,28℃条件下,接种10 d后观测果实发病情况。野生型菌株W16作为阳性对照。每个处理重复3次,用DPS软件进行差异显著性分析。
1.6 生物信息学分析
应用ClustalX对GcAP1β基因序列及其cDNA序列进行比对,以确定GcAP1β基因内含子的有无。应用InterProScan[32-33]对GcAP1β进行保守结构域（Conserved Domain）分析。应用SignalP 4.1[34]对GcAP1β的N端进行信号肽分析。应用ProtComp 9.0（http://www.softberry.com）对GcAP1β进行亚细胞定位分析。应用TMHMM Server v. 2.0（http://www.cbs. dtu.dk/services/TMHMM/）对GcAP1β进行跨膜螺旋结构分析。2 结果
2.1 GcAP1β的克隆与分析
Southern blot分析结果表明,T-DNA插入突变体A7346经过酶切的泳道只杂交出一个条带,推测其表型是由T-DNA以单拷贝插入至基因组所致（图2）。运用hiTail-PCR技术克隆得到的T-DNA右翼序列与围小从壳菌菌株23基因组进行比对,发现T-DNA插入位点位于一个注释但功能尚未解析的基因（e_gw1.11.622.1）的第一个外显子内。InterProScan分析显示该基因蛋白序列含有典型的衔接蛋白AP1复合体β亚基结构域（IPR026739）。因此,将这一基因命名为GcAP1β（GcAP1复合体β亚基基因）。GcAP1β基因全长2 321 bp,含有3个内含子,编码720个氨基酸。SignalP预测该基因蛋白序列无信号肽,表明其可能不是分泌蛋白。
显示原图|下载原图ZIP|生成PPT图 2T-DNA插入突变体A7346 Southern杂交分析
-->Fig. 2Southern hybridization analysis of T-DNA insertional mutant A7346
-->
2.2 Δgcap1β突变体和GcAP1β基因恢复菌株的分子鉴定
特定位点的PCR检测结果（图3-A）显示,Δgcap1β突变体Δgcap1β-1、Δgcap1β-2和Δgcap1β-3基因组中的GcAP1β均被敲除。Southern blot分析结果显示,Δgcap1β-1、Δgcap1β-2和Δgcap1β-3均杂交出一条带,而野生型菌株W16杂交出与敲除突变体长度不同的一条带（图3-B）,表明外源的hph在Δgcap1β突变体基因组中为单位点、单拷贝的插入,也证实GcAP1β被hph替换。RT-PCR检测结果（图3-C）表明,GcAP1β被导入到突变体Δgcap1β-1基因组中,且能在恢复菌株Δgcap1β-GcAP1β中表达。
显示原图|下载原图ZIP|生成PPT图 3Δgcap1β突变体和GcAP1β基因恢复菌株的分子鉴定
-->Fig. 3Molecular confirmation of the Δgcap1β mutants and the GcAP1β complementation strain Δgcap1β-GcAP1β
-->
2.3 GcAP1β缺失对菌丝生长的影响
在PDA平板上,与野生型菌株W16相比,Δgcap1β突变体Δgcap1β-1、Δgcap1β-2和Δgcap1β-3菌落成褶皱状,菌丝生长速度明显减慢（图4、表2）,但分生孢子产量、分生孢子萌发率、附着胞形成率则无显著差异（表2）。
显示原图|下载原图ZIP|生成PPT图4苹果炭疽叶枯病菌野生型菌株及其衍生的转化子菌落形态
-->Fig. 4Colony morphology of G. cingulata strain W16 and its derived transformants
-->
Table 2
表2
表2苹果炭疽叶枯病菌野生型菌株及其衍生的转化子表型分析
Table 2Phenotypic analysis of G. cingulata strain W16 and its derived transformants
| 菌株 Strain | 菌落直径 Colony diameter (cm) | 产孢量 Conidiation (×106) | 萌发率 Germination rate (%) | 附着胞形成率 Appressoria formation rate (%) |
|---|---|---|---|---|
| W16 | 8.6±0.13 | 5.3±1.2 | 96.5±0.5 | 96.3±0.5 |
| Δgcap1β-1 | 3.1±0.09** | 3.5±2.5 | 95.3±0.4 | 94.9±0.9 |
| Δgcap1β-2 | 3.3± 0.12** | 4.7±1.2 | 94.8±0.6 | 94.1±0.6 |
| Δgcap1β-3 | 3.3± 0.11** | 4.4±1.3 | 95.5±0.9 | 94.7±0.3 |
| Δgcap1β-GcAP1β | 8.4±0.15 | 3.6±1.8 | 95.3±0.8 | 95.1±0.6 |
新窗口打开
2.4 GcAP1β的缺失对致病力的影响
致病测验结果显示,Δgcap1β在苹果叶片上引起极小的点状斑,与A7346极为相似,然而同条件下,W16引起明显的坏死斑（图5-A）;Δgcap1β在苹果果实上引起较小的凹陷斑,而W16引起的病斑明显比Δgcap1β大的多（图5-B）。此外,GcAP1β基因恢复菌株Δgcap1β-GcAP1β完全修复了因GcAP1β基因缺失造成的表型缺陷（图5、表2）。这些结果表明GcAP1β的缺失导致苹果炭疽叶枯病菌菌丝生长减慢和致病力降低。
显示原图|下载原图ZIP|生成PPT图5苹果炭疽叶枯病菌野生型菌株及其衍生的转化子致病力测定
-->Fig. 5Pathogenicity assay of G. cingulata strain W16 and its derived transformants
-->
2.5 GcAP1β亚细胞定位
TMHMM Server v. 2.0分析结果表明,GcAP1β不具有跨膜螺旋结构;ProtComp 9.0分析结果表明,GcAP1β没有明确的定位信号。为了确定该蛋白的亚细胞定位,构建了GcAP1β-gfp融合表达的GcAP1β恢复菌株Δgcap1β-GcAP1β。荧光显微镜下检测结果表明,融合蛋白GcAP1β-GFP分布于细胞质中,与T-DNA插入突变体A7346中GFP蛋白（gfp的表达由GAPDH启动子控制）分布特点相似（图6）。
显示原图|下载原图ZIP|生成PPT图6融合蛋白GcAP1β-GFP分布于细胞质中
-->Fig. 6The fused protein GcAP1β-GFP is distributed to cytoplasm
-->
2.6 GcAP1β在苹果炭疽叶枯病菌发育各阶段表达情况
qRT-PCR检测结果表明,GcAP1β在苹果炭疽叶枯病菌各个发育阶段均有表达。相对于菌丝阶段（PDA培养3 d时）的表达量,GcAP1β在分生孢子、芽管和附着胞阶段表达量均有不同程度的降低,而在侵染后（叶片接种48 h）表达量却明显提高（图7）。
显示原图|下载原图ZIP|生成PPT图 7GcAP1β在各发育阶段的表达情况
-->Fig. 7The expression of GcAP1β at different development stages
-->
2.7 GcAP1β的缺失对CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB表达量的影响
由于Δgcap1β突变体导致苹果炭疽叶枯病菌致病力的降低（图7）,笔者假设Δgcap1β突变体中果胶酶相关基因的转录受阻。因此,利用qRT-PCR技术检测了CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB在Δgcap1β突变体中的表达情况。相比于野生型菌株W16,在Δgcap1β突变体中,CgPG1表达量降低至20.3%,CgPG2表达量降低至16.5%,pnl-1表达量降低至8.2%,pnl-2表达量降低至14.4%,pelA表达量降低至4.4%,pelB表达量降低了至0.8%（图8）。这些结果显示GcAP1β正调控基因CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB的表达。
显示原图|下载原图ZIP|生成PPT图 8Δgcap1β突变体CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB的表达量分析
-->Fig. 8The expression analyses of CgPG2, pnl-1, pnl-2, pelA and pelB in the Δgcap1β mutant
-->
3 讨论
真菌菌丝的生长与囊泡（vesicles）的营养运输密切相关[35]。本研究结果表明,GcAP1复合体β亚基的缺失导致苹果炭疽叶枯病菌菌丝生长缓慢;qRT-PCR检测结果表明,GcAP1β在菌丝生长阶段表达量相对较高（图7）。值得注意的是,在哺乳动物中,衔接蛋白AP-1复合体负责运输蛋白至网格蛋白包被的囊泡（clathrin-coated vesicles）中[18,36]。因此,笔者推测真菌菌丝的生长与GcAP1复合体负责运输蛋白至囊泡的功能密切相关。胶胞炭疽菌侵入寄主后,先保持活体营养型生长,然后转为死体营养型生长[30]。苹果炭疽叶枯病菌发病速度特别快,从侵染叶片到发病落叶仅需3—4 d[37]。本研究中,苹果炭疽叶枯病菌叶片接种48 h时已出现病斑。qRT-PCR检测结果显示,GcAP1复合体β亚基基因在苹果炭疽叶枯病菌各个发育阶段均有表达,但是以侵染后（叶片接种48 h时）表达量相对最高（图7）。这些结果说明,GcAP1复合体β亚基基因在苹果炭疽叶枯病菌死体营养型生长阶段中发挥着重要作用。在植物病原真菌中,由附着胞调节的侵入寄主方式在其致病过程起到重要作用。在炭疽菌中,丝裂原活化蛋白激酶（Mitogen-activated protein kinase,MAPK）级联信号途径调控其分生孢子的萌发、附着胞的形成[9]。在酿酒酵母（Saccharomyces cerevisiae）中,衔接蛋白Ste50与G蛋白有关的Cdc42-Ste20激酶复合体互作后,激活MAPKKK基因Ste11的转录,进而激活MAPK级联信号传递链[38]。在胶胞炭疽菌中,MAPK途径基因中的CgMEK1或Cgl-SLT2的缺失,导致萌发的分生孢子不能形成附着胞。本研究中,GcAP1复合体β亚基基因的缺失对分生孢子的萌发、附着孢的形成无影响（表2）。这一结果暗示,GcAP1复合体可能对基因CgMEK1和Cgl-SLT2的转录无直接或间接的影响。
在植物病原菌和寄主相互作用中,病原菌分泌果胶酶降解寄主细胞壁的果胶聚合体,促进病原菌的侵入和定殖[10-14]。此外,病原菌的果胶酶分泌量不仅决定其致病力的强弱,还决定其致害的症状[39]。因为Δgcap1β突变体仅在苹果叶片上引起极小的点状斑以及在苹果果实（有伤接种）上引起较小的凹陷斑（图5）,笔者假设Δgcap1β突变体果胶酶有关基因（CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB）转录可能受阻。qRT-PCR检测结果表明,Δgcap1β突变体果胶酶有关基因的表达量与野性型菌株的相比明显降低（图8）。这些结果表明,Δgcap1β突变体致病力下降的一个重要原因是,GcAP1蛋白β亚基基因的缺失阻遏果胶酶有关基因的表达。但是,GcAP1蛋白是否直接或间接调控果胶酶有关基因的转录,有待进一步研究。
4 结论
衔接蛋白GcAP1复合体分布于细胞质中。GcAP1β在分生孢子、芽管和附着胞阶段均表达,但表达量明显低于菌丝生长阶段,是苹果炭疽叶枯病菌生长发育所需要的。GcAP1调控CgPG1、CgPG2、pnl-1、pnl-2、pelA和pelB的表达,是苹果炭疽叶枯病菌一个重要的毒力因子。The authors have declared that no competing interests exist.
