BtToxin_Digger: a comprehensive and high-throughput pipeline for mining toxin protein genes from Bacillus thuringiensis

Genomic and morphological features of an Amazonian Bacillus thuringiensis with mosquito larvicidal activity

The occurrence of mosquito-borne diseases is increasing, and their geographical range is expanding due to climate change. New control measures are urgently needed to combat these debilitating and, in some cases, fatal diseases. Bacteria of the genus Bacillus are of interest due to the production of bioactive compounds, including those useful for insect control. The discovery and characterization of new species of Bacillus with mosquito larvicidal activity may offer opportunities to develop new products for vector control. In this study, we evaluated larvicidal activity, described morphological characteristics, and sequenced and analyzed the genome of a bacterial strain (GD02.13) isolated from the Amazon region. The metabolites produced by GD02.13 are as effective in killing Aedes aegypti larvae as the commercial product Natular™ DT (Spinosad). Furthermore, the morphological characteristics of the GD02.13 spores and crystal inclusions resemble those previously described for B. thuringiensis. A phylogenetic analysis based on 443 single-copy orthologs indicated that the bacterial strain GD02.13 belongs to the Bacillus thuringiensis species. Its genome, which was assembled and has a size of 6.6 Mb, contains 16 secondary metabolite biosynthetic gene clusters and genes encoding insecticidal proteins, predicted based on sequence similarity. The data obtained in this study support the development of new insecticide products based on the strain GD02.13 of B. thuringiensis.

Assessing the authenticity and purity of a commercial Bacillus thuringiensis bioinsecticide through whole genome sequencing and metagenomics approaches

Biopesticides, biological agents for pest control in plants, are becoming increasingly prevalent in agricultural practices. However, no established methodology currently exists to assess their quality, and there are currently no publicly available authenticity and purity evaluations of commercial products. This lack of data may represent risks because of their widespread dispersal in the environment. We evaluated the potential of whole genome sequencing (WGS) and metagenomics approaches, including nanopore long-read sequencing, to verify both authenticity (i.e., the labeled strain) and biological purity (i.e., the absence of any undesired genetic material) of commercial Bacillus thuringiensis bioinsecticides. Four commercially available bioinsecticidal products containing Bacillus thuringiensis serovar kurstaki strain HD-1 were collected from the European market as a case study. Two sequencing approaches were employed: WGS of isolates and metagenomics sequencing of all genetic material in a product. To assess authenticity, isolate WGS data were compared against the publicly available reference genome of the expected strain. Antimicrobial resistance gene content, insecticidal gene content, and single nucleotide polymorphism differences were characterized to evaluate similarity to the reference genome. To assess purity, metagenomic sequencing data were analyzed using read classification and strain differentiation methods. Additionally, long- and short-read data were used to assess potential large-scale structural variations. Our results confirmed all investigated products to be authentic and pure. With the increasing usage of biopesticides, it is crucial to have adequate quality control methods. Our proposed approach could be adapted for other biopesticides, and similar products, providing a standardized and robust approach to contribute to biopesticide safety.

Genomic Insights into the Bactericidal and Fungicidal Potential of Bacillus mycoides b12.3 Isolated in the Soil of Olkhon Island in Lake Baikal, Russia

The dispersal of plant pathogens is a threat to the global economy and food industry which necessitates the need to discover efficient biocontrol agents such as bacteria, fungi, etc., inhibiting them. Here, we describe the Bacillus mycoides strain b12.3 isolated from the soil of Olkhon Island in Lake Baikal, Russia. By applying the co-cultivation technique, we found that the strain inhibits the growth of plant pathogens, such as the bacteria Xanthomonas campestris, Clavibacter michiganensis, and Pectobacterium atrospecticum, as well as the fungus Alternaria solani. To elucidate the genomic fundament explaining these activities, we leveraged next-generation whole-genome sequencing and obtained a high-quality assembly based on short reads. The isolate bore seven known BGCs (biosynthetic gene clusters), including those responsible for producing bacillibactin, fengycin, and petrobactin. Moreover, the genome contained insecticidal genes encoding for App4Aa1, Tpp78Ba1, and Spp1Aa1 toxins, thus implicating possible pesticidal potential. We compared the genome with the 50 closest assemblies and found that b12.3 is enriched with BGCs. The genomic analysis also revealed that genomic architecture corresponds to the experimentally observed activity spectrum implying that the combination of produced secondary metabolites delineates the range of inhibited phytopathogens Therefore, this study deepens our knowledge of the biology and ecology of B. mycoides residing in the Lake Baikal region.

Whole-Genome Sequencing of Peribacillus frigoritolerans Strain d21.2 Isolated in the Republic of Dagestan, Russia

Pesticide-free agriculture is a fundamental pillar of environmentally friendly agriculture. To this end, there is an active search for new bacterial strains capable of synthesizing secondary metabolites and toxins that protect crops from pathogens and pests. In this study, we isolated a novel strain d21.2 of Peribacillus frigoritolerans from a soil sample collected in the Republic of Dagestan, Russia. Leveraging several bioinformatic approaches on Illumina-based whole-genome assembly, we revealed that the strain harbors certain insecticidal loci (coding for putative homologs of Bmp and Vpa) and also contains multiple BGCs (biosynthetic gene clusters), including paeninodin, koranimine, schizokinen, and fengycin. In total, 21 BGCs were predicted as synthesizing metabolites with bactericidal and/or fungicidal effects. Importantly, by applying a re-scaffolding pipeline, we managed to robustly predict MGEs (mobile genetic elements) associated with BGCs, implying high genetic plasticity. In addition, the d21.2's genome was free from genes encoding for enteric toxins, implying its safety in use. A comparison with available genomes of the Peribacillus frigoritolerans strain revealed that the strain described here contains more functionally important loci than other members of the species. Therefore, strain d21.2 holds potential for use in agriculture due to the probable manifestation of bactericidal, fungicidal, growth-stimulating, and other useful properties. The assembled genome is available in the NCBI GeneBank under ASM4106054v1.

Hybrid de novo whole genome assembly of lipopeptide producing novel Bacillus thuringiensis strain NBAIR BtAr exhibiting antagonistic activity against Sclerotium rolfsii

Bacillus thuringiensis Berliner is recognized as a predominant bioinsecticide but its antifungal potential has been relatively underexplored. A novel B. thuringiensis strain NBAIR BtAr was isolated and morphologically characterized using light and scanning electron microscopy, revealing presence of bipyramidal, cuboidal, and spherical parasporal crystals. The crude form of lipopeptides was extracted from NBAIR BtAr and assessed for its antagonistic activity in vitro, and demonstrated 100 % inhibition of Sclerotium rolfsii Sacc. at a minimum inhibitory concentration of 50 μL of the crude lipopeptide extract per mL of potato dextrose agar. To identify the antagonistic genes responsible, we performed whole genome sequencing of NBAIR BtAr, revealing the presence of circular chromosome of 5,379,913 bp and 175,362 bp plasmid with 36.06 % guanine-cytosine content and 5814 protein-coding sequences. Average nucleotide identity and whole genome phylogenetic analysis delineated the NBAIR BtAr strain as konkukian serovar. Gene ontology analysis revealed associations of 1474, 1323, and 1833 genes with biological processes, molecular function, and cellular components, respectively. Antibiotics & secondary metabolite analysis shell analysis of the whole genome yielded secondary metabolites biosynthetic gene clusters with 100 %, 85 %, 40 %, and 35 % similarity for petrobactin, bacillibactin, fengycin, and paenilamicin, respectively. Also, novel biosynthetic gene clusters, along with antimicrobial genes, including zwittermicin A, chitinase, and phenazines, were identified. Moreover, the presence of eight bacteriophage sequences, 18 genomic islands, insertion sequences, and one CRISPR region indicated prior occurrences of genetic exchange and thus improved competitive fitness of the strain. Overall, the whole genome sequence of NBAIR BtAr is presented, with its taxonomic classification and critical genetic attributes that contribute to its strong antagonistic activity against S. rolfsii.

Insights into the whole genome sequence of Bacillus thuringiensis NBAIR BtPl, a strain toxic to the melon fruit fly, Zeugodacus cucurbitae

Bacillus thuringiensis is the most widely used biopesticide, targets a diversity of insect pests belonging to several orders. However, information regarding the B. thuringiensis strains and toxins targeting Zeugodacus cucurbitae is very limited. Therefore, in the present study, we isolated and identified five indigenous B. thuringiensisstrains toxic to larvae of Z. cucurbitae. However, of five strains NBAIR BtPl displayed the highest mortality (LC50 = 37.3 μg/mL) than reference strain B. thuringiensis var. israelensis (4Q1) (LC50 = 45.41 μg/mL). Therefore, the NBAIR BtPl was considered for whole genome sequencing to identify the cry genes present in it. Whole genome sequencing of our strain revealed genome size of 6.87 Mb with 34.95% GC content. Homology search through the BLAST algorithm revealed that NBAIR BtPl is 99.8% similar to B. thuringiensis serovar tolworthi, and gene prediction through Prokka revealed 7406 genes, 7168 proteins, 5 rRNAs, and 66 tRNAs. BtToxin_Digger analysis of NBAIR BtPl genome revealed four cry gene families: cry1, cry2, cry8Aa1, and cry70Aa1. When tested for the presence of these four cry genes in other indigenous strains, results showed that cry70Aa1 was absent. Thus, the study provided a basis for predicting cry70Aa1 be the possible reason for toxicity. In this study apart from novel genes, we also identified other virulent genes encoding zwittermicin, chitinase, fengycin, and bacillibactin. Thus, the current study aids in predicting potential toxin-encoding genes responsible for toxicity to Z. cucurbitae and thus paves the way for the development of B. thuringiensis-based formulations and transgenic crops for management of dipteran pests.

Lysinibacillus pinottii sp. nov., a novel species with anti-mosquito and anti-mollusk activity

An isolate of a Gram-positive, strictly aerobic, motile, rod-shaped, endospore forming bacterium was originally isolated from soil when screening and bioprospecting for plant beneficial microorganisms. Phylogenetic analysis of the 16S rRNA gene sequences indicated that this strain was closely related to Lysinibacillus fusiformis NRRL NRS-350T (99.7%) and Lysinibacillus sphaericus NRRL B-23268T (99.2%). In phenotypic characterization, the novel strain was found to grow between 10 and 45 °C and tolerate up to 8% (w/v) NaCl. Furthermore, the strain grew in media with pH 5 to 10 (optimal growth at pH 7.0). The predominant cellular fatty acids were observed to be iso-C15: 0 (52.3%), anteiso-C15: 0 (14.8%), C16:1ω7C alcohol (11.2%), and C16: 0 (9.5%). The cell-wall peptidoglycan contained lysine-aspartic acid, the same as congeners. A draft genome was assembled and the DNA G+C content was determined to be 37.1% (mol content). A phylogenomic analysis on the core genome of the new strain and 5 closest type strains of Lysinibacillus revealed this strain formed a distinct monophyletic clade with the nearest neighbor being Lysinibacillus fusiformis. DNA-DNA relatedness studies using in silico DNA-DNA hybridizations (DDH) showed this species was below the species threshold of 70%. Based upon the consensus of phylogenetic and phenotypic analyses, we conclude that this strain represents a novel species within the genus Lysinibacillus, for which the name Lysinibacillus pinottii sp. nov. is proposed, with type strain PB211T (= NRRL B-65672T, = CCUG 77181T).

Comparison of the performance of multiple whole-genome sequence-based tools for the identification of Bacillus cereus sensu stricto biovar Thuringiensis

The Bacillus cereus sensu stricto (s.s.) species comprises strains of biovar Thuringiensis (Bt) known for their bioinsecticidal activity, as well as strains with foodborne pathogenic potential. Bt strains are identified (i) based on the production of insecticidal crystal proteins, also known as Bt toxins, or (ii) based on the presence of cry, cyt, and vip genes, which encode Bt toxins. Multiple bioinformatics tools have been developed for the detection of crystal protein-encoding genes based on whole-genome sequencing (WGS) data. However, the performance of these tools is yet to be evaluated using phenotypic data. Thus, the goal of this study was to assess the performance of four bioinformatics tools for the detection of crystal protein-encoding genes. The accuracy of sequence-based identification of Bt was determined in reference to phenotypic microscope-based screening for the production of crystal proteins. A total of 58 diverse B. cereus sensu lato strains isolated from clinical, food, environmental, and commercial biopesticide products underwent WGS. Isolates were examined for crystal protein production using phase contrast microscopy. Crystal protein-encoding genes were detected using BtToxin_Digger, BTyper3, IDOPS (identification of pesticidal sequences), and Cry_processor. Out of 58 isolates, the phenotypic production of crystal proteins was confirmed for 18 isolates. Specificity and sensitivity of Bt identification based on sequences were 0.85 and 0.94 for BtToxin_Digger, 0.97 and 0.89 for BTyper3, 0.95 and 0.94 for IDOPS, and 0.88 and 1.00 for Cry_processor, respectively. Cry_processor predicted crystal protein production with the highest specificity, and BtToxin_Digger and IDOPS predicted crystal protein production with the highest sensitivity. Three out of four tested bioinformatics tools performed well overall, with IDOPS achieving high sensitivity and specificity (>0.90).IMPORTANCEStrains of Bacillus cereus sensu stricto (s.s.) biovar Thuringiensis (Bt) are used as organic biopesticides. Bt is differentiated from the foodborne pathogen Bacillus cereus s.s. by the production of insecticidal crystal proteins. Thus, reliable genomic identification of biovar Thuringiensis is necessary to ensure food safety and facilitate risk assessment. This study assessed the accuracy of whole-genome sequencing (WGS)-based identification of Bt compared to phenotypic microscopy-based screening for crystal protein production. Multiple bioinformatics tools were compared to assess their performance in predicting crystal protein production. Among them, identification of pesticidal sequences performed best overall at WGS-based Bt identification.

Get to Know Your Neighbors: Characterization of Close Bacillus anthracis Isolates and Toxin Profile Diversity in the Bacillus cereus Group

Unexpected atypical isolates of Bacillus cereus s.l. occasionally challenge conventional microbiology and even the most advanced techniques for anthrax detection. For anticipating and gaining trust, 65 isolates of Bacillus cereus s.l. of diverse origin were sequenced and characterized. The BTyper3 tool was used for assignation to genomospecies B. mosaicus (34), B. cereus s.s (29) and B. toyonensis (2), as well as virulence factors and toxin profiling. None of them carried any capsule or anthrax-toxin genes. All harbored the non-hemolytic toxin nheABC and sphygomyelinase spH genes, whereas 41 (63%), 30 (46%), 11 (17%) and 6 (9%) isolates harbored cytK-2, hblABCD, cesABCD and at least one insecticidal toxin gene, respectively. Matrix-assisted laser desorption ionization-time of flight mass spectrometry confirmed the production of cereulide (ces genes). Phylogeny inferred from single-nucleotide polymorphisms positioned isolates relative to the B. anthracis lineage. One isolate (BC38B) was of particular interest as it appeared to be the closest B. anthracis neighbor described so far. It harbored a large plasmid similar to other previously described B. cereus s.l. megaplasmids and at a lower extent to pXO1. Whereas bacterial collection is enriched, these high-quality public genetic data offer additional knowledge for better risk assessment using future NGS-based technologies of detection.

Genomic–proteomic analysis of a novel Bacillus thuringiensis strain: toxicity against two lepidopteran pests, abundance of Cry1Ac5 toxin, and presence of InhA1 virulence factor

Bacillus thuringiensis (Bt) is a biological alternative to the indiscriminate use of chemical insecticides in agriculture. Due to resistance development on insect pests to Bt crops, isolating novel Bt strains is a strategy for screening new pesticidal proteins or strains containing toxin profile variety that can delay resistance. Besides, the combined genomic and proteomic approaches allow identifying pesticidal proteins and virulence factors accurately. Here, the genome of a novel Bt strain (Bt TOL651) was sequenced, and the proteins from the spore-crystal mixture were identified by proteomic analysis. Toxicity bioassays with the spore-crystal mixture against larvae of Diatraea saccharalis and Anticarsia gemmatalis, key pests of sugarcane and soybean, respectively, were performed. The toxicity of Bt TOL651 varies with the insect; A. gemmatalis (LC₅₀ = 1.45 ng cm^-2) is more susceptible than D. saccharalis (LC₅₀ = 73.77 ng cm^-2). Phylogenetic analysis of the gyrB gene indicates that TOL651 is related to Bt kenyae strains. The genomic analysis revealed the presence of cry1Aa18, cry1Ac5, cry1Ia44, and cry2Aa9 pesticidal genes. Virulence factor genes such as phospholipases (plcA, piplc), metalloproteases (inhA), hemolysins (cytK, hlyIII, hblA, hblC, hblD), and enterotoxins (nheA, nheB, nheC) were also identified. The combined use of the genomic and proteomic data indicated the expression of Cry1Aa18, Cry1Ac5, and Cry2Aa9 proteins, with Cry1Ac5 being the most abundant. InhA1 also was expressed and may contribute to Bt TOL651 pathogenicity. These results provide Bt TOL651 as a new tool for the biocontrol of lepidopteran pests.

Identification of Genetic Markers for the Detection of Bacillus thuringiensis Strains of Interest for Food Safety

Bacillus thuringiensis (Bt), belonging to the Bacillus cereus (Bc) group, is commonly used as a biopesticide worldwide due to its ability to produce insecticidal crystals during sporulation. The use of Bt, especially subspecies aizawai and kurstaki, to control pests such as Lepidoptera, generally involves spraying mixtures containing spores and crystals on crops intended for human consumption. Recent studies have suggested that the consumption of commercial Bt strains may be responsible for foodborne outbreaks (FBOs). However, its genetic proximity to Bc strains has hindered the development of routine tests to discriminate Bt from other Bc, especially Bacillus cereus sensu stricto (Bc ss), well known for its involvement in FBOs. Here, to develop tools for the detection and the discrimination of Bt in food, we carried out a genome-wide association study (GWAS) on 286 complete genomes of Bc group strains to identify and validate in silico new molecular markers specific to different Bt subtypes. The analyses led to the determination and the in silico validation of 128 molecular markers specific to Bt, its subspecies aizawai, kurstaki and four previously described proximity clusters associated with these subspecies. We developed a command line tool based on a 14-marker workflow, to carry out a computational search for Bt-related markers from a putative Bc genome, thereby facilitating the detection of Bt of interest for food safety, especially in the context of FBOs.

Environmental Behaviors of Bacillus thuringiensis (Bt) Insecticidal Proteins and Their Effects on Microbial Ecology

Bt proteins are crystal proteins produced by Bacillus thuringiensis (Bt) in the early stage of spore formation that exhibit highly specific insecticidal activities. The application of Bt proteins primarily includes Bt transgenic plants and Bt biopesticides. Transgenic crops with insect resistance (via Bt)/herbicide tolerance comprise the largest global area of agricultural planting. After artificial modification, Bt insecticidal proteins expressed from Bt can be released into soils through root exudates, pollen, and plant residues. In addition, the construction of Bt recombinant engineered strains through genetic engineering has become a major focus of Bt biopesticides, and the expressed Bt proteins will also remain in soil environments. Bt proteins expressed and released by Bt transgenic plants and Bt recombinant strains are structurally and functionally quite different from Bt prototoxins naturally expressed by B. thuringiensis in soils. The former can thus be regarded as an environmentally exogenous substance with insecticidal toxicity that may have potential ecological risks. Consequently, biosafety evaluations must be conducted before field tests and production of Bt plants or recombinant strains. This review summarizes the adsorption, retention, and degradation behavior of Bt insecticidal proteins in soils, in addition to their impacts on soil physical and chemical properties along with soil microbial diversity. The review provides a scientific framework for evaluating the environmental biosafety of Bt transgenic plants, Bt transgenic microorganisms, and their expression products. In addition, prospective research targets, research methods, and evaluation methods are highlighted based on current research of Bt proteins.

BPPRC database: a web-based tool to access and analyse bacterial pesticidal proteins

Pesticidal proteins derived from the bacterium Bacillus thuringiensis, have provided the bases for a diverse array of pest management tools ranging from natural products used in organic agriculture, to modern biotechnological approaches. With advances in genome sequencing technologies and protein structure determination, an increasing number of pesticidal proteins from myriad bacterial species have been identified. The Bacterial Pesticidal Protein Resource Center (BPPRC) has been established to provide informational and analytical resources on the wide range of pesticidal proteins derived from bacteria that have potential utility for arthropod management. In association with a revised nomenclature for these proteins, BPPRC contains a database that allows users to browse and download sequences. Users can search the database for the best matches to sequences of interest and can incorporate their own sequences into basic informatic analyses. These analyses include the ability to draw and export guide trees from either whole protein sequences or, in the case of the three-domain Cry proteins, from individual domains. The associated website also provides a portal for users to submit protein sequences for naming. The BPPRC provides a single authoritative source of information to which all stakeholders can be referred including academics, government regulatory bodies and research and development personnel in the industrial sector. The database provides information on more than 1060 pesticidal proteins derived from 13 species of bacteria, including insecticidal activities for a subset of these proteins.

Database URL: www.bpprc.org and www.bpprc-db.org/

Draft Genome Sequence of Bacillus thuringiensis ZZQ-130 with Multiple Pesticidal Genes, Isolated from Caka Salt Lake, China

Bacillus thuringiensis is a typical pesticide, with global application for over 40 years. Here, we report the draft genome sequence of B. thuringiensis ZZQ-130 from a salt lake; this strain has 31 pesticidal genes, including five cry genes, one vip gene, two vpa genes, and two vpb genes.

Draft Genome Sequence of Bacillus wiedmannii Biovar thuringiensis ZZQ-138, Isolated from a Saline Lake

Ecologically sound approaches to control mosquitoes like Anopheles stephensi, which are obligatory vectors for malaria transmission, are urgently needed because of increasing insecticide resistance. Bacteria from Bacillus are important resources. Here, we report the whole-genome sequence of Bacillus wiedmannii biovar thuringiensis ZZQ-138, which was isolated from sediment from a saline lake.

IDOPS, a Profile HMM-Based Tool to Detect Pesticidal Sequences and Compare Their Genetic Context

Biopesticide-based crop protection is constantly challenged by insect resistance. Thus, expansion of available biopesticides is crucial for sustainable agriculture. Although Bacillus thuringiensis is the major agent for pesticide bioprotection, the number of bacteria species synthesizing proteins with biopesticidal potential is much higher. The Bacterial Pesticidal Protein Resource Center (BPPRC) offers a database of sequences for the control of insect pests, grouped in structural classes. Here we present IDOPS, a tool that detects novel biopesticidal sequences and analyzes them within their genetic environment. The backbone of the IDOPS detection unit is a curated collection of high-quality hidden Markov models that is in accordance with the BPPRC nomenclature. IDOPS was positively benchmarked with BtToxin_Digger and Cry_Processor. In addition, a scan of the UniProtKB database using the IDOPS models returned an abundance of new pesticidal protein candidates distributed across all of the structural groups. Gene expression depends on the genomic environment, therefore, IDOPS provides a comparative genomics module to investigate the genetic regions surrounding pesticidal genes. This feature enables the investigation of accessory elements and evolutionary traits relevant for optimal toxin expression and functional diversification. IDOPS contributes and expands our current arsenal of pesticidal proteins used for crop protection.