Identification of 29 strains of entomogenous fungi and their toxicity to Bemisia tabaci
-
摘要:目的
对29株虫生真菌菌株进行分类学鉴定并筛选出对烟粉虱Bemisia tabaci具有高毒力的菌株。
方法采用形态学及ITS区、TEF区、Bloc区基因序列相似性分析鉴定菌株，并采用浸渍法测定真菌对烟粉虱的毒力。
结果29株虫生真菌有26株符合白僵菌Beauveria的特征，菌株SB003、SP665和SP670被鉴定属于棒束孢属Isaria；利用TEF区序列、Bloc区序列、TEF-Bloc联合序列构建的系统发育树，结果显示29株菌株共包括24株球孢白僵菌B. bassiana、2株假性球孢白僵菌B. pseudobassiana和3株环链棒束孢I. cateinannulata。烟粉虱2龄若虫的死亡率随孢子浓度的增加而增加，不同菌株致病力不同；用1×108 mL−1的孢子悬浮液处理烟粉虱2龄若虫7 d，致死率最高的菌株为SP433，其次为SB009，分别为87.37%和82.93%。
结论本研究可为虫生真菌的分类提供理论基础，筛选出的高致病力菌株SP433和SB009可为烟粉虱的生物防治提供参考。
Abstract:ObjectiveTo identify 29 strains of entomogenous fungi and select high virulent strains against Bemisia tabaci.
MethodIdentification of fungi were based on morphological analysis and genetic homology analysis of sequences in ITS, TEF and Bloc regions. The fungal toxicity to B. tabaci was evaluated by the immersion method.
ResultTwenty-six of 29 entomogenous fungal strains had the characteristics of Beauveria and the other three strains of SB003, SP665 and SP670 were identified as Isaria. Phylogenetic tree was constructed using TEF region sequence, Bloc region sequence and TEF-Bloc combined sequence, which showed that 29 strains included 24 strains of Beauveria bassiana, 2 strains of B. pseudobassiana, and 3 strains of Isaria cateinannulata. The mortality of the 2nd instar nymph of B. tabaci increased with the increase of spore concentration, and the pathogenicity of different strain was different. The 2nd instar nymphs of B. tabaci were treated with 1×108 mL−1 spore suspension for 7 days, and the SP433 strain had the highest lethality rate of 87.37%, followed by SB009 with a lethality rate of 82.93%.
ConclusionThis study can provide a theoretical basis for the classification of entomophytic fungi. The highly pathogenic strains of SP433 and SB009 were screened out, which can provide a reference for the biological control of B. tabaci.
-
虫生真菌是一类重要的生物资源，也是生态系统的重要组成部分，其寄主广泛，包括半翅目的蚜科、飞虱科、叶蝉科、沫蝉科、粉虱科和盲蝽科，鳞翅目、鞘翅目、双翅目和缨翅目等昆虫[1]。虫生真菌以种类多、安全有效、应用期长、不伤害天敌、不易产生抗性及能快速大量生产等优点著称，在生物防治中有着无可比拟的反复侵染性和生产便利性[2]。白僵菌Beauveria寄主范围广泛，包括重要的农业和森林害虫[3-5]、病媒介昆虫[6-7]和蜱螨类[8]，易培养，且对温血动物和植物无害，另外，不仅对成虫治病力强，还能侵染幼虫、蛹、成虫等虫态，对下一代有持续作用[9-10]。
烟粉虱Bemisia tabaci属于半翅目Homoptera粉虱科Aleyrodidae小粉虱属Bemisia，是重要的农业害虫之一，可通过直接刺吸及有效传播病毒而对作物造成损伤[11-13]。成虫高度杂食性，侵染1 000多种植物，并可传播菜豆金色黄花叶病毒属Begomovirus、毛形病毒属Crinivirus、香石竹潜隐病毒属Carlavirus、番茄托拉多病毒属Torradovirus和甘薯病毒属Ipomovirus等300多种病毒，造成的经济损失每年达数10亿美元[14-15]。烟粉虱若虫拥有蜡质外壳，常固定吸附在作物表面，仅以刺吸式口器危害作物，化学杀虫剂很难渗入其体内，为了达到防治效果往往加大药量，这样不但增加经济成本而且污染环境，还会提高害虫的抗药性[16]。目前烟粉虱对除虫菊酯类、有机磷类等杀虫剂均已产生了较为强烈的抗药性[17]，因此，研究有效的烟粉虱生物防治技术尤为重要。
本研究以华南农业大学生物防治教育部工程研究中心保存的29株菌株为研究对象，对菌株进行形态学及分子生物学分类鉴定，测定不同菌株对烟粉虱的室内致病力，最终筛选出致病力较强的优良菌株，以期为烟粉虱防治提供后备的生产菌株。
1. 材料与方法
1.1 供试材料
虫生真菌共计29株，均为华南农业大学生物防治教育部工程研究中心继代保存的菌株。具体信息详见表1。
表 1 供试菌株信息Table 1. Information of test strains菌株
Strain采样地
Sampling place采样时间
Sampling timeSB003 湖北省神农架林区 Shennongjia forest area, Hubei Province 2007−07 SB004 甘肃省天水市麦草沟自然保护区
Maicaogou nature reserve, Tianshui, Gansu Province2007−08 SB006 甘肃省临夏回族自治州莲花山自然保护区
Lianhuashan nature reserve, Linxia Hui Autonomous Region, Gansu Province2007−08 SB009 四川省成都市青羊区文殊院 Wenshuyuan, Qingyang district, Chengdu, Sichuan Province 2007−10 SB010 云南省昆明市 Kunming, Yunnan Province 2007−09 SB015 河南省鹤壁市 Hebi, Henan Province 2007−10 SB026 四川省雅安市宝兴县蜂桶寨
Fengtong village, Baoxing county, Ya’an, Sichuan Province2007−07 SB032 韩国济州岛天地渊瀑布 Tiandiyuan waterfall, Jeju Island, South Korea 2008−07 SB035 云南省红河州石屏县 Shiping county, Honghe Prefecture, Yunnan Province 2009−05 SB036 内蒙古自治区呼伦贝尔市新巴尔虎左旗
Xinbalhuzuoqi, Hulunbuir, Inner Mongolia Autonomous Region2009−08 SB037 内蒙古自治区兴安盟扎赉特旗
Zhalaitqi, Xing’an League, Inner Mongolia Autonomous Region2009−08 SB038 内蒙古自治区锡林郭勒盟正镶白旗
Zhengxiangbaiqi, Xilinguole League, Inner Mongolia Autonomous Region2009−08 SB039 内蒙古自治区阿拉善盟贺兰山国家级自然保护区
Helan Mountain National Nature Reserve, Alxa League, Inner Mongolia Autonomous Region2009−08 SB041 内蒙古自治区锡林郭勒盟正镶白旗
Zhengxiangbaiqi, Xilinguole League, Inner Mongolia Autonomous Region2009−08 SB043 西藏自治区林芝市波密县古乡
Gu village, Bomi county, Linzhi, Tibet Autonomous Region2009−09 SB050 广东省广州市白云山 Baiyun Mountain, Guangzhou, Guangdong Province 2009−10 SB051 广东省河源市龙川县 Longchuan county, Heyuan, Guangdong Province 2010−05 SB057 贵州省雷山县方祥乡 Fangxiang village, Leishan county, Guizhou Province 2010−07 SB062 内蒙古自治区阿拉善盟古拉本敖包镇
Gulaben Aobao Town, Alxa League, Inner Mongolia Autonomous Region2010−08 SB063 内蒙古自治区阿拉善盟贺兰山脉中部哈拉乌
Harau, central Helan Mountain, Inner Mongolia Autonomous Region2010−08 SB671 云南省镇沅彝族哈尼族拉祜族自治县
Zhenyuan Yi Hani Lahu Autonomous county, Yunnan Province2009−10 SB672 福建省厦门市鼓浪屿 Gulangyu, Xiamen, Fujian Province 2010−09 SB673 云南省迪庆州维西县保和镇 Baohe town, Weixi county, Diqing Prefecture, Yunnan Province 2010−09 SB674 广东省广州市华南农业大学林学与风景园林学院
College of forestry and landscape architecture, South China Agricultural University, Guangzhou,
Guangdong Province2010−10 SP016 湖北省神农架林区老君山 Laojun Mountain, Shennongjia forest area, Hubei Province 2007−08 SP031 甘肃省张掖市肃南县东柳沟 Dongliugou, Sunan county, Zhangye, Gansu Province 2007−08 SP433 内蒙古自治区呼伦贝尔市牙克石市
Yakeshi, Hulunbuir, Inner Mongolia Autonomous Region2009−08 SP665 云南省楚雄彝族自治州峨绿公园西灵宫
Xiling Palace, Erlu Park, Chuxiong Yi Autonomous Prefecture, Yunnan Province2009−09 SP670 湖北省宜昌市蛤蟆泉 Toad Spring, Yichang, Hubei Province 2009−10 棉花“鲁棉研32号”购自山东省潍坊市种子管理站。自然条件下，用经高压灭菌的营养土钵栽培棉花种子，生长期间棉花苗置于干净的养虫笼(60 cm×60 cm×60 cm)内，待苗长至6~8叶期时用于试验。
B型烟粉虱种群2005年采自华南农业大学教学实习农场，寄主植物为番茄。采回后烟粉虱被隔离饲养在生物防治教育部工程研究中心的棉花寄主上继代繁殖。饲养条件为温度(27±1)℃，光周期为14 h光∶10 h暗，相对湿度(75±5)%。
Ezup柱试真菌基因组DNA抽提试剂盒[生工生物工程(上海)股份有限公司]、吐温−80(φ为0.05% )、95%(φ)乙醇溶液、75%(φ)乙醇溶液、葡萄糖、琼脂糖、戊二醛、丙酮、磷酸缓冲液、DL2000 DNA Marker。
1.2 方法
1.2.1 虫生真菌的形态学鉴定
将菌株转接到新的PDA培养皿25 ℃条件下培养3~5 d，取菌丝分布均匀的供试菌一皿，用直径1.5 cm的打孔器打菌饼，将菌饼置于PDA平板中央进行接种，25 ℃光照培养箱内培养。用十字交叉法，分别于第3、5、7、9 和10 天测定各菌株的纵横菌落直径，每株重复测定5~6个皿，以加入无菌水的PDA培养基作为对照。同时，10 d后用相机拍摄菌落正反面。
将上一步骤培养10 d形成的菌落全部置于50 mL 0.05%(φ)的吐温−80溶液中，充分振荡，取200 µL孢子悬浮液，在血球计数板上镜检各菌株的孢子萌发情况。每处理重复3次，以加入无菌水的液体培养基作为空白对照。菌株产孢量计算公式为：
1 cm2 菌株产孢量=(平均每小格孢子数×4×106×稀释倍数)/(3.14×0.09)。
进行孢子的扫描电镜观察。
1.2.2 虫生真菌的分子生物学鉴定
用Ezup柱试真菌基因组DNA抽提试剂盒提取总DNA。
ITS序列、Bloc序列、TEF序列的PCR扩增均采用通用引物，详细引物信息见表2。PCR反应体系为25 μL体系，其中ddH2O 10 μL，2×Taq PCR StarMix 12 μL，上下游引物各1 μL，DNA模板1 μL。ITS扩增条件：94 ℃ 5 min；94 ℃ 1 min，55 ℃ 1 min，72 ℃ 1 min，35个循环；72 ℃延伸10 min。Bloc扩增条件：95 ℃ 5 min；95 ℃ 30 s，54 ℃ 30 s，72 ℃ 80 s，40个循环；72 ℃延伸15 min。TEF扩增条件：TEF采用降落PCR的扩增方式，95 ℃ 5 min；95 ℃ 30 s，66 ℃ 30 s，72 ℃ 1 min，每个循环退火温度降低1 ℃，11个循环；95 ℃ 30 s，56 ℃ 30 s，72 ℃ 1 min，35个循环；72 ℃延伸10 min。将目的条带大小合适样品的纯化回收产物送上海美吉生物医药科技有限公司进行序列测定。利用软件Geneious进行序列拼接，进行Blast序列相似性搜索，采用软件MEGA 7.0构建系统进化树。
表 2 虫生真菌DNA的PCR扩增引物Table 2. PCR amplification primers for entomogenous fungal DNA基因序列 Gene sequence 引物 Primer 引物序列(5′→3′) Primer sequence ITS ITS4F TCCTCCGCTTATTGATATGC ITS5R GGAAGTAAAAGTCGTAACAAGG Bloc B5.1F CGACCCGGCCAACTACTTTGA B3.1R GTCTTCCAGTACCACTACGCC TEF 983F GCYCCYGGHCAYCGTGAYTTYAT 2218R ATGACACCRACRGCRACRGTYTG 1.2.3 虫生真菌对烟粉虱的毒力测定
在PDA平板上培养7 d后，充分产孢的虫生真菌分生孢子用0.05%(φ)吐温−80无菌水洗脱，搅拌，充分振荡，过滤，用血球计数板计数，测量母液的浓度，配制成1×108 mL−1的孢子悬浮液，再稀释成1×107 mL−1的孢子悬浮液备用。
将带有烟粉虱2龄若虫的棉花叶片背面浸于配好的孢子悬浮液中，30 s后取出，待叶片自然晾干后，将叶片背面朝上置于琼脂培养基中，用保鲜膜封口并扎孔通气，将培养皿倾斜放置，置于人工气候箱中，饲养条件为温度(27±1)℃，光周期14 h光∶10 h暗，相对湿度(75±5)%。以0.05%(φ)吐温−80无菌水为空白对照。每个处理重复4次，每个重复100头若虫。于7 d后检查烟粉虱的死亡情况并记录死亡数，计算死亡率。
1.3 数据处理
采用SPSS17.0与Excel软件进行数据处理及统计分析，Duncan’s进行差异显著性检验。
2. 结果与分析
2.1 虫生真菌的形态学观察
2.1.1 虫生真菌的菌落形态
培养第10天分别观察及拍摄PDA培养基上各菌株菌落的正面及背面形态图(图1)。
![]()
图 1 菌株的菌落形态(奇数为正面，偶数为反面)1、2：SB003；3、4：SB004；5、6：SB006；7、8: SB009；；9、10：SB010；11、12：SB015；13、14: SB026；15、16：SB032；17、18：SB035；19、20：SB036；21、22：SB037；23、24：SB038；25、26：SB039；27、28：SB041；29、30：SB043；31、32：SB050；33、34：SB051；35、36：SB057；37、38：SB062；39、40：SB063；41、42：SB671；43、44：SB672；45、46：SB673；47、48：SB674；49、50：SP016；51、52：SP031；53、54：SP433；55、56：SP665；57、58：SP670Figure 1. Colony morphology of strain(the odd number showed the front and the even number showed the back )菌株SB003菌落呈白色，中央气生菌丝突起，向外气生菌丝减少，变薄，背面橘黄色，培养基平板不变色(图1-1、1-2)；菌株SB004菌落呈放射状有褶皱，中间绒毛状，背面橘黄色(图1-3、1-4)；菌株SB006菌落由绒毛状变粉末状，中间为凸起粉末状(图1-5、1-6)。
菌株SB009菌落呈绒毛状，表面白色，背面橘黄色(图1-7、1-8)；菌株SB010菌落由绒毛状变粉末状，背面边缘呈绒毛状，表面白色，背面橘黄色(图1-9、1-10)；菌株SB015菌落由绒毛状变粉末状，边缘无绒毛，背面白色或淡黄色，后期中间凸起厚粉末(图1-11、1-12)。
菌株SB026(图1-13、1-14)菌落形态特征同菌株SB015；菌株SB032(图1-15、1-16)、SB036(图1-19、1-20)、SB037(图1-21、1-22)、SB038(图1-23、1-24)菌落初期为绒毛状，中央气生菌丝突起、疏松，后期中间有凸起厚粉末，边缘呈绒毛状，有褶皱，表面白色或乳白色，培养基平板不变色，此外，菌株SB037、SB038菌落有液滴；菌株SB035(图1-17、1-18)菌落由绒毛状变粉末状，表面白色，背面呈橘黄色。
菌株SB039(图1-25、1-26)、SB041(图1-27、1-28)菌落絮状，表面白色，背面淡黄色；菌株SB043培养10 d的菌落由绒毛变粉末状，边缘呈绒毛状，背面呈淡黄色(图1-29、1-30)。
菌株SB050菌落由绒毛状变为粉末状，表面白色，培养基背面白色或淡黄色(图1-31、1-32)；菌株SB051菌落前期绒毛状，后期粉末状，中间凸起厚粉末，边缘呈绒毛状，菌落放射状，背面淡黄色(图1-33、1-34)；菌株SB057菌落由绒毛状变粉末状，边缘无绒毛，背面白色或淡黄色(图1-35、1-36)。
菌株SB062(图1-37、1-38)、SB673(图1-45、1-46)菌落形态同SB032；菌株SB063(图1-39、1-40)、SB671(图1-41、1-42)、SB672(图1-43、1-44)和SB674(图1-47、1-48)菌落形态同SB050。
菌株SP016(图1-49、1-50)菌落形态同SB043；菌株SP031菌落中间绒毛状，表面白色，背面淡黄色，边缘绒毛状(图1-51、1-52)；菌株SP433(图1-53、1-54)菌落形态同SB032。
菌株SP665(图1-55、1-56)、SP670(图1-57、1-58)菌落绒毛状，中央气生菌丝突起、疏松，向外气生菌丝减少、变薄，中部淡黄色，边缘白色，背面橘黄色或淡黄色，培养基平板不变色。
2.1.2 孢子的扫描电镜观察
菌株SB009的孢子呈卵形和近球形，孢子直径1.8~2.2 µm (图2a)。菌株SB026的孢子呈卵形，孢子直径1.8~2.4 µm (图2b)。菌株SB038产孢细胞的基部膨大，产孢轴“Z”型弯曲，分生孢子有球形和卵形，孢子直径1.7~2.6 µm (图2c)。SB041 孢子形状和大小与SB009比较类似，孢子卵形，部分为近球形，直径1.6~2.5 µm (图2d)。SP031孢子有球形，也有椭圆形，孢子大小1.3~1.5 µm (图2e)。SP665分生孢子梗不规则分枝，长短不一，基部拟椭圆形膨大，变细成明显的颈部。分生孢子球形或椭圆，孢子直径1.4~2.2 µm (图2f)。
2.1.3 菌落生长速率
虫生真菌的各菌落生长速率见图3，菌株培养前期(3 d)，菌株的菌落直径差异不明显，菌株SB010和SP031的直径大于25.0 mm，菌株SB038、SB035、SB039的直径小于20.0 mm，其他菌株的菌落直径为20.8~24.3 mm，菌株间无明显差异；随着培养时间的延长，菌株的菌落直径差异增加，菌株培养5d时，菌株SB010、SB032、SP031、SB006的菌落直径大于35.0 mm，菌株SB010的菌落直径最大，为37.9 mm，菌株SB039与SB035的菌落直径偏小，不到25.0 mm；培养7 d时，菌落直径达40.0 mm以上的菌株占供试菌株的一半，而菌株SB035与SB039的菌落直径却不到30.0 mm；培养9 d时，菌株SP031、SB010、SB006、SB032、SB671、SB062、SB037、SB672的菌落直径达50.0 mm及以上，而菌株SB035的菌落直径仍不到30.0 mm；培养10 d时，菌株SB010与SP031已达60.0 mm以上，分别为62.7 和61.3 mm，生长速率最快；菌株SB041、SB039和SB035的菌落直径不到40.0 mm，菌株SB035的菌落直径仍为各菌株中菌落直径最小的，仅有30.3 mm。
2.1.4 菌株产孢量
菌株产孢量见表3，29株虫生真菌的产孢量存在明显差异。其中，菌株SB032产孢量最大，其次依次分别为SB039、SB038、SB004、SB006、SB050、SP433、SB036和SB043，产孢量均大于108 mL−1。产孢量最少的为SB003和SB041，产孢量不到3.50×107 mL−1。
表 3 虫生真菌的产孢量1)Table 3. Spore productions of entomogenous fungi菌株
Strain产孢量/(×107 mL−1)
Spore production菌株
Strain产孢量/(×107 mL−1)
Spore production菌株
Strain产孢量/(×107 mL−1)
Spore productionSB003 3.25±0.59f SB037 9.86±0.69abcdef SB671 4.64±0.71ef SB004 14.27±2.63abc SB038 14.78±2.06ab SB672 8.35±1.72abcdef SB006 13.65±1.59abcd SB039 16.38±2.83a SB673 5.15±1.10def SB009 9.42±1.19abcdef SB041 3.48±0.47f SB674 6.61±1.75bcdef SB010 6.72±1.92bcdef SB043 11.16±1.48abcdef SP016 4.04±0.49ef SB015 8.29±0.74abcdef SB050 12.64±3.07abcde SP031 9.32±1.73abcdef SB026 9.19±1.73abcdef SB051 6.54±0.47bcdef SP433 12.54±2.36abcde SB032 16.56±1.59a SB057 9.30±0.37abcdef SP665 5.57±0.89cdef SB035 3.84±0.56ef SB062 6.02±1.76bcdef SP670 7.26±0.49bcdef SB036 11.23±2.66abcdef SB063 7.88±1.59abcdef 1)数据后不同小写字母者表示不同菌株间菌落产孢量差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters indicate significant differences in spore production among different strains (P<0.05, Duncan’s test)2.2 虫生真菌的分子生物学鉴定
将ITS区、TEF区、Bloc区的序列拼接后去除ITS4、ITS5、2218R、983F、B5.1F和B3.1R的引物序列，ITS区的片段为550~600 bp，TEF区为1 000 bp左右，Bloc区的片段为1 500 bp左右，在GenBank上进行BLAST比对，比对出的相似性最高的菌株结果见表4，最终比对得出分别属于球孢白僵菌、假球孢白僵菌和环链棒束孢。
表 4 供试菌株ITS区、TEF区和Bloc区的序列比对结果Table 4. Results of comparing the ITS regions, TEF regions and Bloc regions of test strains真菌
Fungus菌株登录号 Accession number of strain ITS区
ITS regionTEF区
TEF regionBloc区
Bloc region球孢白僵菌 Beauveria bassiana KX901309 JQ867178 JQ867116 球孢白僵菌 B. bassiana AY531976 JQ867172 JQ867106 球孢白僵菌 B. bassiana KX901310 JQ867171 JQ867109 球孢白僵菌 B. bassiana GU233705 JQ867157 JQ867118 球孢白僵菌 B. bassiana GU233701 JQ867151 JQ867125 球孢白僵菌 B. bassiana AY883691 AY883788 球孢白僵菌 B. bassiana KU725693 KU725706 假球孢白僵菌 Beauveria pseudobassiana MH259852 HQ880999 HQ880727 假球孢白僵菌 B. pseudobassiana LC163567 LC163595 环链棒束孢 Isaria cateniannulata GU734760 GU734758 将供试的29株菌株的ITS序列、TEF序列和Bloc序列采用最大似然法建树建立系统发育树(图4)。ITS序列构建的系统发育树(图4a)分为2大分支，第1分支为假环链棒束孢ARSEF 6242与供试菌株SB003、SP665和SP670，因此可把它们归为环链棒束孢；第2大分支为白僵菌属，24株供试菌株与白僵菌属的不同球孢白僵菌类聚，为球孢白僵菌，但其中菌株SB039、SB041与Strain E51类聚1分支，为假球孢白僵菌。
![]()
图 4 基于不同序列以最大似然法构建的真菌系统发育树▲为供试菌株Figure 4. Phylogenetic tree of fungi constructed using the maximum likelihood method and based on different sequences▲ indicates test strainsTEF序列构建的系统发育树(图4b)分为3大分支，第1分支为棒束孢所在的分支，菌株SB003、SP665和SP670与菌株ARSEF 6240类聚，为环链棒束孢；第2分支由菌株F884、ARSEF 1855、SB039和SB041组成，为假球孢白僵菌；第3分支由供试的24株菌株和菌株RCEF 0006、RCEF 1276、RCEF 0301、RCEF1491、RCEF0333、ARSEF 1104和CHE-CNRCB 168组成，为球孢白僵菌。
SP665、SP670和SB003的Bloc区未能扩增成功。Bloc序列构建的系统发育树(图4c)分为2大分支，第1分支为假球孢白僵菌ARSEF 1855、F884所在分支，与该种类聚的也是菌株SB039和SB041；第2分支分由3个小分支构成，为球孢白僵菌。各种间的自展支持率基本大于50%。
TEF-Bloc序列构建的系统发育树(图4d)与单基因Bloc序列构建的发育树很相似，自展支持率不一样，但供试菌株与各菌种的类聚位置一致。
2.3 虫生真菌对烟粉虱的毒力测定
不同浓度的各菌株孢子悬浮液对烟粉虱2龄若虫的7 d毒力测定结果见表5，各菌株在同一处理时间下，致死率存在差异。当孢子悬浮液为1×107 mL−1时，菌株SB039的致死率最高(72.77%)，其次是SP433菌株(69.03%)；当孢子悬浮液为1×108 mL−1时，菌株SP433、SB009、SB050、SB036、SB043的致死率达80%以上，菌株SP433的致死率最高(87.37%)，致死率最小的是菌株SP031(52.85%)。同一菌株不同浓度，致死率随着孢子悬浮液浓度的增加而升高，而菌株SB039、SB063、SP031、SP670、SB671、SB035在孢子悬浮液为1×107和1×108 mL−1下处理烟粉虱若虫7 d后得到相似的致死率。
表 5 不同浓度的孢子悬浮液处理烟粉虱2龄若虫7 d的死亡率1）Table 5. Mortality of the 2nd instar nymphs of Bemisia tabaci infected with spore suspensions of different concentrations for 7 days% 菌株 Strain 1×107 mL−1 1×108 mL−1 菌株 Strain 1×107 mL−1 1×108 mL−1 CK 9.43±0.58d 11.73±2.00d SB043 47.67±4.00abc 80.33±3.93ab SB003 56.43±3.45abc 60.06±2.68bc SB050 49.47±1.47abc 81.30±4.86ab SB004 55.23±5.85abc 75.23±5.49abc SB051 59.67±4.17abc 70.87±5.90abc SB006 51.03±0.54abc 72.33±6.90abc SB057 51.67±3.55abc 76.00±2.34abc SB009 67.93±1.68ab 82.93±2.88ab SB062 35.81±3.54c 63.34±3.93abc SB010 47.90±4.42abc 70.20±4.94abc SB063 68.03±5.46ab 71.40±2.03abc SB015 64.30±7.27ab 77.07±4.95abc SB671 59.80±1.12abc 59.50±9.22bc SB026 51.77±4.15abc 67.80±3.0abc SB672 44.67±6.17bc 62.43±7.28bc SB032 56.57±2.95abc 77.69±3.88ab SB673 53.23±5.62abc 70.53±2.57abc SB035 53.07±2.37abc 59.99±2.77bc SB674 67.20±7.04ab 77.20±2.70abc SB036 55.33±6.36abc 80.67±4.41ab SP016 54.47±8.25abc 72.07±4.82abc SB037 50.17±7.69abc 74.00±3.21abc SP031 49.43±2.13abc 52.85±1.00c SB038 48.03±6.08abc 70.87±5.90abc SP433 69.03±3.76ab 87.37±0.43a SB039 72.77±2.14a 73.13±5.20abc SP665 49.70±4.65abc 64.90±3.26abc SB041 50.33±4.36abc 70.33±5.42abc SP670 65.77±0.50ab 72.40±1.56abc 1)相同浓度数据后不同小写字母者表示不同菌株间致病力差异显著(P<0.05，Duncan’s法)
1)Different lowercase letters after data at the same concentration indicate significant differences in pathogenicity among different strains (P<0.05, Duncan’s test)3. 讨论与结论
对本实验室长期保存的29株虫生真菌进行了形态学与分子生物学鉴定。供试的29株真菌菌株间的菌落直径、产孢量、孢子结构存在一定差异，菌株形态与以前描述相符，但菌株的形态特征存在交叉现象，不易区分。同一种分生孢子的大小或菌株形态在培养中会发生一定的变化[18]。在分子鉴定中，用ITS区序列鉴定白僵菌属内各种的效果不理想，ITS区序列的一致性程度很高，用ITS区序列单独建系统发育树时，我们发现供试菌株在系统发育树中不能以种为分支各个分开而是呈现交叉状态，对于需严格区分不同种或不常见的白僵菌的效果不理想，因此需要引入其他辅助基因序列[19]进行研究，而有报道称白僵菌属内有些种的TEF序列对比差异太小会导致难以区分、分类结果错误等问题[20-22]，因此多基因位点分析在近年来有成为研究物种差异主要方法的趋势。张云月等[23]鉴定白僵菌时，用ITS序列建构系统树时也出现了菌属内各种与球孢白僵菌系统交叉的现象，另外，他们还引入了β-tubulin序列，但也没能很好地区分白僵菌内种，验证了ITS序列用于白僵菌属内种分类的效果不理想。芦俊佳等[24]鉴定楚雄腮扁叶蜂寄生真菌时，用ITS序列比对时，发现相似性太高，也不能对菌株的分类做出正确的界定，又对其18 S序列进行了扩增，最终确定了菌株的分类位置。白僵菌属内种的菌株鉴定是否也可以采用ITS和18S序列来解决需要进一步的试验验证。
我们将TEF区、Bloc区的序列及TEF-Bloc联合序列构建系统发育树，运用了最大似然法构建系统发育树，结果显示29株菌株共包括24株球孢白僵菌、2株假性球孢白僵菌和3株环链棒束孢。另外，菌株SB003、SP670和SP655在Bloc区没能扩增出来，出现这种情况的原因还有待进一步的试验验证。分子生物学逐渐成为真菌鉴定的辅助手段，但分子生物学方法也会受基因库的完善程度、基因库中高度同源性系列的多少以及具体物种ITS区的可变程度等因素的影响[25-28]。球孢白僵菌在自然环境中的遗传变异较为广泛，且影响因素很多。从系统发育树来看，与地域环境的影响也有关，如在TEF-Bloc联合构建的系统发育树中，SB036、SB038和SB037聚在同一小支，而它们的采集地都是在内蒙古自治区内；3株从甘肃省采集的菌株SB006、SB004和SP031在ITS、TEF-Bloc的系统发育树都聚在同一支。但来自同一省份采集地却并不都聚在同一小分支，如SB010、SB035、SB671、SB673都采集于云南省，但在TEF-Bloc联合构建的系统发育树上，上述4株菌株各自属于不同的4个小分支。这与方志刚等[29]对不同球孢白僵菌菌株进行地区来源与遗传多样性分析得到的结果类似。
本试验结果显示，菌株SP433和SB009在烟粉虱防治中具有很高的应用潜力。在本研究中，在孢子悬浮液108 mL−1处理7 d，烟粉虱的死亡率为52.85%~87.37%，王慧[30]研究白僵菌对烟粉虱若虫的侵染效果，在孢子悬浮液108 mL−1处理10 d后，烟粉虱的死亡率为52.85%~87.37%；吴乾兴等[31]研究球孢白僵菌对烟粉虱若虫的侵染效果，在孢子悬浮液108 mL−1 处理6 d的烟粉虱死亡率为61.2%~87.9%；以上结果表明，同种的不同菌株可能有不同的致病效果。已有研究通过白僵菌侵染生长动力学试验证实影响毒力差异的一大重要因素是菌体增殖能力不同[32]，另外还可能与采集地点、菌株保存时间、菌株分离对象等有关。
研究白僵菌联合其他杀虫剂或寄生昆虫进行联合防治是研究的热点，银航等[33]使用球孢白僵菌与保幼激素联用，可显著提高对花椒蚜虫的致死率；张爽等[34]使用苏云金杆菌和球孢白僵菌复配来提高对小菜蛾幼虫的杀伤力；杨芷等[35]利用赤眼蜂搭载球孢白僵菌对亚洲玉米螟进行协同防控。未来白僵菌的应用范围将更加广泛、方法将更加多样，如何有效利用我们筛选出的高致病力的菌株还有待进一步的研究。
-
![]()
图 1 菌株的菌落形态(奇数为正面，偶数为反面)
1、2：SB003；3、4：SB004；5、6：SB006；7、8: SB009；；9、10：SB010；11、12：SB015；13、14: SB026；15、16：SB032；17、18：SB035；19、20：SB036；21、22：SB037；23、24：SB038；25、26：SB039；27、28：SB041；29、30：SB043；31、32：SB050；33、34：SB051；35、36：SB057；37、38：SB062；39、40：SB063；41、42：SB671；43、44：SB672；45、46：SB673；47、48：SB674；49、50：SP016；51、52：SP031；53、54：SP433；55、56：SP665；57、58：SP670
Figure 1. Colony morphology of strain(the odd number showed the front and the even number showed the back )
![]()
图 4 基于不同序列以最大似然法构建的真菌系统发育树
▲为供试菌株
Figure 4. Phylogenetic tree of fungi constructed using the maximum likelihood method and based on different sequences
▲ indicates test strains
表 1 供试菌株信息
Table 1 Information of test strains
菌株
Strain采样地
Sampling place采样时间
Sampling timeSB003 湖北省神农架林区 Shennongjia forest area, Hubei Province 2007−07 SB004 甘肃省天水市麦草沟自然保护区
Maicaogou nature reserve, Tianshui, Gansu Province2007−08 SB006 甘肃省临夏回族自治州莲花山自然保护区
Lianhuashan nature reserve, Linxia Hui Autonomous Region, Gansu Province2007−08 SB009 四川省成都市青羊区文殊院 Wenshuyuan, Qingyang district, Chengdu, Sichuan Province 2007−10 SB010 云南省昆明市 Kunming, Yunnan Province 2007−09 SB015 河南省鹤壁市 Hebi, Henan Province 2007−10 SB026 四川省雅安市宝兴县蜂桶寨
Fengtong village, Baoxing county, Ya’an, Sichuan Province2007−07 SB032 韩国济州岛天地渊瀑布 Tiandiyuan waterfall, Jeju Island, South Korea 2008−07 SB035 云南省红河州石屏县 Shiping county, Honghe Prefecture, Yunnan Province 2009−05 SB036 内蒙古自治区呼伦贝尔市新巴尔虎左旗
Xinbalhuzuoqi, Hulunbuir, Inner Mongolia Autonomous Region2009−08 SB037 内蒙古自治区兴安盟扎赉特旗
Zhalaitqi, Xing’an League, Inner Mongolia Autonomous Region2009−08 SB038 内蒙古自治区锡林郭勒盟正镶白旗
Zhengxiangbaiqi, Xilinguole League, Inner Mongolia Autonomous Region2009−08 SB039 内蒙古自治区阿拉善盟贺兰山国家级自然保护区
Helan Mountain National Nature Reserve, Alxa League, Inner Mongolia Autonomous Region2009−08 SB041 内蒙古自治区锡林郭勒盟正镶白旗
Zhengxiangbaiqi, Xilinguole League, Inner Mongolia Autonomous Region2009−08 SB043 西藏自治区林芝市波密县古乡
Gu village, Bomi county, Linzhi, Tibet Autonomous Region2009−09 SB050 广东省广州市白云山 Baiyun Mountain, Guangzhou, Guangdong Province 2009−10 SB051 广东省河源市龙川县 Longchuan county, Heyuan, Guangdong Province 2010−05 SB057 贵州省雷山县方祥乡 Fangxiang village, Leishan county, Guizhou Province 2010−07 SB062 内蒙古自治区阿拉善盟古拉本敖包镇
Gulaben Aobao Town, Alxa League, Inner Mongolia Autonomous Region2010−08 SB063 内蒙古自治区阿拉善盟贺兰山脉中部哈拉乌
Harau, central Helan Mountain, Inner Mongolia Autonomous Region2010−08 SB671 云南省镇沅彝族哈尼族拉祜族自治县
Zhenyuan Yi Hani Lahu Autonomous county, Yunnan Province2009−10 SB672 福建省厦门市鼓浪屿 Gulangyu, Xiamen, Fujian Province 2010−09 SB673 云南省迪庆州维西县保和镇 Baohe town, Weixi county, Diqing Prefecture, Yunnan Province 2010−09 SB674 广东省广州市华南农业大学林学与风景园林学院
College of forestry and landscape architecture, South China Agricultural University, Guangzhou,
Guangdong Province2010−10 SP016 湖北省神农架林区老君山 Laojun Mountain, Shennongjia forest area, Hubei Province 2007−08 SP031 甘肃省张掖市肃南县东柳沟 Dongliugou, Sunan county, Zhangye, Gansu Province 2007−08 SP433 内蒙古自治区呼伦贝尔市牙克石市
Yakeshi, Hulunbuir, Inner Mongolia Autonomous Region2009−08 SP665 云南省楚雄彝族自治州峨绿公园西灵宫
Xiling Palace, Erlu Park, Chuxiong Yi Autonomous Prefecture, Yunnan Province2009−09 SP670 湖北省宜昌市蛤蟆泉 Toad Spring, Yichang, Hubei Province 2009−10 
下载: 导出CSV
表 2 虫生真菌DNA的PCR扩增引物
Table 2 PCR amplification primers for entomogenous fungal DNA
基因序列 Gene sequence 引物 Primer 引物序列(5′→3′) Primer sequence ITS ITS4F TCCTCCGCTTATTGATATGC ITS5R GGAAGTAAAAGTCGTAACAAGG Bloc B5.1F CGACCCGGCCAACTACTTTGA B3.1R GTCTTCCAGTACCACTACGCC TEF 983F GCYCCYGGHCAYCGTGAYTTYAT 2218R ATGACACCRACRGCRACRGTYTG
下载: 导出CSV
表 3 虫生真菌的产孢量1)
Table 3 Spore productions of entomogenous fungi
菌株
Strain产孢量/(×107 mL−1)
Spore production菌株
Strain产孢量/(×107 mL−1)
Spore production菌株
Strain产孢量/(×107 mL−1)
Spore productionSB003 3.25±0.59f SB037 9.86±0.69abcdef SB671 4.64±0.71ef SB004 14.27±2.63abc SB038 14.78±2.06ab SB672 8.35±1.72abcdef SB006 13.65±1.59abcd SB039 16.38±2.83a SB673 5.15±1.10def SB009 9.42±1.19abcdef SB041 3.48±0.47f SB674 6.61±1.75bcdef SB010 6.72±1.92bcdef SB043 11.16±1.48abcdef SP016 4.04±0.49ef SB015 8.29±0.74abcdef SB050 12.64±3.07abcde SP031 9.32±1.73abcdef SB026 9.19±1.73abcdef SB051 6.54±0.47bcdef SP433 12.54±2.36abcde SB032 16.56±1.59a SB057 9.30±0.37abcdef SP665 5.57±0.89cdef SB035 3.84±0.56ef SB062 6.02±1.76bcdef SP670 7.26±0.49bcdef SB036 11.23±2.66abcdef SB063 7.88±1.59abcdef 1)数据后不同小写字母者表示不同菌株间菌落产孢量差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters indicate significant differences in spore production among different strains (P<0.05, Duncan’s test)
下载: 导出CSV
表 4 供试菌株ITS区、TEF区和Bloc区的序列比对结果
Table 4 Results of comparing the ITS regions, TEF regions and Bloc regions of test strains
真菌
Fungus菌株登录号 Accession number of strain ITS区
ITS regionTEF区
TEF regionBloc区
Bloc region球孢白僵菌 Beauveria bassiana KX901309 JQ867178 JQ867116 球孢白僵菌 B. bassiana AY531976 JQ867172 JQ867106 球孢白僵菌 B. bassiana KX901310 JQ867171 JQ867109 球孢白僵菌 B. bassiana GU233705 JQ867157 JQ867118 球孢白僵菌 B. bassiana GU233701 JQ867151 JQ867125 球孢白僵菌 B. bassiana AY883691 AY883788 球孢白僵菌 B. bassiana KU725693 KU725706 假球孢白僵菌 Beauveria pseudobassiana MH259852 HQ880999 HQ880727 假球孢白僵菌 B. pseudobassiana LC163567 LC163595 环链棒束孢 Isaria cateniannulata GU734760 GU734758
下载: 导出CSV
表 5 不同浓度的孢子悬浮液处理烟粉虱2龄若虫7 d的死亡率1）
Table 5 Mortality of the 2nd instar nymphs of Bemisia tabaci infected with spore suspensions of different concentrations for 7 days
% 菌株 Strain 1×107 mL−1 1×108 mL−1 菌株 Strain 1×107 mL−1 1×108 mL−1 CK 9.43±0.58d 11.73±2.00d SB043 47.67±4.00abc 80.33±3.93ab SB003 56.43±3.45abc 60.06±2.68bc SB050 49.47±1.47abc 81.30±4.86ab SB004 55.23±5.85abc 75.23±5.49abc SB051 59.67±4.17abc 70.87±5.90abc SB006 51.03±0.54abc 72.33±6.90abc SB057 51.67±3.55abc 76.00±2.34abc SB009 67.93±1.68ab 82.93±2.88ab SB062 35.81±3.54c 63.34±3.93abc SB010 47.90±4.42abc 70.20±4.94abc SB063 68.03±5.46ab 71.40±2.03abc SB015 64.30±7.27ab 77.07±4.95abc SB671 59.80±1.12abc 59.50±9.22bc SB026 51.77±4.15abc 67.80±3.0abc SB672 44.67±6.17bc 62.43±7.28bc SB032 56.57±2.95abc 77.69±3.88ab SB673 53.23±5.62abc 70.53±2.57abc SB035 53.07±2.37abc 59.99±2.77bc SB674 67.20±7.04ab 77.20±2.70abc SB036 55.33±6.36abc 80.67±4.41ab SP016 54.47±8.25abc 72.07±4.82abc SB037 50.17±7.69abc 74.00±3.21abc SP031 49.43±2.13abc 52.85±1.00c SB038 48.03±6.08abc 70.87±5.90abc SP433 69.03±3.76ab 87.37±0.43a SB039 72.77±2.14a 73.13±5.20abc SP665 49.70±4.65abc 64.90±3.26abc SB041 50.33±4.36abc 70.33±5.42abc SP670 65.77±0.50ab 72.40±1.56abc 1)相同浓度数据后不同小写字母者表示不同菌株间致病力差异显著(P<0.05，Duncan’s法)
1)Different lowercase letters after data at the same concentration indicate significant differences in pathogenicity among different strains (P<0.05, Duncan’s test)
下载: 导出CSV
-
[1] 查岭生, 文庭池, KEVIN D H, et al. 中国虫霉目真菌及其寄主昆虫名录更新[J]. 菌物学报, 2016, 35(6): 666-683. [2] 阙生全, 喻爱林, 刘亚军, 等. 白僵菌应用研究进展[J]. 中国森林病虫, 2019, 38(2): 29-35. doi: 10.3969/j.issn.1671-0886.2019.02.007 [3] 蒲蛰龙, 李增智. 昆虫真菌学[M]. 合肥: 安徽科学技术出版社, 1996: 402-715. [4] WRIGHT M S, LAX A R. Improved mortality of the Formosan subterranean termite by fungi, when amended with cuticle-degrading enzymes or eicosanoid biosynthesis inhibitors[J]. Folia Microbiol, 2016, 61(1): 73-83. doi: 10.1007/s12223-015-0412-0
[5] HEINIG R L, PAAIJMANS K P, HANCOCK P A, et al. The potential for fungal biopesticides to reduce malaria transmission under diverse environmental conditions[J]. J Appl Ecol, 2015, 52(6): 1558-1566. doi: 10.1111/1365-2664.12522
[6] 刘晓晓, 王星冉, 王义勋, 等. 马尾松毛虫高毒力球孢白僵菌菌株的生物学特性[J]. 南方农业学报, 2017, 48(6): 1019-1023. doi: 10.3969/j.issn.2095-1191.2017.06.13 [7] 辛正, 张晓, 刘慧媛, 等. 白僵菌在蟑螂防治中的研究进展[J]. 中华卫生杀虫药械, 2017, 23(5): 401-406. [8] SUN M, REN Q, GUAN G, et al. Eectiveness of Beauveria bassiana sensu lato strains for biological control against Rhipicephalus (Boophilus)microplus (Acari: Ixodidae) in China[J]. Parasitol Int, 2013, 62(5): 412-415. doi: 10.1016/j.parint.2013.04.008
[9] CHEN W, HASEGAWA D D, KAUR N, et al. The draft genome of whitefly Bemisia tabaci MEAM1, a global crop pest, provides novel insights into virus transmission, host adaptation, and insecticide resistance[J]. BMC Biol, 2016, 14(1): 110. doi: 10.1186/s12915-016-0321-y
[10] CHEN M J, HUANG B, LI Z Z, et al. Morphological and genetic characterisation of Beauveria sinensissp. nov. from China[J]. Mycotaxon, 2013, 124(1): 301-308. doi: 10.5248/124.301
[11] KESSIA F C, PANTOJA, KELLY C G, et al. Identification of Capsicum accessions to lerantto Tomato severe rugose virus and resistant to Bemisia tabaci Middle East-Asia Minor1(MEAM1)[J]. Trop Plant Pathol, 2018, 43(2): 138-145. doi: 10.1007/s40858-018-0212-6
[12] 李英梅, 杨苗苗, 刘晨, 等. 杀虫剂对烟粉虱传番茄黄化曲叶病毒病的防控效果[J]. 西北农业学报, 2019, 28(3): 466-474. [13] JURAJ M, JAROSLAV M, JANKA M, et al. Phylogenetic structure and habitat associations of Beauveria species isolated from soils in Slovakia[J]. J Invertebr Pathol, 2016, 140: 46-50. doi: 10.1016/j.jip.2016.08.009
[14] 潘李隆, 迟瑶, 范运运, 等. 烟粉虱对菜豆金黄花叶病毒属病毒的传播:病毒对媒介昆虫的适应[J]. 昆虫学报, 2019, 62(7): 885-894. [15] GILBERTSON R L, BATUMAN O, WEBSTER C G, et al. Role of the insect supervectors Bemisia tabaci and Frankliniella occidentalisin the emergence and global spread of plant viruses[J]. Annu Rev Virol, 2015, 2(1): 67-93. doi: 10.1146/annurev-virology-031413-085410
[16] ABDESSAMAD I, MUZAMMIL H, PAUL M, et al. Entomopathogenic fungus Beauveria: Host specificity, ecology and significance of morpho-molecular characterization in accurate taxonomic classification[J]. J Asia-Pac Entomol, 2017, 20(4): 1204-1212. doi: 10.1016/j.aspen.2017.08.015
[17] 褚栋, 张友军. 近10年我国烟粉虱发生为害及防治研究进展[J]. 植物保护, 2018, 44(5): 51-55. [18] 郭东升, 翟颖妍, 任广伟, 等. 白僵菌属分类研究进展[J]. 西北农业学报, 2019, 28(4): 497-509. doi: 10.7606/j.issn.1004-1389.2019.04.001 [19] 余垚颖, 郭应菊, 杨雪, 等. 我国虫生真菌的研究现状[J]. 四川农业科技, 2018(12): 27-29. doi: 10.3969/j.issn.1004-1028.2018.12.015 [20] REHNER S A, BUCKLEY E P. A Beauveria phylogeny inferred from nuclear ITS and EF1-α sequences: Evidence for cryptic diversification and links to Cordyceps teleomorphs[J]. Mycologia, 2005, 97: 84-98.
[21] REHNER S A, MINNIS A M, SUNG G H, et al. Phylogeny and systematics of the anamorphic, entomopathogenic genus Beauveria[J]. Mycologia, 2011, 103: 1055-1073. doi: 10.3852/10-302
[22] MEYLING N V, PILZ C, KELLER S, et al. Diversity of Beauveria spp. isolates from pollen beetles Meligethes aeneusin Switzerland[J]. J Invertebr Pathol, 2012, 109: 76-82. doi: 10.1016/j.jip.2011.10.001
[23] 张云月, 张军, 李丽莉, 等. 球孢白僵菌有性世代培养及子代毒力稳定性鉴定[J]. 中国生物防治学报, 2017, 33(3): 364-370. [24] 芦俊佳, 夏举飞, 庄飞, 等. 楚雄腮扁叶蜂虫生真菌粉红粘帚霉的鉴定及其生物学特性[J]. 植物保护学报, 2018, 45(3): 622-631. [25] 杨帆, 曹涵铭, 陈园园, 等. AM真菌结构及鉴定方法的研究进展[J]. 现代农业科技, 2019(16): 152-154. doi: 10.3969/j.issn.1007-5739.2019.16.095 [26] 张伟铮, 关文苑, 李松, 等. ITS序列分析与MALDI-TOF MS质谱技术在丝状真菌鉴定中的应用[J]. 菌物学报, 2019, 38(8): 1298-1305. [27] 甘森宁, 陆梅颖, 肖志邦, 等. 植物内生真菌的分类鉴定方法[J]. 安徽农学通报, 2019, 25(10): 36-37. doi: 10.3969/j.issn.1007-7731.2019.10.014 [28] 李美君, 白庆荣, 臧连生, 等. 水稻二化螟病原真菌鉴定及其致病力[J]. 中国生物防治学报, 2019, 35(1): 63-69. [29] 方志刚, 张立钦, 林新春. 不同球孢白僵菌菌株DNA的RAPD分析[J]. 南京林业大学学报(自然科学版), 2001, 25(4): 65-68. doi: 10.3969/j.issn.1000-2006.2001.04.015 [30] 王慧. 绿僵菌和白僵菌不同菌株对烟粉虱若虫的致病力及应用研究[D]. 合肥: 安徽农业大学, 2013. [31] 吴乾兴, 刘勇, 张友军, 等. 烟粉虱常用药剂抗性监测研究[J]. 现代农业科技, 2018(2): 119-120. [32] 王登杰, 王海鸿, 张桃, 等. 三株球孢白僵菌侵染烟粉虱的比较生长动力学及其毒力[J]. 菌物学报, 2016, 35(5): 577-585. [33] 银航, 窦雪绒, 张云霞, 等. 球孢白僵菌、保幼激素联用防治花椒蚜虫[J]. 陕西农业科学, 2019, 65(11): 57-61. doi: 10.3969/j.issn.0488-5368.2019.11.018 [34] 张爽, 杨佳利, 马军, 等. 苏云金杆菌(Bacillus thuringiensis)与球孢白僵菌(Beauveria bassiana)复配对小菜蛾幼虫(Plutella xylostella)的影响研究[J]. 生物灾害科学, 2017, 40(4): 288-291. doi: 10.3969/j.issn.2095-3704.2017.04.63 [35] 杨芷, 路杨, 毛刚, 等. 赤眼蜂携带球孢白僵菌防治亚洲玉米螟技术研究与应用[J/OL]. 中国生物防治学报, 2019. doi: 10.16409/j.cnki.2095-039x.2019.05.011. -
期刊类型引用(4)
1. 詹芳芳，林俨，傅世辉，林华英，林孝春，蔡守平. 侵染松墨天牛幼虫的虫生真菌鉴定及致病性测定. 中国森林病虫. 2024(05): 21-26 . 
百度学术
2. 王露露，王辉，熊焰，周霞，王军，王健华，伍苏然. 虫生真菌防治农作物害虫的研究进展. 热带生物学报. 2022(03): 309-314 . 百度学术
3. 古开平，吴伟，胡琼波. 几种土壤真菌对烟粉虱的生物活性研究. 环境昆虫学报. 2022(04): 1010-1018 . 百度学术
4. 张宣，刘雪莹，胡迪，袁向群，冯克云，康波，张军，罗进仓，南宏宇，李怡萍. 河西走廊棉田主要害虫发生动态及植物源杀虫剂防治技术. 植物保护学报. 2021(05): 1125-1138 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 1056
- HTML全文浏览量: 4
- PDF下载量: 1316
- 被引次数: 9
