摘要：為探明造成養(yǎng)殖沙塘鱧腸道出血的病原菌G-2的特性，對G-2進行全基因組測序和分析?；蚪M組分分析和基因功能注釋結果顯示：G-2的總基因數(shù)為5 134個，其中蛋白編碼基因4 977個，非蛋白編碼基因157個，非編碼基因中tRNA基因107個，rRNA基因42個，sRNA基因8個；分別有5 109、3 339、3 362、2 543個基因在對應的Nr、SwissProt、COG、KEGG數(shù)據(jù)庫中獲得注釋；毒力基因分析表明G-2含有多種腸毒素和溶血素基因、InhA蛋白、磷脂酶及蛋白質毒素，解釋了該菌造成養(yǎng)殖沙塘鱧腸道出血的原因?；贕-2菌株的16S rRNA序列信息構建系統(tǒng)進化樹，確定其為蘇云金芽孢桿菌。

關鍵詞：沙塘鱧；全基因組測序；蘇云金芽孢桿菌

中圖分類號：S943 文獻標識碼：A 文章編號：1006-060X（2024）09-0008-07

Whole Genome Sequencing and Bioinformatics Analysis of a Pathogenic Strain G-2 in the Gut of Odontobutis potamophila

ZHAO Yan-hua1，2，XIA Ai-jun1，2，JIANG Hu-cheng1，2，XUE Hui1，2，LIU Guo-xing1，2

[1. Jiangsu Provincial Freshwater Fisheries Research Institute， Nanjing 210017， PRC; 2. Low-Temperature Germplasm Bank of

Important Economic Fish of Jiangsu Provincial Science and Technology Resources （Agricultural Germplasm Resources）

Coordination Service Platform， Nanjing 210017， PRC]

Abstract： Whole genome sequencing and bioinformatics analysis were carried out to characterize the pathogenic strain G-2 causing gut bleeding of cultured Odontobutis potamophila. The results showed that the genome of strain G-2 contained 5 134 genes， including 4 977 coding genes and 157 non-coding genes （107 tRNA genes， 42 rRNA genes， and 8 sRNA genes）. In Nr， SwissProt， COG， and KEGG， 5 109， 3 339， 3 362， and 2 543 genes， respectively， were annotated. Strain G-2 carried multiple enterotoxin and hemolysin genes and the genes encoding InhA protein， phospholipase， and protein toxins， which explained the gut bleeding caused by this strain. The phylogenetic tree built based on the 16S rRNA gene sequence identified strain G-2 as Bacillus thuringiensis.

Key words： Odontobutis potamophila; whole genome sequencing; Bacillus thuringiensis

全基因組測序（Whole genome sequencing）是對未知基因組序列的物種進行個體的基因組測序，通過將生物的基因組打斷、測序和拼接獲得較為完整的基因組序列，繼而進行序列分析和功能基因注釋，最終從基因水平深入了解生物的各項機理[1-2]。1986年，Renato Dulbecco等科學家最早提出人類基因組測序計劃。1995年，流感嗜血桿菌（Haemophilus influenzae）成為第一個完成全基因組測序的自由生存生物[3]。隨著測序技術的發(fā)展和測序成本的降低，目前全基因組測序已經(jīng)經(jīng)歷了三代技術發(fā)展，對任何生物的基因組進行測序成為可能[4-5]。在病害防控領域，全基因組測序技術可以有效獲取病原菌較為完整的基因組序列，再結合功能基因組學和生物信息學手段，為快速了解致病菌的致病機理和耐藥機制提供了極大的幫助[3]。

蘇云金芽孢桿菌（Bacillus thuringiensis）是一種革蘭氏陽性細菌，可做微生物源殺蟲劑，以胃毒作用為

主[6]。其中，蘇云金芽孢桿菌以色列亞種（Bacillus thuringiensis subsp. israelensis）是一種低毒殺蟲劑，會在芽孢形成過程中產(chǎn)生名為δ-內毒素的殺蟲伴孢晶體蛋白，殺蟲晶體蛋白分為Cry和Cyt兩大類[7]。Cyt蛋白具有溶血和溶細胞特性，對脊椎動物和無脊椎動物均有溶解細胞的毒害作用[8-10]。蘇云金芽孢桿菌也是典型的昆蟲病原菌，對多種昆蟲、線蟲和螨類等具有特異性的毒殺活性，在生物農(nóng)藥領域占據(jù)重要地位，是目前世界上研究最深入、應用最廣泛的微生物殺蟲劑[11-12]。

為探明造成養(yǎng)殖沙塘鱧腸道出血的病原菌G-2的特性，解析G-2的致病機理，該研究結合三代PacBio RS II測序技術和二代Illumina HiSeq 2000測序技術，對病原菌G-2進行全基因組測序分析，并對獲得的數(shù)據(jù)進行Nr（Non-Redundant Protein Sequence

Database）、Swiss-Prot、KEGG（Kyoto Encyclopedia of Genes and Genomes）、GO（Gene Ontology）、COG（Cluster of Orthologous Groups of proteins）、VFDB（Virulence Factors of Pathogenic Bacteria）數(shù)據(jù)庫注釋及分析，明確其基因組遺傳特性，為沙塘鱧的病害防控提供研究基礎。

1 材料與方法

1 供試材料

1.1.1 致病菌 菌株G-2來自江蘇省揚中市沙塘鱧養(yǎng)殖基地，由江蘇省淡水水產(chǎn)研究所病害防治研究室分離、鑒定及保存。

1.1.2 主要試劑和儀器 CTAB試劑（廣州化學試劑廠），Luria-Bertani（LB）培養(yǎng)基、腦心浸液肉湯培養(yǎng)基（廣東環(huán)凱生物科技有限公司），DNA提取試劑盒[生工生物工程（上海）股份有限公司]，瓊脂糖凝膠電泳儀（北京六一生物科技有限公司），NanoDrop微量分光光度計（美國賽默飛世爾科技公司），移液器（德國艾本德股份公司）。

1.2 研究方法

1.2.1 G-2的復蘇及擴大培養(yǎng) 將實驗室凍存的G-2菌株于室溫條件下解凍，并用移液槍吸取100 mL菌液加入腦心浸液肉湯培養(yǎng)基中，置于30℃恒溫培養(yǎng)箱中培養(yǎng)24 h；劃線接種于LB固體培養(yǎng)基，置于同樣培養(yǎng)條件下培養(yǎng)24 h；挑取長勢良好的單菌落接種于LB液體培養(yǎng)基，進行擴大培養(yǎng)，1 000 r/min轉速下離心10 min，收集菌體。

1.2.2 基因組提取、測序及組裝 使用CTAB法提取基因組DNA，將基因組DNA置于-20℃冰箱保存。對提取的DNA進行質量檢測，DNA樣品檢測合格后，進行文庫構建和文庫檢測，庫檢合格后，送至廣州基迪奧科技服務有限公司進行測序。結合三代PacBio RS II測序技術和二代Illumina HiSeq 2000測序技術[13]，對病原菌G-2進行全基因組測序分析。使用三代測序數(shù)據(jù)進行組裝，保證基因組組裝的完整性；使用二代測序數(shù)據(jù)對組裝結果進行校正，保證組裝結果的準確性；使用Fastp軟件對二代測序數(shù)據(jù)進行質控，基于校正后的組裝結果進行基因組組分分析。

1.2.3 基因組組分分析 使用NCBI數(shù)據(jù)庫（National Center for Biotechnology Information）進行基因預測，使用RNAmmer軟件預測基因組中的rRNA基因[14]，使用tRNAscan-SE軟件預測基因組中的tRNA基因，基于Rfam數(shù)據(jù)庫使用CMScan軟件預測基因組中的sRNA基因[15-16]，使用tRNAscan-SE軟件預測tRNA區(qū)域和tRNA的二級結構[17]，使用RepeatMasker軟件預測基因組重復序列，使用TRF軟件預測基因組的串聯(lián)重復序列[18]，使用CRISPRfinder軟件預測基因組的CRISPR區(qū)域，使用在線工具IslandViewer 4對基因組進行基因島預測。

1.2.4 基因功能注釋 使用KEGG數(shù)據(jù)庫（http：//www.genome.jp/kegg/）進行編碼基因功能預測，查找基因所在的生物學通路[19-20]；使用GO數(shù)據(jù)庫（http：//geneontology.org/）對基因和蛋白質功能進行限定和描述[21]；基于COG數(shù)據(jù)庫（http：//www.ncbi.nlm.

nih.gov/COG/）使用DIAMOND軟件進行數(shù)據(jù)比對，進行蛋白功能分類；利用VFDB數(shù)據(jù)庫（http：//www.mgc.ac.cn/）進行BLAST比對，預測毒力因子[22]。

1.2.5 G-2系統(tǒng)發(fā)育樹的構建 基于G-2菌株的16S rRNA序列信息，從GenBank中選取親緣關系較近的10個序列，利用MEGA 4軟件進行序列比對，采用鄰接法分析同源性，構建系統(tǒng)發(fā)育樹。

2 結果與分析

2.1 G-2基因組測序及組裝分析

使用三代測序平臺（PacBio RS II）對G-2進行測序，G-2過濾后的有效reads數(shù)為1 618 433條，平均測序讀長為4 359.5 bp，N50為6 698 bp，組裝后總長度為5 270 034 bp，GC含量為35.36%。

2.2 G-2基因分析及預測

使用二代測序對三代測序數(shù)據(jù)進行校正，最終獲得一條4 400 748 bp的環(huán)狀染色體基因組（圖1）。該基因組的總基因數(shù)為5 134個，其中蛋白編碼基因4 977個，編碼區(qū)域長度為4 327 191 bp，占總長度的98.3%；非蛋白編碼基因157個，其中tRNA基因107個，rRNA基因42個，sRNA基因8個。rRNA基因中23S_rRNA基因14個，16S_rRNA基因14個，5S_rRNA基因14個?；蜷g隔區(qū)域長度為73 557 bp，占總長度的1.7%。菌株G-2基因組的GC-depth分布圖如圖2所示，發(fā)現(xiàn)其散點分布較為集中，無雜菌污染。由G-2的基因長度分布圖（圖3）可知，基因長度小于1 000 bp的基因出現(xiàn)頻率較高，基因長度為2 500～3 000 bp的基因出現(xiàn)的頻率較低。

2.3 基因功能注釋分析

2.3.1 數(shù)據(jù)庫注釋比對 菌株G-2全基因組測序共獲得5 134個編碼基因，具有明確生物學功能的編碼基因5 110個。如圖4所示，將獲得的編碼基因與相應數(shù)據(jù)庫進行比對，分析結果顯示：有5 109個基因在對應的Nr數(shù)據(jù)庫獲得注釋，有3 339個基因在對應的SwissProt數(shù)據(jù)庫獲得注釋，有3 362個基因在對應的COG數(shù)據(jù)庫獲得注釋，有2 543個基因在對應的KEGG數(shù)據(jù)庫獲得注釋。

2.3.2 KEGG代謝通路分析 菌株G-2在KEGG數(shù)據(jù)庫中得到注釋的基因有2 543個。由圖5可知，在代謝類中，有12種通路得到注釋，其中全局總覽的基因個數(shù)最多，有803個；其次為氨基酸代謝通路和碳水化合物代謝通路，分別有270和249個基因；甘氨酸生物合成和代謝通路最少，只有43個基因。在遺傳信息處理類中，有4種通路得到注釋，其中翻譯通路的基因有92個，轉錄通路的基因有5個。在環(huán)境信息處理類中，有2種通路得到注釋，其中膜運輸通路有155個基因。在細胞過程類中，有3種通路得到注釋，其中原核生物基因數(shù)量最多，為107個，細胞生長和死亡通路最少，只有16個基因。2種生物系統(tǒng)類通路得到注釋，環(huán)境適應性和免疫系統(tǒng)的基因分別為4、3個。

2.3.3 GO注釋 GO數(shù)據(jù)庫是用于描述基因和蛋白質功能的標準化分類系統(tǒng)，包括生物過程、細胞組分和分子功能3種基因功能[23]。對蘇云金芽孢桿菌G-2的基因進行GO注釋和分類，判斷其基因參與的主要功能。如圖6所示，共有4 130個基因在GO數(shù)據(jù)庫中得到注釋，在生物過程中，有3 792個基因得到注釋，分類到細胞過程中的基因占比最多；在細胞組分分類中，有2 706個基因得到注釋，分類到細胞和細胞組成中的基因占比最多；在分子功能分類中，有3 767個基因得到注釋，分類到催化活性中的基因占比最多。

2.3.4 COG功能分類 由圖7可知，COG分類的蛋白編碼基因共3 362個，其中，常規(guī)功能預測類基因608個，氨基酸轉運和代謝類基因444個，轉錄類基因365個，無機離子轉運和代謝類基因264個，輔酶轉運和代謝類基因244個。

2.3.5 VFDB注釋 通過VFDB數(shù)據(jù)庫對比，發(fā)現(xiàn)G-2含有多種腸毒素及溶血素基因，如溶血素BL、非溶血性腸毒素基因A、腸毒素C等，還含有InhA蛋白（即免疫抑制因子A）、磷脂酶Phospholipase C precursor和蛋白質毒素Perfringolysin O precursor。溶血素通過溶解紅細胞損害宿主，InhA蛋白通過水解宿主的抗細菌蛋白來增強自身的毒性，磷脂酶則通過水解磷脂酞膽堿、磷脂酞絲氨酸和磷脂酞乙醇破壞宿主腸道，蛋白質毒素Perfringolysin O的編碼基因主要存在于產(chǎn)氣夾膜梭菌中，具有膜穿孔的能力[24]。

2.4 G-2系統(tǒng)發(fā)育樹分析

由圖8可知，G-2與蘇云金芽孢桿菌IAM 12077聚為一支，自展值為61%，因此可以確定菌株G-2為蘇云金芽孢桿菌。

3 討論與結論

隨著水產(chǎn)養(yǎng)殖的集約化發(fā)展，高密度養(yǎng)殖和抗生素的大量使用造成了養(yǎng)殖魚類的細菌性病害日益嚴重。黃志堅等[25]研究發(fā)現(xiàn)蘇云金芽孢桿菌與凡納濱對蝦肝胰腺壞死癥暴發(fā)密切相關，凡納濱對蝦感染蘇云金芽孢桿菌后，其鰓和腸道的Lv-JAK基因和Lv-STAT基因的相對表達量均顯著上調，而JAK基因和STAT基因與天然免疫密切相關。楊灣等[26]研究表明，蘇云金芽孢桿菌能夠感染甲魚，造成甲魚脖子和脾臟腫脹并使其肝臟呈土黃色。

該研究中蘇云金芽孢桿菌G-2會造成養(yǎng)殖沙塘鱧腸道出血，鰓部及頭部潰爛[27]，對G-2進行全基因組測序，獲得了一條4 400 748 bp的環(huán)狀染色體基因組，總基因數(shù)為5 134個，其中蛋白編碼基因

4 977個，非蛋白編碼基因157個，非蛋白編碼基因中tRNA基因107個，rRNA基因42個，sRNA基因8個。rRNA中23S_rRNA基因14個，16S_rRNA基因1 4個，5S_rRNA基因14個。Schnepf等[7]發(fā)現(xiàn)蘇云金芽孢桿菌染色體為環(huán)狀DNA，GC含量一般為32%～35%，基因組大小2.4～5.7 Mb不等，與該研究測序得到的結果相似。

蘇云金芽孢桿菌G-2有2 543個基因在KEGG數(shù)據(jù)庫中得到注釋，代謝類的子條目和基因數(shù)量都是最多的，其中氨基酸代謝通路和碳水化合物代謝通路分別有270和249個基因?？傮w而言，菌株G-2代謝活動比較豐富，需要與外界進行能量交換來維持自身生命活動。GO注釋結果表明，菌株G-2蛋白功能主要集中在生物過程中的代謝過程和分子功能中的催化活性功能。蛋白質COG分析結果表明G-2碳水化合物代謝和氨基酸的轉運活動較多，表明該菌株的基礎代謝較高、生理活動比較豐富。VFDB注釋可以分析不同病原菌的毒力因子構成以及基因組的分布特點，從基因水平上解釋病原菌的入侵機制。通過數(shù)據(jù)庫對比發(fā)現(xiàn)，G-2含有多種腸毒素和溶血素基因，這也解釋了為什么感染該菌株的沙塘鱧的主要癥狀為腸道出血。同時該菌株還含有InhA蛋白、磷脂酶和蛋白毒素；與華中農(nóng)業(yè)大學測序的蘇云金芽孢桿菌YBT-1520比較發(fā)現(xiàn)，YBT-1520和G-2都含有腸毒素、溶血素、磷脂酶和InhA蛋白；但YBT-1520還含有蛋白酶Enhancin和幾丁質酶Chitinase，可以提高殺蟲活性，G-2則含有蛋白質毒素Perfringolysin O，能夠溶解紅細胞。

目前國內針對蘇云金芽孢桿菌的全基因組測序研究相對較少，為探明造成養(yǎng)殖沙塘鱧腸道出血的病原菌蘇云金芽孢桿菌G-2的特性，該研究對G-2進行了全基因組測序，從基因組學上分析了該菌的基因功能，并進行毒力基因分析，解釋了該菌造成養(yǎng)殖沙塘鱧腸道出血的原因，以期為沙塘鱧的病害防控提供科研基礎。

參考文獻：

[1] GOODWIN S，MCPHERSON J D，MCCOMBIE W R. Coming of age：ten years of next-generation sequencing technologies[J]. Nature Reviews Genetics，2016，17（6）：333-351.

[2] GALPERIN M Y，KRISTENSEN D M，MAKAROVA K S，et al. Microbial genome analysis：the COG approach[J]. Briefings in Bioinformatics，2019，20（4）：1063-1070.

[3] 劉貴明. 蘇云金芽孢桿菌YBT-1520全基因組測序和比較基因組學研究[D]. 杭州：浙江大學，2008.

[4] MARGUERAT S，WILHELM B T，B?HLER J. Next-generation sequencing：applications beyond genomes[J]. Biochemical Society Transactions，2008，36（Pt 5）：1091-1096.

[5] 譚壽湖，張文飛，葉大維. 蘇云金芽孢桿菌基因組研究概況[J]. 基因組學與應用生物學，2009，28（1）：202-208.

[6] BEN-DOV E. Bacillus thuringiensis subsp. israelensis and its dipteran-specific toxins[J]. Toxins，2014，6（4）：1222-1243.

[7] SCHNEPF E，CRICKMORE N，VAN RIE J，et al. Bacillus thuringiensis and its pesticidal crystal proteins[J]. Microbiology and Molecular Biology Reviews：MMBR，1998，62（3）：775-806.

[8] WAALWIJK C，DULLEMANS A M，VAN WORKUM M E，et al. Molecular cloning and the nucleotide sequence of the Mr 28000 crystal protein gene of Bacillus thuringiensis subsp. israelensis[J]. Nucleic Acids Research，1985，13（22）：8207-8217.

[9] PROMDONKOY B，ELLAR D J. Membrane pore architecture of a cytolytic toxin from Bacillus thuringiensis[J]. The Biochemical Journal，2000，350（Pt 1）：275-282.

[10] AKIBA T，ABE Y，KITADA S，et al. Crystal structure of the parasporin-2 Bacillus thuringiensis toxin that recognizes cancer cells[J]. Journal of Molecular Biology，2009，386（1）：121-133.

[11] CRICKMORE N，ZEIGLER D R，F(xiàn)EITELSON J，et al. Revision of the nomenclature for the Bacillus thuringiensis pesticidal crystal proteins[J]. Microbiology and Molecular Biology Reviews：MMBR，1998，62（3）：807-813.

[12] 喻子牛，孫明，劉子鐸，等. 蘇云金芽孢桿菌的分類及生物活性蛋白基因[J]. 中國生物防治，1996，12（2）：85-89.

[13] GAN Y Q，ZHANG T，GAN Y Q，et al. Complete genome sequences of two Enterococcus faecium strains and comparative genomic analysis[J]. Experimental and Therapeutic Medicine，2020，19（3）：2019-2028.

[14] LAGESEN K，HALLIN P，R?DLAND E A，et al. RNAmmer：consistent and rapid annotation of ribosomal RNA genes[J]. Nucleic Acids Research，2007，35（9）：3100-3108.

[15] GARDNER P P，DAUB J，TATE J G，et al. Rfam：updates to the RNA families database[J]. Nucleic Acids Research，2009，37：D136-D140.

[16] NAWROCKI E P，KOLBE D L，EDDY S R. Infernal 1.0：inference of RNA alignments[J]. Bioinformatics，2009，25（10）：1335-1337.

[17] LOWE T M，EDDY S R. tRNAscan-SE：a program for improved detection of transfer RNA genes in genomic sequence[J]. Nucleic Acids Research，1997，25（5）：955-964.

[18] BENSON G. Tandem repeats finder：a program to analyze DNA sequences[J]. Nucleic Acids Research，1999，27（2）：573-580.

[19] KANEHISA M，GOTO S，KAWASHIMA S，et al. The KEGG resource for deciphering the genome[J]. Nucleic Acids Research，2004，32（Database issue）：D277-D280.

[20] KANEHISA M，GOTO S，HATTORI M，et al. From genomics to chemical genomics：new developments in KEGG[J]. Nucleic Acids Research，2006，34（Database issue）：D354-D357.

[21] ASHBURNER M，BALL C A，BLAKE J A，et al. Gene ontology：tool for the unification of biology. The Gene Ontology Consortium[J]. Nature Genetics，2000，25（1）：25-29.

[22] CHEN L，YANG J，YU J，et al. VFDB：a reference database for bacterial virulence factors[J]. Nucleic Acids Research，2005，33（suppl_1）：D325-D328.

[23] 王媛. 獼猴桃展青霉素產(chǎn)生菌的識別分析及控制機制研究[D]. 楊凌：西北農(nóng)林科技大學，2017.

[24] HEUCK A P，HOTZE E M，TWETEN R K，et al. Mechanism of

membrane insertion of a multimeric beta-barrel protein：perfringolysin

O creates a pore using ordered and coupled conformational changes[J]. Molecular Cell，2000，6（5）：1233-1242.

[25] 黃志堅，陳勇貴，翁少萍，等. 多種細菌與凡納濱對蝦肝胰腺壞死癥（HPNS）爆發(fā)有關[J]. 中山大學學報（自然科學版），2016，55（1）：1-11.

[26] 楊灣，王琳，朱雅琳，等. 中華鱉源蘇云金芽孢桿菌的分離鑒定及毒力因子檢測[J]. 浙江萬里學院學報，2022，35（6）：85-93.

[27] 趙彥華，史楊白，劉國興，等. 一株養(yǎng)殖沙塘鱧腸道溶血性蘇云金芽孢桿菌的分離鑒定[J]. 水產(chǎn)養(yǎng)殖，2021，42（4）：16-19.

（責任編輯：王婷）