嗜線蟲致病桿菌YZ011對鱗翅目昆蟲殺蟲活性及比較基因組學分析

收稿日期：2024-08-01

基金項目：國家重點研發(fā)計劃項目（2024YFD1400900）；江蘇省農(nóng)業(yè)科技自主創(chuàng)新基金項目[CX（22）1009]；江蘇省國際合作項目（BZ2020039）；揚州市科技計劃項目（YZ2023244）；農(nóng)業(yè)基礎(chǔ)性長期性科技工作項目（NAES069AM04）

作者簡介：林曼曼（1990-），女，江蘇徐州人，博士，助理研究員，主要從事農(nóng)業(yè)害蟲生物防治。（E-mail）lin.man.happy@163.com

通訊作者：徐 健，（E-mail）bio-xj@163.com

摘要： 昆蟲病原線蟲共生菌是一類重要的生物防治資源。為明確本單位分離自小卷蛾斯氏線蟲共生菌YZ011對鱗翅目害蟲的殺蟲活性及機制，本研究在測定菌株YZ011對小菜蛾和草地貪夜蛾血腔注射毒性的基礎(chǔ)上，進一步對其進行全基因組測序和比較基因組學分析。結(jié)果表明，血腔注射100 nL 1.0×10⁶CFU/mL菌株YZ011 24 h后，小菜蛾幼蟲和草地貪夜蛾幼蟲的死亡率分別為75.0%和56.7%，均顯著高于血腔注射磷酸鹽緩沖液（PBS）對照。全基因組測序結(jié)果顯示菌株YZ011為嗜線蟲致病桿菌，其基因組含有19個次級代謝產(chǎn)物基因簇；與其同源性最高的嗜線蟲致病桿菌菌株YL001和菌株SⅡ相比，菌株YZ011含有29個特異基因簇。菌株YZ011基因組中共注釋到586個毒力因子相關(guān)基因。本研究結(jié)果從基因組層面解析了嗜線蟲致病桿菌YZ011的分子特征，對菌株YZ011的綜合利用具有重要意義。

關(guān)鍵詞： 嗜線蟲致病桿菌；殺蟲活性；基因組；次級代謝產(chǎn)物；毒力因子

中圖分類號： S476"" 文獻標識碼： A"" 文章編號： 1000-4440（2025）02-0258-10

Insecticidal activity and comparative genomics analysis of Xenorhabdus nematophila YZ011 against lepidopteran insects

LIN Manman1， LIU Qin1， HAN Guangjie1， HUANG Lixin1， XIA Yang1， LI Chuanming1， ZHANG Nan2， LU Yurong1， XU Bin1， XU Jian^1，2

（1.Institute of Agricultural Sciences of the Lixiahe District in Jiangsu Province/National Agricultural Experimental Station for Agricultural Microbiology in Yangzhou， Yangzhou 225007， China；2.College of Plant Protection， Yangzhou University， Yangzhou 225009， China）

Abstract： Entomopathogenic nematode symbionts are important biological control resources. In order to clarify the insecticidal toxicity and mechanism of the symbiotic bacteria YZ011 isolated from Steinernema carpocapsae on lepidopteran pests， the whole genome sequencing and comparative genomics analysis of strain YZ011 were further carried out on the basis of measuring the hemocoel injection toxicity of strain YZ011 to Plutella xylostella and Spodoptera frugiperda. The results showed that the mortality rates of Plutella xylostella larvae and Spodoptera frugiperda larvae were 75.0% and 56.7%， respectively， after hemocoel injection of 100 nL 1.0×10⁶ CFU/mL strain YZ011 for 24 h， which were significantly higher than those of phosphate buffered saline （PBS） control. The whole genome sequencing results showed that YZ011 was a nematophagous pathogenic bacterium， and its genome contained 19 secondary metabolite gene clusters. Compared with the most homologous strains YL001 and S Ⅱ， YZ011 contained 29 specific gene clusters， and 586 virulence factor-related genes were annotated. The results of this study analyzed the molecular characteristics of Xenorhabdus nematophila YZ011 at the genomic level， which was of great significance for the comprehensive utilization of strain YZ011.

Key words： Xenorhabdus nematophila;insecticidal activity;genome;secondary metabolite;virulence factor

昆蟲病原線蟲是廣泛存在于自然環(huán)境中的一類寄生性線蟲，它與體內(nèi)共生菌互惠共生，具有天敵昆蟲兼病原微生物雙重特點，能通過感染害蟲，并利用其作為寄主進行繁殖和生存，最終導致昆蟲死亡^[1]。昆蟲病原線蟲共生菌存在于昆蟲病原線蟲腸道，屬革蘭氏陰性菌，包括與斯氏線蟲科（Steinernema）線蟲共生的致病桿菌屬（Xenorhabdus）菌和與異小桿線蟲科（Heterorhabditis）線蟲共生的發(fā)光桿菌屬（Photorhabdus）菌^[2-3]。昆蟲病原線蟲共生菌與線蟲專一性共生，線蟲入侵昆蟲血腔后致病桿菌由線蟲肛門排出，發(fā)光桿菌被線蟲反芻從口腔吐出。共生菌進入昆蟲血腔后，產(chǎn)生的一系列次級代謝產(chǎn)物能導致昆蟲死亡，分解昆蟲組織，為線蟲的生長與繁殖提供營養(yǎng)物質(zhì)^[2，4]。作為一種重要的潛在生物殺蟲資源，越來越多具有生物活性的昆蟲病原線蟲共生菌菌株被發(fā)現(xiàn)。截至2024年7月，美國國家生物技術(shù)信息中心（NCBI）分類數(shù)據(jù)庫（https：//www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi）中共收錄了32個分類物種的致病桿菌菌株和24個分類物種的發(fā)光桿菌菌株，目前尚有166個致病桿菌菌株和116個發(fā)光桿菌菌株有待分類^[5]。

共生菌產(chǎn)生的次級代謝產(chǎn)物主要為聚酮合成酶（PKS）、非核糖體肽合成酶（NRPS）和其他類似酶合成的非核糖體衍生肽、聚酮化合物和雜合天然產(chǎn)物，合成這些化合物的基因占共生菌基因組的7.5%^[6-9]。由于共生菌次級代謝產(chǎn)物的抗菌和殺蟲潛力，其在藥品和生物農(nóng)藥的研制中廣受關(guān)注。Palma等^[10]對分離自斯氏線蟲的14株致病桿菌進行生物活性測定和基因組分析，發(fā)現(xiàn)14株致病桿菌均對大蠟螟具有毒殺作用，基因組生物信息分析獲得110種殺蟲蛋白，包括Tc、Txp、Mcf、Pra/Prb和App同系物等，以及其他毒力因子，如殺線蟲蛋白、幾丁質(zhì)酶和用于合成不同生物活性化合物的次級代謝產(chǎn)物。Wang等^[11]分離的122株共生菌對玉米螟幼蟲均有口服殺蟲活性或生長抑制作用，其中，致病桿菌菌株SY5對小菜蛾、甜菜夜蛾、玉米螟、黏蟲和黃粉蟲有較強的殺蟲活性。共生菌產(chǎn)生的具有殺蟲活性的蛋白類毒素根據(jù)其作用方式可分為血腔注射活性和胃毒活性兩類。血腔活性毒素蛋白包括Mcf毒素^[12]、Txp40毒素^[13]、XaxAB毒素^[14]、XhlA毒素^[15]、Pir毒素^[16-17]和RTX毒素^[18]等，胃毒活性殺蟲蛋白主要有Tc毒素^[19]、幾丁質(zhì)酶^[20]和XnGroEL毒素^[21]等。

劉琴等^[22]從罹死二化螟蟲體中分離獲得對二化螟、小菜蛾、草地貪夜蛾等具有廣譜寄生活性的小卷蛾斯氏線蟲及其共生菌YZ011。為明確菌株YZ011功能基因、基因簇及殺蟲活性物質(zhì)，本研究在測定共生菌YZ011對小菜蛾幼蟲和草地貪夜蛾幼蟲殺蟲活性的基礎(chǔ)上，進一步對其進行全基因組測序及生物信息學分析，注釋菌株YZ011基因功能，預測菌株YZ011的次級代謝產(chǎn)物基因簇及毒力因子，以期為菌株YZ011在害蟲生物防治中的應用提供基礎(chǔ)和依據(jù)。

1 材料與方法

1.1 試驗材料

共生菌菌株YZ011分離自罹死二化螟幼蟲體內(nèi)的小卷蛾斯氏線蟲，以甘油-80 ℃保存于江蘇里下河地區(qū)農(nóng)業(yè)科學研究所^[22]。小菜蛾為揚州種群，室內(nèi)用油菜苗飼養(yǎng)繁殖多代^[23]。草地貪夜蛾蟲源引自中國農(nóng)業(yè)科學院植物保護研究所，人工飼料室內(nèi)多代飼養(yǎng)繁殖^[22]。

1.2 菌株YZ011對小菜蛾、草地貪夜蛾的殺蟲毒力

取甘油保存的YZ011菌株10 μL轉(zhuǎn)接至1 mL 胰酪大豆胨液體培養(yǎng)基（TSB）^[22]中，28 ℃、220 r/min振蕩培養(yǎng)16 h后，取500 μL培養(yǎng)液轉(zhuǎn)接至50 mL TSB培養(yǎng)基再在28 ℃、220 r/min振蕩培養(yǎng)24 h，利用平板計數(shù)法測定菌液中菌體含量。菌液離心后用磷酸鹽緩沖液（PBS）重新懸浮菌體至1.0×10⁶CFU/mL。選取4齡小菜蛾幼蟲和草地貪夜蛾幼蟲，于解剖鏡下，利用Nanoject Ⅲ顯微注射器（美國Drummond公司產(chǎn)品）向蟲體倒數(shù)第2～3腹節(jié)間分別注射50 nL和100 nL 1.0×10⁶CFU/mL共生菌，以注射100 nL無菌PBS為對照。小菜蛾飼養(yǎng)于直徑為9.0 cm盛有卷心菜的培養(yǎng)皿中，草地貪夜蛾單頭飼養(yǎng)于盛有人工飼料的24孔細胞培養(yǎng)板內(nèi)，每個處理20頭幼蟲，重復3次，生化培養(yǎng)箱25 ℃培養(yǎng)，光照時間為16 h/d，24 h后進行存活小菜蛾和草地貪夜蛾計數(shù)，計算死亡率。利用GraphPad Prism 9.5.0和IBM SPSS Statistics 21軟件進行數(shù)據(jù)統(tǒng)計和處理間差異顯著性分析（Plt;0.05）。

1.3 菌株YZ011全基因組測序及組分與功能分析

取菌株YZ011培養(yǎng)液50 μL轉(zhuǎn)接至5 mL TSB培養(yǎng)基中，培養(yǎng)16 h離心收集菌體，采用細菌DNA提取試劑盒（美國Omega公司產(chǎn)品）制備DNA樣本，委托南京集思慧遠生物科技有限公司進行菌株的全基因組測序。

采用Unicycler v0.4.8軟件組裝測序數(shù)據(jù)，獲得基因組序列。使用軟件Prokka v1.12軟件進行基因組組分分析，包括編碼基因和非編碼基因的注釋，參數(shù)使用--rfam--addgenes，其余參數(shù)取默認值。使用RepeatMasker 4.1.0軟件進行重復序列預測。采用綜合抗生素耐藥性數(shù)據(jù)庫（CARD）、碳水化合物活性酶數(shù)據(jù)庫（CAZy）、病原與宿主互作數(shù)據(jù)庫（PHI）進行基因功能注釋。利用EffectiveT3軟件預測Ⅲ型分泌系統(tǒng)效應蛋白，利用SignalP 4.1軟件進行分泌蛋白信號肽預測。利用BLAST軟件將預測的基因序列與非冗余蛋白質(zhì)氨基酸序列數(shù)據(jù)庫（NR）、Swiss-Prot蛋白氨基酸序列數(shù)據(jù)庫、京都基因與基因組百科全書（KEGG）、同源蛋白質(zhì)簇（COG）數(shù)據(jù)庫、真核生物同源基因組（KOG）數(shù)據(jù)庫、基因本體（GO）數(shù)據(jù)庫、蛋白質(zhì)家族（Pfam）數(shù)據(jù)庫進行比對，獲得預測基因的注釋信息。利用Circos v0.69-5軟件（https：//circos.ca）及COG注釋繪制菌株YZ011基因組的全基因組圈圖。

1.4 進化分析和比較基因組學分析

上傳菌株YZ011全基因組序列至Type Strain Genome Server（TYGS）平臺（https：//tygs.dsmz.de），分別以16S rDNA和全基因組構(gòu)建系統(tǒng)發(fā)育樹以確定菌株YZ011所屬物種?；谙到y(tǒng)發(fā)育樹分析結(jié)果，從NCBI（https：//www.ncbi.nlm.nih.gov/）獲取具有全基因組序列的嗜線蟲致病桿菌菌株基因組，上傳至TYGS平臺，分別以16S rDNA和全基因組構(gòu)建系統(tǒng)發(fā)育樹，分析獲得同源性高的菌株。利用Geneious 9.0.2和OrthoVenn3軟件（https：//orthovenn3.bioinfotoolkits.net/start/upload）分別對菌株YZ011及其同源性高的菌株進行全基因組基因共線性分析和直系同源核心基因簇及特異基因簇分析。

1.5 次級代謝產(chǎn)物基因簇分析

將菌株YZ011基因組上傳至Anti-SMASH^[24]（https：//antismash.secondarymetabolites.org/#！/start），對其進行次級代謝產(chǎn)物基因簇分析。

1.6 毒力因子分析

采用病原菌毒力因子數(shù)據(jù)庫（Virulence Factor Database， VFDB）（https：//www.mgc.ac.cn/VFs）對菌株YZ011進行毒力因子基因鑒定。利用Blastall 2.2.26軟件對基因序列進行注釋，期望域值（E）設(shè)置為1×10^-10。

2 結(jié)果與分析

2.1 菌株YZ011對小菜蛾和草地貪夜蛾幼蟲的殺蟲活性

菌株YZ011對鱗翅目害蟲小菜蛾和草地貪夜蛾幼蟲均具有感染致死作用，顯微注射50 nL和100 nL的1.0×10⁶CFU/mL共生菌24 h后，小菜蛾幼蟲的死亡率分別為60.0%和75.0%，均顯著高于注射PBS對照；顯微注射50 nL和100 nL 的1.0×10⁶CFU/mL共生菌24 h后，草地貪夜蛾幼蟲的死亡率分別為 35.0%和56.7%，亦顯著高于注射PBS對照（圖1）。且注射100 nL 1.0×10⁶CFU/mL共生菌24 h后，小菜蛾和草地貪夜蛾幼蟲的死亡率均顯著高于注射50 nL 1.0×10⁶CFU/mL共生菌處理。

2.2 菌株YZ011全基因組

菌株YZ011的全基因組圖譜如圖2所示。從圖中可以看出，共生菌菌株YZ011基因組大小為4 408 387 bp（GenBank 登錄號：CP158108），G+C堿基含量為43.9%，預測得到3 985個編碼基因、22個rRNA、76個tRNA、110個misc_RNA和1個tmRNA。

從外到內(nèi)的7個圈依次為基因組大小的標識、正鏈上的編碼序列（CDS）、負鏈上的CDS、正鏈上的rRNA和tRNA、負鏈上的rRNA和tRNA、G+C堿基含量和GC偏移值。

2.3 功能基因數(shù)

基于KEGG、COG、Swiss-Prot、GO和KOG等數(shù)據(jù)庫注釋得到的功能基因較多，分別為3 933個、2 907個、2 871個、2 678個和1 027個，占基因總數(shù)的99.37%、73.45%、72.54%和67.66%和25.95%；基于CARD數(shù)據(jù)庫注釋得到的基因有254個，占基因總數(shù)的6.42%，而基于CAZy數(shù)據(jù)庫注釋得到的基因最少，僅68個，占基因總數(shù)的1.72%。

2.4 菌株YZ011進化和全基因組比對

以菌株YZ011的16S rDNA和全基因組構(gòu)建的系統(tǒng)發(fā)育樹如圖3所示。2個系統(tǒng)發(fā)育樹均顯示，菌株YZ011與嗜線蟲致病桿菌（Xenorhabdus nematophila）處于同一分支，因此，本研究認為共生菌菌株YZ011屬于嗜線蟲致病桿菌。

2.5 比較基因組

從NCBI獲得4株具有全基因組序列的嗜線蟲致病桿菌菌株，菌株SⅡ（Xenorhabdus nematophila，GenBank 登錄號：CP060401）、菌株YL001（Xenorhabdus nematophila，GenBank登錄號：CP032329）、菌株AN6/1（Xenorhabdus nematophila，GenBank登錄號：LN681227）和菌株ATCC 19061（Xenorhabdus nematophila，GenBank 登錄號：FN667742）。將菌株YZ011的基因組序列與上述4株嗜線蟲致病桿菌的基因組序列構(gòu)建系統(tǒng)發(fā)育樹（圖4），從圖中可見菌株YZ011與菌株YL001同源性最高，其次為菌株SⅡ。

菌株YZ011與4株嗜線蟲致病桿菌基因組的DNA同源性值、平均核苷酸一致性、基因組大小和G+C堿基含量見表1。從表中可以看出，菌株YZ011與4個同源菌株的DNA同源性指數(shù)值d0、d4和d6均在88.4%以上，表明菌株YZ011與4株同源菌株基因組的同源性均較高。菌株YZ011與4株同源菌株的平均核苷酸一致性值為99.46%～99.98%，其中菌株YZ011與菌株YL001的平均核苷酸一致性最高（99.98%），與菌株SⅡ的平均核苷酸一致性次之（99.96%）。這一結(jié)果也與圖4結(jié)果相吻合。菌株YZ011與4株同源菌株的基因組大小和G+C堿基含量亦基本接近。

將菌株YZ011與同源性較高的菌株YL001和菌株SⅡ進行共線性分析，結(jié)果如圖5所示。從圖中可以看出，菌株YZ011與菌株YL001和菌株SⅡ有較多大片段的共線性，也存在基因顛倒重排和共線性較低的區(qū)域。

菌株YZ011與菌株YL001和菌株SⅡ的全基因組直系同源基因簇分析結(jié)果如圖6所示。菌株YZ011與菌株YL001和菌株SⅡ的核心基因簇有3 254個、非必須基因簇有269個以及特異基因簇42個（菌株YZ011 29個、菌株YL001 13個）。GO分析發(fā)現(xiàn)菌株YZ011的29個特異基因簇涉及2個條目（GO：0006313和GO：0032196），其中GO：0006313的基因簇數(shù)量為20個，GO：0032196基因簇數(shù)量為9個，主要參與轉(zhuǎn)座和DNA介導。

2.6 菌株YZ011次級代謝產(chǎn)物合成基因簇

菌株YZ011基因組次級代謝產(chǎn)物基因簇分析結(jié)果如表2所示。從表中可以看出，YZ011基因組含有19個次級代謝產(chǎn)物基因簇，包括8個非核糖體合成酶基因簇（NRPS）、2個類非核糖體合成酶基因簇（NRRS-like）、1個I型聚酮化合物合酶基因簇（T1PKS）、1個獨立于NRPS的鐵載體基因簇（NI-siderophore）、1個吩嗪生物合成基因簇（Phenazine）、1個丁內(nèi)酯生物合成基因簇（Butyrolactone）、1個含有β-內(nèi)酯的蛋白酶抑制劑生物合成基因簇（Betalactone）、1個核苷生物合成基因簇（Nucleoside）和3個未知功能的基因簇。Region 11的NRPS基因簇和Region 12未知功能基因簇位于菌株YZ011與菌株YL001和菌株SⅡ共線性較低的區(qū)域。

2.7 菌株YZ011毒力因子

基于VFDB數(shù)據(jù)庫，菌株YZ011基因組序列中共注釋到586個毒力因子相關(guān)基因，如大多數(shù)革蘭氏陰性病原細菌感染宿主的重要分泌系統(tǒng)Ⅲ型分泌系統(tǒng)（T3SS）、Ⅵ型分泌系統(tǒng)（T6SS）、脂多糖、抗菌物質(zhì)、溶血素、毒力增強因子、穿孔毒素、毒素重復序列RTX等類型的基因。部分毒力因子與其他革蘭氏陰性昆蟲病原菌的殺蟲蛋白具有一定的相似性，如溶血素HlyA與嗜水氣單胞菌溶血素AhlA^[25]的相似度為59.6%，穿孔毒素與耶爾森氏菌YaxAB毒素^[26]的相似度為53.8%。

3 討論與結(jié)論

嗜線蟲致病桿菌是與小卷蛾斯氏線蟲共生的昆蟲病原細菌，其產(chǎn)生的次級代謝產(chǎn)物具有抗真菌、抗細菌和殺蟲等活性^[27]。本團隊從二化螟罹死蟲體中分離獲得的小卷蛾斯氏線蟲N-Yz1對二化螟幼蟲和草地貪夜蛾幼蟲均有較高感染致死作用^[28]，進一步從小卷蛾斯氏線蟲體內(nèi)分離獲得了共生菌菌株YZ011，生物活性測定結(jié)果表明菌株YZ011對小菜蛾幼蟲和草地貪夜蛾幼蟲均具有較高的毒力。在全基因組測序基礎(chǔ)上，基于16S rDNA和全基因組序列構(gòu)建的系統(tǒng)進化樹均表明菌株YZ011屬于嗜線蟲致病桿菌，且該菌株基因組大小、G+C堿基含量與其他已報道的嗜線蟲致病桿菌菌株一致。通過與4株已有全基因組序列的嗜線蟲致病桿菌進行基因組比較分析發(fā)現(xiàn)，菌株YZ011與嗜線蟲致病桿菌菌株YL001相似度最高。菌株YL001不僅對馬鈴薯晚疫病病原菌（Phytophthora infestans）、玉米大斑病病原菌（Exserohilum turcicum）、水稻紋枯病病原菌（Rhizoctonia solani）等具有較強抗菌活性^[29]，其發(fā)酵上清液中還含有水溶性苯并芘類化合物Xenocoumacins和酰胺類化合物Nematophin等殺蟲活性物質(zhì)，而在其他嗜線蟲致病桿菌中亦分離得到Xenocoumacins、Nematophin、PAX肽（Peptide-antimicrobial-Xenorhabdus）、Benzylideneacetone、Xenomatides和Xenematides等^{[8，30-31]}殺蟲活性物質(zhì)，但嗜線蟲致病桿菌代謝產(chǎn)物中發(fā)揮殺蟲作用的具體化合物還不明確。

嗜線蟲致病桿菌不能獨立于其宿主長時間存活，因此很難將共生菌直接用作生物農(nóng)藥。基于全基因組測序技術(shù)和生物信息學方法，對致病桿菌的毒力因子進行分析對推進共生菌在害蟲防控領(lǐng)域的應用具有重要意義。本研究對菌株YZ011次級代謝產(chǎn)物合成基因簇及其毒力因子進行了分析，為后續(xù)克隆相關(guān)基因獲得對鱗翅目害蟲具有殺蟲活性的毒素提供基礎(chǔ)。聚酮類化合物（Polyketide， PK）和非核糖體多肽（Nonribosonmal peptide， NRP）是目前研究最多的兩類重要的次級代謝產(chǎn)物，其生物合成需要PK和NRP。PK合成多種聚酮化合物，包括抗生素、抗腫瘤物質(zhì)、免疫抑制劑等，NRP則合成多種天然多肽，如免疫抑制劑、抗生素、抗癌和抗病毒因子、殺蟲劑、鐵載體和生物表面活性劑等^[6-7]。對菌株YZ011預測得到的19個次級代謝產(chǎn)物基因簇中NRP基因簇達9個，表明菌株YZ011中可能存在毒素物質(zhì)合成基因。值得注意的是，次級代謝產(chǎn)物基因簇中Region 11的NRP基因簇位于菌株YZ011與菌株YL001和菌株SⅡ共線性較低的區(qū)域，該共線性較低的區(qū)域還包括Region 12這一未知功能基因簇，說明該區(qū)域可能存在新的潛在活性物質(zhì)基因?；赩FDB數(shù)據(jù)庫，菌株YZ011基因組中注釋到多種毒素、胞外酶和抗菌劑基因，其中一些同源基因在毒殺昆蟲、營養(yǎng)獲取和免疫逃避等方面發(fā)揮重要作用^[32-33]。毒力因子與特定的分泌系統(tǒng)有關(guān)，革蘭氏陰性菌特有的分泌機制可以將蛋白質(zhì)直接從細胞質(zhì)分泌到細胞外或直接進入靶細胞，包括Ⅰ型分泌系統(tǒng)（T1SS）、T3SS和T6SS。其中，T1SS能產(chǎn)生相對分子量較大的毒素蛋白，如RTX蛋白家族的蛋白^[34]。菌株YZ011注釋到9個由T1SS分泌的RTX家族蛋白。T3SS在致病桿菌侵染宿主過程中發(fā)揮著重要作用，病原菌通過T3SS將效應蛋白注射到宿主細胞中，進而調(diào)控其致病力。目前已報道發(fā)光桿菌菌株均編碼至少1個T3SS，而致病桿菌卻不一定存在T3SS^{[9， 35-38]}。另外，致病桿菌通過T6SS分泌系統(tǒng)將毒素分子直接遞送到靶細胞，致病桿菌基因組至少編碼1個T6SS，具有多個T6SS的菌株在與宿主的相互作用中更具有競爭優(yōu)勢^[35，39]。本研究中，菌株YZ011基因組中預測得到5個T3SS效應蛋白基因和32個T6SS效應蛋白基因，說明菌株YZ011在害蟲防控領(lǐng)域具有較高的應用潛力。

不同菌株的分泌系統(tǒng)存在一定差異，菌株分泌的活性物質(zhì)亦變化較大，甚至在相同物種的菌株之間亦是如此^[40]。菌株YZ011毒力因子注釋結(jié)果表明，菌株YZ011中含有多個已報道嗜線蟲致病桿菌毒力因子，同時也存在一定差異。Txp40是一種在昆蟲病原線蟲共生菌中普遍存在且高度保守的殺蟲毒素蛋白，可導致昆蟲中腸腸壁細胞間排列紊亂，對多種鱗翅目和雙翅目昆蟲具有毒殺作用，Kinkar等^[13]研究結(jié)果表明重組Txp40對大蠟螟幼蟲具有致死毒性。本研究在菌株YZ011沒有注釋到Txp40毒素蛋白，但是菌株YZ011對小菜蛾和草地貪夜蛾幼蟲的血腔注射毒性預示其必然存在某些毒力物質(zhì)，對注釋到的其他毒力因子XhlA、Cytolysin等在Cowles等^[15]、Mcquade等^[41]研究中已有報道。共生菌菌株YZ011對小菜蛾幼蟲和草地貪夜蛾幼蟲的殺蟲毒力受哪些毒力基因調(diào)控尚需進一步分析與驗證。

綜上，本研究對分離自小卷蛾斯氏線蟲的嗜線蟲致病桿菌YZ011進行了全基因組序列分析，對其潛在的毒力基因進行了預測，后續(xù)將進一步克隆毒力基因進行異源表達以及對靶標昆蟲生物活性測定，以期為害蟲生物防治開發(fā)新型生物農(nóng)藥提供技術(shù)支撐。

參考文獻：

[1] KOPPENH FER A M， SHAPIRO-ILAN D I， HILTPOLD I. Entomopathogenic nematodes in sustainable food production[J]. Frontiers in Sustainable Food Systems，2020，4：125.

[2] CIMEN H R， TOURAY M， GULSEN S H， et al. Natural products from Photorhabdus and Xenorhabdus：mechanisms and impacts[J]. Applied Microbiology and Biotechnology，2022，106（12）：4387-4399.

[3] 段秀峰，林茂松，李紅梅，等. 伯氏致病桿菌（Xenor habdus bovienii BE）殺蟲活性物質(zhì)的分離和生物測定[J]. 江蘇農(nóng)業(yè)學報，2009，25（1）：59-67.

[4] ABD-ELGAWAD M M M. Xenorhabdus spp. ：an overview of the useful facets of mutualistic bacteria of entomopathogenic nematodes[J]. Life，2022，12（9）：1360.

[5] SCHOCH C L， CIUFO S， DOMRACHEV M， et al. NCBI Taxonomy：a comprehensive update on curation，resources and tools[J]. Database，2020，2020：baaa062.

[6] SHI Y M， BODE H B. Chemical language and warfare of bacterial natural products in bacteria-nematode-insect interactions[J]. Natural Product Reports，2018，35（4）：309-335.

[7] CAI X F， NOWAK S， WESCHE F， et al. Entomopathogenic bacteria use multiple mechanisms for bioactive peptide library design[J]. Nature Chemistry，2017，9（4）：379-386.

[8] TOBIAS N J， WOLFF H， DJAHANSCHIRI B， et al. Natural product diversity associated with the nematode symbionts Photorhabdus and Xenorhabdus[J]. Nature Microbiology，2017，2（12）：1676-1685.

[9] CHASTON J M， SUEN G， TUCKER S L， et al. The entomopathogenic bacterial endosymbionts Xenorhabdus and Photorhabdus：convergent lifestyles from divergent genomes[J]. PLoS One，2011，6（11）：e27909.

[10]PALMA L， FRIZZO L， KAISER S， et al. Genome sequence analysis of native Xenorhabdus strains isolated from entomopathogenic nematodes in Argentina[J]. Toxins，2024，16（2）：108.

[11]WANG H， DONG H， QIAN H T， et al. Isolation，bioassay and characterisation of Xenorhabdus sp. SY5，a highly virulent symbiotic bacterium of an entomopathogenic nematode isolated from China[J]. Nematology，2013，15（2）：153-163.

[12]DABORN P J， WATERFIELD N， SILVA C P， et al. A single Photorhabdus gene，makes caterpillars floppy （mcf），allows Escherichia coli to persist within and kill insects[J]. Proceedings of the National Academy of Sciences of the United States of America，2002，99（16）：10742-10747.

[13]KINKAR O U， PRASHAR A， KUMAR A， et al. Txp40，an insecticidal toxin protein from Xenorhabdus nematophila：purification，toxicity assessment and biophysical characterization[J]. Toxicon，2022，218：40-46.

[14]VIGNEUX F， ZUMBIHL R， JUBELIN G， et al. The xaxAB genes encoding a new apoptotic toxin from the insect pathogen Xenorhabdus nematophila are present in plant and human pathogens[J]. Journal of Biological Chemistry，2007，282（13）：9571-9580.

[15]COWLES K N， GOODRICH-BLAIR H. Expression and activity of a Xenorhabdus nematophila haemolysin required for full virulence towards Manduca sexta insects[J]. Cellular Microbiology，2005，7（2）：209-219.

[16]楊 晴，張 杰，李天慧，等. 嗜線蟲致病桿菌屬Pir毒素基因克隆和生物信息學分析[J]. 中國植保導刊，2017，37（9）：5-11.

[17]AHANTARIG A， CHANTAWAT N， WATERFIELD N R， et al. PirAB toxin from Photorhabdus asymbiotica as a larvicide against dengue vectors[J]. Applied and Environmental Microbiology，2009，75（13）：4627-4629.

[18]SATCHELL K J F. Structure and function of MARTX toxins and other large repetitive RTX proteins[J]. Annual Review of Microbiology，2011，65：71-90.

[19]RODERER D， BR CKER F， SITSEL O， et al. Glycan-dependent cell adhesion mechanism of Tc toxins[J]. Nature Communications，2020，11（1）：2694.

[20]LIU J， NANGONG Z Y， ZHANG J， et al. Expression and characterization of two chitinases with synergistic effect and antifungal activity from Xenorhabdus nematophila[J]. World Journal of Microbiology amp; Biotechnology，2019，35（7）：106.

[21]SHI H X， ZENG H M， YANG X F， et al. An insecticidal protein from Xenorhabdus ehlersii triggers prophenoloxidase activation and hemocyte decrease in Galleria mellonella[J]. Current Microbiology，2012，64（6）：604-610.

[22]劉 琴，黃立鑫，陸玉榮，等. 小卷蛾斯氏線蟲N-Yz1對草地貪夜蛾的致毒作用和感染特性[J]. 中國生物防治學報，2023，39（3）：542-549.

[23]祁建杭，劉 琴，劉雪梅，等. KNO3提高蘇云金芽胞桿菌Bt-8發(fā)酵液殺蟲活性的作用分析[J]. 揚州大學學報（農(nóng)業(yè)與生命科學版），2022，43（5）：136-141.

[24]BLIN K， SHAW S， AUGUSTIJN H E， et al. AntiSMASH 7. 0：new and improved predictions for detection，regulation，chemical structures and visualisation[J]. Nucleic Acids Research，2023，51（W1）：46-50.

[25]WILSON J S， CHURCHILL-ANGUS A M， DAVIES S P， et al. Identification and structural analysis of the tripartite α-pore forming toxin of Aeromonas hydrophila[J]. Nature Communications，2019，10（1）：2900.

[26]WAGNER N J， LIN C P， BORST L B， et al. YaxAB，a Yersinia enterocolitica pore-forming toxin regulated by RovA[J]. Infection and Immunity，2013，81（11）：4208-4219.

[27]張 媛，王 芳，趙紫華. 斯氏線蟲-致病桿菌共生體侵染昆蟲的生物標志物[J]. 植物保護學報，2023，50（5）：1227-1236.

[28]劉 琴，韓光杰，李傳明，等. 昆蟲病原線蟲N-Yz1對水稻二化螟的控制作用[J]. 中國生物防治學報，2024，40（5）：1015-1021.

[29]ZHANG S J， LIU Q， HAN Y F， et al. Nematophin，an antimicrobial dipeptide compound from Xenorhabdus nematophila YL001 as a potent biopesticide for Rhizoctonia solani control[J]. Frontiers in Microbiology，2019，10：1765.

[30]INCEDAYI G， CIMEN H R， ULUG D， et al. Relative potency of a novel acaricidal compound from Xenorhabdus，a bacterial genus mutualistically associated with entomopathogenic nematodes[J]. Scientific Reports，2021，11（1）：11253.

[31]DREYER J， MALAN A P， DICKS L M T. Bacteria of the genus Xenorhabdus，a novel source of bioactive compounds[J]. Frontiers in Microbiology，2018，9：3177.

[32]CALDAS C， CHERQUI A， PEREIRA A， et al. Purification and characterization of an extracellular protease from Xenorhabdus nematophila involved in insect immunosuppression[J]. Applied and Environmental Microbiology，2002，68（3）：1297-1304.

[33]DARSOUEI R， KARIMI J， DUNPHY G B. Functional characterization of outer membrane proteins （OMPs） in Xenorhabdus nematophila and Photorhabdus luminescens through insect immune defense reactions[J]. Insects，2019，10（10）：352.

[34]KANONENBERG K， SCHWARZ C K W， SCHMITT L. Type I secretion systems-a story of appendices[J]. Research in Microbiology，2013，164（6）：596-604.

[35]SILVA C P， WATERFIELD N R， DABORN P J， et al. Bacterial infection of a model insect：Photorhabdus luminescens and Manduca sexta[J]. Cellular Microbiology，2002，4（6）：329-339.

[36]OGIER J C， PAGS S， BISCH G， et al. Attenuated virulence and genomic reductive evolution in the entomopathogenic bacterial symbiont species，Xenorhabdus poinarii[J]. Genome Biology and Evolution，2014，6（6）：1495-1513.

[37]MURFIN K E， WHOOLEY A C， KLASSEN J L， et al. Comparison of Xenorhabdus bovienii bacterial strain genomes reveals diversity in symbiotic functions[J]. BMC Genomics，2015，16：889.

[38]HINCHLIFFE S J. Insecticidal toxins from the Photorhabdus and Xenorhabdus bacteria[J]. The Open Toxinology Journal，2013，3（1）：101-118.

[39]KIM I H， ARYAL S K， AGHAI D T， et al. The insect pathogenic bacterium Xenorhabdus innexi has attenuated virulence in multiple insect model hosts yet encodes a potent mosquitocidal toxin[J]. BMC Genomics，2017，18（1）：927.

[40]BODE E， ASSMANN D， HAPPEL P， et al. EasyPACId，a simple method for induced production，isolation，identification，and testing of natural products from proteobacteria[J]. Bio-protocol，2023，13（13）：e4709.

[41] MCQUADE R， STOCK S P. Secretion systems and secreted proteins in gram-negative entomopathogenic bacteria：their roles in insect virulence and beyond[J]. Insects，2018，9（2）：68.

（責任編輯：石春林）