廣西家蠶僵病流行趨勢及其主要病原遺傳規(guī)律研究

唐亮，胡文娟，蔣滿貴，黃深惠，董桂清，李標，唐名艷，趙燁蕓，何驥，潘志新

（廣西蠶業(yè)技術(shù)推廣站/廣西蠶業(yè)科學研究院，廣西南寧530007）

0 引言

【研究意義】蠶病是在家蠶飼養(yǎng)環(huán)節(jié)中制約蠶繭產(chǎn)量、質(zhì)量及蠶農(nóng)經(jīng)濟效益的主要因素之一（祁廣軍，2013）。其中，僵病是養(yǎng)蠶過程中全齡易發(fā)的真菌蠶病，且病源多樣，常經(jīng)風媒傳播或直接接觸蠶體而感染致病，若防控不及時極易暴發(fā)流行，給蠶桑產(chǎn)業(yè)造成重大經(jīng)濟損失（廖華珍等，2009；唐亮等，2013）。因此，掌握家蠶僵病的發(fā)病規(guī)律及病原特性，建立科學的防控策略和措施，對有效降低家蠶僵病危害、保障蠶業(yè)穩(wěn)定發(fā)展及促進蠶農(nóng)增收均具有重要意義?！厩叭搜芯窟M展】牟志美等（1999）研究報道，常規(guī)年份山東家蠶的僵病危害在5%～8%，與飼養(yǎng)品種、發(fā)育時期及野外昆蟲的發(fā)生密切相關；朱方容（2008）、廖華珍等（2009）調(diào)查發(fā)現(xiàn)廣西2008—2009年連續(xù)2年春季家蠶僵病暴發(fā)流行，造成當?shù)匦Q繭產(chǎn)量銳減，其主要病原可能來源于飛機噴施森林防蟲用的白僵菌。引起家蠶僵病的病原除白僵菌屬外，還有擬青霉屬（朱方容等，1997；苘娜娜等，2015）、曲霉屬、萊氏野村菌屬（王長軍，2000）及蠟蚧屬（苘娜娜等，2015）等多種真菌病原，存在持續(xù)的感染能力（余金勇等，2007），但不同菌株間對宿主的毒力存在明顯差異（Saito and Brownbridge，2016）。Wang等（2005）利用簡單序列重復（ISSR）分子標記研究39個白僵菌屬分離株的遺傳多樣性，結(jié)果表明ISSR分子標記能成功將3個密切相關的物種與球孢白僵菌（）區(qū)分開來，并證實球孢白僵菌分離株與其地理來源存在一定關聯(lián)，但這些分離株與其昆蟲宿主間無明確的相關性。呂思行等（2012）采用微衛(wèi)星分子標記對廣東和廣西蠶區(qū)的白僵菌居群之間和居群之內(nèi)的遺傳多態(tài)性進行分析，結(jié)果發(fā)現(xiàn)兩廣白僵菌居群之間的基因分化系數(shù)（）為0.0590，多態(tài)位點百分率（PPL）為97.73%，Nei's基因多樣性指數(shù)（）為0.1896，Shannon信息多樣性指數(shù)（）為0.3165，表明兩廣家蠶來源的白僵菌居群間遺傳分化較小。唐亮等（2013）通過研究宿主濕度環(huán)境對球孢白僵菌致病力的影響，發(fā)現(xiàn)球孢白僵菌在病蠶體中宿存60 d內(nèi)對家蠶的致死率仍然很高，因此及時清除病蠶尸體并進行尸體病原性滅活，可降低白僵菌的二次感染。陳雪（2014）通過ISSR遺傳多樣性數(shù)據(jù)及菌株對家蠶的致死數(shù)據(jù)研究全國各地球孢白僵菌種群關系，結(jié)果表明家蠶白僵病的發(fā)生與松毛蟲白僵病的自然流行或人工釋放球孢白僵菌殺蟲劑均無直接關系。王婷婷（2019）在柞蠶僵病病原分子鑒定的研究中發(fā)現(xiàn)，野外柞蠶僵蠶、僵蛹及室內(nèi)人工接種僵蠶的內(nèi)轉(zhuǎn)錄間隔區(qū)（ITS）序列相似性存在明顯差異，與已知基因序列（KM205065.1）的相似性在96.01%～100.00%。王露露等（2022）從山東和廣西僵蠶及生防菌中篩選獲得對家蠶具有高毒力的白僵菌，并優(yōu)化其發(fā)酵條件，為制備僵蠶提供了優(yōu)質(zhì)的菌種資源。Wang等（2022）通過對比、、、、及共6個基因序列，以評估中國、泰國15株白僵菌的系統(tǒng)發(fā)育位置，結(jié)果篩選出3株具有生防潛力的白僵菌株?！颈狙芯壳腥朦c】查明廣西家蠶僵病發(fā)生狀況、流行規(guī)律，掌握病原特性及其遺傳規(guī)律等是有效防控該病的基礎，但目前尚缺乏對廣西家蠶僵病的系統(tǒng)性研究。【擬解決的關鍵問題】通過分析2008—2021年廣西主要蠶區(qū)蠶病調(diào)查數(shù)據(jù)，結(jié)合病原菌形態(tài)學與分子生物學研究，旨在掌握廣西主要蠶區(qū)的僵病流行規(guī)律及其病原特性，為精準防控家蠶僵病提供科學依據(jù)。

1 材料與方法

1.1 蠶病數(shù)據(jù)

廣西家蠶僵病發(fā)病原始數(shù)據(jù)由廣西蠶桑病蟲害監(jiān)測網(wǎng)的常規(guī)蠶病監(jiān)測數(shù)據(jù)和相關科技項目蠶病調(diào)查數(shù)據(jù)組成；數(shù)據(jù)匯總年份為2008—2021年，月份為每年3—10月。每批次全齡期蠶病調(diào)查數(shù)據(jù)為三齡起蠶、五齡起蠶、上蔟后第2 d的僵病發(fā)生率總和。

1.2 供試病原

家蠶僵病病原采自桂中、桂南、桂西、桂西北、桂東南等廣西主要桑蠶主產(chǎn)區(qū)的僵蠶，清潔處理后接種于薩氏培養(yǎng)基上，劃線分離純化獲得單菌落，經(jīng)顯微形態(tài)鑒定和產(chǎn)孢測定后，保存于廣西蠶業(yè)技術(shù)推廣站蠶桑病蟲害研究室菌種庫。

1.3 供試菌絲及分生孢子制備

將供試菌株接種至薩氏培養(yǎng)基，培養(yǎng)箱（25±1）℃倒置培養(yǎng)10 d左右，待大量產(chǎn)生菌絲和分生孢子后，刮取表層菌絲及分生孢子，經(jīng)液氮速凍后，分別置于-20和-80℃冰箱保存?zhèn)溆谩?/p>

1.4 病原菌對家蠶致病力測試

采用濾紙接觸法（唐亮等，2013）添毒兩廣二號的二齡起蠶，置于（25±1）℃、相對濕度在90%以上的環(huán)境下飼養(yǎng)，連續(xù)觀察7 d內(nèi)家蠶發(fā)病情況。重復3次。

1.5 病原菌DNA提取

參照朱衡等（1994）的方法提取病原菌總DNA，以1.5%瓊脂糖凝膠電泳檢測DNA質(zhì)量；采用核酸測定儀測定DNA濃度后稀釋至20 ng/μL，-20℃保存?zhèn)溆谩?/p>

1.6 病原菌ITS序列分析

采用真菌通用引物ITS4（5'-TCCTCCGCTTAT TGATATGCITS-3'）/ITS5（5'-GGAAGTAAAAGTCG TAACAAGG-3'）組合進行PCR擴增（陳立杰等，2008），PCR反應體系50.0μL：DNA模板2.0μL，dNTP Mixture（2.5 mmol/L）2.5μL，ITS4/ITS5引物（20μmol/L）各1.5μL，10×ExBuffer（MgPlus）5.0μL，Ex酶（5 U/μL）0.2μL，ddHO補足 至50.0μL。擴增程序：94℃預變性3 min；94℃1 min，55℃1 min，72℃2 min，進行30個循環(huán)；72℃延伸5 min。PCR擴增產(chǎn)物采用1.5%瓊脂糖凝膠電泳進行檢測，回收目的條帶，經(jīng)PCR Purification Kit純化后以擴增引物進行雙向測序，所得測序結(jié)果拼接后輸入GenBank，通過BLAST進行比對分析。

1.7 病原菌ISSR多態(tài)性分析

從測試樣本中選取4個樣本對20對引物進行測試，篩選出10個擴增條帶清晰、多態(tài)性好的引物用于ISSR-PCR反應。反應體系25.0μL：DNA模板約40 ng，10×Bufer 2.5μL，DNA聚合酶1.25 U，引物（10 ng/μL）1.0μL，dNTPs（2.5 mmol/L）0.4μL，MgCl（1.5 mmol/L）2.0μL，ddHO補足至25.0μL。擴增程序：94℃預變性5 min；94℃45 s，55～62℃45 s，72℃1.5 min，進行35個循環(huán)；72℃延伸10 min，4℃保存。PCR擴增產(chǎn)物以2.0%瓊脂糖凝膠電泳進行檢測。

1.8 統(tǒng)計分析

基于ITS4/ITS5引物進行PCR擴增，并將雙向測序結(jié)果輸入GenBank進行BLAST比對分析。ISSRPCR擴增產(chǎn)物在瓊脂糖凝膠電泳的某個相同遷移率位置上有DNA條帶記為1，無DNA條帶則記為0，形成ISSR表型數(shù)據(jù)矩陣，然后使用NTsys計算遺傳相似系數(shù)，并基于遺傳相似系數(shù)進行UPGMA聚類分析。同時，使用PopGene32計算供試材料的遺傳多樣性指數(shù)（安丹丹等，2021；周春寶等，2021），并以SPSS 25.0進行單因素方差分析（One-way ANOVA）。

2 結(jié)果與分析

2.1 不同年份家蠶僵病的發(fā)生規(guī)律

2008—2021年廣西主要蠶區(qū)的家蠶僵病平均發(fā)病情況如圖1所示。廣西家蠶僵病發(fā)病率整體上呈波動下降趨勢，連續(xù)14年的平均發(fā)病率為5.19%。其中，2008年的發(fā)病率為10.08%，顯著高于其他年份（＜0.05，下同）；發(fā)病率最低的年份是2019年，發(fā)病率僅為3.16%。此外，廣西家蠶僵病年均發(fā)病率曲線在2012年呈現(xiàn)一個轉(zhuǎn)折點，2008—2012年整體上呈快速下降趨勢，2012—2021年則現(xiàn)出上下波動，但總體維持在較低水平，年均發(fā)病率僅為4.08%。

2.2 不同月份家蠶僵病的發(fā)生規(guī)律

以月為單元進行分析，結(jié)果顯示廣西家蠶僵病以5月的發(fā)病率最高（圖2），達13.37%，與其他月份的差異均達顯著水平；次高峰出現(xiàn)在4月和10月，對應的發(fā)病率分別為7.91%和6.87%，二者差異不顯著（＞0.05，下同）。廣西家蠶僵病發(fā)病高峰和次高峰時段與廣西溫潤多雨氣候存在高度重疊，推測溫濕度對家蠶僵病的發(fā)生有較大影響。

2.3 不同蠶區(qū)不同地形家蠶僵病的發(fā)生規(guī)律

據(jù)2008—2021年桂中、桂南、桂西、桂西北和桂東南五大蠶區(qū)的調(diào)查數(shù)據(jù)顯示，各蠶區(qū)平緩丘陵地形、山區(qū)地形的家蠶僵病發(fā)生情況如圖3所示。由圖3可知，山區(qū)地形的家蠶僵病發(fā)生率普遍高于平緩丘陵地區(qū)。平緩丘陵蠶區(qū)的家蠶僵病發(fā)病率排序為桂中（13.55%）＞桂西北（12.66%）＞桂南（11.56%）＞桂西（11.05%）＞桂東南（9.53%），可能與桂中、桂西北在春季低溫陰雨天氣偏多，濕度偏高有利于真菌萌發(fā)有關。值得注意的是，5月是各蠶區(qū)家蠶僵病的發(fā)病高峰期，但進入6月中下旬家蠶僵病發(fā)生率迅速下降，可能是平緩丘陵地區(qū)氣溫上升迅速、濕度快速下降，超出了真菌發(fā)育的最適溫濕度而表現(xiàn)出家蠶僵病發(fā)病率急劇下降。

山地蠶區(qū)的家蠶僵病發(fā)病率排序為桂中（18.91%）＞桂西北（17.96%）＞桂西（16.01%）＞桂西（14.90%）＞桂東南（11.61%），可能是桂中、桂西北山區(qū)養(yǎng)蠶多位于河谷地區(qū)，相對日照較短，下午16：00后的太陽照射盲區(qū)及夜間會出現(xiàn)局部相對濕度偏高而有利于真菌萌發(fā)和感染傳播。相對于平緩丘陵地區(qū)5月家蠶僵病發(fā)病高峰過后次月發(fā)病率驟降的情況，山地蠶區(qū)的家蠶僵病發(fā)病率下降趨勢略顯緩慢，尤其是桂中、桂西、桂西北地區(qū)6月的家蠶僵病發(fā)生率仍維持在7.80%～10.28%，提示山地蠶區(qū)的家蠶僵病危害時間更嚴重且更持久。

2.4 病原菌形態(tài)學觀察及其對家蠶毒力測試結(jié)果

本研究收集了桂中、桂南、桂西、桂西北及桂東南5個蠶桑主產(chǎn)區(qū)多年發(fā)生的僵蠶、僵蟲，以及部分白僵菌（農(nóng)藥），以此為基礎分離獲得60株病原單菌落（表1）用于后續(xù)測試。分離獲得的單菌落接種至薩氏培養(yǎng)基上，26℃培養(yǎng)14 d，連續(xù)觀察其生長情況并統(tǒng)計菌落的產(chǎn)孢情況，綜合評價病原菌的危害。然后挑取菌絲和分生孢子制片，在顯微鏡下觀察分生孢子形狀、分生孢子梗分支情況及產(chǎn)孢方式等真菌分類特征，初步鑒定分離菌株即為白僵菌（圖4）。

篩選廣西各蠶區(qū)及不同來源農(nóng)藥白僵菌株產(chǎn)孢量第一或前兩名的菌株，統(tǒng)計分析其中的10株菌株（表2），產(chǎn)孢量排名前3的菌株是崇左白僵菌株（農(nóng)藥）、愛山白僵菌株（農(nóng)藥）和桂南蠶區(qū)白僵菌株，對應的平均產(chǎn)孢量分別為3.4×10、2.8×10和2.6×10分生孢子/mm，其中前兩者是森林防蟲用農(nóng)藥白僵菌，而桂南蠶區(qū)白僵菌是從白僵蠶分離獲得，說明農(nóng)藥源白僵菌在產(chǎn)孢量方面具有明顯優(yōu)勢。其中，桂南蠶區(qū)白僵菌株來源于廣西蠶業(yè)技術(shù)推廣站品種保育過程的僵蠶，由于廣西蠶業(yè)技術(shù)推廣站地處南寧市區(qū)，養(yǎng)蠶環(huán)境獨立、飼養(yǎng)技術(shù)規(guī)范，與其他蠶區(qū)病原基本無交集，故推測該菌株是1株相對獨立的野生白僵菌，但其產(chǎn)孢量在所有測試菌株中最低。

同時測試上述10株白僵菌株對家蠶的致病力，結(jié)果（圖5）顯示，這些白僵菌株對兩廣二號二齡幼蟲的致死速率存在一定差異，但7 d累計致死率均達100.00%，致病力最強的為崇左白僵菌株（農(nóng)藥），其攻毒48 h后即出現(xiàn)家蠶吐水現(xiàn)象，72 h累計致死率達100.00%。72 h累計致死率排名第一梯隊為3種農(nóng)藥源白僵菌株：崇左白僵菌株（農(nóng)藥）（cy34）100.00%、南寧白僵菌株（農(nóng)藥）（cy33）98.89%、愛山白僵菌株（農(nóng)藥）（cy49）97.78%；其次是蠶區(qū)病蠶分離的白僵菌株和松林僵蟲白僵菌株：桂中松林白僵菌株（cy13）94.44%、桂南2白僵菌株（cy28）93.33%、桂中2白僵菌株（cy16）93.33%、桂西北白僵菌株（cy51）92.22%、桂中1白僵菌株（cy24）91.11%、桂東南白僵菌株（cy22）83.33%，即蠶區(qū)病蠶白僵菌株致死率在83.33%～93.33%，低于農(nóng)藥源白僵菌株對家蠶的致死率；最后為獨立的桂南1白僵菌株（cy20）8.89%，明顯低于農(nóng)藥用白僵菌株。經(jīng)顯著差異分析發(fā)現(xiàn)，除愛山白僵菌株（農(nóng)藥）與桂中松林白僵菌株、桂南2白僵菌株、桂中1白僵菌株、桂中2白僵菌株及桂西北白僵菌株的差異不顯著外，南寧白僵菌株（農(nóng)藥）、崇左白僵菌株（農(nóng)藥）均顯著高于蠶區(qū)白僵菌株的致病力，松林白僵菌株顯著高于桂東南蠶區(qū)和桂南獨立蠶區(qū)的白僵菌株。即10株白僵菌株對家蠶的致病力總體上呈現(xiàn)出農(nóng)藥源白僵菌株＞松林白僵菌株＞蠶區(qū)病蠶白僵菌株的趨勢，同時存在桂南、桂中、桂西北蠶區(qū)白僵菌株感染力高于桂東南、桂西蠶區(qū)的白僵菌株?？梢?，桂南、桂中、桂西北蠶區(qū)受白僵菌農(nóng)藥影響的可能性較大，即蠶桑主產(chǎn)區(qū)應慎用白僵菌農(nóng)藥。

2.5 病原菌ITS序列分析結(jié)果

以最初凍存的病原菌絲和分生孢子為樣本，提取60株病原單菌落總DNA，經(jīng)1.5%瓊脂糖凝膠電泳檢測合格（圖4-A）后，利用ITS4/ITS5引物進行目的片段擴增，PCR擴增電泳結(jié)果（圖4-B）顯示60株病原樣本均能擴增獲得預期大小的目的條帶，可用于后續(xù)的ITS序列分析。

將PCR擴增產(chǎn)物送至生工生物工程（上海）股份有限公司進行雙向測序，測序結(jié)果輸入GenBank進行BLAST比對分析，結(jié)果顯示擴增的ITS序列存在高度保守性，有多個序列比對結(jié)果完全一致。在60株病原菌中有56條序列匹配到8個已知序列，序列相似性在97.83%～99.83%（表3）。其中，來自桂中蠶區(qū)2019年的3個樣本和桂南蠶區(qū)2019年的1個樣本，經(jīng)2次重復提取DNA均未能成功測序，故不計入后續(xù)分析。在56株匹配的菌株中，有45株為蠶區(qū)采集病原株，且其中39株為球孢白僵菌（）。可見，導致蠶區(qū)家蠶白僵病的主要病原為球孢白僵菌，占比達86.67%；誘發(fā)家蠶僵病的病原還有曲霉屬、擬青霉屬和蟲草菌屬等，表明蠶區(qū)家蠶僵病病原復雜多樣。

2.6 白僵菌遺傳多樣性及聚類分析結(jié)果

在ITS序列分析的基礎上，選擇45株蠶區(qū)白僵菌株樣本、4株農(nóng)藥源白僵菌株樣本及1株松林白僵菌株樣本共計50株菌株進行多態(tài)性分析。首先選用cy03、cy14、cy40、cy60等4個樣本進行ISSR引物篩選，結(jié)果從20條引物篩選出10條擴增條帶清晰、條帶數(shù)較多、多態(tài)性良好，且易于擴增的引物（表4）。利用篩選出的10條ISSR引物對50株菌株樣本的DNA進行ISSR-PCR擴增，擴增產(chǎn)物以2.0%瓊脂糖凝膠電泳進行檢測，并根據(jù)ISSR-PCR擴增結(jié)果以Quantity One繪制電泳模式圖，結(jié)果顯示，10條ISSR引物從50株菌株樣本中擴增出152條清晰條帶（表5），其中多態(tài)位點147個，多態(tài)百分數(shù)為87.5%～100.0%，平均為96.71%。各引物擴增獲得的條帶數(shù)在13～19條，平均為15.2條。從單條引物擴增得到的多態(tài)位點來看，以引物P8和引物P9擴增得到的多態(tài)位點最多，均為19個。10條引物擴增得到的多態(tài)位點平均為14.7個，表明這10條引物能很好地區(qū)分白僵菌株多樣性。利用PopGene32計算供試材料的遺傳多樣性指數(shù)，結(jié)果顯示，50株菌株樣本的等位基因數(shù)（）為1.8750～2.0000，有效等位基因數(shù)（）為1.2726～1.5625，Nei's基因多樣性指數(shù)（）為0.1740～0.3257，Shannon信息多樣性指數(shù)（）為0.2840～0.4880，均表現(xiàn)出極高的遺傳多樣性，尤其以引物TL-P8的擴增效果最優(yōu)（圖7）。

根據(jù)“0，1”矩陣數(shù)據(jù)，分別計算50個供試材料兩兩間的遺傳距離并構(gòu)建遺傳相似矩陣，結(jié)果顯示，cy13與cy16的遺傳相似系數(shù)（GS）最大，為0.9737，表明二者的遺傳距離非常近；而cy33（白僵屬）與cy25（擬青霉屬）和cy26（擬青霉屬）的GS最小，均為0.4276，表明其遺傳距離較遠。蠶區(qū)來源白僵菌株間的GS波動幅度在0.5921～0.9671，說明蠶區(qū)病原遺傳多樣性較高，同時中位區(qū)間在0.8000～0.8999，占56.41%，預示著超過半數(shù)的蠶區(qū)病原基因相關度較高，GS大于0.8000的比例高達78.00%，則提示區(qū)域內(nèi)的菌株變異較小或區(qū)域間存在頻繁交流。其中，蠶區(qū)病原的GS最高值（0.9737）出現(xiàn)在2013年的2株桂中蠶區(qū)病原菌株間（cy16與cy19），最小值（0.5921）則出現(xiàn)在2013年育種僵蠶菌株cy20與同年桂西北蠶區(qū)的僵蠶病原菌株cy23間，表明相同地域的白僵菌株會表現(xiàn)出更高的親緣關系，而跨區(qū)域菌株間的基因相關度較低。

根據(jù)50株菌株兩兩間的遺傳距離進行UPGMA聚類分析，結(jié)果（圖8）顯示，蠶區(qū)病原菌株與農(nóng)藥源白僵菌株聚類在相鄰位置，提示二者具有較緊密的遺傳關系。如cy28（桂南蠶區(qū)，2017年）與cy33［南寧白僵（農(nóng)藥），2017年］、cy51（桂西北蠶區(qū)，2011年）與cy49［愛山白僵（農(nóng)藥），2011年］、cy22（桂東南蠶區(qū)，2013年）與cy12［南寧白僵（農(nóng)藥），2013年］、cy24（桂中蠶區(qū)，2011年）與cy34［崇左白僵（農(nóng)藥），2011年］間的遺傳距離分別為0.9276、0.9342、0.9605和0.9568，預示著桂南、桂西北、桂東南、桂中等蠶區(qū)的病蠶白僵菌株與農(nóng)藥源白僵菌株存在非常近的親緣關系，進一步證實廣西各蠶區(qū)白僵病的發(fā)生流行與施用白僵菌農(nóng)藥存在密切關系。

3 討論

2008—2021年廣西主要蠶區(qū)的家蠶僵病平均發(fā)病率為5.19%，與牟志美等（1999）研究報道山東家蠶僵病常規(guī)年份的年均發(fā)病率在5%～8%相仿。廣西屬于高溫多濕的亞熱帶季風氣候，山東屬于溫帶季風氣候，兩地溫濕度差異明顯，但兩地的家蠶僵病發(fā)病率相近可能與廣西氣溫偏高，適宜白僵菌繁殖和浸染的時間范圍較窄有關。在2008和2009年廣西蠶區(qū)周邊環(huán)境引入大量森林防控白僵菌（朱方容，2008；廖華珍等，2009）的情形下仍取得良好防控效果，預示著在沒有大量外源病原菌影響的情況下僵病發(fā)病率可能會更低。自2009年廣西蠶區(qū)大面積引入外源性白僵菌后，2010和2011年廣西蠶區(qū)周邊未見大量施用白僵菌的報道，但2010年的家蠶僵病發(fā)病率顯著高于之后的多年平均發(fā)病率（4.08%），廣西家蠶僵病年均發(fā)病率曲線在2012年呈現(xiàn)一個轉(zhuǎn)折點。說明人為引入大量白僵菌株后，對當?shù)丶倚Q僵病的影響可能會維持2年左右，第3年則呈顯著下降趨勢，與余金勇等（2007）在研究白僵菌防松毛蟲持續(xù)致病力得出的結(jié)論相似。本研究通過對廣西主要蠶區(qū)及周邊地區(qū)收集的白僵菌株進行ITS序列分析，進一步明確廣西家蠶僵蠶的病原種類，包括白僵菌屬、曲霉屬、擬青霉屬和蟲草菌屬，表明廣西家蠶僵病病原呈多樣性。其中，擬青霉屬病原與朱方容等（1997）的研究結(jié)果相似。廣西家蠶僵病的主要病原仍是球孢白僵菌（占86.67%），說明廣西蠶區(qū)家蠶僵病的病原高度保守。

本研究對廣西五大蠶桑主產(chǎn)區(qū)家蠶僵病病原、4株生物防控白僵菌（農(nóng)藥）、1株松林白僵分離株和1株研究機構(gòu)育種白僵菌進行毒力測試，結(jié)果顯示生物防控白僵菌株＞松林白僵菌株＞主要蠶區(qū)白僵菌株＞研究機構(gòu)育種白僵菌。生物防控白僵菌株（農(nóng)藥）對家蠶的致病力顯著高于蠶區(qū)白僵菌株，提示其具備區(qū)域性危害的基礎條件，與石美寧和朱方容（2010）、劉吉平等（2011）的研究結(jié)果相似。本研究同時利用ISSR-PCR分析廣西主要蠶區(qū)白僵菌病原間的基因相關度，結(jié)果顯示同一地域相同年份的病原菌株具有非常近的遺傳距離，如cy16與cy19的GS最大（0.9737）；相同區(qū)域經(jīng)過多年后遺傳變異較大，如cy28與cy58的GS最?。?.7895）。廣西桑蠶主產(chǎn)區(qū)家蠶僵病病原的GS中位數(shù)區(qū)間在0.8000～0.9000，約有78.00%的GS大于0.8000，暗示廣西白僵菌居群的遺傳多樣性水平相對較低，與呂思行等（2012）的研究結(jié)果相似。此外，通過對比生物防治用白僵菌株（農(nóng)藥）與廣西部分蠶區(qū)白僵菌株，發(fā)現(xiàn)其遺傳距離非常近，同時蠶區(qū)周邊松林分離白僵株與前兩者存在極高的遺傳關聯(lián)度，預示著松林害蟲可能是白僵菌（農(nóng)藥）感染家蠶的中間宿主（陳雪，2014），即廣西蠶區(qū)白僵病的流行與施用白僵菌農(nóng)藥存在密切聯(lián)系。

4 結(jié)論

廣西家蠶僵病危害高峰期為5月，球孢白僵菌仍是最主要的病原菌，且ITS序列高度保守。不同來源白僵菌對家蠶致病力的排序為農(nóng)藥源白僵菌株＞松林白僵菌株＞蠶區(qū)僵蠶菌株，但三者相互間存在非常高的基因關聯(lián)性，提示廣西家蠶僵病的病原遺傳背景復雜，且蠶桑主產(chǎn)區(qū)應慎用白僵菌農(nóng)藥。