基于FUNGuild的镰刀菌根腐病发病烟株根际真菌群落研究.docx

1、基于FUNGuild的镰刀菌根腐病发病烟株根际真菌群落研究基于FUNGuild的镰刀菌根腐病发病烟株根际真菌群落研究李朋发1, 2，杨龙3，李桂龙1, 2，徐后娟3，王玉军3，李忠佩1, 2【摘 要】为了解镰刀菌根腐病发病烟株与未发病烟株在根际真菌群落结构与功能上的差异，联合应用高通量测序和FUNGuild技术研究了镰刀菌根腐病发病烟株与未发病烟株根际土壤中的真菌群落，结果表明：（1）镰刀菌根腐病发病烟株根际真菌群落的Shannon指数为4.14，而未发病烟株则高达4.71，未发病烟株根际真菌群落多样性显著高于发病烟株，镰刀菌根腐病发病烟株与未发病烟株根际土壤中的真菌群落在组成上有显著差异；（

2、2）镰刀菌根腐病发病烟株根际真菌以病理营养型为主，占54.64%，而未发病烟株根际真菌则以腐生营养型为主，占46.10%；（3）镰刀菌根腐病发病烟株根际的病原菌相对丰度为30.73%，显著高于未发病烟株的1.12%，镰刀菌在病原菌群落中占绝对优势地位。高通量测序与FUNGuild联用可对土传病害感病烟株的根际优势病原真菌群落进行准确鉴定，可为烟草病害的防控提供科学参考。【期刊名称】中国烟草学报【年(卷),期】2019(025)002【总页数】6【关键词】烟草；镰刀菌根腐病；根际；真菌；功能；FUNGuild【文献来源】植物保护引用本文：李朋发，杨龙，李桂龙，等基于FUNGuild的镰刀菌根腐病

3、发病烟株根际真菌群落研究J中国烟草学报，2019，25（2）基金项目：国家自然科学基金项目（41771298）、国家重点基础研究发展计划（973）项目（2015CB150501）网络出版日期：2019-01-31烟草土传病害严重危害我国的烟叶生产1。目前，对烟草土传病害的研究主要集中在三个方面，一是烟草病原菌的致病机理2-4；二是植烟土壤微生物群落组成和多样性5-6；三是烟草病害的综合防治措施7-9。而对于烟草土壤微生物群落功能的研究则鲜有报道。弄清烟草土壤微生物群落功能有助于深入了解烟草发病机理。烟草根际土壤微生物种类繁多，其中真菌是重要的组成部分。受限于技术手段，对烟草根际真菌群落的研究也

4、基本上都集中在群落组成和多样性上10-11。FUNGuild12是一个进行真菌功能比对的数据库，不仅可对真菌的营养类型进行鉴定，还能对真菌进行具体的功能分类。FUNGuild在真菌群落研究中得到了较为广泛的应用13-14。本研究联合使用Illumina高通量测序技术和FUNGuild软件，研究了镰刀菌根腐病烟株与未发病烟株根际土壤真菌群落的组成和功能，旨在深入了解镰刀菌根腐病发生的土壤微生物学机理，为镰刀菌根腐病的控制提供科学参考。1 材料与方法1.1 样品采集采样地位于山东省临沂市沂水县（3597N，11837E），该地区年均气温为13.0，年均降水量875mm。样地为坡地，坡顶为林地，坡中

5、和坡底为多年撂荒的新垦烟地，土壤类型为僵心棕堰土。坡中烟地中的烟株大规模发病，经鉴定为镰刀菌根腐病，平均发病率在95%以上。坡底烟地中的未见镰刀菌根腐病发病烟株。坡中土壤基本理化性质为：pH 6.41，有机碳8.42 gkg-1，总氮 0.89 gkg-1，总磷0.96 gkg-1，总钾17.21 gkg-1，碱解氮 86.65 mgkg-1，速效磷83.38 mg.kg-1，速效钾 103.17 mgkg-1；坡底土壤基本理化性质为：pH 6.38，有机碳8.46 g.kg-1，总氮 0.92 gkg-1，总磷0.97 gkg-1，总钾17.34 gkg-1，碱解氮 84.75 mgkg-1

6、，速效磷85.64 mgkg-1，速效钾 105.83 mg.kg-1。坡中与坡底土壤在理化性质上无显著差异。在坡中、坡底各分别随机选择3块烟地，每个地块分别随机选取5株烟株，取烟株根际土并将同一地块采集的根际土混合。将采集好的土壤样本加冰袋低温运输至实验室。将土样分为两部分，其中一部分立即放入-20冰箱保存，用于DNA提取；另一部分过2 mm筛并风干，研磨后过0.149 mm筛，用于土壤理化性质分析。1.2 土壤DNA提取及PCR扩增土壤DNA使用MP FastDNA试剂盒提取（MP Biomedicals，CA，USA），提取方法参考试剂盒使用说明书。提取的DNA使用PowerClean

7、DNA Clean-up Kit（MoBio，CA，USA）进 行 纯化。DNA浓度和纯度使用NanoDrop ND-1000分光光度计检测，结果显示，经过纯化的DNA浓度 均 在30 ngL-1以 上，OD260/OD280均 在1.8-2.0之间，浓度和纯度均符合ITS扩增要求。使用内转录间隔区1（ITS1）的特异性引物（ITS1-F：CTTGGTCATTTAGAGGAAGTAA15/ ITS2-R：GCTGCGTTCTTCATCGATGC16）对6个DNA样本分别扩增，PCR体系和条件参考Li的方法17。1.3 Illumina测序及序列处理使用Illumina MiSeq PE250测

8、序平台对扩增产物进行双端测序。测序数据经过Flash18和Trimmomatic19进行拼接和质控后，使用微生物生态学定量研究平台（Quantitative Insights into Microbial Ecology ，QIIME）进行处理。使用UCLUST20软件对各序列在97%的相似性水平下聚类操作分类单位（operational taxonomic units，OTUs），各个OTU的代表序列使用UNITE数据库（https:/unite.ut.ee/）进行物种信息比对21。1.4 数据处理与分析真菌的营养型、功能分组使用FUNGuild（1.0）软件进行比对，大体步骤为：（1）将Q

9、IIME输出的OTU表和物种信息表合并，将物种信息附加在OTU表的最后一列并命名为“taxonomy”；（2）将第（1）步得到的表格上传到FUNGuild网站，开始比对；（3）比对完成后下载输出结果到本地。输出结果包含原始的OTU表并依据序列数量进行排序。输出结果还包含物种的分类单位“Taxon”，物种分类级别“Taxon Level”，营养型“Trophic Mode”，功能分组“Guild”，可信度“Confidence”，生长形态“Growth Morphology”，特征“Trait”，备注“Notes”和比对信息来源“Citation/Source”。OTU数、OTU相对丰度等使用R

10、语言（3.4.3）进行聚类计算，主坐标分析（Principal coordinate analysis ，PCoA）、ANOSIM检验也通过R语言（3.4.3）进行。样本的描述性分析、组间差异的显著性分析使用SPSS 19.0中的独立t检验完成。所有图片均通过Origin（9.1）制作。显著性检验水准均为P0.05。2 结果与分析2.1 根际真菌群落结构及多样性比较经过质控，6个样本共得到184,782条高质量序列。在97%的相似性水平下聚类OTU，每个样本按照样本最少序列数抽平到14,644条序列，6个样本共聚类得到600个OTU，其中有444条代表序列通过UNITE数据库比对可鉴定为真菌。

11、在纲水平下对各样本的物种进行分类求和，得到各样本在纲水平下的物种组成，如图1所示。通过图1可以看出，未发病（H）与镰刀菌根腐病（S）烟株的根际真菌群落在物种组成上明显不同。未发病烟株根际真菌群落中，散囊菌纲（Eurotiomycetes）为优势菌纲，平均占比为42.81%，粪壳菌纲（Sordariomycetes）次之，平均占比为13.97%。发病烟株根际真菌群落中则正好相反，粪壳菌纲（Sordariomycetes）平均占比为55.06%，粪壳菌纲（Sordariomycetes）平均占比为14.61%。发病烟株根际真菌群落中的伞菌纲（Agaricomycetes）和鸟纲菌纲（Orbilio

12、mycetes）比例均显著高于未发病烟株（P0.05），但是座囊菌纲（Dothideomycetes）的相对丰度在两组真菌群落中无显著差异（P0.05）。主坐标分析（PCoA）结果显示，未发病和发病烟株根际真菌群落显著分异。进一步通过ANOSIM检验发现，两组烟株根际真菌群落的组间差异显著大于组内差异（P0.05）。分别选取Shannon-Wiener指数、Chao1指数和Observed OTUs指数描述两组烟株根际真菌群落的多样性，结果如表1所示，未发病组的三种多样性指数（Alpha diversities）均显著高于镰刀菌根腐病组（P0.05），说明未发病烟株根际真菌群落的多样性显著高于

13、镰刀菌根腐病发病烟株。2.2 根际真菌群落营养型比较根际FUNGuild数据库的比对结果，对两组烟株根际真菌的营养型（trophic mode）进行分类统计，结果如图3所示，未发病烟株根际真菌的营养型以腐生营养型（saprotroph）为主，平均占比46.10%，病理营养型（pathotroph）次之，平均占比22.85%；镰刀菌根腐病烟株根际真菌的营养型则以病理营养型为主，平均占比54.64%，腐生营养型次之，平均占比12.88%。两组烟株根际的病理营养型真菌和腐生营养型真菌的相对丰度均具有显著性差异（P0.05）。两组烟株根际的共生营养型真菌相对丰度均不足1%，但是未发病烟株根际的共生营养

14、型真菌相对丰度显著低于镰刀菌根腐病烟株（P0.05）。另外，未发病烟株根际有相对更多的真菌营养型尚未鉴别。2.3 根际真菌病原菌群落比较在600个OTU中，共有21个OTU经FUNGuild鉴定为植物病原菌（Plant pathogen），其中有6个OTU的鉴定结果可信度（Confidence Ranking）为“Possible”，因此仅对剩余15个鉴定结果可信度为“Probable”和“Highly probable”的高可信病原菌OTU进行分析。15个高可信度的病原菌OTU分别来自12个菌属，其中，镰刀菌属（Fusarium）、赤霉属（Gibberella）和青霉菌属（Penicilli

15、um）均有2个OTU被鉴定为高可信病原菌且镰刀菌属的两个OTU均为腐皮镰孢菌（F.solina），其它9个菌属均仅有1个OTU被鉴定为高可信病原菌。在所有高可信病原真菌中，镰刀菌占据绝对优势地位。在未发病组中，2个镰刀菌OTU的序列数占病原菌序列总数的75.02%10.67%，而在镰刀菌根腐病组中，2个镰刀菌OTU的序列数占病原菌序列总数的99.10%0.80%。镰刀菌根腐病烟株根际的镰刀菌在序列数和相对丰度上均显著高于未发病组（P0.05）。对两组烟株根际高可信病原菌OTU进行统计，结果如图4所示。未发病烟株根际病原真菌相对丰度平均为1.12%，而镰刀菌根腐病烟株根际病原真菌相对丰度则高达3

16、0.73%，差异显著（P0.05）。3 讨论根际微生物区系被称为植物的第二基因组，对植物的生长发育有决定性的影响22。本研究分析了镰刀菌根腐病发病烟株与未发病烟株根际土壤的ITS测序结果，发现两种根际土壤的真菌群落组成有显著差异，未发病烟株根际真菌多样性显著高于发病烟株，这与李小龙等对青枯病烟株根际真菌群落的研究结果相吻合23。镰刀菌根腐病发病烟株根际真菌群落多样性的降低，很可能是镰刀菌在群落的种间竞争中占据优势，少数病原菌不断积累并逐渐成为优势种，在这个过程中，相应的部分有益菌因为竞争能力弱而被淘汰24，物种数减少，多样性降低。除群落组成和多样性外，镰刀菌根腐病发病烟株与未发病烟株根际真菌的

17、营养型组成也存在显著差异。在未发病烟株根际，真菌以腐生营养型为主，病理营养型次之，而发病烟株根际真菌的营养型则以病理营养型为主，腐生营养型次之。病理营养型表示真菌通过伤害寄主细胞来获取自身所需的营养物质，而腐生营养型则表示真菌通过分解死亡的寄主细胞来获取营养物质25。镰刀菌根腐病发病烟株根际的病理营养型真菌显著高于未发病植株，说明发病植株根际可能有更多的致病真菌，尤其是镰刀菌，通过破坏根系细胞来维持自身生长，导致烟株根系组织遭到破坏，进而使病原菌能更加顺利地从根系对烟株进行侵染，最终导致烟株发病。分析结果显示，发病烟株根际的病原菌相对丰度也显著高于未发病烟株，2个镰刀菌在发病烟株根际的病原菌中

18、占据绝对的优势地位。镰刀菌是烟草致病菌，可引发烟草的镰刀菌根腐病26。其中，尖孢镰刀菌（F.oxysporum）、茄病镰刀菌（F.solani）等是侵染烟草根系并引发根腐病的主要病原菌种26-27，主要通过伤害烟草的维管束组织和根系组织等来侵染烟草27。因此，联合使用高通量测序与FUNGuild软件，可以对土传病害感病烟株的根际优势病原真菌群落进行准确的鉴定，同时也可提供丰富的、非优势的潜在病原真菌信息，这对揭示烟草土传病害发生的土壤微生物学机理具有重要意义，也可为烟草土传病害的防控提供科学参考。4 结论通过高通量测序与FUNGuild的联合使用，对镰刀菌根腐病发病烟株与未发病烟株根际土壤中的

19、真菌群落进行了研究，结果表明：（1）未发病烟株与发病烟株根际的真菌群落在组成上具有显著差异，未发病烟株根际真菌多样性显著高于发病烟株；（2）发病烟株根际真菌的营养型以病理营养型为主，腐生营养型次之，而未发病烟株则正好相反；（3）发病烟株根际的病原菌相对丰度显著高于未发病烟株，镰刀菌在病原菌群落中占绝对优势地位。高通量测序与FUNGuild的联用，可以对土传病害感病烟株的根际优势病原真菌群落进行准确的鉴定，也可为烟草土传病害的防控提供科学参考。参考文献1陈瑞泰, 朱贤朝, 王智发,等.全国16个主产烟省(区)烟草侵染性病害调研报告J.中国烟草科学, 1997, 18(4): 1-7.CHEN R

20、uitai, ZHU Xianchao, WANG Zhifa, et al.A report of investigating and studying tobacco infectious diseases of 16 main tobacco producing provinces (regions)in China J.Chinese Tobacco Science, 1997, 18(4): 1-7.2何永宏, 曾乙心, 刘林,等.温度和品种抗性对烟草青枯病潜育期的影响J.烟草科技, 2017, 50(6): 16-32 HE Yonghong, ZENG Yixin, LIU Li

21、n, et al.Effects of temperature and varietal resistance on latent period of tobacco bacterial wilt J.Tobacco Science & Technology, 2017, 50(6): 16-323潘建菁, 纪成灿, 刘冬霞,等.青枯菌侵染后烟株体内过氧化物酶活性的变化及其与抗病性的关系J.中国烟草科学, 2004, 25(3): 28-30.PANG Jianjing, JI Chengcan, LIU Dongxia, et al.Changes in activities of pero

22、xidase in flue-cured tobacco infected with Ralstonia solanacearum E.F.Smith and their relations to the resistance J.Chinese Tobacco Science, 2004, 25(3): 28-30.4高正良, 周本国, 雷艳丽.黑胫病菌在烟草品种间的侵染差异及I-S关系研究J.中国烟草学报, 2004, 10(4): 31-35.GAO Zhengliang, ZHOU Benguo, LEI Yanli.Study on the Infection difference

23、of tobacco variety to Phytophthora nicotianae and the relationship between incidence and severity J.Acta Tabacaria Sinica, 2004, 10(4): 31-35.5张笑宇, 段宏群, 芦阿虔,等.健康烟田与易感黑胫病烟田烟株不同生长时期根际土壤微生物区系变化规律J.河南农业科学, 2018, 47(3): 63-69.ZHANG Xiaoyu, DUAN Hongqun, LU Aqian, et al.Variation of Rhizosphere Soil Micro

24、bial Flora in Healthy and Black Shank-Susceptible Tobacco Fields at Different Growth Stages J.Journal of Henan Agricultural Sciences, 2018, 47(3): 63-69.6李远远.连作烟田土壤微生物多样性及微生物制剂应用研究D.郑州大学, 2017.LI Yuanyuan.Microbal Diversity in Continuous Cropped Tobacco Field and Application of Microbial Agents D.Zh

25、engzhou University, 2017.7万川, 蒋珍茂, 赵秀兰,等.深翻和施用土壤改良剂对烟草青枯病发生的影响J.烟草科技, 2015(2): 11-26.WAN Chuan, JIANG Zhenmao, ZHAO Xiulan, et al.Effects of Deep-ploughing and Soil Amendment Application on Incidence of obacco Bacterial Wilt J.Tobacco Science & Technology, 2015(2): 11-26.8夏艳, 徐茜, 董瑜,等.烟草青枯病菌拮抗菌的筛选、鉴

26、定及生防特性研究J.中国生态农业学报, 2014, 22(2): 201-207.XIA Yan, XU Qian, DONG Yu, et al.Screening, identification and characterization of antagonistic bacteria against Ralstonia solanacearum J.Chinese Journal of Eco-Agriculture, 2014, 22(2): 201-207.9王晶晶, 蒋士君, 常淑娴,等.两株生防菌对烟草黑胫病的抑制活性及其鉴定J.中国烟草学报, 2011, 17(6): 89-9

27、3.WANG Jingjing, JIANG Shijun, CHANG Shuxian, et al.Antifungal activity of two biocontrol strains against tobacco black shank and their identification J.Acta Tabacaria Sinica, 2011, 17(6): 89-93.10王婷婷, 王仁刚, 李咏梅,等.贵州烟草根围AM真菌多样性的初步研究J.菌物学报, 2014, 33(1): 143-151.WANG Tingting, WANG Rengang, LI Yongmei,

28、 et al.A primary study on diversity of arbuscular mycorrhizal fungi associated with tobacco roots from Guizhou Province J.Mycosystema, 2014, 33(1): 143-151.11王娜, 吕国忠, 孙晓东,等.烟草根际土壤真菌多样性的研究J.菌物学报, 2012, 31(6): 827-836.WANG Na, LV Guozhong, SUN Xiaodong, et al.Fungal diversity in tobacco rhizosphere so

29、il J.Mycosystema, 2012, 31(6): 827-836.12FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild J.Fungal Ecology, 2016, 20(1): 241-248.13Leff J W, Lynch R C, Kane N C, et al.Plant domestication and the assembly of bacterial and fungal communities associated with

30、strains of the common sunflower, Helianthus annuus J.New Phytologist, 2016, 214(1).14Kolakov Z, Kohout P, Krger C, et al.Root-associated fungal communities along a primary succession on a mine spoil: Distinct ecological guilds assemble differently J.Soil Biology & Biochemistry, 2017, 113: 143-152.15

31、Gardes M, Bruns TD.ITS primers with enhanced specificity for basidiomycetes-application to the identification of mycorrhizae and rusts.Molecular Ecology 1993; 2: 113-118.16White T, Bruns T, Lee S, Taylor FJRM, White T, Lee SH, et al.Amplification and direct sequencing of fungal ribosomal RNA genes f

32、or phylogenetics.PCR protocols: a guide to methods and applications.1990.17Li X, Zhang Y, Ding C, et al.Declined soil suppressiveness to Fusarium oxysporum, by rhizosphere microflora of cotton in soil sickness J.Biology & Fertility of Soils, 2015, 51(8): 935-946.18Mago T, Salzberg S L.FLASH: fast length adju