李 铮，王金辉，丁丽丽，等．贝莱斯芽孢杆菌菌株ＮＺ－４生防潜能及基因组学分析［Ｊ］．江苏农业科学，２０２３，５１（２）：１１７－１２５．
ｄｏｉ：１０．１５８８９／ｊ．ｉｓｓｎ．１００２－１３０２．２０２３．０２．０１８
上海好玩的地方推荐贝莱斯芽孢杆菌菌株ＮＺ－４生防潜能及基因组学分析
李 铮１，王金辉１，丁丽丽２，张 岱１，田梦君３，杨志辉１，朱杰华１
（１．河北农业大学植物保护学院／河北省植物病虫害生物防治工程技术研究中心，河北保定０７１０００；
２．河北省农林科学院生物技术与食品科学研究所，河北石家庄０５００５７；３．秦皇岛市植保植检站，河北秦皇岛０６６０００）
本山前妻全部照片Ｚ－４发酵液在不同理化条件下对尖孢镰刀菌的抑菌稳定性；通过鉴定培养基检测菌株ＮＺ－４的抗病促生能力；结合第二代（ＩｌｌｕｍｉｎａＮｏｖａＳｅｑ）及第三代（ＰａｃＢｉｏ）测序数据
拼装出ＮＺ－４基因组完整序列，通过系统发育分析确定该菌株的生物学种类，并将基因组序列与各数据库进行比对分析，预测菌株ＮＺ－４产生的酶类和抗生素相关基因。结果表明，菌株ＮＺ－４无菌发酵滤液耐高温、耐酸碱、耐紫外线，对胰蛋白酶不敏感。该菌株具备能够产生几丁质酶、果胶酶、纤维素酶、蛋白酶、酪蛋白酶并合成铁载体等分解真菌细胞壁和促生能力。Ｎ
Ｚ－４基因组中含与之对应的几丁质酶、葡聚糖酶、纤维素酶、蛋白酶的基因，还含有编码表面活性素、泛革素、杆菌霉素等９种抗生素的基因簇。菌株ＮＺ－４所产生的抑菌活性物质稳定性高，且基因组内存在与之相应的抗生素基因簇和破坏真菌细胞壁的酶类基因。说明该菌株具有很好的生防潜力，适合进一步开发成生防菌剂。 关键词：贝莱斯芽孢杆菌；生物防治；基因组学；抗生素；尖孢镰刀菌
中图分类号：Ｓ１８２ 文献标志码：Ａ 文章编号：１００２－１３０２（２０２３）０２－０１１７－０８
收稿日期：２０２２－０３－３１
基金项目：河北省薯类产业技术体系创新团队专项（编号：ＨＢＣＴ２０１８０８０２０５）；现代农业产业技术体系建设专项（编号：ＣＡＲＳ－０９－Ｐ１８）。
作者简介：李 铮（１９９４—），男，河北邯郸人，硕士研究生，主要从事生物防治研究。Ｅ－ｍａｉｌ：４１５２２８３９３＠ｑｑ．ｃｏｍ。
通信作者：杨志辉，博士，教授，主要从事马铃薯病害研究，Ｅ－ｍａｉｌ：ｂｄｙｚｈ＠ｈｅｂａｕ．ｅｄｕ．ｃｎ；朱杰华，博士，教授，主要从事马铃薯病害研究，Ｅ－ｍａｉｌ：ｚｈｕｊｉｅｈｕａ３５６＠１２６．ｃｏｍ。
暴走大事件 豫章书院 近年来我国马铃薯种植面积逐年增加，但是连作致使马铃薯枯萎病的发生逐年加重，严重影响了马铃薯的产量，成为制约我国马铃薯产业发展的重要因素。该病害由尖孢镰刀菌（Ｆｕｓａｒｉｕｍｏｘｙｓｐｏｒｕｍ）所引起，已成为严重威胁马铃薯产业的
重要土传病害之一［１］。该病原菌抗逆性强，使用化
学杀菌剂防治效果较差，长期使用还会造成环境污
染和农药残留［２］
。而生物防治可作为绿防治手
段，能有效解决化学杀菌剂所带来的环境污染问题。
生物防治具有绿无污染的特点，在多种病害防治过程中得到广泛应用。芽孢杆菌抗逆性强，广泛存在于土壤、植物表面及内部，对人和动物无毒
害的同时还具有广谱抑菌活性［３］
。贝莱斯芽孢杆
菌（
Ｂａｃｉｌｌｕｓｖｅｌｅｚｅｎｓｉｓ）是芽孢杆菌的一种，具有抑菌谱广、抗逆性强、生防机制多样等特点［
４］
，已成功应用于多种植物病害的防治过程，且效果显著。如贝莱斯芽孢杆菌Ｅ６９对稻瘟病的田间防效可达
８５ ９７％［５］
；贝莱斯芽孢杆菌Ｂ－３６在盆栽试验中，
对由尖孢镰刀菌引起的莲根腐病的防效最高可达
７７ １％［６］。贝莱斯芽孢杆菌能够合成分泌多种次级代谢产物实现对病原菌的抑制作用［７－８］，主要包
括由核糖体合成的细菌素和水解酶类，由非核糖体肽合成酶（ｎｏｎｒｉｂｏｓｏｍａｌｐｅｐｔｉｄｅｓｙｎｔｈｅｔａｓｅ，ＮＲＰＳ）合成的脂肽类化合物及多肽类化合物，以及由聚酮合成酶（ｐｏｌｙｋｅｔｉｄｅｓｙｎｔｈａｓｅ，ＰＫＳ）合成的聚酮类化合
物［４，９］
。细菌素、脂肽类物质理化性质稳定、抗菌谱
祚广、不易产生耐药性，是生物防治研究的重点［１０］
。
此外，贝莱斯芽孢杆菌分泌的蛋白酶、几丁质酶和
β
－葡聚糖酶等对抑制病原真菌起重要作用［１１］
。本实验室前期从马铃薯植株根际土壤中分离
得到菌株ＮＺ－４，该菌株对马铃薯枯萎病菌具有良好的抑制作用，抑菌圈直径在２０．０ｍｍ以上。基于ｇｙｒＢ基因的系统发育分析以及形态学特征，将菌株
ＮＺ－４初步鉴定为贝莱斯芽孢杆菌［１２］，但基于单一
基因位点的分析方法难以有效区分包括贝莱斯芽
孢杆菌在内的枯草芽孢杆菌复合［
１３］
。随着测序技术的发展，通过基因组数据分析物种间的亲缘关
系更加行之有效［４］。前期研究表明，菌株ＮＺ－４的抑菌谱较广，对包括马铃薯枯萎病菌、马铃薯病菌（Ｃｏｌｌｅｔｏｔｒｉｃｈｕｍｃｏｃｃｏｄｅｓ）、马铃薯疮痂病菌（Ｓｔｒｅｐｔｏｍｙｃｅｓｓｃａｂｉｅｓ）等有强烈的抑菌作用，且对小麦根腐病菌（Ｃｏｃｈｌｉｏｂｏｌｕ
ｓｓａｔｉｖｕｓ）、梨黑斑病菌（ＡｌｔｅｒｎａｒｉａｋｉｋｕｃｈｉａｎａＴａｎａｋａ）、苹果斑点落叶病菌（Ａｌｔｅｒｎａｒｉａａｌｔｅｒｎａｒｉａｆ．ｓｐ．ｍａｌｉ）、马铃薯早疫病菌（Ａｌｔｅｒｎａｒｉａｓｏｌａｎｉ）等重要农业病害病原菌也有良好的拮抗活性［１２］。此外，本研究室研究表明，该菌株具有促生、增强植株抗性水平、根际定殖能力强等多种应用潜力［１４］。而该菌无菌发酵液（含有抑菌作用的次生代谢产物）对外界环境条件稳定性、抑菌促生相关酶类活性有无及抑菌物质的遗传基础仍不明确。
本研究测定了高温、强酸碱、紫外线和蛋白酶处理下贝莱斯芽孢杆菌ＮＺ－４发酵滤液对引起马铃薯枯萎病的尖孢镰刀菌的抑菌活性稳定性，通过平板试验测定了其产生抑菌活性物质的能力，结合第二代（ＩｌｌｕｍｉｎａＮｏｖａＳｅｑ）和第三代（ＰａｃＢｉｏ）测序技术，获得ＮＺ－４完整的基因组序列，并挖掘了其脂肽类、聚酮类化合物、几丁质酶等主要抑菌活性物质的遗传基础，为生防菌剂开发和机理研究奠定了基础。
张楷依１ 材料与方法
１．１ 试验材料
１．１．１ 供试菌株 贝莱斯芽孢杆菌（Ｂａｃｉｌｌｕｓｖｅｌｅｚｅｎｓｉｓ）菌株ＮＺ－４和尖孢镰刀菌（Ｆｕｓａｒｉｕｍｏｘｙｓｐｏｒｕｍ），均由河北农业大学植物保护学院马铃薯病
害研究实验室提供。
１．１．２ 培养基 ＬＢ培养基、ＰＤＡ培养基、蛋白酶鉴定培养基［１５］、酪蛋白酶鉴定培养基［１６］、几丁质酶鉴定培养基［１７］、纤维素酶鉴定培养基［１８］、果胶酶鉴定培养基［１９］、ＣＡＳ培养基［２０］。
１．１．３ 菌株ＮＺ－４发酵滤液的制备 挑取菌株ＮＺ－４单菌落接种于ＬＢ培养基，１８０ｒ／ｍｉｎ、３７℃振荡培养１２ｈ获得种子液。将种子液接种于ＬＢ培养基中，１８０ｒ／ｍｉｎ、３６℃振荡培养７２ｈ，获得发酵液。将发酵液１００００ｒ／ｍｉｎ离心２０ｍｉｎ，上清液通过０ ２２μｍ的微孔过滤器后得到发酵滤液。
１．２ 菌株ＮＺ－４发酵滤液的抑菌活性稳定性检测１．２．１ 菌株ＮＺ－４发酵滤液处理 （１）酸碱稳定性：将发酵液的ｐＨ值分别调为３．０、４．０、５．０、６．０、７ ０、８．０、９．０、１０．０、１１．０和１２．０，于４℃静置１２ｈ后调节ｐＨ值＝７．０。（２）紫外光稳定性：将菌株ＮＺ－４发酵滤液置于距离紫外灯３０ｃｍ处，照射１、２、４、８、１２、１６、２０ｈ。（３）热稳定性：将发酵液分别在４０、５０、６０、７０、８０、９０、１００℃水浴３０ｍｉｎ或１２１℃处理２０ｍｉｎ后，冷却至室温（２５℃）。（４）蛋白酶稳定性：在发酵液中分别加入蛋白酶Ｋ或胰蛋白酶，使酶终浓度为１ｍｇ／ｍＬ，３７℃反应２ｈ。
１．２．２ 抑菌活性的检测 ＰＤＡ平板中加入１．５ｍＬ发酵滤液，用涂布器充分涂抹均匀后在超净工作台中晾干备用。以尖孢镰刀菌为供试菌、以未经处理的发酵滤液为阳性对照、以清水为阴性对照，２５℃倒置培养５ｄ后，用十字交叉法测量菌落直径，进行２组试验，每组试验设置４次重复。
１．２．３ 数据处理分析 根据下列公式计算菌株ＮＺ－４发酵液对尖孢镰刀菌的相对抑制率：抑制率＝（对照菌落直径－处理菌落直径）／（对照菌落直径－菌饼直径）×１００％；
相对抑菌活性＝处理抑制率／阳性对照抑制率×１００％。
１．３ 菌株ＮＺ－４水解酶的检测
将活化的菌株ＮＺ－４点接到对应的酶鉴定或ＣＡＳ培养基上，２８℃恒温培养３～５ｄ，观察是否产生透明圈或产生凹陷。
１．４ 菌株ＮＺ－４的全基因组测序和分析１．４．１ 菌株ＮＺ－４全基因组测序 收集对数生长期的菌株ＮＺ－４菌体，提取基因组ＤＮＡ。利用ＩｌｌｕｍｉｎａＮｏｖａＳｅｑ和ＰａｃＢｉｏＳｅｑｕｅｌ平台分别测序，对测序数据进行质控分析，利用ＨＧＡＰ［２１］、ＣＡＮＵ［２２］拼装三代数据得到ｃｏｎｔｉｇ序列。将二代的数据用ｐｉｌｏｎ［２３］软件对ｃｏｎｔｉｇ进行校正后获得完整序列。测序和拼装委托上海派森诺生物科技股份有限公司进行。
１．４．２ 基因组系统发育分析 使用ｏｒｔｈｏｆｉｎｄｅｒ［２４］获得与菌株ＮＺ－４近缘的１２个芽孢杆菌的单拷贝同源蛋白序列文件，使用ｍｕｓｃｌｅ［２５］进行多序列比对。后使用ＭｅｇａＸ［２６］默认参数构建ＮＪ（ｎｅｉｇｈｂｏｕｒ－ｊｏｉｎｉｎｇ）发育树，ｂｏｏｔｓｔｒａｐ法重复１０００次检验进化树。
１．４．３ 编码蛋白基因的分析 将菌株ＮＺ－４的完整序列上传至ＮＣＢＩ数据库（ｈｔｔｐｓ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ｇｅｎｏｍｅ／ａｎｎｏｔａｔｉｏｎ＿ｐｒｏｋ／）ＰｒｏｋａｒｙｏｔｉｃＧｅｎｏｍｅＡｎｎｏｔａｔｉｏｎＰｉｐｅｌｉｎｅ（ＰＧＡＰ）进行功能注释。利用ＰＳＩｂｌａｓｔｖ２．９．０［２７］将菌株ＮＺ－４的蛋白序列分别与ＣＯＧ［２８］和ＣＡＺｙ［２９］数据库进行比对，获得功
能注释信息。利用ａｎｔｉＳＭＡＳＨ６．０［３０］
对ＮＺ－４菌
株中次生代谢物合成基因组簇进行预测。２ 结果与分析
２．１ 菌株ＮＺ－４无菌发酵滤液对尖孢镰刀菌的抑菌稳定性
菌株ＮＺ－４发酵液在５０℃处理后抑菌活性为８０．７％，当超过６０℃后，随着温度的升高抑菌
活性逐渐下降（图１－Ａ），１００℃处理后抑菌活性约６０ ０％，１２１℃处理后几乎丧失抑菌活性。菌株
ＮＺ－４发酵滤液抑菌活性的最适ｐＨ值为６～９（图１－Ｂ），ｐＨ值＝３处理后抑菌活性仍在８５％以上，当ｐ
川岛芳子生死之谜大揭秘Ｈ值超过１０后，抑菌活性显著降低。菌株ＮＺ－４的抑菌活性物质对紫外线不敏感，随着紫外线照射时间的延长，其抑菌活性有小幅下降（图１－Ｃ），紫外线照射２０ｈ后，其抑菌活性仍达７６．９％。用蛋白酶Ｋ处理２ｈ后，发酵液抑菌活性显著降低，说明菌株ＮＺ－４发酵液中含有抗菌蛋白（图１－Ｄ）。但发酵液经胰蛋白酶处理后，抑菌活性没有明显改变，表明菌株ＮＺ－４
的抑菌物质对胰蛋白酶不敏感。
２．２ 菌株ＮＺ－４抑菌促生活性物质的测定
参考不同活性酶的选择性鉴定培养基，对菌株ＮＺ－４产酶活性进行检测。结果显示，菌株ＮＺ－４可在几丁质酶鉴定培养基（图２－Ａ）、纤维素酶鉴定培养基（图２－Ｂ）、蛋白酶鉴定培养基（图２－Ｃ）、酪蛋白酶鉴定培养基（图２－Ｄ）和ＣＡＳ培养基（图２－Ｅ）上生长并产生透明圈，能够在果胶酶鉴定培养基上生长，并使培养基产生凹陷（图２－Ｆ）。表明菌株ＮＺ－４能够产生几丁质酶、纤维素酶、蛋白酶、酪蛋白酶、果胶酶并能合成铁载体。该菌株可能通过分泌以上酶类物质，加速土壤中的多种碳水化合物的分解，竞争环境中稀缺的铁元素，发挥其拮抗和促生作用。
２．３ 菌株ＮＺ－４基因组概况和系统发育分析２．３．１ 菌株ＮＺ－４基因组基本概况 菌株ＮＺ－４基因组全长４０８５４０７ｂｐ，为环状染体，其ＧｅｎＢａｎｋ登录号为ＣＰ０７６１１９。ＧＣ含量为４６．３１％，包含４０１４个编码基因，３８９６个蛋白基因序列（ｃｏｄｉｎｇｓｅｑｕｅｎｃｅｓ，ＣＤＳ），６５个假基因（ｐｅｓｕｄｏｇｅｎｅ），２７个核糖体ＲＮＡ（ｒＲＮＡ），８６个转运ＲＮＡ（ｔＲＮＡ）和５个非编码ＲＮＡ（ｎｃＲＮＡ），基因组圈图如图３所示。 菌株ＮＺ－４编码基因的ＣＯＧ注释结果（图３）显示，共有３２１６个基因获得功能注释，注释结果共分为２４类。其中参与次生代谢产物生物合成、运输和分解代谢功能（Ｑ）的基因共有８９个，参与氨基酸转运与代谢（Ｅ）的基因３１２个，参与脂类转运与代
谢（Ｉ）的基因１５７个，参与无机离子转运与代谢（Ｐ）的基因１８７个，参与碳水化合物转运与代谢（Ｇ）的基因２８１个，防御机制（Ｖ）相关基因１１２个。２．３．２ 基于基因组核心基因
的系统发育分析及种类鉴定 基于菌株ＮＺ－４和参考菌株的１３９８个单拷贝核心同源基因构建系统发育树（图４），结果表明，菌株ＮＺ－４与贝莱斯芽孢杆菌ＣＣ０９、ＤＳＹＺ、ＦＺＢ４２、ＣＢＭＢ２０５、ＣＡＵＢ９４６和ＳＣＧＢ１等菌株聚为一个大的分支；而与解淀粉芽孢杆菌（Ｂａｃｉｌｌｕｓａｍｙｌｏｌｉｑｕｅｆａｃｉｅｎｓ）、枯草芽孢杆菌（Ｂａｃｉｌｌｕｓｓｕｂｔｉｌｉｓ）、地衣芽孢杆菌（Ｂａｃｉｌｌｕｓｌｉｃｈｅｎｉｆｏｒｍｉｓ）、短小芽孢杆菌（Ｂａｃｉｌｌｕｓｐｕｍｉｌｕｓ）和蜡样芽孢杆菌（Ｂａｃｉｌｌｕｓｃｅｒｅｕｓ）聚为不同的分支，表明菌株ＮＺ－４在基因组水平上与贝莱斯芽孢杆菌亲缘关系最近。因此基于基因组水平，菌株ＮＺ－４被鉴定为贝莱斯芽孢杆菌（Ｂａｃｉｌｌｕｓｖｅｌｅｚｅｎｓｉｓ）。
２．４ 菌株ＮＺ－４抗生素合成相关基因簇的预测芽孢杆菌抗菌活性与其产生的次生代谢产物密切相关，本研究利用ａｎｔｉＳＭＡＳＨ６．０对菌株ＮＺ－４基因组进行分析，结果显示，菌株ＮＺ－４编码１４种次生代谢产物合成基因簇。其中有９种与ＭｉＢｉＧ中已经描述过的次生代谢产物合成基因簇基本一致（图５），它们分别与表面活性素（ｓｕｒｆａｃｔｉｎ）、大环内酰亚胺Ｈ（ｍａｃｒｏｌａｃｔｉｎＨ）、杆菌素（ｂａｃｉｌｌａｅｎｅ）、杆菌霉素Ｄ（ｂａｃｉｌｌｏｍｙｃｉｎＤ）、泛革素（ｆｅｎｇｙｃｉｎ）、地非西丁（ｄｉｆｆｉｃｉｄｉｎ）、儿茶酚型铁载体（ｂａｃｉｌｌｉｂａｃｔｉｎ）、解淀粉芽孢杆菌素（ａｍｙｌｏｃｙｃｌｉｃｉｎ）和溶杆菌素（ｂａｃｉｌｙｓｉｎ）的合成相关。此外，菌株ＮＺ－４还编码５种没有比对到ＭｉＢｉＧ数据库或相似度比较低、可能为该菌株所特有的次生代谢产物。
这些次生代谢产物中有５种具有抗真菌活性，包括３种脂肽类化合物，分别为表面活性素、泛革素和杆菌霉素，１种二肽溶杆菌素和１种儿茶酚型铁载体。此外，菌株ＮＺ－４还能合成４种对细菌有抑制作用的次生代谢产物，包括３种聚酮类化合物，分别为杆菌素、大环内酰亚胺Ｈ和地非西丁，以及１种环肽解淀粉芽孢杆菌素。菌株ＮＺ－４编码多种抗生素合成基因簇，表明该菌株可能通过这些抗菌产物发挥其抑菌功能。
２．５ 菌株ＮＺ－４碳水化合物活性酶（ＣＡＺｙ）及与生防有关酶类的预测
菌株ＮＺ－４编码１０５个ＣＡＺｙ酶类基因（图６－Ａ），包括４３个糖苷水解酶（ｇｌｙｃｏｓｉｄｅｈｙｄｒｏｌａｓｅ，ＧＨｓ）、３６个糖基转移酶（ｇｌｙｃｏｓｙｌｔｒａｎｓｆｅｒａｓｅ，ＧＴｓ）、１６个碳水化合物酯酶（ｃａｒｂｏｈｙｄｒａｔｅｅｓｔｅｒａｓｅｓ，ＣＥｓ）、６个碳水化合物结合模块（ｃａｒｂｏｈｙｄｒａｔｅ－ｂｉｎｄｉｎｇｍｏｄｕｌｅ，ＣＢＭｓ）（其中４个碳水化合物结合模块又属于糖苷水解酶）、３个多糖裂解酶（ｐｌｏｙｓａｃｃｈａｒｉｄｅｌｙａｓｅｓ，ＰＬｓ）和５个辅助酶类（ａｕｘｉｌｉａｒｙａｃｔｉｖｉｔｉｅｓ，ＡＡｓ）。其中ＧＨｓ家族数目最多，占整个基因家族的４１ ０％，ＧＴｓ家族次之，约占３４．３％（图６－Ｂ）。其中含有信号肽的ＣＡＺｙｍｅ蛋白有２４个（２２．９％），说明它们具有蛋白质跨膜分泌功能。分泌的ＣＡＺｙｍｅ中数量最多的是ＧＨｓ，４３个ＧＨｓ中有１６个含有信号肽，而ＡＡｓ和ＧＴｓ中各有１个成员含有信号肽。在这些基因家族中，ＧＨ５、ＧＨ１１、ＧＨ２６、ＧＨ４３、ＧＨ５１和ＧＨ５３等与纤维素降解有关；ＧＨ１３和ＧＨ１２６与淀粉水解相关；ＧＨ１８、ＧＨ２３、ＣＥ９和ＣＢＭ５０等与几丁质降解相关；ＧＨ４３、ＣＥ４、ＣＥ６和ＣＥ７等与木聚糖降解相关；ＧＨ４３、ＰＬ１、ＰＬ９等与果胶降解酶类相关；ＧＨ２３和ＧＨ７３等与肽聚糖降解相关；ＧＨ３、ＧＨ１６和ＧＨ３０等与葡聚糖酶降解相关。说明菌株ＮＺ－４具有降解纤维素、淀粉、几丁质、果胶、肽聚糖和葡聚糖等物质的潜力。
３ 讨论与结论
芽孢杆菌作为生防菌如贝莱斯芽孢杆菌ＦＺＢ４２等，已在多种作物土传病害的防治过程中取得了显著成效［７，９，３１］。菌株ＮＺ－４是本试验室从马铃薯根际土壤中分离得到的一株贝莱斯芽孢杆菌，对多种重要农作物病原菌有良好的抑制作用［１２］。本试验在此基础上，进一步研究了菌株ＮＺ－４发酵滤液的抑菌活性稳定性，并通过对该菌株的基因组测序和分析，鉴定到该菌株发挥促生和生防作用的潜在物质，为该菌株的开发利用和机理研究奠定了基础。
本研究发现，菌株ＮＺ－４发酵液经过不同温度、ｐＨ值、紫外照射和蛋白酶处理后，其抑菌活性在极端条件下会发生变化，说明该菌株的抑菌物质会受到外界极端环境的影响，该特性与已报道的其他芽孢杆菌相似［３２－３４］。菌株ＮＺ－４的发酵液在常规条件下具有良好的稳定性，但在极端条件下的稳定性变化，可能与发酵液中抗菌物质的种类或活性变化有关。总之，该菌株发酵液稳定性较强，能适应田间不断变化的环境条件，具备在农业生产中防治枯萎病等土传病害的潜力。
利用多种碳源快速繁殖和成功定殖是生防菌发挥生防功能的前提［３５］。菌株ＮＺ－４能够产生并
发布评论